Выравнивание последовательностей
Таблица 1. Характеристики глобального парного выравнивания трёх пар белков
Protein name |
ID 1 |
ID 2 |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
Protein-glutamate methylesterase/protein-glutamine glutaminase |
CHEB_ECOLI |
CHEB_BACSU |
642.5 |
38.9% |
58.2% |
30 |
9 |
2-dehydro-3-deoxygluconokinase |
KDGK_ECOLI |
KDGK_BACSU |
194.5 |
22% |
37.5% |
77 |
15 |
Aminodeoxychorismate lyase |
PABC_ECOLI |
PABC_BACSU |
213.5 |
24.8% |
40.6% |
34 |
6 |
Таблица 2. Характеристики локального парного выравнивания гомологичных белков
Protein name |
ID 1 |
ID 2 |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
Coverage1 |
Coverage2 |
Protein-glutamate methylesterase/protein-glutamine glutaminase |
CHEB_ECOLI |
CHEB_BACSU |
645.5 |
40.9% |
60% |
25 |
6 |
94.8% |
96.4% |
2-dehydro-3-deoxygluconokinase |
KDGK_ECOLI |
KDGK_BACSU |
195.5 |
23.8% |
39.8% |
53 |
12 |
96.8% |
91.4% |
Aminodeoxychorismate lyase |
PABC_ECOLI |
PABC_BACSU |
222 |
26% |
43.8% |
12 |
4 |
94.4% |
88.7% |
Таблица 3. Характеристики глобального парного выравнивания неродственной пары белков
Protein ID 1 andname |
Protein ID 2 and name |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
ALLE_ECOLI | (S)-ureidoglycine aminohydrolase |
EPSN_BACSU | Putative pyridoxal phosphate-dependent aminotransferase EpsN |
22.5 |
3.8% |
5.9% |
505 |
7 |
Таблица 4. Характеристики локального парного выравнивания неродственной пары белков
Protein ID 1 andname |
Protein ID 2 and name |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
Coverage1 |
Coverage2 |
ALLE_ECOLI | (S)-ureidoglycine aminohydrolase |
EPSN_BACSU | Putative pyridoxal phosphate-dependent aminotransferase EpsN |
44 |
20.2% |
32.4% |
75 |
11 |
66% |
25.5% |
Множественное выравнивание белков в Jalview
Для выравнивания я выбрала мнемонику TADA. Всего их 26. Помимо белков кишечной и сенной палочки я выбрала белки гемофильной палочки, бучнеры (эндомимбиот тли), стрептококка и резуховидки.
Выравнивание я делала в Jalview.
Ссылка на выравнивание
Параметры программ needle и water
При запуске программ needle и water без опции -auto запрашивает ввод новых значений для параметров Gap_penalty (штраф за индель) и Extend_penalty (штраф за продолжения инделя).
Для примера я решила изменить условия глобального выравния мнемоники CHEB. Результаты в таблице:
Condition |
ID 1 |
ID 2 |
Score |
%Identity |
%Similarity |
Gaps |
Indels |
-auto |
CHEB_ECOLI |
CHEB_BACSU |
642.5 |
38.9% |
58.2% |
30 |
8 |
Gap_penalty: 100; Extend_penalty: 10 |
CHEB_ECOLI |
CHEB_BACSU |
280 |
32.4% |
49% |
16 |
3 |
При сравнении видно, что рейтинг сильно понизился при повышении штрафа; число инделей уменьшилось, как и их длина.