Static and Dynamic DDI Models as Exposure Approximations

Why the Exposure Approximation Matters

A drug-drug interaction (DDI) can change a victim drug's exposure because the perpetrator changes a process affecting the victim, such as metabolic activity. A prediction must therefore represent both the perpetrator's exposure and the biological response to it. One approach uses a small number of representative values in an algebraic calculation. Another follows concentrations and biological activity through time. These approaches can give similar summary estimates under suitable conditions, but they do not make the same assumptions.

The distinction is not simply “simple model” versus “complex model.” A static calculation compresses time-varying inputs and responses into selected values; a dynamic simulation retains their trajectories and the order in which changes occur. The choice matters when activity changes appreciably during the period over which victim exposure is measured, or when the response to concentration is nonlinear or delayed.

Definition: A static DDI prediction is an algebraic approximation that uses one or more representative exposure or biological-activity values, rather than explicitly tracking their full time courses. A dynamic DDI simulation represents concentrations and relevant biological states as functions of time and uses them to calculate the victim drug's changing disposition and exposure.

A Minimal Model for Comparing Approaches

To make the role of time explicit, consider a deliberately simplified victim drug model. Suppose the victim is given as an intravenous bolus, its initial concentration is \(C_0\), and its only changing disposition process is elimination. Let \(k_0\) be the baseline first-order elimination rate constant, in inverse hours. Let \(q(t)\) be the relative activity of the process affecting elimination: \(q(t)=1\) means baseline activity, and lower values mean reduced activity. In this illustrative model,

\frac{dC(t)}{dt}=-k_0q(t)C(t), \qquad C(0)=C_0.

This model assumes linear elimination at each instant, a fixed volume of distribution, no change in the victim's input after the bolus, and that the interaction changes only the activity multiplier \(q(t)\). It is not a general PBPK model: it omits absorption, distribution among tissues, competing pathways, and changes in other physiological processes. It is useful here because it isolates the difference between holding activity constant and allowing it to vary.

If activity is constant at \(q\) over an interval of length \(T\), then concentration declines exponentially and the area under the concentration-time curve (AUC) over that interval follows directly. Separating variables and integrating from 0 to \(T\) gives:

\frac{dC}{C}=-k_0q\,dt \quad\Longrightarrow\quad C(t)=C_0e^{-k_0qt}, \qquad \mathrm{AUC}_{0-T}=\int_0^T C(t)\,dt =\frac{C_0\left(1-e^{-k_0qT}\right)}{k_0q}.

When \(q=0\), the limiting result is \(\mathrm{AUC}_{0-T}=C_0T\), because concentration does not decline in this idealized model. A static prediction can use the constant-activity expression with a chosen representative \(q\). A dynamic calculation instead applies the differential equation using the activity at each time, then integrates the resulting concentration profile. The static value must be chosen to represent the question of interest; it is not automatically equal to a time average or to the activity evaluated at an average concentration.

Worked Example: A Constant Representative Activity

Suppose, hypothetically, a victim drug has \(C_0=10.0\) mg/L and \(k_0=0.200\) h\(^{-1}\). For an 8.00-hour interval, a static approximation represents the interaction by constant activity \(q=0.750\). Calculate the predicted AUC and compare it with baseline activity.

The constant-activity assumptions allow use of the finite-interval expression. The effective elimination rate is \(k_0q=(0.200)(0.750)=0.150\) h\(^{-1}\), so:

\mathrm{AUC}_{0-8} =\frac{(10.0\ \mathrm{mg/L})\left(1-e^{-(0.150\ \mathrm{h}^{-1})(8.00\ \mathrm{h})}\right)} {0.150\ \mathrm{h}^{-1}} =\frac{10.0(1-e^{-1.20})}{0.150} =46.59\ \mathrm{mg\cdot h/L}.

At baseline, \(q=1\), and the same calculation gives:

