В базе данных Swiss-prot с помощью SRS нашла 1723 белка человека с тем же классом, 959 - с тем же классом и подклассом, 95 - с тремя одинаковыми уровнями классификациями, 52 - полностью классификация совпадает.
С помощью UniProt получила список всех белков мыши (Mus musculus) из Swiss-Prot, у которых класс совпадает с таковым моего фермента. Полученный результат сохранила в виде двух файлов, в одном из них - последовательности в Fasta-формате, в другом — текстовая таблица, в которой против ID каждого белка указан его полный EC-код.
Всего было найдено 1398 белков, у 705 из них совпадает подкласс по EC (второе число), у 95 — также и третье, у 52 — вся классификация.
Пользуясь программой blastp, нашла среди белков, найденных в предыдущем пункте, сходные по последовательности с моим ферментом. Порог по e-value поставила равный 1. Всего нашлось 100 белков. Результаты выдачи программы представила в таблице 1.
Таблица 1. Зависимость ЕС-кода от номера находки.
| Номера находки в выдаче BLASTa | Название белка | EC-код |
| 1 | Q3UFB7 | 2.7.10.1 |
| 2 | Q6VNS1 | 2.7.10.1 |
| 3 | P15209 | 2.7.10.1 |
| 4 | Q61006 | 2.7.10.1 |
| 5 | Q9Z139 | 2.7.10.1 |
| 6 | Q62371 | 2.7.10.1 |
| 7 | Q60751 | 2.7.10.1 |
| 8 | Q9WTL4 | 2.7.10.1 |
| 9 | Q03146 | 2.7.10.1 |
| 10 | Q9Z138 | 2.7.10.1 |
| 11 | P15208 | 2.7.10.1 |
| 12 | P97793 | 2.7.10.1 |
| 13 | P08923 | 2.7.10.1 |
| 14 | P16092 | 2.7.10.1 |
| 15 | P35546 | 2.7.10.1 |
| 16 | Q03145 | 2.7.10.1 |
| 17 | Q03142 | 2.7.10.1 |
| 18 | P21803 | 2.7.10.1 |
| 19 | Q61851 | 2.7.10.1 |
| 20 | Q78DX7 | 2.7.10.1 |
| 21 | Q60629 | 2.7.10.1 |
| 22 | Q4JIM5 | 2.7.10.2 |
| 23 | Q8CBF3 | 2.7.10.1 |
| 24 | P29319 | 2.7.10.1 |
| 25 | P00520 | 2.7.10.2 |
| 26 | O09127 | 2.7.10.1 |
| 27 | Q922K9 | 2.7.10.2 |
| 28 | Q61772 | 2.7.10.1 |
| 29 | Q60750 | 2.7.10.1 |
| 30 | P54763 | 2.7.10.1 |
| 31 | Q03137 | 2.7.10.1 |
| 32 | P97504 | 2.7.10.2 |
| 33 | Q03526 | 2.7.10.2 |
| 34 | P42682 | 2.7.10.2 |
| 35 | P70451 | 2.7.10.2 |
| 36 | Q62270 | 2.7.10.2 |
| 37 | P16056 | 2.7.10.1 |
| 38 | P24604 | 2.7.10.2 |
| 39 | Q62413 | 2.7.10.1 |
| 40 | Q9QVP9 | 2.7.10.2 |
| 41 | P54754 | 2.7.10.1 |
| 42 | P34152 | 2.7.10.2 |
| 43 | Q60805 | 2.7.10.1 |
| 44 | P16879 | 2.7.10.2 |
| 45 | P35991 | 2.7.10.2 |
| 46 | Q00993 | 2.7.10.1 |
| 47 | Q06806 | 2.7.10.1 |
| 48 | Q8BYG9 | 2.7.10.1 |
| 49 | P09581 | 2.7.10.1 |
| 50 | P54761 | 2.7.10.1 |
| 51 | P25911 | 2.7.10.2 |
| 52 | P55144 | 2.7.10.1 |
| 53 | P35969 | 2.7.10.1 |
| 54 | Q02858 | 2.7.10.1 |
| 55 | Q62120 | 2.7.10.2 |
| 56 | P43404 | 2.7.10.2 |
| 57 | P05480 | 2.7.10.2 |
| 58 | P39688 | 2.7.10.2 |
| 59 | P08103 | 2.7.10.2 |
| 60 | P16277 | 2.7.10.2 |
| 61 | Q01887 | 2.7.10.1 |
| 62 | P35917 | 2.7.10.1 |
| 63 | Q62190 | 2.7.10.1 |
| 64 | P06240 | 2.7.10.2 |
| 65 | Q64434 | 2.7.10.2 |
| 66 | Q9R117 | 2.7.10.2 |
| 67 | P41241 | 2.7.10.2 |
| 68 | O54967 | 2.7.10.2 |
| 69 | P14234 | 2.7.10.2 |
| 70 | Q04736 | 2.7.10.2 |
| 71 | Q5XJV6 | 2.7.11.1 |
| 72 | Q61527 | 2.7.10.1 |
| 73 | P41242 | 2.7.10.2 |
| 74 | P35918 | 2.7.10.1 |
| 75 | P52332 | 2.7.10.2 |
| 76 | Q80YE4 | 2.7.11.1 |
| 77 | P70424 | 2.7.10.1 |
| 78 | Q01279 | 2.7.10.1 |
| 79 | Q3TYD6 | 2.7.11.1 |
| 80 | Q62137 | 2.7.10.2 |
| 81 | P26618 | 2.7.10.1 |
| 82 | Q61526 | 2.7.10.1 |
| 83 | P48025 | 2.7.10.2 |
| 84 | P05532 | 2.7.10.1 |
| 85 | P05622 | 2.7.10.1 |
| 86 | Q00342 | 2.7.10.1 |
| 87 | Q99ML2 | 2.7.10.2 |
| 88 | Q3U1V8 | 2.7.11.25 |
| 89 | Q1HKZ5 | 2.7.11.25 |
| 90 | Q6J9G1 | 2.7.10.2 |
| 91 | Q60700 | 2.7.11.25 |
| 92 | Q80XI6 | 2.7.11.25 |
| 93 | Q66L42 | 2.7.11.25 |
| 94 | Q8VDG6 | 2.7.11.25 |
| 95 | Q9ESL4 | 2.7.11.25 |
| 96 | Q99N57 | 2.7.11.1 |
| 97 | P53668 | 2.7.11.1 |
| 98 | Q60855 | 2.7.11.1 |
| 99 | P18654 | 2.7.11.1 |
| 100 | P28028 | 2.7.11.1 |
Из данных, представленных в таблице, видно, что полностью функция с моим ферментом сохранилась у первой 21 находки, следовательно, их мы и будем считать достоверными гомологами. Далее вниз по таблице следуют 50 находок, у которых совпадают три уровня классификации, а затем тем белки, у которых совпадает два уровня классификации. Что интересно, ферменты, имеющие одинаковый ЕС-код с моим ферментом, встречаются в разных местах выборки, но их нельзя считать достоверными гомологами, когда уровень покрытия их меньше 70%.