Документ взят из кэша поисковой машины. Адрес оригинального документа : http://kodomo.cmm.msu.ru/FBB/year_02/science_1/ivanova.doc
Дата изменения: Thu Nov 13 15:44:25 2003
Дата индексирования: Tue Oct 2 18:09:59 2012
Кодировка: koi8-r

Эволюция пространственных структур гомеодоменов.
Иванова Д.П.
Факультет биоинженерии и биоинформатики, Московский государственный
университет им. Ломоносова, 119899
Тьюторы: Алексеевский А.В. , Спирин С.А.


Аннотация
В работе проводилось построение дерева по пространственным структурам
набора гомеодоменов, его исследование и сравнение с деревом по
аминокислотным последовательностям того же набора .


Введение

Одной из важных задач биоинформатики является сравнительная
классификация и изучение эволюции пространственных структур белков, а
также их комплексов с ДНК, РНК, другими лигандами. Это связано с тем , что
функции белков определяются их пространственным строением, и белки,
имеющие сходное строение выполняют, как правило, и сходные функции. При
этом сравнение третичных структур более важно для изучения функций белков,
чем первичных. Так, известен ряд семейств белков, которые выполняют
сходные функции и имеют сходное строение, но аминокислотные
последовательности членов семейства слабогомологичны.
Гомеодомены представляют собой высококонсервативные ДНК-связывающие
домены, входящие в состав многих эукариотических факторов транскрипции.
Белки, содержащие гомеодомен (гомеобелки) участвуют в регуляции
центральных процессов эмбриогенеза и гомеостаза. Замечательно, что
гомеодомены обнаружены во всех таксонах (где проводился поиск), и изучение
их эволюции представляет чрезвычайный интерес.
Гомеодомены принадлежат к большому классу белковых доменов, содержащих
необходимый для контакта и узнавания ДНК структурный элемент спираль-
поворот-спираль (helix-turn-helix, hth).
При этом они обладают очень высокой консервативностью как по первичной,
так и по вторичной, и по третичной структурам. В состав гомеодомена входит
N-концевой наиболее вариабельный по структуре участок («рука») и три (-
спирали. Первая спираль и следующая за ней вторая располагаются
антипараллельно и почти под прямым углом к третьей спирали, которая
собственно является узнающей и взаимодействует со специфической
последовательностью ДНК в большой бороздке. Структурному мотиву HTH
соответствуют вторая и третья спирали. Следует заметить, что при
образовании комплексов с ДНК происходят конформационные перестройки
гомеодомена.
Методами рентгеноструктурного анализа и ЯМР установлена структура ряда
изолированных гомеодоменов и их комплексов с ДНК, с некоторыми из которых
проводилась данная работа.

Материалы и методы.
Гомеодомены отобраны из представленных в PDB и NDB по следующим
критериям: пространственная структура в комплексе с ДНК, разрешение при
рентгеноструктурном анализе ( 2. 5 (. В каждой цепочке рассматривались 408-
456 а.о., кроме 1E3O, где 108-156 а.о. (все последовательности были
приведены к единой нумерации).
Для проведения автоматической классификации пространственных структур
гомеодоменов, представленных в таблице 1, использовался Web-сервис
LifeCore, внедрённый на сервере НИИФХБ им.А.Н.Белозерского МГУ
(http://www.genebee.msu.ru:8080/CoreProc/). На вход программе подается
скрипт, содержащий множественное выравнивание рассматриваемых структур по
ряду определенных атомов. На выходе выдается геометрическое ядро, т.е.
подмножество индексов атомов такое, что соответствующие атомы одинаково
расположены относительно друг друга во всех структурах, и матрица
сравнений структур - матрица размера М(М для М рассматриваемых структур ,
в которой элемент еij вычисляется по правилу: при рассмотрении только i-
го и j-го файлов, их геометрическое ядро будет на еij больше, чем для
всего семейства.
В результате обработки LifeCore данных о пространственных структурах
гомеодоменов были получены данные о геометрическом ядре и матрица
сравнений структур, которые отражают общую картину взаимного родства
структур.
На основе матрицы сравнений была составлена матрица попарных расстояний,
элементы которой получены по формуле dij=1-eij/m , где m - максимальное
число из матрицы сравнений, равное разности количества всех
рассматриваемых атомов гомеодомена и количества атомов геометрического
ядра для всех рассматриваемых структур. По полученной матрице с помощью
алгоритма UPGMA (программа eneighbor пакета EMBOSS) было построено
укорененное дерево.
Затем по аминокислотным последовательностям рассматриваемых гомеодоменов
было сделано множественное выравнивание с помощью команды emma пакета
EMBOSS. По выравниванию было построено дерево аналогично предыдущему (для
построения исходной матрицы по выравниванию исползовалась команда
eprotdist). Оба дерева представлены на рис.1.

Результаты и обсуждение.
Дерево по пространственным структурам гомеодоменов было построено по
матрице, выданной программой Lifecore. В этой матрице для каждой пары
структур фиксируется число С?-атомов в геометрическом ядре. Геометрическое
ядро определяется на основе относительного пространственного положения С?-
атомов в полипептидной цепи. Его размер для каждой пары характеризует
степень сходства пространственных структур, и чем больше число С?-атомов в
ядре, тем более схожи (нередко даже полностью) эти структуры. Как видно,
построение дерева сделано только на основе сходства пространственных
структур и никак не связано с аминокислотными последовательностями.
Дерево по последовательностям строится стандартными методами на основе
выравнивания (см. Материалы и методы).



Сравнение двух деревьев показывает следующее:
1. Все гомеодомены с одинаковой аминокислотной последовательностью в
дереве по пространственным структурам также выделяются в отдельные
ветви. При этом расстояния между листьями таких ветвей маленькие за
исключением белка Pit-1 (1AU7 A, B).
2. Семейства гомеодоменов - POUh , TALE, Antennapedia - выделяются в
отдельные ветви как последовательностям, что естественно, так и по
пространственным структурам.
3. В обоих деревьях одинаково наблюдается ветвь (назовем ее V),
содержащая семейства Eve (1JGG A,B), Antp (9ANT A,B), En (3HDD A) и
MSX (1IG7 A).
4. Основное топологическое различие между деревьями состоит в том, что в
дереве последовательностей Prd кластеризуется с ветвью V, а в дереве
структур - с ветвью POUh.

(соответствия между приведенными названиями семейств и белков и
использованными при построении деревьев PDB-кодами приведены в
таблице 2.)

Тем самым, семейства гомеодоменов могут быть выделены не только на
основе особенностей аминокислотных последовательностей, но и по наборам С?
-атомов, одинаково расположенных в пространстве. Сходство между
определенными обоими способами наборами семейств говорит о зависимости
между последовательностями и тонкими особенностями пространственных
структур.

Выводы.
На примере гомеодоменов показано хорошее соответствие между
филогенетическими деревьями, построенными на основе качественно разных
данных - аминокислотных последовательностей и пространственных структур.
Таким образом показана перспективность подхода, основанного на выделении
геометрических ядер, к изучению эволюции структур.

















Таблица 1.
|Название |Организм |PDB-код |Цепь |Название |
| | | | |семейства |
|Antennapedia(Ant|D.melanogaster |9ANT |A,B | |
|p) C39S | | | | |
| | | | |Antennapedia |
| | | | |(Antp) |
|Ultrabitorax |D.melanogaster |1B8I |A | |
|(Ubx) C39S | | | | |
|HoxB1 |Человек |1B72 |A | |
|Even-skipped(Eve|D.melanogaster |1JGG |A,B |Eve |
|) | | | | |
|Engrailed(En) |D.melanogaster |3HDD |A |En |
|Paired(Prd) |D.melanogaster |1FJL |A,B,C |Prd |
|Mata1 | |1AKH |A |Mata1 |
| |Дрожжи | | | |
| | |1YRN |A | |
| | |1AKH |B |Mat(2 |
|Mat(2 | | | | |
| |Дрожжи | | | |
| | |1MNM |C | |
| | |1YRN |B | |
|Pbx1 |Человек |1B72 |B | |
| | | | |Tale |
|Extradentical(Ex|D.melanogaster |1B8I |B | |
|d) | | | | |
|Oct-1 |Человек |1E3O |C | |
| | | | |POUh |
|Pit-1 |Крыса |1AU7 |A,B | |
|Msx-1 |Мышь |1IG7 |A |Msx-1 |