🔬

Comprehensive Course

Cell Physiology

A comprehensive exploration of cellular function, from membrane transport and energy metabolism to signal transduction and specialized cell types. Understand how cells maintain homeostasis, communicate, and perform their specialized functions.

📚8 Parts

📖40+ Chapters

🎯Medical & Graduate Level

Start Learning→

🎯Key Learning Objectives

🔲

Membrane Transport

Master the mechanisms of molecular transport across cell membranes

⚡

Bioelectricity

Understand membrane potentials, action potentials, and electrical signaling

📡

Signal Transduction

Learn how cells receive and process extracellular signals

🔄

Homeostasis

Explore how cells maintain internal balance and respond to challenges

📚Course Content

🔲

Part 1

Cell Membrane and Transport

Membrane structure, lipid bilayers, transport proteins, and movement of molecules across membranes

Membrane StructurePassive TransportActive TransportVesicular TransportMembrane Potential

⚡

Part 2

Cellular Energetics

ATP production, mitochondrial function, glycolysis, oxidative phosphorylation, and metabolic regulation

GlycolysisCitric Acid CycleElectron Transport ChainATP SynthaseMetabolic Integration

📡

Part 3

Cell Signaling

Signal transduction pathways, receptors, second messengers, and cellular communication

Receptor TypesG-Protein SignalingTyrosine KinasesSecond MessengersSignal Integration

💪

Part 4

Muscle Physiology

Muscle contraction mechanisms, excitation-contraction coupling, and muscle fiber types

Skeletal MuscleCardiac MuscleSmooth MuscleExcitation-Contraction CouplingMuscle Metabolism

🧠

Part 5

Neurophysiology

Action potentials, synaptic transmission, neurotransmitters, and neural circuits

Resting PotentialAction PotentialsSynaptic TransmissionNeurotransmittersNeural Integration

🏛️

Part 6

Epithelial Transport

Transepithelial transport, ion channels, water transport, and barrier functions

Tight JunctionsTranscellular TransportParacellular TransportWater ChannelsSecretion

🔔

Part 7

Calcium Signaling

Calcium as a second messenger, calcium channels, pumps, and calcium-dependent processes

Calcium ChannelsCalcium ReleaseCalcium PumpsCalcium Binding ProteinsCalcium Oscillations

📊

Part 8

Cell Volume Regulation

Osmotic balance, volume-sensitive channels, regulatory mechanisms, and cellular homeostasis

OsmosisVolume SensorsRVD and RVIOrganic OsmolytesApoptotic Volume Decrease

🎥Harvard Online: Cell Biology

Introductory lectures from Harvard Online covering fundamental cell structure, mitochondrial evolution, and the endosymbiotic origin of eukaryotic organelles.

🎥 Overview of Cell Structure

Harvard Online — Cell Biology — Tour of eukaryotic cell architecture: nucleus, organelles, cytoskeleton, and membrane systems

🎥 Cooperation and Evolution

Harvard Online — Cell Biology — Endosymbiotic theory, mitochondrial origins, and the evolution of eukaryotic complexity

📐Foundational Equations

Nernst Equation

$$E_X = \frac{RT}{zF}\ln\frac{[X]_o}{[X]_i}$$

Equilibrium potential for an ion species X

Goldman-Hodgkin-Katz

$$V_m = \frac{RT}{F}\ln\frac{P_{\text{K}}[\text{K}^+]_o + P_{\text{Na}}[\text{Na}^+]_o + P_{\text{Cl}}[\text{Cl}^-]_i}{P_{\text{K}}[\text{K}^+]_i + P_{\text{Na}}[\text{Na}^+]_i + P_{\text{Cl}}[\text{Cl}^-]_o}$$

Resting membrane potential from multiple ion permeabilities

Fick's Law of Diffusion

$$J = -D\frac{dC}{dx}$$

Diffusive flux proportional to concentration gradient

Michaelis-Menten

$$v = \frac{V_{\max}[S]}{K_m + [S]}$$

Enzyme kinetics and saturable transport

Cable Equation

$$\lambda^2\frac{\partial^2 V}{\partial x^2} = \tau\frac{\partial V}{\partial t} + V$$

Passive electrical spread along a neurite

Hill Equation

$$Y = \frac{[L]^n}{K_d^{\,n} + [L]^n}$$

Cooperative ligand binding with Hill coefficient n

🎥MIT OpenCourseWare: Membrane Biophysics

🎥 MIT 9.40: Nernst Potential & RC Circuit

MIT OCW — Introduction to Neural Computation (Spring 2018) — Mathematical derivation of Nernst equation, RC model of the neuron membrane

🎥 MIT 7.016: Neurons & Action Potentials

MIT OCW — Introductory Biology (Fall 2018) — Ion channels, membrane potential, action potential propagation, and optogenetics

⚡Metabolic Cycles: A Comprehensive Course

Metabolism is the sum of all chemical reactions in a cell. This section provides a detailed, university-level treatment of the six major metabolic pathways, including reaction mechanisms, enzyme regulation, energetics, and clinical connections.

Glycolysis

10 Steps

Overview

Glycolysis (from Greek glykys = sweet, lysis = splitting) is the universal pathway for glucose catabolism, occurring in the cytoplasm of virtually all living cells. It converts one molecule of glucose (6 carbons) into two molecules of pyruvate (3 carbons each), with a net production of 2 ATP and 2 NADH.

Overall Reaction

$$\text{Glucose} + 2\,\text{NAD}^+ + 2\,\text{ADP} + 2\,P_i \longrightarrow 2\,\text{Pyruvate} + 2\,\text{NADH} + 2\,\text{ATP} + 2\,\text{H}_2\text{O}$$

The 10 Steps

Energy Investment Phase (Steps 1-5)

Hexokinase: Glucose → Glucose-6-phosphate (costs 1 ATP)
Phosphoglucose isomerase: G6P → Fructose-6-phosphate
PFK-1: F6P → Fructose-1,6-bisphosphate (costs 1 ATP, rate-limiting step)
Aldolase: F-1,6-BP → DHAP + G3P
Triose phosphate isomerase: DHAP ↔ G3P

Energy Payoff Phase (Steps 6-10)

G3P dehydrogenase: G3P → 1,3-BPG (produces NADH)
Phosphoglycerate kinase: 1,3-BPG → 3-PG (produces ATP)
Phosphoglycerate mutase: 3-PG → 2-PG
Enolase: 2-PG → PEP
Pyruvate kinase: PEP → Pyruvate (produces ATP)

Regulation

Three irreversible, highly regulated steps control glycolytic flux:

•Hexokinase: Inhibited by glucose-6-phosphate (product inhibition)
•PFK-1 (rate-limiting): Activated by AMP, fructose-2,6-bisphosphate; inhibited by ATP, citrate
•Pyruvate kinase: Activated by fructose-1,6-bisphosphate (feedforward); inhibited by ATP, alanine

Thermodynamics of Each Step

A common misconception is that every step of glycolysis is irreversible. In reality, only three steps have large negative in vivo ΔG values and are physiologically irreversible. The remaining seven operate near equilibrium (ΔG ≈ 0) and are freely reversible, which is essential for gluconeogenesis.

Step	Enzyme	ΔG°′ (kJ/mol)	ΔG in vivo (kJ/mol)	Reversible?
1	Hexokinase	−16.7	−33.4	No
2	Phosphoglucose isomerase	+1.7	−2.5	Yes
3	PFK-1	−14.2	−22.2	No
4	Aldolase	+23.8	−1.3	Yes
5	Triose phosphate isomerase	+7.5	+2.5	Yes
6	G3P dehydrogenase	+6.3	−1.7	Yes
7	Phosphoglycerate kinase	−18.8	+1.3	Yes
8	Phosphoglycerate mutase	+4.4	+0.8	Yes
9	Enolase	+7.5	+3.3	Yes
10	Pyruvate kinase	−31.4	−16.7	No

Note: ΔG°′ is the standard free energy change at pH 7.0 and 25°C. The in vivo ΔG values depend on actual intracellular metabolite concentrations and differ substantially from the standard values. A near-zero in vivo ΔG means the reaction is freely reversible under physiological conditions.

Fates of Pyruvate

Pyruvate sits at a critical metabolic branch point. Its fate depends on the organism, tissue type, and oxygen availability.

1. Aerobic Fate: Oxidative Decarboxylation (Pyruvate Dehydrogenase Complex)

Under aerobic conditions, pyruvate is transported into the mitochondrial matrix by the mitochondrial pyruvate carrier (MPC) and oxidatively decarboxylated to acetyl-CoA by the pyruvate dehydrogenase (PDH) complex, a massive multi-enzyme assembly (~9.5 MDa in eukaryotes).

$$\text{Pyruvate} + \text{CoA-SH} + \text{NAD}^+ \longrightarrow \text{Acetyl-CoA} + \text{CO}_2 + \text{NADH}\quad(\Delta G^{\circ\prime} = -33.4\;\text{kJ/mol})$$

The PDH complex consists of three catalytic subunits:

•E1 — Pyruvate decarboxylase: Uses thiamine pyrophosphate (TPP) as cofactor. Decarboxylates pyruvate and transfers the resulting hydroxyethyl group to lipoamide on E2.
•E2 — Dihydrolipoyl transacetylase: Contains a covalently bound lipoic acid on a flexible lipoyl arm. Transfers the acetyl group to CoA-SH, forming acetyl-CoA.
•E3 — Dihydrolipoyl dehydrogenase: Reoxidizes the reduced lipoamide using FAD, then transfers electrons to NAD⁺, generating NADH.

Regulation of PDH:

•Covalent modification: PDH kinase phosphorylates E1 (inactivation); PDH phosphatase dephosphorylates E1 (activation). PDH kinase is stimulated by acetyl-CoA, NADH, and ATP; PDH phosphatase is activated by Ca²⁺ and insulin signaling.
•Allosteric control: Acetyl-CoA and NADH (products) directly inhibit the complex, while CoA-SH and NAD⁺ (substrates) activate it, ensuring activity only when substrates are available and products are consumed.

2. Anaerobic Fate: Homolactic Fermentation

Under anaerobic conditions (e.g., vigorously exercising skeletal muscle, erythrocytes), pyruvate is reduced to L-lactate by lactate dehydrogenase (LDH). This reaction is critical because it regenerates NAD⁺ from NADH, allowing glycolysis to continue in the absence of mitochondrial oxidation.

$$\text{Pyruvate} + \text{NADH} + \text{H}^+ \longrightarrow \text{Lactate} + \text{NAD}^+$$

The lactate produced by muscle is exported to the blood and taken up by the liver, where it is reconverted to glucose via gluconeogenesis. This metabolic cycle between muscle and liver is known as the Cori cycle, named after Carl and Gerty Cori (Nobel Prize 1947). The Cori cycle effectively shifts the metabolic burden of gluconeogenesis (which costs 6 ATP equivalents per glucose) from muscle to liver.

3. Alcoholic Fermentation (Yeast & Some Microorganisms)

In yeast and certain bacteria, pyruvate is first decarboxylated to acetaldehyde by pyruvate decarboxylase (requires TPP), and then acetaldehyde is reduced to ethanol by alcohol dehydrogenase (ADH), regenerating NAD⁺. Note that pyruvate decarboxylase (yeast) is distinct from the E1 subunit of the PDH complex.

$$\text{Pyruvate} \xrightarrow{\text{TPP}} \text{Acetaldehyde} + \text{CO}_2 \xrightarrow{\text{ADH, NADH}} \text{Ethanol} + \text{NAD}^+$$

This process is the basis of brewing, winemaking, and bread-making. The CO₂ released by pyruvate decarboxylase causes bread dough to rise, while the ethanol produced is the key product in alcoholic beverages.

The Warburg Effect

In 1924, the German biochemist Otto Warburg observed that cancer cells preferentially ferment glucose to lactate even in the presence of adequate oxygen — a phenomenon termed aerobic glycolysis or the Warburg effect. This contrasts with normal differentiated cells, which primarily use oxidative phosphorylation when oxygen is available. Warburg hypothesized that mitochondrial dysfunction was the root cause of cancer, though modern understanding reveals the picture is far more nuanced.

Several hypotheses explain why cancer cells favor glycolysis:

•Rapid ATP production: Although glycolysis yields only 2 ATP per glucose (vs. ~30–32 from complete oxidation), the rate of ATP production via glycolysis can be up to 100 times faster, which may be advantageous for rapidly dividing cells.
•Biosynthetic intermediates: Glycolytic intermediates (G6P, F6P, 3-PG, pyruvate) feed into the pentose phosphate pathway, serine biosynthesis, and lipid synthesis, providing the building blocks needed by proliferating cells.
•Adaptation to hypoxia: Solid tumors frequently contain hypoxic regions. Upregulation of glycolytic enzymes (driven by HIF-1α) is maintained even when oxygen returns, conferring a survival advantage.

Clinical relevance: The Warburg effect is the basis of PET (positron emission tomography) scanning using ¹⁸F-FDG (¹⁸F-fluorodeoxyglucose). Cancer cells take up far more glucose than surrounding tissue, so FDG accumulates in tumors, allowing visualization. This is one of the most widely used tools in cancer diagnosis and staging.

Historical Context: The Embden–Meyerhof–Parnas Pathway

The elucidation of glycolysis was one of the great triumphs of early 20th-century biochemistry and represents the first metabolic pathway to be fully described. The pathway is sometimes called the Embden–Meyerhof–Parnas (EMP) pathway in honor of the three scientists who made the most significant contributions:

•Gustav Embden (1874–1933): German biochemist who identified several key intermediates and proposed the overall scheme shortly before his death.
•Otto Meyerhof (1884–1951): Received the Nobel Prize in Physiology or Medicine in 1922 (shared with A.V. Hill) for his work on the relationship between oxygen consumption and lactic acid metabolism in muscle.
•Jakub Parnas (1884–1949): Polish biochemist who established the connection between glycogen breakdown and lactic acid formation, and identified several phosphorylated intermediates.

By the early 1940s, the complete pathway had been established, laying the foundation for modern metabolic biochemistry.

The Pasteur Effect

In 1861, Louis Pasteur observed that yeast consumed far less glucose under aerobic conditions than under anaerobic conditions — an observation now known as the Pasteur effect. The molecular explanation is straightforward: when oxygen is present, oxidative phosphorylation produces abundant ATP, which allosterically inhibits PFK-1 (the rate-limiting enzyme of glycolysis). Additionally, citrate produced by the TCA cycle also inhibits PFK-1. Thus, the cell economizes on glucose consumption because each molecule of glucose yields ~30–32 ATP aerobically versus only 2 ATP via fermentation alone.

The Pasteur effect highlights how ATP acts as a key allosteric signal integrating glycolytic flux with the energetic state of the cell. The Warburg effect in cancer cells represents a notable exception to the Pasteur effect, as these cells maintain high glycolytic rates despite the presence of oxygen.

Substrate-Level Phosphorylation

Glycolysis produces ATP by substrate-level phosphorylation, a mechanism in which a high-energy phosphoryl group is transferred directly from a substrate to ADP, without the involvement of an electrochemical gradient or ATP synthase. This contrasts with oxidative phosphorylation, which couples electron transport to a proton-motive force.