\mathrm{AUC}_{0-8,\mathrm{baseline}} =\frac{10.0(1-e^{-(0.200)(8.00)})}{0.200} =\frac{10.0(1-e^{-1.60})}{0.200} =39.91\ \mathrm{mg\cdot h/L}.

The static model therefore predicts an AUC increase of \(46.59-39.91=6.68\) mg·h/L, or about \(16.7\%\) relative to baseline. This is a consequence of the stated constant-\(q\), one-compartment assumptions; it is not a universal prediction for a victim drug exposed to a real perpetrator.

Why the Order of Activity Changes Matters

A dynamic simulation can distinguish activity profiles that have the same average. In the simple model, the solution for changing activity is:

C(t)=C_0\exp\left(-k_0\int_0^t q(s)\,ds\right).

This shows that the concentration at a given time depends on cumulative activity up to that time. But AUC depends on the concentrations throughout the interval, not just on the final concentration. Consequently, two profiles can have the same cumulative activity by the end of an interval and still yield different AUCs if one has lower activity earlier, while concentrations are higher for longer.

Worked Example: Equal Average Activity, Different AUC

Use the hypothetical victim parameters \(C_0=10.0\) mg/L and \(k_0=0.200\) h\(^{-1}\) over 8.00 hours. Compare a dynamic profile with \(q=0.500\) from 0 to 4.00 hours and \(q=1.00\) from 4.00 to 8.00 hours with a static profile having \(q=0.750\) throughout. Both have a time-average activity of 0.750. Assume the stated piecewise-constant activities apply exactly over their intervals.

During the first four hours of the dynamic profile, the elimination rate is \((0.200)(0.500)=0.100\) h\(^{-1}\). The first-interval AUC and the concentration at hour 4 are:

\mathrm{AUC}_{0-4} =\frac{10.0(1-e^{-(0.100)(4.00)})}{0.100} =32.97\ \mathrm{mg\cdot h/L}, \qquad C(4)=10.0e^{-0.400}=6.703\ \mathrm{mg/L}.

During the next four hours, activity returns to baseline, so the rate constant is \(0.200\) h\(^{-1}\). Starting from \(C(4)\), the second-interval AUC is:

\mathrm{AUC}_{4-8} =\frac{6.703(1-e^{-(0.200)(4.00)})}{0.200} =18.46\ \mathrm{mg\cdot h/L}.

Thus the dynamic AUC, calculated from the unrounded component values, is approximately 51.42 mg·h/L. The static prediction for \(q=0.750\), calculated in the preceding example, is 46.59 mg·h/L. The dynamic result is about \(10.4\%\) higher than this static result, even though both profiles have the same average activity. Both profiles produce the same concentration at hour 8, \(10.0e^{-1.20}=3.012\) mg/L, because their cumulative activity is equal. Their AUCs differ because the dynamic profile suppresses elimination earlier, leaving more concentration during much of the interval.

Concentration Averages and Biological Response

A static approximation may choose a representative perpetrator concentration and map it to a representative biological activity. This shortcut needs care. If activity is a nonlinear function of concentration, the activity at the mean concentration generally differs from the mean activity over time. A dynamic model can evaluate the response at each time; a static model must decide which concentration or response summary is appropriate for its intended prediction.

Worked Example: A Nonlinear Response to Concentration

Consider a hypothetical response relationship, already specified for this example, with relative activity \(q(0)=1.00\), \(q(4)=0.400\), and \(q(8)=0.200\), where concentration is in arbitrary units. Suppose a perpetrator concentration spends equal durations at 0 and 8 units. Assume these values are the relevant time-weighted states and use the listed response values without interpolation.

The mean concentration is \((0+8)/2=4\) units. Evaluating the response at that representative concentration gives \(q(4)=0.400\). But the mean of the two time-specific activity values is:

\bar{q}_{\mathrm{time}}=\frac{q(0)+q(8)}{2} =\frac{1.00+0.200}{2}=0.600.

