Расширение люмена крист митохондрий наблюдается при многих заболеваниях и при старении, сопровождается падением производства АТФ [1, 2]. Эти явления связаны с нарушением работы цепи митохондриального электронного транспорта на участке от III к IV дыхательному комплексу, где переносчиком электрона является единственный мобильный в люмене переносчик – цитохром с. Небольшой (около 12 кДа) белок, цитохром c (цитС), диффундирует в межмембранном и люминальном пространстве митохондрий и осуществляет перенос электрона от субъединицы цитохром c1 димера комплекса III к комплексу IV. Эту реакцию можно разбить на 2 стадии: связывание свободно диффундирующего белка цитС с цитС1 комплекса III в предварительный комплекс и электронный транспорт между цитохромами. В данной работе представлена модель броуновской динамики первой стадии этого процесса - образование предварительного (encounter) комплекса водорастворимого белка цитС с кофактором цитС1 комплекса III, фиксированого в мембране. На модели изучали влияние начального расположения цитС относительно комплекса III2 и роль электростатических взаимодействий в образовании предварительного комплекса.

С помощью программы ProKSim [3] нами были построены компьютерные модели диффузии, электростатического взаимодействия и связывания окисленного цитС (PDB ID 3O1Y) с восстановленным цит c1 III2 (PDB ID 1BGY) в растворе и люмене кристы. Для расчета электростатического потенциала вокруг белковых молекул использовали формализм Пуассона-Больцмана. Заряды на атомах белка были расставлены в соответствии с силовым полем CHARMM27, дополненным параметрами для гема и железо-серного кластера. Молекулы белков и электростатические потенциалы визуализировали с помощью PyMol (http://pymol.org/). Оценка параметров модели образования комплекса цитС1-цитС была предварительно выполнена для реакции в растворе. Значения параметров модели оценивали по экспериментальным данным [4], используя условие совпадения модельной и экспериментальной зависимостей константы связывания от ионной силы [5].

III и IV дыхательные комплексы могут образовывать суперкомплекс – респирасому [6], которая создает своеобразный «электростатический канал», направляющий движение окисленного цитС к комплексу III и движение восстановленного цитС к комплексу IV, обеспечивая быстрый перенос электрона между III и IV комплексами. Набухание крист при заболеваниях приводит к разрушению респирасом [7]. На модели мы исследовали влияние начальных положений и ориентации ЦитС относительно дыхательного комплекса III на скорость образования предварительного комплекса ЦитС1(III2)-ЦитС. Мембраны в модели учитывали как геометрические ограничения. Размеры реакционного объема: толщина люмена кристы (ТЛК) 12 и 16 нм, длина мембран 26 и 30 нм. Получены кинетические кривые образования предварительных комплексов при энергии электростатического взаимодействия белков -3.7kT и расстоянии между атомами Fe цитС и субъединицы P III2 менее 3.5 нм для ионной силы 130 мМ, pH 7. Кинетическая кривая связывания белков для каждого начального положения состояла из 20000 точек. Из кривых определяли время полупревращения образования предварительного комплекса белков t1/2. Было проведено по 60000 численных экспериментов диффузии цитС и образования предварительного комплекса с III2.

Исследовали 7 начальных положений цит С, у которых расстояние от атома Fe до Fe гема активного III мономера было 5.6 нм, отстоящих на 0.5 нм от поверхности мембраны. В одной серии численных экспериментов ориентация цитС была случайной относительно фиксированного положения центра масс цитС, в других была фиксирована и соответствовала средней структуре 10000 положений цитС, полученыых в предварительных расчетах диффузии цитС. Эти ориентации цитС может преобрети, если до связывания с cytC1 он долго диффундировал в электростатическом поле респирасомы.

Для 4 начальных положений цит С, лежащих ближе к IV комплексу в респирасоме 5GPN, получено t1/2 0.5-0.9 мкс, для остальных 3 положений t1/2 1.2-1.7·мкс. При увеличении толщины люмена кристы от 12 до 16 нм наблюдается увеличение вероятности нахождения цитС в области неактивной субъединицы цитС1 D III2 и увеличение времени t1/2 на 4-18%. Самое большее отличие t1/2 на 18% было получено для цитС предварительно не ориентированного в поле белков и из положения, лежащего ближе всего в респирасоме 5GPN к IV. Полученный эффект может способствовать уменьшению скорости транспорта электрона цитС при расширении крист и одновременном увеличении расстояния между дыхательными комплексами [7], когда цитС не может предварительно занть ориентацию в электростатическом поле респирасомы. Дальнейшее рассмотрение детальной модели, учитывающей структуру респирасом в кристах, позволит выявить их влияние на перенос электрона в дыхательной цепи.

Исследование выполнено в рамках научного проекта государственного задания МГУ №121032500060-0.



Литература

1. L. Colina-Tenorio et al., 2018 Shaping the mitochondrial inner membrane in health and disease, DOI: 10.1111/joim.13031

2. Siegmund S.E. et al., 2018, DOI: 10.1016/j.isci.2018.07.014

3. Хрущев С.С. и др. 2013, 47-64., DOI: 10.20537/2076-7633-2013-5-1-47-64

4. G. Engstrom, R. Rajagukguk, A. J. Saunders, C. N. Patel, S. Rajagukguk, T. Merbitz-Zahradnik, K. Xiao, G. J. Pielak, B. Trumpower, C. A. Yu, L. Yu, B. Durham, F. Millett,Biochemistry, vol.42, pp. 2816-24 (2003), DOI: 10.1021/bi027213g

5. Abaturova A.M. et. al. // ITM Web of Conferences, V. 31., 2020. DOI: 10.1051/itmconf/20203104001

6. Lapuente-Brun, E.; Moreno-Loshuertos, R.; Acin-Perez, R. et al. (2013).  Science 340 (6140): 1567–1570. DOI:10.1126/science.1230381

7. Sara Cogliati et al., 2013, DOI: 10.1016/j.сell.2013.08.032

Expansion of mitochondrial cristae lumen is observed in many diseases and aging, is accompanied by a decrease in ATP production [1, 2]. These phenomena are associated with disruption of the mitochondrial electron transport chain in the region from III to IV of the respiratory complex, where the electron carrier is the only mobile carrier in the lumen, cytochrome c. A small (about 12 kDa) protein, cytochrome c (cytC), diffuses in the intermembrane and luminal space of mitochondria and carries out electron transfer from the cytochrome c1 subunit of complex III dimer (III2) to complex IV. This reaction can be divided into 2 stages: binding of freely diffusible cytC protein to cytC1 of complex III, and electronic transport between cytochromes. This paper presents a model of the Brownian dynamics of the first stage of this process, i.e., the formation of an encounter complex of the water-soluble cytC protein with the cytC1 cofactor of complex III fixed in the membrane. The model was used to study the influence of the initial position of cytC relative to complex III2 and the role of electrostatic interactions in the formation of the encounter complex.

Using the ProKSim program [3], we constructed computer models of diffusion, electrostatic interaction, and binding of oxidized cytC (PDB ID 3O1Y) with reduced cytC1 III2 (PDB ID 1BGY) in solution and in the lumen of the crista. The Poisson-Boltzmann formalism was used to calculate the electrostatic potential around protein molecules. The charges on the protein atoms were arranged in accordance with the CHARMM27 force field, supplemented with parameters for the heme and iron-sulfur cluster. Protein molecules and electrostatic potentials were visualized using PyMol (http://pymol.org/). The estimation of the parameters of the model for the formation of the cytC1-cytC complex was preliminarily performed for the reaction in solution. The values of the model parameters were estimated from the experimental data [4] using the condition that the model and experimental dependences of the binding constant on the ionic strength coincide [5].

Respiratory complexes III and IV can form a supercomplex, a respirasome [6], which creates a kind of “electrostatic channel” that directs the movement of oxidized cytC to complex III and the movement of reduced cytC to complex IV, providing fast electron transfer between complexes III and IV. Swelling of cristae in diseases leads to the destruction of respirasomes [7]. Using the model, we studied the influence of the initial positions and orientation of CytC relative to the respiratory complex III on the rate of formation of the CytC1(III2)-CytC complex. The membranes in the model were taken into account as geometric constraints. Dimensions of the reaction volume: crista lumen width (CLW) 12 and 16 nm, membrane length 26 and 30 nm. Kinetic curves for the formation of preliminary complexes were obtained at an electrostatic interaction energy of proteins of -3.7kT and a distance between the Fe atoms of cytC and the cytC1 P subunit III2 of less than 3.5 nm for an ionic strength of 130 mM, pH 7. The kinetic curve of protein binding for each initial position consisted of 20000 points. From the curves, the half-life of the formation of the encounter complex of proteins t1/2 was determined. 60000 numerical experiments were carried out on the diffusion of cytC and the formation of a encounter complex with III2.

We studied 7 initial positions of CytC, in which the distance from the Fe atom to the Fe of the heme of the active monomer III was 5.6 nm, 0.5 nm away from the membrane surface. In one series of numerical experiments, the orientation of cytC was random with respect to a fixed position of the center of mass of cytC; in others, it was fixed and corresponded to the average structure of 10000 positions of cytC obtained in preliminary calculations of the diffusion of cytC. CytC acquires such orientations if it has been diffusing for a long time in the electrostatic field of the respirasome before binding with cytC1.

For 4 initial positions of cyt C, which are closer to complex IV in the 5GPN respirasome, t1/2 is 0.5–0.9 mks, for the remaining 3 positions t1/2 is 1.2–1.7 mks. With an increase of CLW from 12 to 16 нм, an increase in the probability of finding cytC in the region of the inactive cytC1 D III2 subunit and an increase in the time t1/2 by 4–18% is observed. The greatest difference in t1/2 was obtained for cytC not previously oriented in the field of proteins and from the position closest to IV in the 5GPN respirasome. The resulting effect may contribute to a decrease in the rate of CytC electron transport during the expansion of the cristae and a simultaneous increase in the distance between the respiratory complexes [7], when cytC can’t obtain a preliminary orientation in the resprasome electrostatic field. Further consideration of a detailed model that takes into account the conformation of the formed complexes and the structure of respirasomes in cristae will reveal their effect on electron transfer in the respiratory chain.

The study was carried out within the framework of the State Budget Project of Lomonosov Moscow State University No. 121032500060-0.



Literature

1. L. Colina-Tenorio et al., 2018 Shaping the mitochondrial inner membrane in health and disease, DOI: 10.1111/joim.13031

2. Siegmund S.E. et al., DOI: 10.1016/j.isci.2018.07.014

3. Khruschev S.S. et al., 2013, 47-64., DOI: 10.20537/2076-7633-2013-5-1-47-64

4. G. Engstrom, R. Rajagukguk, A. J. Saunders, C. N. Patel, S. Rajagukguk, T. Merbitz-Zahradnik, K. Xiao, G. J. Pielak, B. Trumpower, C. A. Yu, L. Yu, B. Durham, F. Millett,Biochemistry, vol.42, pp. 2816-24 (2003), DOI: 10.1021/bi027213g

5. Abaturova A.M. et. al. // ITM Web of Conferences, V. 31., 2020. DOI: 10.1051/itmconf/20203104001

6. Lapuente-Brun, E.; Moreno-Loshuertos, R.; Acin-Perez, R. et al. (2013). Science 340 (6140): 1567–1570. DOI:10.1126/science.1230381

7. Sara Cogliati et al., 2013 DOI: 10.1016/j.cell.2013.08.032