Regulation of Glycolysis

A sprinter's quadriceps can multiply its rate of ATP turnover a hundredfold within seconds of the starting gun. Glycolytic flux must rise to match it almost instantly. This is not accomplished by making more enzyme. There is no time to transcribe or translate anything in the first ten seconds of a 100-meter dash. Instead, it is accomplished by throwing molecular switches on enzymes already sitting in the cytosol. These switches are built from the physical logic of allostery, covalent modification, and the peculiar thermodynamics of reactions that run far from equilibrium. Understanding glycolytic regulation means understanding why only three of the ten glycolytic reactions can serve as control points at all. It also means understanding why the cell chose exactly those three. The reason is buried in Actual Free-Energy Change vs Standard Conditions. Most glycolytic reactions operate near equilibrium inside the cell, with ΔG\Delta G close to zero. A reaction near equilibrium cannot be a valve. Nudging its forward rate up nudges its reverse rate up by nearly the same amount, so net flux barely changes. Only a handful of steps are pulled so far from equilibrium that ΔG\Delta G is large and negative under cellular conditions. These steps are effectively irreversible in vivo. They are the only places where an enzyme's catalytic rate directly sets the net flux through the whole pathway. Hexokinase, phosphofructokinase-1 (PFK-1), and pyruvate kinase are precisely these three nonequilibrium steps. Evolution has decorated each with allosteric sites, competing substrates, and, for pyruvate kinase, a hormonally controlled phosphorylation site. The big picture is a hierarchy of control acting on different timescales. PFK-1 is the master flow-control valve, tuned second-to-second by the energy charge of the cell. It also responds to signals of fuel abundance elsewhere in metabolism. Hexokinase is a local gatekeeper, shut off by its own product. This prevents glucose-6-phosphate from piling up when downstream steps stall. Pyruvate kinase is a committed step near the exit. It is poised to be shut down the moment gluconeogenesis needs to run in the same cell. This reciprocal wiring is exactly what Reciprocal Regulation of Glycolysis and Gluconeogenesis builds upon. This article works through the mechanism of control at each of the three enzymes. It covers the allosteric ligands, the kinetic consequences, and the phosphorylation cascade that inactivates liver pyruvate kinase. It then derives, from first principles of near-equilibrium and far-from-equilibrium thermodynamics, exactly why these three steps are the regulated ones, and no others. It closes with worked numerical examples using real cellular concentrations and ΔG\Delta G values. This way, the abstract claim "these steps are irreversible in the cell" becomes a checkable arithmetic fact rather than a slogan.

Key takeaways
  • Glycolysis has ten reactions but only three control points: hexokinase, PFK-1, and pyruvate kinase.
  • These three catalyze reactions with large negative ΔG\Delta G in vivo; the other seven run near equilibrium (ΔG≈0\Delta G \approx 0) and cannot be flux-controlling.
  • PFK-1 is allosterically inhibited by ATP and citrate, activated by AMP, ADP, and fructose-2,6-bisphosphate — the dominant flux valve.
  • Hexokinase is product-inhibited by glucose-6-phosphate; the liver isozyme glucokinase is not, letting the liver keep absorbing glucose after a meal.
  • Pyruvate kinase is allosterically activated by fructose-1,6-bisphosphate (feedforward) and inactivated by glucagon-driven PKA phosphorylation in liver.
  • Key idea: only reactions far from equilibrium can be rate-controlling, and the cell has placed its allosteric and hormonal switches exactly on those three steps.

Where this lives today

Fructose-2,6-bisphosphate signaling, the PFK-1 control axis derived here, is a direct target of metabolic disease research: the bifunctional enzyme PFKFB (6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase) that makes and degrades F2,6BP is a drug target explored for type 2 diabetes, since boosting hepatic F2,6BP could favor glycolysis over gluconeogenesis and lower fasting blood glucose. The pyruvate kinase story has taken on major clinical importance in oncology: many cancers preferentially express the embryonic PKM2 isoform of pyruvate kinase, whose lower intrinsic activity and unusual allosteric regulation (activated by serine, inhibited by tyrosine-phosphorylated growth-factor signaling) is thought to help redirect glycolytic intermediates into biosynthetic pathways, a molecular underpinning of the Warburg effect exploited by growing tumors; small-molecule PKM2 activators are in clinical development as potential anticancer strategies. Hexokinase regulation is directly clinically relevant in a distinct way: pyruvate kinase deficiency and hexokinase deficiency are both causes of hereditary hemolytic anemia in red blood cells, which depend entirely on glycolysis for ATP since they lack mitochondria, and understanding the control architecture explained here is essential to interpreting these enzymopathies. Finally, the entire near-equilibrium/far-from-equilibrium distinction generalizes far beyond glycolysis — it is the same logic used to identify rate-controlling steps in the urea cycle, the citric acid cycle, and industrial metabolic-engineering efforts to redirect flux in engineered microbial strains for biofuel and pharmaceutical production.

Ask the tutor: How does this same near-equilibrium/far-from-equilibrium logic explain why gluconeogenesis needs its own distinct bypass enzymes at exactly these three steps?

The physical picture — thorough prose

