Lag-Time Models for Delayed Input

Why Model an Absorption Lag?

An oral dose need not begin entering the systemic circulation at the recorded dosing time. A delay can arise before absorption starts, so a model that begins first-order input immediately may predict concentrations too early. A fixed absorption lag represents this delay by shifting the start of input while leaving the subsequent absorption process first order.

The central modeling idea is simple: before the lag has elapsed, the dose contributes no systemic input; after it has elapsed, the usual first-order absorption process begins. This distinction matters because a lag shifts the profile in time, whereas the absorption and disposition rate constants determine its shape once input begins.

Model and Assumptions

Consider one oral dose \(D\) given at time zero. Let \(F\) be the fraction of the dose reaching the systemic circulation, \(V\) the apparent volume of distribution, \(k_a\) the first-order absorption rate constant, and \(k_e\) the first-order elimination rate constant. Both rate constants have units of inverse time; \(D\) is an amount and \(V\) is a volume. Assume a one-compartment disposition model with concentration \(C(t)\) equal to the amount in the central compartment divided by \(V\).

The profile below assumes one dose, no drug in the body at time zero, no input before the lag ends, and first-order absorption after that point. It also assumes linear disposition and constant \(F\), \(V\), \(k_a\), and \(k_e\). The model is not intended to represent a distributed delay or input at a constant rate.

The depot amount remaining after absorption begins is \(D e^{-k_a\tau}\), so the systemic input rate is \(F k_a D e^{-k_a\tau}\). If \(A(t)\) denotes the central-compartment amount, its rate of change after the lag is the systemic input rate minus elimination:

This initial condition is important: at the moment absorption begins, the model has not yet delivered any of the dose to the central compartment. To solve the equation, multiply both sides by the integrating factor \(e^{k_e\tau}\), then integrate from zero to \(\tau\):

Multiplying by \(e^{-k_e\tau}\) gives \(A=F D k_a(e^{-k_e\tau}-e^{-k_a\tau})/(k_a-k_e)\). Dividing by \(V\) gives the concentration formula. If \(k_a=k_e=k\), the integral instead becomes \(\int_0^\tau 1\,ds=\tau\), so \(C(t)=F D k\tau e^{-k\tau}/V\) after the lag.

Calculating a Shifted Profile

For a single dose with the same parameters and disposition, the lagged profile is the zero-lag profile evaluated at \(t-t_{\mathrm{lag}}\), once \(t>t_{\mathrm{lag}}\). Before then, the lagged profile is zero. Consequently, changing only \(t_{\mathrm{lag}}\) translates the post-lag curve horizontally; it does not alter its peak height or its shape as a function of elapsed absorption time.

How the Lag Changes the Peak Time

For \(k_a\ne k_e\), differentiating the post-lag concentration with respect to \(\tau\) gives a factor proportional to \(-k_e e^{-k_e\tau}+k_a e^{-k_a\tau}\). Setting this derivative to zero yields \(k_a e^{-k_a\tau}=k_e e^{-k_e\tau}\). Taking logarithms and solving gives the time from absorption onset to the peak:

For positive, unequal rates, \(\tau_{\max}\) is positive. The lag adds to the time of the peak but does not otherwise enter the expression for \(\tau_{\max}\). This result is a useful check on the model's time shift.

Checks, Alternatives, and Limitations

A useful first diagnostic is whether the model's implied timing is compatible with the observations and the dosing record. If a substantial measured concentration is present before the proposed lag ends, a single-dose model with no baseline concentration and no other input cannot explain it. If early samples are sparse, however, a lag may be weakly informed by the data; absence of an observed early concentration is not by itself proof that a true delay occurred.

The main simple alternative is the same first-order absorption model with \(t_{\mathrm{lag}}=0\). It predicts input immediately after dosing. A fixed lag instead imposes no input during a specified initial interval, then uses the same first-order input. When comparing these formulations, inspect early-time predictions and residual patterns as well as the overall profile. Systematic early prediction errors can indicate that neither assumed onset time is adequate, or that some other model assumption fails.

A fixed lag is a simplified representation. It assumes a sharply defined onset, rather than a gradual or distributed delay, and it does not by itself explain the biological mechanism of delayed input. Further, with few observations near the onset, different combinations of timing and absorption behavior may produce similar observed profiles. The concentration equation describes the consequences of a specified lag; it does not guarantee that the data identify that lag precisely.

Before applying the formula, check that the dose is a single dose at time zero, that the sample is after the lag, that time and rate units agree, and that the unequal-rate formula is being used only when \(k_a\ne k_e\). The special equal-rate limit avoids division by zero when the two rates coincide. Finally, confirm that the calculated coefficient has concentration units: \(F\) is dimensionless, and the dose divided by volume is an amount per volume.