Next Article in Journal
Tumor Cell-Intrinsic BTLA Receptor Inhibits the Proliferation of Tumor Cells via ERK1/2
Next Article in Special Issue
Genetic Inhibition of Mitochondrial Permeability Transition Pore Exacerbates Ryanodine Receptor 2 Dysfunction in Arrhythmic Disease
Previous Article in Journal
Selective Menin Deletion in the Hippocampal CA1 Region Leads to Disruption of Contextual Memory in the MEN1 Conditional Knockout Mouse: Behavioral Restoration and Gain of Function following the Reintroduction of MEN1 Gene
Previous Article in Special Issue
Does Disruption of Optic Atrophy-1 (OPA1) Contribute to Cell Death in HL-1 Cardiomyocytes Subjected to Lethal Ischemia-Reperfusion Injury?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria

1
The Graduate Institute of Biomedical Electronics and Bioinformatics, National Taiwan University, Taipei 106, Taiwan
2
Department of Electrical Engineering, National Taiwan University, Taipei 106, Taiwan
*
Author to whom correspondence should be addressed.
Cells 2022, 11(24), 4020; https://doi.org/10.3390/cells11244020
Submission received: 15 September 2022 / Revised: 5 December 2022 / Accepted: 8 December 2022 / Published: 12 December 2022
(This article belongs to the Special Issue Mitochondrial Dysfunction in Cardiovascular Disease)

Abstract

:
Oxidative phosphorylation (OXPHOS) is an oxygen-dependent process that consumes catabolized nutrients to produce adenosine triphosphate (ATP) to drive energy-dependent biological processes such as excitation-contraction coupling in cardiomyocytes. In addition to in vivo and in vitro experiments, in silico models are valuable for investigating the underlying mechanisms of OXPHOS and predicting its consequences in both physiological and pathological conditions. Here, we compare several prominent kinetic models of OXPHOS in cardiomyocytes. We examine how their mathematical expressions were derived, how their parameters were obtained, the conditions of their experimental counterparts, and the predictions they generated. We aim to explore the general landscape of energy production mechanisms in cardiomyocytes for future in silico models.

1. Introduction

The metabolic process of oxidative phosphorylation (OXPHOS) consumes chemical energy from catabolism to generate adenosine triphosphate (ATP), the energy currency of the cell. OXPHOS complexes sit on the inner mitochondrial membrane (IMM), consisting of electron transport chain (ETC) complexes I, II, III, and IV, as well as ATP synthase (complex V) [1] (Figure 1). The ETC complexes transfer electrons from reducing equivalents (NADH, FADH2) to oxygen molecules in controlled steps and extract the released chemical energy to pump protons out of the mitochondrial matrix, generating a proton electrochemical gradient called the proton motive force (pmf) that drives ATP synthase to produce ATP [2,3,4]. Cells can couple ATP hydrolysis to energy-consuming biological processes, e.g., muscle contraction, ion homeostasis, and macromolecule synthesis. OXPHOS energy generation is more efficient than that of anaerobic glycolysis. Mitochondrial OXPHOS produces approximately 30 ATP molecules for one glucose molecule, compared to two ATP molecules in glycolysis [5]. Thus, mitochondria are famous as “the powerhouse of the cell [6]”.
Oxygen is an excellent electron acceptor, with a redox potential of ~800 mV per electron under physiological conditions [7]. Nonetheless, some oxygen molecules are not reduced completely to water in the ETC. Instead, they become reactive oxygen species (ROS): superoxide, hydrogen peroxide, and hydroxyl radicals [8]. These highly reactive molecules can serve as a redox signaling mechanism from the mitochondria under physiological conditions [9,10,11]. However, excessive ROS can damage biomolecules, such as mitochondrial DNA (mtDNA), mitochondrial proteins for metabolism (e.g., aconitase and complex I), and the peroxidation of membrane lipids crucial for OXPHOS function [12,13]. The ROS scavenging systems include superoxide dismutase (SOD), cytochrome c, catalase, glutathione (GSH), and thioredoxin (Trx) systems [14,15].
Mitochondrial dysfunction is involved in a variety of pathological conditions, such as type 2 diabetes [16,17,18,19,20,21], liver toxicity [22,23], neurodegenerative diseases [24,25,26,27,28], arrhythmia, and heart failure [29,30,31,32,33,34]. In this review, we focus on the kinetic modeling of ATP generation in cardiomyocyte mitochondria because (i) many computational models and experimental data on cardiomyocyte mitochondria are available; (ii) cardiomyocytes rely on mitochondria-derived ATP for action potential firing, muscle contraction, and ion homeostasis in each heartbeat [35]. Mitochondria can take up to 30% of the cardiomyocyte volume [36]. However, cardiomyocytes maintain metabolite stability under varying workloads [37]; (iii) the OXPHOS system is interconnected with other metabolic pathways, such as the citric acid cycle, mitochondrial calcium handling, metabolite transfer, and ROS generation and scavenging. The interactions between these systems are related to cardiomyocyte pathologies such as altered calcium dynamics and electrophysiology, heart failure, ischemic-reperfusion injury, diabetic cardiomyopathy, and chemotherapy-induced cardiomyopathy [38]. Investigating the complex interactions and inner workings of OXPHOS and related metabolic systems using systems biology approaches and computer simulation models can complement the knowledge gathered from in vitro and in vivo studies, identify candidate mechanisms, and corroborate/propose novel hypotheses [39,40,41,42]. For these reasons, we review the mathematical modeling of OXPHOS complexes, focusing on cardiomyocyte bioenergetics models.

2. Description of Enzyme Reaction Rates

2.1. Law of Mass Action

The simplest way to describe a biochemical reaction rate is to assume a sufficiently large number of particles spreading evenly in a solution. Since the reaction rate is proportional to the number of particles colliding with each other, it is also proportional to the product of substrate concentrations to the power of their stoichiometric coefficients. For an elemental reversible reaction aA + bB = cP + dQ, the reaction rate in the PQ-forming direction can be described as:
v = k f ( [ A ] a [ B ] b [ P ] c [ Q ] d / K eq )
where the forward rate constant ( k f ) is often an adjustable parameter determined by the literature and model fitting; the backward rate constant is constrained by the equilibrium constant ( K eq ). The reactants (A, B, P, and Q) can be substrates or different stages of the enzyme catalyzing the reaction.
The law of mass action is the basis of the famous Michaelis-Menten enzyme kinetics and more complex schemes.

2.2. Quasi-Steady State Enzyme Kinetics

Michaelis-Menten kinetic expressions are derived from the law of mass action in enzyme catalysis reactions. Take an irreversible reaction S → P catalyzed by enzyme E as an example, assuming that the first substrate binding reaction is reversible and the second product formation reaction is irreversible.
E + S E S E + P
Suppose we treat the intermediate enzyme-substrate complex (ES) as the fast species and use the conservation relationship (Et = ES + E = constant) to constrain the free enzyme (E). In that case, the quasi-steady-state assumption reduces the two reactions into one by time-scale separation. The product formation rate is then described as a Michaelis-Menten hyperbolic function of the concentration of substrate S [43,44,45,46,47], where V m a x is the maximal reaction rate when the enzyme is saturated, and K m is the Michaelis constant, the concentration of S when the reaction rate reaches half of V m a x .
V = V m a x [ S ] [ S ] + K m
For enzyme reactions involving more than one substrate/product, more combinations, such as a random binding order, a compulsory binding order, ping-pong mechanisms, and dead-end intermediates, are possible [47,48]. The number of parameters and the complexity of the mathematical description increase as more binding and conversion scenarios are considered [43,49]. For simplicity, the enzyme and enzyme-substrate complexes are often treated as fast species, and their transitions are assumed to be in a quasi-steady-state (the rates of change of the transition states are set to zero). The state occupancies and the enzyme turnover rate can be solved systematically by linear algebra or by the Hill–King–Altman diagram method [50,51].
For example, for a reversible enzyme-catalyzed reaction, S + E E S E P E + P the reaction scheme can be drawn as in Figure 2. kij are rate constants for the enzyme state i transitioning to state j. After applying quasi-steady state approximation for the enzyme stages, the state occupancies of each stage can be solved by linear algebra.
[ k 12 S k 13 P k 21   k 31 k 12 S k 21 k 23 k 32 k 13 P k 23 k 31 k 32 1 1 1 ] [ E E S E P ] = [ d E d t d E S d t d E P d t E ] = [ 0 0 0 1 ]
Alternatively, the state occupancies can be solved by the King–Altman–Hill diagram method [50]. This graph-based method first reduces the edges of the diagram to be a set of simple connected graphs without cycles. For example, the reaction cycle of Figure 2a can be broken down into three graphs: Λ, <, and > [50]. The weight of each stage ( w i ) is the sum of the products of the transition rates indicated in Figure 2b.
v = k 23 f 2 k 32 f 3 w 1 = k 23 k 31 + k 32 k 21 + k 31 k 21 w 2 = k 13 k 32 + k 23 k 21 + k 31 k 12 w 3 = k 13 k 23 + k 12 k 23 + k 21 * k 13 k 23 D = w 1 + w 2 + w 3 f 1 = w 1 / D f 2 = w 2 / D f 3 = w 3 / D k 12 = k 12 [ S ] k 13 = k 13 [ P ]
The quasi-steady state reaction flux is then the product of the turnover rate (v) and the enzyme activity. For enzymes with more states and state transitions, a computer program called KAPattern [51] can automatically derive the weights of each state and the overall reaction rate.

3. Lumped Models of the Electron Transport Chain

A simplistic approach to describe reaction rates in the ETC is to treat the ETC complexes (I, II, III, and IV) as a single-unit proton pump.
Korzeniewski’s OXPHOS models [52,53] assumed that complex I and III reactions were close to equilibrium. The overall reaction rate of NADH reducing cytochrome c linearly depended on the thermodynamic span ( Δ G 13 ), the deviation between the current condition and the equilibrium. Δ G 13 was determined by the ratios of NAD/NADH and oxidized/reduced cytochrome c couples, the mitochondrial membrane potential, and the proton motive force.
v C 13 = k redox Δ G 13 Δ G 13 = Δ G C 1 + Δ G C 3 Δ G C 1 = E m , Q E m , N 3 2 Δ p Δ G C 3 = E m , c E m , Q 2 Δ p + 2 Δ Ψ 2 E m , N = E m , N 0 + V T 2 ln NAD + [ NADH ] E m , c = E m , c 0 + V T ln c 3 + c 2 + V T = R T F
Magnus and Keizer [54] built a lumped complex I-III-IV ETC model (Figure 3) based on the six-state proton pump model by Pietrobon and Caplan [55] to study calcium dynamics and bioenergetics in mouse pancreatic beta-cell mitochondria. This model showed that the rates of proton translocation and oxygen consumption depended on the ratio of NADH to NAD and the mitochondrial membrane potential ( Δ Ψ m ).
J H R = 6 ρ r e s f H r a f r r a + r b f Ψ 1 + r 1 f r f B + r 2 + r 3 f r f Ψ J O = 0.5 ρ r e s f H r a + r C 1 f B f A + f Ψ r a + r c 2 f A 1 + r 1 f A f B + r 2 + r 3 f A f Ψ f r = [ NADH ] m / NAD + m f B = exp 6 Δ Ψ B / V T f Ψ = exp 6 g Δ Ψ m / V T f H = 10 6 Δ pH V T = R T / F
where J H R is the proton translocation rate, J O is the oxygen consumption rate, V T is the thermal voltage (~26.7 millivolts at 37 °C), ρ r e s is the concentration of the respiratory complex, Δ Ψ B is the bulk phase boundary potential, g is the correction factor for Δ Ψ m , and r a , r b , r c 1 , r c 2 , r 1 , r 2 , and r 3 are constant parameters for the ETC. This lumped model did not include the ROS generation mechanism or the redox states of ubiquinone and cytochrome c pools.
The mitochondrial energetics model by Cortassa et al. [56] adapted Magnus and Keizer’s model [54] and replaced the dependence on the mitochondrial membrane potential ( Δ Ψ m ) with the proton motive force (pmf, Δ p ). The model added succinate oxidation flux via complexes II, III, and IV, depending on the pmf and the FADH2 to FAD ratio.
V H e = 6 ρ res r a v A r a + r b v H 6 1 + r 1 v A v B 6 + r 2 + r 3 v A v H 6 V H e ( F ) = 4 ρ res ( F ) r a v F r a + r b v H 4 1 + r 1 v F v B 4 + r 2 + r 3 v F v H 4 V O 2 = 0.5 ρ res r a + r c 1 v B 6 + r c 2 v H 6 v A r a v H 6 1 + r 1 v A v B 6 + r 2 + r 3 v A v H 6 v A = K res [ N A D H ] / N A D + v F = K res ( F ) F A D H 2 / [ F A D ] v B = exp Δ Ψ B / V T v H = exp g Δ p / V T Δ p = Δ Ψ m 2.303 Δ pH V T = R T / F
where V H e is the proton translocation rate from complex I-III-IV, V H e ( F ) is the proton translocation rate from complex II-III-IV, V O 2 is the oxygen consumption rate, Δ p is the pmf, ρ r e s ( F ) is the concentration of complex II-III-IV, V T is the thermal voltage, ρ r e s is the concentration of the respiratory complex, Δ Ψ B is the bulk phase boundary potential, g is the correction factor for Δ Ψ m , and r a , r b , r c 1 , r c 2 , r 1 , r 2 , and r 3 are constant parameters for the ETC.

4. Complex I

Mitochondrial complex I, also called NADH: ubiquinone oxidoreductase or type I NADH dehydrogenase, is the largest complex (with 46 subunits in mammalian mitochondria) in the ETC [57,58,59]. Its primary function is to transfer two electrons from one reduced nicotinamide adenine dinucleotide (NADH) molecule to one ubiquinone (coenzyme Q, denoted as Q) and pump four protons from the mitochondrial matrix to the intermembrane space. The chemical reaction is summarized as follows:
NADH + Q + 5 H i n + = NAD + + QH 2 + 4 H o u t +

4.1. The Molecular Structure and Electron Transport in Complex I

The core molecular structure of complex I, sitting in the IMM, consists of one hydrophilic peripheral arm, protruding toward the mitochondrial matrix, and one hydrophobic membrane arm [58]. The peripheral arm conducts electron transfer from NADH to ubiquinone and consists of the NADH dehydrogenase module, the chain of iron-sulfur (FeS) clusters, and the ubiquinone reductase module. The membrane arm contains four parallel proton pumps; each translocates one proton once the bound ubiquinone is fully reduced and triggers a conformational change in the membrane arm. The redox process in complex I is reversible. Under a high pmf and a highly reduced Q pool with a large QH2 to Q ratio, the direction of the reaction favors QH2 oxidation, NAD+ reduction, and proton translocation from the IMS to the matrix. This is known as reverse electron transport (RET) [10]. Conversely, forward electron transport (FET) means that electrons flow ‘normally’ from complex I to ubiquinone.
The forward electron flow in complex I is NADH → FMN → FeS cluster chain (N3 → N1b → N4 → N5 → N6a → N6b → N2) → Q [57,58,59]. The NAD+/NADH couple can transfer only two electrons simultaneously, while the FeS cluster chain can handle only one electron at a time. Thus, FMN acts as a buffer between the NAD+/NADH couple and the FeS cluster chain because FMN can switch between three redox states: fully oxidized FMN, FMNH radical, and fully reduced FMNH2 [60]. Likewise, the ubiquinone molecule bound to the ubiquinone reductase receives electrons from the FeS clusters to form the ubisemiquinone radical (Q), and finally to the ubiquinol molecule (QH2). The rate of electron transport from NADH to FeS cluster N2 is much faster than the overall turnover rate of complex I (a time scale of 100 μs vs. 5 ms), with the FeS clusters mostly in the reduced state under physiological conditions, indicating that ubiquinone binding, reduction, and release, are rate-limiting steps [58]. The molecular structure of complex I also implies that the conformational changes induced by ubiquinone reduction couple the translocation of four protons from the matrix to the IMS [57].

4.2. Reactive Oxygen Species (ROS) Generation Mechanisms

Complex I is considered one of the major ROS producers in the mitochondria [12]. Several possible ROS-producing mechanisms are supported by experimental evidence in complex I [58,61,62,63]. Most studies recognize that the flavin site (IF) in the NADH dehydrogenase module and the quinone site (IQ) in the quinone reductase module are two major sources of ROS generation in complex I [64]. These two sites have different characteristics. ROS generation in the IF site is proportional to fully reduced flavin (FADH2) not bound by NAD/NADH [63]. The FeS cluster N1a can reduce the lifetime of superoxide-forming flavosemiquinone radical (FMNH•) [65]. ROS generation at the IF site is enhanced by blocking electron transport at the IQ site, for instance, by rotenone. ROS generation in the IQ site, on the other hand, is dependent on the fraction of bound semiquinone radicals, which is small during FET, but significant during RET. The inhibition of quinone reductase by rotenone decreases ROS generation during RET. Thus, whether rotenone enhances or inhibits ROS generation in complex I depends on several factors, including different substrates (glutamate/malate vs. succinate), the pmf, and the redox state of the ubiquinone pool, all of which determine the direction of complex I reaction. Kinetic modeling can help us elucidate the complicated dependence of the reaction rate on respiring substrates and the pmf in complex I.

4.3. Kinetic Models of Complex I

