Tutorials › Nonlinear Mixed-Effects Models (NLME) › Why Macroparameters Are Not Physiologic Rate Constants

One-, Two- and Three-Compartment Structural Models · Tutorial 41 of 797

Why Macroparameters Are Not Physiologic Rate Constants

Learn why biexponential macroparameters describe system-wide concentration modes rather than individual physiological transfer or elimination rates.

Beginner 9 min read

What You'll Learn

  • Distinguish two-compartment macroparameters from microconstants that describe individual pathways.
  • Derive how the biexponential rates relate to \(k_{10}\), \(k_{12}\), and \(k_{21}\).
  • Show why, with positive intercompartmental transfer, the terminal rate \(\beta\) is smaller than \(k_{10}\).
  • Calculate hypothetical concentration amplitudes and compare terminal and elimination half-lives.
  • Recognize common interpretation errors and checks for a misleading terminal-phase estimate.

Why the distinction matters

A biexponential concentration curve can be summarized by two rates, \(\alpha\) and \(\beta\). In a two-compartment model, those rates describe the fast and slow exponential modes seen in the central-compartment concentration. They are not, in general, the rates of individual physiological processes. In particular, the terminal slope is \(\beta\), whereas \(k_{10}\) is the first-order elimination rate constant from the central compartment. Treating these as interchangeable can give a substantially incorrect account of elimination.

The distinction is between a description of the observed curve and a description of pathways in a model. A fitted exponential rate describes how quickly a component of the concentration curve decays. A microconstant describes a specified flow process relative to the amount in a specified compartment. The two-compartment equations connect these quantities, but the connection is not an identity.

Definition: In the linear two-compartment model, \(k_{10}\) is the first-order rate constant for elimination from the central compartment, \(k_{12}\) is the first-order transfer rate constant from the central to the peripheral compartment, and \(k_{21}\) is the first-order transfer rate constant from the peripheral to the central compartment. The macroparameters \(\alpha\) and \(\beta\) are the positive exponential rates in the biexponential central concentration curve, conventionally ordered so that \(\alpha>\beta\).

How pathway rates produce exponential rates

Consider a hypothetical intravenous bolus placed in the central compartment. Let \(A_c(t)\) and \(A_p(t)\) be the amounts of drug in the central and peripheral compartments at time \(t\). Assume linear, time-invariant transfer and elimination; positive rate constants; no input after the bolus; and a measured concentration equal to central amount divided by a constant central volume \(V_c\). The amount-balance equations are

$$ \frac{dA_c}{dt}=-(k_{10}+k_{12})A_c+k_{21}A_p,\qquad \frac{dA_p}{dt}=k_{12}A_c-k_{21}A_p. $$

An exponential mode with decay rate \(\lambda\) has amounts proportional to \(e^{-\lambda t}\). Substituting this form into the two amount-balance equations gives a nonzero solution only when the determinant of the resulting coefficient matrix is zero. Equivalently,

$$ \begin{vmatrix} k_{10}+k_{12}-\lambda & -k_{21}\\ -k_{12} & k_{21}-\lambda \end{vmatrix}=0. $$

Expanding the determinant gives

$$ (k_{10}+k_{12}-\lambda)(k_{21}-\lambda)-k_{12}k_{21} =\lambda^2-(k_{10}+k_{12}+k_{21})\lambda+k_{10}k_{21}=0. $$

The two roots of this quadratic are the macroparameters \(\alpha\) and \(\beta\). Comparing the quadratic with \((\lambda-\alpha)(\lambda-\beta)\) yields the key relationships:

$$ \alpha+\beta=k_{10}+k_{12}+k_{21},\qquad \alpha\beta=k_{10}k_{21}. $$

These equations show why neither macroparameter is generally one of the microconstants. Each root depends on the combined behavior of elimination and exchange. A further useful conclusion follows directly. For positive \(k_{10}\), \(k_{12}\), and \(k_{21}\), define \(f(\lambda)=\lambda^2-(k_{10}+k_{12}+k_{21})\lambda+k_{10}k_{21}\). Then \(f(0)=k_{10}k_{21}>0\), while \(f(k_{10})=-k_{10}k_{12}<0\). Since \(f\) is an upward-opening quadratic with roots \(\beta<\alpha\), \(k_{10}\) lies between the roots:

$$ \beta<k_{10}<\alpha. $$

Thus, when both compartments are connected and the rates are positive, the terminal exponential rate \(\beta\) is strictly less than the central elimination rate \(k_{10}\). The inequality is a structural consequence of the model, not a rule that can be applied without checking the model assumptions.

Worked Example: Comparing a Terminal Rate With \(k_{10}\)

