| ||||||||||||||||||||||||||
  |   |
Занятие 9. Мембранные белки, транспортные белки.1. Построение выравнивания заданного белка и белка-прототипа с разметкой трансмембранных сегментов1. Сравнение нумерации остатков белка-прототипа FIEF_ECOLI в UniProt и PDBИз БД PDBsum получим выравнивание последовательности из PDB с идентификатором 3H90 и соответствующей последовательности UniProt. Для этого нужно щелкнуть изображение выравнивания последовательности в структуре, скопировать полученное выравнивание, сохранить его в fasta-формате и импортировать в GeneDoc:![]() Нумерцаия аминоксилотных остатков в UniProt и PDB совпадает, но в PDB представлен лишь фрагмент последовательности из UniProt с 8 по 290 а.о., т.е. обрезанный как с С, так и с N-конца. Файл с выравниванием в формате msf: algn1.msf. 2. Построение полного глобального выравнивания заданного белка FIEF_ERWCT и белка-прототипа FIEF_ECOLIС помощью seqret получим последовательности белков FIEF_ERWCT и FIEF_ECOLI в fasta-формате: fief_erwct.fasta, fief_ecoli.fasta. Затем выровняем их при помощи программы needle пакета EMBOSS, при этом штрафы за гэпы оставим стандартные (10.0 за открытие гэпа, 5.0 за его продление). Получили выравнивание algn2.needle. Импортируем его в GeneDoc для выполнения следующих заданий: algn2.msf.Выравнивание характеризуется малым числом гэпов (2/301), процент идентичности составляет 69,8% (210/301). 3. Создание разметки трансмембранных сегментов в белке-прототипе FIEF_ECOLI по данным БД ОРМВ БД ОРМ по PDB ID найдем описание ТМ-сегментов в белке-прототипе (PDB ID 3h90). Этот белок представляет собой насос, который осуществляет выведение двухвалентного железа из клетки (ferrous-iron efflux pump fieF). Описание этого белка на сайте БД ОРМ таково:Тип: 1. Трансмембранные белки (transmembrane) Класс: 1.1. Альфа-спиральные трансмембранные белки (alpha-helical transmembrane) Надсемейство: 1.1.43. Белки, осуществляющие перенос катионов (cation diffusion facilitator) Семейство: 1.1.43.01. Бактериальные цинковые переносчики (bacterial zinc transporters) На изображении, полученном в Jmol, представлены две субъединицы белка (идентичные цепи): ![]() К полученному в упр.1.2 выравниванию algn2.msf добавим строчку с разметкой трансмембранных сегментов и сохраним в формате *.msf: mark.msf. В этом файле позиции мембранных сегментов отмечены буквой "H", цитоплазматических петель - "+", остальные (внеклеточные сегменты, а также С- и N-концы, отсутствующие в последовательности PDB) отмечены знаком "-". Из выравнивания и изображения Jmol видно, что белок имеет довольно объемную (С-концевую) цитоплазматическую часть, N-конец также находится в цитоплазме. Мембранные сегменты представляют собой 6 альфа-спиралей, расположенных практически перпендикулярно поверхности мембраны. А вот надмембранных участков мало и они небольшого размера. 4. Предсказание топологии заданного белка FIEF_ERWCT с помощью TMHMMС помощью сервера TMHMM получено предсказание топологии белка FIEF_ERWCT. Результаты этого предсказания также прикреплены к файлу mark.msf в виде искусственной последовательности "ТМНММ" с разметкой трансмембранных спиралей в соответствии с предсказанием.2. Сравнение полученного предсказания с данными ОРМДля сравнения результатов предсказания с данными ОРМ был написан скрипт Perl: count.pl. На вход программе подавались файлы с данными, полученными из OPM, и результатами предсказания TМНММ: opm.txt, tmhmm.txt.Поскольку последовательности FIEF_ERWCT и FIEF_ECOLI выровнялись без гэпов практически до конца (все мембранные участки, как предсказанные, так и по данным ОРМ, входят в этот промежуток), можно не менять координаты последовательностей для подсчета. По результатам работы скрипта составим таблицу: Результаты предсказания топологии мембранного белка FIEF_ERWCT
Чувствительность предсказания в данном случае оказалась довольно высокой (90%), то есть были предсказаны практически все мембранные остатки из существующих. Специфичность оказалась чуть ниже, но также на достаточном уровне (около 80%) - большая часть остатков, не погруженных в мембрану, тоже определилась правильно. Точность предсказания была еще ниже - среди остатков, предсказанных как мембранные, правильно угаданы были только ~70%. Недопредсказание в данном случае оказалось совсем незначительным (5,8%), а сверхпредсказание довольно значительным (порядка 30%) - за счет этого не очень велика точность предсказания. ТМНММ предсказала все 6 трансмембранных участков, причем 5 из них (первые четыре и последний по порядку их расположения в последовательности) предсказаны более-менее верно. ТМНММ в основном указывала более широкие границы для этих участков, но все действительно находящиеся в мембране остатки захватывала. Пятый участок (149-158) единственный оказался предсказанным неверно (158-175). Назад |