Tutorials › Nonlinear Mixed-Effects Models (NLME) › Two-Compartment IV Bolus Amplitudes From Initial Conditions

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

Two-Compartment IV Bolus Amplitudes From Initial Conditions

Learn how the bolus concentration at time zero and its initial slope determine the amplitudes of a two-compartment biexponential concentration profile.

Beginner 10 min read

What You'll Learn

  • Define the initial concentration and initial slope for a central IV bolus profile.
  • Derive the two biexponential amplitudes from those initial conditions.
  • Check whether the inferred amplitudes are positive and consistent with the stated rate constants.
  • Calculate amplitudes from hypothetical concentration, slope, and rate data.
  • Recognize when noisy or incompatible initial-condition information makes the calculation unreliable.

Why Initial Conditions Determine the Amplitudes

A two-compartment IV bolus concentration profile has a fast distribution phase and a slower terminal phase. Its two exponential rate constants determine how quickly those phases decline, but they do not by themselves specify the concentration contributed by each phase. The amplitudes are determined by the concentration immediately after the bolus and the initial rate of change of that concentration.

This provides a direct way to connect the biexponential curve to the central compartment at time zero. It also gives a useful consistency check: amplitudes obtained from a standard two-compartment model with positive parameters should be positive and should add to the initial concentration.

Definition: For a central IV bolus, write the concentration profile as \(C(t)=A e^{-\alpha t}+B e^{-\beta t}\), where \(\alpha>\beta>0\). The initial concentration is \(C_0=C(0)\), and the initial slope is \(s_0=\left.\frac{dC(t)}{dt}\right|_{t=0}\). The amplitudes \(A\) and \(B\) have concentration units; \(\alpha\) and \(\beta\) have inverse-time units.

Here \(\alpha\) is the fast rate constant and \(\beta\) the slow rate constant. The terms “fast” and “slow” refer to how rapidly each exponential declines, not to the size of its amplitude. For a positive dose and positive compartment volumes and rate constants, the standard two-compartment model gives positive amplitudes.

Deriving the Amplitude Formulas

At time zero, both exponential terms equal one. Differentiating the profile gives its slope at any time, and evaluating that derivative at zero supplies a second equation. Together, the two initial conditions are sufficient to solve for \(A\) and \(B\).

$$ C_0=A+B,\qquad s_0=-\alpha A-\beta B $$

Solving the first equation for one amplitude, \(B=C_0-A\), and substituting into the slope equation yields

$$ s_0=-\alpha A-\beta(C_0-A) =(\beta-\alpha)A-\beta C_0. $$

Rearranging gives \(A=(s_0+\beta C_0)/(\beta-\alpha)\). Equivalently, define the positive initial fractional decline rate \(r_0=-s_0/C_0\). For a positive initial concentration, this converts the slope equation to \(\alpha A+\beta B=r_0 C_0\). Solving the two linear equations gives the more interpretable formulas:

$$ A=C_0\frac{r_0-\beta}{\alpha-\beta}, \qquad B=C_0\frac{\alpha-r_0}{\alpha-\beta}, \qquad r_0=-\frac{s_0}{C_0}. $$

These formulas are derived from the biexponential form and its initial conditions; they do not require a numerical curve-fitting algorithm. The amplitudes sum to \(C_0\), and their rate-weighted sum equals the magnitude of the initial slope:

$$ A+B=C_0,\qquad \alpha A+\beta B=-s_0=r_0C_0. $$

For both amplitudes to be positive, the initial fractional decline rate must lie between the two exponential rates. This follows directly from the formulas: because \(C_0>0\) and \(\alpha-\beta>0\), \(A>0\) requires \(r_0>\beta\), while \(B>0\) requires \(r_0<\alpha\).

Key condition: In a standard two-compartment IV bolus model with positive parameters, \(\beta<r_0<\alpha\). An initial fractional decline rate outside this interval produces a negative amplitude and signals that the supplied values are incompatible with the assumed positive-amplitude model.

Connection to the Central Compartment

The initial slope has a specific compartmental interpretation. Let \(A_c(t)\) and \(A_p(t)\) be the central and peripheral amounts, and let \(V_c\) be the central volume. For a central IV bolus, the peripheral compartment initially contains no drug. With first-order elimination from the central compartment and first-order exchange between compartments, the central amount balance is

$$ \left.\frac{dA_c}{dt}\right|_{t=0} =-(k_{10}+k_{12})A_c(0)+k_{21}A_p(0) =-(k_{10}+k_{12})A_c(0). $$

Since \(C(t)=A_c(t)/V_c\), the initial concentration is \(C_0=A_c(0)/V_c\), and the corresponding initial slope is \(s_0=-(k_{10}+k_{12})C_0\). Therefore \(r_0=k_{10}+k_{12}\). The initial decline reflects both elimination from the central compartment and distribution out of it; it is not simply the elimination rate constant \(k_{10}\).

Worked Calculations

Worked Example: Calculating Both Amplitudes From the Initial Slope

