How to compute the Hill equation: ligand-receptor cooperative binding model biochemistry pharmacology. Formula: Y = [L]^n / (K_d^n + [L]^n), where Y fractional saturation 0-1 (fraction binding sites occupied), [L] free ligand concentration (mol/L), K_d equilibrium dissociation constant (mol/L), n Hill coefficient cooperativity index (dimensionless). Y = 0.5 when [L] = K_d half-saturation point EC50. For n = 1 non-cooperative simple binding hyperbolic curve Michaelis-Menten enzyme single substrate (myoglobin O₂ binding); for n > 1 positive cooperativity sigmoidal S-curve (hemoglobin O₂ binding n≈2.8 quaternary structure tetramer α₂β₂ allosteric ligand binding facilitates more binding); for n < 1 negative cooperativity (insulin receptor binding initial binding inhibits further). History: Archibald Vivian Hill (UK 1886-1977) Cambridge University Trinity College physiologist 1910 paper "The Possible Effects of the Aggregation of the Molecules of Haemoglobin on its Dissociation Curves" Journal of Physiology 40:iv-vii (Proceedings) original derivation hemoglobin oxygen binding sigmoidal saturation curve, derived empirical equation modeling cooperative interaction multi-subunit protein.

Nobel Prize Physiology Medicine 1922 shared Otto Meyerhof (Germany 1884-1951) 'discoveries relating to the production of heat in muscle' studies muscle physiology lactic acid metabolism not directly cooperative binding but related muscle research. Hemoglobin Hb tetramer α₂β₂ four heme groups iron Fe2+ coordinate O₂ binding cooperative T-state R-state quaternary structure transition Monod-Wyman-Changeux MWC 1965 paper Journal Molecular Biology 12:88 + Koshland-Nemethy-Filmer KNF 1966 paper Biochemistry 5:365 sequential induced fit models. Essential tool for: biochemistry research USA NIH National Institutes Health Bethesda MD 27 institutes founded 1887 NIDDK National Institute Diabetes Digestive Kidney Diseases; Caltech Pasadena CA Linus Pauling Nobel 1954 1962; Harvard Medical School Boston MA founded 1782 oldest US medical school; Johns Hopkins Baltimore MD 1893 medical school; Stanford Medicine CA; UC San Francisco UCSF medical research; Yale School Medicine New Haven CT; UPenn School Medicine Philadelphia 1765 oldest US medical school; pharmacology drug receptor USA FDA Food Drug Administration Silver Spring MD 1906 18000 employees drug approval; ASPET American Society Pharmacology Experimental Therapeutics Bethesda MD 1908 5000 members; USA pharma industry: Pfizer NYC 1849 88000 employees; Merck Rahway NJ 1891 71000 employees; Johnson & Johnson New Brunswick NJ 1886 130000 employees;

Eli Lilly Indianapolis IN 1876 41000 employees diabetes drugs Trulicity Mounjaro; Bristol Myers Squibb NYC 1989 30000 employees; Regeneron Tarrytown NY 1988 11000 employees; Moderna Cambridge MA 2010 mRNA technology COVID-19 vaccine 2020 first commercial mRNA medicine. USA Nobel Physiology Medicine recent: Drew Weissman Katalin Karikó University Pennsylvania Nobel 2023 mRNA modification COVID vaccine technology; James Allison MD Anderson Tasuku Honjo Kyoto Nobel 2018 immune checkpoint inhibitors cancer therapy; William Kaelin Harvard Peter Ratcliffe Oxford Gregg Semenza Johns Hopkins Nobel 2019 oxygen sensing HIF cells.

Hill equation formula

Fundamental equation: Y = [L]^n / (K_d^n + [L]^n). Y fractional saturation 0-1 (occupied sites / total sites). [L] free ligand concentration (M, μM, nM depending system). K_d dissociation constant binding affinity (lower K_d = higher affinity). n Hill coefficient cooperativity (Archibald Hill 1910 derived from data fitting). Equation reduces Michaelis-Menten if n=1: Y = [L]/(K_d + [L]).

Linearization Hill plot: log(Y/(1-Y)) = n × log[L] - n × log K_d. Plot log(Y/(1-Y)) vs log[L]: slope = n Hill coefficient, intercept = -n × log K_d. Linear regression extract n and K_d from saturation data. Validity range 0.1 < Y < 0.9 (saturation extremes degenerate logit transformation).