In Korzeniewski’s rate skeletal muscle OXPHOS models [53,66], the complex I reaction rate ( J C 1 ) depended hyperbolically on its thermodynamic span ( Δ E C 1 ), the deviation of current condition from thermodynamic equilibrium.
J C 1 = k C 1 Δ E C 1 Δ E C 1 + K m G 1
Later, in a rat skeletal muscle OXPHOS model [67] and a rat cardiomyocyte OXPHOS model [68] showed a linear relationship between the complex I reaction rate ( J C 1 ) and its thermodynamic span ( Δ E C 1 ).
J C 1 = k C 1 Δ E C 1
The thermodynamic span ( Δ E C 1 ) can be derived from substrate (NAD/NADH and ubiquinone/ubiquinol couples) concentrations and the pmf ( Δ p ).
Δ E C 1 = E m , Q E m , N 4 2 Δ p E m , N = E m , N 0 + V T 2 ln ( [ NAD + ] [ NADH ] ) E m , Q = E m , Q 0 + V T 2 ln ( [ Q ] [ QH 2 ] ) V T = R T F
Beard’s OXPHOS model [69] also assumed that complex I operated near-equilibrium under physiological conditions and used a simple law of mass action descriptor for complex I reaction flux ( J C 1 ).
J C 1 = X C 1 ( K e q C 1 [ NADH ] [ Q ] [ NAD + ] [ QH 2 ] ) K e q C 1 = [ H + ] in 5 [ H + ] out 4 exp ( Δ r G C 1 0 + 4 F Δ Ψ m R T )
where the enzyme activity ( X C 1 ) is an adjustable parameter fitted by experimental data of porcine cardiomyocyte mitochondria [70], and the equilibrium constant ( K e q C 1 ) depends on the proton concentration, mitochondrial membrane potential ( Δ Ψ m ), and Gibbs free energy change ( Δ r G C 1 0 ) of the reaction. The descriptors for ROS generation were not included.
Heiske et al. [71] built an OXPHOS model, which described the complex I reaction rate ( v C 1 ) based on Henry–Michaelis–Menten kinetics.
v C 1 = k f ϕ N A D H ϕ Q k b ϕ N A D ϕ Q H 2 ( 1 + ϕ N A D H + ϕ N A D ) ( 1 + ϕ Q + ϕ Q H 2 ) ϕ N A D H = [ NADH ] m K m N A D H ,     ϕ N A D = [ NAD ] m K m N A D ϕ Q = [ Q ] K m Q ,     ϕ Q H 2 = [ QH 2 ] K m Q H 2
where k f is the forward reaction constant, k b is the backward reaction constant, and K m x is the apparent Michaelis constant for species x (NADH, NAD, QH2, or Q). The reaction constants and apparent Michaelis constants are thermodynamically constrained by the change in the Gibbs free energy in the reaction.
Markevich and Hoek’s mitochondrial ROS model [63] also included a complex I model using the law of mass action for its 17 reactions. Reactions of NAD/NADH and Q/QH2 couples, the IF (flavin) [60] and IQ (Q-reduction) site, and the FeS clusters of complex I were included to identify the source of ROS formation. They concluded that the majority of ROS production in complex I came from the IF site during FET and the IQ site during RET.
Gauthier et al. [62,72] built a seven-state complex I model (Figure 4) based on the Magnus–Keizer [54,55] six state proton pump model. The enzyme went through a cycle of NADH oxidation, ubiquinone reduction, and proton translocation (state 1→ 2 → 3 → 4 → 7 → 5 → 6 →1). A superoxide-generating pathway (state 4 → 2) was included.
The complex I model by Bazil et al. [61] includes five states by the number of electron(s) the complex holds (Figure 5). They assumed that the electron transfer inside complex I was much faster than reactions with substrates; electrons were distributed in steady-state across the redox centers: FMN, FeS cluster N2, and the Q-site. To eliminate internal state variables, the state transitions were also assumed to be in a quasi-steady state. The enzyme turnover and ROS generation rates were functions of the pmf and substrate concentrations. The kinetic parameters were fitted by the kinetic data and the ROS data from bovine heart mitochondria and constrained by thermodynamics in redox chemistry. They reproduced the reaction rates between NAD/NADH and Q/QH2 couples and the ROS production rates of superoxide and hydrogen peroxide.

5. Complex II

5.1. Molecular Structure and Reaction Mechanism

Mitochondrial complex II is also called succinate dehydrogenase (SDH), or succinate-coenzyme Q reductase (SQR). Complex II couples the citric acid cycle (CAC) with the ETC by transferring electrons from succinate to ubiquinone (Q), forming fumarate and ubiquinol (QH2) [73]:
Succinate + Q = fumarate + QH2
The molecular structure of complex II consists of four subunits: SDHA (succinate dehydrogenase site, with FAD as the prosthetic group), SDHB (FeS clusters), SDHC and SDHD (ubiquinone reductase site) [74]. The electron transfer inside complex II is summarized as follows:
Succinate ↔ FAD ↔ FeS clusters ↔ Q
In the forward direction, succinate is oxidized to fumarate, and Q is reduced to QH2. Nevertheless, the reaction is reversible, depending on the concentrations of succinate, fumarate, Q, and QH2. Significant succinate accumulation is found in ischemic cells and the strong forward reaction upon reoxygenation is responsible for RET-related ROS generation and ischemic-reperfusion injury [9,10,75]. Complex II activity and ROS generation from RET are inhibited by oxaloacetate and malonate [76].

5.2. Kinetic Models of Complex II

In the guinea pig cardiomyocyte mitochondrial model by Cortassa et al. [56], SDH irreversibly oxidized succinate (SUC) into fumarate, and was inhibited by oxaloacetate (OAA) and fumarate (FUM). However, in this model, succinate oxidation was not coupled to proton pumping in the ETC.
V S D H = k c a t E T [ SUC ] [ SUC ] + K m ( 1 + [ OAA ] K i O A A ) ( 1 + [ FUM ] K i F U M )
where k c a t is the catalytic constant, E T is the SDH concentration, K m is the Michaelis constant for succinate, K i O A A is the inhibition constant for OAA, and K i F U M is the inhibition constant for fumarate.
Wei’s acid-base equilibria model [77] (also used in [14]) coupled ETC and CAC reactions by assuming that the FADH2 to FAD ratio in the Cortassa et al. model [56] was determined by the succinate to fumarate ratio. The reaction rate was also inhibited by OAA. However, because complexes II, III, and IV were lumped, there were no reactions of the ubiquinone/ubiquinol couple.
V H e ( F ) = 4 ρ res ( F ) r a v F r a + r b v H 1 + r 1 v F v B + r 2 + r 3 v F v H K i [ O A A ] + K i V O 2 S D H = 0.5 ρ res r a + r C 1 v B + r C 2 v H v A r a v H 1 + r 1 v A v B + r 2 + r 3 v A v H K i [ O A A ] + K i v F = K res SDH P SUC [ Succinate ] / [ Fumarate ] v B = exp 4 Δ Ψ B / V T v H = exp 4 g Δ p / V T V T = R T / F Δ p = Δ Ψ m 2.303 Δ pH
where V H e ( F ) is the proton translocation rate, V O 2 S D H is the oxygen consumption rate, Δ p is the pmf, ρ r e s ( F ) is the concentration of complex II-III-IV, K r e s S D H is the equilibrium constant of succinate oxidation, P S U C is the binding polynomial for unprotonated succinate, K i is the inhibition constant for OAA, Δ Ψ B is the bulk phase boundary potential, g is the correction factor for Δ Ψ m , and r a , r b , r c 1 , r c 2 , r 1 , r 2 , and r 3 are constant parameters for the ETC (complexes II, III, and IV).
Demin’s OXPHOS models [62,78,79] assumed succinate excess. The reaction rate of SDH depended hyperbolically on the fraction of oxidized ubiquinone in the Q pool.
J C 2 = V m a x f Q f Q + K m Q f Q = [ Q ] [ Q ] + [ QH 2 ]
The complex II formulation in the hepatocyte mitochondria model by Mogilevskaya et al. [80] used random order rapid equilibrium Bi-Bi enzyme kinetics. The model was used to study salicylate toxicity to hepatocytes by inhibiting CAC activity.
J S D H = [ SDH ] k f A B k b P Q 1 + A + α B + A B + P + β Q + P Q A = [ Suc ] K E s u c ,     B = [ Q ] K m Q ,     α = K m S U C K E s u c P = [ Fum ] K E f u m ,     Q = [ QH 2 ] K m Q H 2 ,     β = K m F U M K E f u m
The mitochondrial OXPHOS model by Wu et al. [81,82,83] used the Theorell–Chance Bi-Bi mechanism for SDH. Inhibition by oxaloacetate was also included.
V S D H = V f · a · b V b · p · q α i + a + α i b K m A K i A + α i p K m Q K iQ + q + a b + b q K m A K i A + a p K m Q K i Q + p q a = [ SUC ] K i A , b = [ Q ] K m B , p = [ QH 2 ] K m P , q = [ FUM ] K i Q α i = 1 + [ OAA ] K i O A A + [ SUC ] K a S U C + [ FUM ] K a F U M 1 + [ SUC ] K a S U C + [ FUM ] K a F U M
Manhas et al. [84] devised a five-state complex II model to study how substrates and inhibitors affect its reaction rate and ROS generation (Figure 6). The authors assumed that electron transport within complex II is much faster than substrate binding (succinate, fumarate, Q/QH2 couple, and oxygen). The five states represent how many electrons are held by complex II, ranging from zero to five; the state transitions represent reactions between complex II and its substrates. The thermodynamically constrained parameters were fitted against various datasets from bovine, pig, and guinea pig heart mitochondria experiments. The authors reproduced the result of succinate oxidation kinetics with and without inhibitors (atpenin and malonate) and ROS production rates from complex II. The authors concluded that the primary source of ROS is the [3Fe-4S] iron-sulfur complex under most physiological conditions.
Markevich et al. [85] developed a complex II model using mass action kinetics, which included substrate binding and electron transport processes with a total of 31 reactions. The parameters were fitted against experimental data from bovine and rat heart mitochondria. The authors also investigated the difference between intact and disintegrated complex II, where SDHA/SDHB subunits were separated from SDHC/SDHD subunits. The authors concluded that in disintegrated complex II, Q-binding site inhibition, and complex III inhibition had a high-amplitude bell-shaped dependence on succinate in ROS generation. In contrast, intact complex II under physiological conditions had a hyperbolic dependence on succinate in ROS generation.

6. Complex III

Respiratory complex III is also called the cytochrome bc1 complex, coenzyme Q:cytochrome c—oxidoreductase. It transfers electrons from ubiquinol (QH2) to cytochrome c (cytc) by the following reaction:
2 QH 2 ( p ) + Q ( n ) + 2 cytc 3 + + 2 H in + = 2 Q ( p ) + QH 2 ( n ) + 2 cytc 2 + + 4 H o u t +

6.1. Molecular Structure and Reaction Mechanism

Complex III translocates protons outward via Mitchell’s Q-cycle [86]. The first and rate-limiting reaction is ubiquinol oxidation to ubiquinone at the Qo (outer) site. The two electrons from ubiquinol bifurcate toward two paths. One electron goes through the high potential pathway (Qo → iron-sulfur protein (ISP) → cyt c1 → cyt c) to reduce cytochrome c; the other goes through the low potential pathway (Qo → cyt bL → cyt bH → Qi) to reduce ubiquinone at the Qi (inner) site. Two ubiquinol molecules oxidized at the Qo site complete a full cycle and release four protons to the intermembrane space. Meanwhile, one ubiquinone molecule at the Qi site is reduced to ubiquinol, absorbing two protons from the mitochondrial matrix. Overall, the Q-cycle translocates two charges for every ubiquinol oxidized at the Qo site with two electrons transferred [87].
On the inner mitochondrial membrane, complex III monomers form dimers, allowing electron transfer between their two cyt bL centers. However, electron transfer between monomers has been considered nonsignificant; the dimer acts as two independent monomers under physiological conditions [88,89,90].
Commonly used complex III inhibitors in experiments include antimycin A (AA) (blocking the Qi site), stigmatellin (completely blocking the Qo site), and myxothiazol (partially blocking the Qo site) [87,91,92].

6.2. Reactive Oxygen Species Generation Mechanisms

In addition to complex I, complex III is another major source of mitochondrial-derived reactive oxygen species. (mtROS) Semiquinone radicals (SQ) at the Qo site can reduce oxygen to form superoxide. There are two mechanisms for SQ formation, the semiforward and the semireverse mechanisms. The semiforward mechanism states that after ubiquinol (QH2) at the Qo site donates one electron to the iron-sulfur protein (ISP), the remaining SQ intermediate reduces oxygen to form superoxide. [89] In contrast, the semireverse mechanism argues that the SQ intermediate is too unstable to stay and play a significant role in superoxide generation. Instead, QH2 donates all its electrons downstream, one to ISP and the other to bL. [89] Superoxide forms when bL is highly reduced due to high pmf or the downstream pathway is blocked by inhibitors such as antimycin A. An electron can jump from reduced bL to a ubiquinone molecule at the Qo site, forming a transient SQ intermediate, which reduces oxygen to superoxide [87,93]. Both mechanisms of superoxide generation are strongly dependent on the pmf, but the semireverse mechanism implies that the rate of superoxide generation is maximized when the ubiquinone pool is moderately reduced, which is more consistent with experimental findings on antimycin-induced ROS generation from complex III [87]. Complex III releases superoxide to both sides of the IMM [94].

6.3. Kinetic Models of Complex III

In Korzeniewski’s rate skeletal muscle OXPHOS models [53,66], the complex III reaction rate ( J C 3 ) depended hyperbolically on its thermodynamic span ( Δ E C 3 ), the deviation of current condition from thermodynamic equilibrium.
J C 3 = k C 3 Δ E C 3 Δ E C 3 + K m G 3
A rat skeletal muscle OXPHOS model [67] and a rat cardiomyocyte OXPHOS model [68] by Korzeniewski et al. showed a linear relationship between the complex I reaction rate ( J C 3 ) and its thermodynamic span ( Δ E C 3 ).
J C 3 = k C 3 Δ E C 3
The thermodynamic span ( Δ E C 3 ) can be derived from substrate (ubiquinone/ubiquinol and oxidized/reduced cytochrome c couples) concentrations, the mitochondrial membrane potential, and the pH difference across the IMM (embedded in the pmf, Δ p ).
Δ G C 3 = E m , c E m , Q 2 Δ p + 2 Δ Ψ 2 E m , Q = E m , Q 0 + V T 2 ln ( [ Q ] [ QH 2 ] ) E m , c = E m , c 0 + V T ln ( [ c 3 + ] [ c 2 + ] ) V T = R T F
Beard’s OXPHOS model [69] assumed that the complex III reaction operated near equilibrium and used the law of mass action for the reaction rate. Activation from inorganic phosphate ( [ Pi ] m ) was added to explain the data in a mitochondrial experiment [70]. However, a later study [95] found that complex III activity was not enhanced directly by inorganic phosphate. The apparent activation by inorganic phosphate was likely due to associated substrate transport and nonsteady-state experimental conditions.
J C 3 = X C 3 1 + [ Pi ] m K p 1 1 + [ Pi ] m K p 2 ( [ c 2 + ] K e q C 3 [ QH 2 ] [ c 3 + ] [ Q ] ) K e q C 3 = [ H + ] in 2 [ H + ] out 4 exp ( Δ r G C 3 0 + 2 F Δ Ψ m R T )
Heiske et al. [71] built an OXPHOS model, where they described the complex III reaction rate ( J c 3 ) based on Michaelis–Menten kinetics. The apparent activation by inorganic phosphate was inherited from Beard’s OXPHOS model [69].
J C 3 = k f ϕ Q H 2 ϕ c 3 2 k b ϕ Q ϕ c 2 2 ( 1 + ϕ Q H 2 + ϕ Q ) ( 1 + ϕ c 2 + ϕ c 3 ) 2 [ Pi ] [ Pi ] + K A ϕ Q = [ Q ] K m Q ,     ϕ Q H 2 = [ QH 2 ] K m Q H 2 ϕ c 3 = [ c 3 + ] K m c 3 ,     ϕ c 2 = [ c 2 + ] K m c 2
where k f is the forward reaction constant, k b is the backward reaction constant, and K m x is the apparent Michaelis constant for species x (ubiquinone, ubiquinol, reduced and oxidized cytochrome c). The reaction constants and apparent Michaelis constants are thermodynamically constrained by the change in Gibbs free energy in the reaction.
Demin et al. [79] modeled Mitchell’s Q-cycle using the law of mass action. The system included ubiquinone diffusion on the IMM and electron transfer between ubiquinone, complex III redox centers (ISP, cyt c1, bL, and bH), and cyt c. The model used the semiforward mechanism for ROS generation; in the model, superoxide was generated from the oxygen oxidizing Qo site ubisemiquinone radical. The model consisted of two variants: the minimal and channeled models. The minimal model treated the Qo site ubisemiquinone radical as part of the Q pool, not binding to complex III. In the channeled model, on the other hand, Qo site ubisemiquinone radicals were bound to complex III, next to either ISP or heme bL. The superoxide generation rate showed exponential dependence on the mitochondrial membrane potential in the minimal model but sigmoid dependence in the channeled model. The minimal model was later integrated into the redox balance model by Gauthier et al. [62].
Demin et al. [78] also built a newer model by coupling the redox center states of ISP, cyt bL, and cyt bH to ensure electron bifurcation. The semiquinone radical at the Qo site remained bound until bL was oxidized by bH. The locking mechanism prevented both electrons from ubiquinol from going through ISP and bypassing the Q-cycle.
Selivanov’s complex III model [96] included all possible combinations of binding and redox states, including ubiquinone (at both the p-side and n-side), cytochromes bL and bH, ISP, and cytochrome c1, creating a 400-state ODE model. They found that there are two steady-state branches in complex III in simulations and experiments. Increasing the succinate supply or reducing the oxygen concentration switched complex III from a low ROS-generating mode to a high ROS-generating mode, which could explain the increased ROS production from ischemia-reperfusion.
Guillaud et al. [97] constructed a 16-state complex III model in which superoxide was produced via the semireverse mechanism by reduced cytochrome bL and oxidized ubiquinone. Mass action kinetics described the transition between the redox states of ISP, cyt bL, cyt bH, and the Qi site. The model reproduced the bell-shaped dependence of ROS generation on the fraction of reduced quinone in the experiments.
Markevich and Hoek [63] developed a complex III model using mass action kinetics to test competing ROS generation mechanisms under a variety of proton motive forces, substrate levels, and inhibitors. They concluded that the late dissociation of ISP, that is, ISP leaving cyt bL after ubiquinone, is more consistent with experimental data.
Bazil et al. [98] developed an algebraic descriptor of complex III to reproduce experimental results of electron transport, ROS generation, and bistability. The model used a six-state scheme to represent the number of electrons on the Qo and Qi sites. Kinetic rate constants were determined by experimental data, pH, the pmf, and changes in the Gibbs free energy for state transitions. The net turnover flux was derived from state occupancies at steady state, obtained by solving a linear system of equations. Later, a complex III functional dimer model [99] was built upon a similar principle (Figure 7). The dimer model had six-states representing zero to five electrons shared by the redox centers. The dimer model not only reproduced bistability in electron transfer and ROS generation but also showed that low-dose antimycin partially blocking the Qi site stimulates cytochrome c reduction, which can only be explained via a dimer model without using a different parameter set. The dimer model also revealed that ROS were generated predominantly via the semireverse mechanism by reduced cytochrome bL under most conditions. ROS generation from the semiforward mechanism (by semiquinone radical) became significant only with a highly reduced ubiquinone pool under antimycin blocking the Qi site.

