Hemoglobin O₂ Binding: Multiple Occupancy Probabilities

⚡ Key Takeaways

Hemoglobin O₂ Binding: Multiple Occupancy Probabilities is the probabilistic treatment of the four-site cooperative oxygen-binding system in hemoglobin (Hb), in which the equilibrium occupancy distribution over the $2^4 = 16$ microstates (or, by symmetry, the five macrostates $i=0,1,2,3,4$ ) is governed by the Adair equation and its statistical-mechanical (Pauling/MWC/KNF) refinements. Key result: the fractional saturation $\bar{Y}$ and the full occupancy probability vector $P_i(p_{O_2})$ are rational functions of oxygen partial pressure built from successive (Adair) binding constants, and cooperativity manifests as a non-binomial spread of the $P_i$ . Applies in: biophysical chemistry, allosteric enzyme kinetics, respiratory physiology, ligand-binding statistical mechanics. Who/when: Gilbert S. Adair (1925), Linus Pauling (1935), Monod–Wyman–Changeux (1965), Koshland–Némethy–Filmer (1966).

💡 The Core Insight

Treat the tetramer as a small grand-canonical system whose chemical potential is set by dissolved O₂. The binding polynomial (partition function) $Z(x)$ is a degree-4 polynomial in the ligand activity $x \propto p_{O_2}$ ; every occupancy probability is one monomial of $Z$ divided by $Z$ , and every thermodynamic observable is a logarithmic derivative of $Z$ . Cooperativity = curvature of $\ln Z$ beyond the independent-site baseline.

DefinitionBinding Polynomial and Occupancy Probabilities

Let $x = K_{\text{ref}}\,[\text{O}_2]$ be a dimensionless ligand activity (or simply $x=[\text{O}_2]$ with constants carrying units). Define the Adair binding constants $K_1,\dots,K_4$ as the stepwise (macroscopic) association constants:

1

\quad \text{Hb(O}_2)_{i-1} + \text{O}_2 \;\rightleftharpoons\; \text{Hb(O}_2)_i, \qquad K_i = \frac{[\text{Hb(O}_2)_i]}{[\text{Hb(O}_2)_{i-1}]\,[\text{O}_2]}.

Derivation hint: write each step as a mass-action equilibrium; $K_i$ is the ratio of successive species concentrations per unit free ligand. Let $\beta_i = \prod_{j=1}^{i} K_j$ be the cumulative (overall) Adair constants, with $\beta_0 = 1$ . The binding polynomial is

2

\quad Z([\text{O}_2]) = \sum_{i=0}^{4} \beta_i\,[\text{O}_2]^{i} = 1 + \beta_1[\text{O}_2] + \beta_2[\text{O}_2]^2 + \beta_3[\text{O}_2]^3 + \beta_4[\text{O}_2]^4.

Derivation hint: sum the equilibrium concentrations of all macrostates relative to the unliganded reference; $Z$ is the grand partition function with $[\text{O}_2]$ playing the role of fugacity. The multiple occupancy probability of exactly $i$ bound ligands is

3

\quad P_i([\text{O}_2]) = \frac{\beta_i\,[\text{O}_2]^{i}}{Z([\text{O}_2])}, \qquad \sum_{i=0}^{4} P_i = 1.

Derivation hint: each macrostate population is $\beta_i[\text{O}_2]^i$ times the reference; normalize by the total $Z$ . Here $\beta_i$ bundles a statistical degeneracy factor $g_i = \binom{4}{i}$ (the number of microscopic ways to place $i$ ligands on 4 sites) with the intrinsic microscopic affinities. If all sites were independent with intrinsic constant $k$ , then $\beta_i = \binom{4}{i} k^i$ and $P_i$ would be exactly binomial.

TheoremFractional Saturation as a Logarithmic Derivative

The mean fractional saturation $\bar{Y} = \tfrac14\langle i\rangle$ satisfies

4

