Investigating the robustness of the classical enzyme kinetic equations in small intracellular compartments

Classical descriptions of enzyme kinetics ignore the physical nature of the intracellular environment. Main implicit assumptions behind such approaches are that reactions occur in compartment volumes which are large enough so that molecular discreteness can be ignored and that molecular transport occurs via diffusion. Starting from a master equation description of enzyme reaction kinetics and assuming metabolic steady-state conditions, we derive novel mesoscopic rate equations which take into account (i) the intrinsic molecular noise due to the low copy number of molecules in intracellular compartments (ii) the physical nature of the substrate transport process, i.e. diffusion or vesicle-mediated transport. These equations replace the conventional macroscopic and deterministic equations in the context of intracellular kinetics. The latter are recovered in the limit of infinite compartment volumes. We find that deviations from the predictions of classical kinetics are pronounced (hundreds of percent in the estimate for the reaction velocity) for enzyme reactions occurring in compartments which are smaller than approximately 200nm, for the case of substrate transport to the compartment being mediated principally by vesicle or granule transport and in the presence of competitive enzyme inhibitors. This has implications for the common approach of modelling large intracellular reaction networks using ordinary differential equations and also for the calculation of the effective dosage of competitive inhibitor drugs.

💡 Research Summary

The paper challenges the conventional deterministic framework of enzyme kinetics by explicitly incorporating two often‑overlooked aspects of the intracellular environment: (i) the intrinsic stochasticity that arises when reactions take place in sub‑micron compartments containing only a few dozen molecules, and (ii) the physical nature of substrate delivery, which can be either simple diffusion or vesicle‑mediated transport. Starting from a full master‑equation description of the elementary enzymatic steps (E + S ↔ ES → E + P) and assuming a metabolic steady state, the authors perform a system‑size expansion to obtain mesoscopic rate equations that retain terms of order 1/V (V being the compartment volume). These additional terms act as noise‑correction factors that modify the classic Michaelis–Menten velocity law. In the limit V → ∞ the correction vanishes and the traditional deterministic expression is recovered, confirming the consistency of the approach.

The second major contribution is a mechanistic model of substrate influx. For pure diffusion, the influx can be treated as a continuous Poisson process, and the resulting rate equations reduce to the familiar form with a diffusion‑limited substrate concentration. By contrast, vesicle‑ or granule‑mediated transport is modeled as a bursty, jump‑process: a vesicle may carry multiple substrate molecules and deliver them to the compartment in a single stochastic event. This non‑Poissonian arrival pattern dramatically alters the probability distribution of the enzyme–substrate complex, especially when the compartment is small and the number of vesicles arriving per unit time is low. The authors incorporate this burstiness by coupling a compound Poisson process to the master equation, yielding an extra term that depends on the vesicle size, loading efficiency, and transport rate.

A particularly insightful analysis concerns competitive inhibitors. In classical kinetics, the inhibitor I competes with substrate S for the same active site, leading to a simple modification of the Michaelis constant (K_M → K_M(1 +