Proenzyme Activation With Finite Precursor Pools

Why a Finite Precursor Pool Changes the Activation Trajectory

An increase in enzyme activity need not come from faster catalytic function of enzyme that is already active. It can also arise when a proenzyme is converted into active enzyme. If the proenzyme pool is finite, however, activation consumes that pool: the source of new active enzyme weakens over time. At the same time, active enzyme may be cleared. These two processes together determine whether activity rises, falls, or first rises and then falls.

The key distinction from a constant-input model is that formation is not sustained indefinitely. The precursor itself is a state variable, and its depletion must be represented alongside active-enzyme turnover. The resulting trajectory can be derived directly from the coupled differential equations.

Definition: Let $P(t)$ be the amount or concentration of proenzyme at time $t$, and let $E(t)$ be the amount or concentration of active enzyme. A finite-pool activation model assumes first-order conversion of precursor to active enzyme at rate constant $k_{\mathrm{act}}$, and first-order clearance of active enzyme at rate constant $k_{\mathrm{el}}$: $$ \frac{dP}{dt}=-k_{\mathrm{act}}P,\qquad \frac{dE}{dt}=k_{\mathrm{act}}P-k_{\mathrm{el}}E, $$ with initial values $P(0)=P_0$ and $E(0)=E_0$. Both rate constants are nonnegative and have units of inverse time.

Assumptions and Analytical Solution

The model treats activation as a one-to-one conversion: each unit of precursor converted produces one unit of active enzyme. The conversion and clearance rate constants are constant over the time interval considered. The precursor has no replenishment and no loss route other than activation, while active enzyme has no input other than conversion and is removed by first-order clearance. Amounts can be replaced by concentrations if the relevant volume is fixed. If measured catalytic activity is proportional to active-enzyme amount under fixed assay conditions, its time course has the same shape as $E(t)$.

The precursor equation is independent of $E(t)$, so it can be solved first. Substitution of its solution into the active-enzyme equation then gives the activity trajectory.

Solve precursor depletion.
Separating variables and applying $P(0)=P_0$ gives $P(t)=P_0e^{-k_{\mathrm{act}}t}$.

Apply an integrating factor.
Multiply the active-enzyme equation by $e^{k_{\mathrm{el}}t}$. The product rule gives $\frac{d}{dt}(e^{k_{\mathrm{el}}t}E(t))=k_{\mathrm{act}}P_0e^{(k_{\mathrm{el}}-k_{\mathrm{act}})t}$.

Integrate and use the initial condition.
Integrating from 0 to $t$ yields the expressions below. When the two rate constants are equal, the integral is $t$, rather than a quotient with a zero denominator.

E(t)= \begin{cases} E_0e^{-k_{\mathrm{el}}t} +\dfrac{k_{\mathrm{act}}P_0}{k_{\mathrm{el}}-k_{\mathrm{act}}} \left(e^{-k_{\mathrm{act}}t}-e^{-k_{\mathrm{el}}t}\right), & k_{\mathrm{el}}\ne k_{\mathrm{act}},\\[8pt] e^{-k_{\mathrm{act}}t}\left(E_0+k_{\mathrm{act}}P_0t\right), & k_{\mathrm{el}}=k_{\mathrm{act}}. \end{cases}

For unequal rates, the formation term has the difference of two exponentials. It remains nonnegative: the sign of $e^{-k_{\mathrm{act}}t}-e^{-k_{\mathrm{el}}t}$ matches the sign of $k_{\mathrm{el}}-k_{\mathrm{act}}$. The solution therefore respects nonnegative initial amounts and does not imply negative active enzyme. If $k_{\mathrm{act}}=0$, no precursor is converted and $E(t)=E_0e^{-k_{\mathrm{el}}t}$. If $k_{\mathrm{el}}=0$ and $k_{\mathrm{act}}>0$, all precursor is eventually converted and $E(t)$ tends to $E_0+P_0$. When both rates are zero, both states remain at their initial values.

Worked Example: Equal Activation and Clearance Rates

Suppose a hypothetical system begins with $P_0=80$ nmol of proenzyme and $E_0=10$ nmol of active enzyme. Let $k_{\mathrm{act}}=k_{\mathrm{el}}=0.25\ \mathrm{h}^{-1}$. Find both amounts after 4 h.

Because the rate constants are equal, use the equal-rate solution: $$ E(4)=e^{-0.25(4)}\left(10+0.25(80)(4)\right). $$

