Hill equation (biochemistry)

Biochemical binding curves showing the characteristically sigmoidal curves generated by using the Hill equation to model cooperative binding. Each curve depicts a different Hill coefficient: light blue corresponds to n=0.5, red to n=1, gray to n=2, yellow to n=3, and dark blue to n=4. The vertical axis displays the Fraction (Theta) of occupied ligand-binding sites on a protein receptor, equal to ratio of the concentration of Ligand-bound Protein ([PLn_Total]) to the total concentration of Protein receptor ([P_Total]). The horizontal axis is the ratio of the ligand concentration producing half occupation (K_A) to the total Ligand concentration ([L]).

In biochemistry and pharmacology, the binding of a ligand to a macromolecule is often enhanced if there are already other ligands present on the same macromolecule (this is known as cooperative binding). The Hill equation, which was originally formulated by Archibald Hill in 1910 to describe the sigmoidal O₂ binding curve of haemoglobin,^[1] is used to describe the fraction of a macromolecule saturated by ligand as a function of the ligand concentration. The equation is useful for determining the degree of cooperativity of the ligand(s) binding to the enzyme or receptor. The Hill coefficient provides a way to quantify the degree of interaction between ligand binding sites.^[2]

The Hill equation is formally equivalent to the Langmuir isotherm,^[3] and as a result, it is sometimes referred to as the Hill-Langmuir equation.^[4]

Hill Equation (Mathematical Model)

A Hill coefficient of 1 indicates completely independent binding, regardless of how many additional ligands are already bound. Numbers greater than one indicate positive cooperativity, while numbers less than one indicate negative cooperativity. For example, the Hill coefficient of oxygen binding to haemoglobin (an example of positive cooperativity) falls within the range of 1.7-3.2.^[2]

The Hill equation is commonly expressed^[4][5][6] in the following ways:

${\displaystyle \theta ={[{\ce {L}}]^{n} \over K_{d}+[{\ce {L}}]^{n}}={[{\ce {L}}]^{n} \over (K_{A})^{n}+[{\ce {L}}]^{n}}={1 \over \left({K_{A} \over [{\ce {L}}]}\right)^{n}+1}}$,

The equation's terms are defined as follows:

• $\theta$ - Fraction of the ligand-binding sites on the receptor protein which are occupied by the ligand.
• ${\displaystyle {\ce {[L]}}}$ - Free (unbound) ligand concentration
• $K_{d}$ - Apparent dissociation constant derived from the law of mass action (the equilibrium constant for dissociation), which is equal to the ratio of the dissociation rate of the ligand-receptor complex to its association rate (${\textstyle K_{\rm {d}}={k_{\rm {d}} \over k_{\rm {a}}}}$).^[6]
• $K_{A}$ - The ligand concentration producing half occupation (ligand concentration occupying half of the binding sites). Because ${\textstyle K_{A}}$ is defined so that ${\textstyle (K_{A})^{n}=K_{\rm {d}}={k_{\rm {d}} \over k_{\rm {a}}}}$, this is also known as the microscopic dissociation constant. In recent literature, this constant is sometimes referred to as ${\textstyle K_{D}}$.
• $n$ - The Hill coefficient, describing cooperativity (or possibly other biochemical properties, depending on the context in which the Hill equation is being used)

When appropriate, the value of the Hill coefficient describes the cooperativity of ligand binding in the following way:

• $n>1$ - Positively cooperative binding: Once one ligand molecule is bound to the enzyme, its affinity for other ligand molecules increases.
• $n<1$ - Negatively cooperative binding: Once one ligand molecule is bound to the enzyme, its affinity for other ligand molecules decreases.
• $n=1$ - Noncooperative binding: The affinity of the enzyme for a ligand molecule is not dependent on whether or not other ligand molecules are already bound. When n=1, we obtain a model that can be modeled by Michaelis–Menten kinetics,^[7] in which ${\textstyle K_{D}=K_{A}=K_{M}}$, the Michaelis-Menten constant.

Hill Plot

Taking the reciprocal of both sides of the Hill equation, rearranging, and inverting again yields: ${\displaystyle {\theta \over 1-\theta }={[{\ce {L}}]^{n} \over K_{d}}={[{\ce {L}}]^{n} \over (K_{A})^{n}}}$.

