Прямое компьютерное моделирование взаимодействиябелков Прямое компьютерное моделирование взаимодействия белков И.Б. Коваленко, А.М. Абатурова, О.С. Князева, Г.Ю. Ризниченко, А.Б. Рубин МГУ им. М.В.Ломоносова Биологический факультет, кафедра биофизики

План 1. Компьютерная модель взаимодействия белков 2. Моделирование реакции белков пластоцианина и цитохрома f в растворе 3. Моделирование реакции белков пластоцианина и цитохрома f в клетке 4. Предсказание областей связывания белков

Модель 3D кубический реакционный объем Много частиц Pc и cyt f (многочастичная модель) 3D броуновское движение белков Электростатические взаимодействия Модель

Взаимодействие белков 1. Диффузия 2. Электростат. взаимодействие, сближение и взаимная ориентация 3. Предварит. комплекс 4. Конформационные изменения, финальный комплекс Модель

1. Диффузия 2. Электростат. взаимодействие, сближение и взаимная ориентация 3. Предварит. комплекс 4. Конформационные изменения, финальный комплекс Взаимодействие белков: диффузия Модель

Уравнение Ланжевена Большой размер молекул белков по сравнению с молекулами растворителя Влияние молекул растворителя на движение белка через ξ (броуновское движение и сила трения) При больших значениях ξ можно пренебречь инерцией: Диффузия белков Модель

Уравнения Ланжевена для броуновского движения частиц в электростатическом поле Поступательное движение Вращательное движение x – координата, вдоль которой рассматривается движение, ξ t x – коэффициент вязкого трения вдоль этой координаты, f x (t) и F x – проекции случайной и электростатической сил на ось x, соответственно φ- угол поворота вокруг оси, относительно которой рассматривается движение, ξ r x – коэффициент вязкого трения для вращательного движения вокруг оси x, m x (t) и M x – моменты случайной и электростатической сил относительно оси x, соответственно k – постоянная Больцмана, T – температура φ-потенциал Модель (Ermak and McCammon, 1978)

Коэффициенты вязкого трения для эллипсоида вращения поступательное движение вращательное движение b = c -вдоль оси а -вдоль оси b -вокруг оси а -вокруг оси b Модель (Perrin, 1936)

Рс M = 10.5 КДа a=21 Å, b=14 Å Молекулярная масса Оси эллипсоида вращения Cyt f M = 27.9 КДа a=47 Å, b=17 Å Аппроксимация цитохрома f и Рс эллипсоидами вращения для нахождения коэффициентов вязкого трения Модель

Взаимодействие белков: электростатические взаимодействия 1. Диффузия 2. Электростат. взаимодействие, сближение и взаимная ориентация 3. Предварит. комплекс 4. Конформационные изменения, финальный комплекс Модель

Заряженные аминокислотные остатки на поверхности белка Экранирование Уравнение Дебая-Хюккеля При ионной силе 100 мM радиус Дебая равен 8 А На расстоянии 35А за счет экранирования ионами соли потенциал падает на 2 порядка Электростатические взаимодействия белков Модель

Белки представлены как области с низкой диэлектрической проницаемостью Раствор представлен как область с высокой диэлектрической проницаемостью с зарядами (ионами) Белки рассматриваются как твердые тела Модель Электростатическая модель

Расчет потенциала по уравнению Пуассона-Больцмана φ потенциал, ε диэлектрическая проницаемость, ρ prot плотность зарядов в белке, с i bulk концентрация i-иона в растворе, Z i заряд i-иона, e заряд электрона, R молярная газовая постоянная, T температура в К Линеаризованное уравнение Пуассона-Больцмана I ионная сила, κ величина, обратная длине экранирования потенциал, заряд, κ в центре 0-ячейки потенциалы в соседних ячейках диэлектрические проницаемости на границе ячеек из G. M. Ullmann (2004) Модель · ·

уравнение Пуассона-Больцмана ε белка =2 уравнение Дебая-Хюккеля ε белка =80 Сравнение методов расчета потенциала Эквипотенциальные поверхности для восстановленного cytf при ионной силе mM, pH=7, ε р-ра =80; красный цвет -6.5 мВ, синий +6.5 мВ; зеленым цветом обозначены атомы молекулы. Модель

Взаимодействие белков: предварительный комплекс 3. Предварит. комплекс 4. Конформационные изменения, финальный комплекс 1. Диффузия 2. Электростат. взаимодействие, сближение и взаимная ориентация 3. Предварит. комплекс 4. Конформационные изменения, финальный комплекс Модель

Предварительный комплекс Заряды на поверхности Подвижность аминокислотных остатков на поверхности белков Модель (Kleanthous, 2000) (переходное состояние, которое быстро переходит в финальный комплекс) Белок 1 Белок 2 riri

Аппроксимация Рс и cyt f набором сфер для моделирования столкновений Цитохром f Пластоцианин Молекулы показаны зеленым, сферы показаны синим Модель

Взаимодействие белков: финальный комплекс 1. Диффузия 2. Электростат. взаимодействие, сближение и взаимная ориентация 3. Предварит. комплекс 4. Конформационные изменения, финальный комплекс Модель