Glycolysis converts one glucose into two pyruvate through ten enzyme-catalyzed steps, but not all ten steps are created equal thermodynamically. Consider the mass-action ratio Q=[products]/[substrates]Q = [\text{products}]/[\text{substrates}] for each step measured with actual intracellular concentrations, and compare QQ to the equilibrium constant KeqK_{eq}. When Q≈KeqQ \approx K_{eq}, the reaction is near equilibrium: ΔG=RTln⁡(Q/Keq)≈0\Delta G = RT\ln(Q/K_{eq}) \approx 0, and the enzyme, however fast, is not setting the pace — it is simply fast enough to keep the reaction pinned near equilibrium in both directions. When Q≪KeqQ \ll K_{eq}, the reaction is displaced far from equilibrium, ΔG\Delta G is strongly negative, and the reaction proceeds essentially only in the forward direction. Such a step behaves as a one-way valve, and it is only at a one-way valve that turning the valve (activating or inhibiting the enzyme) changes the net downstream flow. In red blood cell and hepatocyte measurements, seven of the ten glycolytic reactions — phosphoglucose isomerase, aldolase, triose phosphate isomerase, glyceraldehyde-3-phosphate dehydrogenase, phosphoglycerate kinase, phosphoglycerate mutase, and enolase — have actual ΔG\Delta G values within a few kJ/mol of zero. The remaining three — hexokinase, PFK-1, and pyruvate kinase — have actual ΔG\Delta G values around −20-20 to −30 kJ/mol-30\ \text{kJ/mol}, an order of magnitude more negative, because each couples a highly exergonic phosphoryl transfer (from ATP, or to ADP in the pyruvate kinase case) to a step that would otherwise be uphill. These three enzymes are the pathway's committed, essentially irreversible reactions, and control is concentrated there because control anywhere else would be thermodynamically futile. The regulatory logic layered onto these three enzymes reads like an energy-accounting ledger. PFK-1 sees ATP (a "too much fuel" signal) as an inhibitor and AMP/ADP (a "fuel is low, burn faster" signal) as an activator — this is the cell's energy charge sensor. It also sees citrate, a signal that the citric acid cycle is saturated with acetyl-CoA, as an inhibitor, coordinating glycolytic flux with downstream capacity. Fructose-2,6-bisphosphate (F2,6BP), synthesized by a bifunctional kinase/phosphatase under hormonal control, is the single most powerful PFK-1 activator and is the molecular handle by which glucagon and insulin steer the liver between glycolysis and gluconeogenesis. Hexokinase, upstream, is shut off by its immediate product glucose-6-phosphate (G6P) so that if PFK-1 is throttled down, G6P does not accumulate indefinitely and glucose is not needlessly trapped as phosphate ester. Pyruvate kinase, at the exit, is activated in a feedforward loop by F1,6BP (the product of PFK-1 reaching downstream to prime the final committed step) and is switched off in liver by glucagon-stimulated phosphorylation, preventing a futile cycle when the liver needs to run gluconeogenesis instead.

Figure
The three nonequilibrium steps of glycolysis and their allosteric controllers
Glycolysis: near-equilibrium steps vs the three control points Glucose 2 Pyruvate Hexokinase ΔG ≈ −33 kJ/mol ⊖ G6P PGI · Aldolase · TPI · GAPDH · PGK · PGM · Enolase (ΔG ≈ 0, near equilibrium) PFK-1 ΔG ≈ −25 kJ/mol ⊖ ATP, citrate ⊕ AMP, ADP, F2,6BP Pyruvate kinase ΔG ≈ −24 kJ/mol ⊕ F1,6BP ⊖ ATP, alanine; ⊖ PKA-P (liver) Only where ΔG is large and negative can enzyme activity set net pathway flux. Q ≪ K_eq at the three red boxes → one-way valves; Q ≈ K_eq elsewhere → freely reversible, non-controlling

The core ideas — the governing equations/relationships

The central quantitative relationship in this entire topic is the connection between the actual free-energy change of a reaction inside the cell and its mass-action ratio:

ΔG=ΔG∘′+RTln⁡Q\Delta G = \Delta G^{\circ\prime} + RT \ln Q

where:

  • ΔG\Delta G — actual free-energy change under cellular conditions (kJ/mol)
  • ΔG∘′\Delta G^{\circ\prime} — standard free-energy change at pH 7, 1 M reactants/products, 298 K
  • R=8.314 J mol−1K−1R = 8.314\ \text{J mol}^{-1}\text{K}^{-1}, TT — absolute temperature (K)
  • QQ — mass-action ratio, the ratio of actual product to substrate concentrations (raised to stoichiometric powers)

Reading: ΔG\Delta G is not a fixed property of a reaction — it depends on how far the actual cellular concentrations (QQ) sit from the equilibrium ratio (KeqK_{eq}, where ΔG∘′=−RTln⁡Keq\Delta G^{\circ\prime} = -RT\ln K_{eq}). A reaction can have a very unfavorable ΔG∘′\Delta G^{\circ\prime} and still run forward in the cell if QQ is kept far below KeqK_{eq} by rapid removal of products. The condition for a step to be a viable control point can be written as a near-equilibrium criterion:

∣ΔGΔG∘′∣≪1(non-controlling, near equilibrium),∣ΔG∣≫RT(controlling, far from equilibrium)\left|\frac{\Delta G}{\Delta G^{\circ\prime}}\right| \ll 1 \quad \text{(non-controlling, near equilibrium)}, \qquad |\Delta G| \gg RT \quad \text{(controlling, far from equilibrium)}

