Учебная страница курса биоинформатики,
год поступления 2018
Практикум 2. Атлас контактов
Исследовательский мини-проект, направленный на изучение контактов биомолекул в составе комплексов, а также контактов с биомолекул с лигандами. В качестве комплексов для изучения можно рассмотреть, например, фрагменты рибосомы, комплексы ВИЧ-интегразы, комплексы биомолекул с антителами.
Для облегчения поиска воспользуйтесь возможностями Advanced Search на сайте PDB. Не стоит выбирать комплексы весом больше 100кДа
Этапы работы:
1. Разбиться на группы по 2-3 человека и записаться в таблицу. Внести ссылку на страницу, на которой будет размещен отчет о вашей деятельности. (Отчет можно разместить только у одного участника команды, но остальные должны поставить ссылку на своих учебных страницах)
2. Время на выполнение: 2 недели. Хотелось бы, конечно, увидеть хотя бы по одному примеру для каждого контакта, но минимум - по 2 контакта (пусть даже одного типа, например, 2 водородных связи) на каждого члена команды.
Формат отчета:
html страница (наличие двух (Ru/En) и более языковых вариантов будет плюсом) с текстовым содержимым и графическими иллюстрациями, полученными при помощи jmol.
Страница с описанием контактов должна содержать следующую информацию:
Введение:
- Описание выбранной биомолекулы (или комплекса)(ее функциольная роль, PDB ID, Uniprot ID отдельных биомолекул и любая уместная информация).
- Со ссылками на литературу
- Описание низкомоллекулярных лигандов (особенно функциональных) в составе комплекса
- Изображение в шариково-стержневой модели с подписанными именами тяжелых атомов + jmol-апплет с соответствующей моделью
В качестве образца описания низкомоллекулярных лигандов можно рассматривать страницы с Википедии. Образец, но не копипасту!
Информация о физико-химических свойствах:
- Название по IUPAC
- Химическая формула
- Брутто-формула
- Молярная масса
Ссылку на идентификатор в базе данных PubChem
Информация о белок-белковых контактах:
Внимательно изучите структуру выбранного вами комплекса и попытайтесь самостоятельно определить ключевые точки межмоллекулярных контактов.
Интересующие нас типы контактов:
- Ковалентные связи
- Водородные связи
Чтобы в PDB файле присутствовали протоны, надо искать структуры с разрешением < 1.2A. Если не получается найти структуры протонированных форм, обсудите вероятность образования данного контакта в зависимости от pKa вовлеченных аминокислот.
- Солевые мостики
- Гидрофобные контакты
Исследуйте плотность упаковки атомов в гидрофобном ядре.
- Выберете гидрофобный остаток, боковая цепь которого упакована в гидрофобное окружениеа. Выбирайте крупный остаток, например тирозин или фенилаланин.
Определите свой остаток в команде define под именем MyResidue. Создайте в скрипте следующие изображения (с раскраской по химическим элементам):
- Изображение вашего остатка в виде ван-дер-ваальсовых радиусов зеленого цвета. Отцентрованное.
b. Предыдущее изображение + 'cpk 20' для соседних остатков (у которых хотя бы один атом находится на расстоянии не более 7 ангстрем от MyResidue). Покрасьте атомы MyResidue зеленым, а остальные атомы - по химическим элементам.
Соседние атомы можно выделить командой 'select within (7.0, MyResidue)', но учтите - сами атомы из MyResidue находятся от самих себя на расстоянии 0 ангстрем, что меньше семи. Так что они тоже выделяются этой командой. Чтобы команда 'cpk 20' не заменила ван-дер-ваальсовы радиусы атомов из MyResidue на маленькие точки, надо придумать более сложную команду select. Используйте логические операторы - AND, OR, NOT.
c. Предыдущее изображение + ван-дер-ваальсовы радиусы для атомов, находящихся на расстоянии не более 1 ангстрем от MyResidue. Цвет - по химическим элементам (cpk)
d. Предыдущее изображение + ван-дер-ваальсовы радиусы для атомов, находящихся на расстоянии не более 2 ангстрем от MyResidue.
e. Предыдущее изображение + ван-дер-ваальсовы радиусы для атомов, находящихся на расстоянии не более 3 ангстрем от MyResidue.
f. Предыдущее изображение + ван-дер-ваальсовы радиусы для атомов, находящихся на расстоянии не более 4 ангстрем от MyResidue.
g. Предыдущее изображение + ван-дер-ваальсовы радиусы для атомов, находящихся на расстоянии не более 5 ангстрем от MyResidue.
h. Предыдущее изображение + ван-дер-ваальсовы радиусы для атомов, находящихся на расстоянии не более 6 ангстрем от MyResidue.
Предыдущее изображение + ван-дер-ваальсовы радиусы для атомов, находящихся на расстоянии не более 7 ангстрем от MyResidue. Для выполнения пп. c-i можно написать много раз одинаковые команды, различающиеся только расстоянием. Однако любители программирования могут использовать цикл for, примерно такого вида
for (var i = 1; i <= 7; i++ ) { # Тут пишем любые команды. # При этом можно использовать переменную i в таком виде: @{i}. # Например, within (@{i}, :a) или echo @{"text " + @{i}} }
На основании пп. c-i сделайте вывод о плотности расположения атомов в гидрофобном ядре. В отчете должны быть ответы на следующие вопросы:- На каком минимальном расстоянии от остатка расположены атомы, которые полностью покрывают его поверхность (так, что практически не видно атомов этого остатка)?
- Какое характерное расстояние между соседними не связанными ковалентно атомами в белке?
- Может ли между соседними атомами поместится еще один? Например, молекула воды, которая по сути - один большой кислород.
При расчете пользуйтесь следующими значениями ван-дер-ваальсовых радиусов атомов (источник: Дж. Эмсли. Элементы. 1993г.): углерод - 1.85 ангстрема, азот - 1.54 ангстрема, кислород - 1.4 ангстрема, сера - 1.85 ангстрема. Учтите, что радиус атома - не диаметр.
Для отображение используйте JMOL-апплет со скриптом:
- который последовательно показывает все перечисленные контакты
- показывает имена тяжелых атомов аминокислот, вовлеченных во взаимодейтсвие
- при помощи команды "echo" выводит строку, содержащую:
- тип контакта
- номера и трехбуквенные коды аминокислот, вовлеченных во взаимодейсвтие.
- Апплет также должен иметь кнопки "Запустить скрипт" и "Сохранить изображение".
- Ссылка на скрипт должна быть размещена на странице отчета.
- Там, где это разумно (например для водородных связей), необходимо отобразить геометрические параметры контактов.
В тексте отчета должна быть кратко описана процедура поиска контактов (какой бы они ни была) и соображения (в свободной форме) о функциональной роли данного контакта. Если вам удастся найти описание выбранного контакта в литературе, то будет просто замечательно.
Информация о лиганд-биомолекулярных контактах:
Если в выбранное структуре присутствуют низкомоллекулярные лиганды, то все аналогично. Обратите внимание, что контакты могут быть не только с белком, но и с нуклеиновой кислотой.
Информация о нк-белковых контактах:
Если в выбранное структуре присутствуют нуклеиновые кислоты, то все аналогично. Достаточно 1-2 примеров.
Личный вклад
В конце отчета добавьте раздел с описанием личного вклада участников. В качестве образца можете взглянуть на пример из Science
Author contributions: IF and SNE conceived the study. The initial dual choice assay was designed by SNE, IF, and RM; and performed by SNE and RM. RI, SRH designed, and SNE performed the two-choice olfactometer assay. RM, A-K B-K and GB designed, performed the GC-MS and analyzed the results. GB also designed the SPME collection of volatiles and designed the release of volatiles for bioassay. RM, PL and SNE designed and performed the respirometry experiments. SNE designed, and SNE and SH performed the mosquito fitness study. SNE performed parasite culturing, mosquito infections as well as the feeding proportion study and the statistical analyses. BL designed the temporal transcriptome study and prepared RNA for Illumina sequencing. SNE, BL, and SRH performed bioinformatics analyses of the sequencing data. BL performed the mosquito total protein assay. BL, SNE constructed figures and tables. IF, SNE and SH designed and SH created the artwork. IF, BL, RI, SRH and SNE wrote the manuscript. All authors proofread and commented on the final draft of the manuscript. The authors declare no competing financial interests.
Ссылки на источники:
Вся использованная литература или ресурсы.