5. Glycolysis

Overview of Glycolysis

Glycolysis (from Greek glykys = sweet, lysis = splitting) is the metabolic pathway that converts glucose ($\text{C}_6\text{H}_{12}\text{O}_6$) into two molecules of pyruvate ($\text{C}_3\text{H}_3\text{O}_3^-$). This pathway occurs in the cytoplasm of virtually all living cells, making it one of the most ancient and conserved metabolic sequences.

The overall equation for glycolysis is:

Glycolysis is divided into two phases:

The Ten Steps of Glycolysis

Energy Investment Phase

Energy Payoff Phase

From this point, each reaction occurs twice per glucose molecule (two GAP molecules enter).

Derivation 1: Thermodynamics of Glycolysis

The distinction between standard free energy change ($\Delta G°'$) and actual cellular free energy change ($\Delta G$) is critical for understanding pathway regulation. The relationship is:

where $Q$ is the mass-action ratio of actual cellular concentrations:

Near-Equilibrium vs Far-from-Equilibrium Steps

At equilibrium, $\Delta G = 0$ and $Q = K_{eq}$. A reaction is near equilibrium when $Q \approx K_{eq}$, meaning $|\Delta G| \ll RT$. These reactions respond rapidly to perturbations and are freely reversible.

A reaction is far from equilibrium when $Q \ll K_{eq}$, giving a large negative $\Delta G$. These steps are effectively irreversible under physiological conditions and serve as regulatory control points.

Example Calculation: PFK-1 (Step 3)

Given $\Delta G°' = -14.2\;\text{kJ/mol}$, at $T = 310\;\text{K}$ (body temperature):

Under cellular conditions with typical concentrations ([F6P] = 0.08 mM, [ATP] = 3.0 mM, [F1,6BP] = 0.02 mM, [ADP] = 0.9 mM):

This is close to the reported $\Delta G = -22.2\;\text{kJ/mol}$, confirming that PFK-1 operates far from equilibrium.

Overall Thermodynamics

Summing the standard free energies:

Under cellular conditions, the actual free energy change is much more negative:

Derivation 2: Allosteric Regulation of PFK-1

Phosphofructokinase-1 (PFK-1) is the committed step and the most important regulatory enzyme in glycolysis. It is a tetrameric enzyme that exhibits sigmoidal kinetics characteristic of allosteric enzymes.

The Hill Equation for Cooperative Binding

The response of an allosteric enzyme to substrate concentration is described by the Hill equation:

where $n_H$ is the Hill coefficient measuring cooperativity, and $K_{0.5}$ is the substrate concentration at half-maximal velocity. Taking the logarithm:

A plot of $\log(v/(V_{\max}-v))$ vs $\log[S]$ (the Hill plot) gives a straight line with slope $n_H$. For PFK-1, $n_H \approx 3\text{-}4$, indicating strong positive cooperativity.

Sensitivity Analysis

The key biological advantage of cooperativity is the switch-like response. The ratio of substrate concentrations needed to go from 10% to 90% activity is:

For a Michaelis-Menten enzyme ($n_H = 1$): $R_s = 81$. For PFK-1 ($n_H = 4$): $R_s = 81^{1/4} = 3.0$. This means PFK-1 switches from 10% to 90% activity over only a 3-fold change in substrate concentration, enabling tight metabolic control.

Allosteric Effectors of PFK-1

Fructose-2,6-Bisphosphate: The Master Regulator

F2,6BP is synthesized and degraded by the bifunctional enzyme PFK-2/FBPase-2. In the liver:

• Fed state (insulin, high F6P): PFK-2 active $\rightarrow$ [F2,6BP]$\uparrow$ $\rightarrow$ glycolysis ON
• Fasted state (glucagon $\rightarrow$ cAMP $\rightarrow$ PKA): FBPase-2 active $\rightarrow$ [F2,6BP]$\downarrow$ $\rightarrow$ glycolysis OFF

Mathematically, F2,6BP modulates the apparent $K_{0.5}$ for F6P. In the presence of F2,6BP:

where $K_a$ is the activation constant for F2,6BP binding ($K_a \approx 0.1\;\mu\text{M}$). Even micromolar concentrations of F2,6BP dramatically shift the dose-response curve to the left, activating PFK-1 at physiological F6P concentrations.

Derivation 3: Energy Yield of Glycolysis

ATP Accounting

Careful bookkeeping of ATP consumed and produced in glycolysis:

NADH Production

Step 6 (GAPDH) produces 1 NADH per GAP. Since two GAP molecules are formed per glucose:

Free Energy Efficiency

The standard free energy of complete glucose oxidation is:

The free energy captured in glycolysis (under standard conditions):

The free energy released in going from glucose to 2 pyruvate is $\Delta G°' = -35.7\;\text{kJ/mol}$ (standard) or $\approx -74\;\text{kJ/mol}$ (cellular). Thus the overall process is strongly exergonic, driving it forward irreversibly under physiological conditions.

Derivation 4: Fermentation Pathways

Under anaerobic conditions, the NADH produced in step 6 cannot be reoxidized by the electron transport chain. To sustain glycolysis, NAD$^+$ must be regenerated by fermentation.

Lactic Acid Fermentation (Mammals)

In anaerobic muscle and erythrocytes, pyruvate is reduced to lactate by lactate dehydrogenase (LDH):

The overall reaction of anaerobic glycolysis becomes:

Note that NAD$^+$ and NADH cancel out — they are recycled. The net yield is only 2 ATP per glucose, with no net change in NAD$^+$/NADH.

Ethanol Fermentation (Yeast)

In yeast under anaerobic conditions, pyruvate is converted to ethanol in two steps:

The overall equation:

NAD$^+$ Regeneration: The Critical Constraint

The fundamental purpose of fermentation is maintaining the NAD$^+$/NADH ratio. The cytoplasmic pool of NAD$^+$ is limited ($\approx 10^{-5}\;\text{M}$), and glycolysis would halt within seconds if NADH were not reoxidized:

In aerobic conditions, NADH is reoxidized via shuttle systems (malate-aspartate shuttle or glycerol-3-phosphate shuttle) that transfer electrons to mitochondrial ETC. Each NADH produces approximately 2.5 ATP via oxidative phosphorylation, making aerobic metabolism ~15-fold more efficient than fermentation alone.

Enzyme Mechanisms in Detail

Hexokinase: Induced Fit Mechanism

Hexokinase undergoes a dramatic conformational change upon glucose binding. The enzyme's two lobes close around the substrate like a clamp, excluding water from the active site. This induced fit mechanism is essential because:

• It prevents the wasteful hydrolysis of ATP by water
• It positions the $\gamma$-phosphate of ATP adjacent to the C-6 hydroxyl of glucose
• It creates a hydrophobic environment that favors phosphoryl transfer

The Mg$^{2+}$-ATP complex is the true substrate. The metal ion coordinates the $\beta$ and $\gamma$ phosphates, neutralizing the negative charge and facilitating nucleophilic attack by the glucose C-6 oxygen:

In the liver, the isozyme glucokinase (hexokinase IV) replaces hexokinase. Glucokinase has a higher $K_m$ ($\approx 10\;\text{mM}$ vs $0.1\;\text{mM}$) and is not inhibited by its product G6P, allowing the liver to act as a glucose buffer.

GAPDH: Coupled Oxidation and Phosphorylation

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) catalyzes a remarkable reaction that couples oxidation of an aldehyde to phosphorylation, preserving the oxidation energy in a high-energy acyl phosphate bond:

