Introduction {#Sec1}
============

The mitochondria are the power plants of the eukaryotic cell. In most species, cell respiration is catalysed by a linear array of three energy-transducing membrane-bound proteins, viz. complexes I, III and IV. Each complex links the catalysed redox reaction (reduction of ubiquinone by NADH, reduction of ferricytochrome *c* by ubiquinol, and reduction of O~2~ by ferrocytochrome *c*, respectively) to net translocation of protons from the inside (mitochondrial matrix) to the outside (intermembrane space) of the inner mitochondrial membrane, thus creating an electrochemical proton gradient (protonmotive force, pmf, Δ*P*, Δ*μ*~H+~) as the primary form of conserved energy. Research over the last 40 years has revealed that the mechanisms of redox-linked proton translocation are fundamentally different for the three complexes^[@CR1]--[@CR5]^. Here, using specific sets of the published data with isolated mitochondria (see below), we address the tightness by which respiratory electron transfer in the three complexes is coupled to transmembrane proton translocation (the "degree of coupling"), and assess the thermodynamic efficiency of this process.

We find that during oxidative phosphorylation (State 3, ref. ^[@CR6]^), when the Δ*P* is submaximal, the average thermodynamic efficiency is ca. 80--90% for all three complexes, and that the degree of coupling of proton translocation to electron transfer is extremely high. When the rate of respiration falls to a resting steady-state level due to a high ATP/ADP^.^P~i~ ratio (State 4), the degree of coupling is significantly lowered and free energy is dissipated as an increase in entropy. This has been ascribed to a proton leak across the inner mitochondrial membrane^[@CR7]--[@CR10]^, and/or to mechanistic "slipping" in the proton pumps themselves^[@CR11]--[@CR13]^. We find that the reduced degree of coupling in State 4 is not uniform among the three respiratory complexes, which suggests control at the mechanistic level of proton pumping.

Results {#Sec2}
=======

The approach {#Sec3}
------------

Conceptually, all three proton-translocating complexes of the respiratory chain (Fig. [1](#Fig1){ref-type="fig"}) function in the same way, although the detailed molecular mechanisms are very different. In all cases, we are dealing with an energetically uphill, endergonic or "driven" reaction of proton translocation from the negatively charged N-side to the positively charged P-side of the inner mitochondrial membrane (defined here as reaction 1). Such a reaction cannot occur spontaneously, but becomes possible when coupled to an energetically downhill, exergonic or "driving" reaction (reaction 2), which is the oxidation of NADH by ubiquinone (catalysed by complex I), oxidation of ubiquinol by ferricytochrome *c* (catalysed by complex III) or oxidation of ferrocytochrome *c* by dioxygen (cataysed by complex IV), or any combination of these three functions (Fig. [1](#Fig1){ref-type="fig"}).Fig. 1The mitochondrial respiratory chain.The protonmotive complexes I, III and IV are depicted in the inner membrane, together with additional components of the chain, complex II, cytochrome *c* and ubiquinone (Q), the hydrogen donors NADH and succinate and the acceptor, O~2~. The combination of N,N,N',N'-tetramethyl-phenylene-diamine (TMPD) and ascorbate (not shown for clarity; see the text) donates electrons directly to cytochrome *c* and thus engages complex IV only. Black arrows depict the hydrogen and electron transfers in the chain, eventually reducing O~2~ to water. Red arrows symbolise the associated protonmotive events where *Z*~n~ is the effective protonmotive H^+^/e^−^ stoichiometry, 2 for complex I and complex IV, 1 for complex III (see "Methods") and n is the name of the complex. Typical average redox potential differences (Δ*E*) across each complex are given. The net free energy change (Δ*G*) for any combination of complexes is given as the arithmetic sum of the relevant redox potential difference and the protonmotive force (Δ*P*) multiplied by the H^+^/e^−^ stoichiometry (*Z*), where Δ*E* has a negative sign and is larger than Δ*P* for the case where redox chemistry drives proton translocation.

Treating this on the general level of thermodynamics of non-equilibrium coupled reactions, one may write (see refs. ^[@CR14]--[@CR16]^)$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{array}{*{20}{c}} {J_1} = {L_{11}X_1 + L_{12}X_2} \\ {J_2} = {L_{21}X_1 + L_{22}X_2} \end{array},$$\end{document}$$where *J*~*1*~ and *J*~*2*~ are the fluxes of the coupled reactions, and *X*~1~ and *X*~2~ refer to the corresponding forces (free energies). In the case of the electron-transport chain complexes, this becomes:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{array}{*{20}{c}} {J_H} = { - L_{11}{\Delta}P - L_{12}{\Delta}E} \\ {J_e} = { - L_{21}{\Delta}P - L_{22}{\Delta}E} \end{array},$$\end{document}$$where *J*~*H*~ is proton flux (negative when protons are pumped from the N to the P-side of the membrane), *J*~*e*~ is the electron flux, Δ*P* is the protonmotive force and Δ*E* is the redox span or difference in redox potential of the redox carriers on the donor and acceptor side of the reaction considered. The cross-coupling coefficients *L*~12~ and *L*~21~ are assumed to be equal, so-called Onsager symmetry, which is met under many conditions of oxidative phosphorylation (see refs. ^[@CR15],[@CR16]^ and below). For a pump, *L*~11~ and *L*~22~ are negative, and *L*~12~ and *L*~21~ are positive.

It is important to note that our treatment of mitochondrial redox-linked proton translocation by equations^[@CR1]^ is a simplification insofar as the proton flux, *J*~*H*~, will be assessed from the proton flux across the F~o~ portion of the ATP synthase leading to ATP synthesis (see below). Equation ([1b](#Equ2){ref-type=""}) does not explicitly include proton flux across the membrane dielectric (proton leak). Such proton leak, as well as incomplete coupling between electrons and protons in the respiratory complexes themselves (i.e. proton or electron "slipping", see refs. ^[@CR7]--[@CR13],[@CR17]^ and Fig. [2](#Fig2){ref-type="fig"}), are instead collectively considered in the "degree of coupling" (*q*) of reactions 1 and 2, as will be defined below.Fig. 2Leaks and slips in the protonmotive function of mitochondria.The cartoon shows a prototype of a protonmotive redox complex of the respiratory chain (left) and a schematic view of the F~1~F~o~ ATP synthase (right), both embedded in the inner mitochondrial phospholipid membrane. Red arrows depict electron transfer from reductant to oxidant, either linked (*J*~*e*~^*P*^) to proton pumping across the membrane (*J*~*H*~^*P*^) or decoupled from it ("electron slip", *J*~*e*~^*S*^). *J*~*H*~^*L*^ depicts an intrinsic proton leak in the pump mechanism ("proton slip"). *J*~*H*~^*M*^ stands for loss of protonmotive force by a proton leak across the membrane. Proton influx across the F~o~ domain of the ATP synthase (*J*~*H*~^*F*^) is assumed to be fully coupled to synthesis of ATP without significant slips (see the text). The stoichiometric aspects (H^+^/e^−^ and H^+^/ATP ratios) are omitted for clarity (but see the text).

Defining the ratios of output and input forces and flows as *x* (=Δ*P*/Δ*E*) and *j* (=*J*~*H*~*/J*~*e*~), respectively, and the mechanistic H^+^/e^−^ stoichiometry, *Z*, as the square root of the ratio between the primary flow/force-coupling coefficients, i.e. *Z* = *√(L*~11~*/L*~22~*)*, one obtains from Eq. (1) ^[@CR14]--[@CR16]^$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\frac{j}{Z} = \frac{{q - Zx}}{{Zxq - 1}},$$\end{document}$$where *q* is given by:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q = \frac{{L_{12}}}{{\sqrt {L_{11}L_{22}} }}.$$\end{document}$$

The factor *q* has been called the degree of coupling between the input and output processes^[@CR14]^ and varies between 0 (no coupling) and 1 (complete coupling). With complete coupling of the two reactions (i.e. *q* = 1), the flux ratio *j* is equal to *−Z*.

The ratio *j*/*Z*, equivalent to (H^+^/e^−^)/*Z*, has been called the "stoichiometric efficiency"^[@CR18]^, which must not be confused with the thermodynamic efficiency (cf. below). *j*/*Z* is also not equivalent to the degree of coupling (*q*), except for the special case where *x* = 0, i.e. when the force of the driven reaction is zero (Δ*P* = 0) and the coupled reaction proceeds without a load (see Eq. ([2](#Equ3){ref-type=""})). This special case has been termed "level flow"^[@CR14]--[@CR16]^, and is a condition usually strived at in most methods of determining the H^+^/e^−^ ratio experimentally (but see below). It is noteworthy that such an experimentally derived H^+^/e^−^ ratio (*j*) at level flow is not equal to the true mechanistic stoichiometry (*Z*), but to that value multiplied by the degree of coupling, i.e. *qZ*.

The focus {#Sec4}
---------

One aim of this work is to assess the thermodynamic efficiency of the redox-driven proton translocation. Since thermodynamic efficiency (*η*) is defined as the ratio between the output and input powers of two coupled reactions, and since power is the product of flux (*J*) and force (*X*) we may write$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\eta = - J_H{\Delta}P/J_e{\Delta}E,$$\end{document}$$where the negative sign signifies that the driven reaction flux (*J*~*H*~) occurs in a direction against its conjugate force (Δ*P*), that is, the net proton flux is dominated by pumping rather than leak. From the earlier definitions of the flow and force ratios (*j* and *x*), this equation simplifies to$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\eta = - jx.$$\end{document}$$

Thus thermodynamic efficiency (*η*) of the whole-electron-transport chain approaches zero in State 4 mitochondria because the proton pumping is equal to the proton leak so that the net proton flux (*J*^*F*^~*H*~, Fig. [2](#Fig2){ref-type="fig"}) approaches zero. However, each individual proton-pumping complex could still be operating at high efficiency in State 4, but the energy transduced by the proton pumps is dissipated by leak through the membrane. In the fully uncoupled state, also called level flow^[@CR14]--[@CR16]^, the protonmotive force (i.e. the output energy, *X*~1~) approaches zero and hence *η* is zero as well. Clearly, the maximum efficiency of the overall reaction is achieved somewhere between static head and level flow (Fig. [3](#Fig3){ref-type="fig"}, for details see below). In a mathematical model of the proton pump of cytochrome *c* oxidase, Kim and Hummer^[@CR19]^ found a very similar bell-shaped dependence of thermodynamic efficiency on the protonmotive force as the theoretical one shown in Fig. [3](#Fig3){ref-type="fig"}.Fig. 3Theoretical dependence of thermodynamic efficiency on force ratio.The figure plots the dependence of efficiency (*η*) on the ratio between output (protonmotive) and input (oxidoreduction) forces normalised by the mechanistic (ideal) H^+^/e^−^ stoichiometry (*Z*) (see also ref. ^[@CR14]^), according to Eq. ([4](#Equ5){ref-type=""}), for four different degrees of coupling (*q*) as defined by Eq. ([7](#Equ8){ref-type=""}). Example curves are given only for the "uphill" case where the redox reaction drives proton translocation. The thin blue lines show the approximate range of force ratios and efficiencies observed in State 3 of isolated rat-liver mitochondria (see Table [1](#Tab1){ref-type="table"}). The area covered by the data is indicated as a blue rectangle.

State 3, where excess of ADP and inorganic phosphate (P~i~) supports a high rate of respiration^[@CR6]^ and concomitant synthesis of ATP driven by the Δ*P*, lies between static head and level flow and is hence the only reasonably well-defined state of isolated mitochondria suitable for determination of thermodynamic efficiency.

Determination of efficiency using Eq. ([5](#Equ6){ref-type=""}) requires knowledge of the force ratio in State 3 (*x*), which can be calculated from literature data where the redox poise of the NAD^+^/NADH, ubiquinone/ubiquinol, ferri-/ferrocytochrome *c* and oxygen/water redox couples have been reported^[@CR20]^. However, the flux ratio (*j*) (i.e. the H^+^/e^−^ ratio) cannot be directly measured in State 3. Fortunately, we have access to a meta-analysis of the amount of ATP formed per electrons transferred (the so-called ATP/2e^−^ or P/O ratio) in the different domains of the mitochondrial respiratory chain in State 3^[@CR21]^. Such ATP/2e^−^ data can be converted to H^+^/2e^−^ ratios (which are twice the flux ratios, j) by the relation$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\mathrm{H}}^ + /2{\mathrm{e}}^ - = {\mathrm{ATP}}/2{\mathrm{e}}^ - \times {\mathrm{H}}^ + /{\mathrm{ATP}}.$$\end{document}$$where H^+^/ATP, i.e. the number of protons translocated from the P-side to the N-side of the inner membrane per synthesised ATP molecule, is 3.67 in animal mitochondria for the case of extramitochondrially produced ATP. Due to the eight c subunits of the membrane ring structure of the ATP synthase in animal mitochondria, eight protons are translocated per full rotation with consequent synthesis of three molecules of ATP^[@CR22]^, yielding 8/3 or 2.67 H^+^/ATP. Further translocation of the formed ATP to the extramitochondrial space, in exchange for ADP and P~i~ in the opposite direction, takes place on the ADP/ATP and phosphate carriers, respectively, in exchange for translocation of one more proton^[@CR22]^.

It is essential to note that the flux ratio, *j*, assessed by this method, only measures the productive proton flux used for ATP synthesis (see Fig. [2](#Fig2){ref-type="fig"} and below).

In order to derive the degree of coupling (*q*) in State 3 mitochondria, it is not sufficient to know the mechanistic stoichiometry (*Z*) and the H^+^/e^−^ ratio (*j*). As mentioned above, *q* = *j*/*Z* only in the special case where there is no load, i.e. where *x* = 0. State 3 mitochondria are not operating at level flow (Δ*P* = 0), but at a considerable force ratio (*x*); Δ*P* is often ca. 170 mV^[@CR7]--[@CR10],[@CR13],[@CR23]^. It is thus necessary to solve Eq. ([3](#Equ4){ref-type=""}) with respect to *q*, yielding$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q = \frac{{Zx - j/Z}}{{1 - jx}}.$$\end{document}$$

Whilst the degree of coupling (*q*) can be considered a constant in this treatment, the H^+^/e^−^ ratio (i.e. the flux ratio *j*) is not, but decreases with increasing load (*x*) from its maximal value of *qZ* at zero load^[@CR14]--[@CR16]^.

Thermodynamic efficiency {#Sec5}
------------------------

Table [1](#Tab1){ref-type="table"} shows the data obtained for the three redox complexes in isolated mitochondria in State 3, either operating together, or assessed separately (see "Methods"). For comparison between the three redox complexes, the force ratio (*x*) is normalised by multiplication with the mechanistic H^+^/e^−^ stoichiometry, *Z*, which is 2, 1 and 2, respectively for complexes I, III and IV (see "Methods"), and the flux ratio (*j*) is normalised by division with *Z*. The data for the individual complexes were obtained as follows (see Fig. [1](#Fig1){ref-type="fig"}). The ATP/2e^−^ for complex I was taken as the difference between the average ratios with NADH-linked substrates and succinate. The ratio for complex IV was based on the average ATP/2e^−^ with TMPD+ ascorbate as substrate. The ratio for complex III was based on the difference between ATP/2e^−^ ratios with succinate and TMPD+ ascorbate as substrates.Table 1Thermodynamic efficiency (*η*) and degree of proton--electron coupling (*q*) in State 3 of isolated rat-liver mitochondria.*Zxj*/*Zηq*Complexes I + III + IV0.82 ± 0.00 (*n* = 3)−0.9779%0.997Complexes III + IV0.82 ± 0.01 (*n* = 4)−0.9679%0.996Complex I0.88 ± 0.01 (*n* = 3)−0.9887%0.999Complex III0.80 ± 0.04 (*n* = 4)−0.9878%0.998Complex IV0.81 ± 0.03 (*n* = 5)−0.9577%0.995(Complex IV)0.70 ± 0.02 (*n* = 5)−0.9566%0.991*x* is the force ratio, between the protonmotive force (assumed at 170 mV) and the redox potential difference across the measured span of the respiratory chain. Extramitochondrial pH is 7.4; matrix pH = 8.0. Midpoint potentials at pH = 7.4 are 230 mV for cytochrome *c*, 66 mV for ubiquinone and 776 mV for O~2~/2H~2~O (at 120 μM O~2~). Midpoint potentials at pH = 8 are 30 mV for ubiquinone and −350 mV for NAD/NADH (see "Methods" and the text). *Z* is the mechanistic H^+^/e^−^ stoichiometry, which is 2, 1 and 2 for complexes I, III and IV, respectively (see "Methods"). *j* is the observed flux ratio (H^+^/e^−^). The values for complex IV as well as for complexes I + III + IV and III + IV are for the cases in which the energy of the irreversible splitting of the O--O bond has been subtracted and the concentration of O~2~ lowered to 12 μM (see the text). The value for complex IV in parenthesis is without subtraction of O--O bond splitting and at a mean O~2~ concentration of 120 μM. Values are expressed as mean ± SD, where *n* is the number of distinct samples.

As shown in Table [1](#Tab1){ref-type="table"}, the calculated thermodynamic efficiency (*η*) varies between \~67 and 87%, the highest being the efficiency of complex I. Two sets of the data are listed for complex IV, which is explained in the next section.

The special case of complex IV (cytochrome *c* oxidase) {#Sec6}
-------------------------------------------------------

Conventional treatment of the complex IV data (Table [1](#Tab1){ref-type="table"}, data within parentheses) indicates that the terminal enzyme of the respiratory chain would have by far the lowest energy transduction efficiency, 67%, as compared with the efficiencies of complex III (78%) and complex I (87%). One may suspect that this is related to the fact that the overall complex IV reaction is irreversible, so that generation of O~2~ from water is not possible as driven by the protonmotive force. However, elucidation of the mechanism of complex IV catalysis^[@CR1],[@CR24]^ has revealed that the irreversibility is the result of an irreversible splitting of the O--O bond in the partial reaction step due to the large Δ*G*^0^, where the oxygen-bound intermediate A is converted into intermediate P (see ref. ^[@CR25]^, Fig. [4](#Fig4){ref-type="fig"}), a reaction that is neither coupled to proton pumping nor to transmembrane electron or proton transfer^[@CR1],[@CR24]^. The binding of O~2~ to intermediate R is fully reversible^[@CR26],[@CR27]^, and most important, all the other four partial reactions of the catalytic cycle are coupled to proton translocation^[@CR1],[@CR24]^, and have indeed been shown to be reversed by protonmotive force^[@CR28]^.Fig. 4Catalytic intermediates of the binuclear centre of complex IV.The binuclear centre is shown as a rectangle, consisting of haem *a*~*3*~ (Fe) the nearby Cu~B~ (Cu) and the conserved tyrosine (tyr-OH, tyrosine; tyr-O^−^, tyrosinate; tyr-O\*, neutral tyrosine radical) covalently linked to one of the three histidine ligands of Cu~B~. The three histidine ligands of Cu~B~ and the proximal histidine of haem *a*~*3*~ are not shown for clarity. The name of each intermediate state is found outside the right corner of each rectangle. The red arrows specify reactions that each involve electron transfer to the site from the low spin haem *a*, net transfer of one substrate proton from the N-side of the membrane to complete oxygen- reduction chemistry, and uptake (from the N-side) and release (to the P-side) of another proton (i.e. proton pumping). Hence, the red arrows indicate partial redox reactions that are reversible at high Δ*P*. The blue arrow depicts the irreversible partial reaction step where the O--O bond of O~2~ is broken in a reaction not linked to energy transduction. Note that the classical "ferrous-oxy" intermediate A^[@CR26]^ is more accurate described as a ferric-superoxide state of haem *a*~*3*~^[@CR5]^.

The irreversibility of the overall cytochrome *c* oxidase reaction should therefore not be taken to suggest loose coupling of the redox chemistry to proton translocation. In order to focus analysis on these coupled reaction steps, we have subtracted the free energy change of the irreversible reaction step A → P (Fig. [4](#Fig4){ref-type="fig"}; ca. 5 kcal/mol^[@CR25]^) from the overall standard driving force. In the basic analysis (Table [1](#Tab1){ref-type="table"}, data in parentheses), we have assumed an average concentration of O~2~ of 120 μM in calculating the driving force. However, as shown by Wilson et al.^[@CR29]^, the redox potential of the donor, cytochrome *c*, is not lowered (cyt. *c* is not getting more reduced) in the steady state until the oxygen concentration is lowered below ca. 10 μM. We have therefore adjusted the redox potential of the O~2~/2H~2~O couple to 12 μM O~2~ to obtain a better estimate of the true efficiency. As shown in Table [1](#Tab1){ref-type="table"}, these two changes raise the thermodynamic efficiency from 67% to 77%, entirely due to the change in the force ratio. We conclude that the average thermodynamic efficiency of the actual proton pump elements of complex IV (red arrows in Fig. [4](#Fig4){ref-type="fig"}) is quite comparable to those of the other complexes, and that the lower overall efficiency of complex IV function is caused by the irreversible step of breaking the O--O bond that is not coupled to generation protonmotive force.

The degree of coupling {#Sec7}
----------------------

Table [1](#Tab1){ref-type="table"} shows the results for the degree of coupling (*q*) calculated for State 3 according to Eq. ([7](#Equ8){ref-type=""}). What is essential here is that the degree of coupling is very high for all three complexes of the respiratory chain; *q* is well over 0.99. Due to this, the normalised output/input force ratio (*Zx*) in State 3 yields a fair approximation of the thermodynamic efficiency (see Eqs. ([2](#Equ3){ref-type=""}) and ([5](#Equ6){ref-type=""})), which gives credence to the use of this ratio as an efficiency indicator (e.g., ref. ^[@CR30]^). The degree of coupling is fairly insensitive to the protonmotive force in State 3, which was assumed here to be 170 mV (see above). If the Δ*P* were higher, the degree of coupling would be even higher than in Table [1](#Tab1){ref-type="table"}. If it were 150 mV, the degree of coupling for succinate oxidation (complexes III + IV), for example, would only fall to 0.992.

The very high degree of coupling (*q* \> 0.99) found here from the data of phosphorylating (State 3) mitochondria might seem to be in sharp contrast to conclusions made earlier on the basis of the dependence of mitochondrial respiratory rates on membrane potential. For example, Hinkle et al.^[@CR23]^ found that the rate of State 4 respiration, measured at the same membrane potential as in the phosphorylating State 3, was a considerable fraction of the rate of respiration in State 3. The energy loss assessed this way was almost 4% with succinate as substrate (complex III + IV activity) and \~7% for complex IV alone (see also "Discussion"). However, the flux ratios (H^+^/e^−^) in State 3 observed here were 0.95--0.98 after normalisation to *Z* (Table [1](#Tab1){ref-type="table"}), suggesting a 2--5% loss to leaks which is consistent with the aforementioned data. It should be emphasised that as long as the degree of coupling is not ideal (i.e. *q* \< 1), the flux ratio is bound to decrease with increasing force ratio^[@CR14]--[@CR16]^. As already concluded above, *j* is equal to *qZ* only at *x* = 0, i.e., at level flow when the protonmotive force is zero.

The coupled reactions run close to equilibrium {#Sec8}
----------------------------------------------

Δ*G* denotes the free energy change of the coupled reaction of oxidoreduction-linked proton translocation for anyone of the three protonmotive complexes or their combination, i.e.$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{array}{*{20}{c}} {{\Delta}G} = { - {\Delta}E + Z{\Delta}P} \\ {} = {X_2 - ZX_1} \end{array}.$$\end{document}$$

The redox span for the different complexes of the respiratory chain (X~2,~ or Δ*E*) may be obtained from the literature as before (Table [2](#Tab2){ref-type="table"}^[@CR20],[@CR31]^). For State 3, we have assumed a pmf of 170 mV (see above), and this yields a free energy change of ca. −44 mV for complexes I and III, which corresponds to \~−1 kcal/mol, showing that even in State 3 these complexes operate quite close to thermodynamic equilibrium. Complex IV appears less reversible, but when the irreversible O--O bond splitting is subtracted and the O~2~ concentration adjusted to physiological levels (see above), the average free energy change of the reactions that actually drive transmembrane proton transfer is only of the order of −1.8 kcal/mol (−78 mV).Table 2Driving force of the redox (Δ*E*) and of the overall redox-coupled proton translocation (Δ*G*) reactions in States 3 and 4.State 3State 4Δ*E* (mV)Δ*G* (mV)Δ*E* (mV)Δ*G* (mV)*q*(*X*~2~)(*ZX*~1~ + *X*~2~)(*X*~2~)(*ZX*~1~ + *X*~2~)Complexes I + III + IV1042 ± 5 (*n* = 3)−1921057 ± 9 (*n* = 3)−1230.880Complexes III + IV621 ± 4 (*n* = 3)−111643 ± 4 (*n* = 3)−820.87Complex I384 ± 6 (*n* = 3)−44379 ± 4 (*n* = 3)−5\>0.99Complex III214 ± 9 (*n* = 4)−44222 ± 6 (*n* = 4)−350.84Complex IV418 ± 15 (*n* = 5)−78434 ± 22 (*n* = 5)−600.86(Complex IV)488 ± 15 (*n* = 5)−148504 ± 22 (*n* = 5)−131All numbers are in millivolts ± standard deviation. Δ*P* = 187 mV in State 4, 170 mV in State 3. The apparent *q* in State 4 is calculated according to Eq. ([9](#Equ10){ref-type=""}). Complex IV data are calculated in two ways (see text and legend to Table [1](#Tab1){ref-type="table"}). Values are expressed as mean ± SD, where *n* is the number of distinct samples.

The values for State 4 are also listed in Table [2](#Tab2){ref-type="table"}. In the absence of an experimental value for Δ*P* in the Muraoka & Slater experiments^[@CR20]^, we have estimated it as follows. We have collected the cases where respiration was studied with succinate and with TMPD + ascorbate as substrates^[@CR20]^, where there is no or very limited net electron flux across complex I. Instead, the generated Δ*P* drives electrons backwards from cytochrome *c* and ubiquinone to NAD^+^ until a steady state is reached, in which the redox reactions of complex I should equilibrate with the Δ*P*. Under those conditions, we found that Δ*E* across complex I was 373 and 374 mV, respectively, i.e. corresponding to a Δ*P* of 187 mV, which seems to be a very good approximate for isolated rat-liver mitochondria in the presence of 10 mM inorganic phosphate^[@CR7],[@CR9]^. Using this value for the Δ*P*, we find that complex I is indeed operating very close to thermodynamic equilibrium in State 4 during oxidation of NAD-linked substrates. By contrast, complexes III and IV deviate more from equilibrium (Table [2](#Tab2){ref-type="table"}; see below).

Is the degree of coupling a constant? {#Sec9}
-------------------------------------

State 4 (static head) is defined as the state where net generation of Δ*P* (*J*~*1*~) vanishes. If *j* is zero in Eq. ([2](#Equ3){ref-type=""}), then the numerator of the right hand side must be zero, from which follows that the force ratio, *x* = Δ*P*/Δ*E*, is a direct function of the degree of coupling^[@CR14]--[@CR16]^, such that$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q = Zx.$$\end{document}$$

This is a simplification because any proton leak across the membrane^[@CR7]--[@CR10]^ will ensure that net generation of Δ*P* (*J*~*1*~) by the pumping complexes still occurs at a rate determined by this leak. However, in well-coupled mitochondria the leak is slow relative to the forward and backward fluxes in the respiratory chain, so that the redox reactions of the respiratory complexes remain fairly close to equilibrium with the protonmotive force in State 4, and even in State 3 (Table [2](#Tab2){ref-type="table"}).

Applying Eq. ([9](#Equ10){ref-type=""}) to the State 4 data (Table [2](#Tab2){ref-type="table"}) suggests that *q* is 0.99 for complex I, but only 0.84 and 0.90 for complexes III and IV, respectively, even though the latter is adjusted for the irreversibility of O--O bond splitting and assessed at 12 μM O~2~, as described above. Values of *q* significantly smaller than 0.99 have been reported before for isolated mitochondria^[@CR15]^, but they have also been based on Eq. ([9](#Equ10){ref-type=""}) and the State 4 situation.

This result suggests a close approach to thermodynamic equilibrium between the oxidoreduction reaction of complex I and its associated proton translocation in State 4, and thus that the degree of coupling is close to unity as it was found to be in State 3 (Table [1](#Tab1){ref-type="table"}). By contrast, the degree of coupling for complexes III and IV appear to be clearly lowered from the value found in State 3, i.e., during active oxidative phoshorylation.

Discussion {#Sec10}
==========

With regard to the requirement of Onsager symmetry in Eq. (1) in order to apply the current treatment, which is based on the principles of near-equilibrium irreversible thermodynamics^[@CR14]^, such symmetry has been ensured experimentally for the reactions of oxidative phosphorylation^[@CR15],[@CR16],[@CR30],[@CR32]^. In addition, it may be pointed out that under physiological conditions complexes I, III and IV turn over at rates that are very slow (see ref. ^[@CR30]^) in comparison with the reaction rates of forward and backward electron and proton transfer within the three complexes. This means that the rates of physiological throughput are only a fraction of the individual forward and backward rates, again suggesting a system optimised to work close to thermodynamic equilibrium.

The thermodynamic efficiency of the three mitochondrial protonmotive redox complexes has not been systematically assessed thus far. Since thermodynamic efficiency is the product of the flux and force ratios of the output (H^+^ translocation) versus input (electron transfer) reactions (*η* = *−jx*; Eq. ([5](#Equ6){ref-type=""})), the flux ratio (*j*), which equals the H^+^/e^−^ ratio, is an essential parameter for such an assessment. A crucial obstacle has been that the thermodynamic efficiency approaches zero in the very experimental conditions strived at when measuring H^+^/e^−^ ratios, viz. conditions where the protonmotive force is zero, or very small ("level flow" condition). In the State 4 condition of maximal protonmotive force, there is an approach to zero net flux that also yields zero thermodynamic energy conservation efficiency. The individual pumping complexes themselves can operate at very high efficiency in State 4, but the overall efficiency of the electron-transport chain driving ATP synthesis is collapsed to zero due to the membrane leak (Fig. [2](#Fig2){ref-type="fig"}).

This problem can be approached by determination of the H^+^/e^−^ ratios of proton translocation by the mitochondrial respiratory chain in the phosphorylating State 3 from measured ATP/2e^−^ ratios, since the H^+^/ATP ratio of the F~1~F~o~ ATP synthase, and the translocation mechanisms of adenine nucleotides and inorganic phosphate are now known for animal mitochondria^[@CR22]^. The implicit assumption here is that in State 3, these processes are "stoichiometric", i.e., that there is no significant slippage in the mechanisms of the ATP synthase or the ADP/ATP and inorganic phosphate carriers. Thus the overall ratio of protons translocated per ATP synthesised extramitochondrially is assumed to be 3.67 (Eq. ([6](#Equ7){ref-type=""})). Some data^[@CR33]^ suggested that the actual H^+^/ATP ratio may be somewhat lower than the one predicted from the structure of the F~o~F~1~ complex, but a later very careful account has suggested that the ratio is precisely determined by the structure^[@CR34]^. If, nevertheless, the H^+^/ATP ratio were lower, the thermodynamic efficiencies reported here would be somewhat lower. More important, if that were the case, the estimates of the degree of coupling would have to be adjusted upwards, which seems unlikely considering that *q* is already higher than 0.99.

State 3 is indeed the only formally relevant mitochondrial condition for determining the thermodynamic efficiency of energy transduction, and was found here to be ca. 75--90% overall (Table [1](#Tab1){ref-type="table"}). As shown in Fig. [3](#Fig3){ref-type="fig"}, such high efficiencies combined with the observed range of force ratios require a degree of coupling of at least 0.99. Rocha and Springett^[@CR30]^ reported thermodynamic efficiencies for complexes I and III in intact isolated cells exceeding 90%, but these were based on the force ratio alone, corresponding to *Zx* in Table [1](#Tab1){ref-type="table"}. Hence they relied on the assumption that the flux ratio equals the mechanistic stoichiometry, and therefore that *q* = 1. As shown in Table [1](#Tab1){ref-type="table"}, this is a reasonable approximation for isolated mitochondria. The higher efficiency observed in cells compared with mitochondria could well be due to the fact that the cells were not studied under true State 3 conditions but under more natural conditions of energy load, which is expected to affect the efficiency^[@CR16]^.

Correcting the experimentally obtained H^+^/e^−^ values for the force ratio in State 3 made it possible to determine the degree of coupling (*q*), which was found to be extremely high for all three complexes. In conditions extrapolated to level flow where Δ*P* = 0, Eq. ([10](#Equ11){ref-type=""}) applies, viz.$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$j = - qZ.$$\end{document}$$

Hence, our finding that *q* ≥ 0.99 means that the flux ratio (*j*) at level flow is very closely equal to the mechanistic H^+^/e^−^ stoichiometry (*Z*), which in turn implies that despite the full redox-driving force, there can be little or no leakage due to "uncoupled" electron transfer ("electron leak").

These findings also corroborate earlier independent conclusions^[@CR8]--[@CR10],[@CR13],[@CR23]^ that slip and leak in the proton pumps of the respiratory chain do not occur significantly in mitochondria during oxidative phosphorylation. A slip in proton pumping is defined as a decrease of the H^+^/e^−^ stoichiometry either due to an electron transfer step being possible without linkage to protons ("electron slip", excluded above for State 3), and pump leak is defined as the possibility of an already translocated proton being able to leak back to its original location by a failure of the pump mechanism ("proton slip"), or by a proton leak of the membrane (see Fig. [2](#Fig2){ref-type="fig"}). The latter two cases of proton leak, driven by Δ*P*, were found to occur to a total maximum of ca. 5% of the total flux in State 3.

In State 4 mitochondria (static head in the nomenclature of the theory of irreversible thermodynamics), Δ*P* is maximised. It is often assumed that there is no net proton pumping and that a near-equilibrium state is reached between the redox and proton translocation reactions. Such a near-equilibrium in State 4 was already concluded by Muraoka and Slater in their seminal work from 1969^[@CR35]^. Although, in reality, a membrane proton leak balances redox-linked proton translocation by the three respiratory chain complexes, this leak is much slower than the forward and backward reactions of redox-linked proton translocation ensuring that a state near-equilibrium is reached. It is well known from the classical work by Nicholls^[@CR7]^, Hinkle^[@CR23]^ and Brand^[@CR8]--[@CR10]^ that when the protonmotive force increases beyond ca. 170 mV, the initially linear (Ohmic) dependence of membrane proton conductance on Δ*P* becomes exponential and the dependence is much steeper. This general proton leak (Fig. [2](#Fig2){ref-type="fig"}), the physiological importance of which has been stressed, has been ascribed to be an unspecific property of the mitochondrial solute carriers of the inner mitochondrial membrane, of which the ADP/ATP carrier is the most abundant^[@CR36]^.

It is therefore not surprising that when the degree of coupling (*q*) is assessed based on Eq. ([9](#Equ10){ref-type=""}), i.e. assuming that redox-linked proton translocation is zero (*j* = *0*), it is found to be lowered from the values observed in State 3 (see Table [2](#Tab2){ref-type="table"}). However, if the flux in State 4 is at least predominantly due to membrane proton leak, then the flux ratio (*j*) cannot be assumed to be zero, but is a function of *q*, *Z* and *x* according to Eq. ([2](#Equ3){ref-type=""}). Interestingly, in State 4 the driving force (Δ*G*) for complex I is very low (Table [2](#Tab2){ref-type="table"}), showing that the redox-driven proton translocation is very close to thermodynamic equilibrium (cf above). However, complexes III and IV appear to remain further from equilibrium (Table [2](#Tab2){ref-type="table"}). It thus seems that in State 4 whilst complex I approaches equilibrium closely, complexes III and IV remain further away. It may be especially revealing that the overall driving forces (ΔG) across complexes I and III are similar in State 3, but that upon the rise of Δ*P* to State 4 values, complex I comes very near to equilibrium, whereas complex III does not. We emphasise that the high absolute value of *q* that we find for complex I in State 4 (Table [2](#Tab2){ref-type="table"}) is not in discord with the fact that there is a significant membrane leak under such conditions^[@CR7]--[@CR10]^. The high State 4 value of *q* in Table [2](#Tab2){ref-type="table"} is due to the simplification by treating State 4 according to Eq. ([9](#Equ10){ref-type=""}), where it is assumed that net proton flux (*J*~*1*~ in Eq. ([1a](#Equ1){ref-type=""}), or *J*~*H*~^*F*^ in Fig. [2](#Fig2){ref-type="fig"}) is zero, when in actual fact it constitutes the sum of the membrane proton leak (*J*~*H*~^*M*^; Fig. [2](#Fig2){ref-type="fig"}) and possible slips in the pumps (*J*~*H*~^*L*^; Fig. [2](#Fig2){ref-type="fig"}).

It is the observation of a different behaviour of complexes III and IV relative to complex I in State 4 that is difficult to explain on the basis of a membrane proton leak, which should affect all three complexes similarly, and this suggests that slips between the oxidoreduction and proton translocation reactions may occur in complexes III and IV driven specifically by the high Δ*P* in that state.

Experiments by Wilson et al.^[@CR37]^ have shown that even a small lowering of the energy state (here considered to be the protonmotive force) can compensate for the lowered rate of respiration due to a decrease in oxygen tension. Thus local hypoxia will increase the level of reduced cytochrome *c*, which causes compensatory enhancement of the turnover of cytochrome *c* oxidase, and hence of the entire respiratory chain. This modulation will be all the more efficient if the linkage to proton pumping is loosened at high protonmotive force, which would otherwise prevent the compensatory increase in turnover. Interestingly, a model study by Krab et al.^[@CR38]^ also showed that a decrease in protonmotive force speeds up the respiratory rate at low O~2~ concentrations, which was observed as a decrease in the apparent K~M~ for O~2~. Clearly, the terminal member of the respiratory chain, complex IV, is a very sensitive detector of tissue anoxia. The departure from equilibrium and partial loss of coupling efficiency of complex IV in State 4 (Table [2](#Tab2){ref-type="table"}) is indicative of an important role in metabolic regulation^[@CR39]--[@CR41]^. The presence of "pulsed" and "resting" forms of complex IV^[@CR42],[@CR43]^ might well be examples of how its energy transduction efficiency may be modulated.

We conclude that the thermodynamic efficiency of the three protonmotive respiratory chain complexes is ca. 80% during oxidative phosphorylation in State 3 of isolated rat-liver mitochondria. Under the same conditions, the degree of coupling between redox activity and proton translocation is more than 99% for all three complexes. However, at a high enough ATP/ADP^.^P~i~ ratio, when net phosphorylation stops and the respiratory rate reaches a minimum (State 4), the degree of coupling is lowered significantly more for complexes III and IV than it is for complex I. Although it is clear that a membrane proton leak is also present in State 4, the different effect on complexes III and IV relative to complex I suggest significant "slipping" of the proton pump machineries of the two former complexes at high protonmotive force. Such slipping is likely to be of physiological importance for at least three reasons. First, it counteracts a blockade of respiratory flux at high energy conditions. Secondly, it minimises accumulation of highly reduced respiratory intermediates and hence their bimolecular reaction with oxygen, producing superoxide. Thirdly, it modulates the reaction of complex IV at low concentrations of oxygen in such a way as to raise the apparent oxygen affinity, an additional security measure that maintains respiratory flux in hypoxia. Finally, our observations suggest that metabolic control mechanisms of cell respiration should primarily be sought from the functions of complexes III and, in particular, complex IV.

Methods {#Sec11}
=======

Experimental data {#Sec12}
-----------------

The basic experimental data (redox states of NAD^+^, ubiquinone (Q) and cytochrome *c*) are derived from the work by Muraoka and Slater^[@CR20]^, who studied isolated rat-liver mitochondria respiring either on one of three NAD-linked substrates, β-hydroxybutyrate, pyruvate and glutamate, on succinate, or on TMPD + ascorbate, at pH = 7.4. The pH in the mitochondrial matrix was assumed to be 8^[@CR23]^. Under such conditions, the Δ*P* in State 4 is estimated to be ca. 187 mV (see "Results"), and is lowered to ca. 170 mV in State 3 (see refs. ^[@CR7],[@CR23]^). The redox potentials of the NAD^+^/NADH, ubiquinone/ubiquinol (Q/QH~2~), and ferric/ferrous cytochrome *c* couples were calculated from the measured extent of reduction^[@CR20]^ using the Nernst equation, and the appropriate midpoint potentials (*E*~m~) relative to the Normal hydrogen electrode (see also Supplementary Data). For complex I, we used the matrix pH of 8.0 to calculate the redox potentials of NAD^+^/NADH and Q/QH~2~ (*E*~m,8~ values of −350 mV and +30 mV, respectively). For complex III, the midpoint redox potential of the donor, Q/QH~2~, was taken to be an *E*~m,7.4~ of +66 mV (note, using the extramitochondrial pH), and for the acceptor (cytochromes *c* + *c*~*1*~) we used +230 mV due to the low ionic strength of the reaction medium.

For complex IV, we used the *E*~m,7.4~ = +230 mV for the donor, as before. For the acceptor, the O~2~/2H~2~O redox couple, we used the extramitochondrial pH (=7.4) as the reference point (see below). The *E*~m,7~ of the O~2~/2H~2~O couple is 815 mV at 1 atm fugacity and 25 °C, corresponding approximately to 1.2 mM of O~2~. In the conventional assessments of the redox potential across complex IV, we have assumed that \[O~2~\] = 120 μM, i.e., at about half-way of the oxygen concentration scale in experiments with isolated mitochondria suspended in air-equilibrated buffer at 25 °C. In the special assessments of the force ratio across the proton-pumping reactions of complex IV (marked separately in Tables [1](#Tab1){ref-type="table"} and [2](#Tab2){ref-type="table"}), we have adjusted the force ratio (and the other parameters) in two ways. First, we have taken into account that one of the partial catalytic reactions of complex IV catalysis is an "irreversible" scission of the O--O bond of O~2~ that is not coupled to the proton pump (Fig. [4](#Fig4){ref-type="fig"}). This roughly amounts to a 5 kcal/mol^[@CR25]^ (\~220 mV) subtraction of the free energy expenditure in the overall O~2~ reduction cycle. Secondly, as the redox potential of cytochrome *c* does not shift towards a more negative value until the oxygen concentration falls below ca. 10 μM^[@CR29]^, we have used an O~2~ concentration of 12 μM.

Thermodynamic H^+^/e^−^ stoichiometries {#Sec13}
---------------------------------------

The reason for using the mitochondrial matrix pH for the donor and acceptor potentials of complex I but the extramitochondrial pH for the potentials of complexes III and IV is explained below. We also address the often confusing thermodynamic H^+^/e^−^ stoichiometries of 2, 1 and 2 for complexes I, III and IV, when the number of protons released (per electron) on the cytoplasmic side of the membrane is different, viz. 2, 2 and 1, respectively. Note, however, that the number of protons taken up from the matrix side coincides with the thermodynamic stoichiometries. For complex I, one proton per electron is taken up from the matrix on reduction of ubiquinone, and one is released back to the matrix on oxidation of the NAD-dependent substrate. Hence, the net substrate proton "count" is zero, even though oxidoreduction of NADH is linked to only 0.5 H^+^/e^−^ at a pH near 7.

Redox energy is transduced into a pH gradient (ΔpH) when protons are taken up from a low proton concentration (matrix side, high pH) and released into a high proton concentration (cytosolic side, low pH). The amount of energy conserved is *n*~*H*~*×*ΔpH, where *n*~*H*~ is the number of protons taken up and released. In contrast, redox energy is transferred into ΔΨ when charge moves across the membrane. The energy conserved is *n*~*Q*~ *×* ΔΨ, where *n*~*Q*~ is the number of translocated charges. The protons taken up or released in the reduction or oxidation of NADH, ubiquinone, or oxygen are "substrate protons" that should basically not be included in the calculation of energy conservation, as they are already included in their redox potentials (but see below, and ref. ^[@CR30]^).

For a pure proton pump, such as complex I, the number of protons taken up and released is equal to the number of charges transferred across the membrane and the energy conserved is *n*~*P*~*(*ΔpH + ΔΨ*)* = *n*~*P*~Δ*P*, where *n*~*P*~ is the number of protons pumped. The Δ*G* (per electron) of the reaction catalysed by complex I (Δ*G*^*I*^*)* is thus:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\Delta}G^I = - \left( {E_h^{{\mathrm{UQ}}} - E_h^{{\mathrm{NADH}}}} \right) + 2{\Delta}P,$$\end{document}$$where *E*~*h*~^UQ^ and *E*~*h*~^NADH^ are the redox potentials of the ubiquinone and NADH pools using the pH in the matrix.

Complex IV pumps one proton per electron, which conserves one *Δ*P of redox energy. In addition, the electron starts on the cytosolic side and moves to the binuclear centre buried in the membrane, and a substrate proton moves from the matrix side to the binuclear centre to complete O~2~ reduction to water. The net effect is to move an additional positive charge across the entire membrane, conserving another ΔΨ of redox energy, giving a total energy conservation of Δ*P* + ΔΨ. The Δ*G* for complex IV is (Δ*G*^*IV*^) is then:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\Delta}G^{IV} = - \left( {E_h^{O,m} - E_h^{{\mathrm{Cyt}}c}} \right) + {\Delta}P + {\Delta}{\Psi},$$\end{document}$$where *E*~*h*~^*O,m*^ is the redox potential of the O~2~/H~2~O couple based on the matrix pH, and *E*~*h*~^Cyt*c*^ is the (pH-independent) redox potential of Cyt*c*. However, if the redox potential of the O~2~/H~2~O couple is instead calculated based on the cytosolic pH (*E*~*h*~^*O,c*^), Eq. ([12](#Equ13){ref-type=""}) becomes:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\Delta}G^{IV} = - \left( {E_h^{O,c} - E_h^{{\mathrm{Cyt}}c}} \right) + 2{\Delta}P.$$\end{document}$$

The situation for complex III is more complicated: for each electron transferred to cytochrome *c*, two protons are released from the Q~o~ centre to the cytosolic side and one proton is taken up from the matrix at the Q~i~ centre. However, this does not result in charge movement across the membrane because the Q~o~ centre is on the cytosolic and the Q~i~ centre on the matrix side of the membrane. Instead, one electron moves across the membrane from the Q~o~ to the Q~i~ centre (via the *b* cytochromes), which is equivalent to one positive charge moving from the matrix to the cytosolic side, and so conserves one ΔΨ of redox energy. The release of protons to the cytosolic side and uptake of protons from the matrix are not basically included in the energy conservation calculation because they are substrate protons of the oxidoreduction of ubiquinone, and hence included in the redox potential of UQ/UQH~2~. However, the redox potential of UQ/UQH~2~ at the Q~o~ site is different from that at the Q~i~ site for the same UQ/UQH~2~ ratio, because the sites use protons from different compartments at different concentrations. Mechanistically, one electron from the oxidation of UQH~2~ at Q~o~ reduces Cyt*c* and one reduces UQ at the Q~i~ site. The Δ*G* for complex III (Δ*G*^*III*^) is then$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\Delta}G^{III} = - \left( {E_h^{{\mathrm{Cyt}}c} + E_h^{{\mathrm{UQi}}} - 2E_h^{{\mathrm{UQo}}}} \right) + {\Delta}{\Psi},$$\end{document}$$where *E*~*h*~^UQi^ and *E*~*h*~^UQo^ are the redox potentials of the UQ/UQH~2~ pools at the Q~i~ and Q~o~ sites, respectively. The difference between *E*~*h*~^UQo^ and *E*~*h*~^UQi^ is ΔpH (in millivolts) (see supplementary data of ref. ^[@CR44]^), so that Eq. ([14](#Equ15){ref-type=""}) can be rewritten as:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\Delta}G^{III} = - \left( {E_h^{{\mathrm{Cyt}}c} - E_h^{{\mathrm{UQo}}}} \right) + {\Delta}P.$$\end{document}$$

Thus, as shown by Eqs. ([13](#Equ14){ref-type=""}) and ([15](#Equ16){ref-type=""}), using the pH of the cytosol to calculate the Δ*G* of complexes III and IV simplifies the thermodynamic analysis.

Statistics and reproducibility {#Sec14}
------------------------------

The P/O ratios required to calculate flux ratios in State 3 (Eq. ([6](#Equ7){ref-type=""}); Table [1](#Tab1){ref-type="table"}) were mean values taken from ref. ^[@CR21]^, omitting the data published before 1960 due to its very high individual variability (Supplementary Tables [1](#MOESM1){ref-type="media"} and [2](#MOESM1){ref-type="media"}). Thus *n*, the number of measurements, was 9, 14 and 5, respectively, for the spans of complexes I + III + IV, III + IV and IV. The redox spans required to determine force ratios (*x*) in Tables [1](#Tab1){ref-type="table"} and [2](#Tab2){ref-type="table"} were obtained from ref. ^[@CR20]^. All the available data were used, i.e., respiration on β-hydroxybutyrate, pyruvate and glutamate (complexes I + III + IV), succinate (complexes III + IV) and ascorbate + TMPD (N,N,N',N'-tetramethyl-p-phenylene diamine; complex IV). The results were taken as the mean ± SD with *n* values of 3 for complex I, 4 for complex III and 5 for complex IV (Supplementary Table [3](#MOESM1){ref-type="media"}).

Reporting summary {#Sec15}
-----------------

Further information on research design is available in the [Nature Research Reporting Summary](#MOESM4){ref-type="media"} linked to this article.

Supplementary information
=========================

 {#Sec16}

###### 

Supplementary Information

###### 

Supplementary Data 1

###### 

Description of Additional Supplementary Files

###### 

Reporting Summary

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Mårten Wikström, Roger Springett.

Supplementary information
=========================

**Supplementary information** is available for this paper at 10.1038/s42003-020-01192-w.

M.W. would like to thank Martin D. Brand, S. Roy Caplan, Antony R. Crofts, Gerhard Hummer, Artur Osyczka, Hagai Rottenberg and John E. Walker for helpful comments. Prof. Olli Silvennoinen (Helsinki Institute of Life Science) is acknowledged for financial support.

The initial concept was advanced by M.W.

All data supporting the findings of this study are available within the paper and its Supplementary Files or available from the corresponding author upon reasonable request. Oxidation states of NAD^+^, UQ and Cyt*c* and P/O ratios that support the findings of this paper were taken from Muraoka and Slater^[@CR20]^ and Hinkle^[@CR21]^, respectively, and reproduced in Supplementary Tables [1](#MOESM1){ref-type="media"}--[3](#MOESM1){ref-type="media"}. The redox potentials and thermodynamic parameters derived from these data are given in the Supplementary Data file.

The authors declare no competing interests.