where:

  • RTRT at 37 °C (310 K) is 8.314×310=2.58 kJ/mol8.314 \times 310 = 2.58\ \text{kJ/mol}, the natural thermal energy scale
  • a ΔG\Delta G of magnitude several RTRT (i.e., tens of kJ/mol) means the reverse reaction is negligible and flux is one-directional

Reading: This is the formal statement of "irreversible in vivo." If ∣ΔG∣|\Delta G| is only comparable to RTRT (a few kJ/mol), both forward and reverse fluxes are substantial and nearly balanced, and modulating the enzyme's total activity barely moves net flux. If ∣ΔG∣|\Delta G| is many multiples of RTRT, essentially all flux is forward, and the enzyme's activity is the net flux. For allosteric regulation itself, the key kinetic relationship is the cooperative (Hill-type) rate law used to describe PFK-1's sigmoidal response to its substrate fructose-6-phosphate (F6P) in the presence of allosteric effectors:

v=Vmax[S]nHK0.5nH+[S]nHv = \frac{V_{max}[S]^{n_H}}{K_{0.5}^{n_H} + [S]^{n_H}}

where:

  • vv — reaction velocity, VmaxV_{max} — maximal velocity
  • [S][S] — substrate concentration (F6P for PFK-1)
  • K0.5K_{0.5} — substrate concentration giving half-maximal velocity (an apparent KmK_m for a cooperative enzyme)
  • nHn_H — Hill coefficient, a measure of cooperativity (nH>1n_H > 1 indicates positive cooperativity)

Reading: Allosteric inhibitors like ATP increase K0.5K_{0.5} (shift the curve right, requiring more F6P for the same velocity) without necessarily changing VmaxV_{max}; activators like AMP or F2,6BP decrease K0.5K_{0.5} (shift left) and can also increase nHn_H sensitivity. This sigmoidal shape means small changes in F6P or effector concentration near K0.5K_{0.5} produce large changes in velocity — an ideal switch, in contrast to the hyperbolic Michaelis–Menten response of a non-regulated near-equilibrium enzyme.

Derivation AWhy only far-from-equilibrium steps can control flux

Goal: show quantitatively that a near-equilibrium enzyme cannot set pathway flux, while a far-from-equilibrium enzyme can, using the relation between net rate and ΔG\Delta G.

Step 1 — net rate as forward minus reverse flux. For any enzyme-catalyzed step A⇌BA \rightleftharpoons B, the net rate is

vnet=vf−vrv_{net} = v_f - v_r

where vfv_f and vrv_r are the (concentration-dependent) forward and reverse unidirectional fluxes catalyzed by the same enzyme population. Basic transition-state theory (or simple mass-action kinetics) gives the ratio of these unidirectional fluxes as

vrvf=QKeq=exp⁡(ΔGRT)\frac{v_r}{v_f} = \frac{Q}{K_{eq}} = \exp\left(\frac{\Delta G}{RT}\right)

This identity follows directly from ΔG=RTln⁡(Q/Keq)\Delta G = RT\ln(Q/K_{eq}) rearranged and exponentiated; it is exact for any elementary or overall reaction obeying mass-action thermodynamics, independent of the enzyme's mechanism.

Step 2 — near-equilibrium regime. Suppose ΔG→0\Delta G \to 0 (a near-equilibrium step, as for phosphoglycerate kinase or triose phosphate isomerase in the cell). Then

vrvf=exp⁡(0)=1\frac{v_r}{v_f} = \exp(0) = 1

so vr≈vfv_r \approx v_f, and the net rate is a small difference of two large, nearly equal unidirectional fluxes:

vnet=vf(1−vrvf)=vf[1−exp⁡ ⁣(ΔGRT)]≈vf(−ΔGRT)(ΔG≪RT)v_{net} = v_f\left(1 - \frac{v_r}{v_f}\right) = v_f\left[1 - \exp\!\left(\frac{\Delta G}{RT}\right)\right] \approx v_f\left(-\frac{\Delta G}{RT}\right) \quad (\Delta G \ll RT)

using the linear expansion ex≈1+xe^x \approx 1+x for small x=ΔG/RTx = \Delta G/RT. Now suppose an allosteric activator doubles the enzyme's total catalytic activity, so both vfv_f and vrv_r double (since the same enzyme catalyzes both directions and QQ, KeqK_{eq} are unchanged by the enzyme). The net rate also just doubles — but QQ will immediately readjust because downstream near-equilibrium reactions rapidly re-equilibrate concentrations, pulling ΔG\Delta G back toward zero and vnetv_{net} back down. In steady state, a near-equilibrium enzyme's total activity is in large excess of what is needed, and its rate constant is not the variable that sets the pathway's steady-state flux; the flux is set elsewhere and this enzyme simply relaxes to whatever QQ the flux demands.

Step 3 — far-from-equilibrium regime. Now suppose ∣ΔG∣≫RT|\Delta G| \gg RT (tens of kJ/mol, several multiples of RT≈2.58 kJ/molRT \approx 2.58\ \text{kJ/mol} at 37 °C), as for PFK-1 in vivo. Then

exp⁡ ⁣(ΔGRT)≈0(ΔG large and negative)\exp\!\left(\frac{\Delta G}{RT}\right) \approx 0 \quad (\Delta G \text{ large and negative})

so vr≈0v_r \approx 0 and

vnet≈vfv_{net} \approx v_f

The net flux through this step is the forward catalytic rate of the enzyme, full stop — there is no comparably-sized reverse flux to cancel against. If an allosteric effector changes the enzyme's forward catalytic capacity (by shifting K0.5K_{0.5} or VmaxV_{max} in the Hill equation), vfv_f changes directly and proportionally, and because vrv_r stays negligible, vnetv_{net} changes by essentially the same factor. Downstream near-equilibrium reactions cannot "erase" this change the way they could in Step 2, because this reaction is not coupled back to its own reverse flux — it is one-way. Hence a change in enzyme activity at a far-from-equilibrium step propagates directly into a change in steady-state pathway flux.

Conclusion: the algebra shows that regulatory sensitivity requires ∣ΔG∣≫RT|\Delta G| \gg RT. This is a thermodynamic, not merely empirical, reason why hexokinase, PFK-1, and pyruvate kinase — and only these three — are the sites of glycolytic control.

Go deeper · Derivation B · 2 minAllosteric shift of PFK-1's apparent K0.5K_{0.5} by ATP and AMP — a Monod–Wyman–Changeux (MWC) sketch
Derivation BAllosteric shift of PFK-1's apparent K0.5K_{0.5} by ATP and AMP — a Monod–Wyman–Changeux (MWC) sketch

Goal: derive how competing allosteric ligands shift the F6P-binding curve of PFK-1, explaining the sigmoidal activation/inhibition pattern quantitatively.

Step 1 — the MWC two-state framework. PFK-1 is a homotetramer existing in equilibrium between a low-affinity "tense" (T) state and a high-affinity "relaxed" (R) state for its substrate F6P:

L=[T0][R0]L = \frac{[T_0]}{[R_0]}

where LL is the allosteric equilibrium constant in the absence of ligand (for PFK-1, LL is large, meaning T dominates and the enzyme is intrinsically sluggish toward F6P unless pushed toward R).

Step 2 — fractional saturation with F6P alone. For nn identical subunits each capable of binding F6P with dissociation constants KRK_R (R state) and KTK_T (T state, KT≫KRK_T \gg K_R), the MWC fractional saturation is

Yˉ=α(1+α)n−1+Lcα(1+cα)n−1(1+α)n+L(1+cα)n\bar{Y} = \frac{\alpha(1+\alpha)^{n-1} + Lc\alpha(1+c\alpha)^{n-1}}{(1+\alpha)^n + L(1+c\alpha)^n}

where α=[F6P]/KR\alpha = [\text{F6P}]/K_R and c=KR/KT≪1c = K_R/K_T \ll 1. This produces the sigmoidal (positively cooperative) F6P-saturation curve characteristic of PFK-1 when LL is large: at low [F6P][\text{F6P}] the enzyme is mostly T (low affinity), and binding of the first F6P molecules progressively pulls the equilibrium toward R (high affinity), steepening the curve — this is the mechanistic origin of the Hill coefficient nH>1n_H > 1 used phenomenologically above.

Step 3 — effect of ATP (allosteric inhibitor) binding preferentially to T. ATP at high concentration binds a distinct allosteric site with much higher affinity for the T conformation than the R conformation. Binding of jj ATP molecules to the T state multiplies the T-state statistical weight by a factor (1+β)n(1+\beta)^n where β=[ATP]/KATP,T\beta = [\text{ATP}]/K_{ATP,T}, effectively increasing the apparent allosteric constant:

Lapp=L (1+β)n(1+βKATP,T/KATP,R)n≈L(1+β)n(KATP,R≫KATP,T)L_{app} = L\,\frac{(1+\beta)^n}{(1+\beta K_{ATP,T}/K_{ATP,R})^n} \approx L(1+\beta)^n \quad (K_{ATP,R} \gg K_{ATP,T})

Because LappL_{app} increases with [ATP][\text{ATP}], the T state is stabilized further, and from Step 2's formula, a larger LL shifts Yˉ\bar Y toward lower values at any fixed [F6P][\text{F6P}] — i.e., it raises the apparent K0.5K_{0.5} for F6P. This is exactly the observed effect: ATP increases the concentration of F6P needed for half-maximal PFK-1 velocity, without changing VmaxV_{max} (since ATP is also the phosphoryl donor substrate at a separate catalytic site, saturating that site).

Step 4 — effect of AMP/ADP and F2,6BP (activators) binding preferentially to R. These ligands bind an allosteric site with higher affinity for R than T, which — by the same algebra with the roles of R and T exchanged — decreases LappL_{app}:

Lapp=L(1+γ)n(γ=[AMP]KAMP,R, KAMP,T≫KAMP,R)L_{app} = \frac{L}{(1+\gamma)^n} \quad \left(\gamma = \frac{[\text{AMP}]}{K_{AMP,R}},\ K_{AMP,T} \gg K_{AMP,R}\right)

A smaller LappL_{app} means less T character intrinsically, so the F6P-saturation curve shifts left (lower apparent K0.5K_{0.5}) and becomes less sigmoidal (approaching hyperbolic as Lapp→1L_{app} \to 1). AMP and F2,6BP therefore both sensitize PFK-1 to physiological F6P concentrations and antagonize ATP inhibition, because ATP and AMP/F2,6BP pull LappL_{app} in opposite directions — this is the molecular basis of the observation that AMP relieves ATP inhibition of PFK-1 even though ATP itself is also PFK-1's phosphoryl-donor substrate.