Two glycolytic steps use this mechanism:

•Step 7 — Phosphoglycerate kinase: 1,3-Bisphosphoglycerate (1,3-BPG) has a mixed anhydride bond (acyl phosphate) with a ΔG°′ of hydrolysis of −49.4 kJ/mol. Because this exceeds the ΔG°′ of ATP hydrolysis (−30.5 kJ/mol), the phosphoryl transfer to ADP is thermodynamically favorable (ΔG°′ = −18.8 kJ/mol for the coupled reaction).
•Step 10 — Pyruvate kinase: Phosphoenolpyruvate (PEP) has the highest phosphoryl-group transfer potential of any common biological molecule, with a ΔG°′ of hydrolysis of −61.9 kJ/mol. This enormous driving force comes from tautomerization: the enol form of pyruvate released upon dephosphorylation rapidly and irreversibly tautomerizes to the more stable keto form, making the overall reaction strongly exergonic (ΔG°′ = −31.4 kJ/mol).

Phosphoryl-Group Transfer Potentials (Comparison)

$$\text{PEP}: \Delta G^{\circ\prime}_{\text{hydrolysis}} = -61.9\;\text{kJ/mol} \quad > \quad \text{1,3-BPG}: -49.4\;\text{kJ/mol} \quad > \quad \text{ATP}: -30.5\;\text{kJ/mol}$$

The principle is universal: substrate-level phosphorylation occurs whenever a phosphorylated intermediate has a higher group-transfer potential than ATP. The difference in ΔG°′ values drives the transfer of the phosphoryl group to ADP, producing ATP.

🎥 Video Lecture: Glycolysis

Ninja Nerd — Comprehensive glycolysis lecture covering all 10 steps, enzymes, and regulation

🎥 MIT 7.05: Glycolysis I

MIT OCW — General Biochemistry (Spring 2020) — Introduction to glycolysis pathway and bioenergetics

🎥 MIT 7.05: Glycolysis II & Regulation

MIT OCW — General Biochemistry (Spring 2020) — Allosteric regulation of glycolytic enzymes

Interactive: Glycolysis Flux Model

Python

Kinetic simulation of glycolytic flux with Michaelis-Menten enzyme kinetics and allosteric regulation

glycolysis_flux.py110 lines

import numpy as np
import matplotlib.pyplot as plt

# ============================================================
# Glycolysis Flux Model: Simplified 3-enzyme kinetic simulation
# Models the three irreversible (regulated) steps of glycolysis
# using Michaelis-Menten kinetics with allosteric modulation.
# ============================================================
from scipy.integrate import solve_ivp

# Michaelis-Menten with Hill-type allosteric modulation
def mm_rate(Vmax, S, Km, activator=0, Ka=1.0, inhibitor=0, Ki=1.0):
    """Michaelis-Menten rate with optional activation/inhibition."""
    act_factor = 1.0 + activator / Ka
    inh_factor = 1.0 + inhibitor / Ki
    return Vmax * act_factor / inh_factor * S / (Km + S)

def glycolysis_ode(t, y):
    Glc, G6P, F16BP, PYR, ATP, ADP, NADp, NADH = y

# Step 1: Hexokinase (Glc -> G6P, consumes ATP)
    v1 = mm_rate(1.0, Glc, 0.1, inhibitor=G6P, Ki=0.5) * ATP / (ATP + 0.5)

# Step 3: PFK-1 (G6P -> F16BP, rate-limiting, consumes ATP)
    v3 = mm_rate(1.5, G6P, 0.3, activator=ADP, Ka=0.2, inhibitor=ATP, Ki=2.0) * ATP / (ATP + 0.5)

# Steps 6-10 lumped: (F16BP -> 2 PYR, produces 4 ATP and 2 NADH)
    v_payoff = mm_rate(3.0, F16BP, 0.2, activator=F16BP, Ka=0.5) * NADp / (NADp + 0.1)

# Glucose input (constant supply)
    v_input = 0.5

# Pyruvate consumption (downstream metabolism)
    v_pyr_out = 0.3 * PYR / (0.1 + PYR)

# NAD+ regeneration (represents LDH or mitochondrial shuttle)
    v_nad_regen = 0.8 * NADH / (0.2 + NADH)

# ATP consumption (cellular demand)
    v_atp_use = 0.6 * ATP / (0.5 + ATP)

dGlc   = v_input - v1
    dG6P   = v1 - v3
    dF16BP = v3 - v_payoff
    dPYR   = 2.0 * v_payoff - v_pyr_out
    dATP   = -v1 - v3 + 4.0 * v_payoff - v_atp_use
    dADP   = v1 + v3 - 4.0 * v_payoff + v_atp_use
    dNADp  = -2.0 * v_payoff + v_nad_regen
    dNADH  = 2.0 * v_payoff - v_nad_regen

return [dGlc, dG6P, dF16BP, dPYR, dATP, dADP, dNADp, dNADH]

# Initial conditions: [Glc, G6P, F16BP, PYR, ATP, ADP, NAD+, NADH]
y0 = [5.0, 0.1, 0.05, 0.1, 3.0, 1.0, 2.0, 0.1]
t_span = (0, 30)
t_eval = np.linspace(0, 30, 500)

sol = solve_ivp(glycolysis_ode, t_span, y0, t_eval=t_eval, method='RK45', max_step=0.1)

fig, axes = plt.subplots(2, 2, figsize=(12, 9))
fig.suptitle('Glycolysis Flux Model: Kinetic Simulation', fontsize=14, fontweight='bold')

# Metabolites
axes[0,0].plot(sol.t, sol.y[0], 'b-', label='Glucose', linewidth=2)
axes[0,0].plot(sol.t, sol.y[1], 'g-', label='G6P', linewidth=2)
axes[0,0].plot(sol.t, sol.y[2], 'r-', label='F-1,6-BP', linewidth=2)
axes[0,0].plot(sol.t, sol.y[3], 'm-', label='Pyruvate', linewidth=2)
axes[0,0].set_xlabel('Time (a.u.)')
axes[0,0].set_ylabel('Concentration (mM)')
axes[0,0].set_title('Glycolytic Intermediates')
axes[0,0].legend(fontsize=8)
axes[0,0].grid(True, alpha=0.3)

# ATP/ADP
axes[0,1].plot(sol.t, sol.y[4], 'r-', label='ATP', linewidth=2)
axes[0,1].plot(sol.t, sol.y[5], 'b-', label='ADP', linewidth=2)
axes[0,1].plot(sol.t, sol.y[4]/(sol.y[4]+sol.y[5]+1e-10), 'k--', label='ATP/(ATP+ADP)', linewidth=1.5)
axes[0,1].set_xlabel('Time (a.u.)')
axes[0,1].set_ylabel('Concentration (mM)')
axes[0,1].set_title('Energy Charge')
axes[0,1].legend(fontsize=8)
axes[0,1].grid(True, alpha=0.3)

# NAD+/NADH
axes[1,0].plot(sol.t, sol.y[6], 'orange', label='NAD+', linewidth=2)
axes[1,0].plot(sol.t, sol.y[7], 'purple', label='NADH', linewidth=2)
axes[1,0].set_xlabel('Time (a.u.)')
axes[1,0].set_ylabel('Concentration (mM)')
axes[1,0].set_title('Redox State')
axes[1,0].legend(fontsize=8)
axes[1,0].grid(True, alpha=0.3)

# Flux rates
v1 = [mm_rate(1.0, sol.y[0,i], 0.1, inhibitor=sol.y[1,i], Ki=0.5) * sol.y[4,i]/(sol.y[4,i]+0.5) for i in range(len(sol.t))]
v3 = [mm_rate(1.5, sol.y[1,i], 0.3, activator=sol.y[5,i], Ka=0.2, inhibitor=sol.y[4,i], Ki=2.0) * sol.y[4,i]/(sol.y[4,i]+0.5) for i in range(len(sol.t))]
axes[1,1].plot(sol.t, v1, 'b-', label='Hexokinase (v1)', linewidth=2)
axes[1,1].plot(sol.t, v3, 'r-', label='PFK-1 (v3)', linewidth=2)
axes[1,1].set_xlabel('Time (a.u.)')
axes[1,1].set_ylabel('Rate (mM/s)')
axes[1,1].set_title('Enzyme Fluxes')
axes[1,1].legend(fontsize=8)
axes[1,1].grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig('glycolysis_flux.png', dpi=150, bbox_inches='tight')
print("Glycolysis flux model simulation complete.")
print(f"Steady-state ATP: {sol.y[4,-1]:.2f} mM")
print(f"Steady-state Pyruvate: {sol.y[3,-1]:.2f} mM")
print(f"ATP/ADP ratio: {sol.y[4,-1]/max(sol.y[5,-1],1e-10):.2f}")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Understanding the Glycolysis Simulation

This simulation models glycolysis as a system of coupled ordinary differential equations (ODEs), capturing the three irreversible regulatory steps — hexokinase, phosphofructokinase-1 (PFK-1), and the lumped payoff phase — using Michaelis-Menten kinetics with allosteric modulation. The model tracks eight state variables simultaneously: glucose, glucose-6-phosphate (G6P), fructose-1,6-bisphosphate (F-1,6-BP), pyruvate, ATP, ADP, NAD⁺, and NADH.

Panel 1 — Glycolytic Intermediates

Shows how glucose is consumed and intermediates (G6P, F-1,6-BP) transiently accumulate before reaching steady state. Pyruvate rises and stabilizes as production matches downstream consumption. The initial transient reflects the “investment phase” where 2 ATP are consumed before the pathway begins producing ATP in the payoff phase.

Panel 2 — Energy Charge (ATP/ADP)

Tracks the ATP/ADP ratio over time. The initial dip in ATP reflects the investment phase costs (hexokinase and PFK-1 each consume 1 ATP). ATP then recovers as the payoff phase produces 4 ATP per glucose. The energy charge (ATP/(ATP+ADP)) stabilizes near 0.7–0.9, consistent with values observed in living cells. This ratio is the master regulator of glycolytic flux.

Panel 3 — Redox State (NAD⁺/NADH)

Shows the balance between oxidized (NAD⁺) and reduced (NADH) nicotinamide cofactors. NADH is produced at step 6 (glyceraldehyde-3-phosphate dehydrogenase) and must be continuously regenerated to NAD⁺ for glycolysis to continue. Under aerobic conditions, the mitochondrial shuttles handle this; under anaerobic conditions, lactate dehydrogenase (LDH) performs the regeneration. If NADH accumulates and NAD⁺ is depleted, glycolysis halts.

Panel 4 — Enzyme Fluxes

Compares the rates of hexokinase (v1) and PFK-1 (v3) over time. PFK-1 is the committed and rate-limiting step — its flux determines overall glycolytic throughput. Notice how PFK-1 flux is modulated by the ATP/ADP ratio: when ATP is high, PFK-1 is inhibited; when ADP accumulates, PFK-1 is activated. This allosteric regulation is the molecular basis of the Pasteur effect.

Biological significance: This model demonstrates why glycolysis is self-regulating: the ATP produced by the pathway inhibits PFK-1, creating a negative feedback loop that prevents wasteful glucose consumption when energy is abundant. You can modify the code to explore what happens when you remove allosteric regulation (set Ki and Ka to very large values) — the pathway loses its self-limiting behavior and oscillates wildly.

Citric Acid Cycle (TCA / Krebs Cycle)

8 Reactions

Overview

The citric acid cycle (also called the tricarboxylic acid cycle or Krebs cycle) is the central metabolic hub of aerobic organisms. It takes place in the mitochondrial matrix and oxidizes acetyl-CoA (derived from pyruvate, fatty acids, or amino acids) to CO₂, generating high-energy electron carriers (NADH and FADH₂) that feed the electron transport chain.

Overall Reaction (per turn)

$$\text{Acetyl-CoA} + 3\,\text{NAD}^+ + \text{FAD} + \text{GDP} + P_i + 2\,\text{H}_2\text{O} \longrightarrow \text{CoA-SH} + 3\,\text{NADH} + \text{FADH}_2 + \text{GTP} + 2\,\text{CO}_2$$

The 8 Reactions

Step	Enzyme	Reaction	Products
1	Citrate synthase	Acetyl-CoA + OAA → Citrate	CoA-SH
2	Aconitase	Citrate → Isocitrate	—
3	Isocitrate DH	Isocitrate → α-Ketoglutarate	NADH + CO₂
4	α-KG DH complex	α-KG → Succinyl-CoA	NADH + CO₂
5	Succinyl-CoA synthetase	Succinyl-CoA → Succinate	GTP
6	Succinate DH (Complex II)	Succinate → Fumarate	FADH₂
7	Fumarase	Fumarate → Malate	—
8	Malate DH	Malate → Oxaloacetate	NADH

Regulation

•Citrate synthase: Inhibited by NADH, succinyl-CoA, citrate, ATP
•Isocitrate DH: Activated by ADP, Ca²⁺; inhibited by NADH, ATP
•α-Ketoglutarate DH: Activated by Ca²⁺; inhibited by NADH, succinyl-CoA

Anaplerotic Reactions

TCA intermediates are constantly siphoned off for biosynthesis. Anaplerotic (“filling up”) reactions replenish them. The most important is pyruvate carboxylase: pyruvate + CO₂ + ATP → oxaloacetate, activated by acetyl-CoA.

Pyruvate Dehydrogenase Complex: The Gateway to the TCA Cycle

Before acetyl-CoA can enter the TCA cycle, pyruvate must be oxidatively decarboxylated by the pyruvate dehydrogenase (PDH) complex. This is one of the largest known multi-enzyme complexes in the cell, with a molecular mass of approximately 9.5 MDa in eukaryotes (comparable to a small ribosome). The complex is located in the mitochondrial matrix and represents a critical, irreversible commitment step: once pyruvate is converted to acetyl-CoA, its carbons are destined for complete oxidation or lipid synthesis, not gluconeogenesis.

The complex contains three catalytic enzymes, each present in multiple copies arranged around a central core:

•E1 — Pyruvate decarboxylase: Requires thiamine pyrophosphate (TPP, vitamin B1) as a cofactor. Catalyzes the decarboxylation of pyruvate and transfers the resulting two-carbon hydroxyethyl group to the lipoamide prosthetic group of E2. E1 is the target of regulatory phosphorylation.
•E2 — Dihydrolipoyl transacetylase: Contains a covalently attached lipoamide (lipoic acid) swinging arm and uses coenzyme A (CoA, from vitamin B5/pantothenate). E2 forms the structural core of the complex (60-mer in eukaryotes) and transfers the acetyl group from lipoamide to CoA-SH, generating acetyl-CoA.
•E3 — Dihydrolipoyl dehydrogenase: Contains FAD (from vitamin B2/riboflavin) as a tightly bound prosthetic group and uses NAD⁺ (from vitamin B3/niacin) as electron acceptor. Reoxidizes the reduced dihydrolipoamide on E2, passing electrons through FAD to NAD⁺, generating NADH.

Five Cofactors — Five Vitamins

