Misplaced Pages

Reversible Michaelis–Menten kinetics

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Enzyme kinetics for reversible reactions

Enzymes are proteins that act as biological catalysts by accelerating chemical reactions. Enzymes act on small molecules called substrates, which an enzyme converts into products. Almost all metabolic processes in the cell need enzyme catalysis in order to occur at rates fast enough to sustain life. The study of how fast an enzyme can transform a substrate into a product is called enzyme kinetics.

The rate of reaction of many chemical reactions shows a linear response as function of the concentration of substrate molecules. Enzymes however display a saturation effect where,, as the substrate concentration is increased the reaction rate reaches a maximum value. Standard approaches to describing this behavior are based on models developed by Michaelis and Menten as well and Briggs and Haldane. Most elementary formulations of these models assume that the enzyme reaction is irreversible, that is product is not converted back to substrate. However, this is unrealistic when describing the kinetics of enzymes in an intact cell because there is product available. Reversible Michaelis–Menten kinetics, using the reversible form of the Michaelis–Menten equation, is therefore important when developing computer models of cellular processes involving enzymes.

In enzyme kinetics, the Michaelis–Menten kinetics kinetic rate law that describes the conversion of one substrate to one product, is often commonly depicted in its irreversible form as:

v = V max s K m + s {\displaystyle v={\frac {V_{\max }s}{K_{\mathrm {m} }+s}}}

where v {\displaystyle v} is the reaction rate, V max {\displaystyle V_{\max }} is the maximum rate when saturating levels of the substrate are present, K m {\displaystyle K_{\mathrm {m} }} is the Michaelis constant and s {\displaystyle s} the substrate concentration.

In practice, this equation is used to predict the rate of reaction when little or no product is present. Such situations arise in enzyme assays. When used to model enzyme rates in vivo , for example, to model a metabolic pathway, this representation is inadequate because under these conditions product is present. As a result, when building computer models of metabolism or other enzymatic processes, it is better to use the reversible form of the Michaelis–Menten equation.

To model the reversible form of the Michaelis–Menten equation, the following reversible mechanism is considered:

E + S k 1 k 1 ES k 2 k 2 E + P {\displaystyle {\ce {{E}+{S}<=>ES<=>{E}+{P}}}}

To derive the rate equation, it is assumed that the concentration of enzyme-substrate complex is at steady-state, that is d e s / d t = 0 {\displaystyle des/dt=0} .

Following current literature convention, we will be using lowercase Roman lettering to indicate concentrations (this avoids cluttering the equations with square brackets). Thus e s {\displaystyle es} indicates the concentration of enzyme-substrate complex, ES.

The net rate of change of product (which is equal to v {\displaystyle v} ) is given by the difference in forward and reverse rates:

v = v f v r = k 2 e s k 2 e   p {\displaystyle v=v_{f}-v_{r}=k_{2}es-k_{-2}e\ p}

The total level of enzyme moiety is the sum total of free enzyme and enzyme-complex, that is e t = e + e s {\displaystyle e_{t}=e+es} . Hence the level of free e {\displaystyle e} is given by the difference in the total enzyme concentration, e t {\displaystyle e_{t}} and the concentration of complex, that is:

e = e t e s {\displaystyle e=e_{t}-es}

Using mass conservation we can compute the rate of change of e s {\displaystyle es} using the balance equation:

d e s d t = k 1 ( e t e s ) s + k 2 ( e t e s ) p ( k 1 + k 2 ) e s = 0 {\displaystyle {\frac {des}{dt}}=k_{1}\left(e_{t}-es\right)s+k_{-2}\left(e_{t}-es\right)p-\left(k_{-1}+k_{2}\right)es=0}

where e {\displaystyle e} has been replaced using e = e t e s {\displaystyle e=e_{t}-es} . This leaves e s {\displaystyle es} as the only unknown. Solving for e s {\displaystyle es} gives:

e s = e t ( k 1 s + k 2   p ) k 1 + k 2 + k 1   s + k 2   p {\displaystyle es={\frac {\mathrm {e_{t}} \left(k_{1}s+k_{-2}\ p\right)}{k_{-1}+k_{2}+k_{1}\ s+k_{-2}\ p}}}

Inserting e s {\displaystyle es} into the rate equation v = k 2 e s k 2 e   p {\displaystyle v=k_{2}es-k_{-2}e\ p} and rearranging gives:

v = e t k 1 k 2 s k 1 k 2 p k 1 + k 2 + k 1 s + k 2 p {\displaystyle v=e_{t}{\frac {k_{1}k_{2}s-k_{-1}k_{-2}p}{k_{-1}+k_{2}+k_{1}s+k_{-2}p}}}

The following substitutions are now made:

k 2 = V max f e t ; K m s = k 1 + k 2 k 1 {\displaystyle k_{2}={\frac {V_{\max }^{f}}{e_{t}}};\quad K_{m}^{s}={\frac {k_{-1}+k_{2}}{k_{1}}}}

and

k 2 = V max r e t ; K m p = k 1 + k 2 k 2 {\displaystyle k_{-2}={\frac {V_{\max }^{r}}{e_{t}}};\quad K_{m}^{p}={\frac {k_{-1}+k_{2}}{k_{-2}}}}

after rearrangement, we obtain the reversible Michaelis–Menten equation in terms of four constants:

v = V max f K m s s V max r K m p p 1 + s K m s + p K m p {\displaystyle v={\frac {{\frac {V_{\max }^{f}}{K_{m}^{s}}}s-{\frac {V_{\max }^{r}}{K_{m}^{p}}}p}{1+{\frac {s}{K_{m}^{s}}}+{\frac {p}{K_{m}^{p}}}}}}

Haldane relationship

This is not the usual form in which the equation is used. Instead, the equation is set to zero, meaning v = 0 {\displaystyle v=0} , indicating we are at equilibrium and the concentrations s {\displaystyle s} and p {\displaystyle p} are now equilibrium concentrations, hence:

0 = V max f s e q / K m s V max r p e q / K m p {\displaystyle 0=V_{\max }^{f}s_{eq}/K_{m}^{s}-V_{\max }^{r}p_{eq}/K_{m}^{p}}

Rearranging this gives the so-called Haldane relationship:

K e q = p e q s e q = V max f K m p V max r K m s {\displaystyle K_{eq}={\frac {p_{eq}}{s_{eq}}}={\frac {V_{\max }^{f}K_{m}^{p}}{V_{\max }^{r}K_{m}^{s}}}}

The advantage of this is that one of the four constants can be eliminated and replaced with the equilibrium constant which is more likely to be known. In addition, it allows one to make a useful interpretation in terms of the thermodynamic and saturation effects (see next section). Most often the reverse maximum rate is eliminated to yield the final equation:

v = V max f / K m S ( s p / K e q ) 1 + s / K m s + p / K m p {\displaystyle v={\frac {V_{\max }^{f}/K_{m}^{S}\left(s-p/K_{eq}\right)}{1+s/K_{m}^{s}+p/K_{m}^{p}}}}

Decomposition of the rate law

The reversible Michaelis–Menten law, as with many enzymatic rate laws, can be decomposed into a capacity term, a thermodynamic term, and an enzyme saturation level. This is more easily seen when we write the reversible rate law as:

v = V max f ( s p / K e q ) 1 1 + s / K m s + p / K m p {\displaystyle v=V_{\max }^{f}\cdot \left(s-p/K_{eq}\right)\cdot {\frac {1}{1+s/K_{m}^{s}+p/K_{m}^{p}}}}

where V max f {\displaystyle V_{\max }^{f}} is the capacity term, ( s p / K e q ) {\displaystyle \left(s-p/K_{eq}\right)} the thermodynamic term and

1 1 + s / K m s + p / K m p {\displaystyle {\frac {1}{1+s/K_{m}^{s}+p/K_{m}^{p}}}}

the saturation term. The separation can be even better appreciated if we look at the elasticity coefficient ε s v {\displaystyle \varepsilon _{s}^{v}} . According to elasticity algebra, the elasticity of a product is the sum of the sub-term elasticities, that is:

ε x a b = ε x a + ε x b {\displaystyle \varepsilon _{x}^{ab}=\varepsilon _{x}^{a}+\varepsilon _{x}^{b}}

Hence the elasticity of the reversible Michaelis–Menten rate law can easily be shown to be:

ε s v = ε s v c a p + ε s v t h e r m o + ε s v s a t {\displaystyle \varepsilon _{s}^{v}=\varepsilon _{s}^{v_{cap}}+\varepsilon _{s}^{v_{thermo}}+\varepsilon _{s}^{v_{sat}}}

Since the capacity term is a constant, the first elasticity is zero. The thermodynamic term can be easily shown to be:

ε s v t h e r m o = 1 1 ρ {\displaystyle \varepsilon _{s}^{v_{thermo}}={\frac {1}{1-\rho }}}

where ρ {\displaystyle \rho } is the disequilibrium ratio and equals Γ / K e q {\displaystyle \Gamma /K_{eq}} and Γ {\displaystyle \Gamma } the mass–action ratio

The saturation term becomes:

ε s v s a t = s / K m s 1 + s / K m s + p / K m p {\displaystyle \varepsilon _{s}^{v_{sat}}={\frac {-s/K_{m}^{s}}{1+s/K_{m}^{s}+p/K_{m}^{p}}}}

References

  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. PMID 9283752.
  2. Briggs, GE; Haldane, JB (1925). "A Note on the Kinetics of Enzyme Action". The Biochemical Journal. 19 (2): 338–9. doi:10.1042/bj0190338. PMC 1259181. PMID 16743508.
  3. Cornish-Bowden, Athel (2012). Fundamentals of enzyme kinetics (4., completely revised and greatly enlarged ed.). Weinheim: Wiley-Blackwell. ISBN 9783527330744.
  4. Rohwer, Johann M.; Hofmeyr, Jan-Hendrik S. (16 December 2010). "Kinetic and Thermodynamic Aspects of Enzyme Control and Regulation". The Journal of Physical Chemistry B. 114 (49): 16280–16289. doi:10.1021/jp108412s. PMID 21028763.
  5. Noor, Elad; Flamholz, Avi; Liebermeister, Wolfram; Bar-Even, Arren; Milo, Ron (2 September 2013). "A note on the kinetics of enzyme action: A decomposition that highlights thermodynamic effects". FEBS Letters. 587 (17): 2772–2777. doi:10.1016/j.febslet.2013.07.028. PMID 23892083. S2CID 41521250.
  6. Woods, John H.; Sauro, Herbert M. (1997). "#Elasticities in Metabolic Control Analysis: algebraic derivation of simplified expressions". Bioinformatics. 13 (2): 123–130. doi:10.1093/bioinformatics/13.2.123. PMID 9146958.
Categories: