Biochem

Cooperative binding: hemoglobin and the Hill equation

12 min read Updated

Michaelis-Menten kinetics assumes each enzyme active site works independently — the binding of a substrate to one enzyme molecule has no effect on any other. Half the interesting biology violates this assumption. Hemoglobin binds four oxygen molecules, but each successive O₂ binds more easily than the previous one; aspartate transcarbamoylase’s catalytic activity switches on across a narrow window of substrate concentration; allosteric enzymes at metabolic branch points respond to their substrate as a sigmoidal switch, not a hyperbolic drift. Cooperativity is the mechanism. The Hill equation is the arithmetic. Together they explain why hemoglobin can pick up oxygen in the lung and release it in tissue with the same 20-mmHg change in P_O₂ that would move a non-cooperative protein less than 15% of the way through its response curve.

Assume familiarity with Michaelis-Menten from why enzymes plateau and pH buffering from physiological buffers.

What cooperative binding looks like

Plot fractional saturation Y (fraction of binding sites occupied) against ligand concentration [L] for a cooperative protein and you get an S-shaped curve, not a hyperbola:

  • At low [L], Y rises slowly — the first ligand binds with difficulty.
  • Over a narrow concentration window, Y rises steeply — subsequent ligands bind with rapidly increasing ease.
  • At high [L], Y plateaus at 1 — full saturation.

Hemoglobin’s oxygen-binding curve (Lehninger 7e Fig 5-12) shifts from 4% saturation at P_O₂ = 5 mmHg (a low-oxygen tissue) to 98% saturation at P_O₂ = 100 mmHg (the alveoli of the lung) — an 86 percentage-point change over a 20× concentration window. Myoglobin, which is not cooperative, is 74% saturated at P_O₂ = 5 mmHg and 99% at 100 mmHg — it hangs onto oxygen so tightly that it never releases to tissue. Hemoglobin’s cooperativity is precisely what makes it a good oxygen carrier rather than a good oxygen sponge.

The Hill equation

Archibald Hill’s 1910 empirical fit to the hemoglobin curve is:

Y=[L]nKdn+[L]nY = \frac{[\text{L}]^n}{K_d^n + [\text{L}]^n}

with:

  • n = Hill coefficient (dimensionless).
  • K_d = apparent binding constant, roughly the [L] where Y = 0.5.
  • [L] = ligand concentration.

At n = 1, Hill collapses to the standard hyperbolic binding curve (equivalent to Michaelis-Menten shape with K_d ≡ K_m). At n > 1, the curve is positively cooperative — steeper than hyperbolic. Larger n → sharper switch. At n < 1, the protein is negatively cooperative, less common but real.

The Hill coefficient n is a lower bound on the number of binding sites. Hemoglobin has four heme groups, so n_max = 4; the experimental n_Hill for hemoglobin is 2.8-3.0, meaning cooperativity is strong but not perfect (a hypothetically perfect all-or-none tetramer would give exactly n = 4). Berg 9e §7.2 gives the range of measured n for common cooperative proteins:

  • Hemoglobin (O₂): n = 2.8-3.0.
  • Aspartate transcarbamoylase (substrate): n ≈ 3.
  • Phosphofructokinase (ATP inhibition): n ≈ 3.8.
  • Calcium-binding proteins (calmodulin): n = 1.5-2 for Ca²⁺.

Deriving n from data: the Hill plot

Rearrange Hill:

log⁡(Y1−Y)=nlog⁡[L]−nlog⁡Kd\log\left(\frac{Y}{1-Y}\right) = n \log[\text{L}] - n \log K_d

A plot of log(Y/(1−Y)) vs log[L] — the Hill plot — is linear with slope n. In practice the plot is only truly linear in the middle of the transition (0.1 < Y < 0.9); at the tails it curves toward the individual K_d values for the first and last binding events. The slope at Y = 0.5 gives the maximum Hill coefficient, which is what “n_Hill” conventionally reports.

The Hill plot is the diagnostic test for cooperativity. If n > 1 at Y = 0.5, cooperativity is positive; n < 1, negative; n = 1, non-cooperative. Voet 5e §10.1 has canonical worked examples using hemoglobin data.

The MWC (Monod-Wyman-Changeux) mechanism

Hill is descriptive but not mechanistic. The mechanism that explains why hemoglobin cooperates is the Monod-Wyman-Changeux (MWC) model (1965): the tetramer exists in equilibrium between two conformational states — T (tense) and R (relaxed) — with all four subunits switching in concert. T has low O₂ affinity; R has high affinity. The equilibrium [R]/[T] shifts toward R as more oxygens bind.