Финальный комплекс Взаимодействия Ван-дер-Ваальса Гидрофобные взаимодействия Солевые мостики Водородные связи Модель

Модели взаимодействия белков пластоцианина (Pc) и цитохрома f (cyt f) в растворе Модель

Электростатический поверхностный потенциал (справа) эквипотенциальные поверхности (слева) (10 мВ) пластоцианина, pH=7, I=100 mM

Электростатический поверхностный потенциал цитохрома f, pH=7, I=100 mM

Изменение потенциала Рс при мутациях модель Ǻ Ǻ Сайты связывания с cyt f WT

модель Изменение потенциала Рс при мутациях Изменения зарядов На аминокислотных остатках Ǻ Ǻ Сайты связывания с cyt f Glu59Lys/Glu60Gln

модель Изменение потенциала Рс при мутациях Ǻ Ǻ Сайты связывания с cyt f Изменения зарядов На аминокислотных остатках Glu59Lys/Glu60Gln/ Glu43Asn

Время Кривая, полученная фиттингом по закону действующих масс Кривая, полученная на модели Концентрация P 1 P 2 Реакция, которую мы моделировали: k V = k[P 1 ][P 2 ] Reaction velocity: После того, как моделирование проведено, нужно вычислить константу скорости реакции белков Определение k фиттингом k Результат моделирования: The model P1+P2 P1P2

результаты Реакция между cyt f и разными мутантными Pc в растворе Зависимость логарифма константы скорости от корня из ионной силы экспериментальные данные из A. Kannt et al.(1996) результаты моделирования Gln88Glu Glu59Lys/Glu60Gln/ Glu43Asn Glu59Lys/Glu60Gln/ Glu43Asn Glu59Lys/Glu60Gln Glu43Gln/Asp44Asn Glu43Lys Glu43Asn Asp42Asn результаты моделирования Kovalenko et al., Phys. Biol. 3, 121–129.

Для вычисления и хранения электрического потенциала небольшого белка необходимо 100 МБ оперативной памяти Для вычисления константы скорости реакции необходимо выполнить 10 7 шагов (шаг с) На каждом временном шаге необходимо вычислить электростатическую и броуновскую силы, переместить и повернуть каждый белок, проверить факты столкновения и реакции белков (тысячи операций для каждого белка) Объемы вычислений P 1 +P 2 P 1 P 2 k

Архитектура хлоропласта

Фрагмент тилакоида. Показаны тилакоидные мембраны и люминальное пространство 100 Å 50 Å

Расположение Cytb 6 f комплекса в мембране и связывание с его субъединицей f (cyt f) белка- подвижного переносчика пластоцианина (Pc) Pc Cyt f PDB ID – 2PCF; 1Q90

Визуализация люмена, ограниченного сверху и снизу тилакоидными мембранами (заштрихованы) и белков пластоцианина и цитохрома f.

Зависимость константы скорости реакции между Pc и cytf в люмене тилакоида при постоянном количестве молекул Pc и cytf от расстояния между мембранами (площадь мембран постоянна =322 х 322 нм 2 ) от длины мембраны (корня из площади мембраны) при постоянном расстоянии между мембранами 16 нм z x x Молекулы цитохрома расположены на мембране Молекулы цитохрома расположены в объеме люмена Молекулы цитохрома расположены в кубе, объем которого равен объему люминального пространства, растворе 16 нм Kovalenko et al., Biophysics 53,

Эквипотенциальные поверхности (6.5 mV) в люмене тилакоида хлоропласта, pH=7, I=100 mM Knyazeva, O.S., Kovalenko, I.B. et al. Biophysics 55, 259–268 (in Russian).

Предсказание областей связывания белковых молекул

Tested protein pairs First proteinSecond protein Barnase +Barstar + Colicin E9 DNase +Im9 + AChE +fasciculin2 – Thrombin + -Thrombomodulin +- Erythropoietin + Epo receptor 1 st domen- Erythropoietin + Epo receptor 2 nd domen+ Colicin E3 RNase -Im3 - Interleukin 4+IL4Binding Protein+ Lipoprotein (pal)Transport pro (tolB) Plastocyanin+Cytochrome f+

barnase barstar Approximation of the proteins barnase and barstar by ellipsoids for calculation of viscous friction coefficients

Electrostatic potential of barnase and barstar barnase barstar

Представление рассчитанной на модели вероятности связывания молекул барназа (справа) и барстар (слева) в виде цветного «глобуса» распределения вероятности при расстоянии между центрами масс белков 40 А; синий цвет ячеек соответствует малой вероятности связывания данной ячейкой, красный цвет – высокой вероятности. Kovalenko et al., Dokl. Biochem. Biophys. 427, 215–217.

Представление рассчитанной на модели вероятности связывания молекул пластоцианина и цитохрома f в виде цветного «глобуса» распределения вероятности при расстоянии между центрами масс белков R=46 A; синий цвет ячеек соответствует малой вероятности связывания данной ячейкой, красный цвет – высокой вероятности. Kovalenko et al., Dokl. Biochem. Biophys. 427, 215–217.

Спасибо за внимание!