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
and
are the forward and backward equilibrium rate constants,
is the reaction rate constant for product release,
is the biomass yield coefficient, and
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
is the biomass mortality rate and
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
;
[4] under the QSS, the kinetic equations describing the MMM problem become
where
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
,
,
,
are explicit function of time
. 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
is chosen as the independent variable and
,
,
and
) are rewritten as functions of
so to obtain
where
has been substituted by
as per mass balance
, with the initial value
when
, and where
has been substituted by
as per the linear relation
expressed by Eq. (4d). The analytic solution to Eq. (5b) is
with the initial biomass concentration
when
. To avoid the solution of a
transcendental function, a polynomial Taylor expansion to the second-order in
is used for
in Eq. (6) as
Substituting Eq. (7) into Eq. (5a} and solving for
with the initial value
, one obtains the implicit solution for
as
with the constants
For any chosen value of
, the biomass concentration can be calculated with Eq. (7) at a time
given by Eq. (8). The corresponding values of
and
can be determined using the mass balances introduced above.
See also
Notes and References
- Michaelis, L.; Menten, M. L. (1913). "Die Kinetik der Invertinwirkung". Biochem Z. 49: 333–369
- Monod J. (1949) The growth of bacterial cultures. Annu. Rev. Microbial. 3, 371–394
- Maggi F. and W. J. Riley, (2010), Mathematical treatment of isotopologue and isotopomer speciation and fractionation in biochemical kinetics, Geochim. Cosmochim. Acta,
- Briggs G.E.; Haldane, J.B.S., "A note on the kinetics of enzyme action", \textit \textbf, \textit, 338–339.
- 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,