Tutorials › Nonlinear Mixed-Effects Models (NLME) › One-Compartment Constant-Rate Infusion

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

One-Compartment Constant-Rate Infusion

Derive and apply the concentration-time equations for constant-rate infusion, including the finite-time approach to steady state and the decline after infusion stops.

Beginner 10 min read

What You'll Learn

  • State the assumptions behind the one-compartment constant-rate infusion model.
  • Derive concentration during infusion from the amount-balance equation.
  • Calculate the concentration at infusion end and after infusion stops.
  • Interpret steady-state concentration as an asymptote and calculate the approach to it.
  • Check units, limiting behavior, and common signs of model mismatch.

Why an Infusion Produces a Rising Concentration

A constant-rate intravenous infusion supplies drug continuously, while elimination removes drug continuously. At the start, the amount entering the body exceeds the amount eliminated, so concentration rises. As concentration increases, the elimination rate increases too. The rise slows as input and elimination approach balance. When the infusion stops, input becomes zero and concentration declines by elimination.

The equations below describe the structural concentration profile for one individual under a one-compartment model. They assume a well-mixed compartment, so the measured concentration is the amount in the compartment divided by its apparent volume. They do not, by themselves, describe between-person differences or measurement error.

Definition: Let \(A(t)\) be the amount of drug in the compartment at time \(t\), \(C(t)=A(t)/V\) its concentration, \(V\) the apparent volume of distribution, \(CL\) clearance, and \(k=CL/V\) the first-order elimination rate constant. Let \(R_0\) be a constant infusion rate, and let \(T\) be the time the infusion stops. During infusion, the amount-balance equation is \(dA/dt=R_0-kA\). After infusion stops, it is \(dA/dt=-kA\).

The model assumes \(V>0\), \(CL>0\), and therefore \(k>0\), with these parameters constant over the period considered. Elimination is first-order: the elimination rate is proportional to the amount present. The drug enters directly into the compartment at a constant rate from time 0 to \(T\), and there is no additional input during that interval. Concentration is continuous when the infusion stops; only the input rate changes.

Deriving the Concentration During Infusion

Because \(A=VC\) and \(V\) is constant, dividing the amount-balance equation by \(V\) gives the concentration equation. Since \(k=CL/V\), the input term \(R_0/V\) can also be written \(kR_0/CL\).

$$ \frac{dC}{dt}=\frac{R_0}{V}-kC $$

For an infusion beginning with concentration \(C(0)=C_0\), multiply by the integrating factor \(e^{kt}\). The product rule gives \(d(e^{kt}C)/dt=(R_0/V)e^{kt}\). Integrating from 0 to \(t\), then solving for \(C(t)\), yields:

$$ e^{kt}C(t)-C_0=\frac{R_0}{V}\frac{e^{kt}-1}{k} \quad\Longrightarrow\quad C(t)=C_0e^{-kt}+\frac{R_0}{CL}\left(1-e^{-kt}\right), \qquad 0\leq t\leq T. $$

For an infusion into a compartment with no drug initially, \(C_0=0\), so the first term vanishes. The concentration rises toward \(R_0/CL\), but for finite \(t\) it remains below that value when the infusion begins at zero. At \(t=0\), the equation gives \(C(0)=0\), as required by the initial condition.

Steady-state concentration: During an indefinitely continued constant-rate infusion, steady state is the equilibrium at which input equals elimination. Setting \(dC/dt=0\) gives \(C_{\mathrm{ss}}=R_0/CL\). This is an asymptote of the concentration curve: with a positive elimination rate constant, the model approaches it but does not reach it at any finite time when starting below it.

The units provide a useful check. If \(R_0\) is in mg/h and \(CL\) is in L/h, \(R_0/CL\) is in mg/L. The exponent \(kt\) is dimensionless when \(k\) is in h\(^{-1}\) and time is in hours.

Worked Example: Concentration During a Six-Hour Infusion

