Совмещение

С помощью сервиса PDBeFold, производящего структурные выравнивания белков с моделями из PDB, был произведен поиск структурных гомологов белка Fis1 (PDB ID: 3o48), использованного в предыдущих практикумах. Сервис выдал для Fis1 127 находок, из которых 12 имели RMSD 0,8-3 Å и N_align 50-90%. Из них для дальнейшей работы были выбраны 4 гомолога (Таблица 1).

Таблица 1. Структурные гомологи 3o48, найденные PDBeFold
PDB НазваниеОрганизмN_matchN_alignN_align/N_queryRMSD
1pc2:AHuman mitochondria fission proteinHomo sapiens1521130,8761.64
1iyg:AFis1p-like and CGI-135 homologous domainMus musculus1331160,8991.82
3sz7:ASgt2 TPR domainAspergillus fumigatus1511060,8222.40
2vgx:BType III Secretion Translocator Chaperone SycDYersinia enterocolitica1451000,7752.36
С данными гомологами было построено множественное структурное выравнивание (fasta-файл). Количество выравненных остатков для данного выравнивания - 99, число выравненных элементов вторичной структуры - 6, общий RMSD - 2.559, общий Q-score - 0.2893.

Файл с совмещёнными структурами: send.pdb

На рисунке 1 показано совмещение выравненных структур. В целом, они достаточно хорошо накладываются друг на друга за исключением петель и более длинных, чем у 3о48, концевых неструктурированных участков отдельных гомологов.

Рисунок 1. Пространственное cовмещение структурных гомологов 3о48 в двух проекциях в изображении cartoons и backbone. Deepteal - 3о48, hotpink - 1iyg:A, purple - 1pc2:А, magenta - 2vgx, lightpink - 3sz7
На рисунке 2 показано структурное выравнивание гомологов, визуализированное с помощью JalView. Рамкой выделена query-структура 3о48. На рисунке 3 - выравнивание тех же гомологов, построенное программой Muscle с дефолтными параметрами.

Выравнивания отличаются достаточно сильно. Визуально, в выравнивании по последовательностям несколько меньше гэпов и невыровненных остатков.

Рисунок 2. Структурное выравнивание гомологов Fis1

Рисунок 3. Выравнивание гомологов Fis1 по последовательностям, построенное Muscle
Рисунок 4. Остаток Gln в 3о48 (бирюзовый), выравненный с
остатками Glu и Asp из 1iyg, 1pc2 и 2vgx (оттенки розового)
по структуре, но не по последовательности. Справа - Leu из
3о48, стоящий в аналогичном положении в выравнивании
последовательностей
Рисунок 5. Остаток Ser 3sz7 (светло-розовый), выравненный с
остатками Ser из 1iyg, 1pc2 и 3о48 по структуре, но не по
последовательности, и остаток Ala 3sz7, выравненный с
остатками Ser из 1iyg, 1pc2 и 3о48 по последовательности, но
не по структуре

Рассмотрим несколько различающихся позиций подробно.

Так, в позиции 40 структурного выравнивания в 3о48 стоит Gln, выравненный с тройкой отрицательно заряженных полярных аминокислот (EED) из других гомологов. В выравнивании последовательностей с этой же тройкой выравнен Leu (позиция 19). Посмотрим, как данный участок выглядит на совмещении структур (Рис. 4). Полярные остатки EED из структур 1iyg, 1pc2 и 2vgx соответственно показаны в оттенках розового, два варианта выровненного с ними остатка (Gln "по версии структурного выравнивания" и Leu "по версии выравнивания последовательностей") показаны бирюзовым. Видно, что Gln, предсказанный структурным выравниванием, довольно хорошо совмещается c остатками из других гомологов и вместе с ними лежит в участке α-спирали, в то время как предсказанный Muscle Leu находится далеко от остальных остатков и лежит в петле. Скорее всего, более правильным в данном случае является выравнивание по структуре.

В позиции 80 структурного выравнивания наблюдается 4 серина в структурах 3о48, 1iyg, 1pc2 и 3sz7. В соответствующей колонке выравнивания по последовательностям на месте Ser 3sz7 стоит аланин. На рисунке 5 представлено изображение данных остатков в совмещенных структурах. В данном случае "альтернативные" остатки являются соседями по цепи. Разумеется, Ser из 3sz7 (светло-розовый), накладывается на все другие Ser лучше, чем Ala из той же структуры, однако по расположению в пространстве Ala отстоит не особенно далеко, к тому же в выравнивании по последовательностям он совпадает с Ala из гомолога 2vgx (на рисунке не показан). Поэтому в этой ситуации сложно однозначно сказать, какое из выравниваний более правильное.










.