Восстановление предкового состояния доменной архитектуры

Упр.1 Восстановление предкового состояния доменной архитектуры

Для работы был выбран домен Integrin_alpha2 (Pfam: PF08441). Этот домен часто входит в состав белков, участвующих в клеточной адгезии либо во взаимодействии клеток с внеклеточным матриксом. Также домен входит в суперсемейство IG-like доменов, третичная структура которых похожа на иммуноглобулин.

Выбранный домен входит в состав 51 архитектуры. В 29 архитектурах из 51 помимо Integrin_alpha2 присутствует домен FG_GAP (часто в повторах). Ещё один домен, часто встречающийся в архитектурах вместе с выбранным доменом (17 из 51) - это VWA (von Willebrand factor type A domain). VWA домены в интегринах опосредуют адгезию, контролируемую ионами металлов. FG-GAP мотивы часто присутствуют на N-концах интегринов, формируя структуру бета-пропеллера. Их предполагаемая функция - связь с ионами кальция.

Было выбрано 3 архитектуры:

• 1. FG-GAP, Integrin_alpha2

  

  Эта архитектура встречается в 181 белке, многие из которых относятся к Uncharacterized protein в Uniprot. Однако выравниваются они с интегрином альфа 5, который представляет из себя рецептор фибронектина.
• 2. VWA, FG-GAP, Integrin_alpha2, Integrin_alpha

  

  Эта архитектура встречается в 104 белках. Помимо ранее описанных доменов в архитектуру входит цитоплазматический домен Integrin_alpha, который отвечает за передачу сигнала внутрь клетки. Белки с этой архитектурой участвуют в контактах между моноцитами, макрофагами и гранулоцитами и опосредуют эндоцитоз.
• 3. FG-GAP x 2, Integrin_alpha2

  

  Архитектура встречается в 80 белках, среди которых есть аннотированные белки, относящиеся к классу интегринов альфа-8/бета-1. Эти белки координируют прикрепление скелетных миобластов и взрослых мышечных волокон к ламинину. Во время дифференцировки он может вызывать изменения в подвижности и форме миобластов.

Сравнение будет проводиться внутри таксона Metazoa (царство животных). В этом таксоне были выбраны подтаксоны Chordata и Arthropoda.

Упр.2 Получение информации об архитектуре всех последовательностей с данным доменом

python swisspfam-to-xls.py -z /srv/databases/pfam/swissprot.gz -p PF08441 -o pfam.xls

python uniprot_to_taxonomy.py -i seqtext -o taxonomy

Получение полного выравнивания для домена и формирование выборки для архитектур

Выравнивание всех последовательностей, содержащих домен (1499) было получено в Jalview путём импорта из PFAM (Fetch sequences -> PFAM -> вводится PFAM ID домена). Файл с проектом - выравнивание в Jalview. Раскрашивание было сделано по аннотации вторичной структуры, поскольку последовательности слишком плохо выравниваются между собой (структурное сходство есть, по последовательности - не очень) и раскраска ClustalX ничего не давала. Для двух интегринов человека - ITA2B и ITAV - с выравниванием были проассоциированы 3D структуры из банка PDB. Также был получен файл с выравниванием в формате .fasta.

Упр.4 Построение филогенетического дерева по выборке доменов из трёх архитектур

С сайта PFAM были вручную получены последовательности доменов из разных организмов (скачать). Выборка была сделана среди доменов интегрина из 3х выбранных архитектур. Все последовательности относятся к одному из двух ранее выбранных таксонов (Arthropoda и Chordata). Они были закодированы идентификаторами вида NumberTaxon_organism. 1, 2, 3 соответствуют номерам архитектур в списке из первого задания, таксоны закодированы буквами A - Arthropoda и C-Chordata. Для 1 и 3 архитектуры были найдены последовательности из обоих таксонов, для 2 архитектуры - только для Chordata. Для второй архитектуры была взята 1 последовательность домена беспозвоночного (не Arthropoda, так как такая структура не представлена у Arthropoda), на дереве она обозначена как 2N_seasqirt. Далее полученные последовательности были выровнены при помощи MUSCLE (скачать выравнивание ). Затем в MEGA было построено дерево методом Maximum Likehood. На рисунке 1 представлен проект Jalview с выровненными доменами, который также можно скачать здесь. Дерево в Newick-формате можно скачать по ссылке.

Рис.1. Общий вид проекта Jalview с раскраской ClustalX по консервативности (порог - 10%)

Рис.2. Дерево, построенное по полученному выравниванию.

Упр.4 Анализ дерева

На дереве чётко выделяется клада, включающая последовательности доменов из одинаковых архитектур. Эта клада показана фиолетовым. Для остальных двух архитектур не наблюдается группировки по архитектурам, но есть выделение общепризнанных таксонов (Invertebrata, Vertebrata). Также можно выделить группу низших и высших эукариот, хотя в этих кладах на дереве есть очевидно не подходящие организмы.

Для дальнейшего анализа возьмём группу доменов из второй архитектуры. Файл с выравниваниями выбранных последовательностей можно скачать здесь. Теперь с помощью пакета HMMER на kodomo построим скрытую марковскую модель с помощью программ:

hmm2build hmm.out 2c_alignment.fasta
hmm2calibrate hmm.out
hmm2search hmm.out all_seq.fasta > hmm2search

В файле hmm2search находятся результаты поиска по всем последовательностям, содержащим этот домен, которые входят в выбранное семейство. С помощью скрипта Дмитрия Пензара построим ROC-кривую. Его скрипту на вход подаём результат работы hmm2search, файл с информацией о доменной архитектуре и таксономической принадлежности, файл, с критериями отношения к нашему семейству. Заведомо положительными мы считаем результаты с выбранной архитекрутой (2с) и таксономией. Данная архитектура включает в себя 4 домена (PF08441, PF00092, PF01839, PF00357). ROC-кривая показана на рисунках 3 и 4.

Рис.3. ROC-кривая 1

Рис.4. ROC-кривая 2