Why Enzyme Turnover Changes the Inactivation Time Course
Irreversible inactivation removes active enzyme molecules, but it does not stop the body from making new enzyme. The observed activity therefore reflects two processes at once: loss of active enzyme through degradation and inhibitor-mediated inactivation, and replacement through enzyme synthesis. A model that omits turnover can mistake gradual replacement for weak inactivation, while a model that assumes enzyme amount instantly reaches equilibrium misses the time needed to change.
The model below treats inhibitor concentration as an input to an enzyme-amount differential equation. Under constant concentration, the equation has a closed-form solution. When concentration changes over time, the instantaneous inactivation rate changes too, and enzyme amount responds with a delay governed by turnover and inactivation together.
Here \(k_{\mathrm{syn}}\) has units of enzyme amount per time, while \(k_{\mathrm{deg}}\), \(k_{\mathrm{inact}}\), and \(k_{\mathrm{obs}}\) have units of inverse time. \(I(t)\) and \(K_I\) must use the same concentration units. This model assumes that inactivation is irreversible over the modeled interval, that the measured activity is proportional to active enzyme amount, and that synthesis and degradation follow the stated forms. It also treats \(I(t)\) as the relevant inhibitor concentration at the enzyme site; if only plasma concentration is available, an exposure or effect-site model may be needed to connect it to that concentration.
In the absence of inhibitor, \(k_{\mathrm{obs}}=0\). If the enzyme has been at its uninhibited steady state, then \(E_{\mathrm{base}}=k_{\mathrm{syn}}/k_{\mathrm{deg}}\). It is often useful to report activity relative to that baseline, \(e(t)=E(t)/E_{\mathrm{base}}\). This normalization is valid when \(E_{\mathrm{base}}>0\), and gives baseline activity \(e=1\).
Constant Inhibitor Concentration: Steady State and Approach to It
If \(I(t)=I\) is constant, then \(k_{\mathrm{obs}}\) is constant as well. The equation becomes a linear differential equation with constant coefficients. Its steady state is the amount at which synthesis exactly balances degradation plus inactivation:
The steady-state fraction of baseline activity is therefore \(E_{\mathrm{ss}}/E_{\mathrm{base}}=k_{\mathrm{deg}}/(k_{\mathrm{deg}}+k_{\mathrm{obs}})\). This shows why an inhibitor's effect cannot be interpreted from \(k_{\mathrm{obs}}\) alone: the same inactivation rate leaves a higher fraction of activity when baseline enzyme turnover is faster.
The time course follows by rewriting the differential equation around its steady state. The following derivation also makes explicit the special case where the initial amount already equals that steady state.
For constant concentration, let \(\lambda=k_{\mathrm{deg}}+k_{\mathrm{obs}}\). Since \(k_{\mathrm{syn}}=\lambda E_{\mathrm{ss}}\), the ODE is \(dE/dt=-\lambda(E-E_{\mathrm{ss}})\).
If \(E(0)\ne E_{\mathrm{ss}}\), separation of variables gives \(\int_{E(0)}^{E(t)}dE/(E-E_{\mathrm{ss}})=-\lambda\int_0^t du\). Thus \(\ln|E(t)-E_{\mathrm{ss}}|-\ln|E(0)-E_{\mathrm{ss}}|=-\lambda t\), and \(E(t)-E_{\mathrm{ss}}=[E(0)-E_{\mathrm{ss}}]e^{-\lambda t}\).
If \(E(0)=E_{\mathrm{ss}}\), the ODE gives \(dE/dt=0\) at that amount. The constant function \(E(t)=E_{\mathrm{ss}}\) solves the equation. Consequently, the same solution formula below applies both when the initial amount differs from steady state and when it equals steady state.
The exponential time constant is \(1/(k_{\mathrm{deg}}+k_{\mathrm{obs}})\), and the half-time for the difference from the inhibitor-dependent steady state to halve is \(\ln(2)/(k_{\mathrm{deg}}+k_{\mathrm{obs}})\). These are properties of the approach to \(E_{\mathrm{ss}}\), not necessarily the enzyme's uninhibited half-life.
Worked Example: Activity During Constant Inhibitor Exposure
Suppose a hypothetical enzyme has an uninhibited amount of 100 nmol, a degradation half-life of 24 h, and is exposed to a constant inhibitor concentration of 2 µmol/L. Let \(K_I=1\) µmol/L and \(k_{\mathrm{inact}}=0.120\ \mathrm{h}^{-1}\). Find activity relative to baseline at 8 h, assuming enzyme amount is initially at its uninhibited steady state.
First convert the half-life to the degradation rate constant and calculate synthesis needed to maintain 100 nmol at baseline:
The concentration-dependent inactivation rate, total loss rate, and inhibitor-dependent steady state are
Using the time-course solution at 8 h:
Relative activity is \(57.27/100=0.5727\), or about 57.3% of baseline at 8 h. It has not yet reached the eventual steady-state fraction, \(26.52/100=26.5\%\), because the enzyme amount changes over time rather than instantaneously.
What \(k_{\mathrm{inact}}\) and \(K_I\) Mean
In this model, \(k_{\mathrm{inact}}\) is the limiting inactivation rate as inhibitor concentration becomes very large: as \(I\) increases without bound, \(k_{\mathrm{obs}}\) approaches \(k_{\mathrm{inact}}\). \(K_I\) is the concentration at which \(k_{\mathrm{obs}}=k_{\mathrm{inact}}/2\). Thus, \(K_I\) describes the concentration scale over which inactivation rate rises, whereas \(k_{\mathrm{inact}}\) sets its maximum. Neither parameter alone determines the activity time course; that also depends on exposure over time and enzyme turnover.
At concentrations much lower than \(K_I\), \(k_{\mathrm{obs}}\) is approximately \((k_{\mathrm{inact}}/K_I)I\), so the ratio \(k_{\mathrm{inact}}/K_I\) governs the low-concentration slope. At concentrations comparable to or greater than \(K_I\), the rate is nonlinear in concentration and approaches saturation. Two compounds with the same ratio can therefore have similar low-concentration behavior but different rates at higher concentrations.
Worked Example: Same Low-Concentration Slope, Different Saturation
Compare two hypothetical inactivation models: A has \(k_{\mathrm{inact}}=0.120\ \mathrm{h}^{-1}\) and \(K_I=3.00\) µmol/L; B has \(k_{\mathrm{inact}}=0.240\ \mathrm{h}^{-1}\) and \(K_I=6.00\) µmol/L. Both have \(k_{\mathrm{inact}}/K_I=0.0400\ \mathrm{L/(µmol\cdot h)}\). Calculate their rates at \(I=1.00\) and \(30.0\) µmol/L.
At 1.00 µmol/L:
At 30.0 µmol/L:
The rates are close at the lower concentration, consistent with the similar low-concentration slopes, but differ substantially at 30.0 µmol/L because the two models have different \(K_I\) and maximal rates. If \(k_{\mathrm{deg}}=0.0500\ \mathrm{h}^{-1}\), their steady-state fractions at 30.0 µmol/L would be \(0.0500/(0.0500+0.1091)=0.314\) for A and \(0.0500/(0.0500+0.2000)=0.200\) for B. Those fractions describe steady states, not activity at a specified finite time.
Changing Concentrations and Interval-by-Interval Calculation
When inhibitor concentration changes, \(k_{\mathrm{obs}}(t)\) changes with it. The enzyme amount cannot generally be calculated by substituting the current concentration into the constant-concentration steady-state formula: that formula gives the steady state for a concentration held constant, not the enzyme amount immediately after a concentration change. For a known time-varying concentration, the linear ODE has the integrating-factor solution below.
This expression accounts for enzyme present at time zero and enzyme synthesized throughout the interval, with each contribution reduced by subsequent loss. For piecewise-constant concentrations, it is simpler to apply the constant-concentration solution successively, carrying the ending enzyme amount in one interval forward as the initial amount in the next.
Worked Example: Two Consecutive Inhibitor Concentrations
Consider a hypothetical enzyme with \(k_{\mathrm{deg}}=0.0500\ \mathrm{h}^{-1}\), baseline amount 100 nmol, and \(k_{\mathrm{syn}}=5.00\ \mathrm{nmol/h}\). Let \(k_{\mathrm{inact}}=0.100\ \mathrm{h}^{-1}\) and \(K_I=2.00\) µmol/L. The enzyme begins at baseline. The inhibitor concentration is \(2.00\) µmol/L for 4 h, then \(6.00\) µmol/L for a further 4 h. Find the amount at 8 h.
In the first interval \(I=K_I\), so \(k_{\mathrm{obs}}=0.100(2.00)/(2.00+2.00)=0.0500\ \mathrm{h}^{-1}\). The combined loss rate is \(0.100\ \mathrm{h}^{-1}\), and the interval-specific steady state is \(5.00/0.100=50.0\) nmol. Therefore,
In the second interval \(I=3K_I\), giving \(k_{\mathrm{obs}}=0.100(6.00)/(2.00+6.00)=0.07500\ \mathrm{h}^{-1}\). The combined loss rate is \(0.1250\ \mathrm{h}^{-1}\), and the new interval-specific steady state is \(5.00/0.1250=40.00\) nmol. Use the amount at 4 h as the initial amount for this interval:
The activity at 8 h is 66.4% of the 100 nmol baseline. It is not the steady state for either interval: the concentration changed before the first interval reached its steady state, and the second interval also lasted a finite time.
Assumptions, Diagnostics, and Common Pitfalls
The model is useful when the measured endpoint tracks active enzyme amount and exposure drives a saturable, irreversible loss process. Its assumptions should be checked against the data and experimental setting. In particular, enzyme synthesis is taken as constant, degradation is first-order, and the inactivation rate follows the specified saturable concentration relationship. If inhibitor concentration varies rapidly, uncertain tissue exposure or delayed access to the enzyme can make a concentration-driven model appear to have a different \(k_{\mathrm{inact}}\) or \(K_I\).
Useful diagnostics include plotting observed activity and model predictions against time, checking whether residuals show systematic early-versus-late deviations, and comparing predicted and observed activity across concentration levels. A persistent mismatch in the time course may indicate that the assumed turnover rate or synthesis is wrong, the concentration at the enzyme site is not represented well, the measured signal is not proportional to active enzyme, or the inactivation mechanism is inadequate. Data at only one concentration may be especially uninformative about both \(k_{\mathrm{inact}}\) and \(K_I\): if all observations lie far below \(K_I\), they primarily inform their ratio; if concentrations do not span a broad range, saturation may be poorly characterized.
A common mistake is to set enzyme amount equal to \(E_{\mathrm{ss}}\) as soon as the inhibitor is given. That confuses a long-run balance with a finite-time response. The correct finite-time amount includes the exponential transient from the starting value. Another mistake is to interpret \(K_I\) as the concentration that halves total enzyme activity. It halves the instantaneous inactivation rate relative to \(k_{\mathrm{inact}}\); the resulting activity also depends on \(k_{\mathrm{deg}}\), synthesis, exposure duration, and initial enzyme amount.
A simpler model with no turnover may be reasonable over a short interval if synthesis and degradation are negligible on that time scale, but it cannot describe long-term replacement of enzyme. At the other extreme, assuming enzyme is always at its concentration-specific steady state discards the time-course transient. Earlier material on slow-binding inhibition versus rapid-equilibrium inhibition concerns reversible binding and equilibration; the present model instead describes loss of active enzyme with synthesis and degradation. This formulation is also limited if synthesis itself changes with treatment, if there are multiple enzyme pools with different turnover, or if inactivation is reversible.
Check Your Understanding
Use the model definitions above. The solutions are included so that each calculation can be checked.
- If \(k_{\mathrm{inact}}=0.0800\ \mathrm{h}^{-1}\), \(K_I=4.00\) µmol/L, and \(I=4.00\) µmol/L, what is \(k_{\mathrm{obs}}\)?
- If \(k_{\mathrm{syn}}=3.00\) nmol/h and \(k_{\mathrm{deg}}=0.0300\ \mathrm{h}^{-1}\), what is baseline enzyme amount? If a constant inhibitor produces \(k_{\mathrm{obs}}=0.0600\ \mathrm{h}^{-1}\), what is the inhibitor-dependent steady-state fraction of baseline?
- For constant concentration, what happens to the solution if \(E(0)=E_{\mathrm{ss}}\), and why is logarithmic separation not needed for this case?
- Why can two inhibitors with the same \(k_{\mathrm{inact}}/K_I\) have different rates at high concentration?
Solutions. 1. \(k_{\mathrm{obs}}=0.0800(4.00)/(4.00+4.00)=0.0400\ \mathrm{h}^{-1}\). 2. \(E_{\mathrm{base}}=3.00/0.0300=100\) nmol. The steady-state fraction is \(k_{\mathrm{deg}}/(k_{\mathrm{deg}}+k_{\mathrm{obs}})=0.0300/(0.0300+0.0600)=0.333\), or 33.3%. 3. The constant solution is \(E(t)=E_{\mathrm{ss}}\) for all \(t\), since the ODE derivative is zero at that amount. The logarithmic derivation divides by \(E-E_{\mathrm{ss}}\), so it applies only when the initial amount differs from steady state. 4. The ratio determines the low-concentration slope, but high-concentration rates depend on saturation and approach each inhibitor's own \(k_{\mathrm{inact}}\), which may differ.