From Static Concentrations to Dynamic Exposure Profiles

From a Dose to a Time-Varying Perpetrator

A static drug–drug interaction (DDI) calculation evaluates an effect at a selected perpetrator concentration. A dynamic calculation instead follows the perpetrator concentration through time and updates inhibition or induction as that concentration changes. The distinction matters when exposure rises and falls substantially within a dosing interval, when inhibition responds quickly, or when induction develops and resolves gradually.

The first step is not an inhibition equation. It is a concentration model consistent with the perpetrator’s dose, route, disposition, and the site where the relevant enzyme is located. In a PBPK model, that concentration may be an unbound concentration in a liver or intestinal compartment. A simpler pharmacokinetic model can be used when its assumptions are appropriate, but a plasma concentration should not automatically be treated as the concentration at an intracellular enzyme.

A useful starting point is a one-compartment model with first-order absorption and elimination after a single oral dose. Let \(D\) be the dose, \(F\) the fraction reaching systemic circulation, \(V\) the apparent volume of distribution, \(k_a\) the absorption rate constant, and \(k_e=CL/V\) the elimination rate constant. Assume the body contains no perpetrator before dosing, absorption and elimination are first order, and the model’s concentration is the concentration used for the effect calculation. Then, for \(k_a\ne k_e\),

This expression describes the rise from absorption and the later decline from elimination. It is not a universal PBPK concentration equation: tissue distribution, transport, nonlinear clearance, and other processes may require a more detailed model. The concentration must have consistent units; for example, if \(D\) is in milligrams and \(V\) is in liters, \(C_p(t)\) is in milligrams per liter.

Mapping Concentration to Reversible Inhibition

For reversible inhibition, assume that the unbound concentration at the enzyme site determines inhibition and that the interaction can be represented by a simple concentration–potency relationship. Let \(K_i\) be the inhibition constant in the same concentration units as \(C_{p,u}(t)\). Under this simplified relationship, the fraction of baseline enzyme activity remaining at time \(t\) is

This equation is a time-dependent version of the usual static reversible-inhibition expression: the chosen concentration is now a function of time. It assumes rapid reversible binding relative to the time scale on which concentrations change, a suitable competitive-inhibition representation, and no other mechanism changing enzyme abundance or activity. The equation does not itself predict the perpetrator concentration. If concentration at the enzyme site is unavailable, substituting total plasma concentration requires an explicit justification rather than an automatic assumption.

Induction Has Its Own Time Scale

Unlike the rapid reversible-inhibition approximation, induction usually requires a change in enzyme abundance. A standard turnover representation makes that delay explicit. Let \(R(t)\) be enzyme abundance relative to its baseline value, so \(R=1\) means baseline abundance. Let \(k_{\mathrm{deg}}\) be the first-order enzyme degradation rate constant, and let \(E_{\max}\) and \(EC_{50}\) describe the maximum induction and concentration producing half of that maximum, respectively. Assuming that perpetrator concentration stimulates enzyme synthesis and that degradation remains first order, one model is

The baseline is consistent with this equation: when \(C_{p,u}(t)=0\), synthesis and degradation balance at \(R=1\). When concentration rises, synthesis increases, but abundance does not jump immediately to its new target. This formulation assumes the chosen concentration drives synthesis through the stated saturable function and that the turnover rate is constant. Other induction mechanisms or models may require different equations.

Putting the Profiles Together

When the same enzyme pathway is affected by both reversible inhibition and induction, a simple combined representation is to multiply enzyme abundance by the reversible-activity fraction. If baseline pathway clearance is \(CL_{\mathrm{path},0}\), then under the assumptions that the induction and inhibition effects act independently in this form,

This is a time-varying pathway clearance, not necessarily total victim clearance. Other pathways may need to be represented separately and added according to the victim model. To predict victim concentrations or exposure, the simulation must use the changing clearance in the victim’s disposition equations and integrate those equations through time. A single concentration or a simple average of the clearance multipliers generally cannot preserve the timing of the interaction.

The concentration and effect equations can be solved together: the dosing input drives perpetrator disposition; \(C_{p,u}(t)\) drives reversible inhibition directly and drives enzyme synthesis in the turnover equation; the resulting pathway activity modifies victim disposition. In a PBPK model, each connection should refer to the appropriate compartment and unbound concentration. The model should also use the same concentration units for \(C_{p,u}\), \(K_i\), and \(EC_{50}\).

Derivation, Diagnostics, and Alternatives

The induction solution used above follows directly from the constant-concentration equation. Write \(R_{\mathrm{target}}=1+E_{\max}C_{p,u}/(EC_{50}+C_{p,u})\). Then \(dR/dt=k_{\mathrm{deg}}(R_{\mathrm{target}}-R)\). Subtracting the constant target gives \(d(R-R_{\mathrm{target}})/dt=-k_{\mathrm{deg}}(R-R_{\mathrm{target}})\), whose solution is \(R(t)-R_{\mathrm{target}}=[R(0)-R_{\mathrm{target}}]e^{-k_{\mathrm{deg}}t}\). Rearranging yields the expression used in the example. Since \(k_{\mathrm{deg}}>0\), the exponential term decays toward zero and \(R(t)\) approaches the target rather than changing instantaneously.

Several checks help identify a broken or implausible implementation:

• Check the concentration profile. Confirm that dose, bioavailability, volume, clearance, absorption, route, and units produce a plausible curve. In the simple oral model, concentration is zero at time zero, rises during absorption, and eventually declines.
• Check the effect bounds. With nonnegative concentration and positive \(K_i\), \(f_{\mathrm{act}}(t)\) must lie between zero and one. Under the stated induction function with nonnegative \(E_{\max}\), the synthesis target must lie between 1 and \(1+E_{\max}\).
• Check timing. Reversible inhibition follows the concentration directly in this approximation. Induction should respond with a delay governed by enzyme turnover; it should not instantaneously attain the concentration-dependent target.
• Check the concentration used. A total plasma concentration substituted for unbound concentration at the site may bias the predicted effect, especially when binding or tissue distribution is important.
• Check numerical resolution. Time steps that are too large can miss a concentration peak or distort the timing of changing effects. Compare results after reducing the step size or tightening the solver tolerance.

A static approach is simpler: it evaluates an effect at one selected concentration, such as a specified maximum or average. That may suit screening when concentration changes little or when a conservative summary is the intended purpose, but it discards the within-interval pattern. A dynamic compartmental or PBPK approach represents concentration and effects over time, at the cost of needing more parameters and careful validation. Between these choices, interpolation of measured concentration data can provide a time-varying input when data are available, but it does not by itself explain concentrations at unmeasured times or predict a new dosing condition.

Dynamic representation has limitations. It cannot compensate for an incorrect perpetrator disposition model, an unsuitable inhibition mechanism, or poorly identified potency and turnover parameters. The enzyme-site concentration may not be directly measured, and different combinations of concentration and potency can yield similar predicted effects. The equations here are illustrative models, not a guarantee that a particular biological system follows simple inhibition and turnover relationships.

Common Mistake

A frequent error is to calculate an induction multiplier from the current concentration and apply it immediately as the enzyme abundance. That treats a synthesis target as if it were the actual abundance at that instant. The correction is to distinguish the target from \(R(t)\) and represent the time required for synthesis and degradation to move abundance toward the target. A related error is to use one concentration value throughout a simulation while describing the result as dynamic; a time-varying effect requires a time-varying concentration input or an explicit reason that concentration is effectively constant.