Suppose a hypothetical two-compartment model has \(k_{10}=0.20\ \mathrm{h}^{-1}\), \(k_{12}=0.60\ \mathrm{h}^{-1}\), and \(k_{21}=0.30\ \mathrm{h}^{-1}\). Find the two macroparameters and compare the slower one with \(k_{10}\).

The characteristic equation is \(\lambda^2-(k_{10}+k_{12}+k_{21})\lambda+k_{10}k_{21}=0\). Substituting the values gives

$$ \lambda^2-1.10\lambda+0.0600=0. $$

Using the quadratic formula,

$$ \alpha,\beta=\frac{1.10\pm\sqrt{1.10^2-4(0.0600)}}{2} =\frac{1.10\pm\sqrt{0.9700}}{2}. $$

Since \(\sqrt{0.9700}\approx0.9849\), the roots are \(\alpha\approx1.042\ \mathrm{h}^{-1}\) and \(\beta\approx0.05756\ \mathrm{h}^{-1}\). The terminal rate is much smaller than \(k_{10}=0.20\ \mathrm{h}^{-1}\); it is not the central elimination rate. The result also satisfies \(\beta<k_{10}<\alpha\), as expected under the stated positive-rate assumptions.

What the macroparameters mean in the concentration curve

For the bolus conditions above, central concentration has the biexponential form \(C(t)=A e^{-\alpha t}+B e^{-\beta t}\), where \(A\) and \(B\) are amplitudes. The rates specify the exponential modes; the amplitudes specify their contributions to the measured concentration. Neither amplitude is, by itself, a pathway rate constant.

