Deriving an Enzyme Abundance Turnover Model

Why Enzyme Induction Needs a Turnover Model

An inducer does not generally make enzyme abundance jump immediately to a new value. Instead, receptor-mediated regulation changes enzyme synthesis, while existing enzyme continues to be degraded. The abundance observed at any time therefore reflects both the current synthesis rate and the history of synthesis and degradation. A turnover model represents these competing processes explicitly.

The Emax model for induction provides a concentration-dependent target for the abundance response. The turnover model adds a dynamic mechanism: inducer concentration controls synthesis, and degradation determines how quickly abundance moves toward that target. The distinction matters whenever inducer exposure changes faster than enzyme abundance can equilibrate.

Definition: Let \(E(t)\) denote enzyme abundance at time \(t\), and let \(C_{\mathrm{ind}}(t)\) be the inducer concentration driving the response. Let \(k_{\mathrm{syn},0}\) be basal enzyme synthesis per unit time and \(k_{\mathrm{deg}}\) be the first-order degradation rate constant. The induction function is \(F(C_{\mathrm{ind}})=1+E_{\max}C_{\mathrm{ind}}/(EC_{50}+C_{\mathrm{ind}})\). Here \(E_{\max}\) is the maximum increase above basal synthesis, expressed as a fraction of basal synthesis, and \(EC_{50}\) is the concentration producing half of that maximal increase. The synthesis rate is \(k_{\mathrm{syn},0}F(C_{\mathrm{ind}})\).

Under this convention, \(F(0)=1\), so the synthesis rate is basal without inducer. At very high inducer concentration, \(F\) approaches \(1+E_{\max}\). Thus if \(E_{\max}=1\), the maximum synthesis rate is twice basal—not one times basal. In some applications, authors define a parameter as the maximum total fold response rather than the maximum increment; the parameterization must be checked before using reported values.

Deriving the Abundance Equation

Assume enzyme is produced at the concentration-dependent synthesis rate and removed by first-order degradation. First-order degradation means that the amount removed per unit time is proportional to the current abundance: the degradation rate is \(k_{\mathrm{deg}}E(t)\). The net rate of change is synthesis minus degradation.

\frac{dE(t)}{dt} = k_{\mathrm{syn},0}F(C_{\mathrm{ind}}(t)) - k_{\mathrm{deg}}E(t)

The terms on the right have units of abundance per time. The first term adds enzyme; the second removes it. In the absence of inducer, a constant basal abundance \(E_0\) is possible when synthesis and degradation balance. Setting the derivative to zero at baseline gives \(k_{\mathrm{syn},0}=k_{\mathrm{deg}}E_0\), or \(E_0=k_{\mathrm{syn},0}/k_{\mathrm{deg}}\). Dividing abundance by \(E_0\) defines the dimensionless relative abundance \(e(t)=E(t)/E_0\), and the model becomes:

\frac{de(t)}{dt} = k_{\mathrm{deg}}\left[F(C_{\mathrm{ind}}(t))-e(t)\right]

This normalized form makes the model's target clear: at a fixed inducer concentration, relative abundance moves toward \(F(C_{\mathrm{ind}})\). The rate constant \(k_{\mathrm{deg}}\) sets how quickly it moves. The model does not say that abundance equals the induction function at every instant.

Solution for a Constant Inducer Concentration

When \(C_{\mathrm{ind}}\) is held constant, \(F\) is constant too. To solve the equation, subtract its equilibrium value \(F\). With \(z(t)=e(t)-F\), the differential equation becomes \(dz/dt=-k_{\mathrm{deg}}z\). Separating variables and integrating gives \(z(t)=z(0)e^{-k_{\mathrm{deg}}t}\). Substituting back yields:

e(t)=F+\left[e(0)-F\right]e^{-k_{\mathrm{deg}}t}

The solution starts at \(e(0)\) and approaches \(F\) exponentially. If abundance begins below the target, it increases; if it begins above the target, it decreases. The equilibrium follows by taking the limit as time increases, since \(e^{-k_{\mathrm{deg}}t}\) tends to zero. In original units, the equilibrium abundance is \(E_0F\).

