Tutorials › Physiologically-Based Pharmacokinetics (PBPK) › When Static and Dynamic Predictions Diverge for AUC and Cmax

Static Versus Dynamic DDI Prediction · Tutorial 484 of 799

When Static and Dynamic Predictions Diverge for AUC and Cmax

Learn why static and dynamic DDI predictions can disagree for AUC and Cmax, and how to trace that disagreement to concentration profiles, enzyme turnover, and victim-drug kinetics.

Advanced 12 min read

What You'll Learn

  • Distinguish AUC ratios from peak-concentration ratios and understand why they need not agree.
  • Derive an AUC result for a periodically varying clearance under specified steady-state conditions.
  • Compare averaging inhibitor concentrations with averaging the resulting nonlinear clearance effects.
  • Quantify how enzyme-turnover delay changes finite-interval AUC and Cmax predictions.
  • Identify how concentration-dependent victim clearance changes the interaction over a concentration profile.
  • Recognize assumptions and diagnostics that matter when interpreting static–dynamic disagreement.

Why the Ratios Can Diverge

A static drug–drug interaction (DDI) prediction represents perpetrator effects using fixed concentrations or fixed enzyme activity. A dynamic prediction allows concentrations, enzyme activity, and victim-drug concentrations to change over time. These approaches need not give the same AUC ratio or Cmax ratio: an AUC summarizes exposure across an interval, whereas Cmax is determined by the highest concentration reached and when it occurs.

The earlier material on static and dynamic DDI models established these as exposure approximations with different time representations. The earlier material on why static models miss DDI onset and offset considered timing. Here the focus is the size and sources of disagreement in the resulting AUC and peak ratios. A difference is not, by itself, evidence that either calculation is incorrect; it may reflect the averaging assumptions, the time interval, or the pharmacokinetic processes represented.

Definition: For a victim drug, the AUC ratio over an interval \([0,T]\) is \(AUCR_{0-T}=AUC_{\mathrm{with\ perpetrator},0-T}/AUC_{\mathrm{control},0-T}\). The Cmax ratio is \(C_{\max,\mathrm{with\ perpetrator}}/C_{\max,\mathrm{control}}\), where each maximum is taken over the same specified interval. A ratio above 1 indicates greater exposure or a higher peak, respectively, with the perpetrator present.

A static prediction is most directly interpretable as an exposure ratio under its assumed fixed conditions, often a steady-state approximation. A dynamic finite-interval ratio instead depends on the simulated concentration-time profiles and the chosen interval. Those quantities coincide only under appropriate assumptions; for example, reaching steady state quickly and maintaining nearly constant interaction effects.

AUC and Cmax Respond Differently

For a linear one-compartment victim model with constant apparent volume \(V\), constant input rate \(R_{\mathrm{in}}\), and time-varying first-order clearance \(CL(t)\), concentration \(C(t)\) satisfies

$$ V\frac{dC(t)}{dt}=R_{\mathrm{in}}-CL(t)C(t). $$

AUC depends on the entire concentration path. Cmax depends on the path's highest point. Consequently, a model can predict a modest AUC change but a more pronounced peak change, or the reverse. For example, rapid inhibition near the time of a victim dose may raise the peak even if the perpetrator effect is brief; a sustained but moderate reduction in clearance may increase AUC without producing an equally large peak change.

One useful result follows by integrating the concentration balance over an interval. If the input is constant and concentration is periodic over a complete cycle, then \(C(T)=C(0)\). This is a specific steady-periodic condition, not a general property of a finite dosing interval.

Derivation: Integrating the balance equation from \(0\) to \(T\) gives \(V[C(T)-C(0)]=R_{\mathrm{in}}T-\int_0^T CL(t)C(t)\,dt\). Under periodic steady state, the left side is zero, so \(\int_0^T CL(t)C(t)\,dt=R_{\mathrm{in}}T\). If clearance is constant, this reduces to \(AUC=R_{\mathrm{in}}T/CL\). If clearance varies, it does not generally reduce to \(R_{\mathrm{in}}T\) divided by an unweighted average clearance: the integral contains the product \(CL(t)C(t)\). In the particular piecewise-constant example below, the clearance-weighted balance can be evaluated from the same mass balance, while Cmax still requires the concentration trajectory.