7. Complex IV

Respiratory complex IV, also called cytochrome c oxidase (CcO), is the terminal enzyme of the ETC, receiving electrons from cytochrome c and catalyzing the reduction of oxygen to water [100]. For every oxygen molecule reduced, four cytochrome c molecules are oxidized, and eight positive charges are translocated outward against the mitochondrial membrane potential.
4 cytc 2 + + O 2 + 8 H in + 4 cytc 3 + + 2 H 2 O + 4 H o u t +

7.1. Molecular Structure and Reaction Mechanism

In complex IV, electrons flow from cytochrome c → CuA → cytochrome a → cytochrome a3-CuB binuclear center (BNC) to oxygen. Because every cytochrome c molecule carries only one electron, it requires four reduced cytochrome c molecules to fully reduce one oxygen molecule to water.
The active site of complex IV to oxygen is the BNC and the catalytic cycle consists of two phases: the eu-oxidase phase and the peroxidase phase [100]. In the eu-oxidase phase, an oxygen molecule binds doubly reduced BNC and is reduced to peroxide. In the peroxidase phase, BNC receives two electrons, fully reducing peroxide to water. Under physiological conditions, complex IV operates continuously and each electron transfer from cytochrome c is coupled to the outward translocation of one proton [100,101,102]. However, slip reactions (electron transfer without proton pumping) have been reported in isolated complex IV [103], with two protons instead of four being pumped per oxygen molecule reduced. It is also hypothesized that under a high pmf, complex IV trades proton pumping efficiency for a higher enzyme turnover via the slip pathway.
The BNC also evens out the redox potential changes for each electron transfer, keeping the midpoint potential for oxygen reduction at approximately −800 mV per electron [102]. The tight interactions between BNC and oxygen ligands also minimize ROS generation during the reaction [12,104].
Complex IV activity is inhibited by nitric oxide (NO) [105,106], cyanide (CN) [107], carbon monoxide (CO) [108], hydrogen sulfide (H2S) [109], azide (N3) [110], and formate [111], the toxicity of which can block the electron transport chain and ATP synthesis. Among complex IV inhibitors, nitric oxide (NO) is the most physiologically relevant because under physiological conditions, a substantial fraction (up to half) of complex IV activity is inhibited by NO [105,106]. NO inhibits complex IV in two modes: competitive inhibition with oxygen binding and uncompetitive inhibition in which NO is reduced to nitrite. NO binding to complex IV is also sensitive to near infrared and red photons, causing the dissociation of NO, and freeing up available complex IV for OXPHOS [106,112].

7.2. Kinetic Models of Complex IV

The early OXPHOS model by Korzeniewski and Froncisz [52] included a three-state (A1, A2, and A3) complex IV reaction cycle [113] that depended on the phosphorylation potential.
v 1 = k 1 A 2 a c 2 + v 1 = k 1 A 3 c 3 + v 2 = k 2 A 3 O 2 v 2 = k 2 A 1 a v 4 = k 4 a A 1 a + k 4 b A 1 b c 2 + A 2 a = A 2 1 1 + g g = c 3 + c 2 + 2 A A P K 5 A 1 a = A 1 1 1 + K 3 A 1 b = A 1 K 3 1 + K 3 K 3 = Q exp E m , a 3 E m , c u F R T E m , a 3 = E m 0 , a 3 R T F · ln 8 × 10 7 + A A P 5 × 10 3 + A A P A A P = [ ATP ] [ ADP ] [ Pi ]
The reaction rate depended on the ratio of oxidized and reduced cytochrome c couple, the oxygen concentration, and the phosphorylation potential term ( A A P ), which is the quotient of ATP by ADP and inorganic phosphate concentrations.
In Korzeniewski’s later OXPHOS models for hepatocytes [114] and myocytes [53], the descriptors for complex IV were simplified and dependent on the pmf. The reaction rate depended on the oxygen concentration, the reduced cytochrome c concentration ( c 2 + ), and the availability of reduced cytochrome a3 ( a 2 + ). a 2 + was assumed in rapid equilibrium and was derived from the ratio of oxidized and reduced cytochrome c couple, the mitochondrial membrane potential ( Δ Ψ ), the pH difference across the inner mitochondrial membrane (embedded in the pmf, Δ p ), and the midpoint potentials of cytochrome c ( E m c ) and cytochrome a3 ( E m a ) in the binuclear center.
v c 4 = k c 4 a 2 + c 2 + [ O 2 ] [ O 2 ] + K m , O 2 V T ln ( a 3 + a 2 + ) = Δ p ( 1 + u ) + E m c + V T ln ( c 3 + c 2 + ) E m a u = Δ Ψ Δ p V T = R T F 26 mV E m c = 250 m V E m a = 540 m V K m , O 2 = 12 μ M
Beard’s OXPHOS model [69] used mass action kinetics for cytochrome c and oxygen and included a hyperbolic dependence for the oxygen concentration.
J C 4 = X C 4 ( [ O 2 ] [ O 2 ] + K O 2 ) ( [ cytc 2 + ] [ cytc ] ) exp ( F Δ Ψ m R T ) ( ( K e q C 4 ) 0.5 [ cytc 2 + ] [ O 2 ] 0.25 [ cytc 3 + ] ) K e q C 4 = [ H + ] in 4 [ H + ] out 2 exp ( Δ r G C 4 0 + 4 F Δ Ψ m R T )
The reaction rate of complex IV ( J C 4 ) depended on the enzyme activity X C 4 , the oxygen concentration ( [ O 2 ] ) and affinity ( K O 2 ), the cytochrome c pool size ( [ cytc ] ), the mitochondrial membrane potential ( Δ Ψ m ), and the apparent equilibrium constant of complex IV ( K e q C 4 ), which depended on the proton concentration, the Δ Ψ m , and the changes in Gibbs free energy ( Δ r G C 4 0 ).
Heiske et al. [71] built an OXPHOS model that described the complex IV reaction rate ( v c 4 ) based on Henry–Michaelis–Menten kinetics.
v C 4 = k f ϕ c 2 ϕ O 2 0.25 k b ϕ c 3 ( 1 + ϕ c 2 + ϕ c 3 ) ( 1 + ϕ O 2 ) 0.25 ϕ c 3 = [ c 3 + ] K m c 3 ,     ϕ c 2 = [ c 2 + ] K m c 2 ϕ O 2 = [ O 2 ] K m O 2
where k f is the forward reaction constant, k b is the backward reaction constant, and K m x is the apparent Michaelis constant for species x (oxygen, reduced and oxidized cytochrome c). The reaction constants and apparent Michaelis constants are thermodynamically constrained by the change in the Gibbs free energy.
Demin et al. [78] built a reaction cycle of four states (Y, Yr, YO, and YOH) for complex IV with mass action kinetics as part of their OXPHOS model (Figure 8). The mitochondrial redox balance model by Gauthier et al. [62] used the King–Altman–Hill diagram method on this complex IV model to eliminate intermediate state variables of Y, Yr, YO, and YOH.
The complex IV model by Krab et al. [115] examined changing Michaelis constants for oxygen under different conditions. The 26-step model represented the redox and oxygen binding states of the redox centers. The steady-state oxygen consumption rate was determined by the King–Altman–Hill diagram method. The authors found that the apparent Michaelis constant for oxygen increased under a high pmf or a highly oxidized cytochrome c pool, consistent with experimental results.
Wilson et al. [116] modeled complex IV around the redox states of the binuclear iron-copper center. Between state transitions, the equilibrium constants came from the midpoint potentials of oxygen and binuclear center reduction reactions; the kinetic constants were fitted from experiments on isolated mitochondria. The steady-state oxygen consumption rate came from mass action kinetics. The authors showed that the reaction rates were highly dependent on the energy state Q, acidity, and redox state of the cytochrome c pool. They also revealed that the apparent Michaelis constant for oxygen increases as the energy state Q increases, similar to the results by Krab et al. [115].
Pannala et al. [117] improved the complex IV model by Wilson et al. [116], by converting the arbitrary energy state Q to the proton motive force (1 Q = 1.5 pmf), adding competitive and uncompetitive inhibition by nitric oxide (NO), adding a proton pumping mechanism coupled with each cyt c oxidation, assuming all steps are reversible to conform to strict thermodynamic constraints, and deriving steady-state reaction rates from the King–Altman–Hill diagram method. This model reproduced turnover rates under a range of pmf and cyt c redox states, as well as NO inhibition under a range of oxygen levels and cyt c redox states.

8. Complex V

OXPHOS complex V, also known as ATP synthase and F1Fo-ATPase, generates magnesium-bound ATP (MgATP) from magnesium-bound ADP (MgADP) and inorganic phosphate (Pi). To make this thermodynamically unfavorable reaction possible, complex V utilizes the proton electrochemical gradient across the inner mitochondrial membrane, the pmf through Peter Mitchell’s chemiosmotic mechanism [3,4,118,119,120].
MgADP + H 2 PO 4 + nH out + MgATP 2 + H 2 O + n H i n +
where n is the number of protons translocated per ATP, also called the H+/ATP ratio, which is approximately three in yeast and mammalian mitochondria [121,122].

8.1. Molecular Structure and Reaction Mechanism

Mitochondrial complex V is a membrane-bound molecular motor on the inner mitochondrial membrane (IMM). Complex V can be separated by water solubility into a soluble globular F1 catalytic sector in the matrix (subunit α, β, γ, δ, and ɛ) and a membrane-bound F0 proton-translocating sector (subunit c-ring, a, b, d, e, f, g, A6L, and oligomycin sensitivity-conferring protein (OSCP)). Mechanically, complex V can be divided into the rotor part (subunit c-ring, γ, δ, and ɛ) and the stator part (the remainder) [120,123,124].
Complex V acts like a turbine generator. The proton gradient across the inner mitochondrial membrane rotates the turbine (oligomeric c ring) and the stalk (gamma subunit) connected to the turbine [124]. In the active site of alpha3/beta3, MgADP and Pi are collected and assembled into MgATP. The mechanical interactions between the rotating stalk and the stator subunits α3/β3 induce a conformational change, altering the binding affinity of MgADP/Pi/MgATP substrates. MgATP is then “pushed” out of complex V, completing the catalytic cycle [120].
Complex V proteins often form a row of V-shaped dimers on the apex of mitochondrial cristae to enhance the efficiency of ATP synthesis. Because proton-pumping ETC complexes concentrate on mitochondrial cristae, the proton gradient is the greatest across the cristae interior and the intermembrane space, allowing efficient ATP synthesis [124,125,126].
Complex V can be inhibited by a range of substances. Oligomycin blocks the proton channels in the c-ring and is commonly used in mitochondrial metabolism studies [123].
ATPase inhibitory factor 1 (IF1) is a regulatory protein [123] that “locks” complex V activity during a low mitochondrial membrane potential or acidification of the matrix, preventing cellular ATP consumption in deenergized mitochondria [127,128,129].
Recent studies by Juhaszova et al. [130,131] demonstrated that complex V can also harness the “potassium motive force” to generate ATP in experiments on proteoliposome (PL)-reconstituted purified ATP synthase and isolated mitochondria. The potassium ion influx is driven by the mitochondrial membrane potential and the potassium concentration gradient.
MgADP + H 2 PO 4 + nK out + MgATP 2 + H 2 O + n K i n +
Potassium ions entering the mitochondrial matrix are translocated back out by the mitochondrial K+/H+ exchanger (mKHE), with a 1:1 stoichiometry of protons entering the mitochondrial matrix [132,133].
K in + + H out + H in + + K out +
Therefore, these study results are fully compatible with Mitchell’s chemiosmotic mechanism.

8.2. Kinetic Models of Complex V

In Korzeniewski’s rate skeletal muscle OXPHOS models [53,66,67], the ATP synthesis rate ( J C 5 ) was nonlinearly dependent on the thermodynamic force ( γ ) of complex V. The thermodynamic span was the displacement from the thermodynamic equilibrium and was determined by the pmf ( Δ p ), the ratio of ADP and inorganic phosphate to ATP, and the free energy change of ATP synthesis ( Δ G p 0 ). The H:ATP ratio ( n A ) was 8/3.
J C 5 = k C 5 γ 1 γ + 1 γ = exp ( n A Δ p F Δ G p 0 R T ) [ ADP ] m [ Pi ] m 1 M [ ATP ] m
Beard’s OXPHOS model [69] assumed that complex V operated close to equilibrium and utilized mass action kinetics for its substrates: ATP, ADP, and inorganic phosphate.
J C 5 = X F 1 ( K e q F 1 [ ADP ] m [ Pi ] m [ ATP ] m ) K e q F 1 = [ H + ] out 3 [ H + ] in 2 P A T P P A D P P P I exp ( Δ r G F 1 0 3 F Δ Ψ m R T )
where X F 1 is the activity of complex V, and K e q F 1 is the apparent equilibrium constant depending on the proton gradient, binding polynomials of substrates ( P A T P , P A D P , P P I ), change in Gibb’s free energy ( Δ r G F 1 0 ), and the mitochondrial membrane potential ( Δ Ψ m ).
Heiske et al. [71] built an OXPHOS model that described the complex V reaction rate ( J c 5 ) based on Henry–Michaelis–Menten kinetics.
J C 5 = k f ϕ D ϕ P k b ϕ T ( 1 + ϕ D + ϕ T ) ( 1 + ϕ P ) ϕ D = [ ADP ] m / K A D P ϕ T = [ ATP ] m / K A T P ϕ P = [ Pi ] m / K P i
where k f is the forward reaction constant, k b is the backwards reaction constant, and K m x is the apparent Michaelis constant for species x (ATP, ADP, inorganic phosphate). The reaction constants and apparent Michaelis constants are thermodynamically constrained by the change in the Gibbs free energy.
The beta cell mitochondria model by Magnus and Keizer [54] assumed that complex V was a six-state proton pump operating in reverse, using the same framework as the ETC in the same model. The ATP synthesis rate and the proton consumption rate strongly depended on both the mitochondrial membrane potential and the phosphorylation potential. The H+/ATP stoichiometry was higher than 3:1 due to the presence of slip reactions.
V ATPase = ρ F 1 10 3 Δ p H p a + p c 1 f B + p c 2 f m f A p a f m 1 + p 1 f A f B + p 2 + p 3 f A f m V H u = 3 ρ F 1 10 3 Δ p H p a f A p a + p b f m 1 + p 1 f A f B + p 2 + p 3 f A f m f A = K e q F 1 [ A T P ] m [ P i ] m [ A D P ] m f B = exp 3 Δ Ψ B / V T f m = exp 3 Δ Ψ m / V T Δ p H = p H i p H m V T = R T F
where V A T P a s e is the ATP generation rate, V H u is the proton flux through complex V, V T is the thermal voltage, ρ F 1 is the concentration of complex V, K e q F 1 is the apparent equilibrium constant, Δ Ψ B is the phase boundary potential, and p a , p b , p c 1 , p c 2 , p 1 , p 2 , and p 3 are the parameters of complex V.
The mitochondrial energetics model for guinea pig ventricular cardiomyocytes by Cortassa et al. [56] replaced the dependence on the mitochondrial membrane potential with the dependence on the pmf in Magnus and Keizer’s model (Figure 9). A factor of 100 was also added to address reversibility under low pmf. Later models [134,135,136,137,138,139] also utilized this schematic for mitochondrial energetics in guinea pig ventricular cardiomyocytes.
V ATPase = ρ F 1 100 p a + p c 1 v B + p c 2 v H v A p a v H 1 + p 1 v A v B + p 2 + p 3 v A v H V H u = 3 ρ F 1 100 p a 1 + v A p a + p b v H 1 + p 1 v A v B + p 2 + p 3 v A v H v A = K e q F 1 [ A T P ] m [ P i ] m [ A D P ] m v B = exp 3 Δ Ψ B / V T v H = exp 3 Δ p / V T Δ p = Δ Ψ m 2.303 V T · Δ p H Δ p H = p H i p H m V T = R T / F
The complex V model by Cortassa et al. in 2022 [140] incorporated ATP generation from the potassium gradient, based on the studies of Juhaszova et al. [130,131]. The mathematical descriptions were based on the model of Magnus and Keizer [54] and Wei et al. [77]. In addition to the proton flux ( J C 5 H ), the model also considered the potassium ( J C 5 K ) and the sodium fluxes ( J C 5 N a ) through complex V. All three fluxes contributed to ATP synthesis ( J C 5 A T P ).
J C 5 A T P = ρ F 1 p a f M C a + p c 1 v B v A p a + p c 2 v A v I / Δ J C 5 H = 3 ρ F 1 10 3 Δ pH p a v A p a + p b v H / Δ J C 5 K = 3 ρ F 1 f K 10 3 Δ pK p a v A p a + p b v K / Δ J C 5 N a = 3 ρ F 1 f N a 10 3 Δ pNa p a v A p a + p b v N a / Δ Δ = 1 + p 1 v A v B + p 2 + p 3 v A v I v B = exp 3 Δ Ψ B / V T v H = exp 3 Δ μ H / V T v K = exp 3 Δ μ K / V T v N a = exp 3 Δ μ N a / V T v I = exp 3 Δ μ H + Δ μ K + Δ μ N a / V T = v H v K v N a v A = K q F 1 MgATP 2 m Pi m [ ADP ] m free f M C a = 10 3 Δ pH + f K 10 3 Δ pK + f N a 10 3 Δ pNa Δ μ H = Δ Ψ m + Z Δ pH Δ μ K = Δ Ψ m + Z Δ pK Δ μ N a = Δ Ψ m + Z Δ pNa Z = 2.303 V T V T = R T / F Δ pH = pH i pH m Δ pK = pK i pK m Δ pNa = pNa i pNa m
The complex V descriptors in Nguyen et al.’s cardiomyocyte OXPHOS model [141] used Hill kinetics for the mitochondrial membrane potential and random-order Bi-Uni kinetics for ADP, ATP, and inorganic phosphate (Pi). The calcium activation of ATP synthesis was also considered. The complex V reaction was assumed to be fully coupled; thus, the H:ATP ratio was exactly 3:1.
V C 5 = V m a x ( f D f P f T ) f v f c 1 + f D + f P + f T + f D f P f D = [ A D P ] / K A D P f T = [ A T P ] / K A T P f P = [ P i ] / K P i f v = Δ Ψ m 8 K v 8 + Δ Ψ m 8 f c = 1 exp ( [ Ca 2 + ] m / K C a )
where V C 5 is the ATP synthesis rate, V m a x is the activity of complex V, K A T P , K A D P , and K P i are the dissociation constants for ATP, ADP, and inorganic phosphate, respectively; K v is the voltage activation factor, and K C a is the calcium activation factor.
The complex V model by Metelkin et al. [79,142] to study the interactions between complex V and the adenine nucleotide translocator (ANT) used a similar scheme and included the concentrations of protons on both sides of the IMM.
V C 5 = V max C 5 H o K H o exp ( χ ϕ ) n N D e n N = f D f P f T K e q H o H i exp ( ϕ ) n Den = 1 + f D f P f H o + f T f H i exp ( n ( 1 χ ) ϕ ) K e q = K e q C 5 K M g T C 5 K M g A T P K M g D C 5 K P i C 5 K M g A D P · 10 7 M 10 7 M + K P , H f D = [ M g A D P ] / K M g D C 5 f T = [ M g A T P ] / K M g T C 5 f P = [ P i ] / K P i C 5 f H o = H o K H o C 5 n f H o = H i K H i C 5 n ϕ = F Δ Ψ m R T