Derivation for a changing inducer: For time-varying \(F(C_{\mathrm{ind}}(t))\), multiply the normalized equation by the integrating factor \(e^{k_{\mathrm{deg}}t}\). Then \(\frac{d}{dt}\left[e(t)e^{k_{\mathrm{deg}}t}\right]=k_{\mathrm{deg}}F(C_{\mathrm{ind}}(t))e^{k_{\mathrm{deg}}t}\). Integrating from \(0\) to \(t\) and rearranging gives \(e(t)=e(0)e^{-k_{\mathrm{deg}}t}+k_{\mathrm{deg}}\int_0^t e^{-k_{\mathrm{deg}}(t-u)}F(C_{\mathrm{ind}}(u))\,du\). This expression shows that current abundance depends on initial abundance and on past induction, with older inputs weighted less because of subsequent degradation.

Worked Examples

Worked Example: The Induced Steady-State Abundance

Suppose the basal abundance of an enzyme is \(100\) pmol/g tissue. The induction parameters are \(E_{\max}=2.0\) and \(EC_{50}=4.0\) micromol/L, and the inducer concentration is held at \(4.0\) micromol/L long enough to reach steady state. Find the predicted abundance.

First calculate the induction function:

F=1+\frac{(2.0)(4.0)}{4.0+4.0}=1+1.0=2.0

At equilibrium, \(e=F=2.0\). Therefore the predicted abundance is \(E=E_0e=(100\ \mathrm{pmol/g})(2.0)=200\) pmol/g tissue. This is a twofold total abundance relative to baseline. At very high concentration, this parameterization would approach \(1+E_{\max}=3.0\) times baseline, not \(2.0\) times baseline.

Worked Example: Approach Toward a New Abundance

Assume an enzyme begins at baseline, \(e(0)=1.0\). A constant inducer concentration produces \(F=1.75\), and \(k_{\mathrm{deg}}=0.150\ \mathrm{h}^{-1}\). Calculate relative abundance after 6.0 h.

The constant-input solution applies because the concentration, and thus \(F\), is unchanged throughout the interval:

e(6.0)=1.75+(1.0-1.75)e^{-(0.150)(6.0)} =1.75-0.75e^{-0.900}

Since \(e^{-0.900}=0.40657\), \(e(6.0)=1.75-(0.75)(0.40657)=1.44507\), or approximately \(1.445\). The abundance has risen above baseline but remains below the eventual target of \(1.75\). If the baseline abundance were \(80.0\) pmol/g, the predicted abundance at 6.0 h would be \((80.0)(1.44507)=115.606\), approximately \(116\) pmol/g.

Worked Example: Two Consecutive Exposure Intervals

Let \(E_{\max}=1.0\), \(k_{\mathrm{deg}}=0.200\ \mathrm{h}^{-1}\), and \(e(0)=1.0\). During the first 3.0 h, the inducer concentration is \(EC_{50}\); during the next 2.0 h, it is zero. Calculate abundance at 5.0 h.

At \(EC_{50}\), the induction increment is half its maximum, so \(F=1+1.0/2=1.50\). Applying the constant-input solution for the first interval gives:

e(3.0)=1.50+(1.00-1.50)e^{-(0.200)(3.0)} =1.50-0.50e^{-0.600} =1.225594\ldots

During the second interval \(F=1.0\). Carry the unrounded abundance from the first interval into the next calculation:

e(5.0)=1.0+(1.225594\ldots-1.0)e^{-(0.200)(2.0)} =1.151220\ldots\approx1.151

Abundance does not return immediately to baseline when inducer concentration falls to zero. It declines toward the basal target, and the predicted value at 5.0 h is about 1.151 times baseline.

Assumptions, Diagnostics, and Alternatives