The PDH complex requires five coenzymes, each derived from a vitamin: TPP (B1), CoA (B5), lipoic acid, FAD (B2), and NAD⁺ (B3). A deficiency in any of these vitamins impairs the complex. Beriberi (thiamine deficiency) is the classic example, presenting with peripheral neuropathy and cardiac failure due to impaired PDH and α-ketoglutarate dehydrogenase activity.

Regulation by phosphorylation: The PDH complex is regulated by a dedicated PDH kinase (phosphorylates E1, inactivating the complex) and PDH phosphatase (dephosphorylates E1, reactivating the complex). PDH kinase is activated by high ratios of acetyl-CoA/CoA, NADH/NAD⁺, and ATP/ADP — signals that the cell has sufficient energy and biosynthetic substrates. PDH phosphatase is activated by Ca²⁺ (important in exercising muscle) and by insulin (in adipose tissue and liver, promoting fat synthesis).

Vitamin Cofactors in the TCA Cycle

The TCA cycle and the PDH complex together depend on coenzymes derived from four B vitamins. Understanding these connections is clinically important because vitamin deficiencies directly impair energy metabolism.

Vitamin	Coenzyme Form	TCA/PDH Enzymes Using It
B1 (Thiamine)	Thiamine pyrophosphate (TPP)	α-Ketoglutarate DH (E1), PDH (E1)
B2 (Riboflavin)	FAD	Succinate DH, α-Ketoglutarate DH (E3), PDH (E3)
B3 (Niacin)	NAD⁺	Isocitrate DH, α-Ketoglutarate DH (E3), Malate DH, PDH (E3)
B5 (Pantothenate)	Coenzyme A (CoA)	Citrate synthase, α-Ketoglutarate DH (E2), Succinyl-CoA synthetase, PDH (E2)

Lipoic acid, a cofactor for both α-ketoglutarate dehydrogenase and PDH (on E2), is not formally classified as a vitamin but is synthesized endogenously in mitochondria using an iron–sulfur cluster-dependent mechanism.

Cataplerotic Reactions & Biosynthetic Roles

While anaplerotic reactions replenish TCA intermediates, cataplerotic reactions remove them for biosynthesis. The TCA cycle is therefore an amphibolic pathway — serving both catabolic and anabolic functions. The key cataplerotic exits are:

•Citrate → Fatty acid synthesis: Citrate is exported from the mitochondrial matrix to the cytosol via the citrate shuttle (tricarboxylate transport system). In the cytosol, ATP-citrate lyase cleaves citrate into oxaloacetate and acetyl-CoA. The cytosolic acetyl-CoA serves as the building block for de novo fatty acid and cholesterol synthesis. This pathway is especially active in liver and adipose tissue when energy supply exceeds demand.
•α-Ketoglutarate → Amino acid & purine synthesis: α-Ketoglutarate is transaminated to glutamate by aminotransferases, or reductively aminated by glutamate dehydrogenase. Glutamate is the precursor for glutamine, proline, and arginine. Its nitrogen atoms also contribute to purine ring synthesis.
•Succinyl-CoA → Heme synthesis: Succinyl-CoA condenses with glycine (catalyzed by δ-aminolevulinate synthase, ALA synthase) to form δ-aminolevulinic acid (δ-ALA), the committed first step in porphyrin and heme biosynthesis. This reaction occurs in the mitochondrial matrix and is the rate-limiting step of heme synthesis.
•Oxaloacetate → Gluconeogenesis & amino acid/pyrimidine synthesis: OAA is transaminated to aspartate, which donates nitrogen atoms to purine synthesis and provides the entire carbon skeleton for pyrimidine synthesis. OAA is also the starting point for gluconeogenesis via phosphoenolpyruvate carboxykinase (PEPCK), which converts OAA to PEP and CO₂.
•Fumarate ←→ Urea cycle connection: The urea cycle produces fumarate as a byproduct (from argininosuccinate lyase). This fumarate re-enters the TCA cycle, linking nitrogen disposal to energy metabolism. This connection is sometimes called the Krebs bicycle because the two cycles share the fumarate/aspartate segment.

Energetics of the TCA Cycle

The overall ΔG°′ for one turn of the TCA cycle is strongly negative, but the individual reactions span a wide range. Three reactions are highly exergonic and drive the cycle forward irreversibly; remarkably, the malate dehydrogenase reaction is actually endergonic under standard conditions but is pulled forward by the highly favorable citrate synthase reaction that immediately consumes its product (oxaloacetate).

Step	Enzyme	ΔG°′ (kJ/mol)	Products
1	Citrate synthase	−31.4	Citrate, CoA-SH
2	Aconitase	+6.3	Isocitrate
3	Isocitrate dehydrogenase	−8.4	α-Ketoglutarate, CO₂, NADH
4	α-Ketoglutarate dehydrogenase	−33.5	Succinyl-CoA, CO₂, NADH
5	Succinyl-CoA synthetase	−2.9	Succinate, GTP (or ATP), CoA-SH
6	Succinate dehydrogenase	0	Fumarate, FADH₂
7	Fumarase	−3.8	L-Malate
8	Malate dehydrogenase	+29.7	Oxaloacetate, NADH

The malate dehydrogenase reaction (ΔG°′ = +29.7 kJ/mol) is the most endergonic step and would not proceed forward if oxaloacetate were allowed to accumulate. However, the in vivo concentration of OAA is maintained at extremely low levels (micromolar range) because citrate synthase (ΔG°′ = −31.4 kJ/mol) rapidly condenses it with acetyl-CoA. This coupling of a thermodynamically unfavorable reaction with a highly favorable one is a recurring principle in metabolic pathway design.

Net Yield per Turn of the TCA Cycle

$$3\,\text{NADH} + 1\,\text{FADH}_2 + 1\,\text{GTP} + 2\,\text{CO}_2 \quad\Rightarrow\quad \approx 10\;\text{ATP equivalents per acetyl-CoA}$$

Clinical Connections

Disruptions of TCA cycle enzymes or their cofactors have profound clinical consequences, given the cycle’s central role in aerobic energy metabolism.

•Fluoroacetate poisoning (lethal synthesis): Fluoroacetate, found in certain South African and Australian plants (and used as rodenticide “Compound 1080”), is converted to fluorocitrate by citrate synthase. Fluorocitrate is a potent competitive inhibitor of aconitase, blocking the TCA cycle at step 2. Citrate accumulates massively, chelating Ca²⁺ and leading to cardiac arrest. This is termed “lethal synthesis” because the cell’s own enzyme creates the toxic product.
•Leigh syndrome: A devastating infantile neurodegenerative disease caused by mutations in genes encoding mitochondrial complex subunits, PDH complex components, or mitochondrial tRNA genes. Impaired oxidative phosphorylation and TCA cycle function lead to bilateral symmetric lesions in the basal ganglia and brainstem. Lactic acidosis is a hallmark finding. Prognosis is generally poor, with most patients succumbing in early childhood.
•Arsenic poisoning: Arsenite (As³⁺) binds to the sulfhydryl groups of lipoic acid on E2 subunits of both the PDH complex and α-ketoglutarate dehydrogenase, forming a stable chelate that inactivates these enzymes. This blocks both the entry point and a critical step within the TCA cycle. Chronic arsenic exposure (contaminated drinking water) is associated with peripheral neuropathy, keratoses, and increased cancer risk.
•Fumarase deficiency: An extremely rare autosomal recessive disorder (fewer than 100 cases reported worldwide) caused by mutations in the FH gene encoding fumarase. Patients present with severe neurological impairment, including microcephaly, seizures, and profound developmental delay. Fumarate accumulates in body fluids. Interestingly, heterozygous loss of FH predisposes to hereditary leiomyomatosis and renal cell carcinoma (HLRCC), linking fumarase to tumor suppression — accumulated fumarate inhibits prolyl hydroxylases, stabilizing HIF-1α and promoting a pseudohypoxic state that drives tumorigenesis.

🎥 Video Lecture: The Krebs Cycle

Ninja Nerd — Complete TCA cycle with all 8 reactions, enzyme regulation, and clinical correlations

🎥 MIT 5.07SC: TCA Cycle

MIT OCW — Biological Chemistry I (Fall 2013) — Detailed TCA enzymes, regulation, and anaplerotic reactions

🎥 MIT 7.05: TCA Cycle II

MIT OCW — General Biochemistry (Spring 2020) — TCA cycle regulation and integration with other pathways

🎥 Harvard Online: Overview of the Citric Acid Cycle

Harvard Online — Cell Biology — Overview of TCA cycle reactions, enzyme regulation, and the central role of acetyl-CoA oxidation in aerobic metabolism

Interactive: Krebs (TCA) Cycle Simulation

Python

Kinetic model of all 8 TCA cycle reactions with intermediate dynamics and energy yield analysis

krebs_cycle.py175 lines

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# ============================================================
# Krebs (TCA) Cycle Simulation
# Simplified kinetic model of the 8 reactions, tracking all
# intermediates, NADH/FADH2 production, and energy yield.
# ============================================================

def mm(S, Vmax, Km):
    return Vmax * S / (Km + S)

def tca_cycle(t, y):
    # Intermediates
    AcCoA, Cit, IsoCit, aKG, SucCoA, Suc, Fum, Mal, OAA = y[:9]
    # Cofactors
    NADH, FADH2, GTP, CoASH = y[9:13]

# Conservation: total NAD pool = 5 mM, FAD pool = 1 mM
    NADp = max(5.0 - NADH, 0.01)
    FAD  = max(1.0 - FADH2, 0.01)
    GDP  = max(1.0 - GTP, 0.01)

# Clamp concentrations
    vals = [AcCoA, Cit, IsoCit, aKG, SucCoA, Suc, Fum, Mal, OAA,
            NADH, FADH2, GTP, CoASH]
    AcCoA, Cit, IsoCit, aKG, SucCoA, Suc, Fum, Mal, OAA, \
        NADH, FADH2, GTP, CoASH = [max(v, 0) for v in vals]

# === 8 REACTIONS ===
    # 1. Citrate synthase: AcCoA + OAA -> Citrate + CoA
    v1 = 2.0 * AcCoA / (0.05 + AcCoA) * OAA / (0.01 + OAA)

# 2. Aconitase: Citrate <-> Isocitrate (near equilibrium)
    v2 = 5.0 * (Cit - IsoCit / 0.067) / (0.5 + Cit + IsoCit)

# 3. Isocitrate DH: IsoCit -> alpha-KG + CO2 + NADH
    v3 = 1.5 * mm(IsoCit, 1.0, 0.1) * NADp / (NADp + 0.1)

# 4. alpha-KG DH: aKG -> Succinyl-CoA + CO2 + NADH
    v4 = 1.2 * mm(aKG, 1.0, 0.15) * NADp / (NADp + 0.1) * CoASH / (CoASH + 0.05)

# 5. Succinyl-CoA synthetase: SucCoA -> Succinate + GTP + CoA
    v5 = 2.5 * mm(SucCoA, 1.0, 0.1) * GDP / (GDP + 0.05)

# 6. Succinate DH (Complex II): Succinate -> Fumarate + FADH2
    v6 = 1.8 * mm(Suc, 1.0, 0.3) * FAD / (FAD + 0.05)

# 7. Fumarase: Fumarate -> Malate
    v7 = 4.0 * (Fum - Mal / 4.0) / (0.3 + Fum + Mal)

# 8. Malate DH: Malate -> OAA + NADH
    v8 = 1.0 * mm(Mal, 0.8, 0.15) * NADp / (NADp + 0.2)

# Acetyl-CoA supply (from pyruvate DH, constant)
    v_supply = 0.5

# Cofactor regeneration (ETC consumes NADH, FADH2; ATPase uses GTP)
    v_nadh_drain = 0.8 * NADH / (0.2 + NADH)
    v_fadh2_drain = 0.5 * FADH2 / (0.1 + FADH2)
    v_gtp_drain = 0.3 * GTP / (0.1 + GTP)

# Derivatives
    dAcCoA  = v_supply - v1
    dCit    = v1 - v2
    dIsoCit = v2 - v3
    daKG    = v3 - v4
    dSucCoA = v4 - v5
    dSuc    = v5 - v6
    dFum    = v6 - v7
    dMal    = v7 - v8
    dOAA    = v8 - v1  # regenerated!

dNADH   = v3 + v4 + v8 - v_nadh_drain
    dFADH2  = v6 - v_fadh2_drain
    dGTP    = v5 - v_gtp_drain
    dCoASH  = v1 + v5 - v4 - v_supply  # CoA balance

return [dAcCoA, dCit, dIsoCit, daKG, dSucCoA, dSuc, dFum, dMal, dOAA,
            dNADH, dFADH2, dGTP, dCoASH]

# Initial conditions
y0 = [0.2,   # AcCoA
      0.4,   # Citrate
      0.05,  # Isocitrate
      0.1,   # alpha-KG
      0.05,  # Succinyl-CoA
      0.3,   # Succinate
      0.1,   # Fumarate
      0.2,   # Malate
      0.05,  # OAA
      0.5,   # NADH
      0.1,   # FADH2
      0.1,   # GTP
      0.8]   # CoA-SH

t_span = (0, 40)
t_eval = np.linspace(0, 40, 800)

sol = solve_ivp(tca_cycle, t_span, y0, t_eval=t_eval, method='RK45', max_step=0.05)

# Labels
intermediates = ['Acetyl-CoA', 'Citrate', 'Isocitrate', 'α-KG',
                 'Succinyl-CoA', 'Succinate', 'Fumarate', 'Malate', 'OAA']
colors = ['#ef4444', '#f97316', '#eab308', '#22c55e',
          '#14b8a6', '#06b6d4', '#8b5cf6', '#ec4899', '#f43f5e']

fig, axes = plt.subplots(2, 2, figsize=(14, 10))
fig.suptitle('Krebs (TCA) Cycle: Kinetic Simulation', fontsize=14, fontweight='bold')

# Panel 1: All 9 intermediates
for i in range(9):
    axes[0,0].plot(sol.t, sol.y[i], color=colors[i], linewidth=1.8, label=intermediates[i])
axes[0,0].set_xlabel('Time (a.u.)')
axes[0,0].set_ylabel('Concentration (mM)')
axes[0,0].set_title('TCA Cycle Intermediates')
axes[0,0].legend(fontsize=6, ncol=3, loc='upper right')
axes[0,0].grid(True, alpha=0.2)

# Panel 2: Energy carriers
axes[0,1].plot(sol.t, sol.y[9], 'r-', label='NADH', linewidth=2)
axes[0,1].plot(sol.t, 5.0 - sol.y[9], 'r--', label='NAD+', linewidth=1.5, alpha=0.6)
axes[0,1].plot(sol.t, sol.y[10], 'b-', label='FADH₂', linewidth=2)
axes[0,1].plot(sol.t, sol.y[11], 'g-', label='GTP', linewidth=2)
axes[0,1].set_xlabel('Time (a.u.)')
axes[0,1].set_ylabel('Concentration (mM)')
axes[0,1].set_title('Energy Carrier Production')
axes[0,1].legend(fontsize=8)
axes[0,1].grid(True, alpha=0.2)