The amplitudes can be checked using the initial conditions. Immediately after a bolus of dose \(D\), \(C(0)=D/V_c\). The initial central amount is \(D\), and the peripheral amount is zero, so the central concentration derivative is \(C'(0)=-(k_{10}+k_{12})D/V_c\). Equating these values to the biexponential curve and its derivative gives \(A+B=D/V_c\) and \(\alpha A+\beta B=(k_{10}+k_{12})D/V_c\). Solving these two equations yields

$$ A=\frac{D}{V_c}\frac{\alpha-k_{21}}{\alpha-\beta},\qquad B=\frac{D}{V_c}\frac{k_{21}-\beta}{\alpha-\beta}. $$

These expressions reinforce the distinction: the curve's amplitudes depend on the model's rates and initial conditions, while the exponents are the system's decay modes. A good fit of the curve does not make an exponent a direct measurement of one physiological pathway.

Worked Example: Checking the Curve Amplitudes

Continue with the hypothetical rate constants in the previous example, and suppose \(D=100\ \mathrm{mg}\) and \(V_c=20\ \mathrm{L}\). Calculate the amplitudes and verify the initial concentration and initial decline.

The initial concentration is \(D/V_c=100/20=5.00\ \mathrm{mg/L}\). Using \(\alpha\approx1.04244\ \mathrm{h}^{-1}\), \(\beta\approx0.05756\ \mathrm{h}^{-1}\), and \(k_{21}=0.30\ \mathrm{h}^{-1}\),

$$ A=5.00\frac{1.04244-0.300}{1.04244-0.05756} \approx3.769\ \mathrm{mg/L}, $$
$$ B=5.00\frac{0.300-0.05756}{1.04244-0.05756} \approx1.231\ \mathrm{mg/L}. $$

The amplitudes sum to \(3.769+1.231=5.000\ \mathrm{mg/L}\), matching \(C(0)\). The initial decline calculated from the curve is

$$ C'(0)=-\alpha A-\beta B \approx-(1.04244)(3.769)-(0.05756)(1.231) \approx-4.00\ \mathrm{mg/L/h}. $$

From the amount-balance equation, the same derivative is \(-(k_{10}+k_{12})D/V_c=-(0.20+0.60)(5.00)=-4.00\ \mathrm{mg/L/h}\). This agreement checks the arithmetic and shows that the initial decline reflects both elimination and distribution out of the central compartment, not \(k_{10}\) alone.

Terminal slope and half-life

At sufficiently late times, the faster exponential term becomes small relative to the slower term, provided its amplitude is nonzero. The concentration then approaches \(B e^{-\beta t}\), and the terminal log-concentration slope is \(-\beta\). The corresponding terminal half-life is \(\ln(2)/\beta\). It describes the decline of the terminal concentration mode; it is not generally the half-life associated with central elimination, \(\ln(2)/k_{10}\).

Worked Example: Same \(k_{10}\), Different Terminal Half-Life

Compare two hypothetical models with the same \(k_{10}=0.20\ \mathrm{h}^{-1}\) and \(k_{21}=0.30\ \mathrm{h}^{-1}\), but different central-to-peripheral transfer rates. Model 1 has \(k_{12}=0.60\ \mathrm{h}^{-1}\), as above; Model 2 has \(k_{12}=0.10\ \mathrm{h}^{-1}\).

For Model 1, the terminal rate is \(\beta\approx0.05756\ \mathrm{h}^{-1}\), so its terminal half-life is

$$ t_{1/2,\beta}=\frac{\ln(2)}{\beta} =\frac{0.6931}{0.05756}\approx12.04\ \mathrm{h}. $$

For Model 2, the sum of the rates in the characteristic equation is \(0.20+0.10+0.30=0.60\ \mathrm{h}^{-1}\), and the product term is \(0.20(0.30)=0.0600\ \mathrm{h}^{-2}\). Therefore,

$$ \lambda^2-0.60\lambda+0.0600=0,\qquad \beta=\frac{0.60-\sqrt{0.60^2-4(0.0600)}}{2} \approx0.1268\ \mathrm{h}^{-1}. $$

Its terminal half-life is \(0.6931/0.1268\approx5.47\ \mathrm{h}\). The two models share the same \(k_{10}\), but their terminal rates and terminal half-lives differ because their distribution rates differ. For comparison, the half-life calculated directly from \(k_{10}=0.20\ \mathrm{h}^{-1}\) is \(0.6931/0.20=3.47\ \mathrm{h}\); that is not either model's terminal half-life.

Interpretation, diagnostics, and alternatives

A fitted biexponential curve is useful for describing concentration-time behavior, including its fast and slow phases. A compartmental model is needed when the question concerns specific pathways such as elimination from the central compartment or exchange between compartments. Interpreting \(\beta\) as \(k_{10}\) silently replaces a model-wide decay mode with a pathway-specific rate, and the numerical comparison above shows why that substitution can fail.

When interpreting a terminal slope, check that the observations actually support a terminal phase: late concentrations should be consistent with an approximately straight segment on a log-concentration plot, and the fitted slope should not be driven by a single point or dominated by measurement noise. Curvature in the proposed terminal region, sensitivity of the slope to which late observations are included, or a fitted rate that changes markedly with the selected time window can indicate that the terminal rate is poorly estimated or that a single exponential is not an adequate description there. These checks assess the fitted curve; they do not establish that the fitted rate equals a physiological microconstant.

The main alternative is to report the biexponential parameters as empirical curve descriptors, without assigning them pathway meanings. If the scientific question is about elimination and intercompartmental exchange, use the compartmental parameterization and its stated assumptions, rather than relabeling the macroparameters. Neither approach makes model assumptions disappear: a biexponential fit does not by itself establish a two-compartment mechanism, and pathway interpretations are conditional on the compartment model being appropriate.

Common mistake

Pitfall: “The terminal log-concentration slope is \(-0.08\ \mathrm{h}^{-1}\), so the physiological elimination rate constant is \(k_{10}=0.08\ \mathrm{h}^{-1}\).” This reasoning equates the terminal macroparameter \(\beta\) with a microconstant. In the positive-rate two-compartment model, \(\beta<k_{10}<\alpha\), so the terminal rate is not \(k_{10}\). Report the fitted terminal slope as \(\beta\), and use the compartmental model to estimate and interpret \(k_{10}\).

Check Your Understanding

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

  1. For positive \(k_{10}\), \(k_{12}\), and \(k_{21}\), where does \(k_{10}\) lie relative to \(\beta\) and \(\alpha\)?
  2. A hypothetical curve has terminal rate \(\beta=0.08\ \mathrm{h}^{-1}\). Is this enough to conclude that \(k_{10}=0.08\ \mathrm{h}^{-1}\)? Why or why not?
  3. For Model 2 in the worked example, what is the terminal half-life to three significant figures?
  4. What does a visibly curved late segment on a log-concentration plot suggest about interpreting one fitted terminal slope?

Solutions. (1) \(k_{10}\) lies strictly between the macroparameters: \(\beta<k_{10}<\alpha\). (2) No. The terminal rate is a macroparameter combining the model's elimination and distribution behavior; it is not generally the central elimination microconstant. (3) Using \(\beta\approx0.1268\ \mathrm{h}^{-1}\), \(t_{1/2,\beta}=0.6931/0.1268\approx5.47\ \mathrm{h}\). (4) The proposed terminal region may not yet follow a single exponential, or the fitted slope may depend on the chosen interval. The terminal rate should not be interpreted physiologically without checking the fit and the compartmental assumptions.