The derivation relies on a defined set of assumptions. First, synthesis is proportional to basal synthesis multiplied by the specified induction function. Second, degradation is first-order with a constant \(k_{\mathrm{deg}}\), so the same rate constant applies across the abundance range and exposure conditions. Third, the measured inducer concentration is an appropriate driver of the synthesis response, with no additional delay or unrepresented compartment between concentration and enzyme regulation. Fourth, the baseline is at steady state when \(E_0=k_{\mathrm{syn},0}/k_{\mathrm{deg}}\) is used. The concentration units in \(C_{\mathrm{ind}}\) and \(EC_{50}\) must match, and time units for \(k_{\mathrm{deg}}\) and the time course must also match.

For a constant concentration, the model predicts a monotonic exponential approach to a single target. A time course with an overshoot, oscillation, or a persistent plateau inconsistent with the calculated target signals that at least one assumption may be inadequate. Other warning signs include different apparent degradation rates across conditions, an estimated negative synthesis rate, or an estimated \(EC_{50}\) outside a plausible concentration range. Before attributing such patterns to a more complex mechanism, check concentration measurements, baseline stability, units, and whether the exposure actually remained constant.

Approach	What it represents	Main limitation or strength
Instantaneous Emax response	Abundance is set directly to the concentration-dependent target at each time.	Simple, but it assumes no delay between inducer concentration and abundance.
Synthesis-and-degradation turnover model	Induction changes synthesis; degradation moves abundance toward its target over time.	Represents enzyme turnover and exposure history, but requires turnover information.
Empirical delayed or multi-process model	Additional steps or response delays are represented explicitly or phenomenologically.	May capture more complex time courses, but adds parameters and requires data to support them.

The turnover formulation is useful when the enzyme's abundance cannot plausibly track concentration instantaneously. Its limitations are equally important: a single well-mixed abundance state cannot represent multiple enzyme pools, delayed transcription and translation as separate processes, or concentration-dependent degradation unless those features are added. The equation predicts abundance, not necessarily catalytic activity. For example, an inhibitor can alter activity without changing abundance, so activity and abundance should not be treated as interchangeable outcomes.

Common Pitfall

A common error is to set \(e(t)=F(C_{\mathrm{ind}}(t))\) at every time and call that the turnover prediction. This is the instantaneous Emax approximation: it discards the degradation process and assumes abundance reaches its target immediately. The turnover equation instead gives \(de/dt=k_{\mathrm{deg}}(F-e)\); unless abundance already equals the target or the system is at equilibrium, the derivative is generally nonzero. Use the constant-input solution or the time-varying integral when predicting abundance over time.

Key takeaway: Inducer concentration sets the synthesis target through \(F\), while enzyme degradation governs how abundance approaches that target. The target is not the time-varying abundance itself.

Check Your Understanding

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

If \(E_{\max}=1.2\) and \(C_{\mathrm{ind}}=EC_{50}\), what induction function \(F\) does this convention imply?
If \(E_{\max}=1.0\), \(EC_{50}=2.0\) micromol/L, and \(C_{\mathrm{ind}}=6.0\) micromol/L, calculate \(F\).
An enzyme starts at \(e(0)=1.20\). With no inducer, \(k_{\mathrm{deg}}=0.250\ \mathrm{h}^{-1}\). What is \(e(2.0\ \mathrm{h})\)?
Why is setting abundance equal to \(F\) immediately after a concentration change generally not the turnover-model prediction?
Name one assumption whose failure could produce a time course that does not approach the model's constant-input target exponentially.

Solutions. 1. At \(EC_{50}\), the induction increment is half of \(E_{\max}\), so \(F=1+1.2/2=1.60\). 2. \(F=1+(1.0)(6.0)/(2.0+6.0)=1.75\). 3. The target is \(1.0\), so \(e(2.0)=1.0+(1.20-1.0)e^{-(0.250)(2.0)}=1+0.20e^{-0.50}=1.121306\ldots\), approximately \(1.121\). 4. Because synthesis changes the target, while existing enzyme is removed at a finite rate; the abundance follows the turnover equation rather than jumping instantly. 5. Examples include non-first-order degradation, an unrepresented delay between inducer concentration and synthesis, or a baseline that was not at steady state.