The response at the mean concentration is therefore 0.400, whereas the mean response is 0.600. Neither number is automatically the correct input for every static prediction: the appropriate summary depends on the model and the exposure metric being approximated. The calculation demonstrates why substituting an average concentration into a nonlinear response relationship cannot be assumed to reproduce the average biological activity.

Assumptions, Diagnostics, and Choice of Approach

Static predictions are most defensible when the representative values reasonably summarize the exposure and biological response over the period relevant to the outcome. That may be plausible when concentrations and activity change little, when the victim's exposure metric is insensitive to the timing of those changes, or when the approximation has been checked against a suitable dynamic model or data. A static calculation can also be valuable as a transparent screen: its assumptions and influential inputs are visible, and it can be evaluated without solving time-dependent equations.

A dynamic simulation is preferable when the question depends on onset, offset, repeated dosing, changing perpetrator concentrations, delayed or reversible biological responses, or a mismatch between the time course of activity and victim exposure. It needs more information, including time-dependent inputs and defensible representations of relevant biological processes. Extra time resolution does not guarantee accuracy: a dynamic model with unsupported parameters or an incorrect structure can be less reliable than a simpler approximation.

Feature	Static approximation	Dynamic simulation
Time representation	One or more representative values	Concentrations and states vary with time
Main strength	Fast, transparent algebraic screening	Represents timing, accumulation, and changing activity
Main assumption	Chosen representative values adequately approximate the relevant response	Time-dependent model and inputs adequately describe the system
Common failure mode	Average or selected values conceal important variation or nonlinear response	Uncertain inputs, parameters, or structure create a precise-looking but unreliable profile

Useful diagnostics begin with the inputs, not merely the final AUC ratio. Plot or tabulate perpetrator concentration and predicted activity over the relevant interval. Check whether a proposed representative value falls within the observed or simulated range, whether the activity is nearly constant, and whether onset or recovery occurs during the victim's exposure. For a dynamic model, inspect concentration and activity trajectories, confirm units and time alignment, and verify that numerical results are stable under a finer time step or solver tolerance. Compare a dynamic result with a static approximation under the same assumptions; a discrepancy is informative, but does not by itself prove which model is correct.

Common pitfall: “The mean activity is 0.750, so using \(q=0.750\) statically must give the same AUC as a dynamic profile whose mean is 0.750.” This is not generally true. Activity affects the concentration trajectory, and AUC integrates that trajectory. As the worked piecewise example showed, the same mean and final concentration can coexist with different intermediate concentrations and different AUCs.

The central limitation of a static result is that it cannot, by itself, describe when an interaction begins or resolves; its summary depends on the representative exposure and response chosen. The central limitation of a dynamic result is dependence on time-resolved data and model structure. The introductory material on drug-drug interaction prediction provides the broad motivation; the present distinction is about how exposure and biological activity are represented. Specific inhibition, induction, and mechanistic equations require their own assumptions and should not be inferred from the generic activity multiplier used here.

Check Your Understanding

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

In the illustrative elimination model, what does \(q(t)=1\) mean, and what does a value below 1 represent?
For the model \(dC/dt=-k_0qC\), why can equal average activity profiles produce different finite-interval AUCs?
A static calculation evaluates a nonlinear activity response at the mean perpetrator concentration. What assumption would be needed to treat that as the mean activity?
Name one situation in which a dynamic simulation is especially useful and one limitation of using it.

Solutions. (1) \(q(t)=1\) represents baseline activity; \(q(t)<1\) represents reduced activity relative to baseline in this simplified model. (2) Earlier activity changes affect concentrations over the remaining interval; AUC integrates the full concentration trajectory, not only the mean activity or the final concentration. (3) The response would need to be approximately linear over the concentration range, or otherwise be known to satisfy that averaging equality for the specific exposure profile. Nonlinearity generally prevents assuming it. (4) A dynamic model is useful when the interaction's onset, offset, or changing activity overlaps the victim's exposure. Its limitation is that it relies on time-resolved inputs and a sufficiently accurate model structure and parameters.