Cooperativity classification: n 1 positive cooperativity sigmoidal (allosteric multi-subunit hemoglobin O₂ n=2.8, ATP synthase, voltage-gated sodium channels). Theoretical maximum n = number of binding sites (n_max = 4 for hemoglobin tetramer 4 hemes).

Dynamic range concentration: from Y=0.1 to Y=0.9 ligand concentration changes by factor (81)^(1/n). For n=1 (non-cooperative): factor 81 (4 orders magnitude). For n=2.8 (hemoglobin): factor (81)^(1/2.8) = 4.86. For n=5: factor (81)^(1/5) = 2.41 (very steep transition narrow range).

Hemoglobin oxygen binding: classic example positive cooperativity n=2.8 P50=26 mmHg (half-saturation oxygen partial pressure). Lung alveoli pO2=100 mmHg arterial saturation 98% O2; tissues pO2=40 mmHg venous saturation 75%; exercise pO2=20 mmHg saturation 30% efficient O2 release. Allosteric effectors: H+ Bohr effect lower pH (lactic acid muscle exercise) shift right curve increase O2 release; CO2 carbamate formation; 2,3-BPG bisphosphoglycerate red blood cells decrease affinity favors O2 release tissues; CO carbon monoxide poisoning competitive inhibitor 200x affinity O2.

History and context

Archibald Vivian Hill (UK 1886-1977) Bristol Cambridge Trinity College student 1905 mathematics natural sciences, Cambridge fellow 1910 paper "The Possible Effects of the Aggregation of the Molecules of Haemoglobin on its Dissociation Curves" Journal of Physiology 40:iv-vii Proceedings supplement original derivation Hill equation hemoglobin oxygen sigmoidal binding curve. Nobel Prize Physiology Medicine 1922 shared Otto Meyerhof (Germany 1884-1951) Heidelberg Berlin 'discoveries relating to the production of heat in muscle' physiologist studies muscle lactic acid metabolism oxidative phosphorylation respiration. Hill mostly known muscle physiology heat production studies, but Hill equation enduring legacy biochemistry pharmacology applies cooperativity ligand binding generally. Hill student Charles Sherrington Cambridge Nobel 1932; Hill professor University College London UCL 1923-1935; National Physical Laboratory NPL Teddington UK 1935-1939; Royal Society Foreign Secretary 1945-1951; Hill participated WWI anti-aircraft research; WWII German refugee scientists relocation England Britain. Hemoglobin structure Max Perutz (Austria-UK 1914-2002) Cambridge MRC Laboratory Molecular Biology 1937-2002 X-ray crystallography hemoglobin Nobel Prize Chemistry 1962 shared John Kendrew (UK 1917-1997) myoglobin structure Nobel 1962. Linus Pauling (USA 1901-1994) Caltech Nobel Chemistry 1954 chemical bond molecular structure hemoglobin alpha helix.

Monod-Wyman-Changeux MWC 1965 paper Journal Molecular Biology 12:88 'On the Nature of Allosteric Transitions' concerted symmetric model T-state tense R-state relaxed equilibrium hemoglobin allosteric transition cooperative binding mechanism. Koshland-Nemethy-Filmer KNF 1966 paper Biochemistry 5:365 sequential induced fit alternative model. USA biochemistry research: NIH National Institutes Health Bethesda MD 27 institutes founded 1887 NIDDK National Institute Diabetes Digestive Kidney Diseases; NCI National Cancer Institute Bethesda MD; NIAID National Institute Allergy Infectious Diseases Anthony Fauci 1984-2022; NHLBI National Heart Lung Blood Institute Bethesda MD; NIMH National Institute Mental Health. USA universities biochemistry Nobel: Caltech Pasadena CA Linus Pauling Nobel 1954 1962 chemical bond molecular structure; Harvard Medical School Boston MA founded 1782 oldest US medical school; Johns Hopkins Baltimore MD 1893 first US research medical school; Stanford Medicine Palo Alto CA founded 1908; UC San Francisco UCSF medical research center 1864 founded; Yale School Medicine New Haven CT 1810; UPenn School Medicine Philadelphia 1765 oldest US medical school Benjamin Franklin founder; UCLA Geffen School Medicine LA; Mount Sinai School Medicine NYC 1968. USA pharma industry: Pfizer NYC 1849 founded 88000 employees Lipitor Viagra COVID vaccine; Merck Rahway NJ 1891 71000 employees diabetes Januvia HPV Gardasil;

