Michaelis-Menten Two-Compartment Model

• The Michaelis–Menten two-compartment model is a framework that describes blood-liver exchange with saturable hepatic metabolism via a nonlinear rate law.
• It is derived by reducing a three-compartment ethanol cascade, treating gut absorption as an external input while assuming well-mixed compartments.
• The model ensures positivity and global stability under proper conditions and is linked to spatial reaction–diffusion reductions in kinetic modeling.

Searching arXiv for the cited papers to ground the article in the relevant literature. The Michaelis–Menten two-compartment model is a compartmental dynamical system in which material is exchanged between two well-mixed compartments while elimination or conversion occurs from one compartment through a saturable Michaelis–Menten rate law. In the ethanol-metabolism setting, the canonical realization is the blood–liver subsystem obtained from a gut–blood–liver cascade by suppressing the gut state or treating gut absorption as an external input to blood. In broader kinetic modeling, closely related two-compartment forms also arise as reduced descriptions of spatial reaction–diffusion systems and bulk–membrane transport problems (Hoang et al., 29 Jun 2026).

1. Canonical compartment structure

In the ethanol framework, the underlying continuous-time model is a three-compartment cascade with state variables $A(t)$, $B(t)$, and $C(t)$, representing gastrointestinal amount, blood amount, and liver amount, respectively. Gut empties into blood at rate $a$, blood and liver exchange ethanol symmetrically at rate $b$, and hepatic metabolism removes ethanol from the liver through a nonlinear rate $f(C)$. The governing equations are

$$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$

The two-compartment specialization is obtained by dropping the gut dynamics, either by setting $A(t)\equiv 0$ or by considering the long-time limit $A(t)=A(0)e^{-at}\to 0$. The resulting blood–liver subsystem is

$$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$

If gut absorption is retained as an external forcing, one writes $B(t)$0 and obtains

$B(t)$1

This formulation rests on standard compartmental assumptions: well-mixed compartments, constant volumes, constant inter-compartment rate(s), constant absorption rate $B(t)$2, and time-invariant parameters. Primary elimination is hepatic; renal elimination and reabsorption are neglected, no endogenous ethanol production is included, and the nonlinear metabolic term $B(t)$3 is assumed continuous, nonnegative, zero only at $B(t)$4, and bounded above, expressing finite metabolic capacity (Hoang et al., 29 Jun 2026).

2. Michaelis–Menten specialization

The Michaelis–Menten specialization enters through the liver compartment. Let $B(t)$5 and $B(t)$6 denote blood and liver volumes, and let $B(t)$7 be liver concentration. The classical hepatic Michaelis–Menten rate in concentration form is

$B(t)$8

When elimination is written in the amount equation, the metabolic term becomes

$B(t)$9

Accordingly, the two-compartment Michaelis–Menten system in amount form is

$C(t)$0

A key structural point is that elimination depends on liver concentration but appears in the amount equation through $C(t)$1; equivalently, the nonlinearity is expressed in the liver amount $C(t)$2.

The same saturating form appears in the generalized ethanol model through

$C(t)$3

with the parameter mapping $C(t)$4 and $C(t)$5. Under this identification, the earlier three-compartment model from Wacker is exactly the Michaelis–Menten case,

$C(t)$6

and the classical two-compartment Michaelis–Menten model is recovered by restricting to the blood–liver subsystem with optional input $C(t)$7. In the terminology of the ethanol paper, this is the Levitt two-compartment model structure specialized to symmetric exchange and Michaelis–Menten hepatic elimination (Hoang et al., 29 Jun 2026).

3. Positivity, boundedness, and equilibrium structure

For nonnegative initial amounts, the continuous-time model preserves nonnegativity. The proof evaluates the vector field on the boundary of the nonnegative orthant and establishes invariance. This ensures that the model is dynamically consistent as a mass-balance system.

Boundedness follows from the total amount

$C(t)$8

which satisfies

$C(t)$9

Hence $a$0 is nonincreasing. In the unforced case $a$1, the compartment amounts are bounded and eventually decrease.

For the two-compartment blood–liver subsystem, the Lyapunov function

$a$2

yields

$a$3

with equality if and only if $a$4. By Lyapunov’s direct method, the equilibrium $a$5 is globally asymptotically stable. Since $a$6 and the full model is cascade-connected, the equilibrium $a$7 is globally asymptotically stable via Seibert–Suarez cascade stabilization. In the Michaelis–Menten case,

$a$8

the derivative becomes

$a$9

so the same conclusion applies directly (Hoang et al., 29 Jun 2026).

The equilibrium structure depends critically on input. Without external input, the only equilibrium is the origin. With constant input $b$0, the steady state satisfies

$b$1

which reduces to

$b$2

For Michaelis–Menten elimination,

$b$3

which exists and is unique if $b$4, with

$b$5

If $b$6, there is no finite steady state. A common misunderstanding is therefore to expect a nonzero equilibrium in the unforced model; in fact, nonzero steady states require sustained input below maximal elimination capacity (Hoang et al., 29 Jun 2026).

4. Discrete-time formulation and numerical behavior

A forward-Euler discretization of the three-compartment system is given by

$b$7

with $b$8. For the two-compartment Michaelis–Menten model with blood input $b$9, the updates are

$f(C)$0

The paper defines

$f(C)$1

lets $f(C)$2, and sets

$f(C)$3

The step-size condition

$f(C)$4

preserves positivity and boundedness and yields local and global stability. In the Michaelis–Menten case,

$f(C)$5

so

$f(C)$6

