Классификация ферментов EC
1.Сбор информации о своем ферменте
DLDH_HUMAN-фермент который будет исследоваться в этой работе.
Соберем все базовые данные о ферменте через запись в банке SwissProt/UniProt.
Full=Dihydrolipoyl dehydrogenase, mitochondrial
Перевод=Митохондриальная дегидрогеназа дигидролипоамида
EC=1.8.1.4
1.---оксиредуктазы
1.8- маленький подкласс, реагирующий с неорганическими и органическими тиолами.
1.8.1 - с НАДФ(+) и НАД(+) в качестве акцепторных ионов.
IntEnz
Функция: Принимает участие в процессах расщепления, как глицина, так и альфа-кетокислотных-дегидрогеназных сложных соединений.
Реакция:
Protein N(6)-(dihydrolipoyl)lysine + NAD(+) <=> protein N(6)-(lipoyl)lysine + NADH
2.Сколько ферментов человека с близкими функциями описано в Swiss-Prot?
Запрос |
1.8.1.4 |
1.8.1 |
1.8 |
1. |
Общее количество |
73 |
547 |
1536 |
37924 |
Совпадений у человека |
1 |
7 |
26 |
571 |
3.Насколько сохраняется функция фермента у белков со сходными последовательностями?
Запрос |
1.8.1.4 |
1.8.1 |
1.8 |
1. |
Совпадений с Mus Musculus |
1 |
7 |
25 |
549 |
После заполнения этой таблицы мы проводим самый последний поиск с сохранением результатов в виде двух файлов:
1)fasta-файл
2)Текстовая таблица
formatdb -i EC.fasta -p T -n base
blastall -p blastp -d base -i EC.fasta -e 1 -o blastp.out
sp|P24815|3BHS1_MOUSE 3 beta-hydroxysteroid dehydrogenase/Delta ... 736 0.0
sp|O35469|3BHS6_MOUSE 3 beta-hydroxysteroid dehydrogenase/Delta ... 626 0.0
sp|P26149|3BHS2_MOUSE 3 beta-hydroxysteroid dehydrogenase/Delta ... 619 0.0
sp|P26150|3BHS3_MOUSE 3 beta-hydroxysteroid dehydrogenase/Delta ... 615 0.0
sp|Q61767|3BHS4_MOUSE 3 beta-hydroxysteroid dehydrogenase type 4; 552 0.0
sp|Q61694|3BHS5_MOUSE 3 beta-hydroxysteroid dehydrogenase type 5; 538 0.0
sp|Q9EQC1|3BHS7_MOUSE 3 beta-hydroxysteroid dehydrogenase type 7; 231 2e-73
sp|Q9D665|D42E1_MOUSE Short-chain dehydrogenase/reductase family... 108 4e-28
sp|Q9R1J0|NSDHL_MOUSE Sterol-4-alpha-carboxylate 3-dehydrogenase... 101 8e-26
sp|Q8K3F7|TDH_MOUSE L-threonine 3-dehydrogenase, mitochondrial; 33 0.005
sp|Q7TNT2|FACR2_MOUSE Fatty acyl-CoA reductase 2; 30 0.064
Полный итоговый файл:blastp.out
ID |
EC |
Уровень совпадения |
3BHS1_MOUSE |
1.1.1.145 |
класс |
3BHS6_MOUSE |
1.1.1.145 |
класс |
3BHS2_MOUSE |
1.1.1.145 |
класс |
3BHS3_MOUSE |
1.1.1.145 |
класс |
3BHS4_MOUSE |
1.1.1.145 |
класс |
3BHS5_MOUSE |
1.1.1.145 |
класс |
3BHS7_MOUSE |
1.1.1.145 |
класс |
D42E1_MOUSE |
1.1.1 |
класс |
NSDHL_MOUSE |
1.1.1.170 |
класс |
TDH_MOUSE |
1.1.1.103 |
класс |
FACR2_MOUSE |
1.2.1 |
класс |
Почти полное отсутствие совпадений по EC кодам было ожидаемо так как в таблице в начале этого задания было указано чрезвычайно низкое количество совпадений даже для первой цифры EC кода. Что удивительно, при этом мы можем наблюдать достаточно неплохие веса выравниваний для выбранных програмой blastp последовательностей. В результате, мы можем считать настоящими гомологами все последовательности, кроме последних двух.
©Melnichuk Anastasia