Part 6: Metabolism & Bioenergetics

ATP & Energy Currency

Every living cell requires a continuous supply of energy to maintain its structure, synthesize macromolecules, transport ions, and perform mechanical work. Adenosine triphosphate (ATP) serves as the universal energy currency — a remarkably versatile molecule that couples energy-releasing (exergonic) reactions to energy-requiring (endergonic) ones. Think of ATP as the "rechargeable battery" of the cell: it is constantly being discharged to ADP and recharged back to ATP, turning over its own weight roughly every 1–2 minutes.

Why Is ATP Hydrolysis Exergonic?

Four factors make ATP hydrolysis thermodynamically favorable:

ATP–ADP Cycle

The ATP–ADP cycle is the central mechanism by which cells transfer energy. Catabolic pathways (glycolysis, the TCA cycle, oxidative phosphorylation) regenerate ATP from ADP + Pᵢ, while biosynthetic reactions, active transport, and muscle contraction consume ATP. A resting human turns over approximately 40–70 kg of ATP per day — roughly equal to body weight — despite having only about 250 g of ATP at any instant.

Other Nucleoside Triphosphates

While ATP is the primary energy currency, cells also use GTP (protein synthesis, signal transduction), UTP (glycogen synthesis), and CTP (lipid synthesis). All are maintained at equilibrium with ATP via nucleoside diphosphate kinase: ATP + NDP ⇌ ADP + NTP (ΔG° ≈ 0).

Cellular Energy Charge

Daniel Atkinson introduced the concept of energy charge in the 1960s to describe the fraction of the total adenine nucleotide pool that carries "high-energy" phosphate bonds:

The energy charge acts as a metabolic "thermostat." Enzymes at critical branch points in metabolism respond to the ratio of ATP/ADP/AMP, ensuring that fuel consumption matches energy demand. For example:

• Phosphofructokinase-1 (PFK-1) — the key glycolytic enzyme — is inhibited by high ATP and activated by AMP
• Isocitrate dehydrogenase in the TCA cycle is activated by ADP and inhibited by ATP
• AMP-activated protein kinase (AMPK) senses rising AMP levels and switches on catabolic pathways while shutting down anabolic ones

import numpy as np
import matplotlib.pyplot as plt

# Calculate energy charge across different ATP:ADP:AMP ratios
# Simulating the energy charge spectrum
atp_fractions = np.linspace(0, 1, 100)
# For simplicity: assume ADP + AMP = 1 - ATP, split equally
adp_fractions = (1 - atp_fractions) * 0.6
amp_fractions = (1 - atp_fractions) * 0.4

total = atp_fractions + adp_fractions + amp_fractions
energy_charge = (atp_fractions + 0.5 * adp_fractions) / total

# Plot energy charge vs ATP fraction
fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 5))

# Left: Energy charge curve
ax1.plot(atp_fractions, energy_charge, 'b-', linewidth=2, label='Energy Charge')
ax1.axhline(y=0.85, color='green', linestyle='--', alpha=0.7, label='Normal Range (~0.85)')
ax1.axhspan(0.80, 0.95, alpha=0.1, color='green', label='Homeostatic Zone')
ax1.set_xlabel('ATP Fraction of Adenine Nucleotide Pool', fontsize=12)
ax1.set_ylabel('Energy Charge', fontsize=12)
ax1.set_title('Cellular Energy Charge', fontsize=14)
ax1.legend(fontsize=10)
ax1.set_ylim(0, 1.05)
ax1.grid(True, alpha=0.3)

# Right: Enzyme response to energy charge
ec_range = np.linspace(0, 1, 200)
# Catabolic enzymes (activated at low EC) — sigmoidal
catabolic_rate = 1 / (1 + np.exp(15 * (ec_range - 0.85)))
# Anabolic enzymes (activated at high EC) — sigmoidal
anabolic_rate = 1 / (1 + np.exp(-15 * (ec_range - 0.85)))

