Why Induction Requires a Turnover Model
An increase in enzyme activity during treatment need not mean that each enzyme molecule has become more active. Exposure can instead increase the rate at which enzyme is synthesized. Existing enzyme does not disappear immediately, and newly synthesized enzyme takes time to accumulate. Consequently, enzyme abundance—and any pharmacodynamic capacity that depends on it—may change more slowly than drug concentration.
A turnover model makes this delay explicit. It links drug concentration to enzyme synthesis, while enzyme degradation removes enzyme continuously. The resulting enzyme time course depends both on the exposure history and on how quickly the enzyme turns over.
The synthesis term has units of enzyme amount per time, while \(k_{\mathrm{deg}}\) has units of inverse time. Thus both terms in the differential equation have units of enzyme amount per time. In this model, concentration changes the input of newly synthesized enzyme; it does not directly multiply the amount of enzyme already present.
The model assumes that baseline synthesis and first-order degradation are adequate descriptions over the modeled range, that induction responds to the specified concentration function, and that other processes do not change enzyme abundance. The concentration-response function is saturable: at low concentration, induction rises approximately in proportion to concentration; at high concentration, it approaches \(I_{\max}\).
Baseline, Steady State, and Delay
With no inducing exposure, \(I(C)=0\). At baseline steady state, synthesis balances degradation:
where \(E_0\) is baseline enzyme abundance. This relation assumes the enzyme was at baseline steady state before the exposure of interest. It is often convenient to normalize enzyme abundance as \(R(t)=E(t)/E_0\). Substituting \(E(t)=E_0R(t)\) and \(k_{\mathrm{syn},0}=k_{\mathrm{deg}}E_0\) gives
For constant concentration, \(I(C)\) is constant. At steady state, the rate of change is zero, so \(R_{\mathrm{ss}}=1+I(C)\). The approach to that steady state follows by separating the deviation from steady state: if \(Z(t)=R(t)-R_{\mathrm{ss}}\), then \(dZ/dt=-k_{\mathrm{deg}}Z\), and therefore \(Z(t)=Z(0)e^{-k_{\mathrm{deg}}t}\). Starting at baseline, \(R(0)=1\), the solution is
This derivation shows why induction is delayed: enzyme abundance approaches its exposure-specific steady state exponentially rather than jumping to it. The enzyme turnover half-life is \(t_{1/2}=\ln(2)/k_{\mathrm{deg}}\). After one half-life at a fixed exposure, the enzyme has completed half the change from its initial value to its new steady state.
Worked Example: Delayed Accumulation During Constant Exposure
Suppose a hypothetical enzyme has \(I_{\max}=1.5\), \(EC_{50}=2\) mg/L, and \(k_{\mathrm{deg}}=0.100\) day\(^{-1}\). Assume it begins at baseline, and concentration is held at \(C=8\) mg/L for 7 days. Calculate the fractional induction, the steady-state enzyme abundance relative to baseline, and the abundance after 7 days.
First calculate the exposure-dependent fractional increase in synthesis:
At this constant concentration, the relative steady-state abundance is \(R_{\mathrm{ss}}=1+I(8)=2.20\). The turnover half-life is \(\ln(2)/0.100=6.93\) days, so 7 days is approximately one half-life. Using the time-course solution:
Thus, after 7 days, enzyme abundance is about 1.60 times baseline, although the exposure-specific steady state is 2.20 times baseline. These are relative abundances, not concentrations in mg/L. If catalytic activity per enzyme molecule is unchanged, enzyme-dependent capacity proportional to enzyme abundance would also be about 1.60 times its baseline value at day 7.
Changing Exposure and Return Toward Baseline
For a concentration that varies over time, the enzyme response reflects prior as well as current exposure. With \(R(0)=1\), the normalized solution is
The exponential factor gives more weight to recent induction than to induction further in the past. Earlier increases in synthesis still contribute, but their contribution diminishes as enzyme is degraded. The expression also makes clear why the current concentration alone does not generally determine current enzyme abundance.
If exposure stops abruptly at time \(t_w\), and there is no induction after that time, then \(I(C)=0\) during washout. The normalized equation reduces to \(dR/dt=k_{\mathrm{deg}}(1-R)\). Starting from \(R(t_w)\), its solution is
The enzyme returns toward baseline; it does not immediately return to baseline when drug concentration falls. This washout equation assumes induction has ceased and degradation remains first order at the same \(k_{\mathrm{deg}}\). If concentration declines gradually rather than stopping abruptly, residual induction continues to contribute, and the time-varying exposure solution is needed.
Worked Example: Enzyme Washout After Exposure Stops
Use the hypothetical results from the preceding example: after 7 days of constant exposure, \(R(7)=1.604\), and \(k_{\mathrm{deg}}=0.100\) day\(^{-1}\). Suppose concentration then falls immediately to zero, so \(I(C)=0\). Calculate relative enzyme abundance 7 days after exposure stops.
The washout interval is 7 days. Substituting into the washout solution gives
The enzyme abundance is therefore about 1.30 times baseline after 7 days of washout. The same half-life that governed the approach toward the induced steady state governs the decay of the excess abundance above baseline in this model. The calculation does not imply that drug concentration itself has a 6.93-day half-life; that is the enzyme turnover half-life specified by \(k_{\mathrm{deg}}\).
Worked Example: Two Exposure Levels in Sequence
Consider a hypothetical enzyme with \(I_{\max}=1.2\), \(EC_{50}=1\) mg/L, and \(k_{\mathrm{deg}}=0.200\) day\(^{-1}\). It begins at baseline. Concentration is 2 mg/L for 3 days, then 0.5 mg/L for 4 days. Calculate its relative abundance at the end of the 7-day sequence.
During the first interval,
Applying the constant-exposure solution for 3 days gives the abundance at the switch:
For the second interval, \(I(0.5)=1.2(0.5)/(1+0.5)=0.400\), so the new steady-state abundance is \(1.400\). Starting from \(R(3)=1.361\), the solution over the next 4 days is
The relative abundance at day 7 is approximately 1.382. It has not reached the second interval’s steady state of 1.400, and it also retains the effect of the first interval. Calculating abundance from the day-7 concentration alone would miss that exposure history.
Connecting Enzyme Abundance to Pharmacodynamic Capacity
If enzyme molecules have unchanged catalytic properties, a simple link to maximal catalytic capacity is \(V_{\max}(t)=k_{\mathrm{cat}}E(t)\), where \(k_{\mathrm{cat}}\) is the catalytic rate per enzyme amount. Under that assumption, \(V_{\max}(t)/V_{\max,0}=E(t)/E_0=R(t)\). This is a consequence of enzyme abundance changing while \(k_{\mathrm{cat}}\) remains constant; it is not a general claim that every measured response changes in direct proportion to enzyme abundance.
A response may also depend on substrate availability, other enzymes, or downstream processes. If those factors vary, the relationship between \(R(t)\) and the measured endpoint needs its own model. In particular, enzyme activity can change without a corresponding change in enzyme abundance, so an activity measurement by itself does not establish that synthesis has increased.
Model Choice, Diagnostics, and Limitations
The induction model is appropriate when the proposed mechanism is increased enzyme synthesis with enzyme turnover. It differs from reversible activation through direct modulation of catalytic capacity: a direct effect can track concentration rapidly, whereas abundance induction generally lags and can persist after concentration falls. It also differs from a model in which drug exposure removes active enzyme; that mechanism changes the enzyme balance in a different direction. Choosing among these mechanisms requires observations that distinguish their expected time courses, not merely a good fit to one endpoint.
Useful checks include whether baseline observations support \(E_0\), whether measured or predicted concentrations plausibly drive the modeled synthesis change, and whether enzyme abundance or activity observations show the expected delay and washout. Compare predictions with data during both rising and falling exposure. Systematic residual patterns, an incorrect predicted lag, or persistent discrepancies during washout can indicate that the concentration-response function, turnover assumption, or mechanism is inadequate. Sparse sampling around exposure changes may make these features difficult to evaluate.
The model has important limits. A saturable \(I(C)\) function is a chosen concentration-response form, not proof that synthesis follows that exact shape. The model assumes one enzyme pool with one degradation rate, and it does not represent delays in the induction signal itself, multiple enzyme states, or time-varying degradation. It also assumes that the concentration supplied to the model is the relevant concentration at the enzyme site. If those assumptions are implausible, a more detailed model may be needed. Additional parameters or compartments should be supported by informative data rather than added solely to improve fit.
Check Your Understanding
Answer the questions, then compare your reasoning with the solutions below.
- At baseline steady state, \(k_{\mathrm{syn},0}=12\) amount/day and \(k_{\mathrm{deg}}=0.3\) day\(^{-1}\). What is \(E_0\), and what are its units?
- An enzyme has \(k_{\mathrm{deg}}=0.1\) day\(^{-1}\). What fraction of the original excess above baseline remains after one half-life of washout?
- If \(I_{\max}=2\) and \(C=EC_{50}\), what is \(I(C)\), and what is the relative induced steady state?
- Why can enzyme abundance remain elevated after concentration has fallen to zero?
Solutions. (1) \(E_0=k_{\mathrm{syn},0}/k_{\mathrm{deg}}=12/0.3=40\) amount, because (amount/day)/(1/day) has units of amount. (2) The half-life is \(\ln(2)/0.1=6.93\) days; after one half-life, \(e^{-k_{\mathrm{deg}}t}=1/2\), so half the excess above baseline remains. (3) \(I(EC_{50})=2EC_{50}/(EC_{50}+EC_{50})=1\), a fractional increase of 1, or 100%; the relative steady state is \(1+I=2\). (4) When concentration reaches zero, the increased synthesis term ends in this model, but the enzyme already present is removed at a finite rate. Its excess therefore decays toward baseline rather than vanishing instantly.