На главную страницу сайта На главную страницу третьего семестра |
---|
Программы пакета BLAST для работы с нуклеотидными последовательностями
Поиск гомологов белка BTUB_ECOLI в геномах родственных бактерий
Сначала был проведен поиск гомологов белка BTUB_ECOLI в геноме Pasteurella multocida. Для этого были созданы индексные
файлы, соответствующие этому геному при помощи программы formatdb:
Выравнивание получено запуском программы blastall, для поиска был выбран алгоритм tblastn:
После анализа результатов при помощи программы entret была получена запись EMBL, соответствующая лучшей находке.
Создание индексных файлов осуществлялось следующим образом:
Выравнивание было получено аналогично предыдущему поиску.
Результаты поиска представлены в таблице:
Увеличение E-value AE006199 при поиске по трем геномам связано с возрастанием объема банка данных. При этом эта находка теперь находится на 36 месте, хотя из Pasteurella multocida она остается лучшей. Таким образом можно сказать, что в Pasteurella multocida нет белков достаточно близких по строению и функции к btuB. По видимому, это связано с наличием другого метода транспорта крупных молекул (в частности, содержащих корриновые и порфириновые кольца), при помощи белков-переносчиков с отличной от btuB структурой. А вот в Xanthomonas campestris и Salmonella typhimurium шанс найти схожие белки гораздо выше. Проведя поиск по базе SwissProt, действительно было обнаружено что белок соответсвующий лучшей находке (AC:AE008893) — это BTUB_SALTY из Salmonella typhimurium, который является непосредственным родственником BTUB_ECOLI и выполняет ту же функцию — транспорт кобаламина. Вторая находка из Xanthomonas campestris (AC:AE012421) соответствует белку Q8P6A3_XANCP, который является рецептором витамина В12. Поиск гомологов гена, кодирующего белок BTUB_ECOLI, с помощью программы BLASTN.По результатам e-value лучшей находки предыдущего поиска увеличился до 4e-68. Это можно объяснить увеличением объема базы данных(с аминокислотной до нуклеотидной). E-value второй и третьей находки гораздо выше (3e-04). Причем большое число хоть сколько-нибудь отдаленных гомологов вообще не были найдены, например AC:AE012421. Таким образом видно, что непосредственных родственников программа BLASTN находит с лёкостью, но когда речь идет о небольших расхождениях, ей не под силу разобрать кто из кандитатов ближе, не говоря уже об очень дальних гомологах. |