Практикум 9. Выравнивание последовательностей

-

Глобальное выравнивание:
Protein NameID_1ID_2ScoreIdentitySimilarityGapsInDels
GTPase DerDER_ECOLIDER_BACSU824.036.1%56.2%606
N-acetylglucosamine-6-phosphate deacetylaseNAGA_ECOLINAGA_BACSU404.529.9%46.5%269
UDP-N-acetylmuramoyl-L-alanyl-D-glutamate--2,6-diaminopimelate ligaseMURE_ECOLIMURE_BACSU809.036.9%53.5%5111
Локальное выравнивание:
Protein NameID_1ID_2ScoreIdentitySimilarityGapsInDelsCoverage_1Coverage_2
GTPase DerDER_ECOLIDER_BACSU831.538.7%60.2%30392.45%97.9%
N-acetylglucosamine-6-phosphate deacetylaseNAGA_ECOLINAGA_BACSU329.025.8%46.2%771396.6%94.44%
UDP-N-acetylmuramoyl-L-alanyl-D-glutamate--2,6-diaminopimelate ligaseMURE_ECOLIMURE_BACSU817.038.3%55.5%37797%95.3%
Выравнивание негомологичных белков:
Alignment typeID_1ID_2ScoreIdentitySimilarityGapsInDelsCoverage_1Coverage_2
LocalDER_ECOLIMURE_BACSU75.520.5%32.7%1342056.5%75.1%
GlobalDER_ECOLIMURE_BACSU63.025.0%46.2%34229--

Количество гэпов составляет около половины длины выравнивания, а общий счет не доходит и до ста. Довольно большой процент сходства в случае глобального выравнивания, но количество гэпов всё равно ужасно. Совпадение 20 и 25% близко к случайному.

Множественное выравнивание

Для UDP-N-acetylmuramoyl-L-alanyl-D-glutamate--2,6-diaminopimelate ligase (УДФ-N-ацетилмурамоил-L-аланил-D-глутамат-2,6-диаминопимелат-лигазы) найдено 592 записи в Swiss-prot. Я выбрала MURE_STAAB, MURE_STRPQ, MURE_MAIZE, MURE_THEMA, MURE_ARATH - просто потому, что были сверху в выдаче. Выравнивание сделано через Jalview. Видны кластеры 245-255, 266-698 (за исключением 300-305, 330-334, 404-408; делеция 446, 470-471, 610-611, 670), 720-786 (делеция 738). Длинная линкерная последовательность на N-конце. Эволюционная консервативность лежит в основе работы Alphafold и заметно отражается на укладке: на месте хорошо выровненных участков выдается высокая достоверность. Таким образом, места высокого сходства этих пяти последовательностей можно экстраполировать на участки сходства все