The exponent is $-1$, and the expression in parentheses is $10+80=90$ nmol. Thus $E(4)=90e^{-1}=33.11$ nmol. The precursor amount is $P(4)=80e^{-0.25(4)}=80e^{-1}=29.43$ nmol.

Initially, the combined amount was $P_0+E_0=90$ nmol; after 4 h it is $29.43+33.11=62.54$ nmol. The reduction is consistent with active-enzyme clearance. Activation transfers material from the precursor state to the active state; it does not, by itself, change their combined amount.

Activity Peaks and Mass Balance

The instantaneous change in active enzyme is the difference between its formation and clearance rates. An interior stationary point therefore has to satisfy rate balance. But this equality identifies only a candidate stationary point; it does not establish that the trajectory has an interior maximum. The initial active-enzyme amount matters.

Peak condition: At any interior time $t_*$ where $E(t)$ has a local maximum or minimum, its derivative is zero, so $k_{\mathrm{act}}P(t_*)=k_{\mathrm{el}}E(t_*)$. For a local maximum, the derivative must also change from nonnegative to nonpositive around $t_*$. A peak need not occur in the observed interval, or at any positive time.

For example, if $E_0=0$, $P_0>0$, and both rate constants are positive, then $dE/dt$ is positive initially. Also $P(t)$ tends to zero, and $E(t)$ tends to zero because active enzyme is cleared. Since the trajectory starts by rising and eventually approaches zero, it has an interior maximum. With a sufficiently large $E_0$, by contrast, clearance can exceed formation from the outset and activity can decline throughout.

A second check follows by adding the two state equations:

\frac{d}{dt}\big(P(t)+E(t)\big) =-k_{\mathrm{el}}E(t).

Thus the combined amount is nonincreasing. It is constant over an interval exactly when $k_{\mathrm{el}}E(t)=0$ throughout that interval—for example, if $k_{\mathrm{el}}=0$, or if $E(t)$ is identically zero there. Otherwise, it decreases over any interval on which active enzyme is positive and clearance is present.

Worked Example: An Interior Peak With a Depleting Pool

Consider a hypothetical system with $P_0=100$ nmol, $E_0=0$ nmol, $k_{\mathrm{act}}=0.40\ \mathrm{h}^{-1}$, and $k_{\mathrm{el}}=0.20\ \mathrm{h}^{-1}$. Find the peak time and active-enzyme amount.

The rates are unequal, so $$ E(t)=\frac{0.40(100)}{0.20-0.40} \left(e^{-0.40t}-e^{-0.20t}\right) =200\left(e^{-0.20t}-e^{-0.40t}\right)\ \mathrm{nmol}. $$

Differentiate to locate a stationary point: $$ \frac{dE}{dt}=200\left(-0.20e^{-0.20t}+0.40e^{-0.40t}\right). $$ Setting this equal to zero gives $0.40e^{-0.40t}=0.20e^{-0.20t}$, or $e^{-0.20t}=0.5$. Hence $t_*=\ln(2)/0.20=3.466$ h. At that time, $e^{-0.40t_*}=0.25$, so $E(t_*)=200(0.5-0.25)=50.00$ nmol.

The derivative is positive before this time and negative afterward, so this stationary point is a maximum. The precursor continues to deplete while active enzyme is cleared; consequently the formation rate eventually becomes too small to offset clearance.

Worked Example: A Declining Trajectory With No Interior Peak

To see why rate balance does not guarantee a peak, take the hypothetical values $P_0=1$ nmol, $E_0=100$ nmol, and $k_{\mathrm{act}}=k_{\mathrm{el}}=1\ \mathrm{h}^{-1}$. The equal-rate solution is $$ E(t)=e^{-t}(100+t)\ \mathrm{nmol}. $$

Differentiating gives $\frac{dE}{dt}=e^{-t}(1-100-t)$ nmol/h. For every $t\geq0$, $1-100-t<0$, so the active-enzyme amount is strictly decreasing and has no interior peak. At $t=1$ h, for example, $E(1)=101e^{-1}=37.16$ nmol.

The formation and clearance rates are equal only if $P(t)=E(t)$. Here that equality does not occur for $t\geq0$: $P(t)=e^{-t}$ nmol, while $E(t)=e^{-t}(100+t)$ nmol. Clearance exceeds formation throughout.