ax2.plot(ec_range, catabolic_rate, 'r-', linewidth=2, label='Catabolic (ATP-generating)')
ax2.plot(ec_range, anabolic_rate, 'b-', linewidth=2, label='Anabolic (ATP-consuming)')
ax2.axvline(x=0.85, color='green', linestyle='--', alpha=0.7, label='Normal EC (~0.85)')
ax2.set_xlabel('Energy Charge', fontsize=12)
ax2.set_ylabel('Relative Enzyme Activity', fontsize=12)
ax2.set_title('Enzyme Response to Energy Charge', fontsize=14)
ax2.legend(fontsize=10)
ax2.grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig('energy_charge_regulation.png', dpi=150, bbox_inches='tight')
plt.show()
print("Energy charge at typical cellular conditions: ~0.85-0.90")
print("ATP turnover in resting adult: ~40-70 kg/day")

NAD+ & FAD Redox Carriers

Metabolism is fundamentally about electron transfer. When fuel molecules (glucose, fatty acids, amino acids) are oxidized, the electrons stripped from carbon–hydrogen bonds don't go directly to oxygen. Instead, they are captured by specialized coenzymes — principally NAD⁺ and FAD — which act as mobile electron carriers, shuttling reducing equivalents from catabolic reactions to the electron transport chain in mitochondria.

Think of NAD⁺ and FAD as "electron taxis": they pick up high-energy electrons at metabolic reaction sites and deliver them to Complex I and Complex II of the electron transport chain, where the electrons' energy is used to pump protons and ultimately drive ATP synthesis.

NAD⁺/NADH: The Primary Electron Carrier

Nicotinamide adenine dinucleotide (NAD⁺) is derived from the vitamin niacin (B3). It accepts two electrons and one proton (a hydride ion, H⁻) at the nicotinamide ring:

NADP⁺/NADPH: The Biosynthetic Electron Carrier

NADP⁺ is structurally identical to NAD⁺ except for an extra phosphate group on the 2'-hydroxyl of the adenine ribose. This seemingly minor difference creates a crucial functional division:

• NADH — primarily donates electrons to the ETC for ATP production (catabolic)
• NADPH — provides reducing power for biosynthetic reactions (anabolic): fatty acid synthesis, cholesterol synthesis, and detoxification by cytochrome P450 enzymes

The pentose phosphate pathway (Section 5) is the primary source of NADPH in most cells.

FAD/FADH2

Flavin adenine dinucleotide (FAD) is derived from the vitamin riboflavin (B2). Unlike NAD⁺, FAD is typically covalently bound to the enzymes that use it (flavoproteins), and it can accept electrons one at a time, passing through a semiquinone radical intermediate:

Coenzyme Q and Cytochrome c

Beyond NAD⁺ and FAD, two additional mobile electron carriers operate within the ETC itself:

• Coenzyme Q (ubiquinone) — a lipid-soluble carrier embedded in the inner mitochondrial membrane; shuttles electrons from Complexes I and II to Complex III
• Cytochrome c — a small, water-soluble protein in the intermembrane space; transfers electrons from Complex III to Complex IV

Glycolysis

Glycolysis (from Greek glykys = sweet, lysis = splitting) is the ancient, universal pathway that splits one molecule of glucose (6 carbons) into two molecules of pyruvate (3 carbons each). It occurs in the cytoplasm of every living cell — from bacteria to neurons — and does not require oxygen, making it the cell's primary anaerobic energy source.

Glycolysis consists of 10 sequential enzymatic reactions divided into two phases: (1) the energy investment phase (steps 1–5), which consumes 2 ATP, and (2) the energy payoff phase (steps 6–10), which generates 4 ATP and 2 NADH. The net yield is 2 ATP + 2 NADH + 2 pyruvate per glucose.

Energy Investment Phase (Steps 1–5)

In the first five steps, the cell "invests" 2 ATP to phosphorylate and rearrange glucose, ultimately splitting it into two molecules of glyceraldehyde-3-phosphate (G3P):

Energy Payoff Phase (Steps 6–10)

In the second phase, each G3P is oxidized and converted to pyruvate, generating ATP and NADH. Since two G3P molecules are produced per glucose, all yields are doubled:

Fates of Pyruvate

Pyruvate sits at a critical metabolic crossroads. Its fate depends on oxygen availability and tissue type:

• Aerobic conditions: Pyruvate enters mitochondria → converted to acetyl-CoA by pyruvate dehydrogenase → enters TCA cycle
• Anaerobic conditions (muscle): Pyruvate → lactate (via lactate dehydrogenase), regenerating NAD⁺ to sustain glycolysis
• Anaerobic conditions (yeast): Pyruvate → ethanol + CO₂ (alcoholic fermentation)

Regulation of Glycolysis