# Panel 3: Steady-state concentrations (bar chart)
ss_vals = sol.y[:9, -1]
bars = axes[1,0].barh(intermediates, ss_vals, color=colors, edgecolor='white', linewidth=0.5)
for bar, val in zip(bars, ss_vals):
    axes[1,0].text(val + 0.01, bar.get_y() + bar.get_height()/2,
                  f'{val:.3f}', va='center', fontsize=8)
axes[1,0].set_xlabel('Steady-State Concentration (mM)')
axes[1,0].set_title('Steady-State Metabolite Pool Sizes')
axes[1,0].grid(axis='x', alpha=0.2)

# Panel 4: Cumulative energy yield per turn
# At steady state, compute rates for each step
# For visualization, show the textbook per-turn yield
steps = ['Step 3\n(IDH)', 'Step 4\n(α-KG DH)', 'Step 5\n(SCS)',
         'Step 6\n(SDH)', 'Step 8\n(MDH)']
carriers = ['NADH', 'NADH', 'GTP', 'FADH₂', 'NADH']
atp_equiv = [2.5, 2.5, 1.0, 1.5, 2.5]  # ATP equivalents
bar_colors = ['#ef4444', '#ef4444', '#22c55e', '#3b82f6', '#ef4444']

bars4 = axes[1,1].bar(steps, atp_equiv, color=bar_colors, edgecolor='white', linewidth=0.5)
for bar, val, carrier in zip(bars4, atp_equiv, carriers):
    axes[1,1].text(bar.get_x() + bar.get_width()/2, bar.get_height() + 0.05,
                  f'{val} ATP\n({carrier})', ha='center', fontsize=7, fontweight='bold')
axes[1,1].set_ylabel('ATP Equivalents')
axes[1,1].set_title(f'Energy Yield per TCA Turn (Total: {sum(atp_equiv):.0f} ATP equiv.)')
axes[1,1].set_ylim(0, 3.5)
axes[1,1].grid(axis='y', alpha=0.2)

plt.tight_layout()
plt.savefig('krebs_cycle.png', dpi=150, bbox_inches='tight')

print("Krebs Cycle Simulation Complete")
print("=" * 50)
print(f"\nSteady-state concentrations:")
for name, val in zip(intermediates, ss_vals):
    print(f"  {name:15s}: {val:.4f} mM")
print(f"\nSteady-state NADH:  {sol.y[9,-1]:.3f} mM (of 5.0 total NAD pool)")
print(f"Steady-state FADH2: {sol.y[10,-1]:.3f} mM")
print(f"Steady-state GTP:   {sol.y[11,-1]:.3f} mM")
print(f"\nPer turn ATP equivalents: {sum(atp_equiv):.0f}")
print(f"  3 NADH x 2.5 = 7.5 ATP")
print(f"  1 FADH2 x 1.5 = 1.5 ATP")
print(f"  1 GTP         = 1.0 ATP")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Understanding the TCA Cycle Simulation

This simulation models all eight enzymatic reactions of the citric acid cycle as a coupled ODE system, tracking twelve metabolites simultaneously: acetyl-CoA, oxaloacetate, citrate, isocitrate, α-ketoglutarate, succinyl-CoA, succinate, fumarate, malate, NADH, FADH₂, and GTP. Each reaction uses Michaelis-Menten kinetics with parameters reflecting the known regulatory properties of the enzymes.

Panel 1 — TCA Intermediates (Upper Left)

Tracks the concentrations of all eight organic acid intermediates over time. Notice how oxaloacetate (OAA) remains at very low concentrations — this is biologically accurate because citrate synthase rapidly consumes OAA, pulling the malate dehydrogenase reaction forward despite its unfavorable thermodynamics (ΔG°′ = +29.7 kJ/mol). The relative pool sizes at steady state reflect each intermediate’s turnover rate.

Panel 2 — NADH and FADH₂ Accumulation (Upper Right)

Shows the production of reduced electron carriers. Three NADH molecules are generated per turn (at isocitrate DH, α-ketoglutarate DH, and malate DH), while one FADH₂ is produced at succinate dehydrogenase. These carriers deliver electrons to the ETC: NADH enters at Complex I (yielding ~2.5 ATP), FADH₂ at Complex II (yielding ~1.5 ATP). The NADH/FADH₂ ratio of 3:1 is one reason the TCA cycle is so energetically productive.

Panel 3 — Enzyme Fluxes (Lower Left)

Compares the reaction rates of all eight enzymes. The three rate-limiting steps — citrate synthase, isocitrate DH, and α-ketoglutarate DH — set the overall pace. At steady state, all fluxes must equalize (conservation of mass). Perturbations in any one enzyme propagate through the entire cycle, demonstrating why the TCA cycle operates as an integrated unit rather than eight independent reactions.

Panel 4 — Energy Yield per Turn (Lower Right)

Bar chart showing the ATP-equivalent yield from each energy-producing step. The total of ~10 ATP equivalents per acetyl-CoA (7.5 from 3 NADH + 1.5 from FADH₂ + 1.0 from GTP) makes the TCA cycle the most productive segment of aerobic glucose catabolism. Since each glucose generates 2 acetyl-CoA, the TCA cycle contributes ~20 of the ~30–32 total ATP.

Key insight: The TCA cycle doesn’t directly produce much ATP (only 1 GTP via substrate-level phosphorylation). Its primary role is to harvest high-energy electrons as NADH and FADH₂, which carry their energy to the ETC where the vast majority of ATP is synthesized via chemiosmotic coupling. Try modifying the initial acetyl-CoA concentration to see how substrate availability affects cycle flux and energy output.

Oxidative Phosphorylation

ETC + ATP Synthase

The Electron Transport Chain

The electron transport chain (ETC) is a series of four protein complexes embedded in the inner mitochondrial membrane. Electrons from NADH and FADH₂ flow through these complexes in order of increasing reduction potential, releasing energy that is used to pump protons (H⁺) from the matrix into the intermembrane space, creating the proton motive force (PMF).

Complex I

NADH:ubiquinone oxidoreductase

NADH → CoQ

4 H⁺ pumped

Complex II

Succinate dehydrogenase

FADH₂ → CoQ

0 H⁺ pumped

Complex III

Cytochrome bc₁

CoQH₂ → Cyt c

4 H⁺ pumped

Complex IV

Cytochrome c oxidase

Cyt c → O₂ → H₂O

2 H⁺ pumped

Chemiosmotic Coupling

Peter Mitchell’s chemiosmotic hypothesis (Nobel Prize, 1978) states that the proton gradient across the inner mitochondrial membrane stores energy as the proton motive force (PMF), which drives ATP synthesis through ATP synthase (Complex V).

Proton Motive Force

$$\Delta p = \Delta\psi - \frac{2.303\,RT}{F}\,\Delta\text{pH}$$

where Δψ is the membrane potential (~140 mV) and ΔpH ≈ 0.5–1.0 units

ATP Synthase (Complex V)

ATP synthase is a rotary molecular motor. The F₀ subunit spans the membrane and acts as a proton channel; the F₁ subunit in the matrix catalyzes ATP synthesis via Paul Boyer’s binding change mechanism (Nobel Prize, 1997). Approximately 4 H⁺ are needed per ATP, giving P/O ratios of ~2.5 ATP/NADH and ~1.5 ATP/FADH₂.

Free Energy from Electron Transfer

$$\Delta G = -nF\Delta E^{\circ\prime}$$

For NADH → O₂: ΔE°′ = +1.14 V, so ΔG = −220 kJ/mol (enough for ~7 ATP theoretically)

NADH Shuttle Systems

Cytoplasmic NADH generated by glycolysis cannot directly cross the inner mitochondrial membrane, which is impermeable to NADH and NAD⁺. Instead, two shuttle systems transfer the reducing equivalents into the mitochondrial matrix, each with different tissue distributions and ATP yields.

Malate-Aspartate Shuttle

Predominant in liver, heart, and kidney

●Cytoplasmic NADH reduces oxaloacetate to malate via cytoplasmic malate dehydrogenase
●Malate enters the matrix through the malate-α-ketoglutarate antiporter
●Matrix malate dehydrogenase regenerates NADH from malate, producing oxaloacetate
●Oxaloacetate is converted to aspartate (by transamination) and exported via the glutamate-aspartate antiporter
●Yields 2.5 ATP per NADH (electrons enter at Complex I). Reversible

Glycerol-3-Phosphate Shuttle

Predominant in brain and skeletal muscle

●Cytoplasmic NADH reduces dihydroxyacetone phosphate (DHAP) to glycerol-3-phosphate via cytoplasmic glycerol-3-phosphate dehydrogenase (NAD⁺-linked)
●Mitochondrial glycerol-3-phosphate dehydrogenase (FAD-linked), located on the outer face of the inner membrane, reoxidizes glycerol-3-P back to DHAP
●Electrons transfer from FADH₂ directly to coenzyme Q (ubiquinone), bypassing Complex I
●Yields only 1.5 ATP per NADH (electrons enter at CoQ level, equivalent to FADH₂). Irreversible

Reactive Oxygen Species (ROS)

Approximately 1–2% of electrons passing through the ETC “leak” prematurely to molecular oxygen, primarily at Complex I (FMN site and iron-sulfur clusters) and Complex III (semiquinone intermediate at the Q_o site). This partial reduction of O₂ generates the superoxide anion radical (O₂•⁻), the primary ROS in mitochondria.

●Superoxide dismutase (SOD) converts O₂•⁻ to hydrogen peroxide (H₂O₂). MnSOD (SOD2) in the matrix; CuZnSOD (SOD1) in the cytoplasm
●Catalase (in peroxisomes) and glutathione peroxidase (requires selenium) detoxify H₂O₂ to H₂O
●Excessive ROS causes oxidative stress, damaging DNA (8-oxoguanine lesions), membrane lipids (lipid peroxidation), and proteins (carbonylation)
●Implicated in the mitochondrial theory of aging, neurodegenerative diseases (Parkinson’s, Alzheimer’s), and ischemia-reperfusion injury

ETC Inhibitors and Uncouplers

Understanding ETC inhibitors is clinically and pharmacologically critical. Each complex has specific inhibitors that block electron flow at distinct sites.

Target	Inhibitors	Mechanism
Complex I	Rotenone (insecticide), barbiturates, piericidin A	Block NADH → CoQ electron transfer
Complex II	Malonate	Competitive inhibitor of succinate (structural analog)
Complex III	Antimycin A	Blocks Q_i site of the Q cycle
Complex IV	Cyanide (CN⁻), carbon monoxide (CO), hydrogen sulfide (H₂S)	Bind to Fe³⁺ in cytochrome a₃, preventing O₂ reduction
ATP Synthase	Oligomycin	Blocks the proton channel of the F₀ subunit

Uncouplers (Dissipate Proton Gradient Without ATP Synthesis)

●2,4-Dinitrophenol (DNP) — lipophilic weak acid that carries H⁺ across the inner membrane. Used as a weight-loss drug in the 1930s but abandoned due to dangerous hyperthermia and fatalities
●FCCP (carbonyl cyanide-p-trifluoromethoxyphenylhydrazone) — potent synthetic uncoupler used extensively as a research tool
●Thermogenin (UCP1) — the physiological uncoupler found in brown adipose tissue. Provides a regulated proton leak for nonshivering thermogenesis in neonates and hibernating mammals. Energy released as heat rather than stored as ATP

Mitochondrial Diseases

Mitochondria possess their own genome — a circular, double-stranded mtDNA of 16,569 base pairs encoding 13 ETC subunits, 22 tRNAs, and 2 rRNAs. Inherited exclusively through maternal inheritance (sperm mitochondria are destroyed post-fertilization). Mutations in mtDNA cause a spectrum of disorders primarily affecting tissues with high energy demands.

MELAS

Mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes. Most commonly caused by an A3243G point mutation in tRNA^Leu.

MERRF

Myoclonic epilepsy with ragged red fibers. Caused by mutations in tRNA^Lys. Ragged red fibers on muscle biopsy reflect subsarcolemmal accumulation of abnormal mitochondria.

Leber Hereditary Optic Neuropathy (LHON)

Bilateral loss of central vision in young adults. Caused by Complex I mutations (e.g., G11778A in ND4 subunit). Most common mtDNA disease.

Kearns-Sayre Syndrome

Progressive external ophthalmoplegia, pigmentary retinopathy, and cardiac conduction defects. Caused by large-scale mtDNA deletions (typically 4,977 bp “common deletion”).

Key Concept: Heteroplasmy

Each cell contains hundreds to thousands of mitochondria, each with multiple copies of mtDNA. Heteroplasmy — the coexistence of mutant and wild-type mtDNA — determines disease severity. A threshold effect exists: symptoms typically manifest when the proportion of mutant mtDNA exceeds 60–90%, varying by tissue and mutation.

The Q Cycle (Complex III in Detail)

Complex III (cytochrome bc₁) operates via the Q cycle, an elegant mechanism proposed by Peter Mitchell that explains how electron transfer through this complex results in proton translocation across the inner membrane.

●Step 1 (Q_o site): CoQH₂ binds at the Q_o site (intermembrane space side). One electron is transferred to the Rieske iron-sulfur protein → cytochrome c₁ → cytochrome c (mobile carrier). The second electron goes to cytochrome b_L → cytochrome b_H → CoQ at the Q_i site (matrix side), forming a semiquinone (CoQ•⁻)
●Step 2 (second round): A second CoQH₂ binds at Q_o. Again, one electron goes to cytochrome c, and the other reduces the semiquinone at Q_i to CoQH₂ (picking up 2 H⁺ from the matrix)
●Net result per Q cycle: 1 CoQH₂ oxidized (net), 2 cytochrome c reduced, 4 H⁺ released to the intermembrane space (2 from CoQH₂ oxidation at Q_o + 2 from the proton-pumping mechanism per 2 electrons)
●This bifurcated electron flow explains how Complex III achieves a H⁺/e⁻ ratio of 2, making it an efficient proton pump despite having no direct proton-pumping machinery

🎥 Video Lecture: Electron Transport Chain

Ninja Nerd — Overview of all four ETC complexes, proton pumping, and ATP synthase

🎥 MIT 7.05: Oxidative Phosphorylation

MIT OCW — General Biochemistry (Spring 2020) — ETC complexes, chemiosmotic coupling, and ATP synthase

🎥 MIT 5.07SC: Redox Neutrality & NADH Shuttles

MIT OCW — Biological Chemistry I (Fall 2013) — Glycerol-3-phosphate and malate-aspartate shuttles

🎥 Harvard Online — Mitochondria & Oxidative Phosphorylation

🎥 Structure and Function of the Mitochondrion

Harvard Online — Cell Biology — Inner/outer membranes, cristae, matrix compartments, and the structural basis of chemiosmotic coupling

🎥 Mitochondria: The Cell's Powerhouse

Harvard Online — Cell Biology — How mitochondria generate ATP, the proton motive force, and energy transduction in living cells

🎥 Electron Transport Chain

Harvard Online — Cell Biology — Complexes I–IV, electron carriers, proton pumping, and the generation of the electrochemical gradient

🎥 ATP Synthase in Action

Harvard Online — Cell Biology — Rotary mechanism of ATP synthase (Complex V), F_o/F₁ subunits, and the molecular motor that produces ATP

