ПОИСК ПРОТЕОМА, СООТВЕТСТВУЮЩЕГО ГЕНОМНОЙ СБОРКЕ
Для работы был рассмотрен организм Natranaerobius thermophilus JW/NM-WN-LF. Поиск геномной сборки проводился в базе NCBI Datasets Genome. Для данного штамма актуальной сборке соответствует RefSeq-идентификатор GCF_000020005.1, а связанный с ней идентификатор INSDC/GenBank — GCA_000020005.1. По имеющимся данным это сборка ASM2000v1. В отчёте использована ссылка именно на страницу сборки в NCBI Datasets Genome.
Далее по идентификатору сборки INSDC был выполнен поиск в UniProt Proteomes. Поисковый запрос:
genome_assembly:GCA_000020005.1
Этот запрос приводит к протеому UP000001683, связанному со сборкой GCA_000020005.1 from ENA/EMBL. На странице протеома указано, что его статус — Reference proteome. Следовательно, данный протеом не является избыточным, не был исключён из базы Proteomes и не требует указания протеома, в пользу которого он был бы удалён. В протеом входит 2847 белковых записей UniProtKB.
ПОИСК И СКАЧИВАНИЕ РЕФЕРЕНСНОГО ПРОТЕОМА
Для поиска ближайшего референсного протеома сначала был определён таксон организма в базе UniProt Taxonomy. Для вида Natranaerobius thermophilus используется TaxID 375929, а для конкретного штамма JW/NM-WN-LF фигурирует штаммовый TaxID 457570. Поиск референсного протеома был начат с уровня вида.
Использованный поисковый запрос в базе Proteomes:
(taxonomy_id:375929) AND (proteome_type:1)
Переходить к более высокому таксону не потребовалось, поскольку референсный протеом уже найден для самого исследуемого организма — это UP000001683.
Для скачивания белковых записей, принадлежащих выбранному протеому, в базе UniProtKB использовался запрос proteome:UP000001683. Для выгрузки записей в формате flat file UniProtKB (txt, соответствующем формату swiss в EMBOSS) в gzip-сжатом виде была использована следующая команда:
curl 'https://rest.uniprot.org/uniprotkb/stream?compressed=true&format=txt&query=%28proteome%3AUP000001683%29' > ~/term2/pr8/UP000001683.swiss.gzВ результате был получен файл UP000001683.swiss.gz.
ОЦЕНКА ЧИСЛА БЕЛКОВ, СОДЕРЖАЩИХ АЛЬФА-СПИРАЛИ
Для оценки числа белков, содержащих альфа-спирали, были проанализированы аннотации поля FT в записях UniProtKB. В описании формата записей UniProtKB ключ HELIX соответствует участкам альфа-спирали, а ключ TRANSMEM — трансмембранным сегментам. Для подсчёта числа записей, а не числа строк, был использован Python-скрипт, читающий gzip-файл напрямую и проверяющий наличие нужных ключей внутри каждой записи до разделителя //:
import gzip
filename = "UP000001683.swiss.gz"
records = 0
with_helix = 0
with_transmem = 0
has_helix = False
has_transmem = False
with gzip.open(filename, "rt") as fh:
for line in fh:
if line.startswith("FT HELIX"):
has_helix = True
elif line.startswith("FT TRANSMEM"):
has_transmem = True
elif line.startswith("//"):
records += 1
if has_helix:
with_helix += 1
if has_transmem:
with_transmem += 1
has_helix = False
has_transmem = False
print("total_records =", records)
print("records_with_HELIX =", with_helix)
print("records_with_TRANSMEM =", with_transmem)При запуске скрипта были получены следующие результаты: total_records = 2847, records_with_HELIX = 0, records_with_TRANSMEM = 699.
Полученный результат показывает, что в данном наборе записей отсутствуют явные аннотации с ключом HELIX, но присутствуют 699 записей с ключом TRANSMEM. Это не означает отсутствие альфа-спиралей в белках протеома. Ключ HELIX относится к явной аннотации вторичной структуры и встречается далеко не во всех записях UniProtKB. Ключ TRANSMEM описывает трансмембранные участки и относится к другой категории локальных особенностей.
Трансмембранные сегменты часто представлены альфа-спиралями, однако сам по себе ключ TRANSMEM не равнозначен ключу HELIX. В данном файле для каждого участка TRANSMEM в примечании указан тип трансмембранного сегмента, например спираль или тяж. Поэтому оценку числа белков с трансмембранными альфа-спиралями следует делать не только по наличию ключа TRANSMEM, а с учётом типа участка в примечании к этому признаку. Следовательно, оценка по HELIX в данном случае является заниженной из-за отсутствия явных аннотаций вторичной структуры, а оценка по TRANSMEM описывает отдельный класс локальных особенностей — трансмембранные сегменты.
ОЦЕНКА КОЛИЧЕСТВА ФЕРМЕНТОВ В ПРОТЕОМЕ
Для оценки количества ферментов в референсном протеоме были использованы два поисковых запроса UniProtKB по разным полям. Первый запрос был основан на наличии EC-номера:
proteome:UP000001683 AND ec:*
Этот запрос дал 619 результатов.
Второй запрос был основан на наличии блока комментария CATALYTIC ACTIVITY:
proteome:UP000001683 AND cc_catalytic_activity:*
Этот запрос дал 482 результата.
Обе оценки связаны с ферментативной активностью, но отражают её по-разному. Запрос по полю ec находит белки, для которых указан код EC, то есть ферментативная активность описана в форме классификации по катализируемой реакции. Запрос по полю cc_catalytic_activity находит записи, где в комментариях явно присутствует блок о каталитической активности.
В данном случае количество находок по запросу ec:* оказалось больше, чем по запросу cc_catalytic_activity:*. Это не является логическим противоречием, поскольку аннотации разных полей в UniProt заполняются не полностью синхронно: часть белков может уже иметь EC-номер, но не иметь отдельного комментария CATALYTIC ACTIVITY, либо такой комментарий может быть оформлен не во всех записях. Следовательно, обе оценки являются приближёнными и скорее задают нижнюю границу числа ферментов.
Запрос по ec в данном случае выглядит более чувствительным, но тоже не гарантирует нахождение всех ферментов, так как не всякая ферментативная активность в базе доведена до присвоения EC-номера. Запрос по cc_catalytic_activity является более строгим и даёт более консервативную оценку. Поэтому реальное число белков с ферментативной активностью в протеоме, вероятно, существенно выше этих значений и может быть заметно больше 619.
ВЫВОДЫ
Для Natranaerobius thermophilus JW/NM-WN-LF был найден соответствующий протеом UniProt UP000001683, связанный со сборкой GCA_000020005.1 и имеющий статус Reference proteome. Именно этот протеом был использован как референсный для дальнейшей работы.
Анализ аннотаций показал, что в протеоме имеется 699 белков с аннотированными трансмембранными сегментами, тогда как явные аннотации HELIX в скачанном файле отсутствуют. Поэтому по ключу HELIX получить осмысленную оценку числа белков с альфа-спиралями не удалось, а данные по TRANSMEM следует интерпретировать отдельно, с учётом типа каждого трансмембранного участка.
Оценка числа ферментов по запросам UniProtKB дала значения 619 и 482. Эти величины следует рассматривать как неполные и, вероятно, заниженные оценки, поскольку реальное число белков с ферментативной активностью в протеоме может быть существенно выше.