This distinction cautions against treating “the average clearance” as a universal shortcut. When clearance and concentration vary together, their association over time matters. Moreover, even an exact AUC calculation supplies no direct information about when the maximum occurs.

Worked Example: Fluctuating Inhibitor Concentrations Under Periodic Input

Consider a hypothetical victim drug given by constant infusion. Its baseline first-order clearance is \(0.10\ \mathrm{L/h}\), and the input is \(0.10\ \mathrm{mg/h}\), so its baseline steady concentration is \(1.00\ \mathrm{mg/L}\). Assume 80% of baseline clearance is through a reversibly inhibited pathway and 20% is unaffected. At inhibitor concentration \(I\), let the inhibited pathway's remaining activity be \(1/(1+I/K_i)\). Assume the inhibitor alternates every 6 hours between \(I=0\) and \(I=9K_i\), and the system has reached a repeating 12-hour cycle.

The relative total clearance is \(r(I)=0.20+0.80/(1+I/K_i)\). Thus \(r(0)=1.00\), while \(r(9K_i)=0.20+0.80/10=0.28\). The dynamic cycle's mean relative clearance is \((1.00+0.28)/2=0.640\). Periodic mass balance fixes the integral of clearance times concentration, not AUC divided by the arithmetic mean clearance, so the AUC ratio must be calculated from the concentration trajectory. Let \(p=e^{-0.60}\), \(q=e^{-0.168}\), and \(A=0.10/0.028\); the periodic starting concentration \(x\) and the concentration \(y\) after the uninhibited segment give:

$$ x=\frac{A(1-q)+q(1-p)}{1-pq}\approx1.742,\qquad y=1+(x-1)p\approx1.407,\qquad \frac{AUC_{\mathrm{dynamic}}}{AUC_{\mathrm{baseline}}} =\frac{6+(x-1)(1-p)/0.10+6A+(y-A)(1-q)/0.028}{12} \approx1.5685. $$

A static calculation using the arithmetic mean inhibitor concentration instead uses \(I_{\mathrm{mean}}=4.5K_i\). Its relative clearance is \(0.20+0.80/(1+4.5)\approx0.3455\), giving a steady-state AUC ratio of \(1/(0.20+0.80/5.5)\approx2.895\). The difference arises because the inhibition function is nonlinear: activity at mean concentration is not mean activity.

The peak also differs. Let \(x\) be the concentration at the start of the uninhibited 6-hour segment and \(y\) the concentration at its end. The constant-clearance solution is \(C(t)=C_{\mathrm{eq}}+[C(0)-C_{\mathrm{eq}}]e^{-kt}\). The uninhibited segment has \(k=0.10\ \mathrm{h^{-1}}\) and \(C_{\mathrm{eq}}=1.00\ \mathrm{mg/L}\); the inhibited segment has \(k=0.028\ \mathrm{h^{-1}}\) and \(C_{\mathrm{eq}}=0.10/0.028=3.571\ \mathrm{mg/L}\). Therefore,

$$ y=1+(x-1)e^{-0.60},\qquad x=3.571+(y-3.571)e^{-0.168}. $$

Solving these two equations under the repeating-cycle condition gives \(x\approx1.742\ \mathrm{mg/L}\) and \(y\approx1.407\ \mathrm{mg/L}\). Concentration falls during the uninhibited segment and rises during the inhibited segment, so the cycle maximum is \(x\), just before the return to the uninhibited segment. The dynamic Cmax ratio is therefore about 1.742. Under the static mean-concentration assumption, the constant steady concentration and Cmax ratio are both about 2.895. This comparison depends on the stated infusion, periodic steady state, and model assumptions; it is not a general rule that dynamic Cmax must be lower.

Delayed Enzyme Turnover Changes Finite-Interval Ratios

For time-dependent inhibition or induction, perpetrator concentration does not necessarily translate immediately into enzyme activity. A dynamic model may represent enzyme abundance \(E(t)\) with a turnover equation. One illustrative form for inactivation is \(dE/dt=k_{\mathrm{deg}}(E_{\mathrm{base}}-E)-k_{\mathrm{inact}}(I)E\), where \(E_{\mathrm{base}}\) is baseline enzyme abundance and \(k_{\mathrm{deg}}\) is the degradation rate constant. Enzyme activity, and therefore victim clearance through that pathway, can lag behind a rapid change in inhibitor concentration.