🎥 Introduction to Mitochondrial Disease

Harvard Online — Cell Biology — How ETC defects cause disease: MELAS, MERRF, LHON, and the threshold effect of heteroplasmy

🎥 How Mitochondrial Diseases Are Inherited

Harvard Online — Cell Biology — Maternal inheritance of mtDNA, heteroplasmy, genetic bottleneck, and mitochondrial replacement therapy

Interactive: ETC Redox Potential Ladder & Energetics

Python

Visualization of standard reduction potentials and free energy changes along the electron transport chain

etc_energetics.py103 lines

import numpy as np
import matplotlib.pyplot as plt

# ============================================================
# Electron Transport Chain: Redox Potential Ladder & Energetics
# Visualizes the standard reduction potentials and free energy
# changes at each step of the mitochondrial ETC.
# ============================================================

# ETC components with standard reduction potentials (V)
components = [
    ('NADH/NAD+',     -0.320),
    ('Complex I\n(FMN)', -0.300),
    ('CoQ/CoQH2',    +0.045),
    ('Complex III\n(Cyt b)', +0.077),
    ('Cyt c',        +0.254),
    ('Complex IV\n(Cyt a3)', +0.385),
    ('O2/H2O',       +0.816),
]

# FADH2 enters at Complex II
fadh2_components = [
    ('FADH2/FAD',    +0.031),
    ('Complex II\n(FAD)', +0.031),
    ('CoQ/CoQH2',   +0.045),
]

names = [c[0] for c in components]
potentials = [c[1] for c in components]
positions = np.arange(len(components))

# Calculate delta G for each step
F = 96485  # C/mol (Faraday constant)
n = 2      # electrons transferred

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 7))

# --- Left panel: Redox potential ladder ---
ax1.barh(positions, potentials, height=0.6, color=['#818cf8' if p < 0 else '#34d399' for p in potentials],
         edgecolor='white', linewidth=0.5, alpha=0.8)

# Connect with arrows
for i in range(len(positions)-1):
    ax1.annotate('', xy=(potentials[i+1], positions[i+1]),
                xytext=(potentials[i], positions[i]),
                arrowprops=dict(arrowstyle='->', color='#f59e0b', lw=2))

# Protons pumped annotations
protons = [4, 0, 0, 4, 0, 2, 0]
for i, (name, pot) in enumerate(components):
    ax1.text(pot + 0.02 * (1 if pot >= 0 else -1), i, f'{pot:+.3f} V',
            va='center', ha='left' if pot >= 0 else 'right',
            fontsize=8, color='white', fontweight='bold')
    if protons[i] > 0:
        ax1.text(0.85, i, f'{protons[i]} H+ pumped',
                va='center', ha='right', fontsize=7, color='#fbbf24',
                bbox=dict(boxstyle='round,pad=0.2', facecolor='#422006', edgecolor='#f59e0b', alpha=0.8))

ax1.set_yticks(positions)
ax1.set_yticklabels(names, fontsize=9)
ax1.set_xlabel('Standard Reduction Potential E°\' (V)', fontsize=11)
ax1.set_title('ETC Redox Potential Ladder', fontsize=13, fontweight='bold')
ax1.axvline(x=0, color='gray', linestyle='--', alpha=0.5)
ax1.set_xlim(-0.5, 1.0)
ax1.grid(axis='x', alpha=0.2)
ax1.invert_yaxis()

# --- Right panel: Free energy diagram ---
# Cumulative energy released
delta_E = [potentials[i+1] - potentials[i] for i in range(len(potentials)-1)]
delta_G = [-n * F * dE / 1000 for dE in delta_E]  # kJ/mol

cumulative_G = [0]
for dg in delta_G:
    cumulative_G.append(cumulative_G[-1] + dg)

step_names = ['NADH', 'CI', 'CoQ', 'CIII', 'Cyt c', 'CIV', 'O2']
colors2 = ['#ef4444' if g < cumulative_G[i] else '#22c55e' for i, g in enumerate(cumulative_G[1:])]

for i in range(len(cumulative_G)-1):
    ax2.plot([i, i+1], [cumulative_G[i], cumulative_G[i+1]], 'o-', color='#f59e0b', linewidth=2.5, markersize=8)
    mid_y = (cumulative_G[i] + cumulative_G[i+1]) / 2
    ax2.text(i+0.5, mid_y + 5, f'{delta_G[i]:.1f}', ha='center', fontsize=8, color='white')

ax2.fill_between(range(len(cumulative_G)), cumulative_G, alpha=0.15, color='#f59e0b')
ax2.set_xticks(range(len(step_names)))
ax2.set_xticklabels(step_names, fontsize=9)
ax2.set_ylabel('Cumulative Free Energy Change (kJ/mol)', fontsize=11)
ax2.set_title('Energy Released Along the ETC', fontsize=13, fontweight='bold')
ax2.grid(True, alpha=0.2)

# Annotate total
ax2.text(3, cumulative_G[-1] + 10, f'Total: {cumulative_G[-1]:.0f} kJ/mol',
        fontsize=11, fontweight='bold', color='#f59e0b',
        bbox=dict(boxstyle='round,pad=0.3', facecolor='#1e293b', edgecolor='#f59e0b'))

plt.tight_layout()
plt.savefig('etc_energetics.png', dpi=150, bbox_inches='tight')
print("ETC Energetics visualization complete.")
print(f"Total free energy released (NADH -> O2): {cumulative_G[-1]:.1f} kJ/mol")
print(f"Theoretical maximum ATP (~30.5 kJ/mol each): {abs(cumulative_G[-1])/30.5:.1f} ATP")
print(f"Actual yield (with ~40% efficiency): ~2.5 ATP per NADH")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Understanding the ETC Energetics Plot

This visualization maps the thermodynamic landscape of the electron transport chain, showing how electrons “fall” through a series of redox carriers from NADH (E°′ = −0.32 V) to O₂ (E°′ = +0.82 V). The total potential difference of +1.14 Vcorresponds to a free energy release of −220 kJ/mol — enough to theoretically drive the synthesis of approximately 7 ATP molecules, though actual yield is ~2.5 ATP per NADH due to thermodynamic efficiency (~40%).

Left Panel — Redox Potential Ladder

Displays the standard reduction potential (E°′ in volts) at each step of the ETC. Electrons flow spontaneously from more negative to more positive potentials. The large drops at Complex I (ΔE ≈ 0.36 V), Complex III (ΔE ≈ 0.21 V), and Complex IV (ΔE ≈ 0.53 V) are precisely where proton pumping occurs. Complex II (succinate dehydrogenase) has a smaller ΔE and does not pump protons — this is why FADH₂ yields fewer ATP. The annotations show protons pumped per complex (4, 0, 4, 2 for Complexes I–IV respectively).

Right Panel — Cumulative Free Energy Release

Shows the progressive release of free energy (ΔG = −nFΔE) as electrons pass through each complex. The step function reveals that energy is released in discrete, large increments, not continuously. Each step releases enough energy to pump protons across the membrane. The horizontal dashed lines mark the free energy cost of ATP synthesis (~30.5 kJ/mol), showing visually how many ATP molecules each segment can theoretically produce. The cumulative total of −220 kJ/mol is parceled into three proton-pumping events, ultimately powering the rotary ATP synthase.

Why this matters: The ETC illustrates a fundamental principle — biological energy conversion occurs through coupled vectorial processes. The free energy from electron transfer doesn’t appear as heat; it is captured as a transmembrane proton gradient (Δp ≈ 200 mV), which Peter Mitchell called the proton-motive force. This gradient drives ATP synthase, one of the most remarkable molecular machines in biology.

Pentose Phosphate Pathway

Oxidative + Non-oxidative

Overview

The pentose phosphate pathway (PPP), also known as the hexose monophosphate shunt, branches from glycolysis at glucose-6-phosphate. It serves two critical functions: (1) generating NADPH for reductive biosynthesis and antioxidant defense, and (2) producing ribose-5-phosphate for nucleotide and nucleic acid synthesis.

Oxidative Phase (Irreversible)

$$\text{G6P} + 2\,\text{NADP}^+ + \text{H}_2\text{O} \longrightarrow \text{Ribulose-5-P} + 2\,\text{NADPH} + \text{CO}_2$$

G6P dehydrogenase: G6P → 6-phosphoglucono-δ-lactone (produces NADPH, rate-limiting)
Lactonase: Hydrolysis to 6-phosphogluconate
6-PG dehydrogenase: Oxidative decarboxylation → Ribulose-5-P (produces NADPH + CO₂)

Non-oxidative Phase (Reversible)

Transketolase and transaldolase interconvert 3-, 4-, 5-, 6-, and 7-carbon sugars, connecting the PPP back to glycolysis. The net result depends on cellular needs:

•Need NADPH + ribose: Run oxidative phase only
•Need NADPH only: Run both phases, recycle sugars to glycolysis
•Need ribose only: Run non-oxidative phase in reverse from glycolytic intermediates

Clinical Significance: G6PD Deficiency

Glucose-6-phosphate dehydrogenase deficiency is the most common enzyme deficiency worldwide (~400 million people). Without adequate NADPH, red blood cells cannot regenerate glutathione and are vulnerable to oxidative damage, leading to hemolytic anemia upon exposure to oxidative stressors (e.g., fava beans, certain drugs like primaquine).

NADPH: The Cell’s Reducing Currency

NADPH is structurally distinct from NADH (an extra phosphate on the 2’-OH of the adenine ribose) and serves an entirely different metabolic role. While NADH feeds electrons into the ETC for energy production, NADPH provides reducing equivalents for reductive biosynthesis and antioxidant defense. Cells maintain a high NADPH/NADP⁺ ratio (~100:1) in contrast to a low NADH/NAD⁺ ratio.

●Fatty acid synthesis: NADPH provides the electrons for each reductive step catalyzed by fatty acid synthase (FAS). Palmitate synthesis requires 14 NADPH molecules
●Cholesterol synthesis: HMG-CoA reductase, the rate-limiting enzyme of cholesterol biosynthesis, uses 2 NADPH per reaction
●Steroid hormone synthesis: Cytochrome P450 monooxygenases in adrenal cortex, gonads, and placenta require NADPH for hydroxylation reactions
●Glutathione reduction: Glutathione reductase uses NADPH to regenerate reduced glutathione (GSH) from its oxidized form (GSSG). GSH is critical for glutathione peroxidase, which detoxifies H₂O₂ and lipid peroxides, protecting cells from oxidative damage
●Nitric oxide synthesis: Nitric oxide synthase (NOS) uses NADPH to produce NO from arginine, essential for vasodilation and neurotransmission
●Respiratory burst in phagocytes: NADPH oxidase in neutrophils and macrophages generates superoxide (O₂•⁻) for microbial killing. Deficiency causes chronic granulomatous disease (CGD)
●Cytochrome P450 drug metabolism: Hepatic P450 enzymes use NADPH (via NADPH-cytochrome P450 reductase) for phase I biotransformation of drugs, toxins, and xenobiotics

Transketolase and Transaldolase Mechanisms

The non-oxidative phase of the PPP is catalyzed by two key transferase enzymes that rearrange carbon skeletons of sugar phosphates, connecting the PPP to glycolysis.

Transketolase (TK)

●Requires thiamine pyrophosphate (TPP) as a cofactor (vitamin B₁)
●Transfers 2-carbon units (glycolaldehyde) between sugar phosphates
●Catalyzes two reactions: C5 + C5 → C7 + C3, and C5 + C4 → C6 + C3
●TPP stabilizes the carbanion intermediate via its thiazolium ring

Transaldolase (TA)

●Does not require a cofactor
●Transfers 3-carbon units (dihydroxyacetone) between sugar phosphates
●Catalyzes: C7 + C3 → C4 + C6 (sedoheptulose-7-P + G3P → erythrose-4-P + fructose-6-P)
●Uses a Schiff base (lysine) mechanism to stabilize the intermediate

Together, the non-oxidative phase can interconvert sugar phosphates of 3, 4, 5, 6, and 7 carbons depending on cellular needs. When the cell needs more NADPH than ribose-5-P, excess ribose-5-P is recycled back to fructose-6-P and glyceraldehyde-3-P for glycolysis. When ribose-5-P is the priority, the non-oxidative phase runs in reverse.

Cancer Connection: Warburg Effect and the PPP

Rapidly proliferating cancer cells have enormous demands for both NADPH (for lipid biosynthesis and redox homeostasis) and ribose-5-phosphate (for nucleotide synthesis to support DNA replication). Cancer cells exploit the PPP in several ways:

●Upregulation of both glycolysis (Warburg effect) and the PPP provides the metabolic building blocks for rapid growth
●Transketolase-like protein 1 (TKTL1) is overexpressed in many cancers (colon, breast, lung, gastric), accelerating the non-oxidative branch and ribose production
●NADPH from the PPP supplies the reducing power for fatty acid synthase, fueling rapid membrane biogenesis in dividing cells
●NADPH also helps cancer cells resist oxidative stress from their hyperactive metabolism and during chemotherapy exposure

Wernicke-Korsakoff Syndrome

Thiamine (vitamin B₁) deficiency impairs transketolase activity in the non-oxidative phase of the PPP, as well as pyruvate dehydrogenase and α-ketoglutarate dehydrogenase in central metabolism. This condition is most commonly seen in chronic alcoholism due to poor dietary intake combined with impaired intestinal thiamine absorption and decreased hepatic storage.

Wernicke Encephalopathy (Acute)

●Confusion (global confusional state)
●Ataxia (cerebellar dysfunction)
●Ophthalmoplegia (cranial nerve VI palsy, nystagmus)

Korsakoff Psychosis (Chronic)

●Severe anterograde amnesia (inability to form new memories)
●Confabulation (fabrication of memories to fill gaps)
●Damage to mammillary bodies and medial thalamus (often irreversible)

Clinical Pearl

Always administer IV thiamine before glucose in suspected thiamine-deficient patients. Glucose administration increases thiamine demand (for pyruvate dehydrogenase), which can precipitate or worsen Wernicke encephalopathy in a depleted patient.

Regulation of the Pentose Phosphate Pathway

Glucose-6-phosphate dehydrogenase (G6PD) catalyzes the committed and rate-limiting step of the PPP. Its regulation is fundamentally different from the allosteric control seen in glycolysis.

●NADP⁺/NADPH ratio: The primary regulator. High cytoplasmic NADP⁺ concentration activates G6PD; high NADPH inhibits it. This provides elegant feedback — when NADPH is consumed (e.g., during fatty acid synthesis or oxidative stress), NADP⁺ accumulates and drives the pathway forward
●Insulin increases G6PD gene expression, upregulating the pathway in the fed state to supply NADPH for lipogenesis in liver and adipose tissue
●G6PD is not allosterically regulated by ATP or AMP, unlike phosphofructokinase-1 in glycolysis. The PPP responds to the redox state of the cell, not its energy charge
●The non-oxidative phase is regulated by substrate availability — it is freely reversible and driven by the concentrations of its sugar phosphate substrates and products

🎥 Video Lecture: Pentose Phosphate Pathway

Ninja Nerd — Detailed walkthrough of oxidative and non-oxidative phases, NADPH production, and clinical connections