The thioester intermediate is key: its hydrolysis has $\Delta G°' = -32\;\text{kJ/mol}$, comparable to ATP hydrolysis. By capturing this energy as an acyl phosphate in 1,3-BPG ($\Delta G°'_{\text{hydrolysis}} = -49.3\;\text{kJ/mol}$), the energy is conserved and can drive ATP synthesis in the next step.

Pyruvate Kinase: Substrate-Level Phosphorylation

The conversion of PEP to pyruvate proceeds via an enolpyruvate intermediate:

The large negative $\Delta G°'$ of this reaction ($-31.4\;\text{kJ/mol}$) is driven by the tautomerization of enolpyruvate to the keto form. The enol-to-keto conversion is thermodynamically very favorable ($\Delta G°' \approx -46\;\text{kJ/mol}$), which is why PEP has such a high phosphoryl transfer potential ($\Delta G°' = -61.9\;\text{kJ/mol}$ for hydrolysis).

Pyruvate kinase requires both Mg$^{2+}$ (coordinates phosphoryl groups) and K$^+$ (stabilizes the transition state). It is allosterically regulated:

• Activated by: F1,6BP (feedforward activation from PFK-1)
• Inhibited by: ATP, alanine (signals amino acid abundance)
• Liver isozyme (L-PK): Inactivated by PKA phosphorylation (glucagon signaling)

Metabolic Integration and Clinical Connections

The Pasteur Effect

Louis Pasteur observed that yeast consume far more glucose under anaerobic conditions than aerobic conditions. This Pasteur effect is explained by the energy yield difference:

Under aerobic conditions, oxidative phosphorylation produces ~32 ATP per glucose (vs. 2 from fermentation alone), so 16-fold less glucose is needed to meet the same ATP demand. The Pasteur effect is mediated by allosteric inhibition of PFK-1 by ATP and citrate when oxidative phosphorylation is active.

The Warburg Effect in Cancer

Otto Warburg (Nobel Prize 1931) observed that cancer cells preferentially use glycolysis even in the presence of oxygen — so-called aerobic glycolysis or the Warburg effect. Cancer cells exhibit:

• Upregulated glucose transporters (GLUT1, GLUT3)
• Increased hexokinase II expression (bound to mitochondrial outer membrane)
• Expression of PKM2 (a less active pyruvate kinase isozyme) that diverts intermediates to biosynthesis
• High lactate production even in normoxic conditions

This seemingly inefficient strategy provides cancer cells with biosynthetic precursors (G6P for pentose phosphate pathway, 3PG for serine, pyruvate for alanine) needed for rapid cell division. The clinical application is FDG-PET imaging, which uses $^{18}$F-fluorodeoxyglucose to detect tumors based on their elevated glucose uptake.

Arsenic Poisoning and Glycolysis

Arsenate (AsO$_4^{3-}$) is structurally analogous to phosphate and can substitute for P$_i$ in the GAPDH reaction. The product, 1-arseno-3-phosphoglycerate, is spontaneously hydrolyzed (unstable), bypassing the PGK step and eliminating one substrate-level phosphorylation:

This arsenolysis reduces the net ATP yield of glycolysis from 2 to 0, effectively uncoupling oxidation from phosphorylation in the payoff phase.

Erythrocyte Glycolysis and 2,3-BPG

Red blood cells lack mitochondria and depend entirely on glycolysis for ATP. They possess a unique bypass called the Rapoport-Luebering pathway:

2,3-BPG is a critical allosteric effector of hemoglobin, decreasing O$_2$ affinity and promoting oxygen release in tissues. This pathway sacrifices one ATP per 1,3-BPG diverted (bypassing PGK) in exchange for regulating oxygen delivery. The concentration of 2,3-BPG in erythrocytes ($\approx 5\;\text{mM}$) is nearly equimolar with hemoglobin.

Genetic Deficiencies

Glycolytic enzyme deficiencies are rare autosomal recessive disorders that predominantly affect erythrocytes (which have no alternative energy source):

Thermodynamic Deep Dive: Mass Action Ratios

For each glycolytic step, we can calculate the mass action ratio ($Q$) and compare it to $K_{eq}$. The ratio $\Gamma = Q/K_{eq}$ reveals how far from equilibrium each reaction operates:

When $\Gamma \approx 1$: the reaction is near equilibrium ($\Delta G \approx 0$), readily reversible, and the enzyme operates near its maximum capacity in both directions.

When $\Gamma \ll 1$: the reaction is far from equilibrium ($\Delta G \ll 0$), effectively irreversible, and the enzyme has reserve capacity — it is a potential regulatory point.

Flux Control Coefficients

Metabolic control analysis (MCA) quantifies each enzyme's contribution to overall pathway flux through flux control coefficients ($C^J_i$):

where $J$ is the pathway flux and $e_i$ is the enzyme concentration. The summation theorem states:

In glycolysis, PFK-1 has the largest flux control coefficient ($C^J_{\text{PFK}} \approx 0.5\text{-}0.8$ depending on conditions), confirming its role as the primary rate-controlling step. However, control is shared — no single enzyme has $C^J = 1$.

Elasticity and Connectivity

The elasticity coefficient ($\varepsilon^v_S$) describes how sensitive an enzyme's rate is to metabolite concentration changes:

For near-equilibrium enzymes, $|\varepsilon| \rightarrow \infty$, meaning they respond dramatically to small perturbations (buffering metabolite concentrations). For far-from-equilibrium enzymes, $\varepsilon$ is smaller and determined by the kinetic properties (allosteric regulation). The connectivity theorem links flux control and elasticity:

This fundamental relationship shows that enzymes with high flux control coefficients must have low elasticities toward shared metabolites, and vice versa — a powerful constraint on metabolic design.

Tissue-Specific Adaptations of Glycolysis

Brain

The brain consumes approximately $120\;\text{g}$ of glucose per day ($\approx 60\%$ of resting glucose consumption), relying almost exclusively on aerobic glycolysis followed by complete oxidation. The brain has limited glycogen stores ($\approx 5\;\text{g}$), making it highly vulnerable to hypoglycemia. Below $\approx 2.5\;\text{mM}$ blood glucose, neurological symptoms appear because:

The brain GLUT3 transporter has $K_m \approx 1.4\;\text{mM}$, so at 5 mM blood glucose the transporter operates at ~78% capacity. During prolonged starvation ($>$3 days), the brain adapts to using ketone bodies (acetoacetate, $\beta$-hydroxybutyrate) for up to 75% of its energy needs.

Erythrocytes

Red blood cells lack mitochondria entirely and depend 100% on anaerobic glycolysis. They produce $\approx 2\;\text{moles}$ of lactate per mole of glucose consumed. The unique Rapoport-Luebering bypass for 2,3-BPG production has already been discussed. Additional features include:

• Hexokinase I (low $K_m$, inhibited by G6P) ensures glucose is phosphorylated even at low blood glucose
• The NADH/NAD$^+$ ratio is maintained by LDH, which is essential for continuous glycolysis
• Pentose phosphate pathway provides NADPH for glutathione reduction, protecting against oxidative damage

Skeletal Muscle

Muscle contains large glycogen stores ($\approx 400\;\text{g}$) that provide fuel for rapid, high-intensity contraction. Key features:

• Hexokinase II has higher $V_{\max}$ than HK-I, matching muscle's higher glycolytic capacity
• PFK-1 in muscle (M-type subunit) is less sensitive to citrate inhibition than the liver (L-type), allowing glycolysis to proceed even during active TCA cycle
• Muscle lacks glucose-6-phosphatase, so muscle glycogen cannot contribute to blood glucose
• The LDH-M$_4$ isozyme favors pyruvate $\rightarrow$ lactate (anaerobic conditions)

Liver

The liver is the master glucose regulator, switching between glucose consumption (glycolysis) and glucose production (gluconeogenesis) depending on the hormonal milieu. Key features:

• Glucokinase (HK-IV) replaces hexokinase, with high $K_m$ ($\approx 10\;\text{mM}$) and no product inhibition
• Contains both PFK-1 and FBPase-1 for bidirectional flux control via F2,6BP
• L-type pyruvate kinase is phosphorylated and inactivated by glucagon (via PKA)
• The LDH-H$_4$ isozyme favors lactate $\rightarrow$ pyruvate, feeding lactate into gluconeogenesis
• Has glucose-6-phosphatase for glucose release into the blood

Glycolysis in Evolutionary Perspective

Glycolysis is considered one of the most ancient metabolic pathways. Several lines of evidence support this:

• Universal distribution: Found in virtually all living organisms, from archaea to mammals
• Cytoplasmic localization: Does not require membrane-bound organelles (consistent with pre-mitochondrial origin)
• Anaerobic operation: Functions without oxygen, consistent with the anoxic early Earth atmosphere
• Conserved enzymes: Glycolytic enzymes show high sequence homology across all domains of life

The pathway likely evolved before the Great Oxidation Event (~2.4 billion years ago). The subsequent rise of atmospheric O$_2$ enabled the evolution of oxidative phosphorylation, increasing ATP yield per glucose from 2 to ~32 — a 16-fold improvement that drove the evolution of complex multicellular life.

Python Simulations