Enzyme 2026-10-13
An enzyme is a biological catalyst: it accelerates a chemical reaction and is regenerated after a successful catalytic cycle. Most enzymes are proteins, but catalytic RNA molecules, called ribozymes, also exist. Many reactions needed by a biological cell would proceed too slowly under its ordinary conditions without enzymes. See Cooper, The Cell: the central role of enzymes as biological catalysts.
The biological beauty is selective speed. A biological cell contains many possible reactants; its enzymes help determine which transformations happen on a useful timescale. In metabolism, one step's reaction product can become another step's substrate, so selective catalysts support an organized flow of material.
The ribosome is a striking biological example: the machinery that builds proteins has RNA at its catalytic centre. The complete ribosome contains both RNA and proteins, but its bond-forming active site is built from RNA. Structural evidence established this distinction in Nissen et al., The structural basis of ribosome activity in peptide bond synthesis, Science 289, 920–930 (2000).
Selection occurs at an active site, the region where substrate binding and catalytic chemistry take place. Shape and chemical interactions favour some substrates over others; binding can also reshape the enzyme, an induced fit. Recognition alone is insufficient: holding a substrate extremely tightly need not make it react quickly. Catalysis requires a favourable route through the transition state, the high-barrier configuration between reactants and reaction products. Positioning reactive groups and stabilizing this configuration relative to the bound substrate can lower that barrier. Specificity is often strong, but need not mean exactly one acceptable substrate. See Cooper, The Cell: mechanisms of enzymatic catalysis.
Figure 1.
Enzymes lower reaction barriers while leaving the overall Gibbs free energy difference unchanged
. Original illustrative profiles, with the free enzyme included in both endpoint states. The horizontal coordinate orders configurations along each route; it is neither elapsed time nor a measured distance. Actual catalytic mechanisms can have several intermediates and barriers. Details.
Read the diagram numerically. Its chosen Gibbs free energies are for the reactants and for the reaction products. The uncatalysed peak is , whereas the catalysed peak is . The forward free-energy barriers are therefore and ; viewed from the reaction products, the reverse barriers are and . Both decrease by , while the endpoint difference stays . These are constructed values, not measurements of a particular enzyme.
This separates chemical kinetics from chemical equilibrium. An enzyme changes how rapidly the system approaches chemical equilibrium; it does not change the equilibrium constant for the same overall chemical reaction under the same conditions. Faster catalysis therefore supplies no extra driving energy. In a biological cell, continuing supply of reactants and removal of reaction products can nevertheless sustain a net flow: the surroundings keep that system away from chemical equilibrium. See IUBMB, Symbolism and terminology in enzyme kinetics: basic definitions.
Some enzymes require a cofactor, an additional nonprotein component such as a metal ion or a small organic molecule. Organic helpers are often called coenzymes; they can assist transfers that the protein alone cannot readily perform. The cofactor may remain bound or exchange with the surroundings, depending on the system. A coenzyme can be chemically changed in one reaction and regenerated by another: regeneration of the enzyme does not mean every helper is unchanged in that individual reaction. Some enzymes need no additional cofactor. See Alberts et al., Molecular Biology of the Cell: small molecules associated with proteins and Alberts et al., activated carriers and their regeneration.
How much faster does adding substrate make a reaction? Enzyme kinetics gives a quantitative answer when a model's assumptions fit. For an illustrative reaction obeying Michaelis-Menten kinetics, use the Michaelis-Menten equation with substrate concentration , limiting reaction rate and Michaelis constant . At , substitution gives . At , it gives ; at , it gives . Thus a hundredfold increase in substrate between the first and last examples yields only a tenfold increase in reaction rate.
Figure 2.
Illustrative Michaelis-Menten kinetics: gives
. Original normalized curve and computed examples, with a logarithmic concentration axis to separate the three values. The dashed upper line is a limiting rate approached at high substrate concentration, not reached at a finite concentration in this model. Details.
For a chosen , those fractions mean reaction rates of , and . Doubling substrate from to raises the normalized rate only from to , a relative increase of about . In this simplest model, the catalytic machinery is already nearly fully occupied. Doubling the amount of active enzyme instead doubles and the predicted reaction rates at fixed free substrate concentration, provided the other conditions are unchanged. This distinction helps decide whether additional substrate or additional catalytic capacity could improve throughput.
The interpretation has limits. The simplest Michaelis-Menten kinetics assumes a single substrate, fixed active enzyme amount, a suitable quasi-steady-state approximation, and initial measurements with negligible reverse reaction from accumulated reaction product. Substantial depletion, interacting binding sites, or inhibition by excess substrate can invalidate this model. Within it, is the substrate concentration giving half the limiting reaction rate; it is not generally the binding dissociation constant . See IUBMB, enzyme kinetics: steady-state approximation and substrate dissociation and IUBMB, Michaelis-Menten and non-Michaelis-Menten kinetics. Measuring activity also requires specified temperature, pH and necessary cofactor availability: a weak signal can reflect unsuitable conditions rather than little enzyme. See Brooks et al., Assay Guidance Manual: basics of enzymatic assays for HTS.
Michaelis constant 2026-10-14
The substrate concentration giving half the limiting reaction rate in the Michaelis-Menten equation. In the elementary mechanism , it is under the quasi-steady-state approximation, whereas the binding dissociation constant is .
The intermediate balance is
Setting this derivative approximately to zero under the quasi-steady-state approximation gives
Both dissociation and reaction with NO remove the intermediate. Its steady concentration is not zero and does not require the first reversible pair to be at equilibrium.
The quasi-steady-state approximation for the nitrate radical balances its formation against its removal:
Equivalently,
The net dinitrogen pentoxide production is . Substituting the steady-state balance gives
The reverse rate has been retained throughout and cancels algebraically. The approximation requires sufficiently rapid radical adjustment and fails if, for example, nitrogen dioxide becomes too scarce to sustain the presumed trapping rate. It does not assume that the radical concentration is zero.
Nitrogen dioxide is consumed in reactions 1 and 2 and regenerated in reaction , so
Using from the nitrate-radical quasi-steady-state approximation gives
The factor of two is essential: an initiating event uses one nitrogen dioxide molecule, and the ensuing net radical-trapping event uses a second one. The exact nitrogen balance before making this approximation is .
Adding the initiation step to the net forward trapping step cancels the nitrate-radical reaction intermediate:
The net stoichiometry agrees with under the quasi-steady-state approximation.
The remaining fraction is . Using the tabulated value rather than replacing it by the fitted value, the constant ozone concentration makes nitrogen dioxide loss a pseudo-first-order reaction:
Integrating with initial concentration and gives
This uses the nitrate-radical quasi-steady-state approximation over the interval, constant temperature and ozone concentration, and no extra nitrogen dioxide sources or sinks. The stated and units cancel directly; no factor of is needed.
For quasi-steady enzyme kinetics at low substrate saturation, the quasi-steady-state approximation is controlled here by , rather than by alone. Even with , taking makes and .
After the initial fast layer , the third equation gives
Since , the last equation gives . Starting from , this implies on bounded intervals of the slow time. The quasi-steady second equation then gives . In the first equation,
Matching to the initial substrate value therefore yields . The relation is an outer approximation: it does not replace the initial layer needed to reconcile with . Likewise the smallness of accumulated inactivated enzyme has been justified for bounded slow times, not arbitrarily long times.
For an outer solution with bounded derivatives as , the final complex equation gives . Working backwards in then gives every , including ; setting only their sum to zero would not suffice. This is the quasi-steady-state approximation, valid after the initial complex-formation layer with the relevant dimensionless rate parameters held fixed.
The recurrence is , where . Thus , . Substitution into enzyme conservation yields
Each state produces reaction product at rate , so
The source's uppercase in the index range is read as the number of binding sites. This approximation preserves the total enzyme population and the different catalytic rates of the occupied states.
The sequential-enzyme Hill limit uses the intended distinguished limit for holds the other rates fixed while , , and tends to a positive finite constant. Then for every , , whereas remains finite and nonzero. Hence intermediate complexes have negligible occupancy relative to the unbound and fully bound states.
In the quasi-steady-state approximation rate, all intermediate numerator and denominator terms disappear. Multiplying by gives
This is the Hill equation with maximal rate and effective half-saturation concentration . Here its power follows from the sequential binding mechanism and the suppression of intermediate states.
The phrase “ finite” needs this scaling interpretation for a nontrivial limiting rate. Merely bounding that product, while permitting unrelated rates to vary, does not ensure that all intermediate weights vanish relative to the terminal one. The derivation explicitly supplies the usual fixed-other-rates, positive-product assumptions. For , the rate in the box is already the exact one-state quasi-steady-state approximation rate, without this two-endpoint limiting argument.
Write , and . Since , the exact equations are
At the prescribed initial data , , rather than zero. The quasi-steady-state approximation would instead demand , giving . This mismatch signals an initial layer.
Set and expand , . At leading order , so , while and . Thus
For a more informative substrate approximation, integration of with gives
On the slow scale and , with matched leading value . Integrating gives . In the overlap , the inner complex tends to and its substrate slope matches the expansion of this outer solution. A leading composite asymptotic expansion is
It satisfies both initial values at leading order and smoothly joins the initial complex build-up to substrate depletion. The displayed inner substrate term records the small binding-induced initial slip if greater accuracy is wanted.
Applying the law of mass action to the sequential two-site enzyme reaction network gives
Adding the appropriate equations gives the conservation laws and . In particular after nondimensionalization.
Define the dimensionless reaction rates , , , and . Substitution, with primes denoting , yields
These right sides are , respectively, with and . When , the form of these ordinary differential equations displays the faster relaxation of the complexes, which motivates a subsequent quasi-steady-state approximation.
On times , the quasi-steady-state approximation sets
The second equation gives
In the first equation the terms then cancel, leaving
Solving these two equations gives
Since ,
At small substrate concentration this allosteric rate is quadratic, , so its graph leaves the origin with zero slope. At large it saturates at . A Michaelis-Menten equation also saturates, but its rate is linear near the origin. The allosteric curve is therefore more sigmoidal:
This is the quasi-steady rate law for a two-substrate allosteric enzyme.
The quasi-steady-state approximation for sets . Substituting into the equation gives
Its equilibrium concentration is therefore
For constant input this is also the exact equilibrium of the full two-variable model: the fast- assumption is needed to replace the transient dynamics, not to calculate a fixed point.
Relative to the direct-activation result, . The extra activation gate suppresses the response at low input and can sharpen its rise. The high-input limit is , which need not reach because saturates at .
For the parameters in part (d), the following Python continuation plots the expression and changes each Hill constant separately:
def z_equilibrium(X, Kxy=1, Kxz=2, Kyz=1):
    Yeq = 2 * hill(X, Kxy, 2)
    return 3 * hill(X, Kxz, 4) * hill(Yeq, Kyz, 4)