Glycolysis is regulated at three irreversible steps, each catalyzed by enzymes with large negative ΔG values:

import numpy as np
import matplotlib.pyplot as plt

# Glycolysis energy profile — ΔG°' for each step
steps = ['Glc', 'G6P', 'F6P', 'F1,6BP', 'DHAP+\nG3P', 'G3P',
         '1,3BPG', '3PG', '2PG', 'PEP', 'Pyr']
# Cumulative standard free energy (kJ/mol) at each intermediate
delta_g_standard = [0, -16.7, -15.0, -29.2, -5.4, -5.4,
                    0.9, -17.6, -13.2, -5.7, -37.1]

# Actual cellular ΔG values (more negative due to concentrations)
delta_g_actual = [0, -33.4, -31.7, -56.8, -33.0, -33.0,
                  -29.7, -47.3, -43.8, -36.3, -73.3]

fig, ax = plt.subplots(figsize=(14, 6))
x = np.arange(len(steps))

ax.plot(x, delta_g_standard, 'bo-', linewidth=2, markersize=8,
        label="Standard ΔG°' (cumulative)")
ax.plot(x, delta_g_actual, 'rs-', linewidth=2, markersize=8,
        label="Actual ΔG (cellular conditions)")

# Highlight irreversible steps (1, 3, 10)
for idx in [1, 3, 10]:
    ax.axvline(x=idx, color='red', linestyle=':', alpha=0.3)
ax.annotate('Hexokinase\n(irreversible)', xy=(1, delta_g_actual[1]),
            xytext=(1.5, -15), fontsize=9, color='red',
            arrowprops=dict(arrowstyle='->', color='red'))
ax.annotate('PFK-1\n(committed step)', xy=(3, delta_g_actual[3]),
            xytext=(3.5, -25), fontsize=9, color='red',
            arrowprops=dict(arrowstyle='->', color='red'))
ax.annotate('Pyruvate kinase\n(irreversible)', xy=(10, delta_g_actual[10]),
            xytext=(8.5, -80), fontsize=9, color='red',
            arrowprops=dict(arrowstyle='->', color='red'))

# Phase labels
ax.axvspan(-0.5, 5.5, alpha=0.05, color='blue', label='Investment Phase')
ax.axvspan(5.5, 10.5, alpha=0.05, color='green', label='Payoff Phase')

ax.set_xticks(x)
ax.set_xticklabels(steps, fontsize=10)
ax.set_ylabel('Cumulative Free Energy Change (kJ/mol)', fontsize=12)
ax.set_title('Glycolysis Energy Profile: Standard vs Cellular Conditions', fontsize=14)
ax.legend(fontsize=10, loc='lower left')
ax.grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig('glycolysis_energy_profile.png', dpi=150, bbox_inches='tight')
plt.show()
print("Net yield per glucose: 2 ATP + 2 NADH + 2 pyruvate")
print("Three irreversible steps: hexokinase, PFK-1, pyruvate kinase")

Gluconeogenesis

Gluconeogenesis (literally "new glucose production") is the synthesis of glucose from non-sugar precursors — primarily lactate, glycerol, and glucogenic amino acids. It occurs mainly in the liver (90%) and kidney cortex (10%), and is essential for maintaining blood glucose during fasting, starvation, and prolonged exercise when glycogen reserves are depleted.

Gluconeogenesis is not simply the reverse of glycolysis. While seven of the ten glycolytic enzymes are shared (their reactions are near-equilibrium and freely reversible), the three irreversible glycolytic steps must be bypassed by different enzymes:

The Four Bypass Reactions

Gluconeogenic Precursors

Three major classes of molecules serve as glucose precursors:

• Lactate — produced by anaerobic glycolysis in muscle and red blood cells; converted to pyruvate by lactate dehydrogenase in the liver
• Glycerol — released from triglyceride lipolysis in adipose tissue; enters at the level of dihydroxyacetone phosphate (DHAP)
• Glucogenic amino acids — 18 of the 20 amino acids can be converted to glucose via their carbon skeletons (e.g., alanine → pyruvate, glutamate → α-ketoglutarate → oxaloacetate)

Cori Cycle

The Cori cycle (named after Carl and Gerty Cori, Nobel Prize 1947) is the metabolic cooperation between muscle and liver that recycles lactate:

Reciprocal Regulation of Glycolysis and Gluconeogenesis

