Statistical Thermodynamics: The Partition Function
How nature counts possibilities and decides probabilities
At the heart of statistical mechanics lies a deceptively simple idea: nature samples from all possible arrangements of a system, but not equally. Each arrangement (or "microstate") gets a weight based on its energy, and the partition function is simply the sum of all these weights.
In the simulation below, ligands float in solution near some receptors. Each receptor can be empty or occupied. Watch how changing the ligand concentration and binding energy shifts the balance between these two states. The partition function Z keeps track of the total "score" of all possibilities.
The Boltzmann Distribution
Statistical mechanics tells us that nature assigns a probability to each microstate based on its energy. States with lower energy are more probable, but how much more? The answer is the Boltzmann distribution:
Probability of microstate i
p(Ei) = e-Ei/kBT / Z
The key insight is the exponential dependence on energy. A state that is just 1 kBT lower in energy is e ≈ 2.7 times more likely. A state 10 kBT lower? Over 22,000 times more likely!
Building the Partition Function
Consider a single receptor that can either be empty or have a ligand bound. Following the lattice model from Physical Biology of the Cell, we assign weights to each state:
Weight(empty) = 1
Weight(bound) = (c/c0) × e-Δε/kBT
where c is concentration, c0 is a reference concentration, and Δε = εbound - εsolution
The partition function Z is the sum of all weights:
Z = 1 + (c/c0) × e-Δε/kBT
And the probability of the receptor being bound is simply the bound weight divided by Z:
pbound = (c/c0) × e-Δε/kBT / Z
This is the famous Langmuir adsorption isotherm (or Hill function with n=1)
Energy vs. Entropy: The Great Competition
Notice that pbound depends on two competing factors:
🔗
ENERGY favors binding
Lower ε → stronger attraction
⚖
🌊
ENTROPY favors freedom
More places to be in solution
At low concentrations, entropy wins—ligands prefer the vast ocean of solution. At high concentrations or with strong binding energy, the energetic payoff of binding overcomes the entropic cost of being "trapped" at the receptor.
The concentration where pbound = 0.5 is called the dissociation constant Kd. It occurs when the two terms in Z are equal, representing the perfect balance between energy and entropy.
Kd = c0 × e+Δε/kBT
The concentration at which half the receptors are occupied