Perutz’s 1970 X-ray crystallography confirmed the atomic-scale reality: T and R are two distinct quaternary structures of the tetramer that differ in the alignment of subunit-subunit contacts. Oxygen binding to one heme triggers a conformational change transmitted through those contacts, pulling nearby subunits into R and raising their O₂ affinity. The mechanism is concerted — the whole tetramer flips at once, not one subunit at a time.

The alternative model (Koshland-Némethy-Filmer, KNF, or “sequential” model) allows individual subunits to switch independently. Real hemoglobin sits between the two — the MWC picture is a good first approximation, with small deviations that KNF captures better.

Regulators of cooperativity: allosteric effectors

Cooperativity is not fixed — it responds to other molecules. Hemoglobin’s O₂ affinity is regulated by four heterotropic effectors:

  • CO₂ / H⁺ (Bohr effect): high CO₂ and low pH (as in exercising muscle) stabilize T. Hemoglobin’s O₂ affinity drops; oxygen is released more readily where it is most needed. The Bohr effect is why exercising tissue gets O₂ efficiently despite ambient P_O₂ being similar to resting tissue.
  • 2,3-BPG (bisphosphoglycerate): binds to the central cavity of the T state and stabilizes it. High-altitude adaptation involves upregulating [2,3-BPG], shifting hemoglobin’s O₂ affinity lower so more O₂ is released to tissue at reduced arterial P_O₂.
  • Temperature: higher T stabilizes R; hyperthermia decreases O₂ affinity, so oxygen is easier to release into hot exercising tissue.
  • Fetal hemoglobin (HbF): contains γ-chains instead of β-chains, binds 2,3-BPG less tightly, and thus has higher O₂ affinity than adult HbA. That difference is what allows the fetus to extract O₂ from the mother’s blood at the placenta.

Each of these effectors changes the position of the sigmoidal curve on the P_O₂ axis without necessarily changing its shape or n_Hill. The clinical shorthand is “P₅₀” — the P_O₂ at which Y = 0.5. Normal adult Hb has P₅₀ = 26 mmHg. HbF has P₅₀ = 19 mmHg. High-altitude adaptation raises P₅₀ to ~30 mmHg through 2,3-BPG upregulation.

Worked example — hemoglobin O2 saturation at rest vs exercise

Worked example — Hb O2 saturation, tissue vs lung, resting vs exercising
Hill coefficient n
2.8 (adult HbA)
P50 (rest)
26 mmHg
P50 (exercising muscle)
34 mmHg (Bohr-shifted by CO2/H+)
Alveolar P_O2
100 mmHg
Resting tissue P_O2
40 mmHg
Exercising muscle P_O2
20 mmHg
  1. Formula
  2. Substitute
  3. Note

    Compute 26^2.8 = 4776. Compute 100^2.8 = 631,000. Y_lung = 631000/(4776+631000) = 0.9925. Nearly fully saturated in the alveoli. At resting tissue P_O2 = 40 mmHg: 40^2.8 = 22150. Y_rest = 22150/(4776+22150) = 0.823 — 82% saturated. In exercising muscle, both the tissue P_O2 (20 mmHg) drops AND P50 rises to 34 mmHg (Bohr effect from CO2 and lactic acid). Compute 34^2.8 = 10570. Compute 20^2.8 = 3193. Y_exercise = 3193/(10570+3193) = 0.232 — 23% saturated. So the delivered oxygen fraction is 0.99 - 0.82 = 0.17 at rest, but 0.99 - 0.23 = 0.76 during exercise. The same hemoglobin protein, without changing its expression or concentration, delivers 4.5x more O2 to exercising muscle simply by responding to the local pH and P_CO2 through the Bohr effect. This is what cooperativity buys: extreme sensitivity to the physiological context, all encoded in the shape of the sigmoidal curve.

  4. Result

When to reach for Hill vs Michaelis-Menten

Practical decision tree:

  • Plotting rate or binding vs [L]: does the curve look hyperbolic or sigmoidal? If sigmoidal, use Hill.
  • Number of binding sites known: if 1, use Michaelis-Menten. If ≥ 2, consider Hill (the protein might still be non-cooperative).
  • Hill plot slope at Y = 0.5: n = 1 means non-cooperative (use MM); n > 1 means positively cooperative (use Hill); n < 1 means negatively cooperative (use Hill, but interpret cautiously).

