СИНТЕЗ АТФ F1-АТФАЗОЙ В СТОХАСТИЧЕСКОЙ МОДЕЛИ

advertisement
СИНТЕЗ АТФ F1-АТФАЗОЙ В СТОХАСТИЧЕСКОЙ
МОДЕЛИ
Погребная А. Ф.
(Россия, Москва)
Данная работа является продолжением цикла работ [1-4] посвященных построению математической модели вращающегося молекулярного мотора F1-АТФазы. В данной работе, в рамках представленной ранее модели рассматривается синтез
АТФ при вращении ротора молекулярного мотора под действием внешней силы.
F1-АТФаза, катализирующая синтез/гидролиз АТФ, — часть
фермента F0F1-АТФсинтазы в энергопреобразующих мембранах
живых клеток. Энергетически выгодный процесс переноса протонов водорода через мембранную часть F0F1-АТФсинтаза — комплекс F0 под действием электрохимического потенциала и синтез
АТФ в F1-АТФазе сопряжены благодаря наличию ротора механически связывающего F0 и F1 комплексы. Протонный ток через F0
вращает вал и является причинной конформационных изменений
в каталитических центрах F1-АТФаза, что приводит к синтезу
АТФ [5]. С другой стороны гидролиз АТФ в F1-АТФазы приводит
к вращению вала в противоположную сторону и создает поток
протонов через мембрану в обратном направлении. Отделенный
от мембранной части комплекс F1-АТФаза сохраняет свою каталитическую способность к гидролизу и синтезу АТФ [6].
Структурные свойства и молекулярный механизм работы F1АТФазы изложены в многочисленных статьях [см., например,
6]. Характерный размер этого молекулярного мотора около 100
ангстрем. Белковый комплекс F1-АТФаза состоит из (αβ)3γ
субъединиц, где центральная субъединица γ окружена цилиндром из α3β3. Кристаллографический анализ показывает, что все
три β каталитические субъединицы находятся в различных кон-
Погребная А. Ф. — МКО — 2005, ч. 3, стр. 944 – 945
Pogrebnaya A. F. — MCE — 2005, vol. 3, p. 944 – 945
формациях: плотно закрытой — (βt), неплотно закрытой — (βl),
открытой — (βe) (рис. 1). При этом активный центр (АЦ) каталитической субъединицы в плотно закрытой конформации (βt)
связан с АТФ, АЦ полузакрытой конформации (βl) связан АДФ
и АЦ открытой конформации (βe) — пуст. Согласованное изменение всех трех конформаций, сопряженное с гидролизом АТФ,
приводит к вращению «ротора» — центральной γ субъединицы
и, наоборот, вращение изогнутой несимметричной γ приводит к
конформационым изменениям в каталитических субъединицах,
что приводит к синтезу АТФ.
А
С
В
Рис. 1. β и γ субъединицы (по данным рентгеноструктурного анализа) в трех
различных конформациях; А — β субъединица связанная с MgATP (закрытая
конформация βt); В — β субъединица со свободным каталитическим центром
(открытая конформация βe); С — β субъединица связанная с MgADP (конформация βl)
В настоящее время вращение и при гидролизе и при синтезе
наглядно экспериментально продемонстрировано. В этих экспериментах к основанию субъединицы γ, выступающему из ансамбля α3β3, прикреплялась нагрузка (фрагмент нити актина или
полистироловая бусинка) с флуоресцирующей меткой. Средняя
скорость вращения ротора в режиме гидролиза была в пределах
0.5–10 Гц (в зависимости от нагрузки), с шагами в 1200, а КПД
мотора близко к 100%. В более поздних экспериментах, выполненных с большим во времени разрешением, было показано, что
945
Раздел 8. Математические методы в биологии, экологии и химии
Part 8. Mathematical methods in biology, ecology and chemistry
повороты ротора на 1200 происходят ступенчато и включают в
себя две фазы: поворот на 800–900, после чего наступает пауза,
затем следует второй поворот на 400–300. В работе [7] представлены результаты по наблюдению синтеза АТФ изолированным
ферментом F1-АТФазы при вращении ротора по часовой стрелке
с помощью внешнего вращающего момента. В этом эксперименте к ротору пришивалась намагниченная бусинка и с помощью
электромагнитов вращалась по часовой стрелке (в сторону соответствующую синтезу АТФ).
Настоящая работа является продолжением наших ранних работ, в которых мы представили достаточно простую динамическую модель F1-АТФазы, описывающую стохастическое движение центральной γ субъединицы в силовом поле, образованном
субъединицами комплекса (αβ)3. Для описания потенциального
рельефа, в котором реализуется движение γ субъединицы, были
определены энергии связи различных субстратов, участвующих в
работе мотора, с активным центром, конформация которого определяется положением γ субъединицы. Наша модель позволяет
описать направленное вращение ротора, как при синтезе, так и
при гидролизе АТФ, а также ответить на вопросы, почему при
работе мотора реализуются паузы именно в положениях 900, 1200,
2100, 2400, и т.д., а также, почему эффективность мотора близка к
100% [3]. В данной работе мы представляем результаты по моделированию синтеза АТФ в рамках нашей модели.
Описание модели. Прежде чем представить уравнения
модели опишем события, происходящие в моторе под действием
внешнего вращающего момента. Изогнутая γ единица вращается
внутри комплекса (αβ)3 и является причиной конформационных
изменений в каталитических субъединицах. Она насильно сгибает подвижную часть β при повороте на 300 и полностью раскрывает активный центр (АЦ) при повороте на 900. Мы предполагаем, что посадка Pi предшествует входу АДФ в карман, но не
приводит к изменению конформации каталитической субъединицы. Далее в открытый АЦ входит АДФ, в каталитической
субъединице происходит переход β e → β l (см рис. 2). Поворачиваясь на 300, γ субъединица переводит βl в βt. Если повороту
946
Погребная А. Ф. — МКО — 2005, ч. 3, стр. 944 – 947
Pogrebnaya A. F. — MCE — 2005, vol. 3, p. 944 – 947
ротора на первые 300 предшествует посадка реагентов (АДФ и
Pi), то после поворота реагенты настолько близки друг к другу,
что становится возможна реакция синтеза АТФ. Во время этого
же поворота другой своей стороной γ субъединица «раскрывает» (переход β t → β e ) предыдущую субъединицу, содержащую в своем АЦ синтезированную АТФ. Таким образом, за поворот γ субъединицы на 1200 происходит сближение реагентов,
реакция синтеза и открытие каталитического центра.
Рис. 2. VATP (θ ) — силовой потенциал образованный каталитической
субъединицей содержащей в АЦ молекулу АТФ, VADP (θ ) — силовой
потенциал образованный каталитической субъединицей содержащей в
АЦ молекулу АДФ; S (θ - θi ) — «открытость АЦ», R(θ - θi ) — реакционная способность АЦ. Более подробно получение и параметры потенциала описаны в работе [3]
Основные уравнения. Для описания движения ротора
мы используем уравнение Ланжевена, которое имеет вид:
ξ
∂V (θ )
∂θ
=− s
+ fl (t ) + Fext ,
∂t
∂θ
где θ — задает угловое положение центральной γ субъединицы; ξ — коэффициент вязкого трения, возникающий при вра947
Раздел 8. Математические методы в биологии, экологии и химии
Part 8. Mathematical methods in biology, ecology and chemistry
щении нагрузки; Vs ( θ ) — силовой потенциал, являющийся аддитивной суммой силовых потенциалов создаваемых всеми
тремя β каталитическими субъединицами, находящимися в положениях θ1 = 00 ,θ 2 = 1200 ,θ3 = 2400 :
Vs (θ ) = V (θ − θ1 ) + V (θ − θ 2 ) + V (θ − θ3 ) ;
Fext — внешний вращающий момент; f1 (t ) — момент вращения
создаваемый силами стохастической природы, обусловленными
тепловыми флуктуациями. При этом
fl (0) fl (τ ) = ξ kT δ (τ ) .
Учитывая, что характерные частоты вращения γ субъединицы
составляют десятки герц и, учитывая значительную вязкость
среды (при η ≈ 10−3 Nsm -2 ξ ≈ 10−21 [3–6]) и сравнительно малую
величину момента инерции ( J = mL2 ≈ 10−32 ), при составлении
d 2θ
уравнения (1) мы опустили член J 2 , отвечающий за инерdt
цию системы.
Изменение деформационного потенциала создаваемого каталитическими центрами, в общем случае определяется присутствием в кармане того или иного комплекса в нем, т.е.:
V (θ - θi ) = ψ ATPVATP (θ ) + ψ ADPVADP (θ )
где ψ ATP ,ψ ADP ∈ {[1,0,] , [ 0,1] , [ 0, 0]} — определяют наличие или
отсутствие того или иного субстрата в активном центре, а
VATP (θ ), VADP (θ ) — определяют потенциалы, возникающие при
присутствии субстрата, соответствующего индексу потенциала
(рис. 2). Эта запись потенциала учитывает наше предположение,
что выход Рi происходит сразу после гидролиза.
«Химический цикл» — посадка, выход и преобразование субстратов в АЦ — мы описываем с помощью системы кинетических уравнений для вероятности заселенностей активных
центров субстратом, которые однозначно определят вероятности ψ ATP ,ψ ADP в уравнении (3):
948
Погребная А. Ф. — МКО — 2005, ч. 3, стр. 944 – 949
Pogrebnaya A. F. — MCE — 2005, vol. 3, p. 944 – 949
dψ ATP ,i
dt
dψ ADP,i
dt
⎛ S (θ − θi ) R(θ − θi ) ⎞
R(θ − θi )
ψ ADP,i + k+ ATP nATP S (θ − θi ),
= −⎜
+
⎟ψ ATP ,i +
τh
τs
(3)
⎝ τ − ATP
⎠
⎛ S (θ − θi ) R(θ − θi ) ⎞
R(θ − θi )
ψ ATP,i + k+ ADP nADP S (θ − θi ),
= −⎜
+
⎟ψ ADP ,i +
τs
τh
⎝ τ − ADP
⎠
где i = 1, 2,3 — номер каталитической субъединицы. Уравнение
3 описывает динамику изменения вероятности заполнения активных центров ψ ATP ( ADP ),i тем или иным субстратом (АТФ или
АДФ), причем, учитывая, что карман может быть заселен только
одной молекулой субстрата, на вероятности ψ ATP ( ADP ),i наложены
ограничения:
0 ≤ ψ ATP ( ADP ),i ≤ 1 ,
ψ ATP ( ADP ),i + ψ ADP ( ATP ),i ≤ 1 .
Константа k ATP ( ADP ) — задает скорость посадки субстрата в полностью открытый центр; nATP и nADP — концентрации АТФ и
АДФ во внешней среде; S (θ − θi ) , R(θ − θi ) — функции, описывающие «открытость» и «реакционную способность» активных
центров в зависимости от положения γ субъединицы (см. рис. 2);
1 τ DT ,1 τ TD — определяют скорость гидролиза и синтеза АТФ
в полностью закрытом центре; 1 τ ATP ( ADP ) — скорость выхода
продуктов реакции из АЦ.
В уравнении (3) первое слагаемое описывает посадку и уход
субстрата из активно центра, скорость такого процесса ограничивается конформацией β субъединицы (функция S (θ − θi ) ).
Второй член описывает реакции гидролиза и синтеза в закрытом
центре. Область конформационных состояний, определяемых
координатой θi , при которых протекает реакция, ограничивается функцией R(θ − θi ) .
Полученная система уравнений позволяет нам полностью
описать динамику работы всего мотора.
Результаты численного моделирования и выводы. При
численном решении системы (1)–(5), мы полагали, что в начальном состоянии АЦ свободны от субстратов. Во всех расчетах ξ = 10−21 Нмс , что соответствует нити филамента актина
949
Раздел 8. Математические методы в биологии, экологии и химии
Part 8. Mathematical methods in biology, ecology and chemistry
длиной
1мкм,
k+ ATP = 4 × 106 c −1 M −1 ,
τ − ADP k+ ATP = 0.25M −1 ,
k+ ADP
тываем
вероятностей
τ − ATP k+ ATP = 0.12 M −1 ,
= 0.5M −1 . При моделировании
k+ ATP
процесса синтеза мы полагали концентрацию АТФ во внешней
среде равной нулю.
Для численного решения уравнения (1) мы воспользовались
простой разностной схемой. Для моделирования шумового воздействия теплового резервуара на вращение центральной субъединицы на каждом шаге по времени ( Δt = 0.01 ÷ 0.1 мс ), мы генерировали случайное число с гауссовским распределением. В
уравнении (4) левые части аппроксимировались разностью и
далее решение уравнений разбивалось на два этапа, т.е. использовался метод разбиения по «физическим» факторам. На первом
шаге Δt / 2 производился расчет вероятности выхода/посадки
субстрата из/в активного центра, а на втором шаге, производился учет химической реакции внутри АЦ.
Подробнее о методике численного расчета можно прочесть в
наших более ранних работах [см., например, 3], отметим только
что, определяя начальные вероятности ψ ATP ( ADP ),i ;t = 0 , мы рассчизначения
на
следующем
шаге
ψ ATP ( ADP ),i ;t =Δt / 2 и ψ ATP ( ADP ),i ;t =Δt , по этим значениям вероятностей
определяем состояние активного центра (заселенность его тем
или иным субстратом). В случае ψ ATP ( ADP ),i ;t =Δt / 2 > μ
(ψ ATP ( ADP ),i ;t =Δt > μ ), где μ ∈ [ 0,1] случайное число с равномерным законом распределения на интервале [ 0,1] , активный центр
считался
заселенным
субстратом,
т.е.
(
)
ψ ATP ( ADP ),i ;t =Δt / 2 ψ ATP ( ADP ),i;t =Δt ←⎯
⎯ 1 (коэффициенты в выраже-
(
)
⎯0.
нии 3), в противном случае ψ ATP ( ADP ),i ;t =Δt / 2 ψ ATP ( ADP ),i ;t =Δt ←⎯
На следующем этапе по заданному состоянию мотора производилось решение уравнения Ланжевена (1), в котором значения
потенциалов определялись исходя из полученных значений
950
Погребная А. Ф. — МКО — 2005, ч. 3, стр. 944 – 951
Pogrebnaya A. F. — MCE — 2005, vol. 3, p. 944 – 951
ψ ATP ( ADP ),i ;t =Δt (0 или 1). На следующем шаге по времени последовательность операций повторялась. Таким образом, на базе
кинетических уравнений, мы могли моделировать поведение
отдельного мотора (одной «реализации»). Для описания усредненных характеристик (например, скорости химических реакций, осуществляемых каталитическими комплексами), необходимо провести моделирование ансамбля моторов, т.е. провести
усреднение по множеству отдельных реализаций.
На рис. 3 представлены полученные нами реализации зависимости угла поворота γ субъединицы от времени при различных концентрациях АДФ, знак минус указывает на направление
вращения по часовой стрелке (в сторону противоположную той,
что при синтезе). Сверху показано количество молекул АТФ,
покидающих один из трех АЦ. Видно, что с ростом концентрации АДФ уменьшается длительность пауз перед каждым 1200
шагом. На рисунке 4 представлена зависимость средней скорости вращения и синтеза от концентрации АДФ в окружающей
среде. Сплошная линия на рисунке является аппроксимацией
точек гиперболой, искомые значения
K m ≈ 0.187 мМ ,
vmax ≈ 22.8c −1 . К сожалению, мы не можем сравнить полученные
значения с экспериментальными данными.
Рис. 3. Реализации зависимости угла
поворота γ субъединицы от времени
при различных концентрациях АДФ
Рис. 4. Зависимость средней скорости вращения и синтеза от
концентрации АДФ в окружающей среде
На рисунке 5 представлены полученные нами реализации зависимости угла поворота γ субъединицы от времени при раз951
Раздел 8. Математические методы в биологии, экологии и химии
Part 8. Mathematical methods in biology, ecology and chemistry
личных величинах вращающего момента ( f ext есть внешний
вращающий момент на единицу угла поворота θ ). Замечу, что
угол наклона силового потенциала, образованного каталитиче∂Vs 20
ским субъединицами (рис. 2) равен
=
≈ 1.667 . Как и
∂θ 120
следовало ожидать, с ростом величины внешнего вращающего
момента скорость вращения растет, но появляются паузы, которые могут быть обусловлены ожиданием посадки АДФ в активный центр. При малых величинах вращающего момента время
перехода от одного центра к другому велико, и молекула АДФ
успевает связаться со следующим активным центром. С ростом
вращающего
момента
время
перехода
уменьшается
( τ step =
ξ
), и паузы связанные с ожиданием посад∂Vs ⎞
⎛
⎜ f ext − ∂θ ⎟ k BT
⎝
⎠
ки АДФ становятся различимы. На рисунке 6 представлена зависимость средней скорости вращения и синтеза от величины
вращающего момента ( f ext ). При больших величинах скорость
синтеза выходит на насыщение, скорость же вращения продолжает практически линейно расти. Это, на наш взгляд, объясняется либо тем, что синтезируемая молекула АТФ не успевает
покинуть АЦ и его «закрывает» γ субъединица, либо молекула
АДФ не успевает связаться с АЦ.
Рис. 5. зависимости угла поворота γ
субъединицы от времени при различных величинах вращающего момента
( f ext )
952
Рис. 6. зависимость средней скорости вращения и синтеза от величины вращающего момента ( f ext )
Погребная А. Ф. — МКО — 2005, ч. 3, стр. 944 – 953
Pogrebnaya A. F. — MCE — 2005, vol. 3, p. 944 – 953
Предложенная модель позволяет смоделировать синтез АТФ
во время вращения ротора F1-АТФазы по часовой стрелке под
действием внешней вращающей силы. Зависимость скорости
синтеза от концентрации АДФ подчиняется кинетике Михаэлиса-Ментена, константа Михаэлис-Ментена K m ≈ 187 мкМ при
используемых нами кинетических параметрах фермента (константы скоростей посадки/выхода субстратов, а также скорости
синтеза/гидролиза АТФ). Зависимость скорости синтеза от величины внешнего вращающего момента носит линейный харакΔV
тер вплоть до значений f ext ≈ 2 s . При больших значениях,
120
скорость вращения продолжает линейно расти, а скорость синтеза остается постоянной, что объясняется тем, что вновь синтезируемая молекула АТФ не успевает покинуть АЦ.
Автор благодарна за помощь в построении модели и всестороннее обсуждение результатов Воронову А.В., Романовскому
Ю.М. и Тихонову А.Н.
Работа выполнена при частичной финансовой поддержке программы «Научные школы России» (2071.2003.4), Российского фонда фундаментальных
исследований (грант РФФИ 05-03-32880), Междисциплинарного научного
проекта МГУ, а так же INTAS (01-0450).
Список литературы:
1. Pogrebnaya A.F., Romanovsky Y.M., Tikhonov A.N. Electrostatic interactions in catalytic centers of F1-ATPase // SPIE. 2003.
V.5068. P.27–35.
2. Тихонов А.Н., Погребная А.Ф., Романовский Ю.М. Электростатические взаимодействия в каталитических центрах F1ATPазы // Биофизика. 2003. Т.48. №6. С.1052–1070.
3. Pogrebnaya A.F., Romanovsky Yu.M., Tikhonov A.N., Stochastic model of molecular motor of ATPase // Complex Dynamics,
Fluctuations, Chaos and Fractals in Biomedical Photonics, Proceedings SPIE. 2004. V.5330. Р.120–131.
4. Walker J.E. (Ed.) The Mechanism of F1F0-ATPase // Biochim.
Biophys. Acta: special issue. 2000. V.1458. №2–3.
5. Senior, Nadinaciva S., Weber J. The molecular mechanism of
ATP synthesis by F1F0-ATP synthase // Biochim. Biophys. Acta.
2002. V.1553. P.188–211.
953
Раздел 8. Математические методы в биологии, экологии и химии
Part 8. Mathematical methods in biology, ecology and chemistry
6. Yasuda R., Noji H., Yoshida M., Kinosita K. Jr. & Itoh H. Resolution of distinct rotational substeps by submillisecond kinetic
analysis of F1-ATPase // Nature. 2001. V.410. №19. Р.898–904.
7. Itoh H., Takahashi A., Adachi K., Noji H., Yasuda R., Yoshida
M., Kinosita K. Mechanically driven ATP synthesis by F1ATPase // Nature. 2004. V.427.
SYNTHESIS OF ATP BY F1-AТPASE IN STOCHASTIC
MODEL
Pogrebnaya A. F.
(Russia, Moscow)
The paper continues our series of papers [1-4] devoted to the development of mathematical model on rotation of F1-AТPase molecular
motor. Here it has been considered the synthesis of ATP induced by
external force applied to the rotor.
954
Download