Получение файлов со структурами, их изучение.
Заданный файл PDB - 1p8k был получен с сайта PDB (единственная выдача -
файл 1p8k.pdb); а биологическая единица
изучаемого комплекса - с сайта Macromolecular
Structure Database EBI (файл 1p8k_mmol.pdb).
Следующий этап - изучение полученных структур средствами RasMol. Результаты: различий между структурами не обнаружено.
Имеются различия в размерах и, соответственно, содержании PDB-файлов: 1p8k.pdb содержит больше информации об "отдельных"
компонентах - ДНК и белке (в т.ч. поля REMARK; о последовательностях ДНК и белка - поле SEQRES; о модифицикации аминокислот - поле MODRES (также в HETATM; есть и брутто-формула - поле FORMUL);
о вторичной структуре белка (шпильки, стренды)- поля HELIX, SHEET и т.д. - не в том смысл, чтобы перечислить все поля файла PDB).
1p8k_mmol.pdb - тут информация, касающаяся комплекса ДНК-белок вцелом (без "идентификации" отдельных молекул).
О функциональных особенностях белка (а это интрон-кодируемая эндонуклеаза I-ANII)- см. п. 4.
В соответствии с заданием, дальнейшая работа ведется с файлом с биологической единицей.
Исследование контактов между молекулами белка и нуклеиновой кислоты.
Как в тексте полученных документов названы нужные атомы:
- атомы остатка сахара - O5*, C5*, C4*, O4*, C3*, O3*, C2*, C1*;
- атомы остатка фосфорной кислоты - O1P, OP2, P;
- атомы азотистых оснований (на примере С) - N1, C2, O2, N3, C4, N4, C5, C6.
Необходимо также определить множества полярных и неполярных атомов ДНК. Критерии: полярными атомами считаются все атомы кислорода и азота, а неполярными - все атомы углерода.
Тогда:
- pol_dr - полярные атомы 2'-дезоксирибозы;
- pol_ph - полярные атомы остатков фосфорной кислоты: " *.O?P "
(это два атома кислорода: O1P и O2P);
- pol_bs_mjg (pol_bs_mig) - полярные атомы азотистых оснований, обращенные в большую (малую)
бороздку;
- npolar_drнеполярные атомы 2'-дезоксирибозы: " *.C?* "
(эквивалентно:
"dna and backbone and carbon", поскольку все атомы углерода остова ДНК
находятся в остатке сахара);
- npolar_bs_mjg (npolar_bs_mig)неполярные атомы ... бороздки: "carbon and ...".
Был выдан скрипт dna.def. В нем определены еще 3 множества:
base - все атомы оснований ДНК
(эквивалентно dna and not backbone);
mjg - атомы оснований ДНК, обращенные в большую бороздку
("major groove");
mig - атомы оснований ДНК, обращенные в малую бороздку
("minor groove").
При поиске полярных контактов использован порог расстояния в 3,5 Å, а при поиске неполярных - 4,5 Å (cоответственно, полярный контакт, или водородная связь,
имеет место, если оба контактирующих атома полярны; неполярный контакт, или гидрофобное взаимодействие, предполагает, что оба атома неполярны).
Для выполнения поставленной задачи - исследования контактов между молекулами белка и ДНК - составлен данный скрипт ("2 в 1": определение
необходимых множеств + выделение исследуемых окрестностей; в поле echo дается "расшифровка" определений).
Результаты оформлены в таблицу "Число атомов белка, участвующих во взаимодействии с ...":
|
Полярные |
Гидрофобные |
Всего |
Контакты белка с ... |
... остатками 2'-дезоксирибозы |
8 |
58 |
66 |
... остатками фосфорной кислоты |
29 |
34 |
63 |
... остатками азотистых оснований со стороны большой
бороздки |
11 |
46 |
57 |
... остатками азотистых оснований со стороны малой бороздки |
0 |
0 |
0 |
Итак:
- число полярных контактов с остатками фосфорной кислоты больше числа полярных контактов с остатками 2'-дезоксирибозы более чем в 3 раза;
это согласуется с функцией белка (режет остов ДНК по фосфатным группам);
- сравнительно велико число гидрофобных взаимодействий "белок - ДНК" (в сумме 92 для сахарофосфатного остова и 46 для азотистых оснований); возможно, это необходимо
для надежной "фиксации" ДНК и белка (необходимость при резке ДНК);
общее число контактов (и полярных, и гидрофобных) весьма велико - 186 (опять же, это связано с функцией - расщепить определенные уч-ки ДНК);
- исследуемый белок взаимодействует с большой бороздкой ДНК (целых 11 полярных контактов в сравнении с 0 для малой бороздки;
46 гидрофобных контактов); вообще не контактирует с остатками азотистых оснований со стороны малой бороздки (ни гидрофобных, ни полярных взаимодействий нет). Понятно,
что эндонуклеаза для узнавания ДНК располагается со стороны большой бороздки (в нее и азотистых оснований, в принципе, обращено больше).
Поиск специфических контактов, обеспечивающих узнавание сайта в молекуле ДНК.
Задание: "Предположить, какие остатки белка отвечают за специфичность узнавания белком ДНК. ".
Для выполнения поставленной задачи необходимо посмотреть, какие аминокислотные остатки контактируют с атомами азотистых оснований большой бороздки ДНК (
образуя, соответственно, полярное и неполярное взаимодействия). Это нетрудно сделать, взяв за основу предложенный ранее скрипт.
Cписок аминокислотных остатков, образующих полярные контакты (в скобках указаны соответствующие атомы) - следовательно, отвечающих за специфичность узнавания белком ДНК:
- Glu 35 (OE1)
- Arg 59 (NH1)
- Arg 61 (NH1, NH2)
- Arg 70 (NH1, NH2)
- Arg 72 (NH1, NH2)
- Tyr 154 (OH)
- Lys 200 (NZ)
- Lys 202 (NZ)
Гидрофобные взаимодействия образуют (число контактирующих атомов обозначено *):
- Glu** 35
- Arg* 59 (NH1)
- Arg* 61 (NH1, NH2)
- Tyr****** 154 (OH)
- Lys* 202 (NZ)
Наиболее интересный, на мой взгляд, контакт Arg 61 (2 полярных взаимодействия + 1 гидрофобное (атом CZ) - "заякоренный а.о.").
Этот контакт достоин рассмотрения, т.к.
Средствами RasMol построено изображение выбранного остатка и контактирующего с ним фрагмента ДНК в шарнирной модели:
![](contact.GIF)
На этой картинке атомы, образующие рассматриваемые полярные взаимодействия, соединены пунктирными линиями. Атомы, образующие гидрофобные контакты -
атомы углерода - изображены серыми шариками.
Скрипт, генерирующий данное изображение, можно увидеть здесь.
В статье Structural and biochemical analyses of DNA and RNA binding by a bifunctional homing endonuclease and group I intron splicing factor (Jill M. Bolduc et al., 2003) я нашла подтверждение
своему предположению: все 4 указанных а.о. arg - Arg 59, Arg 61, Arg 70, Arg 72 - отвечают за специфичность узнавания белком такой пуриновой последовательности ДНК: 5'-GGAGG-3'.
Статья доступна на этом ресурсе.
Описание функций исследованного белка.
Найдено описание заданного белка в UniProt и описание соответствующего ДНК-связывающего домена в Pfam (была выбрана ссылка PF00961 в поле PFAM - это и есть ДНК-связывающий домен).
О белке: он совмещает 2 функции: это интрон-кодируемая ДНК эндонуклеаза I-AniI и мРНК матураза bI1.
- С одной стороны, белок разрезает ДНК в месте соединения 2х экзонов cobA гена (в безинтронных генах), способствуя тем самым вставке интрона, содержащего последовательность
данного белка (группа I интронов; РНК-копии таких участков ДНК являются рибозимами - каталитически активными молекулами РНК; способны к аутосплайсингу
). Этот процесс еще называют "homing'ом" интрона (а эндонуклеазы - 'homing' эндонуклеазами (см. пояснения далее)). Узнавание специфично (последовательности, по которым происходит узнавание ,большие), "резка" происходит по последовательности 5'-GAGGAGGTTTCTCTGTA-3'. Действие - на митохондриальную ДНК. В этом роль белка как эндонуклеазы.
- С другой стороны, белок катализирует процесс вырезания интронов из пре-мРНК и осуществляющий "созревание" мРНК (стабилизация вторичной/четвертичной структур мРНК). В этом его функция как матуразы.
- Cайты узнавания РНК и ДНК и каталитические сайты разобщены (независимы функции эндонуклеазы и матуразы).
- Кофактором является ион Mg2+ - его наличие необходимо для I-AniI эндонуклеазной активности.
- Белок относится к LAGLIDADG семейству эндонуклеаз.
- Гомодимер.
- Дислокация в клетке - в митохондриях.
- Организм - "источник" - Emericella nidulans (Aspergillus nidulans).
ДНК-связывающий домен - домен LAGLIDADG (на полной доменной структуре по данным Pfam изображен красными блоками):
![](ANI1_EMENI.png)
Еще одна "визуализация" (теперь изображены только домены LAGLIDADG):
![](1p8k.png)
![](LAGLIDADG.png)
Эндонуклеазы семейства LAGLIDADG характеризуются способностью "встраивать" в ДНК последовательности ("мобильные элементы" - интроны). 'Homing' - перенос "вставочного" элемента (последовательности) в близкородственную аллель, в которой этого элемента нет.
Процесс этот инициируется в т.ч. сайт-специфичными эндонуклеазами семейства LAGLIDADG, кодируемыми открытыми рамками считывания (ORF) в пределах мобильных элементов. Белки этого семейства имют уникальную стратегию узнавания очень длинных сайтов ДНК (15-30 bp), что позволяет свести к минимум "случайную" резку в пределах гена.
B Pfam было получено множественно выравнивание (оба типа: "сокращенное" - из 91 последовательности и "расширенное" (полное) - из 508 последовательности). Выравнены последовательности доменов, а не всего белка.
Значит номер 61 (по последовательности, см. файл 1p8k_mmol.pdb) аминокислотного остатка Arg 61 не будет соответствовать 61ой позиции выравнивания.
Рассматриваем "сокращенный" вариант выравнивания. И вот что получается: среди консервативных а.о. не оказывается рассмотренного ранее Arg 61. Причем даже не требовалось полной консервативности!!!! Далее открылась очень интересная вещь: когда я попыталась в выравнивании
в строке, соответствующей выравниванию рассматриваемого белка, найти целиком участок VSFRKR (это по нумерации pdb остатки с 56 по 61), потерпела неудачу. Так был найден только участок, отмеченный желтым цветом (в нумерации pdb остатки с 56 по 59). А этого участка - совершенно иная последовательность, нежели все в том же pdb. Тогда я вспомнила об обстоятельстве, смутившем меня
при получении информации из UniProt: в единственном(!) файле-выдаче UniProt содержится информация о белке в 488 аминокислотных остатка - предшественнике интрон-кодируемой ДНК эндонуклеазы I-AniI; в файле содержится цитохром b.
Но вот, опять же, с полной доменной структурой я разобралась выше: cytochrome b соответствует зеленому блоку на схеме, т.е. перед LAGLIDADG доменами (т.о.
его включение не должно повлиять на выравнивание и, тем более, "разобщение" pdb-шных участков 56-59 и 60-61).
Asp61 не оказался консервативным. Это может быть связано только с названной выше проблемой. Хотя... еще одно объяснение (разумное ли?:-)): как я уже говорила, белки данного семейства кодируются ORF'ами,
имеются ограничения на их длину, позицию. М.б. поэтому не консервативен Asp61. Но это как-то маловероятно.
Построение схемы контактов белка с ДНК.
Программа nucplot предназначена для визуализации контактов между ДНК и белком. Запускается на сервере kodomo-count
с использованием команды:
nucplot 1p8k_mmol.pdb
В выдаче получено несколько файлов. Графическое изображение контактов содержится в файле nucplot.ps (всего 4 страницы). Привожу вторую страницу:
![](contact_mmol.jpeg)
Изучив данную схему, видно следующее:
- общее число контактов значительно меньше общего числа контактов, указанных в таблице. Это можно объяснить тем, что алгоритм nucplot предъявляет более жесткие требования к определению взаимодействий (полярное/неполярное), нежели расстояние в 3.5/4.5 ангстрем между одновременно полярными/неполярными атомами.
Особенно меньше число гидрофобных взаимодействий;
- для особенно интересующей меня аминокислоты - Asp61 не оказалось только гидрофобного взаимодействия, в то время как оба полярных контакта есть (для Lys 202, например, полярных взаимодействий программа не выявила); а значит действительно можно предполагать, что она отвечает за специфичность узнавания белком ДНК.
Но, несмотря на превосходные результаты работы программы nucplot, результаты "ручного поиска" не стоит игнорировать: вместе они дают нам целостную
картину взаимодействия белка и ДНК (так, nucplot хорош для определения полярных контактов; гидрофобных же находится меньше (особенности алгоритма).
Но при взаимодействии ДНК и эндонуклеазы их значение велико - см.п.2 "Исследование контактов между молекулами белка и нуклеиновой кислоты").