Suppose a hypothetical patient receives an infusion at \(R_0=30\) mg/h. Let \(CL=6\) L/h and \(V=60\) L, and assume there is no drug in the compartment at the start. Calculate the concentration at 3 hours and at the end of a 6-hour infusion.

First, calculate the elimination rate constant and the steady-state concentration:

$$ k=\frac{CL}{V}=\frac{6\ \mathrm{L/h}}{60\ \mathrm{L}}=0.1000\ \mathrm{h}^{-1}, \qquad C_{\mathrm{ss}}=\frac{R_0}{CL}=\frac{30\ \mathrm{mg/h}}{6\ \mathrm{L/h}}=5.000\ \mathrm{mg/L}. $$

With \(C_0=0\), use \(C(t)=C_{\mathrm{ss}}(1-e^{-kt})\). At 3 hours:

$$ C(3)=5.000\left(1-e^{-(0.1000)(3)}\right) =5.000(1-e^{-0.3000}) =5.000(1-0.740818) \approx 1.296\ \mathrm{mg/L}. $$

At 6 hours:

$$ C(6)=5.000\left(1-e^{-(0.1000)(6)}\right) =5.000(1-e^{-0.6000}) =5.000(1-0.548812) \approx 2.256\ \mathrm{mg/L}. $$

The 6-hour concentration is below the \(5.000\) mg/L steady-state value, as expected: the infusion has not run long enough for the concentration to approach that asymptote closely.

After the Infusion Stops

At \(t=T\), the amount already in the compartment does not disappear when the infusion stops. For \(t>T\), set the input rate to zero and solve \(dC/dt=-kC\) using the concentration at the stop time as the initial value. This gives the post-infusion profile:

$$ C(t)=C(T)e^{-k(t-T)}, \qquad t\geq T. $$

The concentration at infusion end is found from the during-infusion equation, using the initial concentration at \(t=0\). For an initially drug-free compartment, \(C(T)=C_{\mathrm{ss}}(1-e^{-kT})\). The complete profile is therefore the rising infusion curve up to \(T\), followed by exponential decline from that endpoint. The expression after stopping uses elapsed time since stopping, \(t-T\), not time since the infusion began.

Worked Example: Concentration Three Hours After Stopping

Continue the hypothetical infusion above: \(R_0=30\) mg/h, \(CL=6\) L/h, \(V=60\) L, and \(T=6\) h. Find the concentration at \(t=9\) h, three hours after the infusion stops.

The endpoint concentration is \(C(6)\approx2.25594\) mg/L, and \(k=0.1000\) h\(^{-1}\). Apply the post-infusion equation with \(t-T=9-6=3\) h:

$$ C(9)=C(6)e^{-0.1000(9-6)} =2.25594e^{-0.3000} =2.25594(0.740818) \approx 1.671\ \mathrm{mg/L}. $$

The concentration falls from about \(2.256\) mg/L at infusion end to about \(1.671\) mg/L three hours later. The decline is not calculated by subtracting an amount at a constant rate; first-order elimination produces an exponential decline.

Steady State, Duration, and Half-Life

The infusion rate and clearance determine the steady-state concentration, while the elimination rate constant determines how rapidly the profile approaches it. For an initially drug-free compartment, the fraction of steady state reached at time \(t\) is \(1-e^{-kt}\). The half-life is \(t_{1/2}=\ln(2)/k\). After one half-life of infusion, the concentration is half of \(C_{\mathrm{ss}}\); after each additional half-life, the remaining gap between the current concentration and \(C_{\mathrm{ss}}\) is halved.

Worked Example: Three Half-Lives of Infusion

For a hypothetical constant infusion, let \(R_0=20\) mg/h and \(CL=5\) L/h, so \(C_{\mathrm{ss}}=4.000\) mg/L. Suppose the elimination half-life is 8.00 h, and the compartment is initially drug-free. Find the concentration after 24.0 h of infusion and the concentration 8.00 h after the infusion then stops.

Twenty-four hours is three half-lives. Thus \(e^{-k(24)}=(1/2)^3=1/8\). The end-of-infusion concentration is:

$$ C(24)=4.000\left(1-\frac{1}{8}\right) =4.000\left(\frac{7}{8}\right) =3.500\ \mathrm{mg/L}. $$

After stopping, one more half-life elapses, so the concentration halves:

$$ C(32)=C(24)e^{-k(8)} =3.500\left(\frac{1}{2}\right) =1.750\ \mathrm{mg/L}. $$

After three half-lives of infusion, the concentration is 87.5% of steady state, not exactly at steady state. The post-infusion calculation begins from the actual endpoint concentration, not from \(C_{\mathrm{ss}}\).

Assumptions, Checks, and Alternatives

The formulas are exact for the stated model, but the assumptions matter in application. The infusion rate must be constant over the interval being modeled, and the relevant \(CL\) and \(V\) must be constant. The one-compartment description assumes rapid mixing relative to the time scale of interest. The exponential elimination equation assumes first-order elimination, rather than a rate that saturates as concentration rises. If there is drug present before infusion starts, use the general solution with the actual \(C_0\); do not use the drug-free formula.

Several checks can reveal arithmetic or model problems:

  • During infusion from \(C_0=0\), concentration should rise monotonically toward \(R_0/CL\), because \(dC/dt=(R_0/V)e^{-kt}>0\).
  • At \(t=0\), the drug-free formula gives zero; at infusion stop, the during- and post-infusion equations give the same value \(C(T)\).
  • After stopping, concentration should decrease for \(k>0\), and after one half-life it should be half the value at the start of that interval.
  • A negative exponent with time in hours requires \(k\) in h\(^{-1}\). A concentration unit error often signals that the units of infusion rate and clearance do not match.
  • A persistent mismatch between observed profiles and the predicted shape may indicate that the constant-input, first-order, or one-compartment assumptions are inadequate; it is not fixed merely by changing the arithmetic.

The closed-form equation is the direct approach when the input is constant and the one-compartment assumptions are appropriate. For an infusion whose rate changes over time, one can solve the amount-balance equation separately over each constant-rate interval, carrying each interval's ending concentration forward as the next initial condition. More general time-varying inputs or more complex structural models can require numerical solution of differential equations. A bolus input is a different input pattern: it changes the initial amount rather than supplying drug continuously. These approaches are alternatives for different input schedules or model structures, not interchangeable formulas for the same dosing history.

The model has important limits. A good fit to this profile does not establish that the body is truly one well-mixed compartment, and a steady-state calculation is not a guarantee that a finite infusion will reach the target concentration. In practice, an infusion may stop before the asymptote is approached. The equations also describe the structural concentration, not assay variability or all sources of patient-to-patient variability.

Common pitfall: Do not use \(R_0/CL\) as the concentration at every time during infusion. It is the steady-state concentration; starting from zero, the finite-time concentration is \( (R_0/CL)(1-e^{-kt})\). Also do not continue the infusion equation after it stops: for \(t>T\), the concentration is \(C(T)e^{-k(t-T)}\).

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. A drug is infused at 12 mg/h into a one-compartment model with clearance 3 L/h. What is the steady-state concentration?
  2. Starting from zero, what fraction of steady state is reached after two half-lives of infusion?
  3. An infusion stops at \(T=5\) h with \(C(T)=3.0\) mg/L. If \(k=0.20\) h\(^{-1}\), what is the concentration 5 h later?
  4. Why is it incorrect to calculate a finite-time infusion concentration as \(R_0/CL\)?

Solutions. 1. \(C_{\mathrm{ss}}=R_0/CL=(12\ \mathrm{mg/h})/(3\ \mathrm{L/h})=4.0\) mg/L. 2. After two half-lives, the fraction is \(1-(1/2)^2=0.75\), or 75% of steady state. 3. Five hours after stopping, \(C=3.0e^{-(0.20)(5)}=3.0e^{-1}\approx1.104\) mg/L. 4. \(R_0/CL\) is the equilibrium concentration for continued infusion; the concentration starting from zero at finite time \(t\) is \( (R_0/CL)(1-e^{-kt})\), which is lower than steady state for finite positive \(t\).