9. Summaries of the OXPHOS Models

In 1967 Chance proposed the first OXPHOS model using operational flux expressions [143], Bohnensack [144] and Holzhutter [145] later constructed the first kinetic OXPHOS model and used thermodynamic force to drive the OXPHOS reactions [146]. These studies have established the foundation of kinetic modeling of OXPHOS, and more detailed and sophisticated models have been proposed since then [146]. In this section, we summarize the mitochondrial OXPHOS computational modeling mentioned in the previous sections in Table 1. Because many of the studies included two or more respiratory complexes, we provide overall remarks in this section rather than in the respective OXPHOS complex sections above.
The rat hepatocyte [52,114] and skeletal muscle cell OXPHOS models [53,67] by Korzeniewski et al. included OXPHOS complexes I, III, IV, and V, along with proton leak, substrate dehydrogenation, adenine nucleotide translocator (ANT), phosphate carrier (PiC), and adenylate kinase (AdK). Cytochrome c reduction (by complex I and complex III) and ATP synthesis (by complex V) rates were linked to their thermodynamic spans, which is the displacement from thermodynamic equilibria. Complex I, III, and V reactions were assumed to operate near equilibrium [147]. Cytochrome c oxidation (by complex IV) used mass action kinetics. The oxygen-binding reduced heme a3 was assumed to be in rapid equilibrium and derived by oxidized to reduced cytochrome c ratio and the pmf. Korzeniewski et al. used these models to study the metabolic flux control of ATP generation and consumption in different metabolic modes (glycolysis, fatty acid oxidation, and gluconeogenesis) in hepatocytes [114,148], metabolic flux control and the impact of respiratory complex deficiencies in skeletal muscle cells [53], comparing the negative feedback and parallel activation mechanisms of ATP synthesis and consumption in skeletal muscle cells [66,67,149], and the regulation of OXPHOS reactions in cardiomyocytes [68]. They demonstrated that parallel activation to NADH generation, proton pumping by respiratory complexes, ATP synthesis, and ATP consumption could explain why the ATP to ADP ratio barely changes when the work intensity increases. The OXPHOS models by Korzeniewski et al. are simple, thermodynamically consistent (for complexes I, III, V, and ANT), and in good agreement with experimental data of ATP synthesis, proton leak, and oxygen consumption. However, they did not include detailed descriptions of the citric acid cycle, calcium dynamics, and ROS generation/scavenging mechanisms.
Saito et al. [150] built a cardiomyocyte mitochondrial model, which integrated complex I, III, IV, and V models by Korzeniewski et al. [52,53,67,114], citric acid cycle reactions, ion (proton, sodium, calcium, potassium, and phosphate) and respiring substrate (malate and glutamate) dynamics, and high energy phosphate buffering by adenylate kinase and creatine kinase. The authors used this model to study why the metabolite levels of NADH and ATP remain relatively stable under various workloads. They found that inorganic phosphate (Pi) stimulated ATP synthesis more effectively than calcium.
Beard’s OXPHOS model [69] included complexes I, III, IV, and V, as well as ANT, PiC, AdK, potassium leakage, and potassium efflux through mKHE. The reaction rates of the OXPHOS complexes were described by mass action kinetics, assuming they operated near equilibrium. The reaction rate constants were also constrained by the change of Gibbs free energy in each complex and fitted by the Bose et al. dataset [70] of porcine heart and skeletal muscle mitochondria. However, this mitochondrial model did not include explicit citric acid cycle processes, ROS generation/scavenging, or mitochondrial calcium dynamics. This model required complex III to be explicitly activated by inorganic phosphate to fit the experimental data. The simple mass action kinetics descriptors did not follow saturation kinetics in enzyme-catalyzed reactions.
Wu et al. [82] extended Beard’s OXPHOS model [69] and integrated citric acid cycle reactions (including complex II, succinate dehydrogenase) and substrate exchange across the IMM. The parameters were estimated based on the Bose et al. dataset [70] of porcine heart and skeletal muscle mitochondria and the LaNoue et al. dataset [151] of rat heart mitochondria. The authors concluded that ADP and NAD regulated the citric acid cycle fluxes. In contrast, ATP synthesis was regulated by inorganic phosphate, which activated complexes III and V. However, ROS generation/scavenging and calcium dynamics were not included in this model. Bazil et al. [83] extended the mitochondrial model by Wu et al. [82] and integrated sodium, calcium, and potassium influx/efflux reactions as well as mitochondrial matrix volume dynamics due to osmolarity changes caused by the matrix potassium concentration.
Heiske et al. [71] built an OXPHOS model based on the OXPHOS models of Beard [69] and Korzeniewski and Zoladz [67]. The authors used Henry–Michaelis–Menten kinetics instead of mass action kinetics for the rate equations of the OXPHOS complexes (I, III, IV, V) to address saturation in enzyme kinetics. The parameters were constrained by the changes in the Gibbs free energy and fitted by their experimental data [71,152] and the Bose et al. dataset [70]. However, this model did not include citric acid cycle reactions, ROS generation and scavenging, and calcium dynamics.
The complex III monomer [98] and dimer [99] model by Bazil et al., the complex I model [61] by Bazil et al., and the complex II model by Manhas et al. [84] have used quasi-steady state assumptions to describe the overall reaction rates of substrate consumption and ROS formation. Each model had several redox states representing the number of electrons in the redox centers. Redox reactions between the complex and substrates were state transitions. The kinetic parameters were constrained by Gibbs free energy changes of the reactions and fitted by various datasets, primarily from bovine cardiomyocyte mitochondria experiments. By applying the quasi-steady state assumption and setting the rate of change of each state to zero, the overall reaction rates could be solved from a system of linear equations of states and their transition rates. These respiratory complex models have reproduced redox reaction rates under various conditions. They have been used to study the source of ROS generation in ischemia/reperfusion (I/R) injury [81,153] and to test the hypothesis of calcium phosphate inhibiting complex I during calcium overload in a computational-experimental approach [154].
Markevich et al. built complex I [63], III [63], and II [85] models using mass action kinetics. These models have been used to study ROS generation mechanisms and steady-state ROS production rates in various substrate inputs and a range of mitochondrial membrane potentials. Because their mass action kinetics models included detailed reaction steps of substrate binding, proton pumping, enzyme state transitions, and ROS generation, the authors can test different ROS generation mechanisms and simulate complex inhibition by changing the kinetic constants.
The minimal mitochondrial model by Magnus and Keizer [54] included a lumped respiratory chain (complex I-III-IV) and ATP synthase (complex V) for the OXPHOS system as well as adenine ANT, calcium influx via the mitochondrial calcium uniporter (MCU), and calcium efflux via the mitochondrial sodium-calcium exchanger (mNCX). Matrix proton, inorganic phosphate, and sodium concentrations were assumed to be constant, so their dynamics were not included. ROS production and scavenging systems were also absent. The OXPHOS system used the six-state proton pump scheme from Pietrobon and Caplan [55]. The (quasi-)steady-state rates of NADH oxidation, proton translocation, and ATP synthesis were derived by the King–Altman–Hill diagram method [50]. Magnus and Keizer later incorporated this mitochondrial model with cellular glycolysis and ion channels to simulate the mouse beta-cell membrane potential, ATP, and calcium oscillations in the presence of glucose [155,156].
The mitochondrial energetics model by Cortassa et al. [56] integrated the calcium dynamics and OXPHOS from the Magnus–Keizer model [54] and the citric acid cycle (CAC) reactions from the Jafri–Dudycha model [157]. The CAC provided NADH for the ETC and was stimulated by mitochondrial calcium. The model was used to simulate the generation and consumption of NADH and ATP in cardiomyocyte mitochondria under various workloads, demonstrating that mitochondrial calcium can enhance energy supply. However, the succinate oxidation in the respiratory chain did not couple with the succinate dehydrogenase reaction in the CAC. Later the same group integrated the mitochondrial model with electrophysiology and muscle contraction systems and constructed the excitation-contraction coupling/mitochondrial energetics (ECME) model of guinea pig ventricular cardiomyocytes [136]. The simulations were in good agreement with the experimental results of force-frequency relationships, mitochondrial NADH, and mitochondrial calcium changes in response to varying workloads.
Bertram et al. [158] simplified Cortassa et al.’s model [56], known as the BLPS model, to simulate mitochondrial ATP synthesis. The electron transport chain and ATP synthase were described by Michaelis–Menten kinetics for the substrates (NADH, ADP) and sigmoid dependence on the mitochondrial membrane potential. Saa and Siqueira [159] improved and corrected the BLPS model to study the metabolic response under a range of glycolysis fluxes and oscillatory calcium levels.
Nguyen et al. [141] built a cardiomyocyte mitochondrial model that integrated the ETC model from Magnus and Keizer [54], the CAC reactions from the Jafri–Dudycha model [157], and the calcium uniporter and ATP synthase model from the authors’ previous model [160]. The authors added NCX, NHX, and PiC to simulate inorganic phosphate, sodium, and proton dynamics. The models can reproduce mitochondrial calcium concentrations and influx under periodic stimulation, consistent with cardiomyocyte experiments. They also showed that increasing cytosolic sodium levels enhanced calcium extrusion from the mitochondria via the NCX, leading to decreased mitochondrial calcium concentrations and ATP synthesis.
Cortassa et al. [135] built a ROS-induced ROS release (RIRR) model based on their previous model [56]. The authors integrated CAC, OXPHOS, ROS production, ROS-inducible anion channels for ROS efflux, and ROS scavenging by catalase, SOD, and the glutathione (GSH) system. Since there were no mechanistic descriptions for ROS production, it was modeled by taking an adjustable fraction of oxygen consumption to produce superoxide. They simulated the ROS-induced ROS release (RIRR) phenomenon and reproduced mitochondrial membrane potential depolarization and recovery events observed in guinea pig ventricular cardiomyocytes under oxidative stress [161].
Zhou et al. [138] incorporated the ECME model [136] and the RIRR model [135] to study action potential shortening and insensitivity to electrical stimuli due to oxidative stress-induced mitochondrial depolarization, ATP synthase reversal, decreased ATP to ADP ratio, and the opening of ATP-sensitive potassium channels. The ECME-RIRR model has been used to study potential oscillations in a mitochondrial network [137,162], arrhythmogenesis from blocked action potential propagation [139], and abnormal action potential patterns due to ROS disrupting calcium dynamics in the cytosol [163,164].
Wei et al. [77] extended the mitochondrial energetics model by Cortassa et al. [56] by integrating proton, phosphate, sodium, and calcium dynamics. The authors included mathematical descriptors for PiC, MCU, and mNCX. The authors also improved the complex II model by coupling succinate oxidation in the ETC and the CAC. The model results were in good agreement with mitochondrial membrane potential, mitochondrial NADH, and pH changes in guinea pig cardiomyocyte experiments. However, ROS production and scavenging mechanisms were not included in this model.
Kembro et al. [14] extended the mitochondrial model of Wei et al. [77] and the RIRR model of Cortassa et al. [135] to build a two-compartment (cytosol and mitochondria) ROS production, transport, and scavenging model. The mitochondrial energetic-redox (ME-R) model included mitochondrial NADPH generation from transhydrogenase (THD) and NADPH-producing isocitrate dehydrogenase (IDH2), ROS efflux via anion channels (for superoxide) and simple diffusion (for hydrogen peroxide), and ROS scavenging by catalase, SOD, and the GSH and Trx systems [15]. The fatty acid oxidation model by Cortassa et al. [165] added beta oxidation steps for palmitate to the two-compartment ROS generation, transfer, and scavenging model by Kembro et al. [14] The experimental-computational approach showed that ROS excess upon palmitate addition resulted from impaired GSH and Trx ROS scavenging systems and palmitate-induced uncoupling.
Cortassa et al. [140] extended the model of Kembro et al. [14] by including a complex V model utilizing the potassium gradient to generate ATP and an ETC model coupled to mitochondrial matrix volume dynamics, which was controlled by potassium content. The authors reproduced experimental data of changes in matrix volume, potassium uptake, and mitochondrial membrane potential from low to high energy demand.
Gauthier et al. [62] built a mitochondrial model, called the ETC-ROS model, to study the redox balance hypothesis [166], which stated that ROS production is minimal and energy generation is maximal when the cellular environment is moderately reduced. Thus, a mechanism of ROS generation from the respiratory complexes was needed. The complex I model was derived from the six-state proton pump scheme of Pietrobon and Caplan [55]. The complex I model included an oxygen-reducing pathway for superoxide production and an additional state for ubiquinone reduction. It reproduced the guinea pig cardiomyocyte mitochondrial data of Aon et al. [166]. The complex II, III, and IV models were modified from the hepatocyte mitochondrial model by Demin et al. [78,79] to reproduce the macrophage OXPHOS data from Kim et al. [167]. The complex V model was based on the mitochondrial model by Wei et al. [77]. The ROS scavenging systems included the SOD, GSH, and Trx systems. The authors showed that mitochondrial ROS production was minimized when the redox environment (NAD/NADH and NADP/NADPH couples) was moderately reduced, supporting the redox-optimized ROS balance (R-ORB) hypothesis [166].
Later Gauthier et al. reported a new mitochondria model [72] that incorporated the OXPHOS and ROS generation parts from the ETC-ROS model [62] and the ROS scavenging parts from the ME-R model [14] to study ROS production in heart failure. Their model can reproduce respiration rates, NADH to NAD ratios, and ROS production rates at state three (ADP present) and state four (ADP absent). The model also demonstrated that either blocking mitochondrial calcium uniporter or increasing cytosolic sodium concentrations (as in heart failure) suppressed mitochondrial calcium, NADH, NADPH, and GSH levels while ROS production increased. De Oliveira et al. [168] extended this model [72] and included doxorubicin-induced acute and chronic cardiac toxicity. Their model can simulate acute toxicity by inactivating OXPHOS complexes and enhancing ROS production. In contrast, chronic toxicity was simulated by ROS damaging mitochondrial DNA (mtDNA), inhibiting the long-term activity of OXPHOS complexes and leading to a vicious cycle as faulty OXPHOS complexes tend to produce more ROS and damage mtDNA even further. The authors predicted the dose response of both acute and chronic doxorubicin toxicity. They also showed that extending the duration of the iron chelator therapy partially mitigated chronic doxorubicin toxicity and protected cardiomyocytes against mitochondrial dysfunction.
Table 1. List of published OXPHOS kinetic models reviewed in this paper.
Table 1. List of published OXPHOS kinetic models reviewed in this paper.
ModelRespiratory ComplexCell TypeMethodInternal State Variables
Korzeniewski, 1991 [52]I-III
IV
V
Rat hepatocytesTS
MAK
TS
0
2
0
Korzeniewski, 1996 [114]IVRat hepatocytesMAK0
Korzeniewski, 1996 [53]I, IIIRat skeletal muscle cellsTS0
Magnus and Keizer, 1997 [54]I-III-IV
V
Mouse pancreatic beta-cellsKAH
KAH
0
0
Cortassa, 2003 [56]I-III-IV
II-III-IV
V
Guinea pig ventricular cardiomyocytesKAH
KAH
KAH
0
0
0
Saa, 2013 [159]I-III-IV
V
Not specifiedPhenomenal
Phenomenal
0
0
Beard, 2005 [69]I, III, IV, VPorcine heart and skeletal muscle cellsMAK0
Wu, 2007 [82]IIPorcine heart and skeletal muscle cells and rat cardiomyocytesTheorell–Chance Bi-Bi mechanism0
Heiske, 2017 [71]I, III, IV, VBovine cardiomyocytesMM0
Markevich and Hoek, 2015 [63]I
III
Bovine cardiomyocytesMAK12
14
Gauthier, 2013 [62]I
III
Guinea pig ventricular cardiomyocytesKAH
MAK
0
10
Bazil et al., 2014 [61]IBovine cardiomyocytes and rat cardiomyocytesHKA0
Manhas et al., 2020 [84]IIBovine/Pig/Guinea pig cardiomyocytesKAH0
Markevich et al., 2020 [85]IIBovine and rat cardiomyocytesMAK35
Demin et al., 1998 [79]IIINot specifiedMAK9 (minimal)/12 (channeled)
Demin et al., 2001 [78]III
IV
HepatocyteMAK
MAK
12
4
Selivanov et al., 2009 [96]IIIRat brain mitochondriaMAK400
Guillaud et al., 2014 [97]IIIRat skeletal muscle, heart, liver, kidney, and brain cellsMAK16
Bazil et al., 2013 [98]IIIBovine cardiomyocytesKAH0
Krab et al., 2011 [115]IVBovine cardiomyocytesKAH0
Wilson et al., 2014 [169]IVRat hepatocytesQuasi-steady state from MAK0
Pannala et al., 2016 [117]IVRat hepatocytes and bovine cardiomyocytesKAH0
Nguyen et al., 2007 [141]VGuinea pig ventricular cardiomyocytesHill kinetics for MMP,
Random Bi-Uni kinetics for substrate
0
Cortassa et al., 2022 [140]VGuinea pig ventricular cardiomyocytesKAH0
MAK: mass action kinetics. KAH: King–Altman–Hill diagram method. MM: Michaelis–Menten kinetics. TS: thermodynamic span.

10. Other OXPHOS Aspects and Concluding Remarks

In this review, we focus on the mathematical modeling of the OXPHOS complexes, how their reaction mechanisms have been described, and where and how ROS come from. Although we cover the descriptors from complex I to complex V (Table 1), to build an integrated mitochondrial bioenergetics model, additional aspects related to OXPHOS must be addressed.
For instance, ATP generation can be stimulated by mitochondrial calcium, ADP, and inorganic phosphate to meet varying workload demands. Mitochondrial calcium levels increase with cytosolic calcium levels upon excitation-contraction coupling and stimulate dehydrogenases in the CAC. ADP and inorganic phosphate are substrates of complex V to synthesize ATP. Takeuchi et al. [40] conducted a comprehensive review of modeling the tight regulation of ATP levels via calcium, ADP, and inorganic phosphate dynamics.
ROS generation and scavenging regulation is another important topic in mitochondrial OXPHOS modeling. ROS play a crucial role in physiological metabolic signaling and pathological oxidative stress. Pereira et al. [170] and Mazat et al. [171] discussed the mathematical ETC models of ROS production that could be further extended in studies of oxidative stress and signaling. For instance, Kembro et al. [14] built a two-compartment (cytosolic and mitochondrial) model integrating ROS generation from the ETC, ROS scavenging by SOD, catalase, the GSH and the Trx systems, NADPH production, and superoxide-induced mitochondrial depolarization. The authors found that mitochondrial membrane potential oscillation frequencies and amplitudes were dependent on the counteracting forces of ROS production and the antioxidant capacity [172].
Mitochondrial OXPHOS coordinates with other mitochondrial components to regulate ATP production, transportation, and utilization. ATP generated from complex V must be exported from the mitochondrial matrix to the cytosol via ANT, which catalyzes the 1:1 exchange of free ATP to free ADP across the IMM [142,173].
The creatine kinase (CK) system shuttles and buffers ATP synthesis from the OXPHOS and consumption by biological processes, exchanging high energy phosphates between ATP/ADP and creatine/creatine phosphate couples [174]. Mitochondrial CK (mtCK) couples functionally and structurally with ANT to channel ADP recycling, enhance ATP synthesis, and reduce ROS generation [175,176,177]. Mathematical modeling has shown that mtCK can maintain cellular homeostasis and dampen mitochondrial metabolic oscillation [178].
The ETC is fueled by the CAC. CAC-generated NADH is oxidized in complex I and the CAC intermediate succinate is oxidized in complex II. CAC is coupled with carbohydrate (glycolysis), lipid (beta-oxidation) and amino acid metabolism. Some models keep these inputs as constant parameters [56] for simplicity, while others study the impact of varying fuel inputs. For instance, in a glucose-sensing beta cell model by Fridlyand et al. [179], glycolysis flux can enhance CAC, ETC, and mitochondrial ATP synthesis fluxes and induce insulin secretion. A cardiomyocyte mitochondrial model by Cortassa et al. [165] showed that excessive fatty acids can lead to increased ROS production and mitochondrial uncoupling, both computationally and experimentally.
The ETC architecture in vivo could be more complex than the homogeneous assumption in the ODE-based mathematical models. Respiratory complexes I, III dimer, and IV could assemble into a supercomplex, channeling the electron carriers of ubiquinone and cytochrome c between the complexes to increase efficiency, which may deviate from the pool behavior assumed in ODE models [1,180,181,182]. The ETC complexes have also been shown to concentrate on the cristae membrane, the infoldings of the inner mitochondrial membrane, with complex V (ATP synthase) dimers sitting on the apex of the cristae. The membrane potential was shown to be higher across the crista IMM than across the non-crista IMM. Each crista can hold different membrane potentials, acting as independent ATP-generating units [125].
To date, several integrative cellular models have been developed to study the interactions of mitochondrial energetics and the cellular environment. Computational models have been used to study the energy supply and demand in cardiomyocytes [37,136,177,183], glucose sensing and insulin secretion in pancreatic beta cells [156,179], metabolic interactions between neurons and astrocytes in the brain [184,185,186], hepatocellular respiration in the liver [187], and metabolic adaptation during exercise in the skeletal muscle [188,189]. Due to tissue differences in energetics demands, mitochondrial distributions and functions are very different in each tissue [190]. To consider mitochondrial energetics, the model parameters must be adjusted and validated by tissue-specific data. The mitochondria and OXPHOS model could also apply to heart failure studies or iPSC studies in patients. For example, induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) provide a tool for disease modeling, including mitochondrial cardiomyopathy and metabolic disorders with cardiac phenotypes [191]. Cardiac computational modeling [192,193] with bioenergetics details can provide insight into the mechanisms in disease development.
In this review, we focus on deterministic kinetic models of OXPHOS, which emphasize the dynamic, continuous, and time-dependent interactions of each complex and substrate. Stochastic models of the ETC complex that focus on discrete events of electron transfer within a complex are also available [65,90,194,195]. While these models provide more details for the underlying molecular mechanism than deterministic kinetic models, they require more simulation time than deterministic models. They are more challenging to integrate with other modeling components. Another way to build OXPHOS ATP generation models is constraint-based modeling (CBM), which considers biochemical reactions as a whole and uses a stoichiometric matrix of metabolites and reactions to deduce the best plausible steady-state target flux (e.g., ATP synthesis rate) under some constraints in the reaction fluxes [196]. Although CBM requires far less parameters to fit and could handle genome-scale metabolic model (GEMMs) containing thousands of reactions and hundreds of metabolites, CBM cannot capture transient dynamic responses and complex interactions between components in the same way as kinetic models [197].
In conclusion, we review various mathematical kinetic modeling strategies for OXPHOS complexes. The OXPHOS is a highly nonlinear system with complicated time-variant responses on a range of inputs, such as respiring substrates, ions (including protons, phosphate, magnesium, calcium, sodium, and potassium), and the mitochondrial membrane potential. Lumped models require the least number of parameters and intermediate state variables but omit some mechanistic details. Mass action kinetic models give the most detailed inner workings but include many parameters that must be fitted and intermediate state variables that must be computed. Quasi-steady state kinetics models for individual complexes sit in between the two. The OXPHOS system is also interconnected to other mitochondrial bioenergetics systems, as represented mathematically by their inputs: dependence on respiring substrate concentration and the pmf, and by their outputs: generating ATP, pmf, and ROS. We also briefly discuss other bioenergetics topics not covered in detail, such as the CAC, ATP export, the creatine shuttle, ROS production/scavenging, and tissue/species differences. These OXPHOS modeling strategies and other interconnected systems can help us to build an integrative bioenergetics model for cellular metabolism.

Author Contributions