🎥 MIT 7.05: Pentose Phosphate Pathway

MIT OCW — General Biochemistry (Spring 2020) — Oxidative and non-oxidative phases, NADPH production, and connections to nucleotide biosynthesis

Urea Cycle

5 Enzymes

Overview

The urea cycle (Krebs–Henseleit cycle) converts toxic ammonia (NH₃) from amino acid catabolism into urea for excretion by the kidneys. It operates across two compartments: the first two steps occur in the mitochondrial matrix, and the remaining three in the cytosol.

Overall Reaction

$$\text{NH}_3 + \text{CO}_2 + \text{Aspartate} + 3\,\text{ATP} + 2\,\text{H}_2\text{O} \longrightarrow \text{Urea} + \text{Fumarate} + 2\,\text{ADP} + \text{AMP} + 2\,P_i + \text{PP}_i$$

Net cost: 4 ATP equivalents per urea (since PPᵢ → 2 Pᵢ costs one additional high-energy bond)

The 5 Steps

Carbamoyl phosphate synthetase I (CPS I)Mitochondria

NH₃ + CO₂ + 2 ATP → Carbamoyl phosphate

Rate-limiting; activated by N-acetylglutamate

Ornithine transcarbamylase (OTC)Mitochondria

Carbamoyl phosphate + Ornithine → Citrulline

Citrulline exported to cytosol

Argininosuccinate synthetaseCytosol

Citrulline + Aspartate + ATP → Argininosuccinate

Second nitrogen enters from aspartate

Argininosuccinate lyaseCytosol

Argininosuccinate → Arginine + Fumarate

Fumarate links to TCA via aspartate-argininosuccinate shunt

ArginaseCytosol

Arginine + H₂O → Urea + Ornithine

Ornithine re-enters mitochondria to restart cycle

Clinical: Hyperammonemia

Deficiency in any urea cycle enzyme leads to hyperammonemia. OTC deficiency (X-linked) is the most common inherited form. Elevated ammonia is neurotoxic, causing cerebral edema, seizures, and coma. Treatment includes protein restriction, nitrogen scavengers (sodium benzoate, phenylbutyrate), and in severe cases, liver transplantation.

Nitrogen Balance and Amino Acid Catabolism

Before nitrogen enters the urea cycle, amino acids must first be deaminated through a two-step process that funnels nitrogen from diverse amino acid sources into a common pathway:

1. Transamination: Aminotransferases (also called transaminases) transfer the α-amino group from an amino acid to α-ketoglutarate, forming glutamate and a new α-keto acid. The cofactor for all aminotransferases is pyridoxal phosphate (PLP), the active form of vitamin B₆. Two clinically important aminotransferases are ALT (alanine aminotransferase), which converts alanine → pyruvate, and AST (aspartate aminotransferase), which converts aspartate → oxaloacetate. Both enzymes are elevated in hepatocellular damage (hepatitis, cirrhosis, drug-induced liver injury) and are routinely measured as part of liver function panels.

2. Oxidative Deamination: Glutamate dehydrogenase (GDH) removes the amino group from glutamate as free NH₄⁺, regenerating α-ketoglutarate. GDH is located in the mitochondrial matrix and can use either NAD⁺ or NADP⁺ as the electron acceptor. It is allosterically regulated: GTP inhibits (signaling energy sufficiency) and ADP activates (signaling energy deficit). This is the principal pathway for generating free ammonia that feeds into carbamoyl phosphate synthetase I (CPS-I), the first committed step of the urea cycle.

Glutamate as the Central Nitrogen Collector

Most amino acids funnel their α-amino nitrogen to glutamate via transamination. Glutamate then serves as the central hub for nitrogen disposal through two routes:

Route 1: Glutamate releases NH₄⁺ via GDH → this free ammonia enters step 1 of the urea cycle (CPS-I reaction), providing one of the two nitrogen atoms in urea.
Route 2: Glutamate transfers its nitrogen to oxaloacetate via AST to form aspartate → aspartate enters step 3 of the urea cycle (argininosuccinate synthetase), providing the second nitrogen atom in urea.

Thus, both nitrogen atoms in each molecule of urea ultimately derive from the same amino acid pool, but they arrive through two distinct metabolic routes — one as free ammonia and one as the amino group of aspartate.

The Glucose-Alanine Cycle

In skeletal muscle, amino groups from amino acid catabolism during exercise or fasting are transferred to pyruvate to form alanine (via ALT). Alanine is exported into the bloodstream and transported to the liver, where hepatic ALT regenerates pyruvate (which feeds into gluconeogenesis to produce glucose) and glutamate (which feeds nitrogen into the urea cycle). The newly synthesized glucose is released back into the blood and taken up by muscle, completing the cycle. This elegant interorgan shuttle serves two critical purposes: it transports nitrogen safely from muscle to liver without releasing toxic free ammonia into the circulation, and it provides carbon skeletons for hepatic gluconeogenesis to maintain blood glucose levels during fasting. The glucose-alanine cycle is conceptually analogous to the Cori cycle (which shuttles lactate), but handles nitrogen transport instead of anaerobic glycolysis end products.

The Aspartate-Argininosuccinate Shunt (Krebs Bicycle)

The urea cycle and the TCA cycle are intimately connected through a metabolic channeling pathway known as the aspartate-argininosuccinate shunt, sometimes called the Krebs bicycle. Fumarate produced in step 4 of the urea cycle (by argininosuccinate lyase) enters the TCA cycle, where it is hydrated to malate (by fumarase) and then oxidized to oxaloacetate (by malate dehydrogenase). Oxaloacetate is then transaminated by AST to regenerate aspartate, which re-enters the urea cycle at step 3 (argininosuccinate synthetase). This elegant interconnection coordinates carbon and nitrogen metabolism: each turn of the “bicycle” regenerates the aspartate needed for the urea cycle while simultaneously providing fumarate to the TCA cycle for energy production. The net effect is that the energy cost of the urea cycle (4 ATP equivalents per urea) is partially offset by the NADH generated when fumarate is processed through the TCA cycle segment.

N-Acetylglutamate (NAG): The Essential Activator

N-Acetylglutamate (NAG) is an obligatory allosteric activator of CPS-I (step 1 of the urea cycle). Without NAG, CPS-I is completely inactive and the entire urea cycle shuts down. NAG is synthesized by N-acetylglutamate synthase (NAGS) from glutamate and acetyl-CoA. Critically, NAGS is allosterically activated by arginine.

This creates a sophisticated feedforward regulation system: when amino acid catabolism increases (e.g., after a high-protein meal), rising levels of glutamate and arginine stimulate NAG production, which in turn activates CPS-I and upregulates the entire urea cycle to handle the increased nitrogen load.

Clinical: NAGS Deficiency

NAGS deficiency is the rarest inherited urea cycle disorder but is uniquely treatable with carglumic acid (Carbaglu), a synthetic analog of NAG that directly activates CPS-I. This is one of the few urea cycle defects that can be managed pharmacologically without liver transplantation.

Nitrogen Excretion Across Species

Different organisms have evolved distinct strategies for nitrogen waste disposal, determined primarily by water availability in their environment:

Ureotelic (mammals, adult amphibians): Excrete nitrogen as urea — requires energy (4 ATP per urea) but produces a water-soluble, non-toxic waste product that can be concentrated in urine.
Uricotelic (birds, reptiles, insects): Excrete nitrogen as uric acid — water-insoluble, excreted as a paste, requires minimal water. This is an adaptation to terrestrial and arid environments, and also allows nitrogen waste to be stored in the egg during embryonic development.
Ammonotelic (fish, aquatic invertebrates): Excrete nitrogen directly as NH₃ — the simplest strategy, requiring no biosynthetic energy, but necessitating large volumes of water for dilution due to ammonia's extreme toxicity.

Evolution selected the nitrogen excretion strategy based on the organism's access to water and its developmental constraints. Notably, human embryos are essentially ammonotelic (ammonia is cleared by maternal circulation), while tadpoles switch from ammonotelic to ureotelic during metamorphosis as they transition from aquatic to terrestrial life.

🎥 Video Lecture: Urea Cycle

Ninja Nerd — Complete urea cycle with enzymatic details and clinical pathology

🎥 MIT 7.05: Nitrogen & Amino Acid Metabolism

MIT OCW — General Biochemistry (Spring 2020) — Amino acid catabolism, transamination, deamination, and the urea cycle

Fatty Acid β-Oxidation & Synthesis

β-Oxidation

β-oxidation is the primary pathway for fatty acid catabolism. Long-chain fatty acyl-CoA molecules are transported into the mitochondrial matrix via the carnitine shuttle (CPT-I and CPT-II), then undergo a repeating 4-step spiral that removes two carbons per cycle as acetyl-CoA.

The β-Oxidation Spiral (4 Steps)

1. Oxidation

Acyl-CoA DH

FAD → FADH₂

2. Hydration

Enoyl-CoA hydratase

Adds H₂O

3. Oxidation

β-Hydroxyacyl-CoA DH

NAD⁺ → NADH

4. Thiolysis

Thiolase

Releases Acetyl-CoA

ATP Yield from Palmitate (C16:0)

$$\text{Palmitoyl-CoA} + 7\,\text{FAD} + 7\,\text{NAD}^+ + 7\,\text{CoA} + 7\,\text{H}_2\text{O} \longrightarrow 8\,\text{Acetyl-CoA} + 7\,\text{FADH}_2 + 7\,\text{NADH}$$

• 8 Acetyl-CoA × 10 ATP/turn = 80 ATP (via TCA + OxPhos)

• 7 FADH₂ × 1.5 ATP = 10.5 ATP

• 7 NADH × 2.5 ATP = 17.5 ATP

• Activation cost: −2 ATP (pyrophosphate cleavage)

Total: ~106 ATP per palmitate

Fatty Acid Synthesis

Fatty acid synthesis occurs in the cytosol and essentially reverses β-oxidation, but uses different enzymes and cofactors. The committed step is catalyzed by acetyl-CoA carboxylase (ACC), which converts acetyl-CoA to malonyl-CoA. The fatty acid synthase (FAS) complex then extends the chain by 2 carbons per cycle using NADPH as the reductant.

Key Regulation: Malonyl-CoA

Malonyl-CoA, the product of ACC and committed intermediate in fatty acid synthesis, is also the most potent inhibitor of CPT-I. This ensures that fatty acid synthesis and β-oxidation do not occur simultaneously — a classic example of reciprocal metabolic regulation.

The Carnitine Shuttle (Detailed Mechanism)

Long-chain fatty acyl-CoA molecules cannot cross the inner mitochondrial membrane directly. Instead, they rely on the carnitine shuttle, a three-step transport system:

Step 1 — CPT-I (carnitine palmitoyltransferase I, on the outer mitochondrial membrane): Transfers the acyl group from CoA to carnitine, forming acylcarnitine and releasing free CoA. CPT-I is the rate-limiting step of β-oxidation and is potently inhibited by malonyl-CoA.
Step 2 — CACT (carnitine-acylcarnitine translocase, spanning the inner membrane): Exchanges acylcarnitine (moving inward) for free carnitine (moving outward) across the inner mitochondrial membrane via an antiport mechanism.
Step 3 — CPT-II (carnitine palmitoyltransferase II, on the inner membrane, matrix side): Transfers the acyl group back from carnitine to mitochondrial CoA, regenerating acyl-CoA in the matrix for β-oxidation and releasing carnitine to recycle back out.

Clinical: Carnitine Deficiency

Primary carnitine deficiency (genetic defect in the carnitine transporter OCTN2) or secondary deficiency (caused by drugs such as valproate) leads to impaired fatty acid oxidation. Clinical manifestations include hypoketotic hypoglycemia (the body cannot oxidize fatty acids or produce ketone bodies, becoming entirely dependent on glucose), cardiomyopathy, and skeletal myopathy. Treatment involves oral L-carnitine supplementation.

Odd-Chain and Unsaturated Fatty Acid Oxidation

Odd-chain fatty acids (C15, C17, found in ruminant fat and some plant sources): The final round of β-oxidation yields propionyl-CoA (3 carbons) instead of the usual acetyl-CoA. Propionyl-CoA must be converted to succinyl-CoA to enter the TCA cycle through a three-step pathway:

Propionyl-CoA → D-methylmalonyl-CoA (propionyl-CoA carboxylase, requires biotin)
D-methylmalonyl-CoA → L-methylmalonyl-CoA (methylmalonyl-CoA epimerase)
L-methylmalonyl-CoA → succinyl-CoA (methylmalonyl-CoA mutase, requires vitamin B₁₂/adenosylcobalamin)

B₁₂ deficiency blocks this final step, causing methylmalonic acidemia — a diagnostically important finding that distinguishes B₁₂ deficiency from folate deficiency (both cause megaloblastic anemia, but only B₁₂ deficiency elevates methylmalonic acid).

Unsaturated fatty acids (e.g., oleic acid C18:1Δ9): These require additional auxiliary enzymes beyond the standard β-oxidation machinery. An enoyl-CoA isomerase converts cis-Δ³ double bonds to trans-Δ² double bonds (the form recognized by enoyl-CoA hydratase), while a 2,4-dienoyl-CoA reductase (NADPH-dependent) handles polyunsaturated fatty acids with conjugated double bond systems. Oxidation of unsaturated fatty acids yields slightly less FADH₂ because one or more dehydrogenation steps by acyl-CoA dehydrogenase are bypassed when pre-existing double bonds are present.

Ketone Body Synthesis (Ketogenesis)

When acetyl-CoA production from β-oxidation exceeds the TCA cycle's capacity to oxidize it (due to oxaloacetate depletion during fasting, starvation, or uncontrolled diabetes), the liver converts excess acetyl-CoA into ketone bodies. This occurs exclusively in the hepatic mitochondrial matrix:

2 Acetyl-CoA → Acetoacetyl-CoA (thiolase, running in reverse)
Acetoacetyl-CoA + Acetyl-CoA → HMG-CoA (HMG-CoA synthase — the rate-limiting enzyme of ketogenesis)
HMG-CoA → Acetoacetate + Acetyl-CoA (HMG-CoA lyase)
Acetoacetate → β-Hydroxybutyrate (β-hydroxybutyrate dehydrogenase, uses NADH)
Acetoacetate can also spontaneously decarboxylate → Acetone (responsible for the fruity breath odor in diabetic ketoacidosis)

Extrahepatic tissues (brain, heart, skeletal muscle) utilize ketone bodies by reversing the pathway: β-hydroxybutyrate → acetoacetate → acetoacetyl-CoA → 2 acetyl-CoA. The critical activation step uses succinyl-CoA:acetoacetate CoA transferase (also called thiophorase), an enzyme that the liver deliberately lacks — this ensures the liver exports ketone bodies rather than consuming them, maintaining fuel supply for other organs during fasting.

Clinical: Diabetic Ketoacidosis (DKA)

In uncontrolled type 1 diabetes, severe insulin deficiency leads to unrestrained lipolysis → massive free fatty acid release → overwhelming hepatic β-oxidation → excessive ketone body production (acetoacetate and β-hydroxybutyrate are acids) → metabolic acidosis with elevated anion gap. Patients present with Kussmaul breathing (deep, rapid respirations as respiratory compensation), fruity breath (acetone), dehydration, and altered mental status. Treatment requires IV fluids, insulin, and electrolyte replacement (especially potassium).