Diagnostics, Alternatives, and Limitations

The model is useful when an activation mechanism draws on a defined precursor pool and active enzyme is removed on a timescale relevant to the observation period. Its predictions are testable: precursor should decline exponentially under the assumed first-order conversion, and the active-enzyme trajectory should follow the coupled solution. In practice, measurements may be noisy or may report activity rather than enzyme amount, but the state model still specifies what pattern the observations should support.

Check initial conditions and scale. $P_0$ and $E_0$ must be nonnegative and expressed in compatible amount or concentration units. The rate constants must use the inverse of the time unit used for $t$.
Inspect the trajectory, not just a fitted peak time. A fitted interior peak should be supported by observations on both sides and satisfy formation-clearance balance. A parameter estimate that predicts an unobserved peak far outside the sampling window may be weakly identified.
Look for systematic residual patterns. Persistent underprediction of late activity can indicate that the pool is replenished, clearance is not first order, or another source of active enzyme exists. Persistent overprediction of late activity can indicate faster loss or a smaller effective precursor pool than assumed.
Consider what is actually measured. An activity assay can change because of enzyme amount, catalytic function, or assay conditions. Interpreting activity as $E(t)$ requires proportionality under fixed assay conditions.

An unlimited or effectively constant precursor approximation can be reasonable over a short interval if depletion is negligible, but it does not represent eventual exhaustion of a finite pool. Reversible allosteric activation, covered in the earlier material on allosteric activation, changes the function of enzyme already present rather than converting a finite precursor state into active enzyme. These models imply different state transitions and should not be interchanged merely because both can produce increased measured activity.

The finite-pool model is deliberately restrictive. It does not include precursor synthesis, precursor degradation, multiple activation steps, saturable conversion, or active-enzyme loss that depends on concentration. If observations require one of these features, the equations must be extended and supported by informative data; adding parameters without such information can make the mechanism difficult to distinguish from alternatives.

Common Mistakes

Assuming a peak must occur. The equation $k_{\mathrm{act}}P=k_{\mathrm{el}}E$ is necessary at an interior stationary point, but a trajectory can decline from time zero, as the high-$E_0$ example shows.
Treating activation as a permanent input. A finite pool is depleted, so its formation rate $k_{\mathrm{act}}P(t)$ declines over time rather than remaining constant.
Claiming that the total must always fall. The exact derivative is $-k_{\mathrm{el}}E(t)$. The total is constant when this product is zero throughout the interval; it decreases when clearance acts on positive active-enzyme amount.
Using the unequal-rate formula when the rates are equal. Its denominator is zero at equality. Use the limiting equal-rate expression $e^{-k_{\mathrm{act}}t}(E_0+k_{\mathrm{act}}P_0t)$.

Key takeaway: A finite precursor pool makes activation a transient source of active enzyme. The coupled equations predict both precursor depletion and active-enzyme clearance; their balance determines the shape of the activity trajectory, but does not guarantee an interior peak.

Check Your Understanding

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

If $k_{\mathrm{act}}=0$, what are $P(t)$ and $E(t)$?
For $P_0=60$ nmol, $E_0=0$, $k_{\mathrm{act}}=0.30\ \mathrm{h}^{-1}$, and $k_{\mathrm{el}}=0$, calculate $P(2)$ and $E(2)$.
What equation must hold at an interior stationary point of active-enzyme amount, and why is that equation alone not enough to establish a maximum?
For $k_{\mathrm{act}}=k_{\mathrm{el}}=0.50\ \mathrm{h}^{-1}$, $P_0=20$ nmol, and $E_0=0$, calculate $E(2)$.

Solutions. (1) No precursor is converted, so $P(t)=P_0$ and $E(t)=E_0e^{-k_{\mathrm{el}}t}$. (2) $P(2)=60e^{-0.30(2)}=60e^{-0.6}=32.93$ nmol. With no active-enzyme clearance, $E(2)=60-P(2)=27.07$ nmol. (3) The condition is $k_{\mathrm{act}}P(t_*)=k_{\mathrm{el}}E(t_*)$, obtained by setting $dE/dt=0$. It establishes a stationary point only; a maximum additionally requires the derivative to change from nonnegative to nonpositive. (4) Use the equal-rate expression: $E(2)=e^{-0.50(2)}(0+0.50(20)(2))=20e^{-1}=7.358$ nmol.