Taking the logarithm of both sides of the equation leads to an alternative formulation of the Hill equation:

${\displaystyle \log \left({\theta \over 1-\theta }\right)=n\log {[{\ce {L}}]}-\log {K_{d}}=n\log {[{\ce {L}}]}-n\log {K_{A}}}$.

This last form of the Hill equation is advantageous because a plot of ${\textstyle \log \left({\theta \over 1-\theta }\right)}$ versus ${\displaystyle \log {[{\ce {L}}]}}$ yields a linear plot, which is called a Hill plot.^[5][6] Because the slope of a Hill plot is equal to the Hill Coefficient for the biochemical interaction, the slope is denoted by $n_{H}$. A slope greater than one thus indicates positively cooperative binding between the receptor and the ligand, while a slope less than one indicates negatives cooperative binding.

Derivation from Mass Action Kinetics

Consider a protein (${\displaystyle P}$), such as Haemoglobin or a protein receptor, with ${\displaystyle {\mathit {n}}}$ binding sites for ligands (${\displaystyle L}$). The binding of the ligands to the protein can be represented by the chemical equilibrium expression:

${\displaystyle {\ce {{P}+{\mathit {n}}{L}<=>[k_{a}][k_{d}]{P}{L}_{\mathit {n}}}}}$,

where $k_a$ (forward rate, or the rate of association of the protein-ligand complex) and $k_{d}$ (reverse rate, or the complex's rate of dissociation) are the reaction rate constants for the association of the ligands to the protein and their dissociation from the protein, respectively.^[6] From the law of mass action, the apparent dissociation constant $K_{d}$, an equilibrium constant, is given by:

${\displaystyle K_{\rm {d}}={k_{\rm {d}} \over k_{\rm {a}}}={{[{\rm {P}}][{\rm {L}}]^{\mathit {n}}} \over [{\rm {PL_{\mathit {n}}}}]}}$.

At the same time, $\theta$, the fraction of the ligand binding sites on the receptor protein which are occupied by the ligand, is given by:

${\displaystyle \theta ={\mathrm {Occupied\ Binding\ Sites} \over \mathrm {Total\ Binding\ Sites} }={[{\rm {PL_{\mathit {n}}}}] \over {[{\rm {P}}]\ +\ [{\rm {PL_{\mathit {n}}}}]}}}$.

By using the expression obtained earlier for the dissociation constant, we can replace ${\textstyle [{\rm {PL_{\mathit {n}}}}]}$ with ${\textstyle {[{\rm {P}}][{\rm {L}}]^{\mathit {n}} \over K_{\rm {d}}}}$ to yield a simplified expression for ${\textstyle \theta }$:

${\displaystyle \theta ={({[{\rm {P}}][{\rm {L}}]^{\mathit {n}} \over K_{\rm {d}}}) \over {[{\rm {P}}]\ +\ ({[{\rm {P}}][{\rm {L}}]^{\mathit {n}} \over K_{\rm {d}}})}}={{[{\rm {P}}][{\rm {L}}]^{\mathit {n}}} \over {K_{\rm {d}}[{\rm {P}}]\ +\ {[{\rm {P}}][{\rm {L}}]^{\mathit {n}}}}}={{[{\rm {L}}]^{\mathit {n}}} \over {K_{\rm {d}}\ +\ {[{\rm {L}}]^{\mathit {n}}}}}}$,

which is a common formulation of the Hill equation.^[5][6]

Assuming that the protein receptor was initially completely free (unbound) at a concentration ${\textstyle [{\rm {P_{0}}}]}$, then at any time, ${\textstyle {[{\rm {P}}]+[{\rm {PL_{\mathit {n}}}}]}=[{\rm {P_{0}}}]}$ and ${\textstyle \theta ={[{\rm {PL_{\mathit {n}}}}] \over {[{\rm {P_{0}}}]\ }}}$. Consequently, the Hill Equation is also commonly written as an expression for the concentration ${\textstyle [{\rm {PL_{\mathit {n}}}}]}$of bound protein:

${\displaystyle [{\rm {PL_{\mathit {n}}}}]=[{\rm {P_{0}}}]\cdot {{[{\rm {L}}]^{\mathit {n}}} \over {K_{\rm {d}}\ +\ {[{\rm {L}}]^{\mathit {n}}}}}}$ .^[4]

Note that all of these formulations assume that the protein has ${\displaystyle {\mathit {n}}}$ sites to which ligands can bind. In practice, however, the Hill Coefficient ${\displaystyle {\mathit {n}}}$ rarely provides an accurate approximation of the number of ligand binding sites on a protein.^[2][5] Consequently, it has been observed that the Hill Coefficient should instead be interpreted as an "interaction coefficient" describing the cooperativity among ligand binding sites.^[2]

Applications to Mathematical Modeling

Regulation of Gene Transcription

The Hill equation can be applied in modeling the rate at which a gene product is produced when its parent gene is being regulated by transcription factors (activators and/or repressors).^[7] If the production of protein from gene ${\textstyle X}$ is up-regulated (activated) by a transcription factor ${\textstyle Y}$, then the rate of production of protein ${\textstyle X}$ can be modeled as a differential equation in terms of the concentration of activated ${\textstyle Y}$ protein:

${\displaystyle {d[X_{produced}] \over dt}=k\ \cdot {{[{\rm {Y_{active}}}]^{\mathit {n}}} \over {(K_{A})^{n}\ +\ {[{\rm {Y_{active}}}]^{\mathit {n}}}}}}$,

where ${\displaystyle k}$ is the maximal transcription rate of gene ${\textstyle X}$.

Likewise, if the production of protein from gene ${\textstyle Y}$ is down-regulated (repressed) by a transcription factor ${\textstyle Z}$, then the rate of production of protein ${\textstyle Y}$ can be modeled as a differential equation in terms of the concentration of activated ${\textstyle Z}$ protein:

${\displaystyle {d[Y_{produced}] \over dt}={{k} \over {1\ +\ ({[{\rm {Y_{active}}}] \over K_{A}})^{n}}}}$.

Limitations

Because of its assumption that ligand molecules bind to a receptor simultaneously, the Hill equation has been criticized as a physically unrealistic model.^[2] Moreover, the Hill coefficient should not be considered a reliable approximation of the number of cooperative ligand binding sites on a receptor^[2][8] except when the binding of the first and subsequent ligands results in extreme positive cooperativity.^[2]

Unlike more complex models, the relatively simple Hill equation provides little insight into underlying physiological mechanisms of protein-ligand interactions. This simplicity, however, is what makes the Hill equation a useful empirical model, since its use requires little a priori knowledge about the properties of either the protein or ligand being studied.^[4] Nevertheless, other, more complex models of cooperative binding have been proposed.^[5] For more information and examples of such models, see Cooperative binding.

See also

• Cooperative Binding
• Logistic function
• Gompertz curve
• Sigmoid function

References

1. Hill, A. V. (1910-01-22). "The possible effects of the aggregation of the molecules of hæmoglobin on its dissociation curves". J. Physiol. 40 (Suppl): iv–vii. doi:10.1113/jphysiol.1910.sp001386.
2. Weiss, J. N. (1 September 1997). "The Hill equation revisited: uses and misuses.". The FASEB Journal. 11 (11): 835–841. ISSN 0892-6638.
3. Langmuir, Irving (1918). "The adsorption of gases on plane surfaces of glass, mica and platinum.". Journal of the American Chemical Society. 40 (9): 1361–1403. doi:10.1021/ja02242a004.
4. Gesztelyi, Rudolf; Zsuga, Judit; Kemeny-Beke, Adam; Varga, Balazs; Juhasz, Bela; Tosaki, Arpad (31 March 2012). "The Hill equation and the origin of quantitative pharmacology". Archive for History of Exact Sciences. 66 (4): 427–438. doi:10.1007/s00407-012-0098-5. ISSN 0003-9519.
5. Stefan, Melanie I.; Novère, Nicolas Le (27 June 2013). "Cooperative Binding". PLOS Computational Biology. 9 (6): e1003106. doi:10.1371/journal.pcbi.1003106. ISSN 1553-7358.
6. Nelson, David L.; Cox, Michael M. (2013). Lehninger principles of biochemistry (6th ed. ed.). New York: W.H. Freeman. pp. 158–162. ISBN 978-1429234146.
7. Alon, Uri (2007). An Introduction to Systems Biology: Design Principles of Biological Circuits ([Nachdr.] ed.). Boca Raton, FL: Chapman & Hall. ISBN 1-58488-642-0.
8. Monod, Jacque; Wyman, Jeffries; Changeux, Jean-Pierre (1 May 1965). "On the nature of allosteric transitions: A plausible model". Journal of Molecular Biology. 12 (1): 88–118. doi:10.1016/S0022-2836(65)80285-6. Retrieved 5 November 2016.

Further reading

Endrenyi

• Dorland's Illustrated Medical Dictionary
• Nelson, David L.; Cox, Michael M. Lehninger Principles of Biochemistry (4th ed.).
• Coval, ML (December 1970). "Analysis of Hill interaction coefficients and the invalidity of the Kwon and Brown equation". J. Biol. Chem. 245 (23): 6335–6. PMID 5484812.
• d'A Heck, Henry (1971). "Statistical theory of cooperative binding to proteins. Hill equation and the binding potential". J. Am. Chem. Soc. 93 (1): 23–29. doi:10.1021/ja00730a004.
• Atkins, Gordon L. (1973). "A simple digital-computer program for estimating the parameter of the Hill Equation". Eur. J. Biochem. 33 (1): 175–180. doi:10.1111/j.1432-1033.1973.tb02667.x.
• Endrenyi, Laszlo; Kwong, F. H. F.; Fajszi, Csaba (1975). "Evaluation of Hill slopes and Hill coefficients when the saturation binding or velocity is not known". Eur. J. Biochem. 51 (2): 317–328. doi:10.1111/j.1432-1033.1975.tb03931.x.
• Voet, Donald; Voet, Judith G. Biochemistry.
• Weiss, J. N. (1997). "The Hill equation revisited: uses and misuses". FASEB Journal. 11 (11): 835–841. PMID 9285481.
• Kurganov, B. I.; Lobanov, A. V. (2001). "Criterion for Hill equation validity for description of biosensor calibration curves". Anal. Chim. Acta. 427 (1): 11–19. doi:10.1016/S0003-2670(00)01167-3.
• Goutelle, Sylvain; Maurin, Michel; Rougier, Florent; Barbaut, Xavier; Bourguignon, Laurent; Ducher, Michel; Maire, Pascal (2008). "The Hill equation: a review of its capabilities in pharmacological modelling". Fund. Clinic. Pharmac. 22: 633–648. doi:10.1111/j.1472-8206.2008.00633.x.
• Mitavskiy, B.; Chu, D.; Zabet., Radu (2009). "Models of transcription factor binding: Sensitivity of activation functions to model assumptions.". Journal of Theoretical Biology. 257 (3): 419–429. doi:10.1016/j.jtbi.2008.11.026. PMID 19121637.
• Gesztelyi R; Zsuga J; Kemeny-Beke A; Varga B; Juhasz B; Tosaki A (2012). "The Hill equation and the origin of quantitative pharmacology". Archive for History of Exact Sciences. 66: 427–38. doi:10.1007/s00407-012-0098-5.
• Colquhoun D (2006). "The quantitative analysis of drug-receptor interactions: a short history". Trends Pharmacol Sci. 27 (3): 149–57. doi:10.1016/j.tips.2006.01.008. PMID 16483674.
• Rang HP (2006). "The receptor concept: pharmacology's big idea". Br J Pharmacol. 147: S9–16. doi:10.1038/sj.bjp.0706457. PMC 1760743. PMID 16402126.
• Neubig RR; Spedding M; Kenakin T; Christopoulos A (2003). "International Union of Pharmacology Committee on Receptor Nomenclature and Drug Classification. XXXVIII. update on terms and symbols in quantitative pharmacology". Pharmacol Rev. 55: 597–606. doi:10.1124/pr.55.4.4.

External links

• Hill equation calculator