Free Energy, Spontaneity, and Coupled Reactions

A yeast cell manufactures glucose-6-phosphate by spending ATP. This happens even though the phosphorylation reaction alone is thermodynamically uphill. Left to itself, glucose plus inorganic phosphate will never assemble into glucose-6-phosphate in appreciable yield. Yet inside the cell this reaction runs essentially to completion. This happens every second, in every dividing cell on Earth. The trick is not a violation of thermodynamics. It is a masterclass in how thermodynamics actually works when reactions share a common intermediate. This article is about the arithmetic of life's energy economy. It covers how Gibbs free energy determines whether a reaction proceeds. It explains how the standard free energy change relates to the actual free energy change under cellular conditions. It also shows how cells stitch an energetically unfavorable reaction to a favorable one so that the sum runs forward. This is the conceptual machine underneath The Chemical Logic of the Cell. It explains why ATP hydrolysis can "pay for" biosynthesis, ion pumping, and mechanical work. It is also why a pathway's overall thermodynamics can never be judged from a single step in isolation. The big picture is deceptively simple once the algebra is unpacked. ΔG\Delta G is not a fixed property of a reaction. It is a function of concentration. ΔG∘′\Delta G^{\circ\prime} (the biochemist's standard state) is merely the value of ΔG\Delta G at one particular, often unphysiological, concentration point. Because ΔG\Delta G depends logarithmically on the reaction quotient QQ, a reaction with an unfavorable ΔG∘′\Delta G^{\circ\prime} can be pulled forward. This happens simply by keeping product concentration low or substrate concentration high. Coupling two reactions through a shared intermediate is the biological embodiment of this fact. It is thermodynamic bookkeeping, not new chemistry. By the end you will be able to compute ΔG\Delta G from ΔG∘′\Delta G^{\circ\prime} and real concentrations. You will also convert between ΔG∘′\Delta G^{\circ\prime} and Keq′K'_{eq}. You will correctly add coupled-reaction free energies. And you will diagnose why a pathway step that looks impossible in isolation is entirely routine in vivo. This is Thermodynamics of Living Systems made quantitative and operational.

Key takeaways
  • ΔG∘′\Delta G^{\circ\prime} is the free energy change at defined standard conditions (1 M, pH 7, 298 K); it is a constant for a given reaction, but it is not what actually happens in the cell.
  • ΔG\Delta G is the free energy change at actual cellular concentrations; it, not ΔG∘′\Delta G^{\circ\prime}, determines the direction and spontaneity of a reaction as it runs.
  • ΔG∘′\Delta G^{\circ\prime} and Keq′K'_{eq} are two faces of one coin: ΔG∘′=−RTln⁡Keq′\Delta G^{\circ\prime} = -RT\ln K'_{eq}.
  • Free energies of consecutive reactions sharing an intermediate are additive; this lets an unfavorable step be driven by a favorable one.
  • ATP hydrolysis (ΔG∘′≈−30.5 kJ/mol\Delta G^{\circ\prime} \approx -30.5\ \text{kJ/mol}) is the universal energy currency because it couples to otherwise endergonic biosynthetic and transport steps.
  • Cells maintain far-from-equilibrium concentrations (high [ATP]/[ADP][Pi]) precisely so that ΔG\Delta G for ATP hydrolysis is much more negative than ΔG∘′\Delta G^{\circ\prime}.
  • Key idea: spontaneity in the cell is governed by the actual ΔG\Delta G, and coupling through shared intermediates lets a favorable reaction's negative ΔG\Delta G absorb an unfavorable reaction's positive ΔG\Delta G so the net process is spontaneous.

Where this lives today

Every metabolic pathway diagram in modern biochemistry is, at bottom, a free-energy coupling diagram: glycolysis, the citric acid cycle, oxidative phosphorylation, and fatty acid synthesis are all sequences engineered so unfavorable steps are paired with favorable ones through shared ATP/ADP, NAD⁺/NADH, or acyl-CoA intermediates. Pharmaceutical enzyme inhibitor design routinely uses ΔG∘′\Delta G^{\circ\prime} and Keq′K'_{eq} data to predict whether a target reaction is likely to be reversible in vivo (a good drug target is often one sitting at a large, physiologically irreversible ΔG\Delta G, like hexokinase or phosphofructokinase, since blocking it commits the cell rather than merely nudging an equilibrium). Metabolic engineering and synthetic biology use flux-balance and thermodynamic feasibility analysis (software such as OptMDFpathway or the "max-min driving force" method) built directly on the ΔG=ΔG∘′+RTln⁡Q\Delta G = \Delta G^{\circ\prime} + RT\ln Q relationship to redesign microbial pathways for biofuel or pharmaceutical precursor production, checking that every engineered step has ΔG<0\Delta G < 0 under achievable intracellular metabolite concentrations. In exercise physiology and clinical medicine, the actual cellular ΔG\Delta G of ATP hydrolysis (not the textbook −30.5 kJ/mol) is measured noninvasively by ³¹P NMR spectroscopy in muscle, since it governs whether ATP-dependent processes like the sodium-potassium pump or myosin cross-bridge cycling can be sustained under ischemia or intense exertion. Even the chemiosmotic proton-motive force that drives ATP synthase is analyzed with the same ΔG=ΔG∘′+RTln⁡Q\Delta G = \Delta G^{\circ\prime} + RT\ln Q machinery, generalized to include an electrical term for charged species crossing a membrane — the electrochemical potential — making this single equation the quantitative backbone connecting Gibbs's 1876 thermodynamics to the operation of a living mitochondrion.

:::quiz[If [ATP], [ADP], and [Pi] are all set to exactly 1 M, what is ΔG for ATP hydrolysis?]

  • Always zero, since the reaction is at equilibrium
  • Equal to ΔG°′, since Q = 1 and RT ln Q = 0 ✓
  • Undefined, since 1 M is never physiological
  • More negative than ΔG°′, since concentrations are high

Ask the tutor: How does this free-energy coupling framework extend to reactions that move charged species across a membrane, where an electrical potential term must be added to ΔG\Delta G?

The physical picture — free energy is a moving target, not a label

The single most important habit to build is refusing to ask "is this reaction spontaneous?" without immediately asking "at what concentrations?" Gibbs free energy change, ΔG\Delta G, is the criterion for spontaneity at constant temperature and pressure — the conditions of a living cell. A process with ΔG<0\Delta G < 0 is exergonic and proceeds spontaneously as written; ΔG>0\Delta G > 0 is endergonic and requires energy input; ΔG=0\Delta G = 0 means the system sits at equilibrium. But ΔG\Delta G is a function of the instantaneous concentrations of reactants and products, and those concentrations change as the reaction proceeds — indeed ΔG→0\Delta G \to 0 is precisely what happens because the reaction proceeds, since it is running downhill towards equilibrium. ΔG∘′\Delta G^{\circ\prime}, the biochemical standard free energy change, freezes this into a single reference number: it is ΔG\Delta G evaluated when every reactant and product is at 1 M (or, for water, at its pure-liquid activity of 1; for H⁺, at pH 7 rather than 1 M — the biochemist's modification of the chemist's standard state, hence the prime). This is a bookkeeping convenience, not a claim about what is actually in a cell. No cell has 1 M glucose-6-phosphate. ΔG∘′\Delta G^{\circ\prime} is useful precisely because it is a fixed, tabulable, comparable number — but the moment you want to know whether a reaction actually runs forward in vivo, you must convert back to real concentrations. The figure below traces the logic: a reaction's free energy surface as a function of extent of reaction, showing why ΔG∘′\Delta G^{\circ\prime} (the standard-state point) and ΔG\Delta G (the actual operating point) can even have opposite sign.

Figure
Free energy versus extent of reaction: standard-state point vs. actual cellular point
G vs. extent of reaction ξ, A ⇌ B extent of reaction ξ (→ equilibrium) G (total) standard state [A]=[B]=1 M slope here = ΔG°′ (positive: uphill) cellular state [A] high, [B] low minimum: ΔG = 0 (equilibrium) pure A pure B

Notice the essential point encoded in the picture: the slope of the GG-vs-ξ\xi curve at any point is ΔG\Delta G at that composition. At the standard-state point the slope can be positive (reaction unfavorable at 1 M) while at the actual cellular composition — far from 1 M, with substrate kept high and product kept low by downstream consumption — the local slope is negative. Same reaction, same molecules, opposite sign of ΔG\Delta G, because concentration is doing the work. The second organizing idea is coupling through a shared intermediate. If reaction 1 produces species X and reaction 2 consumes the same X, the two reactions can be summed: X cancels, and the overall ΔG\Delta G is the arithmetic sum ΔG1+ΔG2\Delta G_1 + \Delta G_2. This is legal because Gibbs free energy is a state function — path-independent — so free energies of sequential steps genuinely add, exactly like ΔH\Delta H in Hess's law. If reaction 2 is enough exergonic, the sum can be negative even though reaction 1 alone is endergonic. This is not "borrowing energy" in any mystical sense; it is the same algebra as adding two signed numbers.

Figure
Coupling: an endergonic step summed with an exergonic step through a shared intermediate
Energetic coupling via shared intermediate Pi / ATP glucose + Pi → G6P + H₂O ΔG°′ = +13.8 kJ/mol ATP + H₂O → ADP + Pi ΔG°′ = −30.5 kJ/mol shared: Pi Sum: glucose + ATP → G6P + ADP ΔG°′ = 13.8 − 30.5 = −16.7 kJ/mol

The core ideas — the governing equations

The central relationship connects the actual free energy change to the standard free energy change and the instantaneous composition through the reaction quotient QQ:

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

where:

  • ΔG\Delta G = actual free energy change at the current concentrations (kJ/mol) — the true arbiter of spontaneity,
  • ΔG∘′\Delta G^{\circ\prime} = biochemical standard free energy change (1 M solutes, pH 7, 298 K, 1 atm),
  • RR = 8.314 J mol⁻¹K⁻¹, TT = absolute temperature (K),
  • QQ = the reaction quotient, [products]coeffs/[reactants]coeffs[\text{products}]^{\text{coeffs}}/[\text{reactants}]^{\text{coeffs}}, evaluated at the actual instantaneous concentrations (water and H⁺ excluded/folded into the prime convention, solids/pure liquids at activity 1).

Reading: the first term is a fixed offset (the reaction's intrinsic thermodynamic preference at the reference point); the second term is how far actual conditions are displaced from that reference. Because it enters as RTln⁡QRT\ln Q, doubling or halving concentrations produces only a modest, logarithmic shift in ΔG\Delta G — but the shift is unbounded as Q→0Q\to 0 or Q→∞Q\to\infty, meaning sufficiently extreme concentration ratios can force ΔG\Delta G to any sign, however unfavorable ΔG∘′\Delta G^{\circ\prime} is. At equilibrium ΔG=0\Delta G = 0 by definition, and QQ takes its equilibrium value Keq′K'_{eq}, giving the second core relationship:

0=ΔG∘′+RTln⁡Keq′⟹ΔG∘′=−RTln⁡Keq′0 = \Delta G^{\circ\prime} + RT\ln K'_{eq} \quad\Longrightarrow\quad \Delta G^{\circ\prime} = -RT\ln K'_{eq}

where:

  • Keq′K'_{eq} = the equilibrium constant under biochemical standard-state conventions (concentrations at equilibrium, water/H⁺ folded in as above).

Reading: ΔG∘′\Delta G^{\circ\prime} and Keq′K'_{eq} carry exactly the same information — one is the log of the other, scaled by RTRT. A large positive Keq′K'_{eq} (products favored at equilibrium) corresponds to a large negative ΔG∘′\Delta G^{\circ\prime}, and vice versa. This equation is how tables of ΔG∘′\Delta G^{\circ\prime} are actually constructed — nobody measures free energy directly; they measure equilibrium concentrations and back out ΔG∘′\Delta G^{\circ\prime}. The third relationship is additivity for coupled/sequential reactions:

ΔGoverall=∑iΔGiΔGoverall∘′=∑iΔGi∘′\Delta G_{\text{overall}} = \sum_i \Delta G_i \qquad \Delta G^{\circ\prime}_{\text{overall}} = \sum_i \Delta G^{\circ\prime}_i

where:

  • the sum runs over each elementary or overall step in a sequence whose intermediates cancel when the equations are added.

Reading: because GG is a state function, this additivity is exact and requires no correction term — unlike, say, naively adding equilibrium constants (which must be multiplied, not added; Keq,overall′=∏iKeq,i′K'_{eq,\text{overall}} = \prod_i K'_{eq,i}, consistent with the log relationship above).

Derivation AFrom QQ-dependence to the equilibrium-constant relation

Goal. Derive ΔG=ΔG∘′+RTln⁡Q\Delta G = \Delta G^{\circ\prime} + RT\ln Q from the chemical potentials of reactants and products, then obtain ΔG∘′=−RTln⁡Keq′\Delta G^{\circ\prime} = -RT\ln K'_{eq} as the special case ΔG=0\Delta G = 0.

Step 1 — chemical potential of each species. For a reaction aA+bB⇌cC+dDaA + bB \rightleftharpoons cC + dD, each species ii present in dilute solution has a chemical potential

μi=μi∘+RTln⁡ai\mu_i = \mu_i^{\circ} + RT\ln a_i

where μi∘\mu_i^{\circ} is the chemical potential at unit activity (ai=1a_i = 1, i.e., the standard state) and aia_i is the activity (approximated here by concentration in M, relative to the 1 M reference, so aia_i is dimensionless).

Step 2 — free energy of reaction as a sum of chemical potentials. The free energy change for the reaction as written is the weighted difference between product and reactant potentials:

ΔG=[c μC+d μD]−[a μA+b μB]\Delta G = \left[c\,\mu_C + d\,\mu_D\right] - \left[a\,\mu_A + b\,\mu_B\right]

Step 3 — substitute Step 1 into Step 2 and collect standard-state terms. Expanding each μi\mu_i:

ΔG=[c μC∘+d μD∘−a μA∘−b μB∘]+RT[cln⁡aC+dln⁡aD−aln⁡aA−bln⁡aB]\Delta G = \left[c\,\mu_C^{\circ} + d\,\mu_D^{\circ} - a\,\mu_A^{\circ} - b\,\mu_B^{\circ}\right] + RT\left[c\ln a_C + d\ln a_D - a\ln a_A - b\ln a_B\right]

The bracketed standard-state term is by definition ΔG∘′\Delta G^{\circ\prime} (the free energy change when every species sits at unit activity). The logarithmic bracket collapses using ln⁡xn=nln⁡x\ln x^n = n\ln x and ln⁡(xy)=ln⁡x+ln⁡y\ln(xy)=\ln x+\ln y into a single logarithm of a ratio:

ΔG=ΔG∘′+RTln⁡(aC c aD daA a aB b)=ΔG∘′+RTln⁡Q\Delta G = \Delta G^{\circ\prime} + RT\ln\left(\frac{a_C^{\,c}\,a_D^{\,d}}{a_A^{\,a}\,a_B^{\,b}}\right) = \Delta G^{\circ\prime} + RT\ln Q

which is the first core equation, with Q≡aC caD d/(aA aaB b)Q \equiv a_C^{\,c}a_D^{\,d}/(a_A^{\,a}a_B^{\,b}) evaluated at the actual instantaneous concentrations.

Step 4 — impose the equilibrium condition. At equilibrium the system has stopped moving downhill in GG; formally this is the condition ΔG=0\Delta G = 0 (the tangent to the GG-vs-ξ\xi curve is flat — recall the minimum in the earlier figure). Substituting ΔG=0\Delta G = 0 and renaming QQ at equilibrium as Keq′K'_{eq}:

0=ΔG∘′+RTln⁡Keq′0 = \Delta G^{\circ\prime} + RT\ln K'_{eq}

Rearranging isolates the standard free energy change:

ΔG∘′=−RTln⁡Keq′\Delta G^{\circ\prime} = -RT\ln K'_{eq}

This is exact — no approximation was introduced beyond activity ≈ concentration, which is the standard biochemical convention. The sign is essential: a favorable equilibrium (Keq′>1K'_{eq} > 1, products dominate) gives ln⁡Keq′>0\ln K'_{eq} > 0 and hence ΔG∘′<0\Delta G^{\circ\prime} < 0, matching intuition.

Go deeper · Derivation B · 1 minAdditivity of coupled reactions and the free-energy "budget" for ATP-driven biosynthesis
Derivation BAdditivity of coupled reactions and the free-energy "budget" for ATP-driven biosynthesis

Goal. Show why free energies of two reactions sharing an intermediate add algebraically, and use this to derive the condition under which ATP hydrolysis can drive an endergonic biosynthetic step.

Step 1 — write the two half-reactions and their ΔG∘′\Delta G^{\circ\prime}. Consider glucose phosphorylation (unfavorable alone) and ATP hydrolysis (favorable):

Reaction 1: glucose+Pi→glucose-6-phosphate+H2OΔG1∘′=+13.8 kJ/mol\text{Reaction 1: glucose} + P_i \rightarrow \text{glucose-6-phosphate} + H_2O \qquad \Delta G^{\circ\prime}_1 = +13.8\ \text{kJ/mol}
Reaction 2: ATP+H2O→ADP+PiΔG2∘′=−30.5 kJ/mol\text{Reaction 2: ATP} + H_2O \rightarrow \text{ADP} + P_i \qquad \Delta G^{\circ\prime}_2 = -30.5\ \text{kJ/mol}

Step 2 — add the reactions as chemical equations, canceling species common to both sides. PiP_i appears as a product of reaction 2 and a reactant of reaction 1; one H2OH_2O appears as reactant of reaction 1 and product of reaction 2. Adding and cancelling:

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

Step 3 — justify additivity of ΔG∘′\Delta G^{\circ\prime} from Hess's-law logic. Because GG (and hence ΔG∘′\Delta G^{\circ\prime}, a difference of GG values at fixed standard states) is a state function, the free energy change of the sum of two equations equals the sum of the free energy changes of the individual equations — this holds regardless of mechanism or of whether the intermediate (PiP_i) is ever released into free solution, exactly as ΔH\Delta H sums in Hess's law of heat summation. Formally, writing ΔGi∘′=Gproducts,i∘′−Greactants,i∘′\Delta G^{\circ\prime}_i = G^{\circ\prime}_{\text{products},i} - G^{\circ\prime}_{\text{reactants},i} for each step, summing over ii telescopes: intermediate species' G∘′G^{\circ\prime} contributions appear once as a product (positive) and once as a reactant (negative) and cancel, leaving only the net reactants and products of the combined equation.

Step 4 — compute the coupled ΔG∘′\Delta G^{\circ\prime}.

ΔGoverall∘′=ΔG1∘′+ΔG2∘′=(+13.8)+(−30.5)\Delta G^{\circ\prime}_{\text{overall}} = \Delta G^{\circ\prime}_1 + \Delta G^{\circ\prime}_2 = (+13.8) + (-30.5)
ΔGoverall∘′=−16.7 kJ/mol\Delta G^{\circ\prime}_{\text{overall}} = -16.7\ \text{kJ/mol}

Step 5 — state the general coupling condition. An endergonic step with standard free energy +ΔGendergonic∘′+\Delta G^{\circ\prime}_{\text{endergonic}} can in principle be driven by coupling to ATP hydrolysis provided

ΔGATP hydrolysis∘′+ΔGendergonic∘′<0⟺ΔGendergonic∘′<30.5 kJ/mol\Delta G^{\circ\prime}_{\text{ATP hydrolysis}} + \Delta G^{\circ\prime}_{\text{endergonic}} < 0 \quad\Longleftrightarrow\quad \Delta G^{\circ\prime}_{\text{endergonic}} < 30.5\ \text{kJ/mol}

so ATP hydrolysis (−30.5 kJ/mol standard) can in principle offset any single endergonic step whose own standard cost is below about 30.5 kJ/mol — and, crucially, because the actual cellular ΔG\Delta G of ATP hydrolysis is substantially more negative than −30.5 kJ/mol (worked below), the true energetic budget available is even larger. This is why almost all "uphill" biosynthetic steps in metabolism are written as ATP-consuming: the coupling budget is generous, not fine-tuned.

Worked examples

Example 1 — Foundational: computing actual ΔG for ATP hydrolysis inside a real cell. ATP hydrolysis, ATP+H2O→ADP+Pi\text{ATP} + H_2O \rightarrow \text{ADP} + P_i, has ΔG∘′=−30.5 kJ/mol\Delta G^{\circ\prime} = -30.5\ \text{kJ/mol} at 298 K. Typical intracellular concentrations (asymmetric, realistic, non-round numbers, as measured in hepatocytes) are

[ATP]=3.4×10−3 M,[ADP]=1.3×10−4 M,[Pi]=4.8×10−3 M[\text{ATP}] = 3.4\times10^{-3}\ \text{M}, \quad [\text{ADP}] = 1.3\times10^{-4}\ \text{M}, \quad [P_i] = 4.8\times10^{-3}\ \text{M}

Water is the solvent and is folded into the standard-state convention (activity 1), so it does not appear in QQ. The reaction quotient is

Q=[ADP][Pi][ATP]=(1.3×10−4)(4.8×10−3)3.4×10−3Q = \frac{[\text{ADP}][P_i]}{[\text{ATP}]} = \frac{(1.3\times10^{-4})(4.8\times10^{-3})}{3.4\times10^{-3}}

Numerator: (1.3×10−4)(4.8×10−3)=6.24×10−7(1.3\times10^{-4})(4.8\times10^{-3}) = 6.24\times10^{-7}.

Q=6.24×10−73.4×10−3=1.835×10−4Q = \frac{6.24\times10^{-7}}{3.4\times10^{-3}} = 1.835\times10^{-4}

Now apply the master equation with R=8.314 J mol−1K−1R = 8.314\ \text{J mol}^{-1}\text{K}^{-1}, T=310 KT = 310\ \text{K} (physiological body temperature, not 298 K — another asymmetric, deliberately non-round choice):

RTln⁡Q=(8.314)(310)ln⁡(1.835×10−4)RT\ln Q = (8.314)(310)\ln(1.835\times10^{-4})

RT=2577 J/mol=2.577 kJ/molRT = 2577\ \text{J/mol} = 2.577\ \text{kJ/mol}. ln⁡(1.835×10−4)=ln⁡(1.835)+ln⁡(10−4)=0.6070−9.2103=−8.603\ln(1.835\times10^{-4}) = \ln(1.835) + \ln(10^{-4}) = 0.6070 - 9.2103 = -8.603.

RTln⁡Q=2.577×(−8.603)=−22.17 kJ/molRT\ln Q = 2.577\times(-8.603) = -22.17\ \text{kJ/mol}
ΔG=ΔG∘′+RTln⁡Q=−30.5+(−22.17)=−52.7 kJ/mol\Delta G = \Delta G^{\circ\prime} + RT\ln Q = -30.5 + (-22.17) = -52.7\ \text{kJ/mol}

The actual free energy released by ATP hydrolysis under these cellular conditions, ≈−52.7\approx -52.7 kJ/mol, is nearly double the standard value. Limiting-case check: if we artificially set all three concentrations to exactly 1 M, then Q=(1)(1)/1=1Q = (1)(1)/1 = 1, ln⁡Q=0\ln Q = 0, and ΔG=ΔG∘′=−30.5\Delta G = \Delta G^{\circ\prime} = -30.5 kJ/mol exactly — confirming the formula correctly reduces to the definition of the standard state when concentrations are standard. The cell deliberately keeps [ATP] high and [ADP][Pi] comparatively low (via continuous respiration and coupled consumption), which is exactly the displacement from equilibrium that makes Q≪1Q \ll 1 and drives ΔG\Delta G strongly negative — the cell is not at equilibrium, and staying away from equilibrium is the whole point.

Example 2 — Applied: can ATP hydrolysis drive the first step of glycolysis under real hepatocyte concentrations? The hexokinase reaction, glucose + ATP → glucose-6-phosphate + ADP, has ΔG∘′=−16.7\Delta G^{\circ\prime} = -16.7 kJ/mol (computed above by coupling). But we should check the actual ΔG\Delta G, not just the standard value, using realistic concentrations:

[glucose]=5.0×10−3 M,[G6P]=8.3×10−5 M,[ATP]=3.4×10−3 M,[ADP]=1.3×10−4 M[\text{glucose}] = 5.0\times10^{-3}\ \text{M},\quad [\text{G6P}] = 8.3\times10^{-5}\ \text{M},\quad [\text{ATP}] = 3.4\times10^{-3}\ \text{M},\quad [\text{ADP}] = 1.3\times10^{-4}\ \text{M}
Q=[G6P][ADP][glucose][ATP]=(8.3×10−5)(1.3×10−4)(5.0×10−3)(3.4×10−3)Q = \frac{[\text{G6P}][\text{ADP}]}{[\text{glucose}][\text{ATP}]} = \frac{(8.3\times10^{-5})(1.3\times10^{-4})}{(5.0\times10^{-3})(3.4\times10^{-3})}

Numerator: (8.3×10−5)(1.3×10−4)=1.079×10−8(8.3\times10^{-5})(1.3\times10^{-4}) = 1.079\times10^{-8}. Denominator: (5.0×10−3)(3.4×10−3)=1.70×10−5(5.0\times10^{-3})(3.4\times10^{-3}) = 1.70\times10^{-5}.

Q=1.079×10−81.70×10−5=6.35×10−4Q = \frac{1.079\times10^{-8}}{1.70\times10^{-5}} = 6.35\times10^{-4}

At T=310T = 310 K, ln⁡(6.35×10−4)=ln⁡(6.35)+ln⁡(10−4)=1.8485−9.2103=−7.362\ln(6.35\times10^{-4}) = \ln(6.35) + \ln(10^{-4}) = 1.8485 - 9.2103 = -7.362.

RTln⁡Q=2.577×(−7.362)=−18.97 kJ/molRT\ln Q = 2.577\times(-7.362) = -18.97\ \text{kJ/mol}
ΔG=−16.7+(−18.97)=−35.7 kJ/mol\Delta G = -16.7 + (-18.97) = -35.7\ \text{kJ/mol}

The coupled reaction is strongly exergonic in vivo — even more favorable than its already-favorable standard value — confirming that hexokinase runs essentially irreversibly under cellular conditions. This is precisely why hexokinase is a committed, regulated step in glycolysis: its large negative ΔG\Delta G (not merely ΔG∘′\Delta G^{\circ\prime}) makes the reaction physiologically one-directional, a hallmark of the pathway's rate-controlling steps.

Example 3 — Exam/competition level: back-calculating an unknown Keq′K'_{eq} from a multi-step coupled pathway and checking self-consistency. Suppose a hypothetical two-step pathway is:

Step A: X→Y,ΔGA∘′=+21.4 kJ/mol\text{Step A: } X \rightarrow Y, \qquad \Delta G^{\circ\prime}_A = +21.4\ \text{kJ/mol}
Step B: Y→Z,ΔGB∘′=−9.6 kJ/mol\text{Step B: } Y \rightarrow Z, \qquad \Delta G^{\circ\prime}_B = -9.6\ \text{kJ/mol}

(a) Find ΔG∘′\Delta G^{\circ\prime} for the overall process X→ZX \rightarrow Z, and the corresponding Keq′K'_{eq} at T=298T = 298 K. By additivity (justified in the second derivation above, since Y is the shared intermediate and cancels):

ΔGoverall∘′=ΔGA∘′+ΔGB∘′=21.4+(−9.6)=+11.8 kJ/mol\Delta G^{\circ\prime}_{\text{overall}} = \Delta G^{\circ\prime}_A + \Delta G^{\circ\prime}_B = 21.4 + (-9.6) = +11.8\ \text{kJ/mol}

Solve ΔG∘′=−RTln⁡Keq′\Delta G^{\circ\prime} = -RT\ln K'_{eq} for Keq′K'_{eq}:

ln⁡Keq′=−ΔG∘′RT=−11.8×103(8.314)(298)\ln K'_{eq} = -\frac{\Delta G^{\circ\prime}}{RT} = -\frac{11.8\times10^3}{(8.314)(298)}

RT=2478.6RT = 2478.6 J/mol =2.479= 2.479 kJ/mol.

ln⁡Keq′=−11.82.479=−4.761\ln K'_{eq} = -\frac{11.8}{2.479} = -4.761
Keq′=e−4.761=8.55×10−3K'_{eq} = e^{-4.761} = 8.55\times10^{-3}

So at equilibrium, [Z]/[X] ≈0.0086\approx 0.0086 — the overall reaction, taken in isolation at standard concentrations, disfavors Z. This alone would predict the pathway does not run forward.

(b) Now suppose Z is continuously siphoned off by a third, highly exergonic reaction that keeps [Z]/[X]=2.0×10−4[Z]/[X] = 2.0\times10^{-4} in the steady-state cell (well below the equilibrium ratio of 8.55×10−38.55\times10^{-3}). Compute the actual ΔG\Delta G for X→ZX \rightarrow Z under these conditions and interpret.

Q=[Z][X]=2.0×10−4Q = \frac{[Z]}{[X]} = 2.0\times10^{-4}
ln⁡Q=ln⁡(2.0×10−4)=ln⁡(2.0)−4ln⁡(10)=0.6931−9.2103=−8.517\ln Q = \ln(2.0\times10^{-4}) = \ln(2.0) - 4\ln(10) = 0.6931 - 9.2103 = -8.517
RTln⁡Q=2.479×(−8.517)=−21.12 kJ/molRT\ln Q = 2.479\times(-8.517) = -21.12\ \text{kJ/mol}
ΔG=ΔG∘′+RTln⁡Q=11.8+(−21.12)=−9.32 kJ/mol\Delta G = \Delta G^{\circ\prime} + RT\ln Q = 11.8 + (-21.12) = -9.32\ \text{kJ/mol}

Even though ΔG∘′\Delta G^{\circ\prime} for X → Z is positive (+11.8 kJ/mol, thermodynamically unfavorable at standard concentrations), the actual ΔG\Delta G under the cell's maintained steady-state concentrations is negative (−9.3 kJ/mol): the pathway runs forward. Limiting-case check: if the downstream consumption of Z were switched off and [Z]/[X][Z]/[X] were allowed to rise back to the equilibrium ratio Keq′=8.55×10−3K'_{eq} = 8.55\times10^{-3}, then Q=Keq′Q = K'_{eq} and ΔG=ΔG∘′+RTln⁡Keq′=ΔG∘′−ΔG∘′=0\Delta G = \Delta G^{\circ\prime} + RT\ln K'_{eq} = \Delta G^{\circ\prime} - \Delta G^{\circ\prime} = 0 exactly, recovering equilibrium as required. This example is the general mechanism by which "thermodynamically unfavorable" pathway steps are made to run in cells: not by changing ΔG∘′\Delta G^{\circ\prime}, but by never letting QQ approach Keq′K'_{eq}.

Limiting cases, scope & common traps

SituationWhat happens to the equationsInterpretation
Q=1Q = 1 (all species at unit activity)ΔG=ΔG∘′\Delta G = \Delta G^{\circ\prime}Recovers the definition of standard state exactly
Q=Keq′Q = K'_{eq}ΔG=0\Delta G = 0System at equilibrium; no net flux either direction
Q→0Q \to 0 (products depleted)RTln⁡Q→−∞RT\ln Q \to -\infty, ΔG→−∞\Delta G \to -\inftyReaction becomes arbitrarily favorable — cells exploit this by continuous product removal
Q→∞Q \to \infty (products flooded)RTln⁡Q→+∞RT\ln Q \to +\infty, ΔG→+∞\Delta G \to +\inftyReaction is forced backward regardless of ΔG∘′\Delta G^{\circ\prime}
T→0T \to 0RTln⁡Q→0RT\ln Q \to 0ΔG→ΔG∘′\Delta G \to \Delta G^{\circ\prime} regardless of QQ (entropic driving vanishes)
Coupling two reactions, shared intermediateΔGtot∘′=∑ΔGi∘′\Delta G^{\circ\prime}_{\text{tot}} = \sum \Delta G^{\circ\prime}_iValid always for ΔG∘′\Delta G^{\circ\prime} and for actual ΔG\Delta G of the summed reaction, provided the intermediate concentration is the same value in both half-reactions as written
Multi-step Keq′K'_{eq}Keq,tot′=∏iKeq,i′K'_{eq,\text{tot}} = \prod_i K'_{eq,i}Constants multiply (from exponentiating the additive ΔG∘′\Delta G^{\circ\prime} sum) — never add equilibrium constants

The chemist behind it

Josiah Willard Gibbs (1839–1903) spent nearly his entire career at Yale, publishing the foundational papers of chemical thermodynamics — "On the Equilibrium of Heterogeneous Substances" (1876, 1878) — in the relatively obscure Transactions of the Connecticut Academy of Arts and Sciences, so quietly that it took a decade for European scientists like Maxwell and Ostwald to recognize the work's importance and popularize it. Gibbs was trained as an engineer (the first American PhD in engineering, from Yale in 1863) and brought an engineer's insistence on operationally defined, measurable quantities to what had been a physicist's and chemist's somewhat ad hoc subject. His free energy function — what we write today as G=H−TSG = H - TS — provided exactly the single scalar quantity needed to predict, from tabulated data alone, whether any process at constant temperature and pressure would proceed spontaneously, without needing to separately track heat and disorder. Gibbs never applied his function to a living cell; that synthesis, connecting ΔG∘′\Delta G^{\circ\prime}, Keq′K'_{eq}, and the biochemical standard state, was worked out over the following seven decades by biochemists including Fritz Lipmann, who in 1941 crystallized the idea of "high-energy phosphate bonds" and energy coupling that makes ATP the currency it is.

Origins · 1 minJosiah Willard Gibbs and how this was discovered
JW
Josiah Willard Gibbs
1839–1903 · Yale University; founder of chemical thermodynamics and statistical mechanics

Gibbs formalized the free energy function that bears his name and derived the phase rule, giving chemistry its first fully general, quantitative criterion for spontaneity and equilibrium — the mathematical foundation on which every subsequent biochemical energy calculation, including ATP coupling, ultimately rests.

History
- 1873 — Gibbs introduces graphical thermodynamic surfaces relating energy, entropy, and volume.
- 1876/1878 — Gibbs publishes "On the Equilibrium of Heterogeneous Substances," defining the free energy function and chemical potential.
- 1923 — Gilbert N. Lewis and Merle Randall popularize activities and the modern chemical-thermodynamics formalism (ΔG=ΔG∘+RTln⁡Q\Delta G = \Delta G^{\circ} + RT\ln Q) in their landmark textbook.
- 1941 — Fritz Lipmann proposes the concept of the "high-energy phosphate bond" and articulates ATP as the common intermediate coupling catabolism to biosynthesis.
- 1957 — Systematic tables of biochemical standard free energies of hydrolysis (ATP, other phosphorylated compounds) are compiled, enabling routine ΔG∘′\Delta G^{\circ\prime}-based pathway analysis.
Why does biochemistry use ΔG°′ instead of the chemist's ΔG°?

Chemists define the standard state with all species (including H⁺) at 1 M activity. But cellular pH is about 7, meaning [H⁺] ≈ 10⁻⁷ M, not 1 M — using 1 M H⁺ would make every biological standard state absurdly acidic and physiologically meaningless. Biochemists therefore define a modified standard state, denoted with a prime, in which pH is fixed at 7 (and water's concentration, ~55.5 M, is folded into the constant rather than treated as a variable) while all other solutes remain at 1 M. ΔG∘′\Delta G^{\circ\prime} values in biochemistry tables always use this convention; converting between ΔG∘\Delta G^{\circ} and ΔG∘′\Delta G^{\circ\prime} for a reaction involving H⁺ requires knowing the stoichiometric number of protons and applying the corresponding RTln⁡[H+]RT\ln[\text{H}^+] correction.

Is it true that ATP has a special "high-energy bond" that stores energy?

Not in the sense of a uniquely strong or strained covalent bond storing potential energy like a coiled spring. The large negative ΔG∘′\Delta G^{\circ\prime} of ATP hydrolysis arises from several cooperating factors: electrostatic repulsion among the closely spaced negative charges on the triphosphate chain is relieved upon hydrolysis, the products (ADP²⁻ and Pi²⁻) are stabilized by greater resonance delocalization than the triphosphate, and both products are more effectively solvated by water than ATP is. "High-energy bond" is a historically entrenched shorthand (coined by Lipmann) for "a bond whose hydrolysis has a large negative ΔG∘′\Delta G^{\circ\prime}," not a statement about unusual bond strength — in fact the P–O bond being broken is not unusually weak in a bond-dissociation-energy sense.

Can a reaction with a favorable ΔG (negative) still fail to occur at any appreciable rate?

Yes, and this is a crucial distinction from kinetics. ΔG\Delta G tells you the thermodynamic feasibility and direction of a reaction — whether it can happen — but says nothing about the rate at which it happens, which is governed by the activation energy barrier and is the domain of enzyme kinetics, not thermodynamics. The hydrolysis of ATP to ADP and Pi is enormously exergonic under cellular conditions, yet ATP is kinetically stable in water for hours in the absence of a catalyst, because the uncatalyzed reaction has a substantial activation barrier. This kinetic stability is itself biologically essential — ATP would be useless as an energy currency if it decomposed spontaneously the instant it was made; enzymes selectively lower the barrier only when and where the energy is needed.

How is ΔG°′ actually measured or looked up for a given reaction?

In practice, ΔG∘′\Delta G^{\circ\prime} is rarely measured calorimetrically. Instead, the reaction is allowed to reach equilibrium in vitro (often catalyzed by the relevant enzyme to reach equilibrium quickly, since the enzyme changes only rate, not the equilibrium position), the equilibrium concentrations of all reactants and products are measured (by spectrophotometry, HPLC, or coupled enzymatic assays), Keq′K'_{eq} is computed from those concentrations, and ΔG∘′=−RTln⁡Keq′\Delta G^{\circ\prime} = -RT\ln K'_{eq} is then calculated. Comprehensive tables (originally compiled by researchers such as Robert Alberty) list ΔGf∘′\Delta G^{\circ\prime}_f (standard free energy of formation) for individual metabolites, from which the ΔG∘′\Delta G^{\circ\prime} of any reaction connecting them can be computed as a difference of formation energies, analogous to using ΔHf∘\Delta H^{\circ}_f tables in general chemistry.

Does energetic coupling require the two reactions to occur in the same active site?

No — thermodynamic coupling via a shared intermediate works regardless of whether the two half-reactions occur at the same active site, in different active sites of the same enzyme, or in entirely different enzymes and even different cellular compartments, as long as the intermediate species (ATP, NADH, a proton gradient, an acyl-CoA thioester) is genuinely shared between them. However, kinetic efficiency is often far higher when coupling happens within one active site (as in hexokinase, where phosphoryl transfer occurs directly from ATP to glucose without ever releasing free phosphate into solution), since this avoids diffusional loss and unwanted side reactions of a reactive intermediate. Both mechanisms — direct group transfer in one active site, and diffusible shared intermediates connecting separate reactions — are thermodynamically identical in their free-energy bookkeeping.

Why is it wrong to add K'eq values for sequential reactions the way ΔG°′ values are added?

Because ΔG∘′=−RTln⁡Keq′\Delta G^{\circ\prime} = -RT\ln K'_{eq} is a logarithmic relationship: summing ΔG∘′\Delta G^{\circ\prime} values corresponds to summing logarithms of Keq′K'_{eq} values, which by the laws of logarithms means the Keq′K'_{eq} values themselves must be multiplied, not added, to get the overall equilibrium constant. Concretely, Keq,overall′=∏iKeq,i′K'_{eq,\text{overall}} = \prod_i K'_{eq,i} follows directly from exponentiating ΔGoverall∘′=∑iΔGi∘′=−RT∑iln⁡Keq,i′=−RTln⁡(∏iKeq,i′)\Delta G^{\circ\prime}_{\text{overall}} = \sum_i \Delta G^{\circ\prime}_i = -RT\sum_i \ln K'_{eq,i} = -RT\ln\left(\prod_i K'_{eq,i}\right). Adding equilibrium constants directly has no thermodynamic justification and will generally give a badly wrong answer, especially when the constants span many orders of magnitude, as biochemical ones often do.

Before this

Where this leads