Вторичная структура тРНК и ДНК-белковые взаимодействия

Задание 1

Ссылка на блокнот Jupyter

Было проведено предсказание координат стеблей в структуре глутамил-тРНК Thermus thermophilus (PDB ID: 1n78) при помощи программы einverted из пакета EMBOSS и по алгоритму Зукера при помощи ViennaRNA (Рис. 1).

ViennaRNA tRNA prediction
Рисунок 1. Предксазание структуры рассматриваемой глутамил-тРНК при помощи ViennaRNA.

Сравнение результатов предсказания приведено в Таблице 1. Для программы einverted были использованы следующие параметры: -gap 12 -threshold 16 -match 4 -mismatch -9 -maxrepeat 20, вследствие чего (параметр maxrepeat) акцепторный стебель в полученном предсказании отсутствовал.

Таблица 1. Результаты предсказания координат стеблей тРНК при помощи различных методов.
Участок структуры findpair einverted ViennaRNA
Акцепторный стебель 1-7 (66-72) - 1-7 (66-72)
T-стебель 26-32 (38-44) 27-31 (39-43) 27-31 (39-43)
D-стебель 10-14 (21-25) 9-12 (23-26) 9-13 (22-26)
Антикодоновый стебель 49-53 (61-65) 49-53 (61-65) 49-53 (61-65)
Общее число канонических пар 20 14 20

Можно заметить, что предсказания einverted и ViennaRNA на трех из четырех стеблей совпадают и для T- и D-стеблей несколько отличаются от предсказания, полученного при помощи findpair. Однако, общее количество предсказанных канонических пар совпадает для результатов findpair и ViennaRNA, а результат einverted отличается от них только из-за отсутствия акцепторного стебля в предсказании.

При сравнении результатов предсказания было необходимо учитывать тот факт, что нумерация нуклеотидов в тРНК не сквозная, а подчиняется правилам, основанным на выравниваниях и сравнениях структур различных тРНК (Рис. 2). Так в рассматриваемой глутамил-тРНК бактерии Thermus thermophilus присутствует нуклеотид с номером 520A, но отсутствуют с номерами 517 и 547 (в файлах PDB, включающих одновременно структуру тРНК и АРСазы нумерация нуклеотидов тРНК начинается не с 1, а с 501).

tRNA nucleotides numeration
Рисунок 2. Обобщенная структура тРНК. Толстой линией обозначены позиции, нуклеотиды, расположенные на которых, чаще всего инвариантны. Овалами обозначены нуклеотиды, которые могут отсутствовать в тРНК. [1]

Задание 2

Скрипт, определяющий необходимые множества атомов

Скрипт, выводящий изображение всей структуры

Рассматриваемая структура (PDB ID: 1I3J) описывает взаимодействие ДНК-связывающего домена интронной хоуминг-эндонуклеазы I-TevI с дуплексом ДНК. Данный домен содержит мотив "цинковый палец", протяженный участок, взаимодействующий с малой бороздкой, и мотив спираль-поворот-спираль (HTH) [2]. Для этой структуры были рассмотрены полярные и неполярные контакты ДНК-белок (Табл. 2).

Таблица 2. Различные контакты ДНК-белок в рассматриваемой структуре.
Контактирующие атомы ДНК Полярные Неполярные Всего
Отстатки 2'-дезоксирибозы 10 120 130
Остатки фосфорной кислоты 24 27 51
Азотистые основания со стороны большой бороздки 0 25 25
Азотистые основания со стороны малой бороздки 3 10 13

Также была построена карта ДНК-белковых контактов при помощи nucplot (Рис. 3).

Nucplot output
Рисунок 3. Карта контактов ДНК-белок, полученная при помощи программы nucplot.

Согласно карте nucplot наибольшее число контактов с ДНК образует аминокислотный остаток N175, взаимодействующий с нуклеотидными остатками G7 и G8 через малую бороздку (Рис. 4).

DNA-protein interactions
Рисунок 4. Взаимодействия аминокислотного остатка N175 с малой бороздкой ДНК.

Стоит отметить, что nucplot определяет неполярные контакты только в том случае, когда атомы находятся на расстоянии менее 3.35 Å друг от друга, вследствие чего часть неполярных контактов теряется, поэтому для поиска наиболее важных для распознавания последовательности взаимодействий использовались контакты определенные ранее в PyMol вручную.

В представленной структуре обнаруживается малое количество полярных контактов, причем большая их часть связана с сахарофосфатным остовом (Рис. 5), остальные - с малой бороздкой ДНК (Рис. 6). Это, вероятно, является причиной толерантности нуклеазы I-TevI к последовательности [2] и отличает ее от представителей наиболее многочисленного семейства хоуминговых нуклеаз LAGLIDADG (I-TevI относится к семейству GIY-YIG), для которых хаарктерно связывание преимущественно с большой бороздкой.

Protein and DNA backbone interactions
Рисунок 5. Полярные взаимодействия аминокислотного остатка S191 с сахарофосфатным остовом ДНК.
Protein and DNA minor groove interactions
Рисунок 6. Полярные взаимодействия аминокислотного отстатка N201 с нуклеотидами C14 и G40 со стороны малой бороздки.

Однако аминокислотные остатки мотива HTH взаимодействуют с ДНК посредством неполярных взаимодействий, связываясь с 5-метильными группами тиминовых азотистых оснований (Рис. 7), вследствие чего можно предположить, что эти взаимодействия отвечают за специфичность данного домена по отношению к AT-богатым участкам.

Non-polar interactions between DNA and protein
Рисунок 7. Взаимодействие аминокислотных остатков Y231 и L228 с 5-метильными группами азотистых оснований нуклеотидов T33 и T34.

ИСТОЧНИКИ

  1. Mathias Sprinzl, Carsten Horn, Melissa Brown, Anatoli Ioudovitch, Sergey Steinberg, Compilation of tRNA sequences and sequences of tRNA genes, Nucleic Acids Research, Volume 26, Issue 1, 1 January 1998, Pages 148–153, https://doi.org/10.1093/nar/26.1.148

  2. Roey PV, Belfort M, Derbyshire V. Homing Endonuclease I-TevI: An Atypical Zinc Finger with a Novel Function. In: Madame Curie Bioscience Database [Internet]. Austin (TX): Landes Bioscience; 2000-2013. Available from: https://www.ncbi.nlm.nih.gov/books/NBK6328/