X_grid = np.linspace(0, 5, 501)
plt.figure(figsize=(7, 4))
plt.plot(X_grid, 3*hill(X_grid, 2, 4), ":", label="Direct activation")
for parameters, label in [((1, 2, 1), "Baseline"),
                          ((2, 2, 1), "Kxy=2"),
                          ((1, 3, 1), "Kxz=3"),
                          ((1, 2, 2), "Kyz=2")]:
    plt.plot(X_grid, z_equilibrium(X_grid, *parameters), label=label)
plt.xlabel("X concentration")
plt.ylabel("Equilibrium Z")
plt.legend()
plt.tight_layout()
plt.show()
Increasing lowers at fixed finite , delaying opening of the indirect gate. Increasing delays opening of the direct gate. Neither changes the limiting response as , provided the constants remain finite. Increasing makes a given less effective and also lowers the high-input plateau. These changes are visible in the preceding figure. At small , the response is proportional to , while both gates are far below saturation, explaining the increased low-input suppression.
For a persistence detector, retain finite dynamics and use, for example,
with unit loss rates, , , and initially . On switching the input on, and
The time to reach the half-activation threshold of the indirect gate is
An isolated input pulse of duration ends well before this threshold; a duration-four pulse persists long enough to open the gate. The large but finite Hill coefficient makes this a graded filter, not a perfectly sharp cutoff. Once vanishes, its direct gate closes immediately and stops production of , which then decays with time constant one. The protein output itself does not drop discontinuously. This delayed-on, prompt-production-off mechanism is the characteristic persistence detector interpretation of an AND-type coherent feed-forward loop.[ref]
This final Python continuation simulates both pulse durations. Each discontinuity is placed on the time grid and the input is held constant within each integration step:
def persistence_rhs(state, X):
    Y, Z = state
    return np.array([hill(X, 1, 4) - Y,
                     hill(X, 1, 4)*hill(Y, 0.7, 8) - Z])