Conclusion: the MWC algebra shows explicitly how two allosteric ligands acting on the same conformational equilibrium, from the same T/R switch, produce opposite shifts in the apparent K0.5K_{0.5} for the true substrate — turning the tetramer into a graded, reversible, energy-charge–sensing valve.

Worked examples

Example 1 (foundational) — Is the hexokinase step really irreversible in erythrocytes? The hexokinase reaction is

Glucose+ATP→Glucose-6-phosphate+ADP\text{Glucose} + \text{ATP} \rightarrow \text{Glucose-6-phosphate} + \text{ADP}

with ΔG∘′=−16.7 kJ/mol\Delta G^{\circ\prime} = -16.7\ \text{kJ/mol}. Typical human erythrocyte concentrations (in mM) are: glucose =5.0= 5.0, ATP =1.85= 1.85, G6P =0.083= 0.083, ADP =0.14= 0.14. Take T=310 KT = 310\ \text{K}, R=8.314 J mol−1K−1R = 8.314\ \text{J mol}^{-1}\text{K}^{-1}. Compute the mass-action ratio:

Q=[G6P][ADP][Glucose][ATP]=(0.083)(0.14)(5.0)(1.85)Q = \frac{[\text{G6P}][\text{ADP}]}{[\text{Glucose}][\text{ATP}]} = \frac{(0.083)(0.14)}{(5.0)(1.85)}
Q=0.011629.25=1.256×10−3Q = \frac{0.01162}{9.25} = 1.256\times10^{-3}

Then

ΔG=ΔG∘′+RTln⁡Q=−16,700+(8.314)(310)ln⁡(1.256×10−3) J/mol\Delta G = \Delta G^{\circ\prime} + RT\ln Q = -16{,}700 + (8.314)(310)\ln(1.256\times10^{-3})\ \text{J/mol}
ln⁡(1.256×10−3)=−6.680\ln(1.256\times10^{-3}) = -6.680
ΔG=−16,700+(2577)(−6.680)=−16,700−17,214=−33,914 J/mol≈−33.9 kJ/mol\Delta G = -16{,}700 + (2577)(-6.680) = -16{,}700 - 17{,}214 = -33{,}914\ \text{J/mol} \approx -33.9\ \text{kJ/mol}

Check against the near-equilibrium criterion: RT=2.58 kJ/molRT = 2.58\ \text{kJ/mol} at 37 °C, so ∣ΔG∣/RT≈33.9/2.58≈13.1|\Delta G|/RT \approx 33.9/2.58 \approx 13.1 — the actual free-energy change is over thirteen times the thermal energy scale, meaning exp⁡(ΔG/RT)≈e−13.1≈2×10−6\exp(\Delta G/RT) \approx e^{-13.1} \approx 2\times10^{-6}: the reverse flux is negligible by six orders of magnitude relative to the forward flux. Hexokinase is decisively a one-way valve in vivo, fully consistent with it being a control point, and the number (−33.9 kJ/mol-33.9\ \text{kJ/mol}) is more negative than ΔG∘′\Delta G^{\circ\prime} itself because cellular ADP and G6P are both kept well below their equilibrium levels by downstream consumption and by adenylate kinase/other ATP-regenerating reactions.

Limiting case: if downstream consumption of G6P were to stop entirely and G6P accumulated until Q→KeqQ \to K_{eq}, then ΔG→0\Delta G \to 0 and hexokinase would cease to be a valve — this is exactly the physiological function of G6P product inhibition, which throttles hexokinase before QQ ever gets close to that limit.

Example 2 (applied) — PFK-1 sensitivity to the ATP/AMP ratio during muscle contraction At rest, human skeletal muscle typically has approximately ATP =4.0 mM= 4.0\ \text{mM}, ADP =0.013 mM= 0.013\ \text{mM}, and (via the adenylate kinase equilibrium 2ADP⇌ATP+AMP2\text{ADP} \rightleftharpoons \text{ATP} + \text{AMP}) AMP is present at only about 0.0010.001–0.002 mM0.002\ \text{mM} because it is the smallest pool and any ADP produced is rapidly interconverted. The adenylate kinase equilibrium constant is close to 1:

KAK=[ATP][AMP][ADP]2≈1K_{AK} = \frac{[\text{ATP}][\text{AMP}]}{[\text{ADP}]^2} \approx 1

so at rest, with [ADP]=0.013 mM[\text{ADP}] = 0.013\ \text{mM} and [ATP]=4.0 mM[\text{ATP}] = 4.0\ \text{mM}:

[AMP]=[ADP]2[ATP]=(0.013)24.0=1.69×10−44.0=4.2×10−5 mM[\text{AMP}] = \frac{[\text{ADP}]^2}{[\text{ATP}]} = \frac{(0.013)^2}{4.0} = \frac{1.69\times10^{-4}}{4.0} = 4.2\times10^{-5}\ \text{mM}

Now suppose vigorous contraction consumes ATP rapidly enough that ADP rises modestly, only threefold, to 0.039 mM0.039\ \text{mM} (a change that would be almost undetectable if you were only watching ATP itself, since ATP might fall by well under 1% from 4.0 to, say, 3.96 mM). Recompute AMP at the new ADP, holding ATP at 3.96 mM:

[AMP]=(0.039)23.96=1.521×10−33.96=3.84×10−4 mM[\text{AMP}] = \frac{(0.039)^2}{3.96} = \frac{1.521\times10^{-3}}{3.96} = 3.84\times10^{-4}\ \text{mM}

The fold-change in AMP is

3.84×10−44.2×10−5≈9.1\frac{3.84\times10^{-4}}{4.2\times10^{-5}} \approx 9.1

A threefold rise in ADP produces roughly a ninefold rise in AMP, because the adenylate kinase equilibrium makes [AMP][\text{AMP}] scale as [ADP]2/[ATP][\text{ADP}]^2/[\text{ATP}] — the square-law amplification is the physical reason AMP, not ADP or ATP, is the most sensitive internal signal of energy stress, and why PFK-1 (and AMP-activated protein kinase in other contexts) is tuned to respond to AMP specifically. Given PFK-1's Hill-type dependence on activators (Step 4 of the second derivation), this near-order-of-magnitude AMP surge is enough to substantially lower LappL_{app}, shift K0.5K_{0.5} for F6P down, and sharply increase PFK-1 velocity at the essentially unchanged physiological F6P concentration — exactly the flux surge required at the onset of intense exercise.

Limiting case: if ATP were held perfectly constant (no consumption at all), ADP would not rise, AMP would not rise by the square-law amplification, and PFK-1 would receive no activating signal — consistent with the fact that a muscle at rest, with ATP turnover matched by ATP synthesis, shows no PFK-1 activation despite the presence of some basal AMP.

Example 3 (exam/competition level) — Reconciling pyruvate kinase's ΔG∘′\Delta G^{\circ\prime} with observed flux under a phosphorylated (inactivated) state Pyruvate kinase catalyzes

Phosphoenolpyruvate (PEP)+ADP→Pyruvate+ATP,ΔG∘′=−31.4 kJ/mol\text{Phosphoenolpyruvate (PEP)} + \text{ADP} \rightarrow \text{Pyruvate} + \text{ATP}, \qquad \Delta G^{\circ\prime} = -31.4\ \text{kJ/mol}

Suppose in resting hepatocyte cytosol the concentrations are PEP =0.023 mM= 0.023\ \text{mM}, ADP =0.14 mM= 0.14\ \text{mM}, pyruvate =0.051 mM= 0.051\ \text{mM}, ATP =2.5 mM= 2.5\ \text{mM}, at T=310 KT = 310\ \text{K}.

Q=[Pyruvate][ATP][PEP][ADP]=(0.051)(2.5)(0.023)(0.14)=0.12750.00322=39.6Q = \frac{[\text{Pyruvate}][\text{ATP}]}{[\text{PEP}][\text{ADP}]} = \frac{(0.051)(2.5)}{(0.023)(0.14)} = \frac{0.1275}{0.00322} = 39.6
ΔG=−31,400+(8.314)(310)ln⁡(39.6)=−31,400+(2577)(3.679)\Delta G = -31{,}400 + (8.314)(310)\ln(39.6) = -31{,}400 + (2577)(3.679)
ΔG=−31,400+9,481=−21,919 J/mol≈−21.9 kJ/mol\Delta G = -31{,}400 + 9{,}481 = -21{,}919\ \text{J/mol} \approx -21.9\ \text{kJ/mol}

This is still strongly negative (∣ΔG∣/RT≈8.5|\Delta G|/RT \approx 8.5), confirming pyruvate kinase remains a genuine one-way, far-from-equilibrium step even under fasting conditions where flux through it is minimized. Now, glucagon signaling in the fasted liver activates protein kinase A (PKA), which phosphorylates liver pyruvate kinase (L-isozyme) at a regulatory serine. This phosphorylation does not change ΔG∘′\Delta G^{\circ\prime} (a thermodynamic property of the reaction, independent of the enzyme) — it lowers the enzyme's catalytic efficiency toward PEP (raises apparent KmK_m for PEP, and increases sensitivity to allosteric inhibition by ATP and alanine), which lowers vfv_f directly. By the far-from-equilibrium result of the first derivation, vnet≈vfv_{net} \approx v_f here, so a drop in vfv_f from phosphorylation translates almost one-to-one into a drop in net flux through pyruvate kinase — exactly the intended physiological outcome: shutting the glycolytic exit valve so that PEP accumulating from gluconeogenic precursors is not immediately re-consumed by glycolysis, avoiding a futile cycle of PEP →\to pyruvate →\to (reverse gluconeogenic steps) →\to PEP that would burn ATP and GTP for no net synthesis. Notice that the phosphorylated enzyme's QQ would, if flux truly stopped, rise as PEP accumulates upstream — but because pyruvate kinase remains catalytically capable (just less so), a nonzero but sharply reduced vnetv_{net} persists, consistent with observed measurable-but-low hepatic pyruvate kinase flux during fasting.

Limiting case: if [PEP][\text{PEP}] rose enough (as it does during active gluconeogenesis, where PEP is a gluconeogenic intermediate being synthesized) that QQ actually approached KeqK_{eq} for this reaction, ΔG\Delta G would approach zero and pyruvate kinase would cease to be a strict one-way valve at all — this is not observed physiologically because PKA phosphorylation intervenes first, keeping catalytic capacity low enough that PEP never needs to build up that far before flux is throttled.

