A Worked Two-Compartment Parameter Conversion

From Clearance and Volumes to a Concentration Curve

A two-compartment model can be specified using clearances and volumes, or using rate constants that describe transfers between compartments. For a central intravenous bolus, converting between these descriptions gives a direct route from interpretable pharmacokinetic parameters to a biexponential concentration-time curve. The important checks are not only that the rate constants are calculated correctly, but also that the curve starts at the concentration implied by the dose and central volume.

Consider the linear mammillary model: drug is administered to the central compartment, can distribute between central and peripheral compartments, and is eliminated from the central compartment. The calculation below assumes a single instantaneous IV bolus, constant parameters, and first-order linear transfer and elimination. There is no absorption process in this model.

The Conversion and Its Checks

The microconstants follow from clearance divided by the volume of the compartment from which the corresponding flux occurs. The amounts \(A_c\) and \(A_p\) in the central and peripheral compartments obey

The two decay rates are the positive eigenvalues of the negative of this amount-balance system. Equivalently, they are the roots of the characteristic equation

Thus the sum of the decay rates is \(k_{10}+k_{12}+k_{21}\), and their product is \(k_{10}k_{21}\). These are macro-rates for the observed biexponential curve; they are not the individual transfer or elimination microconstants. This distinction is useful when checking a conversion: the rates must satisfy both the sum and product relations.

For a bolus, \(A_c(0)=D\) and \(A_p(0)=0\), so the initial central concentration is \(C(0)=D/V_c\). The initial slope follows directly from the central amount balance: \(C'(0)=-(k_{10}+k_{12})D/V_c\). The amplitudes must match both initial conditions.

A useful sign check follows from the same equation. Write \(s=k_{10}+k_{12}\). The quadratic evaluated at \(s\) is \(-k_{12}k_{21}\), which is negative when all three microconstants are positive. Since the quadratic opens upward, \(s\) lies between its two roots: \(\beta<s<\alpha\). Therefore both amplitude fractions in the result are positive. A negative amplitude in this positive-parameter, central-bolus setting is a signal to recheck the roots, labels, or arithmetic.

Why the Amplitude Calculation Works

The amplitude formulas are obtained by solving two equations, rather than by treating either exponential coefficient as a rate constant. At \(t=0\), the curve must equal \(C_0\). Differentiating the curve gives \(C'(t)=-\alpha A e^{-\alpha t}-\beta B e^{-\beta t}\), so its initial slope must also agree with the compartment model.

Substitute \(B=C_0-A\) into the second equation. Then \((\alpha-\beta)A+\beta C_0=(k_{10}+k_{12})C_0\), so

Subtracting this from \(C_0\) gives \(B=C_0[\alpha-(k_{10}+k_{12})]/(\alpha-\beta)\). This also explains why the two initial-condition checks are informative: a curve can have the right starting concentration but the wrong initial slope if its amplitudes or decay rates were combined incorrectly.

Practical Checks, Alternatives, and Limitations

A reliable hand calculation uses several independent checks rather than relying on a single plausible-looking curve. Confirm that the quadratic’s roots add to \(k_{10}+k_{12}+k_{21}\) and multiply to \(k_{10}k_{21}\). Confirm that \(A+B=D/V_c\), and that \(-\alpha A-\beta B=-(k_{10}+k_{12})D/V_c\). For positive microconstants, both roots and both amplitudes should be positive, with \(\alpha>\beta\).

The direct alternative is to work with the amount-balance differential equations and solve them as a two-state linear system. This is mathematically equivalent: the biexponential rates are determined by the system’s eigenvalues, and the amplitudes are fixed by its initial conditions. For a specified set of \(CL,Q,V_c,V_p\), the clearance-volume route is usually efficient because it makes each microconstant’s source parameter explicit. When checking an independently supplied biexponential curve instead, its sum and product relations can be compared with the microconstants. Do not equate \(\alpha\) or \(\beta\) with a single physiologic transfer rate; each is a system-level decay rate.

This conversion relies on the stated linear, time-invariant, two-compartment structure and an IV bolus into the central compartment. Different input patterns or model structures require different initial conditions or equations. The calculation itself does not establish that this structure is an adequate description of observed data; it converts parameters under the assumed model. Keep enough precision through the roots and amplitude calculations, then round the reported values. Rounding an intermediate discriminant or root too early can make the amplitudes fail the initial checks.

• Wrong: “The slow macro-rate \(\beta\) is the return rate \(k_{21}\).” Correction: \(\beta\) is an eigenvalue of the coupled system; calculate \(k_{21}=Q/V_p\) directly.
• Wrong: “The curve starts at \(A\), so \(A=D/V_c\).” Correction: The two terms are both present at time zero, so \(A+B=D/V_c\).
• Wrong: “A negative amplitude is acceptable because the curve is biexponential.” Correction: With positive parameters in this central-bolus model, the roots straddle \(k_{10}+k_{12}\), making both amplitudes positive. Recheck arithmetic and root ordering.
• Wrong: “The discriminant is small, so its square root can be approximated by the nearest convenient value.” Correction: Evaluate the discriminant explicitly and verify the resulting roots using their sum and product.