Under the same step-size restriction and with $f(C)$7, the discrete-time system is positive and globally asymptotically stable to $f(C)$8.

The numerical experiments reported for nonlinear hepatic rate functions include the monotone saturating law

$f(C)$9

which is a Michaelis–Menten form with $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$0 and $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$1. Representative parameters are $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$2, $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$3, $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$4, and $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$5, hence $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$6 and $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$7. The reported qualitative behavior is that solutions remain positive and decay to zero in the absence of sustained input, while saturation produces a prolonged tail in the decay of liver and blood amounts; nonmonotone rates produce a slower approach to equilibrium than monotone saturating forms (Hoang et al., 29 Jun 2026).

5. Derivations from spatial and reaction–diffusion models

The two-compartment Michaelis–Menten model is not only a pharmacokinetic ansatz; it also emerges from spatial reductions. In the reaction–diffusion setting, Frank, Lax, Walcher, and Wittich analyze irreversible and reversible Michaelis–Menten systems with diffusion under small total enzyme and slow diffusion. With the scaling $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$8, $$\begin{aligned} A'(t) &= -a\,A(t), \ B'(t) &= b\,C(t) - b\,B(t) + a\,A(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$9, $A(t)\equiv 0$0, and slow time $A(t)\equiv 0$1, the complex variable $A(t)\equiv 0$2 is fast and the remaining variables are slow. In the irreversible case, the quasi-steady-state manifold is

$A(t)\equiv 0$3

which yields the familiar local Michaelis–Menten rate

$A(t)\equiv 0$4

Under the same reduction, a two-compartment slow-time model is obtained by coupling well-mixed compartments $A(t)\equiv 0$5 and $A(t)\equiv 0$6 through exchange terms. For the irreversible case,

$A(t)\equiv 0$7

with analogous equations for the total enzyme variables $A(t)\equiv 0$8 and $A(t)\equiv 0$9. If $A(t)=A(0)e^{-at}\to 0$0, enzyme dynamics decouple into simple exchange; if $A(t)=A(0)e^{-at}\to 0$1, an additional coupling through the quasi-steady-state complex appears. The paper reports first-order convergence of the full reaction–diffusion solution to the reduced solution as $A(t)=A(0)e^{-at}\to 0$2 on fixed time intervals, with full and reduced systems visually indistinguishable at $A(t)=A(0)e^{-at}\to 0$3 in the reported simulations (Frank et al., 2017).

A different spatial route to a two-compartment Michaelis–Menten description appears in membrane-receptor modeling. In a spherical-cell geometry, the extracellular region forms a bulk compartment with concentration $A(t)=A(0)e^{-at}\to 0$4, while the membrane acts as a receptor compartment characterized by surface concentration $A(t)=A(0)e^{-at}\to 0$5 and fractional occupancy $A(t)=A(0)e^{-at}\to 0$6. The membrane kinetics are

$A(t)=A(0)e^{-at}\to 0$7

and the diffusive interfacial flux can be written as

$A(t)=A(0)e^{-at}\to 0$8

At steady state, eliminating $A(t)=A(0)e^{-at}\to 0$9 leads to a spatially corrected Michaelis–Menten law

$$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$0

where

$$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$1

Here $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$2, $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$3, and $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$4. The paper states that naive well-mixed calculations significantly overestimate reaction rates in certain biophysical parameter regimes. In the two-compartment interpretation, this overestimation is a consequence of finite interfacial mass transfer and receptor competition, which inflate the apparent half-saturation at low substrate concentration (Cengiz et al., 23 Jan 2025).

6. Identification, scaling, and limitations

For the ethanol two-compartment Michaelis–Menten model, identifiability depends strongly on observability and input design. With only blood measurements, the parameters $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$5 may be only partially identifiable because liver amount $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$6 is unobserved. If gut dynamics are included, $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$7 and the initial gut amount $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$8 are confounded in early-time data. Structural identifiability improves when both blood and liver concentrations are observed or when multiple dosing regimens are used to excite the system. The paper therefore recommends nonlinear least squares or Bayesian inference on $$\begin{aligned} B'(t) &= b\,C(t) - b\,B(t), \ C'(t) &= b\,B(t) - b\,C(t) - f\!\big(C(t)\big). \end{aligned}$$9 using the two-compartment Michaelis–Menten ODEs, with $B(t)$00 treated as a latent state and filtering methods such as extended Kalman or particle filters used to handle measurement noise. It also notes that high-dose data are informative for $B(t)$01, while low-dose data are informative for $B(t)$02 (Hoang et al., 29 Jun 2026).

The same source identifies two operational quantities of interest. Time-to-peak in blood is obtained numerically from the condition $B(t)$03 once $B(t)$04 and $B(t)$05 are known or estimated; early dynamics are dominated by input, whereas later decline reflects exchange and hepatic removal. The area under the blood curve is obtained by integrating $B(t)$06 over time; for repeated dosing one sums AUCs or simulates over long windows.

Several limitations delimit the classical formulation. The assumption of symmetric exchange at rate $B(t)$07 is explicitly described as a simplification, and more realistic pharmacokinetics may require $B(t)$08. Only hepatic elimination is modeled; extrahepatic pathways and renal clearance are neglected. Without sustained input, the only equilibrium is the origin, and although constant-input steady states exist for $B(t)$09, those nonzero equilibria are not analyzed in the stability proofs. A broader implication of the spatial papers is that compartmental Michaelis–Menten models are most reliable when the well-mixed approximation is appropriate; when diffusion or interfacial transport is limiting, spatial corrections or spatially derived two-compartment reductions become necessary (Hoang et al., 29 Jun 2026).