Практикум №10. HMM-профили и эволюционные домены.

Для построения HMM-профиля было выбрано семейство Aminotransferase class I/II (Pfam AC: PF00155). Белки этого семейства являются пиридоксаль-5'-фосфат (PLP)-зависимыми ферментами и катализируют широкий спектр реакций трансаминирования, участвуя в метаболизме аминокислот, биосинтезе вторичных метаболитов и других процессах. Последовательностей в seed и full — 48 и 171 269 соответственно.

Далее было выбрано подсемейство с доменной архитектурой PF03744 - PF00155. Всего было найдено 32 таких белка.

Белковые последовательности этого подсемейства были выравнены алгоритмом MAFFT для создания HMM-профиля.

Выравнивание последовательностей подсемейства

HMM профили

Поиск профилем по белкам семейства

Последовательности белков семейства

Так как белков семейства слишком много (около 310 тысяч), были взяты только reviewed записи, таких оказалось около 1 тысячи.

Таблица результатов

Таблица 1. Численные характеристики выделения подсемейства профилем.
Принадлежность белка подсемейству Классификация по HMM-профилю
Score ≥ порога Score < порога
Белок подсемейства TP = FN =
Белок не из подсемейства FP = TN =