Chemistry:Reversible Hill equation

From HandWiki

The classic Monod–Wyman–Changeux model (MWC) for cooperativity is generally published in an irreversible form. That is, there are no product terms in the rate equation which can be problematic for those wishing to build metabolic models since there are no product inhibition terms.[1] However, a series of publications by Popova and Sel'kov[2] derived the MWC rate equation for the reversible, multi-substrate, multi-product reaction.

The same problem applies to the classic Hill equation which is almost always shown in an irreversible form. Hofmeyr and Cornish-Bowden first published the reversible form of the Hill equation.[1] The equation has since been discussed elsewhere[3][4] and the model has also been used in a number of kinetic models such as a model of Phosphofructokinase and Glycolytic Oscillations in the Pancreatic β-cells[5] or a model of a glucose-xylose co-utilizing S. cerevisiae strain.[6] The model has also been discussed in modern enzyme kinetics textbooks.[7][8]

Derivation

Consider the simpler case where there are two binding sites. See the scheme shown below. Each site is assumed to bind either molecule of substrate S or product P. The catalytic reaction is shown by the two reactions at the base of the scheme triangle, that is S to P and P to S. The model assumes the binding steps are always at equilibrium. The reaction rate is given by:

[math]\displaystyle{ v=k_1\left(E S+2 E S_2+E S P\right)-k_2\left(E P+2 E P_2+E S P\right) }[/math]

Reversible Hill Mechanistic Scheme. k1 and k2 are the catalytic rate constants for the conversion of substrate to product and vice versa. The factor of two is due to the fact that the fully bound dimers have twice the rate of conversion (because there are two possibilities).

Invoking the rapid-equilibrium assumption we can write the various complexes in terms of equilibrium constants to give:

[math]\displaystyle{ v =\frac{V_f \sigma(1-\rho)(\sigma+\pi)}{1+(\sigma+\pi)^2} }[/math]

where [math]\displaystyle{ \rho=\Gamma / K_{eq} }[/math]. The [math]\displaystyle{ \sigma }[/math] and [math]\displaystyle{ \pi }[/math] terms are the ratio of substrate and product to their respective half-saturation constants, namely [math]\displaystyle{ \sigma = S/S_{0.5} }[/math] and [math]\displaystyle{ \pi = P/P_{0.5} }[/math] and

Using the author's own notation, if an enzyme has [math]\displaystyle{ h }[/math] sites that can bind ligand, the form, in the general case, can be shown to be:

[math]\displaystyle{ v = \frac{V_f \sigma(1-\rho)(\sigma+\pi)^{h-1}}{1+(\sigma+\pi)^h} }[/math]

The non-cooperative reversible Michaelis-Menten equation can be seen to emerge when we set the Hill coefficient to one.

If the enzyme is irreversible the equation turns into the simple Michaelis-Menten equation that is irreversible. When setting the equilibrium constant to infinity, the equation can be seen to revert to the simpler case where the product inhibits the reverse step.

A comparison has been made between the MWC and reversible Hill equation.[9]

A modification of the reversible Hill equation was published by Westermark et al[10] where modifiers affected the catalytic properties instead. This variant was shown to provide a much better fit for describing the kinetics of muscle phosphofructokinase.

References

  1. 1.0 1.1 Hofmeyr, Jan-Hendrik S.; Cornish-Bowden, Hofmeyr (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. 
  2. Popova, S.V.; Sel'kov, E.E. (15 May 1975). "Generalization of the model by Monod, Wyman and Changeux for the case of a reversible monosubstrate reaction". FEBS Letters 53 (3): 269–273. doi:10.1016/0014-5793(75)80034-2. 
  3. Saa, Pedro A.; Nielsen, Lars K. (December 2017). "Formulation, construction and analysis of kinetic models of metabolism: A review of modelling frameworks". Biotechnology Advances 35 (8): 981–1003. doi:10.1016/j.biotechadv.2017.09.005. 
  4. Liebermeister, Wolfram; Uhlendorf, Jannis; Klipp, Edda (15 June 2010). "Modular rate laws for enzymatic reactions: thermodynamics, elasticities and implementation". Bioinformatics 26 (12): 1528–1534. doi:10.1093/bioinformatics/btq141. 
  5. Westermark, Pål O.; Lansner, Anders (July 2003). "A Model of Phosphofructokinase and Glycolytic Oscillations in the Pancreatic β-cell". Biophysical Journal 85 (1): 126–139. doi:10.1016/S0006-3495(03)74460-9. PMID 12829470. 
  6. Miskovic, Ljubisa; Béal, Jonas; Moret, Michael; Hatzimanikatis, Vassily (20 August 2019). "Uncertainty reduction in biochemical kinetic models: Enforcing desired model properties". PLOS Computational Biology 15 (8): e1007242. doi:10.1371/journal.pcbi.1007242. PMID 31430276. 
  7. Cornish-Bowden, Athel (2012). Fundamentals of enzyme kinetics (4., completely revised and greatly enlarged ed.). Weinheim: Wiley-Blackwell. ISBN 978-3-527-33074-4. 
  8. Sauro, Herbert M. (2013). Enzyme kinetics for systems biology (2. ed.). New York, NY: Ambrosius Publishing. ISBN 978-0982477335. 
  9. Olivier, B.G.; Rohwer, J.M.; Snoep, J.L.; Hofmeyr, J.-H.S. (2006). "Comparing the regulatory behaviour of two cooperative, reversible enzyme mechanisms". IEE Proceedings - Systems Biology 153 (5): 335. doi:10.1049/ip-syb:20060020. 
  10. Westermark, Pål; Hällgren Kotaleski, Jeanette; Lansner, Anders (2004). "Derivation of a reversible Hill equation with modifiers affecting catalytic properties". WSEAS Transactions on Biology and Biomedicine 1: 91–98. http://kth.diva-portal.org/smash/record.jsf?pid=diva2%3A10651&dswid=-448.