Extrapolating AUC to Infinity

Why Extend AUC Beyond the Last Sample?

A concentration–time profile is observed only at finite sampling times, but exposure is often summarized as the area from time zero to infinity, \(AUC_{0-\infty}\). The unobserved part after the last quantifiable concentration is not measured directly. It is estimated by extending a model of the terminal decline beyond the final observation. The resulting total therefore combines an area calculated from observed data with an area inferred from a terminal-phase model.

The extrapolation uses two quantities: the concentration at the last quantifiable time and the positive terminal rate constant, \(\lambda_z\). The earlier material on selecting the terminal-phase slope addresses how \(\lambda_z\) is chosen. Here, the focus is what that rate implies for the remaining area and how to combine that area with the observed AUC.

Deriving the Residual-Area Formula

Assume that after \(t_{\text{last}}\), concentration follows a monoexponential terminal decline with the same positive rate constant \(\lambda_z\). This means concentration at a time \(t\) after the last sample is modeled as \(C_{\text{last}}\) multiplied by an exponential decay from that last concentration. The residual area is the integral of this modeled concentration from the last sampling time to infinity.

To integrate, set \(u=t-t_{\text{last}}\). When \(t=t_{\text{last}}\), \(u=0\); as \(t\) tends to infinity, so does \(u\). The change of variable gives the following integral and evaluation:

The final equality requires \(\lambda_z>0\), so that \(e^{-\lambda_z u}\) tends to zero as \(u\) tends to infinity. The total AUC estimate is then the observed area plus this residual area:

The units provide a useful check. If concentration is measured in mg/L and time in hours, then \(\lambda_z\) has units of h\(^{-1}\). Thus \(C_{\text{last}}/\lambda_z\) has units of mg·h/L, the same as AUC. The formula also shows the roles of its inputs: a larger last concentration increases the estimated tail area, while a larger terminal rate constant decreases it.

Calculating the Tail and Total AUC

The observed area, \(AUC_{0-t_{\text{last}}}\), is calculated from the concentration data through the last quantifiable time using the appropriate numerical integration method. The extrapolation does not replace that observed-data calculation: it adds an estimated terminal tail to it. Before applying the formula, confirm that \(C_{\text{last}}\) is the concentration at the stated final quantifiable time and that \(\lambda_z\) is expressed in the reciprocal of the time unit used for AUC.

A useful descriptive quantity is the percentage of the total AUC that is extrapolated. It is calculated using the residual area in the numerator and the estimated total AUC in the denominator:

For the preceding example, this is \(100\times 6.67/486.67=1.37\%\). It describes the share of the estimated total contributed by the extrapolation; it does not by itself prove that the terminal model is appropriate.

Combining Observed Segments With the Extrapolation

The next example shows the components together. The observed AUC is calculated from the stated concentrations; the tail is calculated from the final concentration and terminal rate. For a declining interval with positive concentrations, the log-trapezoidal area is \((C_1-C_2)/\ln(C_1/C_2)\) multiplied by the time interval. The first interval below rises from zero, so its area is calculated with a linear trapezoid.

Assumptions, Interpretation, and Alternatives

The residual-area formula is exact for the assumed exponential tail, but the tail itself is a model-based continuation. Its use requires a positive terminal rate and a terminal decline that is reasonably represented by that rate beyond the last quantifiable sample. The last concentration must be quantifiable and correctly paired with its sampling time. The formula does not require the entire concentration–time profile to be monoexponential; it applies to the modeled terminal portion after \(t_{\text{last}}\).

A terminal rate used for this extrapolation is not automatically the drug’s elimination rate in every setting. For example, after extravascular dosing, the terminal decline can reflect absorption as well as elimination. The formula still computes the area under the assumed terminal concentration curve, but interpretation of \(\lambda_z\) as a specific physiological elimination constant requires additional justification.

A basic diagnostic check is whether the estimated residual area is positive and has the same units as the observed AUC. A zero or negative \(\lambda_z\) cannot yield a finite, decaying exponential tail under this model. Also check that the chosen terminal rate corresponds to the terminal phase and that the last concentration is not a substituted value such as zero. AUC extrapolation can be sensitive to the rate: a small positive \(\lambda_z\) produces a large residual area. More detailed evaluation of whether an extrapolation is reliable is a separate question.

There are two common alternatives in reporting or analysis. One may report \(AUC_{0-t_{\text{last}}}\) alone, which is based on observed data but does not represent the full area to infinity. Alternatively, a compartmental or other model-based analysis can integrate a fitted concentration curve, including its tail; that approach depends on the chosen model and its fit. The terminal-rate formula is a direct noncompartmental extrapolation and should not be mistaken for a separately observed portion of exposure.