Glycolysis and gluconeogenesis are reciprocally regulated to prevent futile cycling — the wasteful simultaneous operation of both pathways. The same allosteric effectors and hormonal signals that activate one pathway inhibit the other:

Pentose Phosphate Pathway

The pentose phosphate pathway (PPP) — also called the hexose monophosphate shunt — is an alternative route for glucose-6-phosphate oxidation that serves two critical purposes: (1) production of NADPH for reductive biosynthesis and antioxidant defense, and (2) synthesis of ribose-5-phosphate for nucleotide and nucleic acid biosynthesis. The PPP operates in the cytoplasm and is especially active in tissues with high biosynthetic demands (liver, adipose tissue, adrenal cortex, lactating mammary gland, red blood cells).

Oxidative Phase: NADPH Production

The oxidative phase is irreversible and generates NADPH by oxidizing glucose-6-phosphate to ribulose-5-phosphate:

Non-oxidative Phase

The non-oxidative phase is a series of reversible sugar interconversions catalyzed by transketolase (requires thiamine pyrophosphate, vitamin B1) and transaldolase. This phase connects the PPP to glycolysis by interconverting 3-, 4-, 5-, 6-, and 7-carbon sugars:

Metabolic Modes of the PPP

The cell can operate the PPP in different modes depending on which products it needs most:

import numpy as np
import matplotlib.pyplot as plt

# PPP flux distribution across different cell types
cell_types = ['Red Blood\nCells', 'Adipose\nTissue', 'Liver\n(Fed)', 'Bone\nMarrow', 'Adrenal\nCortex', 'Muscle']
nadph_demand = [90, 85, 75, 40, 95, 10]  # relative NADPH demand
r5p_demand = [5, 10, 20, 90, 5, 5]       # relative R5P demand
ppp_flux = [80, 70, 60, 55, 85, 3]       # % glucose through PPP

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

# Left: Stacked bar — NADPH vs R5P demand
x = np.arange(len(cell_types))
width = 0.35
ax1.bar(x - width/2, nadph_demand, width, label='NADPH Demand',
        color='#3B9797', alpha=0.8)
ax1.bar(x + width/2, r5p_demand, width, label='Ribose-5-P Demand',
        color='#BF092F', alpha=0.8)
ax1.set_xticks(x)
ax1.set_xticklabels(cell_types, fontsize=9)
ax1.set_ylabel('Relative Demand (%)', fontsize=12)
ax1.set_title('NADPH vs Ribose-5-P Demand by Tissue', fontsize=13)
ax1.legend(fontsize=10)
ax1.grid(axis='y', alpha=0.3)

# Right: PPP flux as % of glucose metabolism
colors = ['#BF092F' if f > 50 else '#16476A' for f in ppp_flux]
bars = ax2.barh(cell_types, ppp_flux, color=colors, alpha=0.8, edgecolor='white')
ax2.set_xlabel('% Glucose Metabolized via PPP', fontsize=12)
ax2.set_title('Pentose Phosphate Pathway Activity by Tissue', fontsize=13)
for bar, val in zip(bars, ppp_flux):
    ax2.text(bar.get_width() + 1, bar.get_y() + bar.get_height()/2,
             f'{val}%', va='center', fontsize=10, fontweight='bold')
ax2.set_xlim(0, 105)
ax2.grid(axis='x', alpha=0.3)

plt.tight_layout()
plt.savefig('ppp_tissue_distribution.png', dpi=150, bbox_inches='tight')
plt.show()
print("Oxidative phase: 2 NADPH + 1 CO2 per glucose-6-phosphate")
print("G6PD deficiency affects ~400 million people worldwide")

Metabolic Regulation

Metabolism is not a collection of independent reactions but a tightly coordinated network regulated at multiple levels to match energy supply with demand. Cells employ four principal regulatory strategies — substrate availability, allosteric regulation, covalent modification, and transcriptional control — operating on timescales from milliseconds to hours.

Levels of Metabolic Control

Key Allosteric Regulators

Several small molecules serve as universal metabolic "sensors" whose concentrations reflect the cell's energy state:

• ATP — signals energy sufficiency; inhibits catabolic enzymes (PFK-1, isocitrate dehydrogenase, pyruvate dehydrogenase)
• AMP — signals energy deficit; activates AMPK, PFK-1, and glycogen phosphorylase
• Citrate — signals biosynthetic precursor abundance; inhibits PFK-1, activates acetyl-CoA carboxylase
• Acetyl-CoA — signals fatty acid oxidation status; activates pyruvate carboxylase, inhibits pyruvate dehydrogenase
• NADH/NAD⁺ ratio — high ratio inhibits TCA cycle dehydrogenases and pyruvate dehydrogenase

Hormonal Regulation

Systemic metabolic coordination is achieved primarily through two opposing pancreatic hormones:

The Fructose-2,6-Bisphosphate Switch

The bifunctional enzyme PFK-2/FBPase-2 is a masterstroke of metabolic engineering. This single protein has two catalytic domains: a kinase domain (PFK-2) that makes F2,6BP, and a phosphatase domain (FBPase-2) that degrades it. The balance between these activities is controlled by phosphorylation:

• Insulin (fed) → dephosphorylation → PFK-2 active → ↑ F2,6BP → ↑ glycolysis
• Glucagon (fasted) → PKA phosphorylation → FBPase-2 active → ↓ F2,6BP → ↑ gluconeogenesis

This elegant mechanism allows a single phosphorylation event to simultaneously accelerate one pathway and brake the other — the quintessential reciprocal regulatory switch.

import numpy as np
import matplotlib.pyplot as plt

# Metabolic pathway regulation under fed vs fasted states
pathways = ['Glycolysis', 'Gluconeogenesis', 'Glycogen\nSynthesis',
            'Glycogenolysis', 'Lipogenesis', 'Lipolysis',
            'β-Oxidation', 'Ketogenesis']
fed_activity = [90, 10, 85, 15, 80, 10, 20, 5]      # % activity in fed state
fasted_activity = [20, 85, 10, 90, 5, 85, 90, 70]    # % activity in fasted state

fig, ax = plt.subplots(figsize=(12, 6))
x = np.arange(len(pathways))
width = 0.35

bars1 = ax.bar(x - width/2, fed_activity, width, label='Fed State (Insulin)',
               color='#3B9797', alpha=0.85, edgecolor='white')
bars2 = ax.bar(x + width/2, fasted_activity, width, label='Fasted State (Glucagon)',
               color='#BF092F', alpha=0.85, edgecolor='white')

ax.set_xticks(x)
ax.set_xticklabels(pathways, fontsize=10, rotation=15, ha='right')
ax.set_ylabel('Relative Pathway Activity (%)', fontsize=12)
ax.set_title('Hormonal Control of Metabolism: Fed vs Fasted States', fontsize=14)
ax.legend(fontsize=11, loc='upper right')
ax.set_ylim(0, 105)
ax.grid(axis='y', alpha=0.3)

# Add value labels
for bar in bars1:
    ax.text(bar.get_x() + bar.get_width()/2, bar.get_height() + 1,
            f'{int(bar.get_height())}%', ha='center', fontsize=8, fontweight='bold')
for bar in bars2:
    ax.text(bar.get_x() + bar.get_width()/2, bar.get_height() + 1,
            f'{int(bar.get_height())}%', ha='center', fontsize=8, fontweight='bold')

plt.tight_layout()
plt.savefig('metabolic_regulation_fed_fasted.png', dpi=150, bbox_inches='tight')
plt.show()
print("Insulin: store fuel (anabolic)")
print("Glucagon: mobilize fuel (catabolic)")
print("F2,6BP is the most potent regulator of glycolysis vs gluconeogenesis")

Practice Exercises

Metabolism & Bioenergetics Worksheet

Conclusion & Next Steps

In this article, we explored the foundational machinery of cellular energy metabolism — from ATP as the universal energy currency and NAD⁺/FAD as electron shuttles, through the complete glycolysis pathway, the reverse pathway of gluconeogenesis, and the versatile pentose phosphate pathway. We also examined how cells regulate metabolic flux through allosteric effectors, covalent modification, and hormonal signaling via insulin and glucagon.

Key takeaways include: (1) ATP is continuously recycled at extraordinary rates (~40-70 kg/day in a resting adult), (2) glycolysis is regulated at three irreversible steps with PFK-1 as the primary gatekeeper, (3) glycolysis and gluconeogenesis are reciprocally regulated to prevent futile cycling, (4) the PPP provides NADPH and ribose-5-phosphate with tissue-specific flux modes, and (5) the insulin/glucagon axis coordinates metabolism across the entire organism.