Классификация ферметов (EC)
Сбор информации о ферменте O95336
Для фермента O95336 (6-phosphogluconolactonase) определим EC-код. Расшифруем EC-код и переведем на русский активность, соответствующую каждому из четырех чисел кода, воспользовавшись сайтом "Enzyme Nomenclature" (NC-IUBMB).
EC код для 6-фосфоглюконолактоназы: EC=3.1.1.31
Расшифруем: что означают отдельные цифры в коде:
EC 3. - Hydrolases (Гидролазы - это класс ферментов, катализирующий гидролиз ковалентной связи). Общий вид реакции,
катализируемой гидролазой выглядит следующим образом:
AB + H2O > AOH + BH
EC 3.1 - Acting on Ester Bonds (Действует на сложноэфирные связи)
EC 3.1.1 - Carboxylic Ester Hydrolases (Карбоксэфирные гидролазы)
EC 3.1.1.31 - 6-phosphogluconolactonase (6-фосфоглюконолактоназа)
Сколько ферментов человека с близкими функциями описано в Swiss-Prot?
С помощью SRS выясним, сколько белков человека являются ферментами с тем же классом по EC, сколько — с тем же классом и подклассом, сколько имеют общие с ферментом O95336 три уровня классификации и сколько — все четыре. Эти данные представлены в табл.1.
Таблица 1. Количество ферментов с близкими функциями описанные в Swiss-Prot
EC 3. |
EC 3.1 |
EC 3.1.1 |
EC 3.1.1.31 |
|
Общий уровень классификации | Класс | Подкласс | Подподкласс | Порядковый номер |
Количество белков | 168297 | 15293 | 2644 | 127 |
Насколько сохраняется функция фермента у белков со сходными последовательностями?
С помощью UniProt получим список всех белков мыши (Mus musculus) из Swiss-Prot, у которых класс (т.е., первое число) классификации EC совпадает с таковыми фермента O95336. Для этого, в поле "Query" напишем
Полученные данные и дальнейшая с ними работа приведена на рис.2.
Для того, что бы получить только нужную информацию из всего многообразия полученнных данных, зайдем в раздел "Customize". |
Для того, что бы получить данные для ID белка и EC оставим в правой колонке только запись EC number. Добавим ее кнопкой "show" из левого столбца, а все остальные строки из правого столбца скроем кнопкой "Hide". |
Полученный результат сохраним в виде двух файлов, в одном последовательности в Fasta-формате, в другом — текстовая таблица, в которой против ID каждого белка указан его полный EC-код.
Всего нашлось 1,719 белков.
Таблица 2. Количество белков M.musculus в Swiss-Prot c конкретным EC
EC 3. |
EC 3.1 |
EC 3.1.1 |
EC 3.1.1.31 |
|
Общий уровень классификации | Класс | Подкласс | Подподкласс | Порядковый номер |
Количество | 1719 | 624 | 94 | 2 |
Теперь с помощью программы blastp, среди найденных белков, определим сходные по последовательности с ферментом O95336. При этом в качестве порога по e-value будем использовать значение равное 1.
-outfmt 7 ( 7 = tabular with comment lines)
На выходе мы получим файл - 6PGL_HUMAN_homologs.out.
Таблица 3. Результаты поиска схожих с ферментом P08319 последовательностей в Mus musculus
Позиция в выдаче |
ID белка |
Общий уровень |
E-value |
EC-код |
%idenity |
Можно считать |
1 | Q9CQ60 | Вся классификация | 4e-158 | 3.1.1.31 | 88.76 | да |
2 | Q8CFX1 | Вся классификация | 4e-22 | 1.1.1.47; 3.1.1.31 | 29.74 | нет |
3 | Q9CRC9 | Класс | 5e-06 | 3.5.99.6 | 23.39 | нет |
4 | O88958 | Класс | 2e-04 | 3.5.99.6 | 23.98 | нет |
5 | Q9EQF5 | Класс | 0.31 | 3.5.2.2 | 37.50 | нет |
6 | D3Z7P3 | Класс | 0.77 | 3.5.1.2 | 43.33 | нет |
Табл.3 отражает насколько сохраняется функция (т.е. подкласс, подподкласс и порядковый номер EC) в зависимости от позиции находки в выдаче blastp.
© Nuzhdina Ekaterina, 2012