Mapping Two-Compartment Microconstants to Macroparameters

From Compartment Rates to a Concentration Curve

A two-compartment IV-bolus model describes drug moving between a central compartment, where concentration is measured, and a peripheral compartment. Its rate constants describe elimination and exchange pathways. A measured concentration-time curve, by contrast, is often expressed using two exponential terms. The mapping between these descriptions matters because the exponential rates and amplitudes are not individual compartment rate constants: each is determined by the coupled system.

Here the dose \(D\) enters the central compartment instantaneously at time zero. The central and peripheral amounts are \(X_c(t)\) and \(X_p(t)\), and the measured concentration is \(C(t)=X_c(t)/V_c\), where \(V_c\) is the central volume. The rate constants \(k_{10}\), \(k_{12}\), and \(k_{21}\) have units of inverse time: elimination from the central compartment, transfer from central to peripheral, and transfer from peripheral to central, respectively.

For positive rates, the central concentration can be written as \(C(t)=A e^{-\alpha t}+B e^{-\beta t}\), with \(\alpha>\beta>0\). The parameters \(\alpha\) and \(\beta\) are exponential rates, in inverse-time units; \(A\) and \(B\) are concentration amplitudes, in concentration units. This tutorial derives how all four are determined by the microconstants and the bolus initial conditions.

Deriving the Exponential Rates

The coupled amount equations have a two-by-two system matrix. An exponential mode with time dependence \(e^{rt}\) is possible when \(r\) satisfies the characteristic equation of that matrix. Equivalently, substituting the matrix into its characteristic determinant gives

The stable modes have \(r=-\alpha\) and \(r=-\beta\). Therefore the corresponding positive exponential rates are the roots of the polynomial obtained by changing the sign of the mode variable:

The discriminant is positive when all three rates are positive. To see why, write \(S=k_{10}+k_{12}+k_{21}\) and \(P=k_{10}k_{21}\). The discriminant also equals \((k_{10}+k_{12}-k_{21})^2+4k_{12}k_{21}\), which is strictly positive. The two roots are real and distinct. Their sum and product are positive, so both roots are positive.

A further useful ordering follows by evaluating \(f(x)=x^2-Sx+P\) at \(k_{21}\). This gives \(f(k_{21})=-k_{12}k_{21}<0\). A quadratic with positive leading coefficient is negative between its two real roots, so \(\beta<k_{21}<\alpha\). This ordering will also show that both amplitudes are positive for a positive bolus.

Deriving the IV-Bolus Amplitudes

The exponential rates describe the system's modes, but the initial conditions determine how strongly each mode contributes to central concentration. Let \(C_0=D/V_c\), the concentration immediately after the bolus. At time zero, the macroparameter form gives \(A+B=C_0\). From the central amount equation and \(X_p(0)=0\), the initial concentration slope is \(C'(0)=-(k_{10}+k_{12})C_0\). Differentiating the biexponential form gives \(C'(0)=-\alpha A-\beta B\). Consequently,

Solving these two equations gives \(A=C_0(k_{10}+k_{12}-\beta)/(\alpha-\beta)\) and \(B=C_0(\alpha-k_{10}-k_{12})/(\alpha-\beta)\). Using \(\alpha+\beta=k_{10}+k_{12}+k_{21}\), these expressions can be written directly in terms of \(k_{21}\):

This is the amplitude mapping for the specified central IV bolus. The amplitude attached to the faster exponential is \(A\), and the amplitude attached to the slower exponential is \(B\). Since \(\beta<k_{21}<\alpha\), both fractions are positive. They sum to one, so \(A+B=D/V_c\), as required by the initial concentration.

Checks, Alternatives, and Limitations

There are two complementary ways to organize the mapping. The characteristic-polynomial approach gives the exponential rates through their sum and product, then the initial concentration and slope give the amplitudes. Alternatively, one can solve the coupled differential equations directly by representing the solution as a sum of exponential modes. Both approaches give the same mapping; the polynomial-and-slope route makes the roles of the rates and bolus initial conditions especially transparent.

Several checks help identify arithmetic or modeling errors:

• Ordering: with positive exchange and elimination rates, \(\alpha>k_{21}>\beta>0\). A result that violates this ordering suggests a calculation or labeling error.
• Sum and product: verify \(\alpha+\beta=k_{10}+k_{12}+k_{21}\) and \(\alpha\beta=k_{10}k_{21}\), allowing for rounding.
• Initial concentration: verify \(A+B=D/V_c\); both sides must have concentration units.
• Initial slope: verify \(\alpha A+\beta B=(k_{10}+k_{12})D/V_c\); both sides must have concentration-per-time units.
• Units: \(k_{10}, k_{12}, k_{21}, \alpha,\beta\) must use the same inverse-time unit. The amplitudes use the same units as \(D/V_c\).

The mapping also depends on the stated structural assumptions. A peripheral elimination pathway, time-varying rates, nonlinear processes, or a non-bolus input changes the equations or initial conditions, and the formulas for this bolus model need not apply unchanged. With no exchange (\(k_{12}=0\)), one mode may have zero amplitude: the mathematical system has two rates, but a central bolus does not necessarily excite both modes. For example, if \(k_{10}=0.25\ \text{h}^{-1}\), \(k_{21}=0.60\ \text{h}^{-1}\), and \(k_{12}=0\), the roots are \(0.60\) and \(0.25\ \text{h}^{-1}\), but the amplitude for the \(0.60\ \text{h}^{-1}\) mode is zero. The central concentration is just \(C_0e^{-0.25t}\).

A common mistake is to identify \(\alpha\) with \(k_{10}+k_{12}\), or \(\beta\) with \(k_{21}\). Neither equality generally holds: the observed rates are the roots of a coupled-system polynomial. Another mistake is to calculate the roots correctly but assign the larger rate to \(B\). By convention here, \(\alpha\) is the larger, faster rate and has amplitude \(A\); \(\beta\) is the smaller, slower rate and has amplitude \(B\). The formulas map a specified microparameter set to a curve; they do not imply that the macroparameters themselves are individual physiologic rate constants.