Для выполнения данного задания я выбрала 1 из белков, рассматриваемых в 1 задании:
FMT_*, где * - принадлежность к какому-либо организму (полное имя белка см. в задании1). Через базу
UniProt KB и ее расширенный поиск я посмотрела, у каких организмов встречается данный белок. Учитывая банк
Swiss-Prot, запрос выдал 696 результатов, из которых я выбрала 7 следующих:
FMT_HUMAN, FMT_ECOLI, FMT_YEAST, FMT_BOVIN, FMT_SALTY, FMT_SCHPO, FMT_BACSU;
соответственно принадлежащие организмам:
-Homo sapiens (Human), Escherichia coli (strain K12), Saccharomyces cerevisiae, Bos taurus (Bovine), Salmonella typhimurium (strain LT2), Schizosaccharomyces pombe (strain 972), Bacillus subtilis (strain 168).
Затем я создала файл
FMT_txt, строками которого являются
sw:FMT_*. Далее командами
seqret @FMT.txt FMT.fastа
muscle -in FMT.fasta -out FMT_alignment.fasta
формируется нужный файл в формате fasta. Полученный файл был открыт в
Jalview, ссылка на который представлена ниже:
Исходя из результатов, сложно судить о гомологии всех 7 выбранных белков, так как визуально можно оценить большое количество гэпов, например блоки [347-358], [283-295], крупный блок [1-49], а также малую протяженность консервативных блоков. Наряду с ними можно увидеть схожие блоки: [52-72], [85-92], [121-138], но самое заметное сходство наблюдается в блоке [160-262], в котором практически нет даже единичных гэпов и заметны консервативные колонки.
Также можно рассмотреть белки попарно. В упомянутом блоке с 1 по 49 аминокислотный остаток наблюдается практически абсолютное сходство между белковыми последовательностями
FMT_HUMAN и
FMT_BOVIN, принадлежащими таксиномически близким организмам, поэтому между ними наблюдается наибольшее сходство, что говорит о гомологии данных белков. Подобную ситуацию можно также наблюдать при сравнении последовательностей
FMT_ECOLI и
FMT_SALTY имеющих одинаковые длину последовательноти,консервативные участки и даже количество гэпов.