Предсказание вторичной структуры тРНК

Задание 1

Мной была проанализированна тРНК 1FLU. С помощью нескольких программ (find_pair, einverted, ViennaRNA) сгенерированна вторичная структура тРНК на основе ее последовательности и создана сравнительная таблица.
Таблица 1. Сравнение выдачи программ, предсказывающих вторичную структуру тРНК

Участок структуры Позиции в структуре (по результатам find_pair ) Результаты предсказания с помощью einverted Результаты предсказания по алгоритму Зукера
Акцепторный стебель (1-7) - (72-66) (1-7) - (72-66) (1-5) - (72-68)
D-стебель (10-13) - (25-22) (6-12) - (17-23)
T-стебель (49-53) - (65-61) (25-28) - (64-67)
Антикодоновый стебель (39-44) - (31-26) (36-40) - (52-56)
Общее число канонических пар нуклеотидов 18 7 27

Из приведенной таблицы можно сделать вывод, что наиболее верный и точный результат выдает программа find_pair,так как она работает с проверенной pdb структурой, программа einverted предсказала лишь координаты акцепторного стебелька, а при изменении параметров (были в итоге применены: gap=3, thershold=3, match=1, mismatch=-1) выдавал длинные последовательности с большим количеством гэпов.
Программа ViennaRNA успешно предсказала акцепторный стебелек и D-петлю, однако в антикодоновом- и T-стебле неверно были построены водородные связи, более правильным были бы пары нуклеотидов (25-28) - (40-36) для антикодоновой и (52-56) - (67-64) для T-стебля.
Также был найден набор субптиальных структур, всего таковых получилось 234, одной из структур сильнее похожей на реальную вторичную структуру была выдача №215.

2. Поиск ДНК-белковых контактов в заданной структуре

Скрипт в JMol (с выделением 3х групп атомов).

Скрипт в JMol (дает последовательное изображение всей структуры только ДНК в проволочной модели с выделенными множествами).

Таблица 2. Предсказанные контакты белка с разными участками ДНК.

Контакты атомов белка с Полярные Неполярные Всего
остатками 2'-дезоксирибозы 8 39 47
остатками фосфорной кислоты 88 19 107
остатками азотистых оснований со стороны большой бороздки 21 20 41
остатками азотистых оснований со стороны малой бороздки 12 0 12
остатками азотистых оснований (неуточнённые по бороздке) 0 4 4

С помощью nucplot было получено изображение контактов ДНК-белок.

Рисунок 5

Рис. 1. Схема, полученная с помощью nucplot

Аминокислотный остаток с наибольшим числом контактов - Arg84 (2 контакта с остатками фосфорной кислоты).

Рисунок 5
Рис. 2. Аминокислота Arg84 с самым большим количеством контактов (справа видны контакты с ДНК).

Arg69 наиболее важный для распознавания последовательности ДНК, так как он взаимодействует именно с азотистым основанием (гуанином).

Рисунок 5
Рис. 3. Связь Arg69 с гуанином.