Использовалась программа einverted из пакета EMBOSS. Лучшее предсказание показано ниже.
SEQUENCE: Score 72: 24/24 (100%) matches, 20 gaps 1 ggggt------a-tcgccaag---cggtaag-gcaccgg-attc 32 ||||| | ||| || ||| || | |||| | || 69 ccccatgctcctaagc--tt-ggagcc-tt-ac--ggcctt-ag 34
Кака видно, нормально посчитан только акцепторный стебель, в других стеблях имеются вставленные гэпы или они не совпадают вовсе.
Был использован онлайн сервис RNAfold. Результаты представлены в таблице.
Участок структуры | Позиции в структуре (по результатам find_pair) | Результаты предсказания с помощью einverted | Результаты предсказания по алгоритму Зукера |
---|---|---|---|
Акцепторный стебель | 2-7:66-71 (6 пар) | 1-6:64-69 (6 пар) | 1-6:64-69 (6 пар) |
D-стебель | 10-12:23-25 | 10-13:53-56 | 9-11:21-23 |
T-стебель | 49-52:62-65 | - | 49-53:59-63 |
Антикодоновый стебель | 26-30:40-44 (5 пар) | - | 25-29:39-43 |
Общее число канонических пар нуклеотидов | 19 | 9 | 19 |
Из таблицы видно, что алгоритм Зукера является более эфективным в предсказании. Тогда как einverted дало малоудовлетворительные результаты.
Созданные скрипты:
script.spt define.sptБудем считать полярными атомы кислорода и азота, а неполярными атомы углерода, фосфора и серы. Назовем полярным контактом ситуацию, в которой расстояние между полярным атомом белка и полярным атомом ДНК меньше 3.5Å. Аналогично, неполярным контактом будем считать пару неполярных атомов на расстоянии меньше 4.5Å.
Требуемые атомы ДНК были получены в JMol с помошью команды define. Затем, пользуясь командой select within (расстояние, множество) были получены атомы, контактирующие с данными атомами ДНК.
Контакты атомов белка с | Полярные | Неполярные | Всего |
---|---|---|---|
Остатками 2'-дезоксирибозы | 10 | 5 | 15 |
Остатками фосфорной кислоты | 57 | 15 | 72 |
Остатками азотистых оснований со стороны большой бороздки | 8 | 3 | 11 |
Остатками азотистых оснований со стороны малой бороздки | 5 | 10 | 15 |
Схема ДНК-белковых контактов на рисунке:
По рисунку видно, что аспарагин 81, аргенин 33, лизин 79 и лизин 25 контактируют с белком по два раза. Однако, остаток аспарагин 81 важнее, так как он взаимодействует с двуия отнованиями, когда лизин25 реагирует с одним и сахарофосфатным остовом, а остальные реагируют только с сахарофосфатным остовом.
Term 3 Main page