Hill equation (biochemistry)

Last updated
Binding curves showing the characteristically sigmoidal curves generated by using the Hill equation to model cooperative binding. Each curve corresponds to a different Hill coefficient, labeled to the curve's right. The vertical axis displays the proportion of the total number of receptors that have been bound by a ligand. The horizontal axis is the concentration of the ligand. As the Hill coefficient is increased, the saturation curve becomes steeper. Hill-Langmuir equation.svg
Binding curves showing the characteristically sigmoidal curves generated by using the Hill equation to model cooperative binding. Each curve corresponds to a different Hill coefficient, labeled to the curve's right. The vertical axis displays the proportion of the total number of receptors that have been bound by a ligand. The horizontal axis is the concentration of the ligand. As the Hill coefficient is increased, the saturation curve becomes steeper.

In biochemistry and pharmacology, the Hill equation refers to two closely related equations that reflect the binding of ligands to macromolecules, as a function of the ligand concentration. A ligand is "a substance that forms a complex with a biomolecule to serve a biological purpose" (ligand definition), and a macromolecule is a very large molecule, such as a protein, with a complex structure of components (macromolecule definition). Protein-ligand binding typically changes the structure of the target protein, thereby changing its function in a cell.

Contents

The distinction between the two Hill equations is whether they measure occupancy or response. The Hill equation reflects the occupancy of macromolecules: the fraction that is saturated or bound by the ligand. [1] [2] [nb 1] This equation is formally equivalent to the Langmuir isotherm. [3] Conversely, the Hill equation proper reflects the cellular or tissue response to the ligand: the physiological output of the system, such as muscle contraction.

The Hill equation was originally formulated by Archibald Hill in 1910 to describe the sigmoidal O2 binding curve of haemoglobin. [4]

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 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. [5]

The Hill equation (for response) is important in the construction of dose-response curves.

Proportion of ligand-bound receptors

Plot of the % saturation of oxygen binding to haemoglobin, as a function of the amount of oxygen present (expressed as an oxygen pressure). Data (red circles) and Hill equation fit (black curve) from original 1910 paper of Hill. Plot of %25 saturation of O2 binding to haemoglobin as a function of O2 pressure - from 1910 Hill paper.png
Plot of the % saturation of oxygen binding to haemoglobin, as a function of the amount of oxygen present (expressed as an oxygen pressure). Data (red circles) and Hill equation fit (black curve) from original 1910 paper of Hill.

The Hill equation is commonly expressed in the following ways. [2] [7] [8]

,

where:

The special case where is a Monod equation.

Constants

In pharmacology, is often written as , where is the ligand, equivalent to L, and is the receptor. can be expressed in terms of the total amount of receptor and ligand-bound receptor concentrations: . is equal to the ratio of the dissociation rate of the ligand-receptor complex to its association rate (). [8] Kd is the equilibrium constant for dissociation. is defined so that , this is also known as the microscopic dissociation constant and is the ligand concentration occupying half of the binding sites. In recent literature, this constant is sometimes referred to as . [8]

Gaddum equation

The Gaddum equation is a further generalisation of the Hill-equation, incorporating the presence of a reversible competitive antagonist. [1] The Gaddum equation is derived similarly to the Hill-equation but with 2 equilibria: both the ligand with the receptor and the antagonist with the receptor. Hence, the Gaddum equation has 2 constants: the equilibrium constants of the ligand and that of the antagonist

Hill plot

A Hill plot, where the x-axis is the logarithm of the ligand concentration and the y-axis is the transformed receptor occupancy. X represents L and Y represents theta. Hill Plot.png
A Hill plot, where the x-axis is the logarithm of the ligand concentration and the y-axis is the transformed receptor occupancy. X represents L and Y represents theta.

The Hill plot is the rearrangement of the Hill equation into a straight line.

Taking the reciprocal of both sides of the Hill equation, rearranging, and inverting again yields: . Taking the logarithm of both sides of the equation leads to an alternative formulation of the Hill-Langmuir equation:

.

This last form of the Hill equation is advantageous because a plot of versus yields a linear plot, which is called a Hill plot. [7] [8] Because the slope of a Hill plot is equal to the Hill coefficient for the biochemical interaction, the slope is denoted by . A slope greater than one thus indicates positively cooperative binding between the receptor and the ligand, while a slope less than one indicates negatively cooperative binding.

