The Hill equation is a mathematical formula that describes how molecules bind to proteins in a cooperative way, where the binding of one molecule makes it easier (or harder) for the next one to bind. It was originally developed by Archibald Hill in 1910 to describe how oxygen binds to hemoglobin, and it remains one of the most widely used equations in biochemistry, pharmacology, and modern systems biology.
At its core, the equation captures something intuitive: some biological systems don’t respond gradually to increasing concentrations of a substance. Instead, they stay quiet for a while, then switch on rapidly over a narrow concentration range. The Hill equation puts a number on how sharp that switch is.
The Formula and Its Variables
The standard form of the Hill equation is:
Θ = [L]ⁿ / (K_D ⁿ + [L]ⁿ)
Each variable represents something concrete:
- Θ (theta) is the fractional response, ranging from 0 to 1. It could represent the fraction of binding sites occupied on a protein, or the fraction of maximum biological response achieved. Think of it as “how much of the total possible effect is happening right now.”
- [L] is the concentration of the ligand, meaning whatever molecule is doing the binding (oxygen, a drug, a signaling molecule).
- K_D is the dissociation constant, the ligand concentration at which half the binding sites are occupied. In pharmacology, this same value is called the EC₅₀ (the dose producing 50% of the maximum effect) or IC₅₀ (the dose producing 50% inhibition).
- n is the Hill coefficient (often written n_H), the number that captures cooperativity and controls the steepness of the response curve.
When n equals 1, the equation simplifies to the same form as the Michaelis-Menten equation, which describes simple, non-cooperative binding. The response curve in that case is a gentle hyperbola. As n increases above 1, the curve becomes sigmoidal (S-shaped), meaning the system resists responding at low concentrations but then switches on sharply once a threshold is crossed.
What the Hill Coefficient Tells You
The Hill coefficient is the most informative part of the equation. It acts as a measure of cooperativity, telling you how the binding of one molecule influences the binding of the next.
- n = 1: No cooperativity. Each binding event is independent. The system responds proportionally to increasing ligand concentration.
- n > 1: Positive cooperativity. Binding one molecule makes additional binding easier. The response curve becomes steeper and more switch-like. For a protein with n subunits, the theoretical maximum Hill coefficient equals n.
- n < 1: Negative cooperativity. Binding one molecule makes subsequent binding harder. The response curve is shallower than a simple hyperbola, meaning the system responds more gradually across a wider range of concentrations.
A real example of negative cooperativity: epidermal growth factor (EGF) binding to its receptor (EGFR) produces a Hill coefficient less than 1. This allows the receptor system to respond more gradually across a wider range of growth factor concentrations, rather than flipping on like a switch. Positive cooperativity does the opposite, compressing the responsive range and creating a sharper transition.
The Classic Example: Oxygen and Hemoglobin
The equation was born from the puzzle of hemoglobin. Each hemoglobin molecule has four binding sites for oxygen, and Hill noticed that the binding curve wasn’t a simple hyperbola. Instead, it was S-shaped: hemoglobin resists picking up the first oxygen molecule, but once one binds, the remaining sites fill much more easily.
Hill’s equation captured this behavior with a coefficient of roughly 2.8 for hemoglobin. This is less than the theoretical maximum of 4 (one for each binding site), which tells you the cooperativity is strong but not perfectly concerted. The practical consequence is profound: hemoglobin loads up almost completely with oxygen in the lungs, where oxygen is abundant, and then releases it efficiently in tissues where oxygen is scarce. That steep S-shaped curve makes hemoglobin a far better oxygen transporter than it would be with simple, non-cooperative binding.
Hill himself acknowledged the equation was an empirical approximation. It slightly underestimates saturation below about 30%, a region that rarely matters physiologically. More complex models, like Adair’s four-step binding equation, offer better accuracy across the full range, but Hill’s two-parameter formula remains dominant because of its simplicity.
How It’s Used in Pharmacology
In drug development, the Hill equation is the standard tool for analyzing dose-response curves. When you plot drug concentration on the x-axis and biological response on the y-axis, the data almost always fit a Hill-type curve. The two key parameters extracted from this fit are the EC₅₀ (or IC₅₀), which tells you how potent the drug is, and the Hill coefficient, which tells you how steeply the response changes around that midpoint.
A drug with a high Hill coefficient produces a narrow window between “almost no effect” and “near-maximum effect.” This can be desirable when you want a clean on/off response, but it also means small changes in dose can produce large swings in effect, making dosing less forgiving. A drug with a Hill coefficient near 1 produces a more gradual response across a wider dose range.
The equation also helps researchers determine whether drug combinations are additive, synergistic, or antagonistic. By comparing the observed dose-response curve of a combination against what the Hill equation predicts for each drug alone, you can quantify whether the drugs are enhancing or interfering with each other’s effects.
The Hill Plot: Measuring Cooperativity Experimentally
To determine the Hill coefficient from experimental data, researchers use a logarithmic transformation called a Hill plot. The equation rearranges to:
log[Θ / (1 – Θ)] = n_H × log[L] – n_H × log K_D
This has the form of a straight line (y = mx + b), where the slope equals the Hill coefficient. By plotting log[Θ / (1 – Θ)] against log[L], you get a line whose slope between 10% and 90% saturation gives you n_H. The x-intercept, where the line crosses zero, gives you log K_D.
In practice, real binding data don’t always produce a perfectly straight Hill plot. The slope often changes at the extremes of saturation, which is one reason the Hill equation is considered an approximation rather than a mechanistic model. It tells you cooperativity exists and how strong it is, but it doesn’t explain the underlying structural mechanism that produces it.
Beyond Binding: Gene Regulation and Biological Switches
The Hill equation has found a second life in systems biology, far from its origins in oxygen binding. When a transcription factor binds to a gene’s promoter region, the relationship between transcription factor concentration and gene output often follows a Hill function. In this context, the Hill coefficient measures “ultrasensitivity,” meaning how sharply a gene switches from off to on as its regulatory signal changes.
A Hill coefficient greater than 1 creates switch-like behavior. This is essential for biological circuits that need to make clean yes-or-no decisions: a cell differentiating into one type or another, a signaling pathway activating in response to a threshold stimulus, or an oscillator cycling between two states. Hill functions with high coefficients underlie the emergence of multiple steady states (where a cell commits to one of two fates) and limit cycle oscillations (where a system produces rhythmic behavior).
A 2024 paper in the Proceedings of the National Academy of Sciences showed that Hill functions represent a fundamental upper limit on the sharpness of input-output responses in biochemical networks. No matter how complex the underlying molecular machinery, the sharpness of the response can always be described and bounded by a Hill function. This makes the century-old equation not just a convenient approximation, but a universal reference point for understanding biological sensitivity.