Практикум 9
Глобальное парное выравнивание гомологичных белков
Табл.1 Характеристики глобального парного выравнивания трёх пар белков
| Protein Name |
ID 1 |
ID 2 |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
| Adenine phosphoribosyltransferase |
APT_ECOLI |
APT_BACSU |
441.5 |
50.3% |
61.2% |
13 |
3 |
| Flagellar motor switch protein |
FLIM_ECOLI |
FLIM_BACSU |
369.0 |
28.8% |
50.9% |
14 |
8 |
| DNA polymerase III subunit delta |
HOLA_ECOLI |
HOLA_BACSU |
73.0 |
18.2% |
33.8% |
110 |
19 |
Локальное парное выравнивание гомологичных белков
Табл.2 Характеристики локального парного выравнивания трёх пар белков
| Protein Name |
ID 1 |
ID 2 |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
Coverage 1 |
Coverage 2 |
| Adenine phosphoribosyltransferase |
APT_ECOLI |
APT_BACSU |
450.0 |
56.2% |
66.7% |
0 |
0 |
88.5% |
95.3% |
| Flagellar motor switch protein |
FLIM_ECOLI |
FLIM_BACSU |
374.5 |
29.4% |
52.0% |
9 |
6 |
94.9% |
96.4% |
| DNA polymerase III subunit delta |
HOLA_ECOLI |
HOLA_BACSU |
92.0 |
19.5% |
37.5% |
66 |
14 |
87.8% |
89.6% |
Значения Score, Identity и Similarity практически совпадают, однако при выравнивании water значения несколько больше. Значения в Gaps и Indels различаются довольно сильно, needle для первого выравнивания выдал 13 гэпов, а water - 0, для второго 14 и 9, для третьего 110 и 66 соответственно. Таким образом, water выравнивает последовательности так, чтобы уменьшить количество гэпов, в результате повышается значения Score, Identity и Similarity. Если были выравнены меньшие по длине участки последовательности или начало выравнивания изменилось, то пары букв, сопоставленных в локальном выравнивании, оказывались не сопоставлены в глобальном.
Adenine phosphoribosyltransferase
При глобальном выравнивании Identity: 50.3%, при локальном - 56.2% это говорит о том, что белки не гомологичны по всей длине, но довольно высокие значения Identity и Similarity подтверждают наличие гомологичных участков. Локальное выравнивание информативно, оно подтверждает гомологичность, но также указывает на небольшие несовпадения в отдельных областях.
Flagellar motor switch protein
Значения Identity и Similarity у белков низкие. Локальное выравнивание показывает высокое покрытие (94.9% и 96.4%), что может указывать на наличие гомологичных участков, но белки не являются гомологичными по всей длине (вероятно, есть консервативные домены, разделённые негомологичными участками). В данном случае локальное выравнивание так же более информативно, так как оно показывает, что гомологичные участки занимают почти всю длину.
DNA polymerase III subunit delta
Как в глобальном, так и в локальном выравниваниях низкие значения Identity и Similarity, поэтому белки не являются гомологичными по всей длине. Довольно высокие значения покрытия подтверждают наличие гомологичных участков. Локальное выравнивание выявляет гомологичные участки, но общее сходство остаётся низким.
Результат применения программ выравнивания к неродственным белкам
Глобальное выравнивание
Табл.3 Характеристики глобального выравнивания пары неродственных белков
| Protein Name |
ID 1 |
ID 2 |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
Carbamoyl phosphate synthase small chain /
Magnesium transport protein CorA |
CARA_ECOLI |
CORA_BACSU |
32.0 |
9.7% |
17.4% |
349 |
17 |
Локальное выравнивание
Табл.4 Характеристики локального выравнивания пары неродственных белков
| Protein Name |
ID 1 |
ID 2 |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
Coverage 1 |
Coverage 2 |
Carbamoyl phosphate synthase small chain /
Magnesium transport protein CorA |
CARA_ECOLI |
CORA_BACSU |
54.5 |
19.4% |
32.6% |
98 |
11 |
47.4% |
55.2% |
Количество гэпов в локальном выравнивании снова оказалось значительно меньше, чем в глобальном, что повлияло на Score. Значения Identity и Similarity низкие, то есть белки не являются гомологичными по всей длине. Coverage меньше 60% это означает, что большая часть последовательностей не совпадает. Глобальное выравнивание подтверждает, что белки не родственны по общей структуре (низкие Score и %Identity, огромное число гэпов). Локальное выравнивание показывает некоторое слабое сходство на отдельных участках, что может быть результатом случайного совпадения. Все эти данные указывают на отсутствие эволюционного родства белков.
Множественное выравнивание белков и импорт в Jalview
Мнемоника: HOLA, 11 записей в Swiss-Prot
Полное имя: DNA polymerase III subunit delta
Выравнивание было сделано с помощью Jalview
Последовательности выровнены хорошо, с равномерными промежутками и правильно сформированными столбцами. Хотя белки и имеют общие функции нельзя сказать, что они гомологичны по всей длине. Есть несколько гомологичных участков в области столбцов: 30, 80, 190, 210, 273, 351, 362. Однако также есть области со значительным расхождением, что указывает на вариативность этих участков у гомологов. В выравнивании встречается относительно невысокое количество гэпов.
График консервативности (под выравниванием) показывает уровень консервативности каждого остатка в последовательностях. Сильно консервативные позиции отмечены высокими столбцами (окрашенны в коричневый и жёлтый), которые указывают на функциональную или структурную важность этих остатков. На графике показаны участки как с высокой, так и с низкой консервативностью, что говорит о наличии консервативных доменов.