Diffusion influence on Michaelis–Menten kinetics

Effective reaction rates for diffusion-limited reaction cycles

Biological signals in cells are transmitted with the use of reaction cycles, such as the phosphorylation-dephosphorylation cycle, in which substrate is modified by antagonistic enzymes. An appreciable share of such reactions takes place in crowded environments of two-dimensional structures, such as plasma membrane or intracellular membranes, and is expected to be diffusion-controlled. In this work, starting from the microscopic bimolecular reaction rate constants and using estimates of the mean first-passage time for an enzyme-substrate encounter, we derive diffusion-dependent effective macroscopic reaction rate coefficients (EMRRC) for a generic reaction cycle. Each EMRRC was found to be half of the harmonic average of the microscopic rate constant (phosphorylation c or dephosphorylation d), and the effective (crowding-dependent) motility divided by a slowly decreasing logarithmic function of the sum of the enzyme concentrations. This implies that when c and d differ, the two EMRRCs scale differently with the motility, rendering the steady-state fraction of phosphorylated substrate molecules diffusion-dependent. Analytical predictions are verified using kinetic Monte Carlo simulations on the two-dimensional triangular lattice at the single-molecule resolution. It is demonstrated that the proposed formulas estimate the steady-state concentrations and effective reaction rates for different sets of microscopic reaction rates and concentrations of reactants, including a non-trivial example where with increasing diffusivity the fraction of phosphorylated substrate molecules changes from 10% to 90%.

Reversible Stochastically Gated Diffusion-Influenced Reactions

An approximate but accurate theory is developed for the kinetics of reversible binding of a ligand to a macromolecule when either can stochastically fluctuate between reactive and unreactive conformations. The theory is based on a set of reaction-diffusion equations for the deviations of the pair distributions from their bulk values. The concentrations are shown to satisfy non-Markovian rate equations with memory kernels that are obtained by solving an irreversible geminate (i.e., two-particle) problem. The relaxation to equilibrium is not exponential but rather a power law. In the Markovian limit, the theory reduces to a set of ordinary rate equations with renormalized rate constants.

Diffusional correlations among multiple active sites in a single enzyme

Simulations of the enzymatic dynamics of a model enzyme containing multiple substrate binding sites indicate the existence of diffusional correlations in the chemical reactivity of the active sites. A coarse-grain, particle-based, mesoscopic description of the system, comprising the enzyme, the substrate, the product and solvent, is constructed to study these effects. The reactive and non-reactive dynamics is followed using a hybrid scheme that combines molecular dynamics for the enzyme, substrate and product molecules with multiparticle collision dynamics for the solvent. It is found that the reactivity of an individual active site in the multiple-active-site enzyme is reduced substantially, and this effect is analyzed and attributed to diffusive competition for the substrate among the different active sites in the enzyme.

Effective reaction rates in diffusion-limited phosphorylation-dephosphorylation cycles

We investigate the kinetics of the ubiquitous phosphorylation-dephosphorylation cycle on biological membranes by means of kinetic Monte Carlo simulations on the triangular lattice. We establish the dependence of effective macroscopic reaction rate coefficients as well as the steady-state phosphorylated substrate fraction on the diffusion coefficient and concentrations of opposing enzymes: kinases and phosphatases. In the limits of zero and infinite diffusion, the numerical results agree with analytical predictions; these two limits give the lower and the upper bound for the macroscopic rate coefficients, respectively. In the zero-diffusion limit, which is important in the analysis of dense systems, phosphorylation and dephosphorylation reactions can convert only these substrates which remain in contact with opposing enzymes. In the most studied regime of nonzero but small diffusion, a contribution linearly proportional to the diffusion coefficient appears in the reaction rate. In this regime, the presence of opposing enzymes creates inhomogeneities in the (de)phosphorylated substrate distributions: The spatial correlation function shows that enzymes are surrounded by clouds of converted substrates. This effect becomes important at low enzyme concentrations, substantially lowering effective reaction rates. Effective reaction rates decrease with decreasing diffusion and this dependence is more pronounced for the less-abundant enzyme. Consequently, the steady-state fraction of phosphorylated substrates can increase or decrease with diffusion, depending on relative concentrations of both enzymes. Additionally, steady states are controlled by molecular crowders which, mostly by lowering the effective diffusion of reactants, favor the more abundant enzyme.

Role of diffusion in the kinetics of reversible enzyme-catalyzed reactions

The accurate expression for the steady-state velocity of an irreversible enzyme-catalyzed reaction obtained by Shin and co-workers is generalized to allow for the rebinding of the product. The amplitude of the power-law (t(-1/2)) relaxation of the free- and bound-enzyme concentrations to steady-state values is expressed in terms of the steady-state velocity and the intrinsic (chemical) rate constants. This result is conjectured to be exact, even though our expression for the steady-state velocity in terms of microscopic parameters is only approximate.