Johnson & Johnson New Brunswick NJ 1886 130000 employees medical devices consumer pharmaceuticals; Eli Lilly Indianapolis IN 1876 41000 employees diabetes drugs Trulicity Mounjaro insulin Humalog; Bristol Myers Squibb NYC 1989 30000 employees; Regeneron Tarrytown NY 1988 11000 employees biotechnology; Moderna Cambridge MA 2010 mRNA technology COVID-19 vaccine 2020 first commercial mRNA medicine; AbbVie North Chicago IL 30000 employees Humira Skyrizi; Gilead Sciences Foster City CA 1987 antivirals; Amgen Thousand Oaks CA 1980 8000 employees biotechnology Epogen Procrit Neulasta. USA universities biochemistry Caltech: Linus Pauling Caltech Nobel Chemistry 1954 chemical bond molecular structure 'The Nature of the Chemical Bond' 1939 textbook foundational; Nobel Peace Prize 1962 anti-nuclear weapons activist double Nobel laureate only person two unshared Nobel prizes. USA hemoglobin sickle cell research: James Herrick (USA 1861-1954) Chicago 1910 first description sickle cell anemia; Linus Pauling 1949 first molecular disease description Hemoglobin S Glu6Val; Vernon Ingram 1957 amino acid substitution mechanism. USA Nobel Physiology Medicine recent: Drew Weissman Katalin Karikó University Pennsylvania Nobel 2023 mRNA modification; James Allison MD Anderson Tasuku Honjo Kyoto Nobel 2018 immune checkpoint inhibitors PD-1 CTLA-4 cancer therapy; William Kaelin Harvard Peter Ratcliffe Oxford Gregg Semenza Johns Hopkins Nobel 2019 oxygen sensing HIF cells;

Charles Rice Rockefeller Michael Houghton Alberta Harvey Alter NIH Nobel 2020 hepatitis C virus discovery. ASPET American Society Pharmacology Experimental Therapeutics Bethesda MD 1908 founded 5000 members. AAPS American Association Pharmaceutical Scientists Arlington VA 1986 founded 8000 members. ASBMB American Society Biochemistry Molecular Biology Rockville MD 1906 founded 12000 members. FDA Food Drug Administration Silver Spring MD 1906 founded 18000 employees drug approval CDER Center Drug Evaluation Research. Software pharmacology dose-response: GraphPad Prism San Diego CA dose-response curve fitting Hill equation 4-parameter logistic model; SigmaPlot Systat Software San Jose CA; OriginLab Northampton MA scientific graphing fitting; MATLAB Statistics Toolbox; R packages drc dose-response curves; Python scipy.optimize curve_fit; FDA Phoenix WinNonlin pharmacokinetic modeling; Monolix Lyon France population PK/PD.

Hill coefficient examples table

• Myoglobin Mb O2 binding: n=1.0 simple non-cooperative
• Hemoglobin Hb O2 binding: n=2.8 positive cooperativity tetramer 4 hemes
• Hemoglobin fetal HbF: n=2.5 lower cooperativity higher O2 affinity
• Hemoglobin sickle cell HbS: n=2.6 mutation Glu6Val polymerization
• Insulin receptor IR: n=0.5-0.7 negative cooperativity
• Acetylcholine receptor nicotinic: n=1.5-2 cooperative
• Glycogen phosphorylase: n=1.6-2 cooperative ATP allosteric
• Tetrameric calmodulin Ca2+: n=2 binding cooperative
• Hemoglobin Bohr effect: pH dependence O2 affinity
• P50 hemoglobin adult: 26 mmHg arterial
• P50 hemoglobin fetal HbF: 18 mmHg higher affinity
• Lung alveoli pO2: 100 mmHg saturation 98%
• Tissues venous pO2: 40 mmHg saturation 75%
• Exercise muscle pO2: 20 mmHg saturation 30%
• Carbon monoxide CO affinity: 200x hemoglobin O2
• Oxygen Hb molecules: 4 binding sites tetramer α2β2
• Drug ED50: dose 50% effect therapeutic dose
• Drug LD50: lethal dose 50% mortality
• Therapeutic index TI: LD50/ED50 safety margin
• Hill plot slope: n Hill coefficient extracted
• Linearization log(Y/(1-Y)): vs log[L] slope n
• Validity range Hill plot: 0.1<Y<0.9
• Dynamic range Y 10-90%: factor (81)^(1/n)
• n=1 dynamic range: 81 factor (4 orders magnitude)
• n=2.8 hemoglobin range: 4.86 factor
• n=5 ultracoop range: 2.41 factor
• Hill 1910 paper: Journal Physiology 40:iv-vii Proceedings
• Hill Nobel 1922: shared Meyerhof muscle heat production
• MWC 1965: Monod-Wyman-Changeux concerted symmetric
• KNF 1966: Koshland-Nemethy-Filmer sequential induced fit
• Linus Pauling Nobel 1954: Caltech chemical bond
• USA Drew Weissman Katalin Karikó 2023: Penn mRNA

