Комплексы ДНК-белок
Поиск ДНК-белковых контактов в структуре 1DDN
В JMol с помощью команды define
можно задавать множества атомов. Определим следующие множества:
Таблица 1. Контакты разного типа в комплексе 1DDN.pdb
Контакты атомов белка с |
Полярные |
Неполярные |
Всего |
остатками 2'-дезоксирибозы |
7 |
31 |
38 |
остатками фосфорной кислоты |
46 |
26 |
72 |
остатками азотистых оснований со стороны большой бороздки |
7 |
1 |
8 |
остатками азотистых оснований со стороны малой бороздки |
0 |
1 |
1 |
Программа nucplot, предназначенная для визуализации контактов между ДНК и белком, запускается на сервере kodomo. К сожалению, программа работает только со старым форматом PDB, поэтому используем программу remediator для перевода.
В результате выполенения программы будут созданы файлы:
Также программа выводит количество цепей ДНК(2) и белка(8) в структуре, при запуске в самом терминале.
Но нас интересует файл nucplot.ps (ps = формат Postscript).
Откроем его ассоциированной программой
GSview. Мы получим
изображение популярной схемы ДНК-белковых контактов.
P.S. Удобнее открывать формат ps оказалось онлайн-вьювером.
Рис.1 Популярная схема ДНК-белковых контактов, полученная с помощью программы nucplot
![]() |
![]() |
![]() |
![]() |
Из рис.1 видно, что аминокислотный остаток с наибольшим числом указанных на схеме контактов с ДНК - Arg29 (3 контакта)
Аминокислотный остаток наиболее важный для распознавания последовательности ДНК должен взаимодействовать
с азотистым основанием. Gln43:С образует целых два таких взаимодействия.
Предсказание вторичной структуры 1G59 тРНК
Программа einverted из пакета EMBOSS позволяет найти инвертированные участки в нуклеотидных последовательностях.
На вход программе дается файл с последовательностью тРНК в fasta-формате.
Найдем возможные комплементарные участки в последовательности структуры 1G59 тРНК.
Сравним результаты с данными, полученными
ранее с помощью find_pair
(см. 3 семестр -> A- и В- формы ДНК. Структура РНК).
Результаты сравнения занесем в таблицу 2, приведенную ниже.
Примечание: * - это вводимый пользвателем параметр. Если не вводить свои значения, а просто нажимать Enter
то программа идет со стандартными настройками.
При стандартных установках программа не находила комплементарных пар, но задав Minimum score threshold = 15,
был получен файл seq.inv. Из файла видно, что найден только акцепторный стебель,
причем не весь и немного смещенный, по сравнению с результатами, найденнными программой find pair.
Найдены спаренные основания 3-7, 65-69, а программой find pair найденны 1-7, 66-72.
Программа mfold реализует алгоритм Зукера. Воспользуемся web вариантом
ссылка.
Параметр P - единственный, который имеет смысл менять. Он указывает, на сколько процентов выдаваемое предсказание структуры
может отличаться по своей вычисленной энергии от оптимального.
Чем больше значение этого параметра, тем больше вариантов предсказания будет выдано.
Были получены результаты при p = 10.
Таблица 2. Контакты разного типа в комплексе 1DDN.pdb
Участок структуры |
Позиции в структуре (по результатам find_pair) |
Результаты предсказания с помощью einverted |
Результаты предсказания по Zucker алгоритму |
Акцепторный стебель |
5'-501 - 507-3' |
предсказано 5 пар из 7 реальных |
предсказано 7 пар из 7 реальных |
D-стебель |
5'-510 - 512-3' |
предсказано 0 пар из 3 реальных |
предсказано 5 пар из 3 реальных |
T-стебель |
5'-549 - 553-3' |
предсказано 0 пар из 5 реальных |
предсказано 5 пар из 5 реальных |
Антикодоновый стебель |
5'-538 - 544-3' |
предсказано 0 пар из 7 реальных |
предсказано 5 пар из 7 реальных |
Общее число канонических пар нуклеотидов |
19 |
5 |
19 |
© Nuzhdina Ekaterina, 2013