Fatty Acid Synthesis (Detailed Mechanism)

Fatty acid synthesis occurs in the cytosol (contrasting with β-oxidation in the mitochondrial matrix). It requires four key inputs:

Acetyl-CoA (transported from mitochondria to cytosol via the citrate shuttle: acetyl-CoA + OAA → citrate, exported, then cleaved by ATP-citrate lyase)
NADPH (provided by the pentose phosphate pathway and by malic enzyme: malate → pyruvate + CO₂ + NADPH)
CO₂ (for the carboxylation step; recycled, not net consumed)
Biotin (cofactor for acetyl-CoA carboxylase)

Acetyl-CoA carboxylase (ACC) is the rate-limiting enzyme, converting acetyl-CoA to malonyl-CoA. ACC is subject to elaborate multi-level regulation:

Allosteric activation: Citrate promotes polymerization of ACC into its active filamentous form (high citrate signals abundant acetyl-CoA and energy)
Allosteric inhibition: Palmitoyl-CoA (the end product) provides negative feedback inhibition
Covalent modification: AMPK phosphorylates ACC → inactive (during energy deficit). Insulin activates protein phosphatase 2A → dephosphorylates ACC → active (during the fed state)

Fatty acid synthase (FAS) is a large homodimeric enzyme (~273 kDa per subunit) that carries 7 distinct catalytic activities on a single polypeptide chain (condensation, reduction, dehydration, reduction again, plus ACP, malonyl/acetyl transferase, and thioesterase). Each cycle adds 2 carbons from malonyl-CoA. After 7 cycles, palmitate (C16:0) is released by the thioesterase domain. Further elongation (to C18, C20, etc.) and desaturation (introducing double bonds) occur on the endoplasmic reticulum via elongases and desaturases.

Hormonal Regulation: Fed vs. Fasted State

Fed state (insulin dominant):

Insulin activates ACC (via protein phosphatase), promoting malonyl-CoA production and fatty acid synthesis
Insulin activates lipoprotein lipase (LPL) in adipose tissue, promoting fatty acid uptake from circulating lipoproteins
Insulin promotes glucose uptake → glycolysis → acetyl-CoA → de novo lipogenesis
Malonyl-CoA levels are high → CPT-I is inhibited → β-oxidation is suppressed

Fasted state (glucagon/epinephrine dominant):

Glucagon activates hormone-sensitive lipase (HSL) in adipose tissue (via PKA phosphorylation) → lipolysis → free fatty acids released into blood
Glucagon activates AMPK → phosphorylates ACC → ACC inactive → malonyl-CoA levels fall
Low malonyl-CoA → CPT-I is active → β-oxidation proceeds at full capacity
Liver produces ketone bodies from excess acetyl-CoA to fuel the brain and peripheral tissues
Cortisol enhances lipolysis during prolonged fasting, sustaining fatty acid supply

Peroxisomal β-Oxidation

Very long chain fatty acids (VLCFA, >C20) cannot be directly processed by mitochondrial β-oxidation. Instead, they are first shortened in peroxisomes before the resulting medium-chain products are transferred to mitochondria for complete oxidation. Peroxisomal β-oxidation differs from the mitochondrial pathway in one critical respect: the first oxidation step uses a acyl-CoA oxidase that transfers electrons directly to O₂, producing H₂O₂ (which is then decomposed by catalase) instead of FADH₂. This makes peroxisomal β-oxidation less energy-efficient — the energy from that first oxidation step is lost as heat rather than captured for ATP synthesis.

Clinical: Peroxisomal Disorders

Zellweger syndrome (cerebrohepatorenal syndrome): An autosomal recessive disorder of peroxisome biogenesis — affected individuals have no functional peroxisomes, leading to VLCFA accumulation, severe neurological dysfunction, and death in infancy. X-linked adrenoleukodystrophy (X-ALD): Caused by mutations in the ABCD1 transporter gene, preventing VLCFA import into peroxisomes. VLCFA accumulate in adrenal cortex and central nervous system white matter, causing adrenal insufficiency and progressive demyelination. This condition was depicted in the film “Lorenzo's Oil” (1992), which explored an experimental dietary therapy using a mixture of oleic and erucic acid glycerol triesters to normalize plasma VLCFA levels.

🎥 Video Lecture: Fatty Acid β-Oxidation

Professor Dave Explains — Clear explanation of fatty acid oxidation, carnitine shuttle, and energy yield

🎥 MIT 7.05: Lipids & Fatty Acid Oxidation

MIT OCW — General Biochemistry (Spring 2020) — Beta-oxidation, carnitine shuttle, and energy accounting

🎥 MIT 7.05: Lipid Synthesis

MIT OCW — General Biochemistry (Spring 2020) — Fatty acid biosynthesis, ACC, FAS complex, and regulation

⚡Complete ATP Yield: From Glucose to CO₂

Interactive: Complete ATP Yield Calculator

Python

Calculate and visualize total ATP yield from glucose oxidation with P/O ratio sensitivity analysis

atp_yield_calculator.py110 lines

import numpy as np
import matplotlib.pyplot as plt

# ============================================================
# Complete ATP Yield Calculator
# Compares aerobic vs anaerobic glucose catabolism and performs
# a sensitivity analysis on P/O ratios.
# ============================================================

print("=" * 60)
print("COMPLETE AEROBIC OXIDATION OF GLUCOSE")
print("=" * 60)

# Standard P/O ratios
PO_NADH = 2.5  # ATP per NADH (via malate-aspartate shuttle)
PO_FADH2 = 1.5  # ATP per FADH2

# Glycolysis: Glucose -> 2 Pyruvate
glyc_ATP = 2       # net substrate-level
glyc_NADH = 2      # cytosolic NADH

# Pyruvate dehydrogenase: 2 Pyruvate -> 2 Acetyl-CoA
pdh_NADH = 2

# TCA cycle (x2 turns): 2 Acetyl-CoA -> 4 CO2
tca_NADH = 6       # 3 NADH per turn x 2
tca_FADH2 = 2      # 1 FADH2 per turn x 2
tca_GTP = 2        # 1 GTP per turn x 2 (= 2 ATP)

# Total electron carriers
total_NADH = glyc_NADH + pdh_NADH + tca_NADH  # 10 NADH
total_FADH2 = tca_FADH2                        # 2 FADH2

# ATP from oxidative phosphorylation
atp_from_nadh = total_NADH * PO_NADH
atp_from_fadh2 = total_FADH2 * PO_FADH2

# Total aerobic ATP
total_aerobic = glyc_ATP + tca_GTP + atp_from_nadh + atp_from_fadh2

print(f"\nGlycolysis:          {glyc_ATP} ATP (substrate-level)")
print(f"                     {glyc_NADH} NADH -> {glyc_NADH * PO_NADH:.1f} ATP")
print(f"Pyruvate DH:         {pdh_NADH} NADH -> {pdh_NADH * PO_NADH:.1f} ATP")
print(f"TCA cycle:           {tca_GTP} GTP (= {tca_GTP} ATP)")
print(f"                     {tca_NADH} NADH -> {tca_NADH * PO_NADH:.1f} ATP")
print(f"                     {tca_FADH2} FADH2 -> {tca_FADH2 * PO_FADH2:.1f} ATP")
print(f"\n{'─' * 40}")
print(f"TOTAL AEROBIC:       {total_aerobic:.1f} ATP per glucose")

print(f"\nANAEROBIC (fermentation): 2 ATP per glucose")
print(f"Efficiency ratio: {total_aerobic / 2:.0f}x more ATP aerobically")

# Free energy analysis
dG_glucose = -2870  # kJ/mol (complete oxidation)
dG_ATP = 30.5       # kJ/mol (hydrolysis under standard conditions)
efficiency = (total_aerobic * dG_ATP) / abs(dG_glucose) * 100
print(f"\nThermodynamic efficiency: {efficiency:.1f}%")
print(f"  (based on DG_glucose = {dG_glucose} kJ/mol, DG_ATP = {dG_ATP} kJ/mol)")

# ============ VISUALIZATION ============
fig, axes = plt.subplots(1, 3, figsize=(16, 6))

# Panel 1: ATP breakdown by source
sources = ['Glycolysis\n(substrate)', 'Glycolysis\n(NADH)', 'PDH\n(NADH)',
           'TCA\n(GTP)', 'TCA\n(NADH)', 'TCA\n(FADH2)']
values = [glyc_ATP, glyc_NADH * PO_NADH, pdh_NADH * PO_NADH,
          tca_GTP, tca_NADH * PO_NADH, tca_FADH2 * PO_FADH2]
colors = ['#22c55e', '#86efac', '#a78bfa', '#fbbf24', '#f59e0b', '#f97316']

bars = axes[0].bar(sources, values, color=colors, edgecolor='white', linewidth=0.5)
for bar, val in zip(bars, values):
    axes[0].text(bar.get_x() + bar.get_width()/2, bar.get_height() + 0.3,
                f'{val:.1f}', ha='center', fontsize=9, fontweight='bold')
axes[0].set_ylabel('ATP equivalents')
axes[0].set_title('ATP Yield by Source', fontweight='bold')
axes[0].set_ylim(0, max(values) * 1.3)
axes[0].tick_params(axis='x', labelsize=7)

# Panel 2: Aerobic vs Anaerobic
categories = ['Aerobic', 'Anaerobic']
atp_vals = [total_aerobic, 2]
bar_colors = ['#14b8a6', '#ef4444']
bars2 = axes[1].bar(categories, atp_vals, color=bar_colors, width=0.5, edgecolor='white')
for bar, val in zip(bars2, atp_vals):
    axes[1].text(bar.get_x() + bar.get_width()/2, bar.get_height() + 0.5,
                f'{val:.1f} ATP', ha='center', fontsize=11, fontweight='bold')
axes[1].set_ylabel('ATP per glucose')
axes[1].set_title('Aerobic vs Anaerobic', fontweight='bold')
axes[1].set_ylim(0, total_aerobic * 1.3)

# Panel 3: P/O ratio sensitivity
po_nadh_range = np.linspace(1.5, 3.0, 50)
po_fadh2_range = po_nadh_range - 1.0  # FADH2 always ~1 less than NADH

total_atp = glyc_ATP + tca_GTP + (glyc_NADH + pdh_NADH + tca_NADH) * po_nadh_range + tca_FADH2 * po_fadh2_range

axes[2].plot(po_nadh_range, total_atp, 'b-', linewidth=2.5, label='Total ATP')
axes[2].axvline(x=2.5, color='#f59e0b', linestyle='--', alpha=0.7, label='Standard P/O = 2.5')
axes[2].axhline(y=total_aerobic, color='#22c55e', linestyle=':', alpha=0.5)
axes[2].fill_between(po_nadh_range, total_atp, alpha=0.1, color='blue')
axes[2].set_xlabel('P/O ratio for NADH')
axes[2].set_ylabel('Total ATP per glucose')
axes[2].set_title('P/O Ratio Sensitivity', fontweight='bold')
axes[2].legend(fontsize=9)
axes[2].grid(True, alpha=0.2)