Conceptualization, W.-W.T. and A.-C.W.; investigation, W.-W.T.; writing—original draft preparation, W.-W.T. and A.-C.W.; writing—review and editing, W.-W.T. and A.-C.W.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology in Taiwan, grant number MOST 111-2636-B-002-018, and the National Health Research Institute, grant number NHRI-EX-111-11121SC.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Letts, J.A.; Sazanov, L.A. Clarifying the supercomplex: The higher-order organization of the mitochondrial electron transport chain. Nat. Struct. Mol. Biol. 2017, 24, 800–808. [Google Scholar] [CrossRef]
  2. Junge, W.; Nelson, N. ATP synthase. Annu. Rev. Biochem. 2015, 84, 631–657. [Google Scholar] [CrossRef] [Green Version]
  3. Boyer, P.D.; Chance, B.; Ernster, L.; Mitchell, P.; Racker, E.; Slater, E.C. Oxidative phosphorylation and photophosphorylation. Annu. Rev. Biochem. 1977, 46, 955–966. [Google Scholar] [CrossRef]
  4. Mitchell, P. Keilin’s respiratory chain concept and its chemiosmotic consequences. Science 1979, 206, 1148–1159. [Google Scholar] [CrossRef] [Green Version]
  5. Hinkle, P.C. P/O ratios of mitochondrial oxidative phosphorylation. Biochim. Biophys. Acta 2005, 1706, 1–11. [Google Scholar] [CrossRef] [Green Version]
  6. Siekevitz, P. Powerhouse of the cell. Sci. Am. 1957, 197, 131–144. [Google Scholar] [CrossRef]
  7. Wood, P.M. The potential diagram for oxygen at pH 7. Biochem. J. 1988, 253, 287–289. [Google Scholar] [CrossRef]
  8. Jastroch, M.; Divakaruni, A.S.; Mookerjee, S.; Treberg, J.R.; Brand, M.D. Mitochondrial proton and electron leaks. Essays Biochem. 2010, 47, 53–67. [Google Scholar] [CrossRef] [Green Version]
  9. Scialò, F.; Sriram, A.; Fernández-Ayala, D.; Gubina, N.; Lõhmus, M.; Nelson, G.; Logan, A.; Cooper, H.M.; Navas, P.; Enríquez, J.A.; et al. Mitochondrial ROS produced via reverse electron transport extend animal lifespan. Cell Metab. 2016, 23, 725–734. [Google Scholar] [CrossRef] [Green Version]
  10. Scialò, F.; Fernández-Ayala, D.J.; Sanz, A. Role of mitochondrial reverse electron transport in ROS signaling: Potential roles in health and disease. Front. Physiol. 2017, 8, 428. [Google Scholar] [CrossRef]
  11. Brand, M.D. Mitochondrial generation of superoxide and hydrogen peroxide as the source of mitochondrial redox signaling. Free Radic. Biol. Med. 2016, 100, 14–31. [Google Scholar] [CrossRef]
  12. Zhao, R.-Z.; Jiang, S.; Zhang, L.; Yu, Z.-B. Mitochondrial electron transport chain, ROS generation and uncoupling (Review). Int. J. Mol. Med. 2019, 44, 3–15. [Google Scholar] [CrossRef] [Green Version]
  13. Orrenius, S.; Gogvadze, V.; Zhivotovsky, B. Mitochondrial oxidative stress: Implications for cell death. Annu. Rev. Pharmacol. Toxicol. 2007, 47, 143–183. [Google Scholar] [CrossRef]
  14. Kembro, J.M.; Aon, M.A.; Winslow, R.L.; O’Rourke, B.; Cortassa, S. Integrating mitochondrial energetics, redox and ROS metabolic networks: A two-compartment model. Biophys. J. 2013, 104, 332–343. [Google Scholar] [CrossRef] [Green Version]
  15. Aon, M.A.; Stanley, B.A.; Sivakumaran, V.; Kembro, J.M.; O’Rourke, B.; Paolocci, N.; Cortassa, S. Glutathione/thioredoxin systems modulate mitochondrial H2O2 emission: An experimental-computational study. J. Gen. Physiol. 2012, 139, 479–491. [Google Scholar] [CrossRef] [Green Version]
  16. Szendroedi, J.; Phielix, E.; Roden, M. The role of mitochondria in insulin resistance and type 2 diabetes mellitus. Nat. Rev. Endocrinol. 2011, 8, 92–103. [Google Scholar] [CrossRef]
  17. Haythorne, E.; Rohm, M.; van de Bunt, M.; Brereton, M.F.; Tarasov, A.I.; Blacker, T.S.; Sachse, G.; Silva Dos Santos, M.; Terron Exposito, R.; Davis, S.; et al. Diabetes causes marked inhibition of mitochondrial metabolism in pancreatic β-cells. Nat. Commun. 2019, 10, 2474. [Google Scholar] [CrossRef] [Green Version]
  18. Sivitz, W.I.; Yorek, M.A. Mitochondrial dysfunction in diabetes: From molecular mechanisms to functional significance and therapeutic opportunities. Antioxid. Redox Signal. 2010, 12, 537–577. [Google Scholar] [CrossRef] [Green Version]
  19. Rovira-Llopis, S.; Bañuls, C.; Diaz-Morales, N.; Hernandez-Mijares, A.; Rocha, M.; Victor, V.M. Mitochondrial dynamics in type 2 diabetes: Pathophysiological implications. Redox Biol. 2017, 11, 637–645. [Google Scholar] [CrossRef]
  20. Anello, M.; Lupi, R.; Spampinato, D.; Piro, S.; Masini, M.; Boggi, U.; Del Prato, S.; Rabuazzo, A.M.; Purrello, F.; Marchetti, P. Functional and morphological alterations of mitochondria in pancreatic beta cells from type 2 diabetic patients. Diabetologia 2005, 48, 282–289. [Google Scholar] [CrossRef]
  21. Fex, M.; Nicholas, L.M.; Vishnu, N.; Medina, A.; Sharoyko, V.V.; Nicholls, D.G.; Spégel, P.; Mulder, H. The pathogenetic role of β-cell mitochondria in type 2 diabetes. J. Endocrinol. 2018, 236, R145–R159. [Google Scholar] [CrossRef] [Green Version]
  22. Mihajlovic, M.; Vinken, M. Mitochondria as the Target of Hepatotoxicity and Drug-Induced Liver Injury: Molecular Mechanisms and Detection Methods. Int. J. Mol. Sci. 2022, 23, 3315. [Google Scholar] [CrossRef]
  23. Nadanaciva, S.; Will, Y. Current concepts in drug-induced mitochondrial toxicity. Curr. Protoc. Toxicol. 2009, 40, 2–15. [Google Scholar] [CrossRef]
  24. Koopman, W.J.H.; Distelmaier, F.; Smeitink, J.A.M.; Willems, P.H.G.M. OXPHOS mutations and neurodegeneration. EMBO J. 2013, 32, 9–29. [Google Scholar] [CrossRef]
  25. Schon, E.A.; Przedborski, S. Mitochondria: The next (neurode)generation. Neuron 2011, 70, 1033–1053. [Google Scholar] [CrossRef] [Green Version]
  26. Procaccini, C.; Santopaolo, M.; Faicchia, D.; Colamatteo, A.; Formisano, L.; de Candia, P.; Galgani, M.; De Rosa, V.; Matarese, G. Role of metabolism in neurodegenerative disorders. Metab. Clin. Exp. 2016, 65, 1376–1390. [Google Scholar] [CrossRef]
  27. Golpich, M.; Amini, E.; Mohamed, Z.; Azman Ali, R.; Mohamed Ibrahim, N.; Ahmadiani, A. Mitochondrial Dysfunction and Biogenesis in Neurodegenerative diseases: Pathogenesis and Treatment. CNS Neurosci. Ther. 2017, 23, 5–22. [Google Scholar] [CrossRef]
  28. Reddy, P.H. Role of mitochondria in neurodegenerative diseases: Mitochondria as a therapeutic target in Alzheimer’s disease. CNS Spectr. 2009, 14, 8–13; discussion 16. [Google Scholar] [CrossRef]
  29. Fillmore, N.; Mori, J.; Lopaschuk, G.D. Mitochondrial fatty acid oxidation alterations in heart failure, ischaemic heart disease and diabetic cardiomyopathy. Br. J. Pharmacol. 2014, 171, 2080–2090. [Google Scholar] [CrossRef] [Green Version]
  30. Bayeva, M.; Gheorghiade, M.; Ardehali, H. Mitochondria as a therapeutic target in heart failure. J. Am. Coll. Cardiol. 2013, 61, 599–610. [Google Scholar] [CrossRef]
  31. El-Hattab, A.W.; Scaglia, F. Mitochondrial Cardiomyopathies. Front. Cardiovasc. Med. 2016, 3, 25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Murphy, E.; Ardehali, H.; Balaban, R.S.; DiLisa, F.; Dorn, G.W.; Kitsis, R.N.; Otsu, K.; Ping, P.; Rizzuto, R.; Sack, M.N.; et al. American Heart Association Council on Basic Cardiovascular Sciences, Council on Clinical Cardiology, and Council on Functional Genomics and Translational Biology Mitochondrial function, biology, and role in disease: A scientific statement from the american heart association. Circ. Res. 2016, 118, 1960–1991. [Google Scholar] [CrossRef] [PubMed]
  33. Varga, Z.V.; Ferdinandy, P.; Liaudet, L.; Pacher, P. Drug-induced mitochondrial dysfunction and cardiotoxicity. Am. J. Physiol. Heart Circ. Physiol. 2015, 309, H1453–H1467. [Google Scholar] [CrossRef] [Green Version]
  34. Yang, K.-C.; Bonini, M.G.; Dudley, S.C. Mitochondria and arrhythmias. Free Radic. Biol. Med. 2014, 71, 351–361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhou, B.; Tian, R. Mitochondrial dysfunction in pathophysiology of heart failure. J. Clin. Investig. 2018, 128, 3716–3726. [Google Scholar] [CrossRef] [Green Version]
  36. Piquereau, J.; Caffin, F.; Novotova, M.; Lemaire, C.; Veksler, V.; Garnier, A.; Ventura-Clapier, R.; Joubert, F. Mitochondrial dynamics in the adult cardiomyocytes: Which roles for a highly specialized cell? Front. Physiol. 2013, 4, 102. [Google Scholar] [CrossRef] [Green Version]
  37. Tran, K.; Loiselle, D.S.; Crampin, E.J. Regulation of cardiac cellular bioenergetics: Mechanisms and consequences. Physiol. Rep. 2015, 3, e12464. [Google Scholar] [CrossRef]
  38. D’Oria, R.; Schipani, R.; Leonardini, A.; Natalicchio, A.; Perrini, S.; Cignarelli, A.; Laviola, L.; Giorgino, F. The role of oxidative stress in cardiac disease: From physiological response to injury factor. Oxid. Med. Cell. Longev. 2020, 2020, 5732956. [Google Scholar] [CrossRef] [PubMed]
  39. Cortassa, S.; Sollott, S.J.; Aon, M.A. Computational Modeling of Mitochondrial Function from a Systems Biology Perspective. Methods Mol. Biol. 2018, 1782, 249–265. [Google Scholar] [CrossRef]
  40. Takeuchi, A.; Matsuoka, S. Integration of mitochondrial energetics in heart with mathematical modelling. J. Physiol. 2020, 598, 1443–1457. [Google Scholar] [CrossRef]
  41. Winslow, R.L.; Cortassa, S.; Greenstein, J.L. Using models of the myocyte for functional interpretation of cardiac proteomic data. J. Physiol. 2005, 563, 73–81. [Google Scholar] [CrossRef]
  42. Cortassa, S.; Aon, M.A. Computational modeling of mitochondrial function. Methods Mol. Biol. 2012, 810, 311–326. [Google Scholar] [CrossRef] [Green Version]
  43. Imperial, S.; Centelles, J.J. Enzyme kinetic equations of irreversible and reversible reactions in metabolism. JBM 2014, 2, 24–29. [Google Scholar] [CrossRef]
  44. Alberty, R.A. Relations between biochemical thermodynamics and biochemical kinetics. Biophys. Chem. 2006, 124, 11–17. [Google Scholar] [CrossRef]
  45. Gunawardena, J. A linear framework for time-scale separation in nonlinear biochemical systems. PLoS ONE 2012, 7, e36321. [Google Scholar] [CrossRef] [Green Version]
  46. Srinivasan, B. A guide to the Michaelis-Menten equation: Steady state and beyond. FEBS J. 2021, 289, 6086–6098. [Google Scholar] [CrossRef]
  47. Einav, T.; Mazutis, L.; Phillips, R. Statistical mechanics of allosteric enzymes. J. Phys. Chem. B 2016, 120, 6021–6037. [Google Scholar] [CrossRef] [Green Version]
  48. Leskovac, V.; Trivić, S.; Pericin, D.; Kandrac, J. Deriving the rate equations for product inhibition patterns in bisubstrate enzyme reactions. J. Enzyme Inhib. Med. Chem. 2006, 21, 617–634. [Google Scholar] [CrossRef]
  49. Cornish-Bowden, A.; Hofmeyr, J.-H.S. Enzymes in context: Kinetic characterization of enzymes for systems biology. Biochemist 2005, 27, 11–14. [Google Scholar] [CrossRef] [Green Version]
  50. King, E.L.; Altman, C. A Schematic Method of Deriving the Rate Laws for Enzyme-Catalyzed Reactions. J. Phys. Chem. 1956, 60, 1375–1378. [Google Scholar] [CrossRef]
  51. Qi, F.; Dash, R.K.; Han, Y.; Beard, D.A. Generating rate equations for complex enzyme systems by a computer-assisted systematic method. BMC Bioinform. 2009, 10, 238. [Google Scholar] [CrossRef] [Green Version]
  52. Korzeniewski, B.; Froncisz, W. An extended dynamic model of oxidative phosphorylation. Biochim. Biophys. Acta BBA Bioenerg. 1991, 1060, 210–223. [Google Scholar] [CrossRef]
  53. Korzeniewski, B.; Mazat, J.P. Theoretical studies on the control of oxidative phosphorylation in muscle mitochondria: Application to mitochondrial deficiencies. Biochem. J. 1996, 319 Pt 1, 143–148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Magnus, G.; Keizer, J. Minimal model of beta-cell mitochondrial Ca2+ handling. Am. J. Physiol. 1997, 273, C717–C733. [Google Scholar] [CrossRef]
  55. Pietrobon, D.; Caplan, S.R. Flow-force relationships for a six-state proton pump model: Intrinsic uncoupling, kinetic equivalence of input and output forces, and domain of approximate linearity. Biochemistry 1985, 24, 5764–5776. [Google Scholar] [CrossRef]
  56. Cortassa, S.; Aon, M.A.; Marbán, E.; Winslow, R.L.; O’Rourke, B. An integrated model of cardiac mitochondrial energy metabolism and calcium dynamics. Biophys. J. 2003, 84, 2734–2755. [Google Scholar] [CrossRef] [Green Version]
  57. Baradaran, R.; Berrisford, J.M.; Minhas, G.S.; Sazanov, L.A. Crystal structure of the entire respiratory complex I. Nature 2013, 494, 443–448. [Google Scholar] [CrossRef] [Green Version]
  58. Sazanov, L.A. A giant molecular proton pump: Structure and mechanism of respiratory complex I. Nat. Rev. Mol. Cell Biol. 2015, 16, 375–388. [Google Scholar] [CrossRef]
  59. Vinogradov, A.D.; Grivennikova, V.G. Oxidation of NADH and ROS production by respiratory complex I. Biochim. Biophys. Acta 2016, 1857, 863–871. [Google Scholar] [CrossRef]
  60. Kussmaul, L.; Hirst, J. The mechanism of superoxide production by NADH:ubiquinone oxidoreductase (complex I) from bovine heart mitochondria. Proc. Natl. Acad. Sci. USA 2006, 103, 7607–7612. [Google Scholar] [CrossRef]
  61. Bazil, J.N.; Pannala, V.R.; Dash, R.K.; Beard, D.A. Determining the origins of superoxide and hydrogen peroxide in the mammalian NADH:ubiquinone oxidoreductase. Free Radic. Biol. Med. 2014, 77, 121–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Gauthier, L.D.; Greenstein, J.L.; Cortassa, S.; O’Rourke, B.; Winslow, R.L. A computational model of reactive oxygen species and redox balance in cardiac mitochondria. Biophys. J. 2013, 105, 1045–1056. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Markevich, N.I.; Hoek, J.B. Computational modeling analysis of mitochondrial superoxide production under varying substrate conditions and upon inhibition of different segments of the electron transport chain. Biochim. Biophys. Acta 2015, 1847, 656–679. [Google Scholar] [CrossRef] [Green Version]
  64. Brand, M.D. The sites and topology of mitochondrial superoxide production. Exp. Gerontol. 2010, 45, 466–472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Ransac, S.; Arnarez, C.; Mazat, J.-P. The flitting of electrons in complex I: A stochastic approach. Biochim. Biophys. Acta 2010, 1797, 641–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Korzeniewski, B. Regulation of ATP supply during muscle contraction: Theoretical studies. Biochem. J. 1998, 330 Pt 3, 1189–1195. [Google Scholar] [CrossRef] [Green Version]
  67. Korzeniewski, B.; Zoladz, J.A. A model of oxidative phosphorylation in mammalian skeletal muscle. Biophys. Chem. 2001, 92, 17–34. [Google Scholar] [CrossRef]
  68. Korzeniewski, B.; Noma, A.; Matsuoka, S. Regulation of oxidative phosphorylation in intact mammalian heart in vivo. Biophys. Chem. 2005, 116, 145–157. [Google Scholar] [CrossRef]
  69. Beard, D.A. A biophysical model of the mitochondrial respiratory system and oxidative phosphorylation. PLoS Comput. Biol. 2005, 1, e36. [Google Scholar] [CrossRef] [Green Version]
  70. Bose, S.; French, S.; Evans, F.J.; Joubert, F.; Balaban, R.S. Metabolic network control of oxidative phosphorylation: Multiple roles of inorganic phosphate. J. Biol. Chem. 2003, 278, 39155–39165. [Google Scholar] [CrossRef]
  71. Heiske, M.; Letellier, T.; Klipp, E. Comprehensive mathematical model of oxidative phosphorylation valid for physiological and pathological conditions. FEBS J. 2017, 284, 2802–2828. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Gauthier, L.D.; Greenstein, J.L.; O’Rourke, B.; Winslow, R.L. An integrated mitochondrial ROS production and scavenging model: Implications for heart failure. Biophys. J. 2013, 105, 2832–2842. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Dhingra, R.; Kirshenbaum, L.A. Succinate dehydrogenase/complex II activity obligatorily links mitochondrial reserve respiratory capacity to cell survival in cardiac myocytes. Cell Death Dis. 2015, 6, e1956. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Sun, F.; Huo, X.; Zhai, Y.; Wang, A.; Xu, J.; Su, D.; Bartlam, M.; Rao, Z. Crystal structure of mitochondrial respiratory membrane protein complex II. Cell 2005, 121, 1043–1057. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Chouchani, E.T.; Pell, V.R.; Gaude, E.; Aksentijević, D.; Sundier, S.Y.; Robb, E.L.; Logan, A.; Nadtochiy, S.M.; Ord, E.N.J.; Smith, A.C.; et al. Ischaemic accumulation of succinate controls reperfusion injury through mitochondrial ROS. Nature 2014, 515, 431–435. [Google Scholar] [CrossRef] [Green Version]
  76. Stepanova, A.; Shurubor, Y.; Valsecchi, F.; Manfredi, G.; Galkin, A. Differential susceptibility of mitochondrial complex II to inhibition by oxaloacetate in brain and heart. Biochim. Biophys. Acta 2016, 1857, 1561–1568. [Google Scholar] [CrossRef] [Green Version]
  77. Wei, A.-C.; Aon, M.A.; O’Rourke, B.; Winslow, R.L.; Cortassa, S. Mitochondrial energetics, pH regulation, and ion dynamics: A computational-experimental approach. Biophys. J. 2011, 100, 2894–2903. [Google Scholar] [CrossRef] [Green Version]
  78. Demin, O.V.; Gorianin, I.I.; Kholodenko, B.N.; Westerhoff, H.V. Kinetic modeling of energy metabolism and generation of active forms of oxygen in hepatocyte mitochondria. Mol. Biol. 2001, 35, 1095–1104. [Google Scholar] [CrossRef]
  79. Demin, O.V.; Kholodenko, B.N.; Skulachev, V.P. A model of O2.-generation in the complex III of the electron transport chain. Mol. Cell. Biochem. 1998, 184, 21–33. [Google Scholar] [CrossRef]
  80. Mogilevskaya, E.; Demin, O.; Goryanin, I. Kinetic model of mitochondrial Krebs cycle: Unraveling the mechanism of salicylate hepatotoxic effects. J. Biol. Phys. 2006, 32, 245–271. [Google Scholar] [CrossRef]
  81. Bazil, J.N.; Beard, D.A.; Vinnakota, K.C. Catalytic coupling of oxidative phosphorylation, ATP demand, and reactive oxygen species generation. Biophys. J. 2016, 110, 962–971. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Wu, F.; Yang, F.; Vinnakota, K.C.; Beard, D.A. Computer modeling of mitochondrial tricarboxylic acid cycle, oxidative phosphorylation, metabolite transport, and electrophysiology. J. Biol. Chem. 2007, 282, 24525–24537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Bazil, J.N.; Buzzard, G.T.; Rundell, A.E. Modeling mitochondrial bioenergetics with integrated volume dynamics. PLoS Comput. Biol. 2010, 6, e1000632. [Google Scholar] [CrossRef] [Green Version]
  84. Manhas, N.; Duong, Q.V.; Lee, P.; Richardson, J.D.; Robertson, J.D.; Moxley, M.A.; Bazil, J.N. Computationally modeling mammalian succinate dehydrogenase kinetics identifies the origins and primary determinants of ROS production. J. Biol. Chem. 2020, 295, 15262–15279. [Google Scholar] [CrossRef] [PubMed]
  85. Markevich, N.I.; Markevich, L.N.; Hoek, J.B. Computational modeling analysis of generation of reactive oxygen species by mitochondrial assembled and disintegrated complex II. Front. Physiol. 2020, 11, 557721. [Google Scholar] [CrossRef] [PubMed]
  86. Mitchell, P. Protonmotive redox mechanism of the cytochrome b-c 1 complex in the respiratory chain: Protonmotive ubiquinone cycle. FEBS Lett. 1975, 56, 1–6. [Google Scholar] [CrossRef] [Green Version]
  87. Bleier, L.; Dröse, S. Superoxide generation by complex III: From mechanistic rationales to functional consequences. Biochim. Biophys. Acta 2013, 1827, 1320–1331. [Google Scholar] [CrossRef] [Green Version]
  88. Crofts, A.R.; Holland, J.T.; Victoria, D.; Kolling, D.R.J.; Dikanov, S.A.; Gilbreth, R.; Lhee, S.; Kuras, R.; Kuras, M.G. The Q-cycle reviewed: How well does a monomeric mechanism of the bc(1) complex account for the function of a dimeric complex? Biochim. Biophys. Acta 2008, 1777, 1001–1019. [Google Scholar] [CrossRef] [Green Version]
  89. Crofts, A.R.; Hong, S.; Wilson, C.; Burton, R.; Victoria, D.; Harrison, C.; Schulten, K. The mechanism of ubihydroquinone oxidation at the Qo-site of the cytochrome bc1 complex. Biochim. Biophys. Acta 2013, 1827, 1362–1377. [Google Scholar] [CrossRef] [Green Version]
  90. Ransac, S.; Parisey, N.; Mazat, J.-P. The loneliness of the electrons in the bc1 complex. Biochim. Biophys. Acta 2008, 1777, 1053–1059. [Google Scholar] [CrossRef]
  91. Quinlan, C.L.; Gerencser, A.A.; Treberg, J.R.; Brand, M.D. The mechanism of superoxide production by the antimycin-inhibited mitochondrial Q-cycle. J. Biol. Chem. 2011, 286, 31361–31372. [Google Scholar] [CrossRef] [Green Version]
  92. Chan, K.; Truong, D.; Shangari, N.; O’Brien, P.J. Drug-induced mitochondrial toxicity. Expert Opin. Drug Metab. Toxicol. 2005, 1, 655–669. [Google Scholar] [CrossRef] [PubMed]
  93. Dröse, S.; Brandt, U. The mechanism of mitochondrial superoxide production by the cytochrome bc1 complex. J. Biol. Chem. 2008, 283, 21649–21654. [Google Scholar] [CrossRef] [Green Version]
  94. Muller, F.L.; Liu, Y.; Van Remmen, H. Complex III releases superoxide to both sides of the inner mitochondrial membrane. J. Biol. Chem. 2004, 279, 49064–49073. [Google Scholar] [CrossRef] [Green Version]
  95. Vinnakota, K.C.; Bazil, J.N.; Van den Bergh, F.; Wiseman, R.W.; Beard, D.A. Feedback regulation and time hierarchy of oxidative phosphorylation in cardiac mitochondria. Biophys. J. 2016, 110, 972–980. [Google Scholar] [CrossRef] [Green Version]
  96. Selivanov, V.A.; Votyakova, T.V.; Zeak, J.A.; Trucco, M.; Roca, J.; Cascante, M. Bistability of mitochondrial respiration underlies paradoxical reactive oxygen species generation induced by anoxia. PLoS Comput. Biol. 2009, 5, e1000619. [Google Scholar] [CrossRef] [Green Version]
  97. Guillaud, F.; Dröse, S.; Kowald, A.; Brandt, U.; Klipp, E. Superoxide production by cytochrome bc1 complex: A mathematical model. Biochim. Biophys. Acta 2014, 1837, 1643–1652. [Google Scholar] [CrossRef] [Green Version]
  98. Bazil, J.N.; Vinnakota, K.C.; Wu, F.; Beard, D.A. Analysis of the kinetics and bistability of ubiquinol:cytochrome c oxidoreductase. Biophys. J. 2013, 105, 343–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Bazil, J.N. Analysis of a functional dimer model of ubiquinol cytochrome c oxidoreductase. Biophys. J. 2017, 113, 1599–1612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Konstantinov, A.A. Cytochrome c oxidase: Intermediates of the catalytic cycle and their energy-coupled interconversion. FEBS Lett. 2012, 586, 630–639. [Google Scholar] [CrossRef]
  101. Yoshikawa, S.; Shimada, A. Reaction mechanism of cytochrome c oxidase. Chem. Rev. 2015, 115, 1936–1989. [Google Scholar] [CrossRef]
  102. Wikström, M.; Krab, K.; Sharma, V. Oxygen activation and energy conservation by cytochrome c oxidase. Chem. Rev. 2018, 118, 2469–2490. [Google Scholar] [CrossRef] [Green Version]
  103. Bloch, D.; Belevich, I.; Jasaitis, A.; Ribacka, C.; Puustinen, A.; Verkhovsky, M.I.; Wikström, M. The catalytic cycle of cytochrome c oxidase is not the sum of its two halves. Proc. Natl. Acad. Sci. USA 2004, 101, 529–533. [Google Scholar] [CrossRef] [Green Version]
  104. Muramoto, K.; Ohta, K.; Shinzawa-Itoh, K.; Kanda, K.; Taniguchi, M.; Nabekura, H.; Yamashita, E.; Tsukihara, T.; Yoshikawa, S. Bovine cytochrome c oxidase structures enable O2 reduction with minimization of reactive oxygens and provide a proton-pumping gate. Proc. Natl. Acad. Sci. USA 2010, 107, 7740–7745. [Google Scholar] [CrossRef] [Green Version]
  105. Mason, M.G.; Nicholls, P.; Wilson, M.T.; Cooper, C.E. Nitric oxide inhibition of respiration involves both competitive (heme) and noncompetitive (copper) binding to cytochrome c oxidase. Proc. Natl. Acad. Sci. USA 2006, 103, 708–713. [Google Scholar] [CrossRef] [Green Version]
  106. Sarti, P.; Forte, E.; Giuffrè, A.; Mastronicola, D.; Magnifico, M.C.; Arese, M. The Chemical Interplay between Nitric Oxide and Mitochondrial Cytochrome c Oxidase: Reactions, Effectors and Pathophysiology. Int. J. Cell Biol. 2012, 2012, 571067. [Google Scholar] [CrossRef] [Green Version]
  107. Jensen, P.; Wilson, M.T.; Aasa, R.; Malmström, B.G. Cyanide inhibition of cytochrome c oxidase. A rapid-freeze e.p.r. investigation. Biochem. J. 1984, 224, 829–837. [Google Scholar] [CrossRef] [Green Version]
  108. Alonso, J.-R.; Cardellach, F.; López, S.; Casademont, J.; Miró, O. Carbon monoxide specifically inhibits cytochrome c oxidase of human mitochondrial respiratory chain. Pharmacol. Toxicol. 2003, 93, 142–146. [Google Scholar] [CrossRef]
  109. Borisov, V.B.; Forte, E. Impact of hydrogen sulfide on mitochondrial and bacterial bioenergetics. Int. J. Mol. Sci. 2021, 22, 12688. [Google Scholar] [CrossRef]
  110. Yoshikawa, S.; Caughey, W.S. Infrared evidence of azide binding to iron, copper, and non-metal sites in heart cytochrome c oxidase. J. Biol. Chem. 1992, 267, 9757–9766. [Google Scholar] [CrossRef]
  111. Liesivuori, J.; Savolainen, H. Methanol and formic acid toxicity: Biochemical mechanisms. Pharmacol. Toxicol. 1991, 69, 157–163. [Google Scholar] [CrossRef]
  112. Hashmi, J.T.; Huang, Y.-Y.; Osmani, B.Z.; Sharma, S.K.; Naeser, M.A.; Hamblin, M.R. Role of low-level laser therapy in neurorehabilitation. PM R 2010, 2, S292–S305. [Google Scholar] [CrossRef] [Green Version]
  113. Antonini, E.; Brunori, M.; Colosimo, A.; Greenwood, C.; Wilson, M.T. Oxygen “pulsed” cytochrome c oxidase: Functional properties and catalytic relevance. Proc. Natl. Acad. Sci. USA 1977, 74, 3128–3132. [Google Scholar] [CrossRef] [Green Version]
  114. Korzeniewski, B. Simulation of oxidative phosphorylation in hepatocytes. Biophys. Chem. 1996, 58, 215–224. [Google Scholar] [CrossRef]
  115. Krab, K.; Kempe, H.; Wikström, M. Explaining the enigmatic K(M) for oxygen in cytochrome c oxidase: A kinetic model. Biochim. Biophys. Acta 2011, 1807, 348–358. [Google Scholar] [CrossRef] [Green Version]
  116. Wilson, D.F.; Harrison, D.K.; Vinogradov, A. Mitochondrial cytochrome c oxidase and control of energy metabolism: Measurements in suspensions of isolated mitochondria. J. Appl. Physiol. 2014, 117, 1424–1430. [Google Scholar] [CrossRef] [Green Version]
  117. Pannala, V.R.; Camara, A.K.S.; Dash, R.K. Modeling the detailed kinetics of mitochondrial cytochrome c oxidase: Catalytic mechanism and nitric oxide inhibition. J. Appl. Physiol. 2016, 121, 1196–1207. [Google Scholar] [CrossRef] [Green Version]
  118. Weber, J.; Senior, A.E. ATP synthase: What we know about ATP hydrolysis and what we do not know about ATP synthesis. Biochim. Biophys. Acta BBA Bioenerg. 2000, 1458, 300–309. [Google Scholar] [CrossRef] [Green Version]
  119. Gout, E.; Rébeillé, F.; Douce, R.; Bligny, R. Interplay of Mg2+, ADP, and ATP in the cytosol and mitochondria: Unravelling the role of Mg2+ in cell respiration. Proc. Natl. Acad. Sci. USA 2014, 111, E4560-7. [Google Scholar] [CrossRef] [Green Version]
  120. Senior, A.E.; Nadanaciva, S.; Weber, J. The molecular mechanism of ATP synthesis by F1F0-ATP synthase. Biochim. Biophys. Acta BBA Bioenerg. 2002, 1553, 188–211. [Google Scholar] [CrossRef]
  121. Petersen, J.; Förster, K.; Turina, P.; Gräber, P. Comparison of the H+/ATP ratios of the H+-ATP synthases from yeast and from chloroplast. Proc. Natl. Acad. Sci. USA 2012, 109, 11150–11155. [Google Scholar] [CrossRef] [Green Version]
  122. Alexandre, A.; Reynafarje, B.; Lehninger, A.L. Stoichiometry of vectorial H+ movements coupled to electron transport and to ATP synthesis in mitochondria. Proc. Natl. Acad. Sci. USA 1978, 75, 5296–5300. [Google Scholar] [CrossRef] [Green Version]
  123. Hong, S.; Pedersen, P.L. ATP synthase and the actions of inhibitors utilized to study its roles in human health, disease, and other scientific areas. Microbiol. Mol. Biol. Rev. 2008, 72, 590–641. [Google Scholar] [CrossRef] [Green Version]
  124. Jonckheere, A.I.; Smeitink, J.A.M.; Rodenburg, R.J.T. Mitochondrial ATP synthase: Architecture, function and pathology. J. Inherit. Metab. Dis. 2012, 35, 211–225. [Google Scholar] [CrossRef] [Green Version]
  125. Wolf, D.M.; Segawa, M.; Kondadi, A.K.; Anand, R.; Bailey, S.T.; Reichert, A.S.; van der Bliek, A.M.; Shackelford, D.B.; Liesa, M.; Shirihai, O.S. Individual cristae within the same mitochondrion display different membrane potentials and are functionally independent. EMBO J. 2019, 38, e101056. [Google Scholar] [CrossRef]
  126. Blum, T.B.; Hahn, A.; Meier, T.; Davies, K.M.; Kühlbrandt, W. Dimers of mitochondrial ATP synthase induce membrane curvature and self-assemble into rows. Proc. Natl. Acad. Sci. USA 2019, 116, 4250–4255. [Google Scholar] [CrossRef] [Green Version]
  127. Campanella, M.; Parker, N.; Tan, C.H.; Hall, A.M.; Duchen, M.R. IF(1): Setting the pace of the F(1)F(o)-ATP synthase. Trends Biochem. Sci. 2009, 34, 343–350. [Google Scholar] [CrossRef]
  128. García-Bermúdez, J.; Cuezva, J.M. The ATPase Inhibitory Factor 1 (IF1): A master regulator of energy metabolism and of cell survival. Biochim. Biophys. Acta 2016, 1857, 1167–1182. [Google Scholar] [CrossRef]
  129. Chinopoulos, C. Mitochondrial consumption of cytosolic ATP: Not so fast. FEBS Lett. 2011, 585, 1255–1259. [Google Scholar] [CrossRef] [Green Version]
  130. Juhaszova, M.; Kobrinsky, E.; Zorov, D.B.; Nuss, H.B.; Yaniv, Y.; Fishbein, K.W.; de Cabo, R.; Montoliu, L.; Gabelli, S.B.; Aon, M.A.; et al. ATP Synthase K+- and H+-Fluxes Drive ATP Synthesis and Enable Mitochondrial K+-”Uniporter” Function: I. Characterization of Ion Fluxes. Function 2022, 3, zqab065. [Google Scholar] [CrossRef]
  131. Juhaszova, M.; Kobrinsky, E.; Zorov, D.B.; Nuss, H.B.; Yaniv, Y.; Fishbein, K.W.; de Cabo, R.; Montoliu, L.; Gabelli, S.B.; Aon, M.A.; et al. ATP Synthase K+- and H+-fluxes Drive ATP Synthesis and Enable Mitochondrial K+-”Uniporter” Function: II. Ion and ATP Synthase Flux Regulation. Function 2022, 3, zqac001. [Google Scholar] [CrossRef] [PubMed]
  132. Garlid, K.D.; Paucek, P. Mitochondrial potassium transport: The K+ cycle. Biochim. Biophys. Acta BBA Bioenerg. 2003, 1606, 23–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  133. Mitchell, P.; Moyle, J. Translocation of some anions cations and acids in rat liver mitochondria. Eur. J. Biochem. 1969, 9, 149–155. [Google Scholar] [CrossRef] [PubMed]
  134. Wei, A.-C.; Liu, T.; Cortassa, S.; Winslow, R.L.; O’Rourke, B. Mitochondrial Ca2+ influx and efflux rates in guinea pig cardiac mitochondria: Low and high affinity effects of cyclosporine A. Biochim. Biophys. Acta 2011, 1813, 1373–1381. [Google Scholar] [CrossRef] [Green Version]
  135. Cortassa, S.; Aon, M.A.; Winslow, R.L.; O’Rourke, B. A mitochondrial oscillator dependent on reactive oxygen species. Biophys. J. 2004, 87, 2060–2073. [Google Scholar] [CrossRef] [Green Version]
  136. Cortassa, S.; Aon, M.A.; O’Rourke, B.; Jacques, R.; Tseng, H.-J.; Marbán, E.; Winslow, R.L. A computational model integrating electrophysiology, contraction, and mitochondrial bioenergetics in the ventricular myocyte. Biophys. J. 2006, 91, 1564–1589. [Google Scholar] [CrossRef] [Green Version]
  137. Zhou, L.; O’Rourke, B. Cardiac mitochondrial network excitability: Insights from computational analysis. Am. J. Physiol. Heart Circ. Physiol. 2012, 302, H2178–H2189. [Google Scholar] [CrossRef] [Green Version]
  138. Zhou, L.; Cortassa, S.; Wei, A.-C.; Aon, M.A.; Winslow, R.L.; O’Rourke, B. Modeling cardiac action potential shortening driven by oxidative stress-induced mitochondrial oscillations in guinea pig cardiomyocytes. Biophys. J. 2009, 97, 1843–1852. [Google Scholar] [CrossRef] [Green Version]
  139. Zhou, L.; Solhjoo, S.; Millare, B.; Plank, G.; Abraham, M.R.; Cortassa, S.; Trayanova, N.; O’Rourke, B. Effects of regional mitochondrial depolarization on electrical propagation: Implications for arrhythmogenesis. Circ. Arrhythm. Electrophysiol. 2014, 7, 143–151. [Google Scholar] [CrossRef] [Green Version]
  140. Cortassa, S.; Aon, M.A.; Juhaszova, M.; Kobrinsky, E.; Zorov, D.B.; Sollott, S.J. Computational modeling of mitochondrial K+- and H+-driven ATP synthesis. J. Mol. Cell. Cardiol. 2022, 165, 9–18. [Google Scholar] [CrossRef]
  141. Nguyen, M.-H.T.; Dudycha, S.J.; Jafri, M.S. Effect of Ca2+ on cardiac mitochondrial energy production is modulated by Na+ and H+ dynamics. Am. J. Physiol. Cell Physiol. 2007, 292, C2004–C2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  142. Metelkin, E.; Demin, O.; Kovács, Z.; Chinopoulos, C. Modeling of ATP-ADP steady-state exchange rate mediated by the adenine nucleotide translocase in isolated mitochondria. FEBS J. 2009, 276, 6942–6955. [Google Scholar] [CrossRef] [PubMed]
  143. Chance, E.M. A computer simulation of oxidative phosphorylation. Comput. Biomed. Res. 1967, 1, 251–264. [Google Scholar] [CrossRef] [PubMed]
  144. Bohnensack, R. Control of energy transformation of mitochondria. Analysis by a quantitative model. Biochim. Biophys. Acta 1981, 634, 203–218. [Google Scholar] [CrossRef] [PubMed]
  145. Holzhütter, H.G.; Henke, W.; Dubiel, W.; Gerber, G. A mathematical model to study short-term regulation of mitochondrial energy transduction. Biochim. Biophys. Acta 1985, 810, 252–268. [Google Scholar] [CrossRef] [PubMed]
  146. Schmitz, J.P.J.; Vanlier, J.; van Riel, N.A.W.; Jeneson, J.A.L. Computational modeling of mitochondrial energy transduction. Crit. Rev. Biomed. Eng. 2011, 39, 363–377. [Google Scholar] [CrossRef] [Green Version]
  147. Mazat, J.-P.; Ransac, S.; Heiske, M.; Devin, A.; Rigoulet, M. Mitochondrial energetic metabolism-some general principles. IUBMB Life 2013, 65, 171–179. [Google Scholar] [CrossRef]
  148. Korzeniewski, B.; Froncisz, W. Theoretical studies on the control of the oxidative phosphorylation system. Biochim. Biophys. Acta 1992, 1102, 67–75. [Google Scholar] [CrossRef]
  149. Korzeniewski, B. The modeling of oxidative phosphorylation in skeletal muscle. Jpn. J. Physiol. 2004, 54, 511–516. [Google Scholar] [CrossRef] [Green Version]
  150. Saito, R.; Takeuchi, A.; Himeno, Y.; Inagaki, N.; Matsuoka, S. A simulation study on the constancy of cardiac energy metabolites during workload transition. J. Physiol. 2016, 594, 6929–6945. [Google Scholar] [CrossRef]
  151. LaNoue, K.; Nicklas, W.J.; Williamson, J.R. Control of citric acid cycle activity in rat heart mitochondria. J. Biol. Chem. 1970, 245, 102–111. [Google Scholar] [CrossRef]
  152. Heiske, M.; Nazaret, C.; Mazat, J.-P. Modeling the respiratory chain complexes with biothermokinetic equations—The case of complex I. Biochim. Biophys. Acta 2014, 1837, 1707–1716. [Google Scholar] [CrossRef]
  153. Duong, Q.V.; Levitsky, Y.; Dessinger, M.J.; Strubbe-Rivera, J.O.; Bazil, J.N. Identifying Site-Specific Superoxide and Hydrogen Peroxide Production Rates from the Mitochondrial Electron Transport System Using a Computational Strategy. Function 2021, 2, zqab050. [Google Scholar] [CrossRef]
  154. Malyala, S.; Zhang, Y.; Strubbe, J.O.; Bazil, J.N. Calcium phosphate precipitation inhibits mitochondrial energy metabolism. PLoS Comput. Biol. 2019, 15, e1006719. [Google Scholar] [CrossRef] [Green Version]
  155. Magnus, G.; Keizer, J. Model of β-cell mitochondrial calcium handling and electrical activity. I. Cytoplasmic variables. Am. J. Physiol. Cell Physiol. 1998, 274, C1158–C1173. [Google Scholar] [CrossRef]
  156. Magnus, G.; Keizer, J. Model of β-cell mitochondrial calcium handling and electrical activity. II. Mitochondrial variables. Am. J. Physiol. Cell Physiol. 1998, 274, C1174–C1184. [Google Scholar] [CrossRef]
  157. Jafri, M.S.; Dudycha, S.J. Increases in mitochondrial Ca can result in increases in NADH production. Biophys. J. 2001, 80, 589. [Google Scholar]
  158. Bertram, R.; Gram Pedersen, M.; Luciani, D.S.; Sherman, A. A simplified model for mitochondrial ATP production. J. Theor. Biol. 2006, 243, 575–586. [Google Scholar] [CrossRef] [Green Version]
  159. Saa, A.; Siqueira, K.M. Modeling the ATP production in mitochondria. Bull. Math. Biol. 2013, 75, 1636–1651. [Google Scholar] [CrossRef] [Green Version]
  160. Nguyen, M.-H.T.; Jafri, M.S. Mitochondrial calcium signaling and energy metabolism. Ann. N. Y. Acad. Sci. 2005, 1047, 127–137. [Google Scholar] [CrossRef]
  161. Aon, M.A.; Cortassa, S.; Marbán, E.; O’Rourke, B. Synchronized whole cell oscillations in mitochondrial metabolism triggered by a local release of reactive oxygen species in cardiac myocytes. J. Biol. Chem. 2003, 278, 44735–44744. [Google Scholar] [CrossRef] [PubMed]
  162. Zhou, L.; Aon, M.A.; Almas, T.; Cortassa, S.; Winslow, R.L.; O’Rourke, B. A reaction-diffusion model of ROS-induced ROS release in a mitochondrial network. PLoS Comput. Biol. 2010, 6, e1000657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  163. Yang, R.; Ernst, P.; Song, J.; Liu, X.M.; Huke, S.; Wang, S.; Zhang, J.J.; Zhou, L. Mitochondrial-Mediated Oxidative Ca2+/Calmodulin-Dependent Kinase II Activation Induces Early Afterdepolarizations in Guinea Pig Cardiomyocytes: An in Silico Study. J. Am. Heart Assoc. 2018, 7, e008939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  164. Li, Q.; Su, D.; O’Rourke, B.; Pogwizd, S.M.; Zhou, L. Mitochondria-derived ROS bursts disturb Ca2+ cycling and induce abnormal automaticity in guinea pig cardiomyocytes: A theoretical study. Am. J. Physiol. Heart Circ. Physiol. 2015, 308, H623–H636. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  165. Cortassa, S.; Sollott, S.J.; Aon, M.A. Mitochondrial respiration and ROS emission during β-oxidation in the heart: An experimental-computational study. PLoS Comput. Biol. 2017, 13, e1005588. [Google Scholar] [CrossRef] [Green Version]
  166. Aon, M.A.; Cortassa, S.; O’Rourke, B. Redox-optimized ROS balance: A unifying hypothesis. Biochim. Biophys. Acta 2010, 1797, 865–877. [Google Scholar] [CrossRef] [Green Version]
  167. Kim, N.; Ripple, M.O.; Springett, R. Measurement of the mitochondrial membrane potential and pH gradient from the redox poise of the hemes of the bc1 complex. Biophys. J. 2012, 102, 1194–1203. [Google Scholar] [CrossRef] [Green Version]
  168. De Oliveira, B.L.; Niederer, S. A biophysical systems approach to identifying the pathways of acute and chronic doxorubicin mitochondrial cardiotoxicity. PLoS Comput. Biol. 2016, 12, e1005214. [Google Scholar] [CrossRef] [Green Version]
  169. Wilson, D.F.; Vinogradov, S.A. Mitochondrial cytochrome c oxidase: Mechanism of action and role in regulating oxidative phosphorylation. J. Appl. Physiol. 2014, 117, 1431–1439. [Google Scholar] [CrossRef] [Green Version]
  170. Pereira, E.J.; Smolko, C.M.; Janes, K.A. Computational Models of Reactive Oxygen Species as Metabolic Byproducts and Signal-Transduction Modulators. Front. Pharmacol. 2016, 7, 457. [Google Scholar] [CrossRef] [Green Version]
  171. Mazat, J.-P.; Devin, A.; Ransac, S. Modelling mitochondrial ROS production by the respiratory chain. Cell. Mol. Life Sci. 2020, 77, 455–465. [Google Scholar] [CrossRef] [PubMed]
  172. Kembro, J.M.; Cortassa, S.; Aon, M.A. Complex oscillatory redox dynamics with signaling potential at the edge between normal and pathological mitochondrial function. Front. Physiol. 2014, 5, 257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  173. Metelkin, E.; Goryanin, I.; Demin, O. Mathematical modeling of mitochondrial adenine nucleotide translocase. Biophys. J. 2006, 90, 423–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  174. Saks, V.A.; Kongas, O.; Vendelin, M.; Kay, L. Role of the creatine/phosphocreatine system in the regulation of mitochondrial respiration. Acta Physiol. Scand. 2000, 168, 635–641. [Google Scholar] [CrossRef]
  175. Meyer, L.E.; Machado, L.B.; Santiago, A.P.S.A.; da-Silva, W.S.; De Felice, F.G.; Holub, O.; Oliveira, M.F.; Galina, A. Mitochondrial creatine kinase activity prevents reactive oxygen species generation: Antioxidant role of mitochondrial kinase-dependent ADP re-cycling activity. J. Biol. Chem. 2006, 281, 37361–37371. [Google Scholar] [CrossRef] [Green Version]
  176. Schlattner, U.; Tokarska-Schlattner, M.; Wallimann, T. Mitochondrial creatine kinase in human health and disease. Biochim. Biophys. Acta 2006, 1762, 164–180. [Google Scholar] [CrossRef]
  177. Vendelin, M.; Lemba, M.; Saks, V.A. Analysis of functional coupling: Mitochondrial creatine kinase and adenine nucleotide translocase. Biophys. J. 2004, 87, 696–713. [Google Scholar] [CrossRef] [Green Version]
  178. Hettling, H.; van Beek, J.H.G.M. Analyzing the functional properties of the creatine kinase system with multiscale “sloppy” modeling. PLoS Comput. Biol. 2011, 7, e1002130. [Google Scholar] [CrossRef]
  179. Fridlyand, L.E.; Philipson, L.H. Glucose sensing in the pancreatic beta cell: A computational systems analysis. Theor. Biol. Med. Model. 2010, 7, 15. [Google Scholar] [CrossRef] [Green Version]
  180. Lenaz, G.; Genova, M.L. Mobility and function of coenzyme Q (ubiquinone) in the mitochondrial respiratory chain. Biochim. Biophys. Acta 2009, 1787, 563–573. [Google Scholar] [CrossRef] [Green Version]
  181. Letts, J.A.; Fiedorczuk, K.; Sazanov, L.A. The architecture of respiratory supercomplexes. Nature 2016, 537, 644–648. [Google Scholar] [CrossRef] [PubMed]
  182. Lenaz, G.; Genova, M.L. Kinetics of integrated electron transfer in the mitochondrial respiratory chain: Random collisions vs. solid state electron channeling. Am. J. Physiol. Cell Physiol. 2007, 292, C1221–C1239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  183. Hatano, A.; Okada, J.-I.; Washio, T.; Hisada, T.; Sugiura, S. Distinct functional roles of cardiac mitochondrial subpopulations revealed by a 3D simulation model. Biophys. J. 2015, 108, 2732–2739. [Google Scholar] [CrossRef] [Green Version]
  184. Occhipinti, R.; Somersalo, E.; Calvetti, D. Astrocytes as the glucose shunt for glutamatergic neurons at high activity: An in silico study. J. Neurophysiol. 2009, 101, 2528–2538. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  185. Cloutier, M.; Bolger, F.B.; Lowry, J.P.; Wellstead, P. An integrative dynamic model of brain energy metabolism using in vivo neurochemical measurements. J. Comput. Neurosci. 2009, 27, 391–414. [Google Scholar] [CrossRef] [Green Version]
  186. Cakir, T.; Alsan, S.; Saybaşili, H.; Akin, A.; Ulgen, K.O. Reconstruction and flux analysis of coupling between metabolic pathways of astrocytes and neurons: Application to cerebral hypoxia. Theor. Biol. Med. Model. 2007, 4, 48. [Google Scholar] [CrossRef] [Green Version]
  187. Yang, Y.; Nadanaciva, S.; Will, Y.; Woodhead, J.L.; Howell, B.A.; Watkins, P.B.; Siler, S.Q. Mitosym®: A mechanistic, mathematical model of hepatocellular respiration and bioenergetics. Pharm. Res. 2015, 32, 1975–1992. [Google Scholar] [CrossRef] [Green Version]
  188. Li, Y.; Lai, N.; Kirwan, J.P.; Saidel, G.M. Computational Model of Cellular Metabolic Dynamics in Skeletal Muscle Fibers during Moderate Intensity Exercise. Cell. Mol. Bioeng. 2012, 5, 92–112. [Google Scholar] [CrossRef] [Green Version]
  189. Dash, R.K.; Dibella, J.A.; Cabrera, M.E. A computational model of skeletal muscle metabolism linking cellular adaptations induced by altered loading states to metabolic responses during exercise. Biomed. Eng. Online 2007, 6, 14. [Google Scholar] [CrossRef] [Green Version]
  190. Benard, G.; Faustin, B.; Passerieux, E.; Galinier, A.; Rocher, C.; Bellance, N.; Delage, J.P.; Casteilla, L.; Letellier, T.; Rossignol, R. Physiological diversity of mitochondrial oxidative phosphorylation. Am. J. Physiol. Cell Physiol. 2006, 291, C1172–C1182. [Google Scholar] [CrossRef] [Green Version]
  191. Parrotta, E.I.; Lucchino, V.; Scaramuzzino, L.; Scalise, S.; Cuda, G. Modeling Cardiac Disease Mechanisms Using Induced Pluripotent Stem Cell-Derived Cardiomyocytes: Progress, Promises and Challenges. Int. J. Mol. Sci. 2020, 21, 4354. [Google Scholar] [CrossRef] [PubMed]
  192. Niederer, S.A.; Lumens, J.; Trayanova, N.A. Computational models in cardiology. Nat. Rev. Cardiol. 2019, 16, 100–111. [Google Scholar] [CrossRef] [PubMed]
  193. Kernik, D.C.; Morotti, S.; Wu, H.; Garg, P.; Duff, H.J.; Kurokawa, J.; Jalife, J.; Wu, J.C.; Grandi, E.; Clancy, C.E. A computational model of induced pluripotent stem-cell derived cardiomyocytes incorporating experimental variability from multiple data sources. J. Physiol. 2019, 597, 4533–4564. [Google Scholar] [CrossRef] [Green Version]
  194. Moser, C.C.; Farid, T.A.; Chobot, S.E.; Dutton, P.L. Electron tunneling chains of mitochondria. Biochim. Biophys. Acta 2006, 1757, 1096–1109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  195. Kim, Y.C.; Hummer, G. Proton-pumping mechanism of cytochrome c oxidase: A kinetic master-equation approach. Biochim. Biophys. Acta 2012, 1817, 526–536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  196. Maarleveld, T.R.; Khandelwal, R.A.; Olivier, B.G.; Teusink, B.; Bruggeman, F.J. Basic concepts and principles of stoichiometric modeling of metabolic networks. Biotechnol. J. 2013, 8, 997–1008. [Google Scholar] [CrossRef] [Green Version]
  197. Maldonado, E.M.; Taha, F.; Rahman, J.; Rahman, S. Systems biology approaches toward understanding primary mitochondrial diseases. Front. Genet. 2019, 10, 19. [Google Scholar] [CrossRef]
