Tutorials › Nonlinear Mixed-Effects Models (NLME) › Repeated-Dose Accumulation With Delayed and Extended Absorption

Absorption Models · Tutorial 84 of 797

Repeated-Dose Accumulation With Delayed and Extended Absorption

Calculate repeated-dose concentrations by adding each dose’s contribution when absorption is lagged, extended, or distributed over time.

Intermediate 10 min read

What You'll Learn

  • State when repeated-dose contributions can be added under linear pharmacokinetics.
  • Calculate concentrations from lagged first-order absorption after multiple doses.
  • Calculate concentrations during overlapping zero-order input.
  • Use a transit-time input distribution to calculate a repeated-dose concentration.
  • Distinguish ongoing absorption from accumulation due to elimination.
  • Recognize when simple accumulation shortcuts do not apply.

Why Repeated Doses Can Overlap in More Than One Way

A new dose does not necessarily arrive only after the previous dose has finished entering the body. A lag time can postpone absorption, a zero-order input can continue at a fixed rate for hours, and transit compartments can spread absorption over time. With repeated dosing, these inputs may overlap. The concentration at a sampling time must then reflect both drug already in the body and drug still entering from one or more earlier doses.

The calculations below use linear one-compartment disposition: drug entering the systemic circulation is distributed into a fixed apparent volume \(V\), and is eliminated with first-order rate constant \(k_e\). The central practical result is to calculate the contribution from each dose at the time of interest, then add those contributions. This addition is valid only under the stated linearity assumptions.

Definition: A dose’s concentration contribution at time \(t\) is the concentration that dose would produce by itself under the specified absorption and disposition model. Under linear disposition, the total concentration is the sum of contributions from all doses administered by time \(t\).

Conditions for Adding Dose Contributions

Assume each dose has the same bioavailability \(F\), apparent volume \(V\), and absorption and elimination parameters, unless a worked example says otherwise. For a dose \(D\), let \(R_{\mathrm{in}}(s)\) denote its systemic input rate at time \(s\), measured in amount per time. The concentration contribution at time \(t\) is

$$ C(t)=\frac{1}{V}\int_0^t R_{\mathrm{in}}(s)e^{-k_e(t-s)}\,ds. $$

The exponential factor is the fraction of drug entering at time \(s\) that remains in the body at time \(t\), under first-order elimination. If the total input is the sum of inputs from doses at times \(t_j\), then the integral of that sum equals the sum of the integrals. Therefore,

$$ C_{\mathrm{total}}(t)=\sum_{j:t_j\leq t} C_j(t). $$

This superposition requires linear disposition: \(V\) and \(k_e\) do not change with concentration, and one dose does not alter another dose’s absorption or disposition. It also assumes that the specified input profiles correctly represent each dose. Saturable elimination, concentration-dependent bioavailability, or dose-dependent absorption can invalidate simple addition.

Key result: Under linear one-compartment disposition, calculate each dose’s input-driven concentration contribution at the observation time and add the contributions. Ongoing absorption from an earlier dose is included through that dose’s input profile; it is not a separate dose.

Lagged First-Order Absorption

For a dose administered at time \(t_j\), suppose absorption begins after lag time \(T_{\mathrm{lag}}\) and is then first order with rate constant \(k_a\). Its absorption start time is \(\tau_j=t_j+T_{\mathrm{lag}}\). For \(t\geq\tau_j\), the contribution is

$$ C_j(t)=\frac{F D k_a}{V(k_a-k_e)} \left[e^{-k_e(t-\tau_j)}-e^{-k_a(t-\tau_j)}\right], \qquad k_a\ne k_e. $$

Before \(\tau_j\), this dose contributes zero concentration in this model. After absorption begins, the expression accounts for both continuing absorption and elimination. It is not appropriate to treat a lagged dose as if the full dose entered the body instantaneously at \(\tau_j\).

Worked Example: Two Lagged Doses at a Sampling Time

Consider hypothetical doses of \(100\) mg every \(8\) hours, with \(F=0.80\), \(V=20\) L, \(k_a=1.2\ \mathrm{h}^{-1}\), \(k_e=0.20\ \mathrm{h}^{-1}\), and \(T_{\mathrm{lag}}=2\) h. Find the concentration contribution at \(t=18\) h from doses given at \(t=0\), \(8\), and \(16\) h.

The dose at \(16\) h has not begun absorption at \(18\) h: its absorption start time is \(16+2=18\) h, so its contribution at that instant is zero. For the doses at \(0\) and \(8\) h, the elapsed times since absorption began are \(a=18-(0+2)=16\) h and \(a=18-(8+2)=8\) h. The shared coefficient is

$$ \frac{F D k_a}{V(k_a-k_e)} =\frac{(0.80)(100)(1.2)}{20(1.2-0.20)} =4.8\ \mathrm{mg/L}. $$

Thus the first two contributions are

$$ C_0(18)=4.8(e^{-0.20(16)}-e^{-1.2(16)}) \approx 0.1957\ \mathrm{mg/L}, $$
$$ C_8(18)=4.8(e^{-0.20(8)}-e^{-1.2(8)}) =4.8(0.201897-0.0000677) \approx 0.9688\ \mathrm{mg/L}. $$

Using the unrounded contributions, \(C_{\mathrm{total}}(18)\approx1.1644\ \mathrm{mg/L}\). The dose at \(16\) h is at its absorption start time, not an instantaneous systemic input, so it adds no concentration at that exact time under the lagged first-order model.

Extended Zero-Order Absorption

Suppose each dose’s bioavailable amount \(FD\) enters at a constant rate over duration \(T\). Its input rate is \(FD/T\) from dose time \(t_j\) through \(t_j+T\), and zero outside that interval. Let \(a=t-t_j\) be elapsed time since that dose. While input is ongoing, \(0\leq a\leq T\), integration gives

$$ C_j(t)=\frac{F D}{V T k_e}\left(1-e^{-k_e a}\right). $$

If input has ended, \(a>T\), the concentration at the end of input is reduced by elimination over the remaining time:

$$ C_j(t)=\frac{F D}{V T k_e} \left(1-e^{-k_e T}\right)e^{-k_e(a-T)}. $$

These expressions assume the full bioavailable amount is delivered uniformly across \(T\). If the next dose is given before \(T\) has elapsed, the input rates from the two doses overlap, even though their concentration contributions are still computed separately and added.

Worked Example: A New Dose Arrives Before Zero-Order Input Ends

Suppose hypothetical \(120\)-mg doses are given every \(8\) h. Let \(F=0.75\), \(V=30\) L, \(k_e=0.10\ \mathrm{h}^{-1}\), and let each dose enter over \(T=10\) h. Find the concentration at \(t=9\) h.

Each dose supplies \(FD=(0.75)(120)=90\) mg systemically over \(10\) h. Since the next dose was administered at \(8\) h, both the dose at \(0\) h and the dose at \(8\) h are still being absorbed at \(t=9\) h. For an elapsed input time \(a\), the coefficient is

$$ \frac{F D}{V T k_e}=\frac{90}{(30)(10)(0.10)}=3.00\ \mathrm{mg/L}. $$

For the first dose, \(a=9\) h; for the second, \(a=1\) h. Both are less than \(T\), so

$$ C_0(9)=3.00(1-e^{-0.10(9)}) =3.00(1-0.406570) \approx1.780\ \mathrm{mg/L}, $$
$$ C_8(9)=3.00(1-e^{-0.10(1)}) =3.00(1-0.904837) \approx0.2855\ \mathrm{mg/L}. $$

The total is approximately \(1.780+0.2855=2.066\ \mathrm{mg/L}\). During the interval from \(8\) to \(10\) h, the two input rates are each \(90/10=9.00\) mg/h, for a combined input rate of \(18.0\) mg/h. The concentration calculation still uses the two dose-specific contributions because each has a different history of elimination.

Transit Absorption and Repeated Doses

Transit absorption distributes the time to systemic input rather than assigning one fixed delay. For an Erlang transit-time distribution with \(N\) transit compartments and transit rate \(k_{\mathrm{tr}}\), the absorption-time density for \(u\geq0\) is

$$ f(u)=\frac{k_{\mathrm{tr}}^N u^{N-1}e^{-k_{\mathrm{tr}}u}}{(N-1)!}. $$

This density integrates to one. A dose’s input rate is \(FD f(u)\) at elapsed time \(u\), and its concentration contribution follows by integrating that input against the elimination kernel. The earlier material on transit compartments describes this distributed delay; here the important repeated-dose step is to evaluate each dose’s contribution at its own elapsed time and sum.

For \(N=2\), \(f(u)=k_{\mathrm{tr}}^2u e^{-k_{\mathrm{tr}}u}\). With \(a=t-t_j\geq0\), the contribution is

$$ C_j(t)=\frac{F D}{V} k_{\mathrm{tr}}^2 e^{-k_e a} \int_0^a u e^{-(k_{\mathrm{tr}}-k_e)u}\,du. $$

When \(k_{\mathrm{tr}}\ne k_e\), integrating \(u e^{-xu}\) gives \(\{1-e^{-xa}(1+xa)\}/x^2\), with \(x=k_{\mathrm{tr}}-k_e\). This yields a directly evaluable expression for this two-compartment transit input.

Worked Example: Summing Two-Dose Transit Contributions

For hypothetical \(100\)-mg doses at \(0\) and \(4\) h, take \(F=1\), \(V=20\) L, \(N=2\), \(k_{\mathrm{tr}}=0.50\ \mathrm{h}^{-1}\), and \(k_e=0.20\ \mathrm{h}^{-1}\). Calculate concentration at \(t=5\) h. The dose-specific elapsed times are \(a=5\) h and \(a=1\) h, and \(FD/V=5.00\ \mathrm{mg/L}\). Here \(x=0.50-0.20=0.30\ \mathrm{h}^{-1}\), so

$$ C_j(t)=\frac{(5.00)(0.50)^2 e^{-0.20a}}{(0.30)^2} \left[1-e^{-0.30a}(1+0.30a)\right]. $$

For the dose at \(0\) h,

$$ C_0(5)=13.8889e^{-1} \left[1-e^{-1.5}(2.5)\right] \approx 2.259\ \mathrm{mg/L}. $$

For the dose at \(4\) h,

$$ C_4(5)=13.8889e^{-0.2} \left[1-e^{-0.3}(1.3)\right] \approx0.4200\ \mathrm{mg/L}. $$

Therefore \(C_{\mathrm{total}}(5)\approx2.259+0.4200=2.679\ \mathrm{mg/L}\). The second dose is still being absorbed at the sampling time; its concentration contribution is not equal to the fraction of its dose that has entered, because elimination acts on drug as it enters.

Derivation, Diagnostics, and Common Pitfalls

The justification for superposition follows directly from the concentration integral. Write the total input as \(R_{\mathrm{in}}(s)=\sum_j R_j(s)\). Then

$$ C_{\mathrm{total}}(t) =\frac{1}{V}\int_0^t \left(\sum_j R_j(s)\right)e^{-k_e(t-s)}\,ds =\sum_j \frac{1}{V}\int_0^t R_j(s)e^{-k_e(t-s)}\,ds =\sum_j C_j(t). $$

The equality is valid because integration is linear and the same elimination kernel applies to every dose. It does not imply that absorption profiles are identical at a given clock time: each is shifted by its dose time and has its own elapsed absorption time.

When checking a calculation, verify that each dose has been included only after its administration, and use elapsed time relative to that dose’s lag or input start. For finite-duration zero-order input, use the ongoing-input expression only when \(0\leq a\leq T\); after \(T\), use the post-input expression. For transit input, the contribution must be zero before the dose and nonnegative afterward. If the model parameters are positive, a computed negative concentration signals an algebraic, timing, or numerical error.

A frequent mistake is to assume a new dose resets or replaces absorption from the earlier dose. It does not under the models above: each input continues according to its own profile. Another mistake is to multiply a single-dose concentration by a dose count as a shortcut. That is justified only in special circumstances, such as identical contributions at the time of interest; generally, doses have different elapsed times and different elimination histories. At periodic steady state, repeated profiles may be represented by a sum over prior doses, but the same dose-specific timing logic still applies.

The alternatives differ in their representation of input, not in the superposition principle. Lagged first-order absorption uses a waiting period followed by an exponentially declining input rate. Zero-order absorption uses a constant rate over a defined duration. Transit absorption spreads input according to a distribution of transit times. The choice should follow the intended absorption mechanism and data support; these formulas do not determine which mechanism is appropriate. They also do not cover nonlinear disposition or changing bioavailability.

Common pitfall: “The next dose was given, so the previous dose’s absorption is over.” This is wrong whenever the previous input duration extends past the next dose time, or when a distributed absorption process still has input. Check each dose’s absorption window separately, calculate its contribution at the observation time, and then add the contributions.

Check Your Understanding

Answer each question, then compare your reasoning with the brief solutions.

  1. Under what key disposition condition can concentration contributions from different doses be added?
  2. A zero-order input lasts \(12\) h, and doses are given every \(8\) h. For how long do two consecutive doses’ inputs overlap?
  3. For lagged first-order absorption, a dose is given at \(t=5\) h and has a lag of \(3\) h. At what time does its modeled absorption begin?
  4. Why is multiplying one dose’s concentration contribution by the number of doses generally incorrect?

Solutions. (1) Addition requires linear disposition with a common concentration-independent elimination process and volume, so that superposition holds. (2) The overlap lasts \(12-8=4\) h, from the second dose at \(8\) h until the first dose’s input ends at \(12\) h. (3) Absorption begins at \(5+3=8\) h. (4) Contributions generally differ because each dose has a different elapsed absorption time and has undergone elimination for a different duration; calculate them separately and sum them.