Detailed numerical examples

Example 1 - Hemoglobin O2 binding lung alveoli. Hb tetramer α2β2 4 sites O2 cooperative binding. n=2.8, Kd=P50=26 mmHg, [L]=pO2 alveoli 100 mmHg arterial. Y = 100^2.8 / (26^2.8 + 100^2.8) = 5.0e5 / (1.16e4 + 5.0e5) = 0.977. Saturation 97.7% O2 efficient lung uptake. In vein pO2=40 mmHg: Y = 40^2.8 / (26^2.8 + 40^2.8) = 30000/(11600+30000) = 0.72 = 72%. Tissue O2 release 25% di 100% = 25 mL/dL blood arterial-venous difference. USA: cardiovascular research AHA Dallas TX 32000 employees American Heart Association.

Example 2 - Myoglobin O2 binding muscle simple. Mb monomer 1 site O2 simple binding. n=1, Kd=1 mmHg, [L]=3 mmHg. Y = 3 / (1 + 3) = 0.75 = 75% saturation. Hyperbolic curve no cooperativity. Mb O2 storage muscle high affinity O2 reservoir tissue extreme low pO2 conditions exercise diving mammals whales. USA: research Caltech Pasadena CA Linus Pauling Nobel 1954 1962 hemoglobin alpha helix structure.

Example 3 - Insulin receptor negative cooperativity. IR insulin receptor heterotetramer α2β2 transmembrane RTK receptor tyrosine kinase. n=0.5 negative cooperativity (binding insulin first inhibits second). Kd=1 nM, [L]=1 nM. Y = 1^0.5 / (1^0.5 + 1^0.5) = 1/(1+1) = 0.5 half saturation. n<1 negative cooperativity binding curve flatter than hyperbolic. USA: diabetes research ADA American Diabetes Association Arlington VA 23000 members 1940 founded; Eli Lilly Indianapolis IN 1876 insulin Humalog largest insulin manufacturer USA.

Example 4 - Drug receptor binding pharmacology. Drug Kd=50 nM, n=1.2 weak positive cooperativity, Bmax=100 receptors, [L]=100 nM. Y = 100^1.2 / (50^1.2 + 100^1.2) = 251/(123+251) = 0.67. Bound = 100 × 0.67 = 67 receptors occupied. EC50=50 nM half-effective concentration. Therapeutic index TI = LD50/ED50 typical drugs 5-100 safety margin. USA FDA Silver Spring MD drug approval CDER; ASPET Bethesda MD 1908 5000 members; AAPS Arlington VA 1986 8000 members.

Example 5 - Extreme cooperativity allosteric n=5. Voltage-gated sodium channel Nav1.5 cardiac, Hill coefficient n=5 extreme cooperativity 4 transmembrane domains S4 voltage sensors. Kd=5 (membrane potential mV threshold), [L]=5 mV. Y = 5^5/(5^5+5^5) = 0.5 half-activation Vh. Steep transition voltage 5 mV change activates channel from 10% to 90% open over Δ=2.41 mV (factor (81)^(1/5)=2.41). Cardiac action potential rapid depolarization lithium-ion channel block lidocaine arrhythmia treatment.

Limits and considerations

1. Empirical phenomenological Hill equation: not derived from physical first principles, MWC 1965 KNF 1966 detailed mechanistic models needed.
2. n max = number binding sites theoretical: hemoglobin n_max=4 but observed n=2.8 indicates partial cooperativity not perfect concerted.
3. Single ligand simple model: real systems multiple ligands competitive inhibitors allosteric effectors complex.
4. Equilibrium thermodynamic assumption: kinetic non-equilibrium dynamics important transient processes.
5. Linear Hill plot range 0.1<Y<0.9: outside range logit transformation breaks down extreme saturation.
6. Cooperativity sign reversal: some systems n>1 first phase, n<1 second phase, complex curves.
7. Temperature dependence: Hill coefficient n weakly temperature dependent, K_d strongly T dependent van t Hoff equation.
8. Ionic strength pH effects: Bohr effect hemoglobin allosteric H+ binding modify cooperativity.
9. Numerical precision: very high n (>10) extreme cooperativity edge case computational stability.
10. Physical constraints: n cannot exceed number sites (4 for tetramer), apparent n higher indicates model limitation.