Figure 1. Schematic of OXPHOS, including the electron transport chain (ETC) and ATP synthase. The figure was adapted and created from “Electron Transport Chain”, by BioRender.com (accessed on 20 August 2022). Retrieved from https://app.biorender.com/biorender-templates (accessed on 20 August 2022).
Figure 1. Schematic of OXPHOS, including the electron transport chain (ETC) and ATP synthase. The figure was adapted and created from “Electron Transport Chain”, by BioRender.com (accessed on 20 August 2022). Retrieved from https://app.biorender.com/biorender-templates (accessed on 20 August 2022).
Cells 11 04020 g001
Figure 2. (a) Schematic of an enzyme-catalyzed reaction with three stages: the free enzyme (E, stage 1), the substrate-bound form (ES, stage 2), and the product-bound form (EP, stage 3). kij are rate constants transitioning from stage i to stage j. (b) The weight of each stage is the sum of the products of the transitional rates in the three connected graphs: Λ, <, and >.
Figure 2. (a) Schematic of an enzyme-catalyzed reaction with three stages: the free enzyme (E, stage 1), the substrate-bound form (ES, stage 2), and the product-bound form (EP, stage 3). kij are rate constants transitioning from stage i to stage j. (b) The weight of each stage is the sum of the products of the transitional rates in the three connected graphs: Λ, <, and >.
Cells 11 04020 g002
Figure 3. Schematic of the lumped complex I-III-IV model by Magnus and Keizer [54]. The six-state proton pump consists of proton translocation (5 → 6 → 1 → 2), NADH oxidation (2 → 3 → 4), and a slip reaction (5 → 2). Each cycle represents one electron transfer from NADH to oxygen, coupled with the translocation of 6 protons from the matrix to the intermembrane space (IMS). The arrows represent the forward reaction direction, although all reactions are reversible. Reprinted/adapted with permission from Ref. [54].
Figure 3. Schematic of the lumped complex I-III-IV model by Magnus and Keizer [54]. The six-state proton pump consists of proton translocation (5 → 6 → 1 → 2), NADH oxidation (2 → 3 → 4), and a slip reaction (5 → 2). Each cycle represents one electron transfer from NADH to oxygen, coupled with the translocation of 6 protons from the matrix to the intermembrane space (IMS). The arrows represent the forward reaction direction, although all reactions are reversible. Reprinted/adapted with permission from Ref. [54].
Cells 11 04020 g003
Figure 4. Schematic of the complex I model by Gauthier et al. [72]. The seven-state proton pump is adapted from the six-state proton pump model of Magnus and Keizer [54,55]. The model includes proton pumping (5 → 6 → 1→ 2), NADH oxidation (2 → 3 → 4), ubiquinone reduction (4 → 7 → 5), and superoxide generation (4 → 2). Each cycle represents one electron transfer from NADH to ubiquinone or oxygen. The arrows represent the forward reaction direction, although all reactions are reversible. Reprinted/adapted with permission from Ref. [62].
Figure 4. Schematic of the complex I model by Gauthier et al. [72]. The seven-state proton pump is adapted from the six-state proton pump model of Magnus and Keizer [54,55]. The model includes proton pumping (5 → 6 → 1→ 2), NADH oxidation (2 → 3 → 4), ubiquinone reduction (4 → 7 → 5), and superoxide generation (4 → 2). Each cycle represents one electron transfer from NADH to ubiquinone or oxygen. The arrows represent the forward reaction direction, although all reactions are reversible. Reprinted/adapted with permission from Ref. [62].
Cells 11 04020 g004
Figure 5. Schematic of the complex I model by Bazil et al. [61]. The 5 states represent the number of electrons held by complex I. The model includes peroxide generation reactions (green arrows) taking two electrons from complex I to oxygen, superoxide generation reactions (orange arrows) taking one electron from complex I to oxygen, NADH oxidation reactions (blue arrows) taking two electrons from NADH to complex I, and finally ubiquinone reduction reactions (black arrows) taking two electrons from complex I to ubiquinone, coupled with 4-proton translocation. The figure was reprinted/adapted with permission from Ref. [61].
Figure 5. Schematic of the complex I model by Bazil et al. [61]. The 5 states represent the number of electrons held by complex I. The model includes peroxide generation reactions (green arrows) taking two electrons from complex I to oxygen, superoxide generation reactions (orange arrows) taking one electron from complex I to oxygen, NADH oxidation reactions (blue arrows) taking two electrons from NADH to complex I, and finally ubiquinone reduction reactions (black arrows) taking two electrons from complex I to ubiquinone, coupled with 4-proton translocation. The figure was reprinted/adapted with permission from Ref. [61].
Cells 11 04020 g005
Figure 6. Schematic of the complex II model by Manhas et al. [84]. The 5 states represent the number of electrons in complex II. The model includes peroxide generation reactions (green arrows) taking two electrons from complex II to oxygen, superoxide generation reactions (orange arrows) taking one electron from complex II to oxygen, succinate oxidation reactions (blue arrows) taking two electrons from NADH to complex II, and finally ubiquinone reduction reactions (black arrows) taking one electron from complex II to ubiquinone. The figure was reprinted/adapted from Ref. [84].
Figure 6. Schematic of the complex II model by Manhas et al. [84]. The 5 states represent the number of electrons in complex II. The model includes peroxide generation reactions (green arrows) taking two electrons from complex II to oxygen, superoxide generation reactions (orange arrows) taking one electron from complex II to oxygen, succinate oxidation reactions (blue arrows) taking two electrons from NADH to complex II, and finally ubiquinone reduction reactions (black arrows) taking one electron from complex II to ubiquinone. The figure was reprinted/adapted from Ref. [84].
Cells 11 04020 g006
Figure 7. Schematic of the dimeric complex III model by Bazil et al. [99]. The 6 states represent the number of electrons in complex III. The model includes one-electron transfer reactions (orange arrows): ubiquinol oxidation at the Qo site (one electron goes to heme bL, the other goes through ISP, cytochrome c1, and finally to cytochrome c) as well as oxygen reduction to superoxide. The model also has a two-electron transfer reaction: Q reduction at the Qi site. The figure was reprinted/adapted from Ref. [99].
Figure 7. Schematic of the dimeric complex III model by Bazil et al. [99]. The 6 states represent the number of electrons in complex III. The model includes one-electron transfer reactions (orange arrows): ubiquinol oxidation at the Qo site (one electron goes to heme bL, the other goes through ISP, cytochrome c1, and finally to cytochrome c) as well as oxygen reduction to superoxide. The model also has a two-electron transfer reaction: Q reduction at the Qi site. The figure was reprinted/adapted from Ref. [99].
Cells 11 04020 g007
Figure 8. Schematic of the complex IV model by Demin et al. [78] with 4 stages (Y, Yr, YO, and YOH). The figure was reprinted/adapted with permission from Ref. [62].
Figure 8. Schematic of the complex IV model by Demin et al. [78] with 4 stages (Y, Yr, YO, and YOH). The figure was reprinted/adapted with permission from Ref. [62].
Cells 11 04020 g008
Figure 9. Schematic of the complex V model by Magnus and Keizer [54], assuming complex V is a 6-state proton pump operating in reverse. Not accounting for slip reactions (5 → 2), one ATP molecule is generated for every 3 protons translocated from the intermembrane space (IMS) to the matrix. The figure was reprinted/adapted from Magnus and Keizer [54].
Figure 9. Schematic of the complex V model by Magnus and Keizer [54], assuming complex V is a 6-state proton pump operating in reverse. Not accounting for slip reactions (5 → 2), one ATP molecule is generated for every 3 protons translocated from the intermembrane space (IMS) to the matrix. The figure was reprinted/adapted from Magnus and Keizer [54].
Cells 11 04020 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tseng, W.-W.; Wei, A.-C. Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria. Cells 2022, 11, 4020. https://doi.org/10.3390/cells11244020

AMA Style

Tseng W-W, Wei A-C. Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria. Cells. 2022; 11(24):4020. https://doi.org/10.3390/cells11244020

Chicago/Turabian Style

Tseng, Wen-Wei, and An-Chi Wei. 2022. "Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria" Cells 11, no. 24: 4020. https://doi.org/10.3390/cells11244020

APA Style

Tseng, W. -W., & Wei, A. -C. (2022). Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria. Cells, 11(24), 4020. https://doi.org/10.3390/cells11244020

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop