The Michaelis-Menten Elimination Rate Law in PBPK Models

Why a Saturable Rate Law Is Used

A PBPK model needs a rule connecting the concentration available to an eliminating enzyme with the rate at which the enzyme removes compound. At low concentrations, a rate proportional to concentration may be an adequate approximation. But an enzyme has finite capacity: as concentration rises, its rate approaches a maximum rather than increasing without bound. The Michaelis-Menten rate law represents this saturable relationship with two parameters.

Here \(C\) means the substrate concentration at the enzyme site, or the concentration variable the PBPK model uses to represent it. It may be an unbound concentration, a total concentration, or another defined concentration measure, depending on how the kinetic parameters were obtained. The concentration definition used with \(C\) must match the definition used for \(K_m\).

Definition: The Michaelis-Menten elimination rate law is \(v=V_{\max}C/(K_m+C)\), where \(v\) is the elimination rate, \(C\) is the relevant substrate concentration, \(V_{\max}\) is the maximum rate represented by the model, and \(K_m\) is the concentration at which the rate is half of \(V_{\max}\). \(C\) and \(K_m\) have the same concentration units. \(v\) and \(V_{\max}\) have the same rate units.

For example, if \(C\) and \(K_m\) are measured in micromoles per litre and \(v\) is expressed in nanomoles per minute per milligram of microsomal protein, then \(V_{\max}\) must also be in nanomoles per minute per milligram of microsomal protein. The denominator \(K_m+C\) has concentration units, so the fraction \(C/(K_m+C)\) is dimensionless. Multiplying this fraction by \(V_{\max}\) therefore gives a rate in the same units as \(V_{\max}\).

A PBPK implementation may use a rate normalized to tissue volume, such as micromoles per litre of tissue per minute, or a total rate for an entire organ, such as micromoles per minute. The choice does not change the equation, but it does change what \(V_{\max}\) represents. In a tissue-volume-normalized formulation, converting \(v\) to an organ amount per time requires multiplication by the relevant tissue volume.

Where the Equation Comes From

The standard derivation considers an enzyme \(E\) that binds substrate \(S\) to form an enzyme-substrate complex \(ES\), which can then form product \(P\). It is useful to show the derivation because it clarifies what the two parameters represent. Let \(C\) denote the substrate concentration in this derivation, and let \(E_T\) be the total enzyme concentration.

E+S \xrightleftharpoons[k_{-1}]{k_1} ES \xrightarrow{k_{\mathrm{cat}}} E+P

The usual initial-rate derivation assumes one effective substrate-enzyme system, a well-mixed reaction, constant total enzyme, negligible product accumulation and reverse product reaction during the measurement, and a steady concentration of \(ES\) over the rate measurement. The steady-state assumption means that the rate of forming \(ES\) equals the rate of its breakdown. Since free enzyme concentration is \(E_T-[ES]\), this gives the following balance.

\frac{d[ES]}{dt} = k_1 C(E_T-[ES])-(k_{-1}+k_{\mathrm{cat}})[ES]=0

Rearranging the balance and defining \(K_m=(k_{-1}+k_{\mathrm{cat}})/k_1\) gives the concentration of enzyme-substrate complex. Product forms at rate \(k_{\mathrm{cat}}[ES]\). Substituting the complex concentration into that rate produces the Michaelis-Menten form, with \(V_{\max}=k_{\mathrm{cat}}E_T\).

[ES]=\frac{E_T C}{K_m+C}, \qquad v=k_{\mathrm{cat}}[ES] =\frac{k_{\mathrm{cat}}E_T C}{K_m+C} =\frac{V_{\max}C}{K_m+C}

This derivation describes a simplified enzyme system. In a PBPK model, the same equation can also serve as a compact representation of measured or scaled metabolic activity. In that use, the parameters should be interpreted on the basis and under the conditions for which they were determined; the equation alone does not establish that every molecular assumption of the derivation is exactly satisfied in vivo.

Interpreting the Parameters

\(V_{\max}\) sets the rate scale. For nonnegative \(C\), the fraction \(C/(K_m+C)\) is between zero and one, so the predicted rate does not exceed \(V_{\max}\). As \(C\) becomes large relative to \(K_m\), the fraction approaches one and the rate approaches \(V_{\max}\). Thus, \(V_{\max}\) is the maximum rate represented by this pathway under the model’s specified enzyme abundance and rate basis.

\(K_m\) identifies a particular point on the rate curve. Substituting \(C=K_m\) into the law gives \(v=V_{\max}/2\). For example, if \(K_m=20\) micromoles per litre and \(V_{\max}=0.90\) micromoles per minute, then at \(C=20\) micromoles per litre the predicted rate is \(0.45\) micromoles per minute. This interpretation is part of the rate law; it does not mean that every biological system has an abrupt switch from unsaturated to saturated at \(K_m\).

The parameter values must also correspond to the enzyme capacity being modelled. An enzyme-normalized \(V_{\max}\), such as a rate per milligram of protein, is not automatically a whole-organ rate. It needs an appropriate scaling to the enzyme amount or tissue basis in the PBPK model. Likewise, \(K_m\) is a concentration and cannot be combined directly with an amount unless the model converts that amount to a concentration using a specified volume.

Worked Example: Rate at the Half-Maximal Concentration

Suppose a hypothetical enzyme assay gives \(V_{\max}=120\) nanomoles per minute per milligram of protein and \(K_m=30\) micromoles per litre. Calculate the rate at \(C=30\) micromoles per litre.

The rate law is \(v=V_{\max}C/(K_m+C)\). The concentration units match, and the rate will have the same units as \(V_{\max}\).

v =\frac{(120\ \mathrm{nmol\,min^{-1}\,mg^{-1}})(30\ \mathrm{\mu mol\,L^{-1}})} {30\ \mathrm{\mu mol\,L^{-1}}+30\ \mathrm{\mu mol\,L^{-1}}} =120\left(\frac{30}{60}\right) =60\ \mathrm{nmol\,min^{-1}\,mg^{-1}}

At \(C=K_m\), the predicted rate is \(60\) nanomoles per minute per milligram of protein, exactly half of \(V_{\max}=120\) nanomoles per minute per milligram of protein.

Worked Example: Rates Across Three Concentrations

For a hypothetical pathway, let \(V_{\max}=0.90\) micromoles per minute and \(K_m=20\) micromoles per litre. Calculate rates at \(C=5\), \(20\), and \(80\) micromoles per litre.

At \(C=5\), substitution gives \(v=0.90(5)/(20+5)=0.18\) micromoles per minute. At \(C=20\), \(v=0.90(20)/(20+20)=0.45\) micromoles per minute. At \(C=80\), \(v=0.90(80)/(20+80)=0.72\) micromoles per minute.

\begin{aligned} C=5\ \mathrm{\mu mol\,L^{-1}} &: \quad v=0.90\left(\frac{5}{25}\right)=0.18\ \mathrm{\mu mol\,min^{-1}}\\ C=20\ \mathrm{\mu mol\,L^{-1}} &: \quad v=0.90\left(\frac{20}{40}\right)=0.45\ \mathrm{\mu mol\,min^{-1}}\\ C=80\ \mathrm{\mu mol\,L^{-1}} &: \quad v=0.90\left(\frac{80}{100}\right)=0.72\ \mathrm{\mu mol\,min^{-1}} \end{aligned}

The rates increase as concentration increases, but they remain below \(0.90\) micromoles per minute. At \(C=80\), the rate is \(0.80V_{\max}\), since \(80/(20+80)=0.80\). These are predictions from the specified rate law, not measured values from a real study.

Using the Rate on a PBPK Tissue Basis

A PBPK model often tracks amounts in compartments, while kinetic parameters may be expressed per unit tissue volume. The rate law can be evaluated on the concentration basis first; if its result is a rate per unit volume, multiplying by the compartment volume gives the total amount eliminated per time. The concentration and volume must refer to compatible definitions and bases.

Worked Example: Converting a Tissue-Normalized Rate

Consider a hypothetical, well-mixed tissue compartment of \(1.5\) litres. Suppose its modelled enzyme rate has \(V_{\max}=0.60\) micromoles per litre of tissue per minute and \(K_m=12\) micromoles per litre. At a tissue concentration of \(C=6\) micromoles per litre, calculate both the rate per litre of tissue and the total compartment rate.

