Введение в биоинформатику. Современное положение. Задачи

advertisement
Сворачивание белка (folding)
“Principles that govern the folding of protein chain” Anfinsen, C. (1973) Science 181, 223230.
“the native conformation is determined by the totality of inter-atomic
interactions and hence by the amino acid sequences, in a given environment.
(solvent, pH, ionic strength, chemicals,etc)
Проблема сворачивания белка
• Функции белков определяются исключительно их 3D
структурой, конформацией
• Вопрос: можно ли предсказать 3D структуру белка, исходя
из данной аминокислотной последовательности?
• Ответ: в общем случае – нет!
• Но методы, которые позволяют разрешить частично 3D
структуры, бывают полезны.
Сворачивание – почему это так сложно?
• Линейные молекулы белков
предопределённые 3D структуры.
очень
быстро
сворачиваются
в
• Свойства любого белка определяются его 3D структурой.
• Белки могут денатурировать под воздействием химических веществ или
теплоты, но затем они сворачиваются вновь в исходную структуру.
Так почему же так сложно разрешить сворачивание?
• Структура белка может быть определена экспериментально (X-Rays or NMR)
но эта процедура не всегда возможна и не всегда даёт хорошие результаты.
• 3D структура последовательности определяется последовательностью
допустимых углов поворота, в которой каждый угол «состоит» из 2-х планарных
углов.
• Эта проблема может быть решена путём дискретизации (при некоторой потере
точности) путём ограничения количества возможных путей достижения каждой
из точек (координат атомов).
Если у вас нет структуры…
Подходы в предсказании
структуры
• Предсказание 1D:
–вторичные структуры
–доступность для растворителя
–трансмембранальные спирали
• Предскание 2D:
–контакты между аминокислотами/нитями бета.
• Предсказание 3D:
–моделирование гомологов
–распознавание фолда (e.g. via threading)
–ab initio предсказание (e.g. via молекулярная динамика)
Задачи
• Сравнение всех известных структур друг с другом
• Классификация и организация всех известных структур
• Поиск общих структурных шаблонов и мотивов
• Определение эволюционных расстояний между структурами белков
• Докинг – изучение взаимодействия между структурами
• Предсказание структур на основе последовательности
• Дизайн новых лекарств
Зачем?
• Первый шаг к предсказанию третичной структуры
• Один из основных элементов в распознавании фолда (для
моделирования далеких в эволюционном плане белков)
Предсказание вторичной структуры
http://www.new-science-press.com/
Table II of Williams, R.W. et al.: Biochimica
et Biophysica Acta 1987, 916:200-204.
Предсказание вторичной структуры
– предсказание для каждой аминокислоты в выбранном окне соседних
аминокислот (13-21)
– скоринг, обучение модели и предсказание 2D структуры (маппирование
элемента вторичной структуры на окно)
Методы
I. Chou-Fasman / GOR метод
II. Модели нейронных сетей
III. Методы «ближайшего соседа»
Метод Chou-Fasman (1974)
• Разработан Chou & Fasman в 1974 -1978
• База – известные 3D структуры глобулярных белков
– Частоты аминокислот в -спиралях
– Частоты аминокислот в -листах
– Частоты аминокислот в -поворотах
– Правила образования -спиралей и -листов
• Основан на растворимых, глобулярных белках – начальная
база 15 белков
Метод Chou-Fasman (1974)
Развитие Chou-Fasman
1. Присвоение каждой аминокислоте определенного пула
параметров
2. Идентификация a-helix и b-sheet. Удлинение этих областей в
обоих направлениях.
3. При перекрытии – сравнение P(H) и P(E) и скоринг.
Развитие Chou-Fasman
1.
Вероятности
T
P(H)
69
P(E)
147
P(turn) 114
S
P
T
A
E
L
M
R
S
T
G
77
75
143
57
55
152
69
147
114
142
83
66
151
37
74
121
130
59
145
105
60
98
93
95
77
75
143
69
147
114
57
75
156
P(H)=[(#H in helix)/(#H)]/(fraction helix {all})
Поиск a-спирали
2.
Поиск областей, где 4 из 6 аминокислот имеют P(H) >100 (“ ядро aспирали”)
P(H)
P(H)
T
S
P
T
A
E
L
M
R
S
T
G
69
77
57
69
142
151
121
145
98
77
69
57
T
S
P
T
A
E
L
M
R
S
T
G
69
77
57
69
142
151
121
145
98
77
69
57
Удлинение ядра a-спирали
3.
Расширение области ядра, пока 4 аминокислоты имеют среднее P(H)
>100.
P(H)
T
S
P
T
A
E
L
M
R
S
T
G
69
77
57
69
142
151
121
145
98
77
69
57
Поиск β-листа
4. Поиск областей, где 3 из 5 аминокислот имеют P(E) >100 “ядро β-листа”
5. Удлинение ядра до тех пор, пока 4 соседних аминокислоты имеют
среднее P(E) > 100
6. Если score области > 105 и среднее P(E) > среднее P(H), значит эта
область - β-лист
GCG Programs
•
PepPlot
– Plot on parallel panels
– -cff option, text output
•
PeptideStructure
– text output (Most useful for detail)
•
PlotStructure
– two outputs
• squiggles “protein-like”
• parallel panels
GOR III (Garnier-Osguthorpe-Robson. Gibrat J.F., J.Mol.Biol, 1987)
Модели нейронных сетей
-
Машинное обучение
Сет структур (e.g. a-спирали, не a-спирали)
Обучение распознавать шаблоны, структуры в известных белках
Эффективность ~ 70 –75%
NPS@
сервер
Rost B. Sander C. Prediction of Protein Secondary Structure at
Better then 70% Accuracy. J.Mol.Biol., 1993, vol. 232. 584-599.
Сервисы и алгоритмы
Eva
Предсказание вторичной структуры
Предсказание вторичной структуры
Predict Protein (Mega) - secondary structure ( PHDsec, and PROFsec)
PSI-pred (PSI-BLAST profiles used for prediction; David Jones, Warwick)
PHD - Rost & Sander, EMBL, Germany
ASPSSP server Raghava, INDIA
DSC - King & Sternberg (this server)
PREDATOR - Frischman & Argos (EMBL)
ZPRED server Zvelebil et al., Ludwig, U.K.
nnPredict Cohen et al., UCSF, USA.
BMERC PSA Server Boston University, USA
SSP (Nearest-neighbor) Solovyev and Salamov, Baylor College, USA.
•
•
JPRED Consensus prediction (Cuff & Barton, EBI)
NPS@
Предсказание вторичной структуры
Предсказание вторичной структуры
Предсказание функции
Еще одна важная задача протеомики — анализ и
предсказание функции белка. Известно, что функция
белка определяется его активными сайтами, поэтому
накопление и систематизация информации об активных
сайтах белков чрезвычайно актуальна. В. Иванисенко,
Д. Григоровичем и С. Пинтусом разработана
компьютерная база данных PDBSite, которая содержит
информацию о более чем 12 тысячах активных сайтов
белков. Источником информации служат хорошо
документированные пространственные структуры
белков.
Основные методы в
биоинформатике:
• парное выравнивание;
• множественное выравнивание;
• поиск гомологов, threading;
• структурное выравнивание.
CASP
Critical Assessment of Techniques for Protein Structure Prediction
CASP1 (1994) CASP2 CASP3 CASP4 CASP5…..CASP9 (2010)
• Comparative modeling (CM)
• Fold-recognition (FR)
• CAFASP meta-server ver. 3
• New folds (NF)
• Ten most wanted sec. struct. contacts, protein-protein docking,
and disordered predictions.
About CASP: CASP is a blind study/experiment that aims at establishing the current state of the art
in protein structure prediction; identifying what progress has been made; and highlighting where
future effort may be most productively focused (Every two years).
This blind study is held over an ~8 month time period and ends in a meeting held every two years,
in Asilomar, CA, starting from 1994. For the procedure of the experiment, CASP participants are first
provided target sequences (around May) via the Protein Structure Prediction Center at Lawrence
Livermore National Laboratory. The participants have a few months to determine the template
structure, alignment, model structure and evaluate their results.
The sequence targets are categorized by homology and difficulty for predicting their structure. The
fairly simple targets have med. sequence homology (>30% seq. identity) are considered comparative
modeling (CM) predictions; the med. difficulty targets have med.-to-low sequence homology (~1030% seq. identity) are considered fold-recognition (FR) predictions; and the difficult targets have low
seq. homology and usually require an ab initio methods are considered new folds (NF).
During the prediction time (~May-Oct.), researchers (structural biologist in x-ray or NMR) work
on solving the experimental structure of each of the target sequences and they hold back the structure
coordinate information from the predictors. By Nov., all participants submit their models (as
coordinates) to the Livermore Center and the researchers (who solve the target structure) finalize and
post their results. Finally, in Dec., all participants and the CASP organizers meet to evaluate the
results of the experiment comparing each model with the experimental structure and discussing the
methodologies used.
The goal of CAFASP is to evaluate the performance of fully automatic structure prediction servers
available to the community. In contrast to the normal CASP procedure, CAFASP aims to answer the
question of how well servers do without any intervention of experts, i.e. how well ANY user using
only automated methods can predict protein structure. CAFASP assesses the performance of methods
without the user intervention allowed in CASP.
CASP
CASP отчеты, результаты
Предсказание сворачивания белка vs
предсказание структуры
Престказание процесса фолдинга белка связано с процессом приобретения
белком его 3D формы, очертаний – физико-химические принципы.
Предсказание структуры – используются любые статистические,
теоретические и эмпирические данные.
4 подхода:
Моделирование гомологов (Homology Modeling)
Распознавание фолда (Sequence-Structure Threading (secondary structure
prediction)):
• Dynamic programming
• Knowledge-based potentials
Предсказание Ab initio
Docking and Drug Design
Техники фолдинга белка
• Моделирование гомологов (homology
modeling)
• Ab initio предсказание
• Распознавание сворачивания “Threading'‘
• Докинг
Сравнительное моделирование гомологов
Для последовательностей с гомологичностью > 25-30%
использовать известную PDB структуру как отправной пункт
для
создания
3D
модели
структуры
неизвестной
последовательности.
Нужно
использовать
координаты
основной
цепи
(N-Cα-C) гомологичной структуры как шаблон для модели
70% и более гомологичности – очень высокое качество
модели, даже положения боковых цепей могут быть предсказаны
с высокой точностью.
40%-65% - средняя точность предсказания. Могут быть
серьёзные ошибки даже в положении основной цепи, особенно в
областях петель\изгибов.
Лекарственные средства, разработанные с использованием
неверных представлений о структуре белка, могут быть
токсичны или обладать неучтёнными побочными эффектами.
Для эффективности этого метода требуется по меньшей мере
3,000 уникальных, совершенно точно определённых
структур. На конец 2001 года имелось только 1,000
уникальных структур среди 16973 в PDB. 2008 год – 53000
структур в PDB, homo sapience ~1500.
Сравнительное моделирование гомологов
1. Последовательность-цель – первичная структура
белка, 3D структуру которого следует определить
2. Шаблон – белок, чья 3D структура ясна
3. Выравнивание последовательностей 1 и 2
Желательно также иметь
Биохимическую и структурную информацию
(литература)
Дополнительные последовательности
гомологов с известной структурой
Моделирование гомологов
1. Fragment-based modeling:
Выравнивание с целью идентификации структурнопостоянных областей (SCR): а) области без вставокделеций и в) области с четко определяемой вторичной
структурой. VR – области между SCR. Composer (Sybil),
Homology (InsightII)
2. Restraint-based modeling:
получение score-функции путём комбинирования
«ограничений» - расстояний между Сα, торсионных
углов и т.д. Оценка результатов MD данной scoreфункцией. Modeller
Swiss-PDBViewer
• Моделирование гомологов
• Предсказание Ab initio
• Распознавание сворачивания “Threading'‘
• Докинг
Предсказание Ab initio
Применяется, когда неизвестны гомологи, нет структуры,
которую можно было бы использовать как шаблон
Есть только одна последовательность. Предсказание 3D
основано на «базовых» принципах, таких, как энергетические
и статистические законы и правила.
Это – симуляция физических сил и процессов, которые
могут привести развёрнутый белок в нативную (стабильную,
присутствующую в природе) конформацию на компьютере
Стабильность с точки зрения термодинамики: нативная
конформация белка есть его глобальный минимум свободной
энергии. Белок должен сворачиваться так самостоятельно.
Предсказание Ab initio
Электростатические
Ван-Дер-Ваальс
Водородные связи
Энергия торсионных связей
Предсказание Ab initio - сворачивание
Полный расчёт энергий – очень затратный с точки
зрения вычислений процесс.
Поэтому требуется разработка неких эвристических
энергетических функций, которые бы надёжно
различали «правильную» и «неправильную» структуры
и лучше «понимали» бы силы, которые управляют
сворачиванием белка.
Folding
139 минимумов
Тетрапептид аланина
Фолдинг. Предсказание Ab initio
Protein Folding: A Perspective from Theory and Experiment
Christopher M. Dobson,* Andrej Sœ ali, and Martin Karplus*
Предсказание Ab initio
Сравнение расчётной и экспериментальной модели для белка миоглобина и
использованием refined potential function. Рассчитанная структура является 3D
структурой, полученной в результате 3-х разных расчётов с дальнейшей
кластеризацией и выбором структуры с наименьшей энергией. Общее время
симуляции на кластере из 16 машин CM-5 massively parallel computer составило 60
часов, в течении которых было генерировано порядка 5 миллионов структур. RMS
составляет 6.2 Å.
Парадокс Левенталя
Время, за которое белок скручивается,
(принимает конечное 3D состояние) на много
порядков меньше времени перебора всех
возможных конфигураций.
Допустим, в белке 100 атомов, каждый из которых принимает 3 положения:
3 100 = 5 × 10 47 конформаций.
Наибыстрейшее движение – 10- 15 с. Перебор всех конформаций займёт
5 × 10 32 с или 1.6 × 10 25 лет (возраст Вселенной ~ 13,75 × 109)
• Homology Modeling
• Ab initio prediction
• Fold Recognition or “Threading'‘
Распознавание сворачивания (“Threading”)
Напоминает метод моделирования гомологов, но не требует
структур с высокой степенью идентичности.
Интересующая нас последовательность «протягивается»
через все возможные позиции основной цепи во всех
известных белковых структурах в PDB, и для каждой итерации
рассчитывается её свободная энергия.
Структура, которая даст лучший показатель энергии
принимается за «шаблон» и дальнейший процесс напоминает
моделирование гомологов
Threading не может быть применён для тех белков, для
которых в базе PDB нет похожих структур.
Распознавание сворачивания (“Threading”)
Profiles-3D scoring function: оценка
локального структурного
выравнивания (укладки) каждой
аминокислоты в последовательности без учета попарного взаимодействия аминокислот+склонность к
H/E/L структурам+полярность
(solvent exposure)
Из «Methods in Molecular Biology, vol 143, Methods and ProtocolMethods and Protocols.
Protein Structure Prediction, еdited by David M. Webster»
Распознавание сворачивания (“Threading”)
Рисунок из R. Lathrop et al, “Analysis and Algorithms for Protein Sequence-Structure Alignment” in
Computational Methods in Molecular Biology, Salzberg et al. editors, 1998.
Fold Recognition – The Fold
PDB
Groups clustered
by a common
resemblance
Calculated
Folds
Сколько
всего
фолдов?
Genome Sequencing
Homology
Structure Conservation
Количество фолдов ~ 4000
БД из 930 фолдов ~ 90% семейств белков
Fold Recognition – недостатки
Этот метод редко приводит к тому качеству структурного
выравнивания, которое предоставляет моделирование гомологов.
Серверы
PredictProtein Server•
ModBase (a database of three-dimensional protein models •
calculated by comparative modeling)
3D PSSM & ModBase
3D-PSSM предсказание 3D структуры по последовательности и вероятность
этой структуры
ModBase – база данных 3D структур, построенных на основе сравнительного
моделирования
Download