Практикум 7. Трансмембранные белки

1. Выбранный белок

Для работы был выбран белок флиппаза липида II MurJ (Lipid II flippase MurJ), структура которого имеет идентификатор PDB 6NC9. Идентификатор UniProt - B7IE18 (MURJ_THEAB). Белок состоит из 475 аминокислот и происходит из бактерии Thermosipho africanus strain TCF52B.

MurJ находится во внутренней (цитоплазматической) мембране бактерии и является многопроходным α-спиральным мембранным белком. Он участвует в синтезе пептидогликана: переносит липид-связанный предшественник клеточной стенки с внутренней стороны мембраны на внешнюю. Структура 6NC9 соответствует outward-facing конформации белка.

2. Трансмембранные участки по данным OPM

По данным OPM белок содержит 14 трансмембранных α-спиралей. Для цепи A указаны следующие координаты:

  1. 3–28
  2. 41–65
  3. 80–102
  4. 124–146
  5. 150–171
  6. 178–200
  7. 211–238
  8. 253–273
  9. 290–313
  10. 331–355
  11. 359–381
  12. 388–409
  13. 414–434
  14. 440–462

3. Предсказание DeepTMHMM

Для анализа использовалась аминокислотная последовательность MurJ, взятая из записи UniProt B7IE18. Последовательность была запущена в DeepTMHMM. Программа определила белок как alpha TM и также предсказала 14 трансмембранных участков.

Координаты, предсказанные DeepTMHMM: 8–29, 38-60, 80–100, 121–143, 152–171, 178–200, 218–241, 246–269, 291–314, 329–352, 363–385, 387–407, 414–434, 441–462.

Текстовая выдача DeepTMHMM в формате GFF3.

Результат DeepTMHMM для белка MurJ 6NC9
Рис. 1. Результат предсказания трансмембранной топологии MurJ с помощью DeepTMHMM. На верхнем графике показана наиболее вероятная топология: красным обозначены трансмембранные участки, розовым — участки с внутренней стороны мембраны, синим - с внешней. На нижнем графике показаны вероятности принадлежности каждого остатка к этим состояниям.

4. Сравнение OPM и DeepTMHMM

Таблица 1. Сравнение координат трансмембранных участков MurJ.
Участок OPM DeepTMHMM Сравнение
TM13–288–29перекрываются
TM241–6538–60перекрываются
TM380–10280–100перекрываются
TM4124–146121–143перекрываются
TM5150–171152–171перекрываются
TM6178–200178–200полное совпадение
TM7211–238218–241перекрываются
TM8253–273246–269перекрываются
TM9290–313291–314перекрываются
TM10331–355329–352перекрываются
TM11359–381363–385перекрываются
TM12388–409387–407перекрываются
TM13414–434414–434полное совпадение
TM14440–462441–462перекрываются

Как видно из таблицы 1, результаты двух методов совпадают очень хорошо. OPM и DeepTMHMM нашли одинаковое число трансмембранных участков - 14, и каждый участок одного метода перекрывается с соответствующим участком другого. TM6 (178–200) и TM13 (414–434) совпали полностью. Для большинства остальных спиралей границы отличаются только на несколько аминокислот; максимальное различие границы составляет 7 аминокислот для TM7 и TM8.

Небольшие различия ожидаемы из-за разных подходов. DeepTMHMM делает предсказание только по аминокислотной последовательности, а OPM использует трехмерную структуру белка и его положение относительно мембраны. Поэтому точные границы одной и той же спирали могут немного отличаться.

Вывод: DeepTMHMM очень хорошо воспроизвел трансмембранную топологию MurJ: ни одна из 14 α-спиралей не была пропущена, а координаты всех предсказанных участков хорошо перекрываются с данными OPM.