Transformations of equations into linear forms such as this were very useful before the widespread use of computers, as they allowed researchers to determine parameters by fitting lines to data. However, these transformations affect error propagation, and this may result in undue weight to error in data points near 0 or 1. [nb 2] This impacts the parameters of linear regression lines fitted to the data. Furthermore, the use of computers enables more robust analysis involving nonlinear regression.

Tissue response

A trio of dose response curves Dose response antagonist.jpg
A trio of dose response curves

A distinction should be made between quantification of drugs binding to receptors and drugs producing responses. There may not necessarily be a linear relationship between the two values. In contrast to this article's previous definition of the Hill equation, the IUPHAR defines the Hill equation in terms of the tissue response , as [1]

where is the drug concentration, is the Hill coefficient, and is the drug concentration that produces a 50% maximal response. Dissociation constants (in the previous section) relate to ligand binding, while reflects tissue response.

This form of the equation can reflect tissue/cell/population responses to drugs and can be used to generate dose response curves. The relationship between and EC50 may be quite complex as a biological response will be the sum of myriad factors; a drug will have a different biological effect if more receptors are present, regardless of its affinity.

The Del-Castillo Katz model is used to relate the Hill equation to receptor activation by including a second equilibrium of the ligand-bound receptor to an activated form of the ligand-bound receptor.

Statistical analysis of response as a function of stimulus may be performed by regression methods such as the probit model or logit model, or other methods such as the Spearman–Kärber method. [9] Empirical models based on nonlinear regression are usually preferred over the use of some transformation of the data that linearizes the dose-response relationship. [10]

Hill coefficient

The Hill coefficient is a measure of ultrasensitivity (i.e. how steep is the response curve).

The Hill coefficient, or , may describe cooperativity (or possibly other biochemical properties, depending on the context in which the Hill equation is being used). When appropriate,[ clarification needed ] the value of the Hill coefficient describes the cooperativity of ligand binding in the following way:

The Hill coefficient can be calculated approximately in terms of the cooperativity index of Taketa and Pogell [12] as follows: [13]

.

where and are the input values needed to produce the 10% and 90% of the maximal response, respectively.


Reversible form

The most common form of the Hill equation is its irreversible form. However, when building computational models a reversible form is often required in order to model product inhibition. For this reason, Hofmeyr and Cornish-Bowden devised the reversible Hill equation. [14]

Relationship to the elasticity coefficients

The Hill coefficient is also intimately connected to the elasticity coefficient where the Hill coefficient can be shown to equal:

where is the fractional saturation, , and the elasticity coefficient.

This is derived by taking the slope of the Hill equation:

and expanding the slope using the quotient rule. The result shows that the elasticity can never exceed since the equation above can be rearranged to:

Applications

The Hill equation is used extensively in pharmacology to quantify the functional parameters of a drug[ citation needed ] and are also used in other areas of biochemistry.

The Hill equation can be used to describe dose-response relationships, for example ion channel open-probability (P-open) vs. ligand concentration. [15]

Regulation of gene transcription

The Hill equation can be applied in modelling the rate at which a gene product is produced when its parent gene is being regulated by transcription factors (e.g., activators and/or repressors). [11] Doing so is appropriate when a gene is regulated by multiple binding sites for transcription factors, in which case the transcription factors may bind the DNA in a cooperative fashion. [16]

If the production of protein from gene X is up-regulated (activated) by a transcription factor Y, then the rate of production of protein X can be modeled as a differential equation in terms of the concentration of activated Y protein:

,

where k is the maximal transcription rate of gene X.

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

,

where k is the maximal transcription rate of gene Y.

Limitations

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

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. [2] Nevertheless, other, more complex models of cooperative binding have been proposed. [7] For more information and examples of such models, see Cooperative binding.

Global sensitivity measure such as Hill coefficient do not characterise the local behaviours of the s-shaped curves. Instead, these features are well captured by the response coefficient measure. [18]

There is a link between Hill Coefficient and Response coefficient, as follows. Altszyler et al. (2017) have shown that these ultrasensitivity measures can be linked. [13]

See also

Notes

  1. For clarity, this article will use the International Union of Basic and Clinical Pharmacology convention of distinguishing between the Hill-Langmuir equation (for receptor saturation) and Hill equation (for tissue response)
  2. See Propagation of uncertainty. The function propagates errors in as . Hence errors in values of near or are given far more weight than those for

Related Research Articles

In chemistry, biochemistry, and pharmacology, a dissociation constant (KD) is a specific type of equilibrium constant that measures the propensity of a larger object to separate (dissociate) reversibly into smaller components, as when a complex falls apart into its component molecules, or when a salt splits up into its component ions. The dissociation constant is the inverse of the association constant. In the special case of salts, the dissociation constant can also be called an ionization constant. For a general reaction:

Cooperative binding occurs in molecular binding systems containing more than one type, or species, of molecule and in which one of the partners is not mono-valent and can bind more than one molecule of the other species. In general, molecular binding is an interaction between molecules that results in a stable physical association between those molecules.

Cooperativity is a phenomenon displayed by systems involving identical or near-identical elements, which act dependently of each other, relative to a hypothetical standard non-interacting system in which the individual elements are acting independently. One manifestation of this is enzymes or receptors that have multiple binding sites where the affinity of the binding sites for a ligand is apparently increased, positive cooperativity, or decreased, negative cooperativity, upon the binding of a ligand to a binding site. For example, when an oxygen atom binds to one of hemoglobin's four binding sites, the affinity to oxygen of the three remaining available binding sites increases; i.e. oxygen is more likely to bind to a hemoglobin bound to one oxygen than to an unbound hemoglobin. This is referred to as cooperative binding.

<span class="mw-page-title-main">Receptor (biochemistry)</span> Protein molecule receiving signals for a cell

In biochemistry and pharmacology, receptors are chemical structures, composed of protein, that receive and transduce signals that may be integrated into biological systems. These signals are typically chemical messengers which bind to a receptor and produce physiological responses such as change in the electrical activity of a cell. For example, GABA, an inhibitory neurotransmitter, inhibits electrical activity of neurons by binding to GABAA receptors. There are three main ways the action of the receptor can be classified: relay of signal, amplification, or integration. Relaying sends the signal onward, amplification increases the effect of a single ligand, and integration allows the signal to be incorporated into another biochemical pathway.

<span class="mw-page-title-main">Pharmacodynamics</span> Area of Academic Study

Pharmacodynamics (PD) is the study of the biochemical and physiologic effects of drugs. The effects can include those manifested within animals, microorganisms, or combinations of organisms.

IC<sub>50</sub> Half maximal inhibitory concentration

Half maximal inhibitory concentration (IC50) is a measure of the potency of a substance in inhibiting a specific biological or biochemical function. IC50 is a quantitative measure that indicates how much of a particular inhibitory substance (e.g. drug) is needed to inhibit, in vitro, a given biological process or biological component by 50%. The biological component could be an enzyme, cell, cell receptor or microorganism. IC50 values are typically expressed as molar concentration.

<span class="mw-page-title-main">Monod–Wyman–Changeux model</span> Biochemical model of protein transitions

In biochemistry, the Monod–Wyman–Changeux model describes allosteric transitions of proteins made up of identical subunits. It was proposed by Jean-Pierre Changeux in his PhD thesis, and described by Jacques Monod, Jeffries Wyman, and Jean-Pierre Changeux. It contrasts with the sequential model.

In statistics, Poisson regression is a generalized linear model form of regression analysis used to model count data and contingency tables. Poisson regression assumes the response variable Y has a Poisson distribution, and assumes the logarithm of its expected value can be modeled by a linear combination of unknown parameters. A Poisson regression model is sometimes known as a log-linear model, especially when used to model contingency tables.

The Scatchard equation is an equation used in molecular biology to calculate the affinity and number of binding sites of a receptor for a ligand. It is named after the American chemist George Scatchard.

<span class="mw-page-title-main">Enzyme kinetics</span> Study of biochemical reaction rates catalysed by an enzyme

Enzyme kinetics is the study of the rates of enzyme-catalysed chemical reactions. In enzyme kinetics, the reaction rate is measured and the effects of varying the conditions of the reaction are investigated. Studying an enzyme's kinetics in this way can reveal the catalytic mechanism of this enzyme, its role in metabolism, how its activity is controlled, and how a drug or a modifier might affect the rate.

EC<sub>50</sub> Concentration of a compound where 50% of its maximal effect is observed

Half maximal effective concentration (EC50) is a measure of the concentration of a drug, antibody or toxicant which induces a biological response halfway between the baseline and maximum after a specified exposure time. More simply, EC50 can be defined as the concentration required to obtain a 50% [...] effect and may be also written as [A]50. It is commonly used as a measure of a drug's potency, although the use of EC50 is preferred over that of 'potency', which has been criticised for its vagueness. EC50 is a measure of concentration, expressed in molar units (M), where 1 M is equivalent to 1 mol/L.

<span class="mw-page-title-main">Dose–response relationship</span> Measure of organism response to stimulus

The dose–response relationship, or exposure–response relationship, describes the magnitude of the response of an organism, as a function of exposure to a stimulus or stressor after a certain exposure time. Dose–response relationships can be described by dose–response curves. This is explained further in the following sections. A stimulus response function or stimulus response curve is defined more broadly as the response from any type of stimulus, not limited to chemicals.

In biochemistry, receptor–ligand kinetics is a branch of chemical kinetics in which the kinetic species are defined by different non-covalent bindings and/or conformations of the molecules involved, which are denoted as receptor(s) and ligand(s). Receptor–ligand binding kinetics also involves the on- and off-rates of binding.

<span class="mw-page-title-main">Schild equation</span>

In pharmacology, Schild regression analysis, based upon the Schild equation, both named for Heinz Otto Schild, are tools for studying the effects of agonists and antagonists on the response caused by the receptor or on ligand-receptor binding.

<span class="mw-page-title-main">Rotational diffusion</span>

Rotational diffusion is the rotational movement which acts upon any object such as particles, molecules, atoms when present in a fluid, by random changes in their orientations. Whilst the directions and intensities of these changes are statistically random, they do not arise randomly and are instead the result of interactions between particles. One example occurs in colloids, where relatively large insoluble particles are suspended in a greater amount of fluid. The changes in orientation occur from collisions between the particle and the many molecules forming the fluid surrounding the particle, which each transfer kinetic energy to the particle, and as such can be considered random due to the varied speeds and amounts of fluid molecules incident on each individual particle at any given time.

Receptor theory is the application of receptor models to explain drug behavior. Pharmacological receptor models preceded accurate knowledge of receptors by many years. John Newport Langley and Paul Ehrlich introduced the concept that receptors can mediate drug action at the beginning of the 20th century. Alfred Joseph Clark was the first to quantify drug-induced biological responses. So far, nearly all of the quantitative theoretical modelling of receptor function has centred on ligand-gated ion channels and G protein-coupled receptors.

<span class="mw-page-title-main">Langmuir adsorption model</span> Model describing the adsorption of a mono-layer of gas molecules on an ideal flat surface

The Langmuir adsorption model explains adsorption by assuming an adsorbate behaves as an ideal gas at isothermal conditions. According to the model, adsorption and desorption are reversible processes. This model even explains the effect of pressure i.e. at these conditions the adsorbate's partial pressure, , is related to the volume of it, V, adsorbed onto a solid adsorbent. The adsorbent, as indicated in the figure, is assumed to be an ideal solid surface composed of a series of distinct sites capable of binding the adsorbate. The adsorbate binding is treated as a chemical reaction between the adsorbate gaseous molecule and an empty sorption site, S. This reaction yields an adsorbed species with an associated equilibrium constant :

In chemistry, binding selectivity is defined with respect to the binding of ligands to a substrate forming a complex. Binding selectivity describes how a ligand may bind more preferentially to one receptor than another. A selectivity coefficient is the equilibrium constant for the reaction of displacement by one ligand of another ligand in a complex with the substrate. Binding selectivity is of major importance in biochemistry and in chemical separation processes.

<span class="mw-page-title-main">Ultrasensitivity</span>

In molecular biology, ultrasensitivity describes an output response that is more sensitive to stimulus change than the hyperbolic Michaelis-Menten response. Ultrasensitivity is one of the biochemical switches in the cell cycle and has been implicated in a number of important cellular events, including exiting G2 cell cycle arrests in Xenopus laevis oocytes, a stage to which the cell or organism would not want to return.

<span class="mw-page-title-main">Competitive inhibition</span> Interruption of a chemical pathway

Competitive inhibition is interruption of a chemical pathway owing to one chemical substance inhibiting the effect of another by competing with it for binding or bonding. Any metabolic or chemical messenger system can potentially be affected by this principle, but several classes of competitive inhibition are especially important in biochemistry and medicine, including the competitive form of enzyme inhibition, the competitive form of receptor antagonism, the competitive form of antimetabolite activity, and the competitive form of poisoning.