First evaluate the rate law. The numerator and denominator use the same concentration units, and the resulting rate is on a per-litre-of-tissue basis.

v =\frac{(0.60\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}) (6\ \mathrm{\mu mol\,L^{-1}})} {12\ \mathrm{\mu mol\,L^{-1}}+6\ \mathrm{\mu mol\,L^{-1}}} =0.60\left(\frac{6}{18}\right) =0.20\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}

The whole-compartment rate is the rate per litre multiplied by the tissue volume: \(0.20\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\times1.5\ \mathrm{L} =0.30\ \mathrm{\mu mol\,min^{-1}}\). Thus, the model predicts elimination of \(0.30\) micromoles per minute from this compartment at that concentration, under the stated assumptions.

Assumptions, Checks, and Alternatives

Before applying the law, check that its concentration input represents the substrate available to the modelled enzyme and that its kinetic parameters use the same concentration definition. Confirm that \(V_{\max}\) is on the intended basis—such as per milligram of protein, per litre of tissue, or for an entire organ—and that the model does not multiply or scale that capacity a second time. If enzyme abundance changes, a fixed \(V_{\max}\) may no longer represent the changed capacity.

The one-enzyme Michaelis-Menten form may be inadequate when other processes substantially affect the rate, such as multiple unresolved pathways, substrate interactions, inhibition, or transport into the enzyme’s location. A mismatch between observed and predicted rates across the concentration range can signal that the assumed form or the input data are unsuitable. Also check the recorded concentration range: measurements confined to a narrow range may not clearly distinguish the curvature or the limiting rate. Estimating the parameters from rate data is a separate task; this tutorial uses specified parameters to calculate and interpret rates.

The main simple alternatives make different assumptions about how rate responds to concentration. A first-order rate law, \(v=kC\), predicts a rate proportional to concentration and does not impose a finite maximum rate. A constant zero-order rate does not vary with concentration over the range in which it is used. The Michaelis-Menten law provides a concentration-dependent rate with a finite \(V_{\max}\), and it reduces to the half-maximal rate at \(C=K_m\). These forms are not interchangeable: the choice should reflect the intended mechanism and the concentration range the model must represent.

Common pitfall: It is incorrect to say that \(K_m\) is the maximum concentration the enzyme can handle, or to treat it as a rate. \(K_m\) has concentration units; it is the concentration at which this rate law predicts \(v=V_{\max}/2\). The rate approaches \(V_{\max}\) as concentration increases, and \(V_{\max}\) has rate units.

Key takeaway: Keep three checks together: \(C\) and \(K_m\) must share concentration units and concentration definitions; \(v\) and \(V_{\max}\) must share rate units and a consistent tissue or enzyme basis; and at \(C=K_m\), the equation must give \(v=V_{\max}/2\).

Check Your Understanding

Try each question, then compare your reasoning with the solutions below.

If \(V_{\max}=50\) nanomoles per minute and \(K_m=10\) micromoles per litre, what is the rate at \(C=10\) micromoles per litre?
For the same parameters, calculate the rate at \(C=30\) micromoles per litre.
A model uses \(C\) in micromoles per litre but \(K_m\) in milligrams per litre. What must be addressed before evaluating the equation?
A pathway has a \(V_{\max}\) expressed per milligram of enzyme preparation. Can this value be used directly as the total rate for a whole organ? Explain.
What do \(V_{\max}\) and \(K_m\) represent, and what are their units relative to \(v\) and \(C\)?

Solutions. (1) At \(C=K_m\), \(v=50(10)/(10+10)=25\) nanomoles per minute. (2) \(v=50(30)/(10+30)=37.5\) nanomoles per minute. (3) Convert one concentration so \(C\) and \(K_m\) use the same units and compatible concentration definitions. (4) No. The enzyme-normalized capacity must be scaled to an appropriate enzyme amount or tissue basis before it represents a whole-organ rate. (5) \(V_{\max}\) is the maximum rate represented by the law and has the same units as \(v\); \(K_m\) is the concentration at which \(v\) is half of \(V_{\max}\) and has the same units as \(C\).