Hill treats cooperativity phenomenologically; MWC and KNF give the mechanism. For teaching or curve-fitting, Hill is enough. For enzyme kinetics near a regulatory setpoint or for structure-function studies, upgrade to MWC.

Practice problems

Problem 1 Easy An allosteric enzyme has Hill coefficient n = 3 and P50 = K = 1 mM. What is the fractional activity at [S] = 0.5 mM?

Answer: Y = 0.5^3 / (1^3 + 0.5^3) = 0.125 / 1.125 = 0.111 — 11% active. Compare to Michaelis-Menten (n = 1) at the same [S]/Km = 0.5: v/Vmax = 0.5 / 1.5 = 0.333 — 33% active. The cooperative enzyme is more turned off at sub-Km substrate, sharper switch as [S] rises past K.

Problem 2 Intermediate At high altitude, [2,3-BPG] in red cells doubles. If this shifts P50 from 26 to 32 mmHg (leaving n = 2.8 unchanged), what is Y at typical high-altitude arterial P_O2 = 60 mmHg? Compare to sea-level (P_O2 = 100, P50 = 26).

Answer: At altitude: 60^2.8 / (32^2.8 + 60^2.8) = 5.51e4 / (1.71e4 + 5.51e4) = 0.763. Sea level: 100^2.8 / (26^2.8 + 100^2.8) = 6.31e5 / (0.48e4 + 6.31e5) = 0.992. Altitude Hb still delivers more O2 to tissue than sea-level Hb at the same tissue P_O2, but arterial Hb is only 76% saturated instead of 99% — the compromise is fewer O2 molecules per pass, but each is more releasable. That is the physiological rationale for polycythemia (increased red cell count) as a complementary adaptation.

Problem 3 Hard You measure the initial rate of an allosteric enzyme at [S] = 0.1, 0.5, 1.0, 2.0, 5.0 mM and get v/Vmax = 0.01, 0.14, 0.50, 0.86, 0.99. Estimate the Hill coefficient.

Answer: Compute log(Y/(1-Y)) at each: log(0.01/0.99) = -2.00; log(0.14/0.86) = -0.79; log(0.50/0.50) = 0; log(0.86/0.14) = 0.79; log(0.99/0.01) = 2.00. Compute log[S]: log(0.1) = -1; log(0.5) = -0.30; log(1) = 0; log(2) = 0.30; log(5) = 0.70. Slope of Hill plot at Y = 0.5: use the two nearby points around Y = 0.5 (0.14 and 0.86). Rise = 0.79 - (-0.79) = 1.58. Run = 0.30 - (-0.30) = 0.60. Slope = 1.58/0.60 = 2.63. Hill coefficient n ≈ 2.6, consistent with a strongly cooperative 3-4-subunit enzyme. K_d = 1 mM (where Y = 0.5).

Where to go next

  • The Michaelis-Menten calculator is not directly usable for cooperative enzymes but is a good comparison — plot v/Vmax for the same K and see how much more gradual the MM curve is than the Hill curve at n = 3.
  • The buffer calculator models the pH environment that drives the Bohr effect; simulating pH shifts under exercise loads gives an intuition for how physiological pH couples to O₂ delivery.
  • Beyond hemoglobin, cooperative binding underlies gene regulation (transcription factors on multiple binding sites), calcium signaling (calmodulin’s four EF hands), and drug-receptor pharmacology (agonists with n > 1). The Hill equation shows up throughout biology.

Cooperative binding is what lets a protein act like a switch instead of a dimmer. Michaelis-Menten describes proteins whose response is monotonic and gradual; Hill describes proteins whose response is sigmoidal and sharp. Every regulatory decision point in biology — muscle contraction, gene expression, hemoglobin loading and unloading, allosteric enzyme control — reaches for cooperativity to convert a gradual chemical signal into a decisive functional output.

Sources

  • Nelson & Cox, Lehninger Principles of Biochemistry, 7th ed. (2017), §5.1 — Oxygen-binding proteins
  • Berg, Tymoczko, Gatto & Stryer, Biochemistry, 9th ed. (2019), §7.2 — Cooperative binding
  • Voet, Voet & Pratt, Biochemistry, 5th ed. (2016), §10.1
  • M. F. Perutz (1970), Nature 228, 726 — Stereochemistry of cooperative effects in hemoglobin
  • A. V. Hill (1910), J. Physiol. 40, iv-vii — Original Hill equation for hemoglobin
  • J. Monod, J. Wyman & J.-P. Changeux (1965), J. Mol. Biol. 12, 88 — MWC model