def pulse_trajectory(duration, dt=0.005):
    times = np.arange(round(8/dt) + 1) * dt
    out = np.zeros((len(times), 2))
    for j in range(len(times)-1):
        X = 2.0 if 1 <= times[j] + dt/2 < 1 + duration else 0.0
        # Each step is integrated with its constant input.
        out[j+1] = integrate(lambda s: persistence_rhs(s, X),
                             out[j], [0, dt])[-1]
    return times, out

fig, axes = plt.subplots(1, 2, figsize=(10, 4))
for duration in [0.25, 4]:
    times, out = pulse_trajectory(duration)
    axes[0].plot(times, out[:, 0], label=f"Y, duration {duration}")
    axes[1].plot(times, out[:, 1], label=f"Z, duration {duration}")
    print(duration, "maximum output", out[:, 1].max())
axes[0].axhline(0.7, color="black", linestyle=":")
for ax in axes:
    ax.set_xlabel("Time")
    ax.set_ylabel("Concentration")
    ax.legend()
plt.tight_layout()
plt.show()
The chosen example produces peak of about for the short pulse and for the long pulse. These are original illustrative simulations, not fitted experimental measurements.
Figure 1.
A coherent feed-forward loop distinguishes short and persistent input pulses
. Both inputs switch from zero to X=2 at time one. One ends after 0.25 time units and the other after four. The intermediate Y crosses its activation threshold only for the sustained input, producing a much larger target response. Details.
Finally, the strict instantaneous- reduction used for the equilibrium calculation has no separate intermediate accumulation delay. Its remaining one-variable dynamics can still attenuate short pulses through ordinary output relaxation, but it cannot demonstrate the additional delay caused by finite buildup. The pulse test therefore uses the full two-variable model.
Persistence detector Created 2026-10-20 Updated 2026-10-25
A persistence detector responds strongly to a sustained input and weakly to short isolated inputs. In an activating coherent feed-forward loop, an intermediate protein accumulates toward a threshold while the target also requires the direct input. Eliminating that intermediate by an instantaneous quasi-steady-state approximation removes its accumulation delay.
A quasi-steady-state approximation for enzyme-substrate complexes can remain useful when total enzyme and initial substrate concentrations are comparable. The relevant parameter can involve the sum , rather than alone. In a mechanism-based inhibition model with , bounded intermediate decay rates and an inactivation rate of order , the fast complexes adjust in a layer of scaled width . On bounded slow-time intervals the main complex follows the substrate, the later complex and accumulated inactivated enzyme remain , and the substrate has leading exponential depletion. An outer algebraic complex concentration does not satisfy the unbound initial condition inside that initial layer.
Reversible enzyme inhibition Created 2026-10-18 Updated 2026-10-20
Construct a minimal model with one substrate and no appreciable accumulated reaction product. Only the enzyme-substrate complex produces product, with catalytic rate constant ; and are inactive. Write and for free molar concentrations, and for total enzyme molar concentration. Assume a suitable quasi-steady-state approximation, rapid inhibitor binding and release relative to turnover, and negligible depletion of free inhibitor. In this construction the population relations are
Here is the uninhibited Michaelis constant, and are the inhibition constants for the two indicated binding branches. The symbol denotes a free concentration, not an amount added to the tube. Define . Conservation of enzyme gives an original short derivation:
Thus fitting a Michaelis-Menten equation at a fixed would return the apparent parameters
“Apparent” means conditional on that inhibitor concentration. The uninhibited Michaelis constant has not been redefined. Dimensional analysis checks the denominator: both and have units of molar concentration; the resulting has units of reaction rate.
Figure 1. . Original calculations, not measured rates or fits to trypsin. Every panel uses , and , with free . Competitive inhibition uses ; pure noncompetitive inhibition uses ; the unequal-constant mixed inhibition example uses . Dots mark . Dotted horizontal lines are limiting rates, approached rather than reached at finite . Details.
The shared numerical choices make the mechanisms comparable. At and , all three models give . At , the respective reaction rates are , , and ; at they are , , and . These are independently chosen teaching models, not measured inhibition constants for a named enzyme.
For sequential enzyme binding, a quasi-steady-state approximation gives occupied-state weights , , relative to free enzyme. If all intermediate weights tend to zero while the terminal weight remains nonzero, enzyme conservation leaves only unbound and fully bound states, and the production rate becomes the displayed Hill equation. One realization holds other rates fixed while the first binding rate tends to zero and the last tends to infinity with positive finite product. A statement merely bounding that product does not by itself specify this nontrivial distinguished limit.
Sequential two-site enzyme reaction network Created 2026-10-14 Updated 2026-10-20
An enzyme may bind one molecule of a substrate to form , then a second to form , with product-release reactions and . The law of mass action gives conservation laws and . Scaling substrate by , complexes by and time by exposes the small ratio : the substrate evolves on the slow scale while the complex equations have the form . This is a starting point for quasi-steady-state approximation when .
Steady state (chemistry) 2026-10-13
A steady state has time-independent concentrations because formation and removal rates balance for the species considered. Nonzero reaction fluxes can continue, so a steady state need not be chemical equilibrium. A quasi-steady-state approximation imposes this balance approximately for a fast intermediate while slower species evolve.