О незначимости энергетической оценки при моделировании РНК
Это описание эксперимента: его планирование и полученные результаты. Автор: Сергей Яковлев |
Введение
правитьВ статье «О детерминированном моделировании РНК-петли» было показано каким образом можно перейти от стохастического моделирования РНК к детерминированному. Так же там был упомянут метод „X-тюнинг“, который целенаправленно конструирует конечное состояние (нативное) РНК. При этом он использует для получения грубой структуры РНК специальным образом видоизмененную по сравнению с работой [1] оценку энергии [2].
В данном исследовании была сделана попытка при конструировании РНК, основываться на такую энергетическую оценку, в которой учтено только отсутствие ковалентных связей и близость водородных связей. Для определения водородных связей должна быть известна вторичная структура РНК.
Модельный организм
правитьВ качестве модельного организма взят вироидный рибозим NC_003540 с классической структурой типа “головки молотка” (Hamerhead ribozyme). Первичная структура фрагмента RZ+ следующая: aagaggucggcaccugacgucgguguccugaugaagauccaugacaggaucgaaaccucuu. Вторичная представлена на рисунке, а третичная неизвестна и её нужно получить.
Ход эксперимента
правитьРибозим разделяется на фрагменты в соответствии со вторичной структурой L1 (ggcaccugacgucgguguc), L2 (gauccaugacaggauc) и то что остается. Фрагменты L1 и L2 собираются поотдельности, вначале формируя петлю, а затем достраивая последующий водородные связи.
Вначале получаем грубую структуру РНК-петли используя метод „Быстрое охлаждение“. В методе используется придельно упрощенная энергетическая оценка, основанная только на расчете близости водородных связей и отсутствии ковалентных связей между атомами разных нуклеотидов (естественно, за исключением связей фосфора, образующих цепь РНК). Затем применяется метод „X-тюнинг“.
Результаты
правитьВыводы
правитьВо многом форму стабильной РНК определяют фрагменты, где идут подряд две (или более) водородные связи вида g-c + g-c.
Для отдельных односпиральных фрагментов (таких как L1 и L2) - применение сокращенной энергетической оценки оказалось достаточно. Но для многоспиральных фрагментов ориентирование только на образование водородных связей и отсутствие ковалентных недостаточно, так как в противном случае сворачивание идет по неправильному случайному пути. В результате могут быть получены различные псевдо-третичные структуры, которые не соответствуют реальной структуре РНК фрагмента.