Common errors

1. Confusing Hill equation with Michaelis-Menten: Hill general n cooperativity, MM specific n=1.
2. Neglecting units consistency [L] vs Kd: same units required (M vs mM vs μM).
3. Pretending n = number sites exact: theoretical max only, real data lower partial cooperativity.
4. Confusing Kd with EC50: in Hill model EC50 = Kd half-saturation, but for enzyme kinetics MM Km different.
5. Ignoring allosteric effectors: pH CO2 2,3-BPG modify hemoglobin n curve.
6. JavaScript Math.pow precision: very large or small bases L^n n high accumulate float64 errors.
7. Confusing positive vs negative cooperativity: n>1 positive (sigmoidal), n<1 negative (sub-hyperbolic).
8. Neglecting Bohr effect physiological: arterial-venous pH differ slightly modify hemoglobin curve.
9. Ignoring carbon monoxide CO competitive: 200x hemoglobin O2 affinity, dangerous low concentration.
10. Pretending Hill applicable enzyme kinetics: MM Henri-Michaelis-Menten 1913 different model substrate-product.

Numbered applications

1. USA NIH biochemistry: National Institutes Health Bethesda MD 27 institutes founded 1887; NIDDK National Institute Diabetes Digestive Kidney Diseases; NCI National Cancer Institute; NIAID Anthony Fauci 1984-2022.
2. USA pharma industry: Pfizer NYC 1849 88000 employees Lipitor Viagra COVID; Merck Rahway NJ 1891 71000 employees Januvia HPV Gardasil; Johnson & Johnson NJ 1886 130000 employees; Eli Lilly Indianapolis 1876 41000.
3. USA biotech mRNA: Moderna Cambridge MA 2010 mRNA technology COVID-19 vaccine 2020 first commercial mRNA medicine; Pfizer-BioNTech Mainz Germany Pennsylvania Drew Weissman Katalin Karikó Nobel 2023.
4. USA universities biochemistry Nobel: Caltech Pasadena CA Linus Pauling Nobel 1954 1962 hemoglobin alpha helix; Harvard Medical 1782 oldest; Johns Hopkins Baltimore MD 1893; Stanford Medicine 1908; UCSF medical 1864.
5. USA hemoglobin sickle cell research: James Herrick Chicago 1910 first description; Linus Pauling 1949 first molecular disease HbS Glu6Val; Vernon Ingram 1957 amino acid substitution mechanism.
6. USA Nobel Physiology Medicine recent: Drew Weissman Katalin Karikó Penn Nobel 2023 mRNA; James Allison MD Anderson Tasuku Honjo Kyoto Nobel 2018 immune checkpoint PD-1 CTLA-4; Kaelin Ratcliffe Semenza Nobel 2019 oxygen sensing HIF.
7. USA professional societies: ASPET American Society Pharmacology Bethesda MD 1908 5000 members; AAPS American Association Pharmaceutical Scientists Arlington VA 1986 8000; ASBMB American Society Biochemistry Molecular Biology Rockville MD 1906 12000.
8. USA FDA pharma regulator: Food Drug Administration Silver Spring MD 1906 18000 employees drug approval CDER Center Drug Evaluation Research clinical trials Phase I-IV.
9. MWC KNF allosteric models: Monod-Wyman-Changeux 1965 concerted symmetric T-state R-state; Koshland-Nemethy-Filmer 1966 sequential induced fit alternative; foundation modern allosteric protein theory.
10. Software pharmacology dose-response: GraphPad Prism San Diego CA Hill equation 4-parameter logistic; SigmaPlot Systat San Jose CA; OriginLab Northampton MA; R drc; Python scipy.optimize curve_fit; Phoenix WinNonlin pharmacokinetic.

For support

For technical questions: info@calcoi.com. Educational biochemistry pharmacology tool. Applications Hill AV 1910 'The Possible Effects of the Aggregation of the Molecules of Haemoglobin on its Dissociation Curves' Journal of Physiology 40:iv-vii, Monod J Wyman J Changeux JP 1965 'On the Nature of Allosteric Transitions: A Plausible Model' Journal Molecular Biology 12:88, Koshland DE Nemethy G Filmer D 1966 'Comparison of Experimental Binding Data and Theoretical Models in Proteins Containing Subunits' Biochemistry 5:365, Perutz MF 1970 'Stereochemistry of Cooperative Effects in Haemoglobin' Nature 228:726.

Updated July 2026