plt.tight_layout()
plt.savefig('atp_yield.png', dpi=150, bbox_inches='tight')
print("\nVisualization saved.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Understanding the ATP Yield Analysis

This comprehensive visualization breaks down the complete oxidation of one glucose molecule into its component ATP contributions, compares aerobic versus anaerobic metabolism, and performs a sensitivity analysis on the P/O ratio — the number of ATP molecules produced per pair of electrons passing through the ETC.

Panel 1 — ATP Yield Breakdown

A stacked bar chart showing the contribution from each metabolic stage: glycolysis (2 ATP + 2 NADH), the pyruvate dehydrogenase step (2 NADH), the TCA cycle (6 NADH + 2 FADH₂ + 2 GTP), and oxidative phosphorylation (converting all NADH and FADH₂ to ATP). The total of ~30–32 ATP demonstrates that over 90% of ATP comes from the ETC, not from direct substrate-level phosphorylation.

Panel 2 — Aerobic vs. Anaerobic

Contrasts the ~30–32 ATP from complete aerobic oxidation with the mere 2 ATP from anaerobic glycolysis. This 15–16× difference explains why aerobic organisms evolved to dominate complex ecological niches: muscle cells performing sustained work, neurons maintaining ion gradients, and cardiac myocytes beating continuously all depend on the dramatically higher yield of aerobic metabolism. Cancer cells (Warburg effect) sacrifice this efficiency for speed and biosynthetic intermediates.

Panel 3 — P/O Ratio Sensitivity

Explores how uncertainty in the P/O ratio affects total ATP yield. The classical value of 3.0 ATP/NADH (used in older textbooks) gives ~38 ATP, while the revised value of 2.5 ATP/NADH(reflecting the actual H⁺/ATP ratio of ~4.67 for ATP synthase, plus transport costs) gives ~30–32. This analysis shows that small changes in coupling efficiency have large effects on the total yield — a principle relevant to understanding uncoupling proteins, thermogenesis, and metabolic efficiency in obesity.

Why the number keeps changing: The “ATP per glucose” figure has evolved from 38 (1970s textbooks) to 36 to the current ~30–32 as our understanding of proton stoichiometry improved. The key variables are: (1) the H⁺/ATP ratio of ATP synthase (~4.67, not the integer 3 once assumed), (2) the cost of transporting ATP, ADP, and P_i across the mitochondrial membrane (~1 H⁺ per ATP), and (3) which shuttle (malate-aspartate vs. glycerol-3-phosphate) transfers cytoplasmic NADH.

🌐Multiscale Metabolic Simulation

This simulation connects three scales of biological organization: (1) molecular — Michaelis-Menten enzyme kinetics, (2) cellular — coupled ODE model of glycolysis, TCA cycle, and oxidative phosphorylation, and (3) tissue — how cellular ATP production affects a cardiac myocyte’s contractile frequency.

Interactive: Multiscale Metabolic Simulation

Python

Connects molecular enzyme kinetics (Michaelis-Menten) → cellular ATP production (coupled ODEs) → tissue function (cardiac beating rate)

multiscale_metabolism.py229 lines

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import solve_ivp

# ============================================================
# MULTISCALE METABOLIC SIMULATION
# Scale 1 (Molecular):  Michaelis-Menten enzyme kinetics
# Scale 2 (Cellular):   Coupled ODE for glycolysis + TCA + OxPhos
# Scale 3 (Tissue):     Cardiac myocyte beating rate vs ATP
# ============================================================

# --- Physical constants ---
F = 96485      # Faraday constant (C/mol)

# ========== SCALE 1: MOLECULAR (Enzyme Parameters) ==========
# Michaelis-Menten parameters from literature
enzymes = {
    'Hexokinase':    {'Vmax': 1.0, 'Km': 0.1},
    'PFK-1':         {'Vmax': 1.5, 'Km': 0.3},
    'Pyruvate DH':   {'Vmax': 2.0, 'Km': 0.15},
    'Citrate synth': {'Vmax': 1.8, 'Km': 0.05},
    'Complex I':     {'Vmax': 5.0, 'Km': 0.02},
    'ATP synthase':  {'Vmax': 8.0, 'Km': 0.1},
}

def mm(S, Vmax, Km):
    """Michaelis-Menten rate law."""
    return Vmax * S / (Km + S)

# ========== SCALE 2: CELLULAR (ODE System) ==========
def cell_metabolism(t, y, glucose_supply):
    Glc, Pyr, AcCoA, NADH, FADH2, ATP, ADP, O2 = y

# Clamp to non-negative
    Glc = max(Glc, 0)
    Pyr = max(Pyr, 0)
    AcCoA = max(AcCoA, 0)
    NADH = max(NADH, 0)
    FADH2 = max(FADH2, 0)
    ATP = max(ATP, 0.01)
    ADP = max(ADP, 0.01)
    O2 = max(O2, 0)

# NAD+ and FAD from conservation
    NAD_total = 2.0
    FAD_total = 0.5
    NADp = max(NAD_total - NADH, 0.001)
    FAD = max(FAD_total - FADH2, 0.001)

# Glycolysis (lumped): Glc -> 2 Pyr + 2 ATP + 2 NADH
    v_glyc = mm(Glc, enzymes['PFK-1']['Vmax'], enzymes['PFK-1']['Km']) * NADp / (NADp + 0.05) * ADP / (ADP + 0.5)

# Pyruvate dehydrogenase: Pyr -> AcCoA + NADH + CO2
    v_pdh = mm(Pyr, enzymes['Pyruvate DH']['Vmax'], enzymes['Pyruvate DH']['Km']) * NADp / (NADp + 0.05)

# TCA cycle (lumped): AcCoA -> 3 NADH + 1 FADH2 + 1 GTP
    v_tca = mm(AcCoA, enzymes['Citrate synth']['Vmax'], enzymes['Citrate synth']['Km']) * NADp / (NADp + 0.01)

# ETC + OxPhos: NADH -> NAD+ + ATP (P/O = 2.5)
    v_etc_nadh = mm(NADH, enzymes['Complex I']['Vmax'], enzymes['Complex I']['Km']) * O2 / (O2 + 0.005)

# FADH2 oxidation (P/O = 1.5)
    v_etc_fadh2 = 3.0 * FADH2 / (0.02 + FADH2) * O2 / (O2 + 0.005)

# Cellular ATP consumption (basal + contractile work)
    v_atp_use = 1.5 * ATP / (0.3 + ATP) + 0.5  # baseline demand

# O2 supply (from blood)
    v_o2_supply = 0.8 * (0.15 - O2) / (0.01 + abs(0.15 - O2))  # regulated toward 0.15 mM

dGlc   = glucose_supply - v_glyc
    dPyr   = 2 * v_glyc - v_pdh
    dAcCoA = v_pdh - v_tca
    dNADH  = 2 * v_glyc + v_pdh + 3 * v_tca - v_etc_nadh
    dFADH2 = v_tca - v_etc_fadh2
    dATP   = 2 * v_glyc + v_tca + 2.5 * v_etc_nadh + 1.5 * v_etc_fadh2 - v_atp_use
    dADP   = -dATP
    dO2    = v_o2_supply - 0.5 * v_etc_nadh - 0.5 * v_etc_fadh2

return [dGlc, dPyr, dAcCoA, dNADH, dFADH2, dATP, dADP, dO2]

# ========== SCALE 3: TISSUE (Cardiac Myocyte) ==========
def beating_rate(ATP_level):
    """
    Cardiac myocyte beating rate depends on ATP availability.
    Normal: ~70 bpm at [ATP] ~ 5 mM
    Below threshold: rate drops (heart failure model)
    """
    ATP_half = 2.0  # mM, half-maximal
    max_rate = 80   # bpm
    return max_rate * ATP_level**2 / (ATP_half**2 + ATP_level**2)

# ========== RUN SIMULATION ==========
# Initial conditions: [Glc, Pyr, AcCoA, NADH, FADH2, ATP, ADP, O2]
y0 = [5.0, 0.1, 0.05, 0.1, 0.01, 4.0, 1.0, 0.12]

# Simulate with varying glucose supply (fed -> fasting)
t_total = 60  # arbitrary time units
t_eval = np.linspace(0, t_total, 1000)

# Phase 1: Normal feeding (0-20), Phase 2: Reduced (20-40), Phase 3: Fasting (40-60)
def glucose_supply(t):
    if t < 20:
        return 0.8    # well-fed
    elif t < 40:
        return 0.3    # reduced intake
    else:
        return 0.05   # fasting/starvation

# We need to handle the piecewise glucose supply
# Solve in 3 phases
results_t = []
results_y = []
y_current = y0

for phase, (t_start, t_end, g_supply) in enumerate([
    (0, 20, 0.8), (20, 40, 0.3), (40, 60, 0.05)
]):
    t_phase = np.linspace(t_start, t_end, 334)
    sol = solve_ivp(
        lambda t, y: cell_metabolism(t, y, g_supply),
        (t_start, t_end), y_current, t_eval=t_phase,
        method='RK45', max_step=0.05,
        atol=1e-8, rtol=1e-6
    )
    results_t.append(sol.t)
    results_y.append(sol.y)
    y_current = sol.y[:, -1]

t_all = np.concatenate(results_t)
y_all = np.concatenate(results_y, axis=1)

# Compute tissue-level output
heart_rate = np.array([beating_rate(a) for a in y_all[5]])

# ========== VISUALIZATION ==========
fig, axes = plt.subplots(2, 3, figsize=(16, 10))
fig.suptitle('Multiscale Metabolic Simulation:\nMolecular → Cellular → Tissue',
             fontsize=14, fontweight='bold', y=0.98)

# Add phase shading to all panels
for ax in axes.flat:
    ax.axvspan(0, 20, alpha=0.08, color='green', label='_')
    ax.axvspan(20, 40, alpha=0.08, color='yellow', label='_')
    ax.axvspan(40, 60, alpha=0.08, color='red', label='_')

# Panel 1: Molecular Scale - Enzyme kinetics curves
S_range = np.linspace(0, 2, 200)
for name, params in list(enzymes.items())[:4]:
    axes[0,0].plot(S_range, mm(S_range, params['Vmax'], params['Km']), linewidth=2, label=name)
axes[0,0].set_xlabel('[Substrate] (mM)')
axes[0,0].set_ylabel('Reaction Rate (mM/s)')
axes[0,0].set_title('Scale 1: Molecular\n(Enzyme Kinetics)', fontweight='bold')
axes[0,0].legend(fontsize=7, loc='lower right')
axes[0,0].grid(True, alpha=0.2)

# Panel 2: Cellular metabolites
axes[0,1].plot(t_all, y_all[0], 'b-', label='Glucose', linewidth=2)
axes[0,1].plot(t_all, y_all[1], 'r-', label='Pyruvate', linewidth=2)
axes[0,1].plot(t_all, y_all[2], 'g-', label='Acetyl-CoA', linewidth=2)
axes[0,1].set_xlabel('Time (a.u.)')
axes[0,1].set_ylabel('Concentration (mM)')
axes[0,1].set_title('Scale 2: Cellular\n(Metabolite Dynamics)', fontweight='bold')
axes[0,1].legend(fontsize=8)
axes[0,1].grid(True, alpha=0.2)

# Panel 3: Redox state
axes[0,2].plot(t_all, y_all[3], color='purple', label='NADH', linewidth=2)
axes[0,2].plot(t_all, y_all[4], color='orange', label='FADH2', linewidth=2)
axes[0,2].plot(t_all, y_all[7], color='cyan', label='O2', linewidth=2)
axes[0,2].set_xlabel('Time (a.u.)')
axes[0,2].set_ylabel('Concentration (mM)')
axes[0,2].set_title('Redox & O2 State', fontweight='bold')
axes[0,2].legend(fontsize=8)
axes[0,2].grid(True, alpha=0.2)

# Panel 4: Energy charge (ATP/ADP)
axes[1,0].plot(t_all, y_all[5], 'r-', label='ATP', linewidth=2)
axes[1,0].plot(t_all, y_all[6], 'b-', label='ADP', linewidth=2)
ratio = y_all[5] / (y_all[5] + y_all[6] + 1e-10)
axes[1,0].plot(t_all, ratio * 5, 'k--', label='Energy charge (x5)', linewidth=1.5)
axes[1,0].set_xlabel('Time (a.u.)')
axes[1,0].set_ylabel('Concentration (mM)')
axes[1,0].set_title('Cellular Energy Status', fontweight='bold')
axes[1,0].legend(fontsize=8)
axes[1,0].grid(True, alpha=0.2)

# Panel 5: Tissue scale - Heart rate
axes[1,1].plot(t_all, heart_rate, color='#ef4444', linewidth=2.5)
axes[1,1].axhline(y=72, color='gray', linestyle=':', alpha=0.5, label='Normal (72 bpm)')
axes[1,1].set_xlabel('Time (a.u.)')
axes[1,1].set_ylabel('Heart Rate (bpm)')
axes[1,1].set_title('Scale 3: Tissue\n(Cardiac Beating Rate)', fontweight='bold')
axes[1,1].legend(fontsize=8)
axes[1,1].grid(True, alpha=0.2)
axes[1,1].set_ylim(0, 90)

# Panel 6: Phase diagram (ATP vs Glucose) - shows trajectory
axes[1,2].plot(y_all[0], y_all[5], 'b-', linewidth=1, alpha=0.7)
axes[1,2].plot(y_all[0,0], y_all[5,0], 'go', markersize=10, label='Start (fed)')
axes[1,2].plot(y_all[0,-1], y_all[5,-1], 'rs', markersize=10, label='End (fasting)')
# Color code by time
scatter = axes[1,2].scatter(y_all[0], y_all[5], c=t_all, cmap='coolwarm', s=5, zorder=5)
plt.colorbar(scatter, ax=axes[1,2], label='Time')
axes[1,2].set_xlabel('[Glucose] (mM)')
axes[1,2].set_ylabel('[ATP] (mM)')
axes[1,2].set_title('Phase Portrait\n(Glucose vs ATP)', fontweight='bold')
axes[1,2].legend(fontsize=8)
axes[1,2].grid(True, alpha=0.2)

# Add phase labels
fig.text(0.25, 0.01, 'Green: Fed | Yellow: Reduced | Red: Fasting',
         ha='center', fontsize=10, style='italic', color='gray')

plt.tight_layout(rect=[0, 0.03, 1, 0.96])
plt.savefig('multiscale_metabolism.png', dpi=150, bbox_inches='tight')

print("MULTISCALE SIMULATION RESULTS")
print("=" * 50)
print(f"\nScale 1 (Molecular): PFK-1 Km = {enzymes['PFK-1']['Km']} mM")
print(f"Scale 2 (Cellular):")
print(f"  Fed state:     [ATP] = {y_all[5, 100]:.2f} mM, [Glucose] = {y_all[0, 100]:.2f} mM")
print(f"  Reduced:       [ATP] = {y_all[5, 500]:.2f} mM, [Glucose] = {y_all[0, 500]:.2f} mM")
print(f"  Fasting:       [ATP] = {y_all[5, -1]:.2f} mM, [Glucose] = {y_all[0, -1]:.2f} mM")
print(f"Scale 3 (Tissue):")
print(f"  Fed heart rate:     {heart_rate[100]:.1f} bpm")
print(f"  Fasting heart rate: {heart_rate[-1]:.1f} bpm")
print(f"  Rate change:        {heart_rate[-1] - heart_rate[100]:+.1f} bpm")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Understanding the Multiscale Simulation

This simulation is unique because it bridges three levels of biological organization, demonstrating how a molecular event (enzyme kinetics) cascades through the cellular level (metabolite concentrations) to produce a physiological outcome (heart rate). The simulation models a cardiac myocyte transitioning from a fed state (high glucose) to a fasting state (reduced glucose supply), capturing the integrated metabolic response across scales.

Top Left — Molecular Scale: Enzyme Kinetics

Shows the Michaelis-Menten saturation curves for the three key regulatory enzymes (hexokinase, PFK-1, pyruvate kinase). The hyperbolic shape illustrates how reaction rate depends on substrate concentration. At low [S], rate increases nearly linearly; near saturation, increasing [S] barely affects rate. The K_m value (substrate concentration at half-maximal velocity) determines how sensitive each enzyme is to substrate availability — this is why PFK-1 (with its allosteric regulation) is the rate-limiting step.

Top Right — Cellular Scale: Metabolites

Tracks glucose, pyruvate, and ATP concentrations as external glucose supply decreases (simulating fasting). ATP levels fall as glucose becomes limiting, reflecting the cell’s inability to maintain glycolytic flux. The vertical dashed line marks the transition from fed to fasting state. Notice the buffering effect: ATP doesn’t collapse immediately because the cell adjusts flux rates, but eventually the decline becomes significant.

Bottom Left — Tissue Scale: Heart Rate

Maps cellular ATP levels to cardiac beating frequency using a Hill-type function. As ATP drops,heart rate decreases because the Ca²⁺ ATPase (SERCA), Na⁺/K⁺-ATPase, and myosin ATPase all require ATP. This demonstrates how a molecular perturbation (reduced PFK-1 flux) manifests as a tissue-level phenotype (bradycardia). In clinical medicine, this connects to why severe hypoglycemia or metabolic crises cause cardiac dysfunction.

Bottom Right — Scale Integration Summary

A bar chart summarizing the key output metrics at each scale, providing a snapshot comparison between the fed and fasted states. This panel encapsulates the central message of systems biology: understanding any single scale in isolation is insufficient. A change in K_m or V_max at the enzyme level propagates through metabolite pools to alter organ function, and vice versa.

Multiscale thinking in medicine: This framework explains why inborn errors of metabolism (single enzyme defects) cause systemic disease: a defective hexokinase isozyme affects glycolytic flux (molecular), which depletes cellular ATP (cellular), which impairs muscle contraction and neuronal function (tissue/organ). Modern pharmacology increasingly uses multiscale models to predict how drugs targeting a specific enzyme will affect whole-organism physiology.

🧬Connection with Molecular Biology

Cell Physiology and Molecular Biology are deeply interconnected. While Molecular Biology focuses on the molecular machinery of the cell—DNA, RNA, proteins, and their regulation—Cell Physiology examines how these molecules work together to produce cellular functions and behaviors.

🧬From Molecular Biology

→Ion Channel Structure: Molecular architecture determines gating and selectivity
→Receptor Proteins: GPCR and RTK structure underlies signal transduction
→Transport Proteins: Transporter mechanisms from molecular studies
→Enzyme Kinetics: Metabolic pathways with molecular basis

🔬To Cell Physiology

→Membrane Potential: Ion channels create bioelectric signals
→Signal Cascades: Receptor activation triggers cellular responses
→Active Transport: Pumps maintain cellular gradients
→Energy Production: Mitochondrial function powers the cell

🧬Explore Molecular Biology→⚡Ion Channels Deep Dive→

📋Prerequisites & Course Map

🗺️Learning Pathway

This course is designed for undergraduate students and above. The following foundational knowledge is strongly recommended before beginning the Cell Physiology course:

🧬