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.
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\).
Solving the first equation for one amplitude, \(B=C_0-A\), and substituting into the slope equation yields
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:
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:
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\).
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
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:
The rate difference is \(\alpha-\beta=1.284-0.01557=1.26843\ \mathrm{h^{-1}}\). Using the displayed rates in both amplitude calculations,
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,
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:
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.
Use the central concentration after a bolus and identify \(\alpha\) as the larger rate and \(\beta\) as the smaller rate.
Calculate \(r_0=-s_0/C_0\), checking that it has inverse-time units.
For positive amplitudes, verify \(\beta<r_0<\alpha\).
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.
Check Your Understanding
Answer each question, then compare with the solutions below.
- 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\).
- Why does the initial central concentration slope involve both \(k_{10}\) and \(k_{12}\) after a bolus?
- If \(r_0=\beta\), what does the amplitude formula imply for \(A\), and does this satisfy the strict positive-amplitude condition?
- 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.