Free Energy and the Direction of Biochemical Reactions

Every claim a biochemist makes about whether a pathway "can" run in a given direction — glycolysis versus gluconeogenesis sharing most enzymes but running opposite net directions under different hormonal and concentration conditions, the citric acid cycle's near-irreversible steps acting as metabolic ratchets, or why the cell needs seven distinct enzymatic steps rather than one to reverse glycolysis — is a direct application of comparing ΔG (not ΔG°′) across physiological concentration ranges. Clinical metabolic disease panels, when they report elevated lactate or altered ATP/ADP ratios, are implicitly reporting a shift in cellular Q away from the values that keep central pathways thermodynamically favorable; profound hypoxia, for instance, collapses the ATP/(ADP·Pi) ratio toward equilibrium and correspondingly collapses the free energy available for every ATP-coupled process in the cell, from ion pumping to muscle contraction. Pharmaceutical enzyme inhibitor design routinely targets the single most thermodynamically committed (large negative ΔG, functionally irreversible) step of a pathway — such as phosphofructokinase-1 in glycolysis — because that is the step whose forward flux is hardest for the cell to reverse by mass action alone, making it the most reliable control point. Bioengineered metabolic pathways, in synthetic biology, are explicitly redesigned around free-energy profiles: pathway engineers compute the ΔG of every proposed step at anticipated intracellular concentrations before ever cloning a gene, specifically to avoid designing a pathway with a thermodynamically "uphill" bottleneck that no enzyme, however good a catalyst, could ever push forward. The same logic underlies why NAD⁺/NADH and ATP/ADP ratios are measured as key biomarkers in cancer metabolism and neurodegeneration research — they are direct proxies for the Q that determines whether central bioenergetic reactions still carry a favorable ΔG in a diseased cell.

Ask the tutor: How does the biochemical standard state (ΔG°′) get defined precisely, and why does fixing pH at 7 change the numeric value of ΔG° for proton-releasing reactions?

Think of a reaction A ⇌ B as a ball on a landscape whose height is the Gibbs free energy of the system, and whose horizontal position is the extent of reaction (how far the mixture has progressed from pure A toward pure B). The ball rolls downhill — that is the Second Law in miniature — and it stops rolling when it reaches the bottom of the local valley, which is the equilibrium composition, the point where the free energy of the mixture is at a minimum. ΔG is the slope of this landscape at the current composition, not the depth of the valley. A steeply negative slope (large negative ΔG) means the reaction is far from equilibrium and strongly favored to proceed forward; a slope of zero (ΔG = 0) means you are sitting exactly at the equilibrium composition, and there is no net driving force in either direction. Crucially, the shape of this landscape — where the valley floor sits — is fixed by the chemistry (bond energies, standard free energies of formation), captured in ΔG°′. But the position of the ball on that landscape is fixed by the actual concentrations of A and B present right now, captured in the reaction quotient Q = [B]/[A]. A cell can hold the ball artificially high on the uphill side of the valley by continuously removing B (product) as fast as it forms and continuously resupplying A (substrate) — that is exactly what a metabolic pathway does when the next enzyme immediately consumes the product of the previous one. This is why an individual step with an unfavorable ΔG°′ can still have a favorable (negative) ΔG in vivo: the cell keeps Q artificially small compared to K′eq by mass action, coupling, and enzymatic pull-through. The second diagram below makes the ΔG–ΔG°′–Q relationship quantitative and shows how the same reaction can sit on either side of ΔG = 0 depending on where the concentrations put it, while the third figure traces the ATP hydrolysis free-energy landscape which recurs constantly as the "downhill" partner in coupled reactions.

:::svg[ΔG of ATP hydrolysis as a function of the actual [ATP]/([ADP][Pi]) ratio in the cell]

The general relation between the actual free-energy change of a reaction and its standard value is:

where:

• $\Delta G$ — the actual Gibbs free-energy change under the concentrations present (kJ/mol), the quantity that determines spontaneity right now
• $\Delta G^{\circ\prime}$ — the biochemical standard free-energy change (1 M reactants/products, pH 7, 298 K, 1 atm)
• $R$ — gas constant, 8.315 J mol⁻¹ K⁻¹
• $T$ — absolute temperature (K)
• $Q$ — the reaction quotient, the ratio of actual product to substrate concentrations (each raised to its stoichiometric coefficient), evaluated at the current (non-equilibrium) concentrations

Reading: the sign of ΔG, not ΔG°′, tells you which way the reaction runs. If Q is much smaller than its equilibrium value (products scarce, substrates abundant), RT ln Q is a large negative number and can push ΔG negative even when ΔG°′ is positive. At equilibrium, by definition, ΔG = 0 and Q takes the special value K′eq, giving the bridge equation:

where:

• $K'_{eq}$ — the equilibrium constant measured under biochemical standard-state conventions (pH 7 fixed, water and H⁺ activities folded in)

Reading: this equation converts a measurable equilibrium constant (from concentrations at equilibrium) into a standard free-energy change, and vice versa — it is how the ΔG°′ values in textbook tables are actually obtained. Combining the two, the working formula used constantly in practice for a reaction A + B ⇌ C + D is:

where each bracketed quantity is the actual cellular (or in vitro) concentration at the moment of interest, not the equilibrium concentration.

Reading: plug in real numbers for what's actually in the cytosol, and you get the real thermodynamic driving force — this is the calculation that separates "textbook favorable" from "physiologically happening."

Example 1 — Foundational: computing ΔG from ΔG°′ and non-standard concentrations. Consider the isomerization catalyzed by triose phosphate isomerase:

with $\Delta G^{\circ\prime} = +7.5$ kJ/mol (unfavorable at standard state — DHAP is thermodynamically preferred at equilibrium). Suppose in a red blood cell the steady-state concentrations are [DHAP] = 0.030 mM and [G3P] = 0.0056 mM (asymmetric, non-round, realistic cellular values — G3P is kept low because it's rapidly consumed by the next glycolytic enzyme). Take $T = 310$ K. The reaction quotient (product over substrate, one-to-one stoichiometry):

Interpretation: even with G3P held well below DHAP, ΔG is still slightly positive here — this single step, in isolation, is only weakly reversible and sits close to equilibrium, consistent with triose phosphate isomerase being one of the most catalytically perfect (diffusion-limited) enzymes known: it must equilibrate its substrates rapidly in both directions because the thermodynamics gives it very little slack. The reaction is pulled forward overall in glycolysis only because G3P is continuously drained by the next irreversible step (glyceraldehyde-3-phosphate dehydrogenase), which keeps [G3P] low enough, on average, for the net flux through the pathway to run forward.

Limiting-case check: if [G3P] = [DHAP] exactly (Q = 1), then ln Q = 0 and ΔG = ΔG°′ = +7.5 kJ/mol exactly — recovering the standard-state definition, as required, since Q = 1 corresponds precisely to the 1 M (or here, equal-concentration) standard-state condition.

Example 2 — Applied: ATP hydrolysis under real cytosolic conditions. The textbook value for ATP hydrolysis,

is $\Delta G^{\circ\prime} = -30.5$ kJ/mol at 25 °C, pH 7. But cells never operate at 1 M ATP, 1 M ADP, 1 M Pi — actual hepatocyte cytosolic concentrations are approximately [ATP] = 3.4 mM, [ADP] = 1.3 mM, [Pi] = 4.8 mM (asymmetric, physiologically realistic, non-round values). Take T = 310 K (37 °C, so we also adjust slightly for temperature, though the dominant effect here is concentration).

Interpretation: the actual free energy released by ATP hydrolysis in a living hepatocyte (about −47 kJ/mol) is substantially more negative than the standard value (−30.5 kJ/mol), because the cell maintains a high [ATP]/([ADP][Pi]) ratio — this is exactly the mechanism that makes ATP such a reliable energy currency: cells actively work (via oxidative phosphorylation and substrate-level phosphorylation) to keep Q far below K′eq, which amplifies the usable free energy per hydrolysis event well beyond the textbook number. This is the number that should be used whenever asking "how much work can this cell get from one ATP," not −30.5 kJ/mol.

Limiting-case check: as the cell's phosphorylation potential collapses toward equilibrium (e.g., in ischemia or severe hypoxia, where ATP is depleted and ADP/Pi accumulate), Q rises toward K′eq, RT ln Q → 0, and ΔG → ΔG°′ = −30.5 kJ/mol — the driving force shrinks toward the standard value and, in the true equilibrium limit, to zero, which is exactly the pathological state of a cell that can no longer do useful work.

Example 3 — Exam/competition level: coupling an unfavorable reaction to ATP hydrolysis and finding the overall equilibrium constant. Consider glutamine synthesis from glutamate and ammonia:

This is unfavorable at standard state. The cell couples it (via the enzyme glutamine synthetase, through a phospho-glutamate intermediate) to ATP hydrolysis:

Because free energies are additive for sequential reactions sharing intermediates (a result you will see formalized separately for coupled reactions), the overall reaction:

has

Now find the overall equilibrium constant at T = 298 K, using the bridge equation:

Interpretation: coupling to ATP hydrolysis flips the equilibrium constant from strongly unfavorable (for glutamine synthesis alone, $K'_{eq,1} = \exp(-14200/2477.9) = \exp(-5.73) \approx 3.2\times10^{-3}$, meaning at equilibrium glutamate would outnumber glutamine roughly 300 to 1) to strongly favorable (K′eq ≈ 720 for the coupled reaction, glutamine now favored roughly 700 to 1 at equilibrium). This is the quantitative essence of "energetic coupling": the two reactions don't chemically merge into one mechanism at the level of orbitals in this case (they proceed through a real phosphorylated intermediate, γ-glutamyl phosphate, so mechanistically they do share an intermediate) but thermodynamically their ΔG°′ values simply add, and the equilibrium constants multiply — a direct algebraic consequence of ΔG°′ = −RT ln K′eq applied to each step, since adding exponents multiplies the exponentials:

(matching the direct calculation above up to rounding, since $K'_{eq,2} = \exp(30500/2477.9) = \exp(12.31) \approx 2.2\times10^{5}$) — confirming the additivity of ΔG°′ is equivalent to multiplicativity of K′eq, a useful cross-check on any coupled-reaction calculation.

Josiah Willard Gibbs (1839–1903) spent almost his entire working life at Yale, where he held an unpaid professorship for nearly a decade before the university finally offered him a salary — a reflection of how little contemporary American science establishment grasped the significance of what he was doing. Working largely alone and publishing in the relatively obscure Transactions of the Connecticut Academy of Arts and Sciences, Gibbs built the entire modern edifice of chemical thermodynamics between 1873 and 1878 in a single extraordinary memoir, "On the Equilibrium of Heterogeneous Substances." In it he introduced the free energy function that now bears his name, defined chemical potential rigorously for the first time, and derived the phase rule — essentially inventing the mathematical machinery by which every subsequent chemist, biochemist, and chemical engineer would decide whether a process could occur spontaneously. His work was so densely mathematical and so far ahead of the experimental chemistry of his era that it was barely read in the United States; it took European physical chemists — Ostwald, van der Waals, and later Duhem — translating and championing it before its full importance was recognized. The application of Gibbs's free energy function specifically to metabolic reactions, and the biochemical standard-state convention (fixing pH at 7, folding water and H⁺ activity into ΔG°′) that makes it usable for real cellular chemistry, is a twentieth-century biochemical development built directly on Gibbs's nineteenth-century thermodynamic foundation — a striking case of pure mathematical physics becoming, decades later, the operating language of cell biology.