Suppose a hypothetical central IV bolus profile has \(C_0=10.0\ \mathrm{mg/L}\), initial slope \(s_0=-12.0\ \mathrm{mg/L/h}\), fast rate \(\alpha=1.284\ \mathrm{h^{-1}}\), and slow rate \(\beta=0.01557\ \mathrm{h^{-1}}\). Find the amplitudes.

First calculate the initial fractional decline rate:

$$ r_0=-\frac{-12.0\ \mathrm{mg/L/h}}{10.0\ \mathrm{mg/L}} =1.20\ \mathrm{h^{-1}}. $$

The rate difference is \(\alpha-\beta=1.284-0.01557=1.26843\ \mathrm{h^{-1}}\). Using the displayed rates in both amplitude calculations,

$$ A=10.0\frac{1.20-0.01557}{1.26843} =\frac{11.8443}{1.26843} \approx 9.3378\ \mathrm{mg/L}, $$
$$ B=10.0\frac{1.284-1.20}{1.26843} =\frac{0.8400}{1.26843} \approx 0.6622\ \mathrm{mg/L}. $$

The checks agree with the inputs: \(A+B\approx9.3378+0.6622=10.0000\ \mathrm{mg/L}\). Also, \(\alpha A+\beta B\approx1.284(9.3378)+0.01557(0.6622)\approx12.000\ \mathrm{mg/L/h}\), so the initial slope is approximately \(-12.0\ \mathrm{mg/L/h}\). Small discrepancies in the final digits are due to rounding.

Worked Example: Checking the Result Against the Compartmental Slope

For a second hypothetical profile, suppose \(C_0=8.00\ \mathrm{mg/L}\), \(k_{10}=0.15\ \mathrm{h^{-1}}\), and \(k_{12}=0.60\ \mathrm{h^{-1}}\). A fitted biexponential has \(\alpha=0.95\ \mathrm{h^{-1}}\) and \(\beta=0.05\ \mathrm{h^{-1}}\). Calculate the initial slope and amplitudes, then check their consistency.

At time zero the peripheral amount is zero, so \(r_0=k_{10}+k_{12}=0.75\ \mathrm{h^{-1}}\), and \(s_0=-r_0C_0=-0.75(8.00)=-6.00\ \mathrm{mg/L/h}\). The rate difference is \(0.95-0.05=0.90\ \mathrm{h^{-1}}\). Therefore,

$$ A=8.00\frac{0.75-0.05}{0.90} =8.00\frac{0.70}{0.90} \approx6.222\ \mathrm{mg/L}, $$
$$ B=8.00\frac{0.95-0.75}{0.90} =8.00\frac{0.20}{0.90} \approx1.778\ \mathrm{mg/L}. $$

The amplitudes sum to \(6.222+1.778=8.000\ \mathrm{mg/L}\). Using the unrounded amplitude values, the slope check is \(-0.95(6.222\ldots)-0.05(1.777\ldots)=-6.00\ \mathrm{mg/L/h}\). The rate condition also holds: \(0.05<0.75<0.95\ \mathrm{h^{-1}}\), so both amplitudes are positive.

Worked Example: Detecting Incompatible Initial Information

Suppose hypothetical inputs are \(C_0=10.0\ \mathrm{mg/L}\), \(s_0=-9.0\ \mathrm{mg/L/h}\), \(\alpha=0.80\ \mathrm{h^{-1}}\), and \(\beta=0.10\ \mathrm{h^{-1}}\). Determine whether they yield positive amplitudes under the standard model.

The implied fractional decline rate is \(r_0=9.0/10.0=0.90\ \mathrm{h^{-1}}\), which is greater than the fast rate \(0.80\ \mathrm{h^{-1}}\). Substitution makes the incompatibility explicit:

$$ A=10.0\frac{0.90-0.10}{0.80-0.10} =\frac{8.0}{0.70} \approx11.43\ \mathrm{mg/L}, \qquad B=10.0\frac{0.80-0.90}{0.80-0.10} =\frac{-1.0}{0.70} \approx-1.429\ \mathrm{mg/L}. $$

The amplitudes still sum to \(10.0\ \mathrm{mg/L}\), and their weighted sum reproduces the supplied slope algebraically. But the negative slow amplitude is inconsistent with the standard positive-parameter two-compartment IV bolus profile. The correct conclusion is not to replace the negative value with zero: the inputs, or the assumed model, need to be checked.

Using and Checking Initial-Condition Estimates

In an ideal mathematical model, \(C_0\) and \(s_0\) are exact values. In data analysis, however, the slope is usually not observed directly. It may be estimated from very early measurements or inferred from an estimated compartmental model. A difference between two measured concentrations is not automatically the slope at time zero: it is a finite-interval average change, and the biexponential curve is generally changing slope over that interval.

The initial-condition calculation is most useful when the quantities refer to the same central concentration profile and use consistent time and concentration units. For example, if the rate constants are in \(\mathrm{h^{-1}}\), the slope must be in concentration per hour. An immediate post-bolus concentration is also a model-based value at \(t=0\); an actual sample taken later may already reflect distribution.

