Module 3: Oxidative Phosphorylation
The respiratory chain couples the oxidation of NADH and FADH2 to molecular oxygen to the synthesis of ATP—not directly, but through an electrochemical proton gradient across the inner mitochondrial membrane. Peter Mitchell’s chemiosmotic hypothesis (1961, Nature) revolutionized bioenergetics by explaining this coupling in thermodynamic rather than chemical-intermediate terms, winning the 1978 Nobel Prize. Four multi-subunit enzyme complexes (I–IV), plus two mobile carriers (ubiquinone and cytochrome c), assemble into higher-order respirasomes that achieve electron transfer with ~40–50% thermodynamic efficiency.
1. Mitchell 1961: The Chemiosmotic Hypothesis
Before 1961, bioenergetics was dominated by the search for a “high-energy intermediate” X∼I that would link oxidation to phosphorylation. This mythical molecule was never found. Peter Mitchell, working almost alone at Glynn House in Cornwall, proposed that the intermediate was not a chemical bond but a vectorial proton gradient maintained across a topologically closed membrane (Mitchell 1961, Nature, 191, 144–148).
The central claims of the chemiosmotic hypothesis are:
- Electron-transport complexes are vectorial pumps that translocate protons from the matrix to the intermembrane space (IMS).
- The inner mitochondrial membrane is impermeable to protons, so the gradient persists.
- ATP synthase (F0F1) is a reversible molecular motor that uses the gradient to synthesize ATP.
- Coupling is obligatory but indirect: any process that discharges the gradient (uncoupling, leak, ionophores) blocks ATP synthesis without blocking electron transport.
The hypothesis was met with intense skepticism. Paul Boyer and others argued for a conformational coupling model. Experimental vindication came gradually through the 1960s and 1970s: reconstitution of purified bacteriorhodopsin into vesicles with F0F1(Racker & Stoeckenius 1974) showed that a light-driven proton pump alone was sufficient to drive ATP synthesis—proof that the “intermediate” was indeed just a proton gradient. Mitchell received the Nobel Prize in Chemistry in 1978.
\[ \Delta p \;=\; \Delta\psi \;-\; \frac{2.303\,RT}{F}\,\Delta\text{pH} \]
The total proton-motive force \(\Delta p\) (in mV) decomposes into an electrical component \(\Delta\psi\) and a chemical (pH) component. At 37 °C, \(2.303\,RT/F \approx 61.5\) mV per pH unit.
For a typical mammalian mitochondrion with \(\Delta\psi \approx 140\) mV (negative inside the matrix) and \(\Delta\text{pH} \approx 1.0\)(alkaline inside), the total PMF is\(\Delta p \approx 140 + 61.5 \times 1.0 \approx 200\) mV. In plants and some bacteria the balance is different: chloroplasts at pH 8 matrix vs. pH 5 stroma have almost all of their PMF in the \(\Delta\text{pH}\) term.
2. Thermodynamics of the Proton Gradient
The free-energy change for moving one proton from the IMS (high [H+], low pH) to the matrix (low [H+], high pH) is
\[ \Delta G_{H^+} \;=\; -F\,\Delta p \;=\; -F\,\Delta\psi \;+\; 2.303\,RT\,\Delta\text{pH} \]
With \(F = 96{,}485\) C/mol and \(\Delta p = 0.2\) V, each proton stores \(\sim 19.3\) kJ/mol of free energy. If ATP synthase couples ~3 H+ to one ATP (Module 4 will refine to the c-ring stoichiometry), the input is \(\sim 58\) kJ/mol—comfortably exceeding the\(\Delta G_{\text{ATP}} \approx +50\) kJ/mol required under cellular conditions ([ATP]/[ADP][Pi] far from equilibrium).
Full oxidation of one NADH by O2 releases\(\Delta G^\circ = -220\) kJ/mol. The ETC pumps approximately 10 H+ per NADH (4 at C1, 2 at C3, 4 at C4), capturing\(10 \times 19.3 \approx 193\) kJ/mol in the gradient, for a thermodynamic efficiency of \(\sim 88\%\). Downstream, ATP synthase converts the gradient to ATP at \(\sim 65\%\) efficiency, yielding an overall oxidation-to-ATP efficiency of ~40%.
Bioenergetic partitioning of the PMF:
- ATP synthesis (~65%) at F0F1
- Pi/H+ symport and ATP/ADP antiport (adenine nucleotide translocase)
- Ca2+ uptake via MCU (~1 Ca2+ per 2 H+)
- Protein import via TOM/TIM (requires \(\Delta\psi\))
- Proton leak (10–25% under basal conditions)—heat generation
3. Complex I: NADH:Ubiquinone Oxidoreductase
Complex I is the entry point for electrons from NADH into the chain and the largest respiratory enzyme. Mammalian C1 contains 45 subunits in an L-shape with a hydrophilic arm protruding into the matrix (which hosts FMN and 8 iron-sulfur clusters) and a membrane arm embedded in the IMM (which contains the proton translocation machinery).
The reaction:
\[ \text{NADH} + \text{Q} + 5\,\text{H}^+_{\text{matrix}} \;\to\; \text{NAD}^+ + \text{QH}_2 + 4\,\text{H}^+_{\text{IMS}} \]
Two electrons flow from NADH through FMN and the iron-sulfur chain to ubiquinone. 4 H+ are pumped per 2 e-.
The electron path: NADH donates a hydride to FMN (flavin mononucleotide), producing FMNH2. Electrons then tunnel through eight Fe-S clusters (N1a, N3, N1b, N4, N5, N6a, N6b, N2) spaced <14 Å apart (Moser-Dutton criterion for biological tunneling)—a physiological electron wire ~90 Å long from FMN to the Q-binding site. The final cluster N2 transfers electrons to bound quinone, which after two-electron reduction and protonation leaves as QH2.
The proton pumping mechanism is uniquely long-range: four antiporter-like subunits (ND2, ND4, ND5, and the swapped-symmetry fourth channel in ND6/ND4L) are coupled to the Q-binding site by a ~100 Å lateral helix (HL) on the matrix side of ND5. This helix acts as a mechanical piston: Q reduction at the hydrophilic arm triggers conformational waves that propagate the length of the membrane arm, driving sequential protonation/deprotonation in each of the four half-channels. The mechanism was first proposed from the bacterial cryo-EM structure by Baradaran et al. (2013, Nature), refined in ovine C1 by Zhu et al. (2016, Nature), and in active/deactive conformations by Agip et al. (2020, Nat. Struct. Mol. Biol.).
Seven of the 45 subunits are mitochondrially encoded (ND1–ND6, ND4L) and constitute the proton-pumping core. This explains why Complex I is the most common site of inherited mitochondrial disease (LHON, MELAS-associated mutations, Leigh syndrome). The inhibitor rotenone binds the Q-binding pocket and is used experimentally to dissect electron flow; it was historically a fish poison used by indigenous peoples.
Complex I also has an active/deactive (A/D) transition: when oxygen is absent, the enzyme adopts an inactive “D” state that is slow to reactivate. On reperfusion, a large fraction of C1 is blocked in D-form, producing a lag that contributes to ischemia-reperfusion injury (Chouchani et al. 2013). During the reactivation phase, C1 runs in reverse as described below.
Reverse electron transport (RET):
At high PMF and a reduced Q-pool (typical after succinate oxidation with blocked C3/C4), electrons run backward from QH2 to NAD+through C1, producing superoxide at the flavin site. RET-generated ROS is a major contributor to oxidative damage in ischemia-reperfusion (Chouchani 2014, Nature—covered in Module 5) and is exceptionally sensitive to succinate accumulation.
4. Complex II: Succinate Dehydrogenase
Complex II is the respiratory chain’s simplest enzyme and the only one that participates dually in the TCA cycle. It catalyzes the oxidation of succinate to fumarate and donates the resulting electrons directly to the ubiquinone pool. Unlike C1, C3, and C4, C2 does not pump protons—a thermodynamic consequence of the small redox drop from FADH2to ubiquinone (\(\Delta E \approx 30\) mV).
The mammalian enzyme has four subunits:
- SDHA: catalytic flavoprotein, carries FAD covalently linked to histidine
- SDHB: iron-sulfur subunit with three Fe-S clusters [2Fe-2S], [4Fe-4S], [3Fe-4S]
- SDHC, SDHD: membrane anchors, cytochrome b560 heme, Q-binding pocket
All four subunits are nuclear-encoded. Mutations in SDHB, SDHC, SDHD cause hereditary paraganglioma and pheochromocytoma (Baysal et al. 2000, Science). Loss of SDH activity causes succinate accumulation, which stabilizes HIF-1α by inhibiting prolyl hydroxylases—linking mitochondrial dysfunction to oncometabolism.
Malonate is the canonical competitive inhibitor of succinate dehydrogenase (structural analog of succinate). It has re-emerged as a therapeutic candidate in ischemia-reperfusion: preemptive malonate blocks the succinate buildup that would otherwise drive RET-dependent ROS (Valls-Lacalle et al. 2016).
\[ \text{succinate} + \text{Q} \;\to\; \text{fumarate} + \text{QH}_2 \qquad (\text{no H}^+\text{ pumping}) \]
5. Complex III: The Q-Cycle
Complex III (cytochrome bc1) transfers electrons from QH2 to cytochrome c and translocates 2 H+/2 e-. It is a dimer of 11 subunits (mammalian) with three redox-active centers per monomer: cytochrome b (two hemes: low-potential bL and high-potential bH), cytochrome c1 (one heme c type), and the Rieske iron-sulfur protein(ISP, [2Fe-2S] cluster with coordinating histidines).
The mechanism is the famous Q-cycle proposed by Mitchell himself (1975, FEBS Lett.) and later refined with cryo-EM. The electrons from QH2 bifurcate at the outer Q-binding pocket (Qo site, near the IMS):
- First electron: from QH2 at Qo to Rieske ISP, then to cytochrome c1, then to soluble cytochrome c. Two protons released into IMS.
- Second electron: from the semiquinone Q•- intermediate at Qo to heme bL, tunneling through the membrane to bH, reducing a second Q at the inner Qi site to semiquinone.
- A second round repeats and fully reduces the Q at Qi to QH2. Two protons taken up from matrix.
\[ \text{QH}_2 + 2\,\text{cyt}\,c^{3+} + 2\,\text{H}^+_\text{matrix} \;\to\; \text{Q} + 2\,\text{cyt}\,c^{2+} + 4\,\text{H}^+_\text{IMS} \]
Net: 2 H+ translocated per 2 electrons. Two additional protons appear in IMS from QH2 oxidation—but these are “scalar” and counted as part of the net arithmetic.
The Rieske ISP executes a dramatic ~60° mechanical rotation during catalysis, carrying its Fe-S cluster between contact with QH2 at Qo and delivery to cytochrome c1. This “headpiece motion” was first revealed in yeast by Iwata et al. (1998, Science) and is the slowest step in the cycle.
Two classes of C3 inhibitor illustrate the bifurcation:
- Myxothiazol, stigmatellin: block Qo site, prevent first-electron transfer
- Antimycin A: blocks Qi site, traps semiquinone at Qo → ROS generation
The antimycin-induced semiquinone at Qo is the major physiological source of superoxide from C3 (covered in Module 5). The site is exposed to both leaflets of the IMM, so ROS appear on both sides—an important distinction from C1 which releases superoxide only into the matrix.
6. Complex IV: Cytochrome c Oxidase
Complex IV (cytochrome c oxidase, COX) is the terminal electron acceptor of the chain. It accepts electrons from cytochrome c and uses them to reduce molecular oxygen to water. The mammalian enzyme has 13 subunits; the three catalytic subunits (COX1, COX2, COX3) are mitochondrially encoded, while the remaining ten are nuclear.
The catalytic center is a binuclear heme a3-CuB site buried in the COX1 subunit. Heme a3and CuB are only ~4.5 Å apart and coordinate O2 as a peroxide bridge during turnover. Electrons enter from cytochrome c at a second copper site, CuA (dinuclear, in COX2), pass to heme a (in COX1), and then to the binuclear center for O2 reduction.
The overall reaction (per 4 electrons, 1 O2):
\[ 4\,\text{cyt}\,c^{2+} + \text{O}_2 + 8\,\text{H}^+_\text{matrix} \;\to\; 4\,\text{cyt}\,c^{3+} + 2\,\text{H}_2\text{O} + 4\,\text{H}^+_\text{IMS} \]
4 H+ pumped per 4 e- (= 2 H+/2 e-); an additional 4 “substrate” protons are consumed on the matrix side to make water, effectively adding to the charge separation.
The reaction cycle passes through four spectroscopically distinct intermediates (R, A, PR, F, OH) that cleave the O=O bond without releasing harmful intermediates—a triumph of enzyme design, since free hydroxyl or superoxide released here would wreak oxidative havoc. The structure of bovine COX was solved at atomic resolution by Tsukihara et al. (1996, Science).
Cyanide, azide, and carbon monoxide are classic COX inhibitors. Cyanide binds heme a3 as an axial ligand and brings respiration to a near-immediate halt, causing rapid death by energy starvation. COX has extraordinarily high O2 affinity (\(K_m \approx 0.1\) µM), which is why tissue pO2 can fall to extremely low values before respiration limits.
7. Mobile Carriers: Ubiquinone and Cytochrome c
Electrons cross from complex to complex via two small mobile carriers:
- Ubiquinone (Q): a hydrophobic isoprenoid quinone (coenzyme Q10 in humans) that diffuses within the lipid bilayer. Q/QH2 redox potential \(\approx +60\) mV; 50–100 Q per complex in the pool. As a two-electron carrier, it shuttles 2 e- + 2 H+per turnover, translating reducing equivalents efficiently across the membrane.
- Cytochrome c: a small (12 kDa, 104 residues) water-soluble heme protein that diffuses in the IMS along the outer surface of the IMM. Redox potential \(\approx +250\) mV; carries one electron per turnover. Also plays a second role: when released to the cytosol during apoptosis, it triggers caspase-9 activation and programmed cell death.
Cytochrome c is covalently tethered to its heme via two thioether bonds to cysteines in a conserved CXXCH motif—a rare modification found in few other proteins. The bovine holo-cytochrome c structure (Bushnell, Louie & Brayer 1990) shows a compact fold that presents a hydrophobic patch for docking to cytochrome c1 and COX via surface electrostatic complementarity.
Fun facts about cytochrome c:
- One of the first proteins sequenced (Margoliash 1962)—a classical molecular-evolution tool.
- Deeply conserved: human and yeast cytochromes c share ~50% identity and are functionally interchangeable.
- Released through BAX/BAK pores in the outer membrane in intrinsic apoptosis.
- Peroxidase activity when bound to cardiolipin—generates oxidized lipids that release it in apoptosis (Kagan et al. 2005).
8. Supercomplexes and the Respirasome
For decades the “fluid mosaic” view held that complexes I–IV diffuse independently in the IMM and encounter substrates by random Q/c collision. Blue-native PAGE (Schägger & Pfeiffer 2000, EMBO J.) changed this: mild solubilization preserves intact supercomplexeswith defined stoichiometries:
- I1III2: minimal respirasome core
- I1III2IV1: the canonical mammalian respirasome
- I1III2IV2: megacomplex in some tissues
- III2IV1-2: supercomplex without C1 (bacteria, plants)
- ATP synthase dimers: V-shape rows along cristae ridges
Cryo-EM structures of the mammalian respirasome by Letts et al. (2016, Nature), Sousa et al. (2016), and Gu et al. (2016) revealed the architecture: C1’s membrane arm docks against the C3 dimer, which in turn contacts C4. The assembly creates a protected micro-environment where Q and cyt c can diffuse short distances without being “lost” to the bulk membrane.
The adaptive value of supercomplex assembly is debated. Substrate channeling (Lapuente-Brun et al. 2013, Science) would speed electron transfer and reduce ROS escape. Alternative views emphasize structural stabilization without kinetic channeling (Blaza et al. 2014). A compromise view: supercomplexes partition the Q pool into accessible and sequestered fractions that respond differently to substrate and inhibitor challenges.
The supercomplex assembly factor SCAFI(COX7A2L) has been implicated in respirasome biogenesis; its expression varies by mouse strain and tissue, and correlates with the fraction of C1 in supercomplex form. This introduces an intriguing form of respiratory-capacity regulation at the organelle level.
9. Arithmetic of Respiration: Protons and the P/O Ratio
Accounting for every proton in and out yields the P/O ratio: moles of ATP synthesized per gram-atom of oxygen consumed. Classical biochemistry texts used 3 for NADH and 2 for FADH2, but careful modern measurements give smaller values:
Proton stoichiometries per 2 e-:
- Complex I: 4 H+ pumped
- Complex II: 0 H+
- Complex III: 2 H+ pumped (+2 scalar)
- Complex IV: 2 H+ pumped (+2 substrate)
Total pumped H+ per NADH-to-O2 oxidation:\(\approx 10\). Per FADH2-to-O2:\(\approx 6\).
Now the ATP synthesis side (details in Module 4): the c-ring of mammalian F0has 8 c-subunits, each carrying one H+; three ATPs are made per full revolution (120° per ATP, three beta catalytic sites). This gives\(8/3 \approx 2.67\) H+/ATP. Adding the cost of Pi/H+ symport into the matrix (+1 H+ per ATP) and ATP/ADP antiport (neutral but draws on \(\Delta\psi\)), the physiological total is ~4 H+ per cytosolic ATP.
\[ \text{P/O (NADH)} \;=\; \frac{10\ \text{H}^+}{4\ \text{H}^+/\text{ATP}} \;\approx\; 2.5 \]
\[ \text{P/O (FADH}_2\text{)} \;=\; \frac{6\ \text{H}^+}{4\ \text{H}^+/\text{ATP}} \;\approx\; 1.5 \]
Hinkle (2005, BBA) re-analyzed the experimental literature and recommends P/O = 2.5 (NADH) and 1.5 (FADH2) as best-current values.
Over a human lifetime at basal metabolic rate, ATP synthase spins at roughly 100 revolutions/s per enzyme, and the body synthesizes roughly its own weight in ATP each day (Module 4 will quantify this).
10. Pharmacology of the Respiratory Chain
The respiratory chain is the target of many natural products and drugs:
Site-selective inhibitors:
- Rotenone (C1, Q-binding pocket): fish poison, Parkinson’s animal model
- Metformin (C1, mild inhibitor; type 2 diabetes)
- Malonate (C2, competitive with succinate); ischemia-reperfusion therapy (Valls-Lacalle 2016)
- Myxothiazol / stigmatellin (C3, Qo site)
- Antimycin A (C3, Qi site) — classical ROS-inducer
- Cyanide, azide (C4, heme a3)
- Oligomycin (F0 c-ring, blocks ATP synthase proton channel)
Uncouplers (dissipate PMF):
- 2,4-Dinitrophenol (DNP): weight-loss drug in 1930s, removed after hyperthermia deaths
- FCCP, CCCP: research-use protonophores
- UCP1 (thermogenin): physiological uncoupler in brown adipose tissue; activated by free fatty acids during cold exposure
- Salicylate, aspirin (high dose): mild uncoupling partly explains toxicity
Uncouplers illustrate Mitchell’s argument beautifully: in their presence, the ETC keeps running (electrons still flow to O2), the PMF collapses, and ATP synthesis stops. Oxygen consumption actually increases (“state 3u”) because the brake of thermodynamic back-pressure is released.
11. Schematic: The Respiratory Chain
Simulation 1: Multi-Complex ETC Kinetics with Inhibitor Titration
Kinetic simulation of the respiratory chain modeled as a pool-balance system (following Beard 2005 and Wu et al. 2007). We integrate the ODEs governing NADH, QH2, cytochrome c redox, and PMF, with inhibitor parameters for rotenone (C1), malonate (C2), antimycin A (C3), and cyanide (C4). The output shows the site-specific effects: rotenone collapses PMF only when succinate is unavailable (succinate bypass via C2); antimycin dramatically over-reduces the Q pool; cyanide halts downstream electron flow at the terminal step. The rotenone dose-response illustrates non-linear control by C1.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Simulation 2: PMF Buildup and Uncoupler Titration
We simulate the dynamics of PMF formation after substrate addition to de-energized mitochondria, decomposing the total \(\Delta p\) into its\(\Delta\psi\) and \(\Delta\text{pH}\) components. Uncouplers (DNP, FCCP, CCCP) are modeled as proton leaks that short-circuit the gradient. The third panel shows the classical uncoupler dose-response on steady-state PMF and ATP synthesis, illustrating the catastrophic collapse of ATP flux when \(\Delta p\)falls below the ~150 mV threshold for ATP synthase. The fourth panel simulates UCP1-mediated thermogenesis in brown adipose tissue, where intentional uncoupling generates heat for thermoregulation during cold exposure.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Key References
• Mitchell, P. (1961). “Coupling of phosphorylation to electron and hydrogen transfer by a chemi-osmotic type of mechanism.” Nature, 191, 144–148.
• Mitchell, P. (1975). “The protonmotive Q cycle: a general formulation.” FEBS Lett., 59, 137–139.
• Racker, E. & Stoeckenius, W. (1974). “Reconstitution of purple membrane vesicles catalyzing light-driven proton uptake and adenosine triphosphate formation.” J. Biol. Chem., 249, 662–663.
• Baradaran, R., Berrisford, J.M., Minhas, G.S., & Sazanov, L.A. (2013). “Crystal structure of the entire respiratory complex I.” Nature, 494, 443–448.
• Zhu, J., Vinothkumar, K.R., & Hirst, J. (2016). “Structure of mammalian respiratory complex I.” Nature, 536, 354–358.
• Agip, A.-N.A. et al. (2018 & 2020). “Cryo-EM structures of complex I from mouse heart mitochondria in two biochemically defined states.” Nat. Struct. Mol. Biol., 25, 548–556.
• Iwata, S., Lee, J.W., Okada, K. et al. (1998). “Complete structure of the 11-subunit bovine mitochondrial cytochrome bc1 complex.” Science, 281, 64–71.
• Tsukihara, T. et al. (1996). “The whole structure of the 13-subunit oxidized cytochrome c oxidase at 2.8 Å.” Science, 272, 1136–1144.
• Schägger, H. & Pfeiffer, K. (2000). “Supercomplexes in the respiratory chains of yeast and mammalian mitochondria.” EMBO J., 19, 1777–1783.
• Letts, J.A., Fiedorczuk, K., & Sazanov, L.A. (2016). “The architecture of respiratory supercomplexes.” Nature, 537, 644–648.
• Gu, J. et al. (2016). “The architecture of the mammalian respirasome.” Nature, 537, 639–643.
• Sousa, J.S., Mills, D.J., Vonck, J., & Kühlbrandt, W. (2016). “Functional asymmetry and electron flow in the bovine respirasome.” eLife, 5, e21290.
• Lapuente-Brun, E. et al. (2013). “Supercomplex assembly determines electron flux in the mitochondrial electron transport chain.” Science, 340, 1567–1570.
• Blaza, J.N., Serreli, R., Jones, A.J.Y., Mohammed, K., & Hirst, J. (2014). “Kinetic evidence against partitioning of the ubiquinone pool.” PNAS, 111, 15735–15740.
• Hinkle, P.C. (2005). “P/O ratios of mitochondrial oxidative phosphorylation.” Biochim. Biophys. Acta, 1706, 1–11.
• Baysal, B.E. et al. (2000). “Mutations in SDHD, a mitochondrial complex II gene, in hereditary paraganglioma.” Science, 287, 848–851.
• Valls-Lacalle, L. et al. (2016). “Succinate dehydrogenase inhibition with malonate during reperfusion reduces infarct size.” Cardiovasc. Res., 109, 374–384.
• Nicholls, D.G. & Ferguson, S.J. (2013). Bioenergetics 4. Academic Press.