Michaelis–Menten–Monod kinetics explained

For Michaelis–Menten–Monod (MMM) kinetics it is intended the coupling of an enzyme-driven chemical reaction of the Michaelis–Menten type[1] with the Monod growth of an organisms that performs the chemical reaction.[2] The enzyme-driven reaction can be conceptualized as the binding of an enzyme E with the substrate S to form an intermediate complex C, which releases the reaction product P and the unchanged enzyme E. During the metabolic consumption of S, biomass B is produced, which synthesizes the enzyme, thus feeding back to the chemical reaction. The two processes can be expressed as

where

k1

and

k-1

are the forward and backward equilibrium rate constants,

k

is the reaction rate constant for product release,

Y

is the biomass yield coefficient, and

z

is the enzyme yield coefficient.

Transient kinetics

The kinetic equations describing the reactions above can be derived from the GEBIK equations[3] and are written as

where

\muB

is the biomass mortality rate and

\muE

is the enzyme degradation rate. These equations describe the full transient kinetics, but cannot be normally constrained to experiments because the complex C is difficult to measure and there is no clear consensus on whether it actually exists.

Quasi-steady-state kinetics

Equations 3 can be simplified by using the quasi-steady-state (QSS) approximation, that is, for

d[C]
dt

=0

;[4] under the QSS, the kinetic equations describing the MMM problem become

where

K=(k-1+k)/k1

is the Michaelis–Menten constant (also known as the half-saturation concentration and affinity).

Implicit analytic solution

If one hypothesizes that the enzyme is produced at a rate proportional to the biomass production and degrades at a rate proportional to the biomass mortality, then Eqs. 4 can be rewritten as

where

S

,

P

,

E

,

B

are explicit function of time

t

. Note that Eq. (4b) and (4d) are linearly dependent on Eqs. (4a) and (4c), which are the two differential equations that can be used to solve the MMM problem. An implicit analytic solution[5] can be obtained if

P

is chosen as the independent variable and

t(P)

,

S(P)

,

E(P)

and

B(P

) are rewritten as functions of

P

so to obtain

where

S(t)

has been substituted by

S(P)=S0-P

as per mass balance

S0+P0=S+P

, with the initial value

S0=S(P)

when

P=P0=0

, and where

E(t)

has been substituted by

zB(P)

as per the linear relation

E=zB

expressed by Eq. (4d). The analytic solution to Eq. (5b) is

with the initial biomass concentration

B0=B(P)

when

P=0

. To avoid the solution of a transcendental function, a polynomial Taylor expansion to the second-order in

P

is used for

B(P)

in Eq. (6) as

Substituting Eq. (7) into Eq. (5a} and solving for

t(P)

with the initial value

t(P=0)=0

, one obtains the implicit solution for

t(P)

as

with the constants

For any chosen value of

P

, the biomass concentration can be calculated with Eq. (7) at a time

t(P)

given by Eq. (8). The corresponding values of

S(P)=S0-P

and

E(P)=zB(P)

can be determined using the mass balances introduced above.

See also

Notes and References

  1. Michaelis, L.; Menten, M. L. (1913). "Die Kinetik der Invertinwirkung". Biochem Z. 49: 333–369
  2. Monod J. (1949) The growth of bacterial cultures. Annu. Rev. Microbial. 3, 371–394
  3. Maggi F. and W. J. Riley, (2010), Mathematical treatment of isotopologue and isotopomer speciation and fractionation in biochemical kinetics, Geochim. Cosmochim. Acta,
  4. Briggs G.E.; Haldane, J.B.S., "A note on the kinetics of enzyme action", \textit \textbf, \textit, 338–339.
  5. Maggi F. and La Cecilia D., (2016), "An implicit analytic solution of Michaelis–Menten–Monod kinetics", American Chemical Society, ACS Omega 2016, 1, 894−898,