A static prediction based on the eventual enzyme effect can therefore overstate the effect over an early observation interval. Conversely, a static estimate based on a brief measurement of inhibitor concentration may understate a later effect if enzyme activity accumulates or recovers slowly. The finite-interval AUC and peak must be calculated from the concentration profile for that interval; they are not automatically equal to the eventual steady-state ratio.

Worked Example: A Delayed Effect During a 24-Hour Interval

Suppose a hypothetical static steady-state calculation predicts an AUC ratio and Cmax ratio of 2.00 after full enzyme suppression. A dynamic simulation begins at baseline and, because enzyme activity declines gradually, gives victim concentrations relative to the control profile of 1.00, 1.15, 1.35, 1.55, and 1.70 at 0, 6, 12, 18, and 24 hours. Assume the control concentration is constant at \(1.00\ \mathrm{mg/L}\) over this illustrative interval. Use trapezoidal integration for the dynamic AUC.

The control AUC is \(24\ \mathrm{h}\times1.00\ \mathrm{mg/L}=24.0\ \mathrm{mg\cdot h/L}\). The dynamic AUC is the sum of four trapezoids:

$$ AUC_{\mathrm{dynamic}} =6\left(\frac{1.00+1.15}{2} +\frac{1.15+1.35}{2} +\frac{1.35+1.55}{2} +\frac{1.55+1.70}{2}\right) =32.4\ \mathrm{mg\cdot h/L}. $$

Thus \(AUCR_{0-24}=32.4/24.0=1.35\). Since the control is constant at \(1.00\ \mathrm{mg/L}\) and the largest dynamic value in the interval is \(1.70\), the Cmax ratio is \(1.70/1.00=1.70\). Both dynamic ratios are below the static value of 2.00, but they differ from each other and from the static estimate. The listed profile is an illustrative hypothetical simulation output, not a claim that every turnover model produces these values.

A useful diagnostic is to inspect enzyme activity and victim concentration together rather than comparing only final ratios. If enzyme activity changes slowly but a simulated victim concentration responds immediately, check whether the model has inadvertently coupled clearance directly to perpetrator concentration instead of to the turnover state. Also check that the AUC interval is the same in both conditions and that Cmax is obtained from the full profile, including times between reporting points when appropriate.

Changing Victim Concentrations Can Change the Interaction

A static interaction effect may also fail to represent a profile in which victim concentrations vary substantially. With saturable metabolism, the enzyme-mediated rate is not proportional to concentration over the entire range. For example, with a competing linear clearance pathway, the total elimination rate can be written \(v(C)=CL_{\mathrm{other}}C+V_{\max}C/(K_m+C)\). The apparent clearance \(v(C)/C\) then depends on \(C\). If an inhibitor reduces \(V_{\max}\), the size of the clearance change depends on the victim concentration at which it is evaluated.

Worked Example: Concentration-Dependent Local Clearance Ratios

For a hypothetical victim, let \(CL_{\mathrm{other}}=1.00\ \mathrm{L/h}\), \(V_{\max}=9.00\ \mathrm{mg/h}\), and \(K_m=1.00\ \mathrm{mg/L}\). Assume inhibition reduces \(V_{\max}\) by 80%, to \(1.80\ \mathrm{mg/h}\), while leaving the other pathway unchanged. At \(C=1.00\ \mathrm{mg/L}\), baseline apparent clearance is \(1.00+9.00/(1.00+1.00)=5.50\ \mathrm{L/h}\). Under inhibition it is \(1.00+1.80/2.00=1.90\ \mathrm{L/h}\); the local baseline-to-inhibited clearance ratio is \(5.50/1.90=2.895\).

At \(C=9.00\ \mathrm{mg/L}\), baseline apparent clearance is \(1.00+9.00/10.00=1.90\ \mathrm{L/h}\). Inhibited apparent clearance is \(1.00+1.80/10.00=1.18\ \mathrm{L/h}\), so the local ratio is \(1.90/1.18=1.610\). The same 80% reduction in \(V_{\max}\) therefore produces different local clearance ratios at these two concentrations. These are local comparisons, not AUC ratios: predicting AUC and Cmax requires solving the concentration-time model with its input and initial conditions. A static calculation at one concentration cannot automatically represent a profile that moves between such ranges.

Comparing Approaches and Checking for Failure

Static and dynamic calculations answer related but distinct questions. A static calculation is compact and can be useful when concentrations and enzyme effects are adequately represented by fixed values for the intended exposure comparison. A dynamic simulation is needed to represent changing perpetrator concentrations, enzyme turnover, dosing schedules, or concentration-dependent victim disposition explicitly. A dynamic model is not inherently more reliable: an incorrect turnover rate, clearance pathway, or dosing input can produce a detailed but misleading profile.

FeatureStatic approximationDynamic simulation
Time representationFixed concentrations or effectsTime-varying concentrations and states
AUC interpretationOften a steady-state or representative-condition ratioRatio over the specified simulated interval
Cmax interpretationMay represent a steady concentration or simplified peak estimateMaximum of the simulated concentration-time profile
Common sensitivityChoice of representative concentration and nonlinear averagingInput schedules, initial conditions, time steps, and model structure

When results disagree, first align the definitions: use the same victim, control condition, dosing regimen, and AUC interval. Then identify what each calculation averages. In particular, do not substitute an average inhibitor concentration into a nonlinear inhibition function and assume it equals the average enzyme activity. For dynamic outputs, check the concentration-time profile and numerical integration, and verify that the maximum was searched over the whole interval. If changing victim concentrations are important, inspect whether the clearance law is linear over the range encountered.

  • Pitfall: treating a steady-state static ratio as the ratio for every interval. A static prediction based on a fully developed effect need not describe the first hours or days of exposure. Specify the interval and account for turnover or changing concentrations if they matter.
  • Pitfall: assuming the AUC ratio determines the Cmax ratio. AUC is an integrated quantity; the timing and shape of the profile determine Cmax. Calculate and interpret the peak separately.
  • Pitfall: averaging before applying a nonlinear function. In general, \(f(\bar I)\ne\overline{f(I)}\). Calculate the effect over the changing concentrations, or state and justify the approximation.
  • Pitfall: calling a clearance ratio an AUC ratio without checking assumptions. The inverse-clearance relationship requires specified linear, steady-state conditions; it does not generally apply unchanged to time-varying or saturable clearance.

The main limitation of any comparison is that disagreement alone does not identify which representation is appropriate for a particular use. Static predictions can conceal timing and variability, while dynamic predictions depend on the quality of their input profiles and parameterization. The key interpretive discipline is to state the quantity, interval, and assumptions behind each ratio, then trace any difference to the mechanisms actually represented.

Key takeaway: AUC and Cmax are distinct exposure summaries. Differences between static and dynamic ratios can arise from nonlinear averaging, delayed enzyme turnover, and victim concentration-dependent clearance. Compare like intervals and conditions, and use the model's concentration-time profile—not a single ratio—to explain the disagreement.

Check Your Understanding

Answer each question, then compare with the concise solutions below.

  1. Why can applying a nonlinear inhibition function to mean inhibitor concentration differ from averaging the inhibition effect over time?
  2. Under what condition did the periodic mass-balance argument in the first worked example set the net change in victim amount to zero?
  3. In the delayed-turnover example, why is the 24-hour dynamic AUC ratio not automatically 2.00?
  4. Can an AUC ratio alone determine the Cmax ratio? Explain.
  5. In a saturable clearance model, why might the local clearance ratio differ at two victim concentrations even if the fractional reduction in \(V_{\max}\) is unchanged?

Solutions: (1) A nonlinear function does not generally commute with averaging, so \(f(\bar I)\) need not equal the mean of \(f(I)\). (2) The system had reached periodic steady state and the interval covered a complete cycle, so \(C(T)=C(0)\) and the net amount change was zero. (3) The static value represented the eventual full effect, whereas enzyme activity and victim concentration changed gradually during the finite interval; integration of that profile gave 1.35. (4) No. AUC integrates the profile, while Cmax depends on its highest value and timing. (5) Saturable enzyme clearance per unit concentration depends on \(C\), so the affected pathway contributes a different fraction of total apparent clearance at different concentrations.