References

  1. 1 2 3 Neubig, Richard R. (2003). "International Union of Pharmacology Committee on Receptor Nomenclature and Drug Classification. XXXVIII. Update on Terms and Symbols in Quantitative Pharmacology" (PDF). Pharmacological Reviews. 55 (4): 597–606. doi:10.1124/pr.55.4.4. PMID   14657418. S2CID   1729572.
  2. 1 2 3 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. S2CID   122929930.
  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. Hill, A. V. (1910-01-22). "The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves". J. Physiol. 40 (Suppl): iv–vii. doi:10.1113/jphysiol.1910.sp001386. S2CID   222195613.
  5. 1 2 3 4 5 Weiss, J. N. (1 September 1997). "The Hill equation revisited: uses and misuses". The FASEB Journal. 11 (11): 835–841. doi: 10.1096/fasebj.11.11.9285481 . ISSN   0892-6638. PMID   9285481. S2CID   827335.
  6. "Proceedings of the Physiological Society: January 22, 1910". The Journal of Physiology. 40 (suppl): i–vii. 1910. doi:10.1113/jphysiol.1910.sp001386. ISSN   1469-7793. S2CID   222195613.
  7. 1 2 3 Stefan, Melanie I.; Novère, Nicolas Le (27 June 2013). "Cooperative Binding". PLOS Computational Biology. 9 (6): e1003106. Bibcode:2013PLSCB...9E3106S. doi: 10.1371/journal.pcbi.1003106 . ISSN   1553-7358. PMC   3699289 . PMID   23843752.
  8. 1 2 3 4 Nelson, David L.; Cox, Michael M. (2013). Lehninger principles of biochemistry (6th ed.). New York: W.H. Freeman. pp. 158–162. ISBN   978-1429234146.
  9. Hamilton, MA; Russo, RC; Thurston, RV (1977). "Trimmed Spearman–Karber method for estimating median lethal concentrations in toxicity bioassays". Environmental Science & Technology . 11 (7): 714–9. Bibcode:1977EnST...11..714H. doi:10.1021/es60130a004.
  10. Bates, Douglas M.; Watts, Donald G. (1988). Nonlinear Regression Analysis and its Applications . Wiley. p.  365. ISBN   9780471816430.
  11. 1 2 Alon, Uri (2007). An Introduction to Systems Biology: Design Principles of Biological Circuits ([Nachdr.] ed.). Boca Raton, FL: Chapman & Hall. ISBN   978-1-58488-642-6.
  12. Taketa, K.; Pogell, B. M. (1965). "Allosteric inhibition of rat-liver fructose 1,6-diphosphatase by adenosine 5'-monophosphate". J. Biol. Chem. 240 (2): 651–662. doi: 10.1016/S0021-9258(17)45224-0 . PMID   14275118.
  13. 1 2 Altszyler, E; Ventura, A. C.; Colman-Lerner, A.; Chernomoretz, A. (2017). "Ultrasensitivity in signaling cascades revisited: Linking local and global ultrasensitivity estimations". PLOS ONE. 12 (6): e0180083. arXiv: 1608.08007 . Bibcode:2017PLoSO..1280083A. doi: 10.1371/journal.pone.0180083 . PMC   5491127 . PMID   28662096.
  14. Hofmeyr, Jan-Hendrik S.; Cornish-Bowden, Athel (1997). "The reversible Hill equation: how to incorporate cooperative enzymes into metabolic models". Bioinformatics. 13 (4): 377–385. doi: 10.1093/bioinformatics/13.4.377 . PMID   9283752.
  15. Ding, S; Sachs, F (1999). "Single Channel Properties of P2X2 Purinoceptors". J. Gen. Physiol. 113 (5). The Rockefeller University Press: 695–720. doi:10.1085/jgp.113.5.695. PMC   2222910 . PMID   10228183.
  16. Chu, Dominique; Zabet, Nicolae Radu; Mitavskiy, Boris (2009-04-07). "Models of transcription factor binding: Sensitivity of activation functions to model assumptions" (PDF). Journal of Theoretical Biology. 257 (3): 419–429. Bibcode:2009JThBi.257..419C. doi:10.1016/j.jtbi.2008.11.026. PMID   19121637.
  17. Monod, Jacques; 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. PMID   14343300.
  18. Kholodenko, Boris N.; et al. (1997). "Quantification of information transfer via cellular signal transduction pathways". FEBS Letters. 414 (2): 430–434. doi: 10.1016/S0014-5793(97)01018-1 . PMID   9315734. S2CID   19466336.

Further reading