1
Specify the profile.
Use the central concentration after a bolus and identify \(\alpha\) as the larger rate and \(\beta\) as the smaller rate.
2
Form the initial decline rate.
Calculate \(r_0=-s_0/C_0\), checking that it has inverse-time units.
3
Check compatibility.
For positive amplitudes, verify \(\beta<r_0<\alpha\).
4
Calculate and verify.
Use the amplitude formulas, then check \(A+B=C_0\) and \(\alpha A+\beta B=-s_0\), allowing for rounding.

There are alternatives to calculating amplitudes from initial conditions. A biexponential curve can be fitted directly, estimating amplitudes and rates together; this uses the full set of observations rather than just the initial concentration and slope. Alternatively, amplitudes can be derived from a compartmental parameterization. The initial-condition approach is an algebraic relationship and a consistency check, not a substitute for fitting data when the initial slope is uncertain.

Each approach has limitations. Direct fitting can be sensitive to sparse or noisy data, while an initial-condition calculation can be especially sensitive to slope error because it uses a local rate at time zero. The algebra also becomes numerically unstable when \(\alpha\) and \(\beta\) are very close, since the denominator \(\alpha-\beta\) is then small. A negative calculated amplitude can indicate measurement noise, mismatched units, imprecise early-time information, or a mismatch between the assumed two-compartment model and the data. It should prompt investigation rather than be interpreted as a positive compartmental contribution.

Common Mistakes

  • Using a positive slope sign. A declining concentration has \(s_0<0\). Define \(r_0=-s_0/C_0\) to use a positive decline rate, and retain the negative sign in \(s_0=-r_0C_0\).
  • Confusing \(r_0\) with \(k_{10}\). Initially, distribution out of the central compartment also contributes to the decline, so \(r_0=k_{10}+k_{12}\), not just \(k_{10}\).
  • Switching the fast and slow rates. Keep \(\alpha>\beta\); the amplitude paired with \(\alpha\) is \(A\), and the one paired with \(\beta\) is \(B\).
  • Rounding inconsistently. Substituting rounded rates may not reproduce amplitudes calculated from higher-precision rates. Use the same displayed inputs throughout and report an appropriate number of significant figures.
  • Treating any early finite difference as the initial slope. A measurement interval is not generally infinitesimal. State how the slope was obtained and recognize that its uncertainty affects the amplitudes.
Key takeaway: For \(C(t)=A e^{-\alpha t}+B e^{-\beta t}\), the initial concentration and slope imply \(A+B=C_0\) and \(\alpha A+\beta B=-s_0\). With \(r_0=-s_0/C_0\), the amplitudes are \(A=C_0(r_0-\beta)/(\alpha-\beta)\) and \(B=C_0(\alpha-r_0)/(\alpha-\beta)\). Positive amplitudes require \(\beta<r_0<\alpha\).

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. For \(C_0=12.0\ \mathrm{mg/L}\), \(s_0=-6.0\ \mathrm{mg/L/h}\), \(\alpha=0.70\ \mathrm{h^{-1}}\), and \(\beta=0.10\ \mathrm{h^{-1}}\), calculate \(r_0\), \(A\), and \(B\).
  2. Why does the initial central concentration slope involve both \(k_{10}\) and \(k_{12}\) after a bolus?
  3. If \(r_0=\beta\), what does the amplitude formula imply for \(A\), and does this satisfy the strict positive-amplitude condition?
  4. For \(C_0=5.0\ \mathrm{mg/L}\), \(r_0=0.40\ \mathrm{h^{-1}}\), \(\alpha=0.60\ \mathrm{h^{-1}}\), and \(\beta=0.05\ \mathrm{h^{-1}}\), calculate the amplitudes and verify their sum.

Solutions. (1) \(r_0=6.0/12.0=0.50\ \mathrm{h^{-1}}\). The denominator is \(0.70-0.10=0.60\ \mathrm{h^{-1}}\), so \(A=12.0(0.50-0.10)/0.60=8.00\ \mathrm{mg/L}\) and \(B=12.0(0.70-0.50)/0.60=4.00\ \mathrm{mg/L}\). (2) At time zero the peripheral amount is zero, but central drug is both eliminated at rate \(k_{10}\) and distributed out at rate \(k_{12}\); thus \(s_0=-(k_{10}+k_{12})C_0\). (3) If \(r_0=\beta\), then \(A=C_0(\beta-\beta)/(\alpha-\beta)=0\), so it does not meet the strict condition for two positive amplitudes. (4) \(A=5.0(0.40-0.05)/(0.60-0.05)=3.182\ \mathrm{mg/L}\), and \(B=5.0(0.60-0.40)/(0.60-0.05)=1.818\ \mathrm{mg/L}\). Their sum is \(5.000\ \mathrm{mg/L}\); also \(0.05<0.40<0.60\ \mathrm{h^{-1}}\), so both are positive.