Figure
MWC allosteric shift of PFK-1: ATP vs AMP effect on the F6P saturation curve
PFK-1 velocity vs [F6P]: allosteric effector shifts [F6P] (mM) v / V_max +AMP, +F2,6BP (low K₀.₅) basal +ATP, +citrate (high K₀.₅) physiological [F6P]

Limiting cases, scope & common traps

SituationWhat happensWhy
ΔG→0\Delta G \to 0 at a step (e.g., TPI, PGK in vivo)Enzyme cannot control net flux even if very abundant/fastvr≈vfv_r \approx v_f; net rate is a small residual, set by downstream demand pulling QQ
∥ΔG∥≫RT\|\Delta G\| \gg RT (hexokinase, PFK-1, pyruvate kinase)Enzyme activity directly sets net fluxvr≈0v_r \approx 0; vnet≈vfv_{net} \approx v_f, no cancellation
Liver at high blood glucoseGlucokinase (hexokinase D/IV) not inhibited by G6P, high KmK_mLiver keeps phosphorylating glucose even as G6P accumulates, buffering blood glucose
Fasting liver, glucagon highPyruvate kinase phosphorylated (inactive), PFK-1 low F2,6BP (inactive)Blocks glycolysis so gluconeogenic flux is not futile-cycled away
Fed liver, insulin highPyruvate kinase dephosphorylated (active), PFK-1 high F2,6BP (active)Glycolysis flux favored, glucose committed to catabolism/storage
Exercising muscleAMP surges nonlinearly (square-law via adenylate kinase)PFK-1 activated sharply even with small ATP drop
Cancer cells (Warburg effect)PFK-1 and pyruvate kinase (PKM2 isoform) regulation alteredSustained high glycolytic flux even with oxygen present

The chemist behind it

The modern framework for glycolytic control rests heavily on the work of Earl Sutherland, who discovered cyclic AMP and hormone-triggered phosphorylation cascades — the molecular chain that lets glucagon, via PKA, phosphorylate and inactivate liver pyruvate kinase — and on Edwin Krebs and Edmond Fischer, who established reversible protein phosphorylation as a general regulatory mechanism, work for which they shared the 1992 Nobel Prize in Physiology or Medicine. Equally central is Jacques Monod, whose collaboration with Jeffries Wyman and Jean-Pierre Changeux produced the concerted (MWC) model of allosteric transitions used above to derive PFK-1's sigmoidal kinetics; Monod's parallel insight — that regulatory and catalytic functions of a protein can be physically separated onto distinct sites without direct steric overlap — is the conceptual key that makes it possible for ATP to be simultaneously PFK-1's substrate at one site and its inhibitor at another.

Origins · 1 minJacques Monod and how this was discovered
JM
Jacques Monod
1910–1976 · Institut Pasteur, Paris

Monod's studies of bacterial enzyme induction and allosteric feedback inhibition, developed alongside François Jacob and formalized mathematically with Jeffries Wyman and Jean-Pierre Changeux in the 1965 MWC model, gave biochemistry its first rigorous, testable theory of how a protein's function can be regulated by ligands binding away from the catalytic site. The 1965 paper "On the Nature of Allosteric Transitions: A Plausible Model" derived the two-state (T/R) partition-function treatment used directly in modeling PFK-1's response to ATP and AMP. Monod shared the 1965 Nobel Prize in Physiology or Medicine with Jacob and Lwoff for work on genetic regulation of enzyme synthesis, and his allosteric framework became the universal language for describing PFK-1, hemoglobin, aspartate transcarbamoylase, and virtually every cooperative regulatory enzyme discovered since.

History
- 1943 — Otto Meyerhof and colleagues complete the elucidation of the glycolytic pathway's ten enzymatic steps.
- 1956 — Earl Sutherland discovers cyclic AMP, establishing the second-messenger basis for glucagon/epinephrine signaling to glycolytic enzymes.
- 1965 — Monod, Wyman, and Changeux publish the concerted (MWC) allosteric model, later applied directly to PFK-1 kinetics.
- 1968–1980 — Edwin Krebs, Edmond Fischer, and colleagues establish reversible phosphorylation as the mechanism inactivating liver pyruvate kinase.
- 1980 — Simon Pilkis, Henri-Géry Hers and colleagues identify fructose-2,6-bisphosphate as the potent PFK-1 activator and key hormonal control signal.
Why are there exactly three control points in glycolysis and not more or fewer?

The number is dictated by thermodynamics, not convention. Of the ten glycolytic reactions, only three — hexokinase, PFK-1, and pyruvate kinase — have actual intracellular ΔG\Delta G values large and negative enough (tens of kJ/mol, many multiples of RTRT) that the reverse reaction is negligible. As derived from vr/vf=exp⁡(ΔG/RT)v_r/v_f = \exp(\Delta G/RT), only at such steps does a change in enzyme catalytic activity translate directly into a change in net pathway flux; at the other seven near-equilibrium steps, forward and reverse fluxes are comparable and nearly cancel, so modulating the enzyme barely moves the net rate. Evolution places regulatory machinery (allosteric sites, phosphorylation sites) where it can actually be effective, which is exactly at these three thermodynamically irreversible reactions.

How can ATP be both a substrate and an inhibitor of PFK-1 at the same time?

