- Создание профиля для доменов белков группы специфичности frur
Для начала надо добавить веса в выравнивании доменов белков. Для этого воспользуемся командой:
pfw -m frur.fasta > out.w.fasta
Выходной файл: out.w.fasta
Затем дадим на вход програме pfmake полученный файл для построения профиля:
pfmake -m out.w.fasta /usr/share/pftools23/blosum45.cmp > profile.prf
Выходной файл: profile.prf
Затем дадим получившийся профиль программе autoscale для нормировки относительно случайной базы:
autoscale -m profile.prf > profile.scale.prf
Выходной файл: profile.scale.prf
- Поиск белков заданной группы специфичности (frur) в протеоме Marinomonas sp. (strain MWYL1)
Для этого скачаем файл с протеомом Marinomonas sp. (strain MWYL1).
Затем проведем в протеоме поиск белков заданной группы специфичности (frur).
Для этого воспользуемся программой pfsearch:
pfsearch -f profile.scale.prf proteom.fasta > pan.search;
Выходной файл: pan.search
pfsearch -x -f profile.scale.prf proteom > pan_s.fasta
Выходной файл: pan_s.fasta
Далее используем программу ClustalW2 (на сервере kodomo-count запускается командой: /home/preps/golovin/progs/bin/clustalw2)
для выравнивания найденных последовательностей под выравнивание семейства. На выходе получаем
выравнивание
и изображение:
Зеленым выделены ДНК-связывающие домены группы специфичности frur.
Снизу расположены последовательности 12 белковых находок.
Данные находки вряд ли принадлежат к группе специфичности frur, т.к. они, как минимум, не содержат специфических позиций, характерных для данной группы специфичности, т.е. уже на этом этапе они не проходят отбор.