\quad \bar{Y} = \frac{1}{4}\sum_{i=0}^{4} i\,P_i = \frac{1}{4}\,\frac{[\text{O}_2]\,Z'([\text{O}_2])}{Z([\text{O}_2])} = \frac{1}{4}\,\frac{\partial \ln Z}{\partial \ln [\text{O}_2]}.

Derivation hint: $\langle i\rangle = \sum i P_i$ ; recognize $\sum i\beta_i[\text{O}_2]^i = [\text{O}_2] Z'$ , then convert to a log-derivative. Explicitly, the Adair equation:

5

\quad \bar{Y} = \frac{\beta_1[\text{O}_2] + 2\beta_2[\text{O}_2]^2 + 3\beta_3[\text{O}_2]^3 + 4\beta_4[\text{O}_2]^4}{4\left(1 + \beta_1[\text{O}_2] + \beta_2[\text{O}_2]^2 + \beta_3[\text{O}_2]^3 + \beta_4[\text{O}_2]^4\right)}.

Derivation hint: substitute $Z$ and $[\text{O}_2]Z'$ into [4] and factor the leading $\tfrac14$ .

ProofFrom Mass Action to the Occupancy Distribution

Step 1 — Enumerate macrostates. A tetramer with four equivalent sites has macrostates indexed by $i\in\{0,1,2,3,4\}$ . Choose Hb (fully deoxy) as reference, $[\text{Hb}]_0 \equiv [\text{Hb}(\text{O}_2)_0]$ .

Step 2 — Express each population. Iterating [1], $[\text{Hb}(\text{O}_2)_i] = K_i[\text{O}_2][\text{Hb}(\text{O}_2)_{i-1}] = \big(\prod_{j\le i}K_j\big)[\text{O}_2]^i[\text{Hb}]_0 = \beta_i[\text{O}_2]^i[\text{Hb}]_0$ .

Step 3 — Total tetramer concentration. Summing, $[\text{Hb}]_{\text{tot}} = [\text{Hb}]_0 \sum_{i=0}^4 \beta_i[\text{O}_2]^i = [\text{Hb}]_0\,Z$ .

Step 4 — Probabilities. $P_i = [\text{Hb}(\text{O}_2)_i]/[\text{Hb}]_{\text{tot}} = \beta_i[\text{O}_2]^i / Z$ , which is [3] and automatically normalized.

Step 5 — Mean occupancy. Bound O₂ per tetramer $= \sum_i i[\text{Hb}(\text{O}_2)_i]/[\text{Hb}]_{\text{tot}} = \sum_i iP_i$ . Since $\sum_i i\beta_i[\text{O}_2]^i = [\text{O}_2]\,\partial_{[\text{O}_2]}Z$ , divide by $Z$ to get [4]; divide by 4 sites for fractional saturation [5].

Step 6 — Microscopic decomposition. Each $\beta_i$ factors as $g_i \kappa_i$ where $g_i = \binom{4}{i}$ is the number of distinct site patterns and $\kappa_i$ is the product of intrinsic association constants of an average pattern. The independent-site null model is $\kappa_i = k^i$ , giving the binomial $P_i = \binom{4}{i}(k[\text{O}_2])^i/(1+k[\text{O}_2])^4$ . Deviations of $\kappa_i$ from $k^i$ are exactly the thermodynamic signature of cooperativity. $\blacksquare$

∎

LP

Linus Pauling

1901–1994

In 1935 Pauling produced the first statistical-mechanical theory of hemoglobin cooperativity, modeling the four hemes as interacting sites with a nearest-neighbor "interaction free energy" that boosted the affinity of empty sites once a neighbor was occupied. His treatment recast the Adair phenomenology as a microscopic lattice problem and anticipated, by three decades, the allosteric formalisms of Monod–Wyman–Changeux and Koshland–Némethy–Filmer. Pauling later identified the molecular basis of sickle-cell disease (1949) as a hemoglobin defect — coining "molecular disease" — connecting the same molecule's occupancy thermodynamics to clinical pathology.

🧠 Intuition

Picture four mailboxes on one pole; an O₂ molecule is a letter. Independent sites: each mailbox fills at random, so the count of filled boxes is binomial — broad and symmetric near half-saturation. Positive cooperativity: the first letter "unlocks" the others (the deoxy → oxy quaternary switch makes remaining sites hungrier). Now the system prefers to be either nearly empty or nearly full, and dislikes the half-filled states. The probability mass on $P_1, P_2, P_3$ is depleted while $P_0$ and $P_4$ swell — the distribution becomes bimodal-leaning and the saturation curve sigmoidal. Mathematically, cooperativity = the variance [6] exceeds the binomial variance near $\bar Y = \tfrac12$ , and the Hill coefficient $n_H = 1 + \big(\partial \ln(\bar Y/(1-\bar Y))/\partial\ln[\text{O}_2] - 1\big)$ rises above 1.

Figure

Occupancy probabilities

P_i

versus

\log p_{O_2}

for cooperative hemoglobin. Intermediate states

P_1

–

P_3

are suppressed;

P_0

and

P_4

dominate; the saturation curve

\bar Y

is sigmoidal.

ExampleComputing the Full Occupancy Vector

Given: Adair cumulative constants (illustrative, units of mmHg $^{-i}$ ): $\beta_1 = 0.10$ , $\beta_2 = 0.012$ , $\beta_3 = 0.0011$ , $\beta_4 = 0.00060$ , at $p_{O_2} = 26$ mmHg ( $\approx$ p50 region).

Step 1 — Monomials $\beta_i p^i$ :

$\beta_0 = 1$
$\beta_1 p = 0.10 \times 26 = 2.60$
$\beta_2 p^2 = 0.012 \times 676 = 8.112$
$\beta_3 p^3 = 0.0011 \times 17576 = 19.334$
$\beta_4 p^4 = 0.00060 \times 456976 = 274.19$

Step 2 — Partition function: $Z = 1 + 2.60 + 8.112 + 19.334 + 274.19 = 305.24$ .

Step 3 — Probabilities $P_i = \beta_i p^i/Z$ :

$P_0 = 1/305.24 = 0.00328$
$P_1 = 2.60/305.24 = 0.00852$
$P_2 = 8.112/305.24 = 0.02658$
$P_3 = 19.334/305.24 = 0.06334$
$P_4 = 274.19/305.24 = 0.89829$

Step 4 — Saturation: $\langle i\rangle = 0(0.00328)+1(0.00852)+2(0.02658)+3(0.06334)+4(0.89829) = 4(0.89829)+\dots = 3.844$ ; $\bar Y = 3.844/4$ .

\boxed{\;P = (0.0033,\,0.0085,\,0.0266,\,0.0633,\,0.8983),\quad \bar Y = 0.961\;}

Note the strong suppression of $P_1,P_2,P_3$ relative to $P_4$ — the hallmark of cooperativity at high $p_{O_2}$ .

ExampleVariance and Hill Slope at Half-Saturation

Given: A symmetric cooperative model with $\beta_1 = \beta_3 = 4c\,L$ small and $\beta_4 = L$ chosen so that $\bar Y = \tfrac12$ at $p = p_{50} = 26$ mmHg, with measured $P = (0.18, 0.05, 0.04, 0.05, 0.18)$ (intermediates strongly suppressed; note these sum with the rest of mass; rescaled here for illustration so $\sum P_i = 0.50$ on the empty side — treat as normalized so $P_0+\dots+P_4 = 1$ with $P_0=P_4=0.36, P_1=P_3=0.10, P_2=0.08$ ).

Step 1 — Mean: by symmetry $\langle i\rangle = 2$ , $\bar Y = 0.5$ . ✓

Step 2 — Second moment: $\langle i^2\rangle = 0^2(0.36)+1^2(0.10)+2^2(0.08)+3^2(0.10)+4^2(0.36)$ $= 0 + 0.10 + 0.32 + 0.90 + 5.76 = 7.08$ .

Step 3 — Variance: $\operatorname{Var}(i) = 7.08 - 2^2 = 3.08$ .

Step 4 — Compare to independent-site (binomial) baseline at $\bar Y=\tfrac12$ : $\operatorname{Var}_{\text{indep}} = 4\cdot\tfrac12\cdot\tfrac12 = 1.0$ .

Step 5 — Cooperativity ratio: $\operatorname{Var}/\operatorname{Var}_{\text{indep}} = 3.08$ , confirming a strongly super-binomial (positively cooperative) distribution. Using $n_H \approx \operatorname{Var}/[\,4\bar Y(1-\bar Y)\,]$ at half-saturation:

\boxed{\;\operatorname{Var}(i) = 3.08,\quad n_H \approx 3.08/(4\cdot 0.5\cdot 0.5)=3.08\;}

(Real human Hb gives $n_H \approx 2.8$ –3.0; the maximum possible for a tetramer is 4.)

⚠️ Common mistake

The Hill coefficient is not the number of binding sites. $n_H \le 4$ for hemoglobin but typically $\approx 2.8$ –3.0; it is a cooperativity index (the maximum log-log slope of $\bar Y/(1-\bar Y)$ ), not a stoichiometry. Only in the unphysical limit of infinite cooperativity ("all-or-none" binding, where $P_1=P_2=P_3=0$ ) does $n_H$ reach the site number 4. Conversely $n_H=1$ means independent sites, not one site.

⚠️ Common mistake

Adair constants are macroscopic and cannot, alone, recover microscopic site affinities. Oxygen-equilibrium curves determine only the four $K_i$ (or $\beta_i$ ). But there are eight distinct microscopic species among the singly/doubly/triply liganded states (e.g. two ways to be doubly liganded on a tetramer with two $\alpha$ and two $\beta$ chains in distinct geometries). The map from microscopic constants to macroscopic Adair constants is many-to-one; the inverse problem is underdetermined without additional data (NMR, kinetics, isotope-edited measurements).

CounterexampleIndependent Sites Mimicking Cooperativity Is Impossible, but the Reverse Trap Exists

A naive fit might report $n_H > 1$ and conclude positive cooperativity, yet heterogeneity of independent sites produces $n_H < 1$ (apparent negative cooperativity), never $n_H>1$ . Consider two non-interacting site classes with affinities $k_a \ne k_b$ (e.g. distinct $\alpha$ and $\beta$ chains, 2 each), each binding independently: $\bar Y = \tfrac12\big[\frac{k_a p}{1+k_a p} + \frac{k_b p}{1+k_b p}\big]$ . This curve is flatter than a single hyperbola, giving $n_H < 1$ . Thus a measured $n_H>1$ genuinely requires site–site interaction, but a measured $n_H<1$ is ambiguous between true negative cooperativity and mere affinity heterogeneity — a classic identifiability pitfall.

connection

The occupancy framework links several neighboring theories. The Monod–Wyman–Changeux Allosteric Model derives the $\beta_i$ from a two-state (T/R) equilibrium with allosteric constant $L$ and ratio $c=K_R/K_T$ , giving a concerted microscopic origin for the Adair constants. The Koshland–Némethy–Filmer Sequential Model instead derives them from ligand-induced conformational change propagating through subunit interfaces. Both reduce, at the level of observables, to the Adair polynomial [2], illustrating that macroscopic occupancy data underdetermine the microscopic mechanism. The variance relation [6] connects to the Fluctuation-Dissipation Theorem (occupancy fluctuations ↔ susceptibility to ligand). The binding polynomial is the Grand Canonical Partition Function of a four-site system, tying the whole subject to Statistical Mechanics of Adsorption (the Langmuir/Ising correspondence) and to Ising Model on Small Lattices via Pauling's nearest-neighbor interaction picture.

🔬 Where this appears

⏳ History

- 1904 — Christian Bohr, Hasselbalch, Krogh: report the sigmoidal oxygen dissociation curve and the Bohr effect (pH dependence), establishing that Hb O₂ binding is not hyperbolic.

- 1910 — Archibald Hill: introduces the empirical Hill equation

\bar Y = (Kp)^n/(1+(Kp)^n)

to fit the sigmoid;

n

becomes the cooperativity index (though Hill knew it was not a literal stoichiometry).

- 1925 — Gilbert S. Adair: measures hemoglobin's molecular weight (≈68 kDa) implying four hemes, and formulates the four-constant Adair equation [5] as the exact intermediate-state description.

- 1935 — Linus Pauling: first statistical-mechanical model of cooperativity via heme–heme interaction energy, microscopically grounding the Adair constants.

- 1949 — Pauling, Itano, Singer, Wells: identify sickle-cell anemia as a "molecular disease" of hemoglobin.

- 1959–1960 — Max Perutz: solves the crystal structure of hemoglobin (Nobel 1962), revealing the deoxy (T) and oxy (R) quaternary structures that physically underlie cooperativity.

- 1965 — Monod, Wyman, Changeux: publish the concerted (MWC) allosteric model, deriving the occupancy distribution from a two-state T⇌R equilibrium.

- 1966 — Koshland, Némethy, Filmer: publish the sequential (KNF) induced-fit model, an alternative microscopic origin for the same Adair phenomenology.

- 1970s–1990s — Ackers, Gary Gill, et al.: develop linkage thermodynamics and precise resolution of microscopic versus macroscopic constants, establishing the identifiability limits and the role of the binding polynomial as partition function.

open-problem

Microscopic resolution of intermediate states. Even with high-precision oxygen-equilibrium data, only four macroscopic Adair constants are accessible, whereas the tetramer ( $\alpha_2\beta_2$ ) has up to ten distinct microscopic ligation species (counting $\alpha$ / $\beta$ asymmetry and dimer–dimer interface states). Disentangling which intermediate microstates dominate — and whether the T→R quaternary switch occurs predominantly after the first, second, or third ligand — remains experimentally contentious and model-dependent (MWC predicts a sharp switch governed by $L$ ; KNF a graded one). Related open issues: (i) the precise statistical-thermodynamic role of the $\alpha_1\beta_2$ dimer interface and tetramer–dimer dissociation coupling to ligation (the "Ackers tertiary two-state" debates); (ii) whether a single equilibrium binding polynomial suffices, or whether kinetic intermediates and conformational dynamics require a non-equilibrium occupancy treatment; (iii) extending the framework to combined O₂/CO₂/H⁺/2,3-BPG multi-ligand binding polynomials, where the partition function becomes multivariate and the occupancy "vector" becomes a high-dimensional tensor whose identifiability is still incompletely characterized.

Hemoglobin O₂ Binding: Multiple Occupancy Probabilities

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also