PFK-1 has two distinct binding sites for ATP: a high-affinity catalytic site where ATP is the phosphoryl donor for converting fructose-6-phosphate to fructose-1,6-bisphosphate, and a separate, lower-affinity allosteric site that, when occupied at high cellular ATP concentrations, stabilizes the low-affinity T conformation of the enzyme. Because the allosteric site only becomes significantly occupied when ATP is in excess of what is needed to saturate the catalytic site, PFK-1 is fully functional at normal-to-low ATP levels (catalytic site engaged, allosteric site mostly empty) but becomes progressively throttled as ATP rises further, exactly the "fuel is abundant, slow down" signal the cell needs. This dual-site architecture, formalized by the Monod–Wyman–Changeux model, is what allows one molecule to play opposite catalytic and regulatory roles simultaneously.

What is fructose-2,6-bisphosphate and why is it so important?

Fructose-2,6-bisphosphate (F2,6BP) is not a glycolytic intermediate itself — it never appears on the main glycolytic path from glucose to pyruvate — but is a dedicated regulatory molecule synthesized from fructose-6-phosphate by the kinase activity of the bifunctional enzyme PFK-2/FBPase-2, and degraded by the phosphatase activity of the same bifunctional protein. Its concentration is set by the phosphorylation state of PFK-2/FBPase-2, which is itself controlled by glucagon (via PKA phosphorylation, favoring the phosphatase, lowering F2,6BP, dampening PFK-1) and insulin (favoring dephosphorylation, favoring the kinase, raising F2,6BP, activating PFK-1). Because F2,6BP is PFK-1's most potent allosteric activator and also inhibits the opposing gluconeogenic enzyme fructose-1,6-bisphosphatase, it acts as a single molecular switch that reciprocally controls glycolysis and gluconeogenesis in the liver.

Does regulation at hexokinase, PFK-1, and pyruvate kinase matter in tissues without hormonal control, like exercising muscle?

Yes, though the balance of mechanisms differs. Skeletal muscle glycolysis is controlled predominantly by local allosteric signals — rising AMP and ADP (via the adenylate kinase equilibrium) and falling ATP/citrate during contraction directly activate PFK-1 within milliseconds to seconds, without needing any hormone. Muscle pyruvate kinase is largely the M1 isozyme, which (unlike the liver L isozyme) is not a PKA substrate and is not switched off by glucagon, consistent with muscle's need to sustain glycolytic flux during exercise regardless of hormonal state. Hexokinase in muscle is still subject to G6P product inhibition, preventing glucose-6-phosphate buildup if downstream flux transiently outpaces glucose uptake. So the same three control points operate everywhere, but the dominant regulatory inputs (allosteric energy-charge signals versus hormonal phosphorylation) differ by tissue and isozyme.

How is the near-equilibrium/far-from-equilibrium classification actually measured experimentally?

Researchers measure the intracellular concentrations of substrates and products for each glycolytic step, typically using rapid freeze-quenching of tissue followed by enzymatic or mass-spectrometric assay of metabolite levels, and independently measure (or take from the literature) each reaction's standard free-energy change ΔG∘′\Delta G^{\circ\prime}. Plugging measured concentrations into the mass-action ratio QQ and applying ΔG=ΔG∘′+RTln⁡Q\Delta G = \Delta G^{\circ\prime} + RT\ln Q gives the actual in vivo ΔG\Delta G for each step. This kind of metabolite-profiling analysis, first done systematically for erythrocytes and liver in the mid-20th century, is what originally identified hexokinase, PFK-1, and pyruvate kinase as the outliers with large negative ΔG\Delta G, and it remains the standard method (now often done with LC-MS/MS metabolomics) for identifying rate-controlling steps in any pathway.

Why does citrate inhibit PFK-1, given that citrate is a citric acid cycle intermediate, not a glycolytic one

Citrate is exported from mitochondria when the citric acid cycle is running at capacity relative to the demand for oxidative phosphorylation — a signal that acetyl-CoA and ATP are already abundant. In the cytosol, elevated citrate binds an allosteric site on PFK-1 (distinct from the ATP site but acting through the same T-state-stabilizing mechanism) and synergizes with ATP inhibition, raising the apparent K0.5K_{0.5} for F6P even further. This is a feedback loop that spans two pathways: it prevents glycolysis from continuing to feed pyruvate (and hence acetyl-CoA) into a citric acid cycle that is already saturated, avoiding a wasteful pileup of intermediates and coordinating cytosolic carbon flux with mitochondrial oxidative capacity.

Is the "three control points" picture the complete story, or are there other layers of regulation?

It captures the core enzymatic logic, but real cells layer several additional mechanisms on top. Substrate channeling and enzyme compartmentalization can modulate local effective concentrations of intermediates; transcriptional control (for instance, insulin- and glucose-driven induction of glucokinase, PFK-1, and pyruvate kinase gene expression in liver) sets the total enzyme amount available over hours to days, a slower layer beneath the fast allosteric/covalent switches discussed here; and cross-talk with signaling pathways (AMPK, PI3K/Akt) further modulates glycolytic enzyme activity in many cell types, including cancer cells exploiting the Warburg effect. The three-control-point framework nonetheless remains the essential first-order picture for understanding how flux is set on a second-to-minute timescale.

Before this

Where this leads