AUC and Peak-Concentration PTA

Why Calculate PTA for AUC and Peak Concentration?

The same simulated concentration profile can support different exposure summaries. AUC measures accumulated exposure over a specified interval, whereas peak concentration measures the highest concentration reached in that interval. A regimen can therefore attain an AUC target often but a peak target less often, or vice versa. Which index matters depends on the exposure-response evidence and the clinical target being used for the dose-selection estimand.

This tutorial treats PTA as a population probability calculated from simulated profiles for a defined regimen, interval, and target population. The introductory material on probability of target attainment establishes the general idea; here the key task is to define and calculate AUC- and peak-based events consistently. A PTA is not a property of a dose alone: it depends on the selected index, target, population, dosing conditions, and model used to generate the profiles.

Definition: For a simulated individual \(i\) under candidate regimen \(d\), let \(C_i(t;d)\) denote the concentration at time \(t\) in the specified analysis interval \([0,\tau]\). The interval AUC is \(AUC_i(d)=\int_0^\tau C_i(t;d)\,dt\). The peak concentration is \(C_{\max,i}(d)=\max_{0\leq t\leq\tau}C_i(t;d)\). A lower-bound AUC PTA for target \(AUC_{\mathrm{target}}\) is \(P(AUC_i(d)\geq AUC_{\mathrm{target}})\). A peak-target PTA can be defined for a lower bound, \(P(C_{\max,i}(d)\geq C_{\max,\mathrm{target}})\), or for an acceptable range, \(P(C_{\max,\mathrm{low}}\leq C_{\max,i}(d)\leq C_{\max,\mathrm{high}})\). The target event must be specified before calculating PTA.

Define the Index and Target Before Counting Attainment

The interval in the AUC definition must match the intended estimand. For a steady-state regimen dosed every \(\tau\) hours, one common choice is AUC over one dosing interval at steady state. For first-dose exposure, the interval and profile must instead correspond to the first-dose question. These are not interchangeable summaries. Similarly, a peak target must state whether any concentration above a minimum is desirable or whether concentrations above an upper bound are unacceptable.

AUC has concentration multiplied by time as its units, such as mg·h/L; \(C_{\max}\) has concentration units, such as mg/L. Check that target units, profile units, and time units agree. If an interval is expressed in hours, an AUC computed from concentrations in mg/L is in mg·h/L. Rescaling time from hours to days changes the numerical AUC unless its units are rescaled too.

When simulated profiles are available only at discrete times, AUC is commonly approximated by numerical integration. For adjacent times \(t_j\) and \(t_{j+1}\), the linear trapezoidal contribution is the mean of the two endpoint concentrations multiplied by the time increment. The sum approximates the integral. For \(C_{\max}\), taking the largest stored concentration is adequate only if the profile grid captures the peak sufficiently well. A rapidly changing profile may peak between stored times; use a sufficiently fine grid or a numerical maximum of the simulated concentration-time function.

Worked Example: Calculating AUC and Peak From One Profile

Suppose a hypothetical simulated profile over a 12-hour interval is piecewise linear, with concentrations 0, 6, 12, 8, and 3 mg/L at times 0, 2, 4, 8, and 12 hours, respectively. Calculate interval AUC and \(C_{\max}\), then assess an AUC target of at least 80 mg·h/L and a peak target range of 10 to 14 mg/L.

For the linear trapezoidal calculation, add the area from each adjacent pair of observations:

\begin{aligned} AUC &\approx \sum_j \frac{C(t_j)+C(t_{j+1})}{2}(t_{j+1}-t_j)\\ &= \frac{0+6}{2}(2-0) +\frac{6+12}{2}(4-2) +\frac{12+8}{2}(8-4) +\frac{8+3}{2}(12-8)\\ &= 6+18+40+22\\ &=86\ \text{mg·h/L}. \end{aligned}

The largest concentration is 12 mg/L, so \(C_{\max}=12\) mg/L. The profile meets the AUC target because \(86\geq80\), and it meets the peak-range target because \(10\leq12\leq14\). The piecewise-linear assumption makes the trapezoidal area exact for this hypothetical profile. With a model-generated curved profile, the same calculation is an approximation whose accuracy depends on the time grid.

Estimate PTA From Simulated Individuals

Suppose \(N\) simulated individuals represent the population and regimen specified by the estimand. For a target event \(E_i(d)\), such as \(AUC_i(d)\geq AUC_{\mathrm{target}}\), the Monte Carlo estimator is the fraction of simulated individuals meeting that event. The same calculation applies to a peak threshold or range, with the corresponding event substituted.

\widehat{\mathrm{PTA}}(d)=\frac{1}{N}\sum_{i=1}^{N} I\{E_i(d)\}, \qquad SE_{\mathrm{MC}}\{\widehat{\mathrm{PTA}}(d)\} \approx \sqrt{\frac{\widehat{\mathrm{PTA}}(d)[1-\widehat{\mathrm{PTA}}(d)]}{N}}.

Here \(I\{E_i(d)\}\) equals 1 when simulated individual \(i\) meets the target and 0 otherwise. If the simulated individuals are independent draws from the intended predictive population, the estimator is unbiased for the target-event probability, and the displayed standard error estimates its Monte Carlo variability. It does not measure all uncertainty about the model, target, or clinical relevance of the threshold.

Why the Monte Carlo Formula Works

For one simulated individual, define \(Y_i=I\{E_i(d)\}\). Under independent simulation from a fixed target population, \(Y_i\) is Bernoulli with probability \(p=P(E_i(d))\): it is 1 with probability \(p\), and 0 with probability \(1-p\). Therefore \(E(Y_i)=p\) and \(\operatorname{Var}(Y_i)=p(1-p)\). Since \(\widehat{\mathrm{PTA}}=\bar{Y}\), independence gives \(\operatorname{Var}(\bar{Y})=p(1-p)/N\). Replacing the unknown \(p\) with \(\widehat{\mathrm{PTA}}\) yields the estimated Monte Carlo standard error above.

This justification relies on independent draws for the event indicators and a well-defined simulation distribution. If profiles share parameter draws, are generated in correlated batches, or arise from a dependent simulation scheme, treating every profile as an independent draw can understate Monte Carlo error. The error calculation must reflect the actual simulation design. In all cases, a small Monte Carlo standard error indicates numerical precision for the specified simulation setup; it does not establish that the model or target is correct.

Worked Example: Estimating Separate AUC and Peak PTAs

Consider 10 hypothetical simulated individuals under one regimen. Their interval AUCs, in mg·h/L, are 92, 104, 118, 99, 130, 107, 88, 115, 101, and 96. Their peak concentrations, in mg/L, are 11.2, 12.4, 14.1, 10.8, 13.0, 11.9, 12.8, 13.6, 9.7, and 12.1. Use an AUC target of at least 100 mg·h/L and a peak target of at least 12 mg/L.

Six of the 10 AUC values meet the target: 104, 118, 130, 107, 115, and 101 mg·h/L. Thus:

\widehat{\mathrm{PTA}}_{\mathrm{AUC}}=\frac{6}{10}=0.600, \qquad SE_{\mathrm{MC}}\approx\sqrt{\frac{0.600(1-0.600)}{10}}=0.155.

Six peak values meet or exceed 12 mg/L: 12.4, 14.1, 13.0, 12.8, 13.6, and 12.1 mg/L. Therefore:

\widehat{\mathrm{PTA}}_{C_{\max}}=\frac{6}{10}=0.600, \qquad SE_{\mathrm{MC}}\approx\sqrt{\frac{0.600(1-0.600)}{10}}=0.155.

Both estimated PTAs are 60.0%, but that does not mean the same individuals met both targets; the marginal counts alone do not show which events occurred together. With only 10 simulated individuals, each Monte Carlo standard error is about 15.5 percentage points, so the estimates are imprecise as Monte Carlo calculations. A larger simulation may improve numerical precision, while leaving model and target limitations unchanged.

Different Indices Can Rank Regimens Differently

AUC and \(C_{\max}\) summarize different features of a concentration-time profile. AUC reflects total exposure over the selected interval; peak concentration reflects the highest level. A regimen with a relatively concentrated exposure profile can have a higher peak without a proportionate increase in AUC. A different schedule can produce a flatter profile and a lower peak while maintaining or changing total exposure. The simulated model and regimen determine the actual pattern; no general rule says that the regimen preferred by one index must also be preferred by the other.

Consequently, compare each candidate against the target event that is relevant to the dose-selection question. Do not substitute mean AUC or mean \(C_{\max}\) for the probability of individual attainment: a population mean can meet a threshold even when many individuals do not, or miss it when many individuals do meet it. Report the target, interval, simulated population, PTA estimate, and Monte Carlo precision together.

Worked Example: AUC and Peak Targets Rank Candidates Differently

Suppose two hypothetical candidate regimens are evaluated using 100 simulated individuals per regimen from the same target population and under the same analysis conditions. Regimen A has 88 individuals with AUC at least 100 mg·h/L; regimen B has 72. For a peak target range of 10 to 16 mg/L, regimen A has 58 individuals in range and regimen B has 81.

For the AUC target, the estimated PTAs are:

\widehat{\mathrm{PTA}}_{\mathrm{AUC}}(A)=\frac{88}{100}=0.880, \qquad \widehat{\mathrm{PTA}}_{\mathrm{AUC}}(B)=\frac{72}{100}=0.720.

Regimen A has the higher AUC PTA, by 0.160 or 16.0 percentage points. For the peak range, the estimates are:

\widehat{\mathrm{PTA}}_{C_{\max}\ \mathrm{range}}(A)=\frac{58}{100}=0.580, \qquad \widehat{\mathrm{PTA}}_{C_{\max}\ \mathrm{range}}(B)=\frac{81}{100}=0.810.

Regimen B has the higher peak-range PTA, by 0.230 or 23.0 percentage points. These hypothetical results show a ranking reversal: the AUC criterion favors A, while the specified peak range favors B. This is not a contradiction. The targets describe different exposure properties, and the simulated profiles can differ in their total exposure and their peak distribution. Which candidate is preferred requires a decision rule grounded in the relevant clinical target; these marginal calculations alone do not combine the criteria into a simultaneous requirement.

Diagnostics, Alternatives, and Limitations

Before interpreting an AUC- or peak-based PTA, verify that the exposure summaries actually represent the intended simulation. The following checks are especially useful:

Confirm the interval and dosing condition. Ensure that every profile uses the specified first-dose or steady-state condition and the same intended interval. AUC values from different intervals are not directly comparable without an explicit reason.
Check numerical integration and peak capture. Inspect representative profiles, time grids, and calculated summaries. Sparse time points can miss a narrow peak or distort the AUC, particularly where concentration changes rapidly.
Check units and target boundaries. Confirm time and concentration units, and check how values exactly equal to a threshold are classified. For a lower-bound event, equality is included when the definition uses \(\geq\); for a range, specify whether each boundary is included.
Assess simulation precision and integrity. Report the number of valid simulations and Monte Carlo error. Investigate failures, implausible profiles, or patterns in failed draws rather than silently excluding them in a way that changes the target population.
Look beyond a single threshold. A small change in a target can change the attainment fraction, especially when many simulated values lie near the boundary. Examine the exposure distribution or calculate PTA at clinically justified nearby thresholds as a sensitivity check.

AUC PTA and peak PTA are not the only useful summaries. Mean exposure describes an average, but not the fraction of individuals meeting a target. Time-above-target PTA uses the duration or proportion of an interval above a concentration threshold and answers a different exposure question. Choosing among these indices requires exposure-response justification rather than selecting whichever produces the most favorable dose ranking.

PTA remains conditional on the simulated model, population, dosing assumptions, and target definition. A large number of simulations can make the Monte Carlo error small without resolving model misspecification, weak exposure-response evidence, or an unsuitable target. Also, monotonicity across doses should not be assumed automatically: nonlinear pharmacokinetics, changing regimens, or a bounded peak target can produce non-monotone attainment patterns. Evaluate each candidate with the same definitions and simulation procedures.

Common pitfall: “Regimen A has the higher mean AUC, so it must have the higher AUC PTA and the better peak attainment.” This reasoning is not valid. A mean does not determine the fraction above an AUC threshold, and an AUC summary does not determine the distribution of peak concentrations. Calculate the target-event indicator for each simulated profile, estimate each PTA separately, and interpret each result against its own justified target.

Check Your Understanding

Answer each question using the definitions and calculations above.

A simulated profile has AUC 95 mg·h/L and \(C_{\max}=15\) mg/L. Does it meet an AUC target of at least 90 mg·h/L and a peak range of 10 to 14 mg/L?
In 200 independent simulated individuals, 150 meet a lower-bound AUC target. Calculate the estimated PTA and its Monte Carlo standard error.
Why can two regimens have different rankings for AUC PTA and peak-range PTA?
What is one reason that the maximum concentration on a sparse time grid might not equal the profile’s true \(C_{\max}\)?

Solutions. (1) It meets the AUC target because \(95\geq90\), but it does not meet the peak range because 15 mg/L exceeds the upper boundary of 14 mg/L. (2) The estimate is \(150/200=0.750\). Its Monte Carlo standard error is \(\sqrt{0.750(1-0.750)/200}=\sqrt{0.0009375}=0.0306\), or about 3.06 percentage points. (3) AUC describes accumulated exposure, whereas a peak-range event describes the maximum concentration; candidates can have different exposure profiles and therefore different attainment fractions for the two events. (4) The peak may occur between stored time points, especially if concentration changes rapidly, so a sparse grid can miss it.

What You'll Learn