Time-Varying Inhibitor and Inducer Exposure

Why Exposure Time Course Matters

A perpetrator drug may inhibit or induce an enzyme, but its concentration at the enzyme site is rarely constant during repeated dosing. Concentrations rise after each dose, decline between doses, and may accumulate across dosing intervals. The enzyme response does not necessarily follow those changes at the same speed: reversible inhibition can track concentration relatively quickly, whereas enzyme inactivation and induction are governed by enzyme turnover and therefore can persist after concentrations fall.

The central modeling task is to use the time-varying unbound perpetrator concentration, \(C_u(t)\), as an input to one or more enzyme-response equations. The resulting activity profile depends both on the shape of the exposure profile and on the response mechanism. The concentration profile must represent the relevant site and be on the time scale of the enzyme response; a nominal dose or average concentration alone may not capture peak inhibition, delayed induction, or recovery.

Representing Repeated-Dose Exposure

A simulation needs a concentration-time profile rather than just a dosing interval. That profile may be supplied as measured or simulated values, or generated by a drug disposition model. Here, treat \(C_u(t)\) as a known input. A repeated-dose profile typically has a peak after each dose and a trough before the next one, but accumulation can raise both peaks and troughs across successive intervals.

For a simple reversible competitive inhibitor at fixed substrate concentration \(S\), with Michaelis–Menten constant \(K_m\) and inhibition constant \(K_i\), the rate relative to the rate without inhibitor is

Thus, \(K_{i,\mathrm{app}}\) is the apparent inhibition constant under the stated fixed-substrate conditions. It is not generally equal to \(K_i\). If substrate concentration changes over time, the first expression should use the contemporaneous \(S(t)\); a single fixed apparent constant would no longer represent that changing substrate condition. Mixed and noncompetitive mechanisms require their corresponding inhibition expressions rather than this competitive formula.

Slower Changes in Active Enzyme Amount

Let \(E(t)\) denote active enzyme amount relative to its baseline amount, so \(E=1\) is the baseline. Let \(k_{\mathrm{deg}}\) be the first-order enzyme degradation rate constant, in \(\mathrm{h^{-1}}\). In the absence of modulation, synthesis balances degradation and \(E\) remains at 1. With mechanism-based inactivation, define the concentration-dependent inactivation rate as

Here \(k_{\mathrm{inact}}\) has units of \(\mathrm{h^{-1}}\), and \(K_I\) has concentration units. For induction, a common concentration-response representation of enzyme synthesis is \(g_{\mathrm{ind}}(t)=1+E_{\max}C_u(t)/(EC_{50}+C_u(t))\). The factor \(g_{\mathrm{ind}}\) is dimensionless; \(E_{\max}\) is the maximum fractional increase in synthesis above baseline, and \(EC_{50}\) has concentration units.

A simple model allowing inactivation and induction to operate together is

The first term represents synthesis, increased by induction. The loss term includes ordinary enzyme degradation and additional inactivation of active enzyme. In the absence of both effects, \(g_{\mathrm{ind}}=1\) and \(k_{\mathrm{obs}}=0\), so \(E=1\) is the steady state. If only induction is present, the enzyme amount moves toward the increased synthesis target. If only inactivation is present, it falls during exposure and recovers through baseline synthesis after exposure ends.

Deriving the Piecewise-Constant Update

For a time interval over which \(C_u\), and therefore \(g_{\mathrm{ind}}\) and \(k_{\mathrm{obs}}\), are constant, write \(a=k_{\mathrm{deg}}+k_{\mathrm{obs}}\). The abundance equation becomes \(dE/dt=k_{\mathrm{deg}}g_{\mathrm{ind}}-aE\). Its equilibrium value on that interval is \(E_{\mathrm{eq}}=k_{\mathrm{deg}}g_{\mathrm{ind}}/a\). Subtracting this equilibrium from both sides gives \(d(E-E_{\mathrm{eq}})/dt=-a(E-E_{\mathrm{eq}})\). Separating variables and integrating from the start of the interval to its end gives the update:

This expression shows why enzyme response has memory: the amount at the end of an interval depends on the amount at its beginning, not only on the current concentration. For a concentration profile that changes continuously, the same differential equation can be integrated numerically with a time-varying input.

Simulation Choices, Checks, and Limitations

For repeated dosing, divide the input profile into intervals fine enough to represent concentration changes and update the enzyme state at each interval, or solve the differential equation with a numerical integrator. The concentration profile and enzyme state must use compatible time units. For example, if rate constants are in \(\mathrm{h^{-1}}\), interval durations must be in hours. For simultaneous reversible inhibition, calculate \(q_I(t)\) from the appropriate mechanism expression and combine it with \(E(t)\) only when the model assumes that reversible binding changes activity without separately changing enzyme abundance.

A useful diagnostic is to inspect concentration, inhibition factor, and enzyme amount on the same time axis. Reversible inhibition should rise and fall with concentration under the assumed rapid-equilibrium model. Enzyme abundance should evolve smoothly according to its turnover rate; it should not jump at a dose unless the model explicitly includes an instantaneous change. When concentration returns to zero, reversible inhibition disappears in this model, but inactivated enzyme recovers gradually and induced abundance declines gradually. Implausible jumps, negative enzyme amounts, or results that change markedly when the numerical time step is reduced may signal an implementation or time-scale problem.

The major alternatives differ in what they retain about exposure. A single average concentration can be convenient for a steady-state approximation, but generally loses peak and trough information. A constant-inhibitor equilibrium calculation can estimate an instantaneous inhibition factor at a specified concentration, but does not represent enzyme turnover or the history of repeated exposure. A time-resolved model preserves this history at the cost of requiring a defensible concentration profile and additional parameters. It remains a simplification: the equations assume specified concentration-response relationships, first-order enzyme turnover, and adequate representation of the enzyme-site concentration. They do not by themselves establish that the chosen exposure profile or mechanism is correct.