Nonequilibrium distribution function formalism for diffusion-influenced bimolecular reactions: Beyond the superposition approximation

Interplay of reactive interference and crowding effects in the diffusion-influenced reaction kinetics

We investigate the interplay of reactive interference and crowding effects in the irreversible diffusion-influenced bimolecular reactions of the type A+B→P+B by using the Brownian dynamics simulation method. It is known that the presence of nonreactive crowding agents retards the reaction rate when the volume fraction of the crowding agents is large enough. On the other hand, a high concentration of B is known to increase the reaction rate more than expected from the mass action law, although the B's may also act as crowders. Therefore, it would be interesting to see which effect dominates when the number density of B as well as the number density of the crowders increases. We will present an approximate theory that provides a reasonable account for the Brownian dynamics simulation results.

Multisite reversible association in membranes and solutions: From non-Markovian to Markovian kinetics

The role of diffusion on the kinetics of reversible association to a macromolecule with two inequivalent sites is studied. Previously, we found that, in the simplest possible description, it is not sufficient to just renormalize the rate constants of chemical kinetics, but one must introduce direct transitions between the bound states in the kinetic scheme. The physical reason for this is that a molecule that just dissociated from one site can directly rebind to the other rather than diffuse away into the bulk. Such a simple description is not valid in two dimensions because reactants can never diffuse away into the bulk. In this work, we consider a variety of more sophisticated implementations of our recent general theory that are valid in both two and three dimensions. We compare the predicted time dependence of the concentrations for a wide range of parameters and establish the range of validity of various levels of the general theory.

Surface reaction-diffusion kinetics on lattice at the microscopic scale

Microscopic models of reaction-diffusion processes on the cell membrane can link local spatiotemporal effects to macroscopic self-organized patterns often observed on the membrane. Simulation schemes based on the microscopic lattice method (MLM) can model these processes at the microscopic scale by tracking individual molecules, represented as hard spheres, on fine lattice voxels. Although MLM is simple to implement and is generally less computationally demanding than off-lattice approaches, its accuracy and consistency in modeling surface reactions have not been fully verified. Using the Spatiocyte scheme, we study the accuracy of MLM in diffusion-influenced surface reactions. We derive the lattice-based bimolecular association rates for two-dimensional (2D) surface-surface reaction and one-dimensional (1D) volume-surface adsorption according to the Smoluchowski-Collins-Kimball model and random walk theory. We match the time-dependent rates on lattice with off-lattice counterparts to obtain the correct expressions for MLM parameters in terms of physical constants. The expressions indicate that the voxel size needs to be at least 0.6% larger than the molecule to accurately simulate surface reactions on triangular lattice. On square lattice, the minimum voxel size should be even larger, at 5%. We also demonstrate the ability of MLM-based schemes such as Spatiocyte to simulate a reaction-diffusion model that involves all dimensions: three-dimensional (3D) diffusion in the cytoplasm, 2D diffusion on the cell membrane, and 1D cytoplasm-membrane adsorption. With the model, we examine the contribution of the 2D reaction pathway to the overall reaction rate at different reactant diffusivity, reactivity, and concentrations.

Theory of Diffusion-Influenced Reaction Networks

A formalism is developed to describe how diffusion alters the kinetics of coupled reversible association-dissociation reactions in the presence of conformational changes that can modify the reactivity. The major difficulty in constructing a general theory is that, even to the lowest order, diffusion can change the structure of the rate equations of chemical kinetics by introducing new reaction channels (i.e., modifies the kinetic scheme). Therefore, the right formalism must be found that allows the influence of diffusion to be described in a concise and elegant way for networks of arbitrary complexity. Our key result is a set of non-Markovian rate equations involving stoichiometric matrices and net reaction rates (fluxes), in which these rates are coupled by a time-dependent pair association flux matrix, whose elements have a simple physical interpretation. Specifically, each element is the probability density that an isolated pair of reactants irreversibly associates at time t via one reaction channel on the condition that it started out with the dissociation products of another (or the same) channel. In the Markovian limit, the coupling of the chemical rates is described by committors (or splitting/capture probabilities). The committor is the probability that an isolated pair of reactants formed by dissociation at one site will irreversibly associate at another site rather than diffuse apart. We illustrate the use of our formalism by considering three reversible reaction schemes: (1) binding to a single site, (2) binding to two inequivalent sites, and (3) binding to a site whose reactivity fluctuates. In the first example, we recover the results published earlier, while in the second one we show that a new reaction channel appears, which directly connects the two bound states. The third example is particularly interesting because all species become coupled and an exchange-type bimolecular reaction appears. In the Markovian limit, some of the diffusion-modified rate constants that describe new transitions become negative, indicating that memory effects cannot be ignored.

Theoretical description of spin-selective reactions of radical pairs diffusing in spherical 2D and 3D microreactors

In this work, we treat spin-selective recombination of a geminate radical pair (RP) in a spherical "microreactor," i.e., of a RP confined in a micelle, vesicle, or liposome. We consider the microreactor model proposed earlier, in which one of the radicals is located at the center of the micelle and the other one undergoes three-dimensional diffusion inside the micelle. In addition, we suggest a two-dimensional model, in which one of the radicals is located at the "pole" of the sphere, while the other one diffuses on the spherical surface. For this model, we have obtained a general analytical expression for the RP recombination yield in terms of the free Green function of two-dimensional diffusion motion. In turn, this Green function is expressed via the Legendre functions and thus takes account of diffusion over a restricted spherical surface and its curvature. The obtained expression allows one to calculate the RP recombination efficiency at an arbitrary magnetic field strength. We performed a comparison of the two models taking the same geometric parameters (i.e., the microreactor radius and the closest approach distance of the radicals), chemical reactivity, magnetic interactions in the RP and diffusion coefficient. Significant difference between the predictions of the two models is found, which is thus originating solely from the dimensionality effect: for different dimensionality of space, the statistics of diffusional contacts of radicals becomes different altering the reaction yield. We have calculated the magnetic field dependence of the RP reaction yield and chemically induced dynamic nuclear polarization of the reaction products at different sizes of the microreactor, exchange interaction, and spin relaxation rates. Interestingly, due to the intricate interplay of diffusional contacts of reactants and spin dynamics, the dependence of the reaction yield on the microreactor radius is non-monotonous. Our results are of importance for (i) interpreting experimental data for magnetic field effects on RP recombination in confined space and (ii) for describing kinetics of chemical reactions, which occur predominantly on the surfaces of biomembranes, i.e., lipid peroxidation reactions.

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.

Concentration effects on the rates of irreversible diffusion-influenced reactions

We formulate a new theory of the effects of like-particle interactions on the irreversible diffusion-influenced bimolecular reactions of the type A + B → P + B by considering the evolution equation of the triplet ABB number density field explicitly. The solution to the evolution equation is aided by a recently proposed method for solving the Fredholm integral equation of the second kind. We evaluate the theory by comparing its predictions with the results of extensive computer simulations. The present theory provides a reasonable explanation of the simulation results.

General theory of multistage geminate reactions of isolated pairs of reactants. I. Kinetic equations

General matrix approach to the consideration of multistage geminate reactions of isolated pairs of reactants depending on reactant mobility is formulated on the basis of the concept of "effective" particles. Various elementary reactions (stages of multistage reaction including physicochemical processes of internal quantum state changes) proceeding with the participation of isolated pairs of reactants (or isolated reactants) are taken into account. Investigation has been made in terms of kinetic approach implying the derivation of general (matrix) kinetic equations for local and mean probabilities of finding any of the reaction species in the sample under study (or for local and mean concentrations). The recipes for the calculation of kinetic coefficients of the equations for mean quantities in terms of relative coordinates of reactants have been formulated in the general case of inhomogeneous reacting systems. Important specific case of homogeneous reacting systems is considered.

A new multicompartmental reaction-diffusion modeling method links transient membrane attachment of E. coli MinE to E-ring formation

Many important cellular processes are regulated by reaction-diffusion (RD) of molecules that takes place both in the cytoplasm and on the membrane. To model and analyze such multicompartmental processes, we developed a lattice-based Monte Carlo method, Spatiocyte that supports RD in volume and surface compartments at single molecule resolution. Stochasticity in RD and the excluded volume effect brought by intracellular molecular crowding, both of which can significantly affect RD and thus, cellular processes, are also supported. We verified the method by comparing simulation results of diffusion, irreversible and reversible reactions with the predicted analytical and best available numerical solutions. Moreover, to directly compare the localization patterns of molecules in fluorescence microscopy images with simulation, we devised a visualization method that mimics the microphotography process by showing the trajectory of simulated molecules averaged according to the camera exposure time. In the rod-shaped bacterium Escherichia coli, the division site is suppressed at the cell poles by periodic pole-to-pole oscillations of the Min proteins (MinC, MinD and MinE) arising from carefully orchestrated RD in both cytoplasm and membrane compartments. Using Spatiocyte we could model and reproduce the in vivo MinDE localization dynamics by accounting for the previously reported properties of MinE. Our results suggest that the MinE ring, which is essential in preventing polar septation, is largely composed of MinE that is transiently attached to the membrane independently after recruited by MinD. Overall, Spatiocyte allows simulation and visualization of complex spatial and reaction-diffusion mediated cellular processes in volumes and surfaces. As we showed, it can potentially provide mechanistic insights otherwise difficult to obtain experimentally.

Rate kernel theory for pseudo-first-order kinetics of diffusion-influenced reactions and application to fluorescence quenching kinetics

Theoretical foundation of rate kernel equation approaches for diffusion-influenced chemical reactions is presented and applied to explain the kinetics of fluorescence quenching reactions. A many-body master equation is constructed by introducing stochastic terms, which characterize the rates of chemical reactions, into the many-body Smoluchowski equation. A Langevin-type of memory equation for the density fields of reactants evolving under the influence of time-independent perturbation is derived. This equation should be useful in predicting the time evolution of reactant concentrations approaching the steady state attained by the perturbation as well as the steady-state concentrations. The dynamics of fluctuation occurring in equilibrium state can be predicted by the memory equation by turning the perturbation off and consequently may be useful in obtaining the linear response to a time-dependent perturbation. It is found that unimolecular decay processes including the time-independent perturbation can be incorporated into bimolecular reaction kinetics as a Laplace transform variable. As a result, a theory for bimolecular reactions along with the unimolecular process turned off is sufficient to predict overall reaction kinetics including the effects of unimolecular reactions and perturbation. As the present formulation is applied to steady-state kinetics of fluorescence quenching reactions, the exact relation between fluorophore concentrations and the intensity of excitation light is derived.

Diffusion-influenced excited-state reversible transfer reactions, A*+B⇌C*+D, with two different lifetimes: Theories and simulations

We report accurate Brownian simulation results for the kinetics of the pseudo-first-order diffusion-influenced excited-state reversible transfer reaction A(*) + Bright harpoon over left harpoonC(*) + D with two different lifetimes using two different propagation algorithms. The results are used to test approximate solutions for this many-particle problem. Available theories fail when one of the two reactions or (decay) rate constants is large. To remedy this situation, we develop two uniform approximations, which are based on introducing a generalized Smoluchowski term into the relaxation-time approximation. The best of these is the extended unified theory of reversible target reactions, which reduces correctly in all limits and exhibits superior agreement with simulations.

Elementary Steps in Excited-State Proton Transfer

The absorption of a photon by a hydroxy-aromatic photoacid triggers a cascade of events contributing to the overall phenomenon of intermolecular excited-state proton transfer. The fundamental steps involved were studied over the last 20 years using a combination of theoretical and experimental techniques. They are surveyed in this sequel in sequential order, from fast to slow. The excitation triggers an intramolecular charge transfer to the ring system, which is more prominent for the anionic base than the acid. The charge redistribution, in turn, triggers changes in hydrogen-bond strengths that set the stage for the proton-transfer step itself. This step is strongly influenced by the solvent, resulting in unusual dependence of the dissociation rate coefficient on water content, temperature, and isotopic substitution. The photolyzed proton can diffuse in the aqueous solution in a mechanism that involves collective changes in hydrogen-bonding. On longer times, it may recombine adiabatically with the excited base or quench it. The theory for these diffusion-influenced geminate reactions has been developed, showing nice agreement with experiment. Finally, the effect of inert salts, bases, and acids on these reactions is analyzed.

The integral encounter theory of multistage reactions containing association-dissociation reaction stages. III. Taking account of quantum states of reactants

The formalism developed in Part I [K. L. Ivanov, N. N. Lukzen, A. A. Kipriyanov, and A. B. Doktorov, Phys. Chem. Chem. Phys. 6, 1706 (2004)] of the present contribution is extended to treat the reacting particles with internal quantum states. Initial spatial correlations of reactants are considered in the framework of this formalism as well.

Influence of diffusion on the kinetics of excited-state association–dissociation reactions: Comparison of theory and simulation

Several recent theories of the kinetics of diffusion influenced excited-state association--dissociation reactions are tested against accurate Brownian dynamics simulation results for a wide range of parameters. The theories include the relaxation time approximation (RTA), multiparticle kernel decoupling approximations and the so-called kinetic theory. In the irreversible limit, none of these theories reduce to the Smoluchowski result. For the pseudo-first-order target problem, we show how the RTA can be modified so that the resulting formalism does reduce correctly in the irreversible limit. We call this the unified Smoluchowski approximation, because it unites modern theories of reversible reactions with Smoluchowski's theory of irreversible reactions.

Exact dynamics of a continuous time random walker in the presence of a boundary: beyond the intuitive boundary condition approach

We derive the exact dynamics of a random walker with arbitrary non-Markovian transport and reaction rate distribution at a boundary, and present exact solutions in the continuum limit. We find that the ultimate escape probability of the particle is independent of the transport mechanism in contradiction to the long-standing belief based on the conventional approach. We also find a phase transition in the relaxation kinetics associated with the heterogeneity of the transport media.

Green's function approach to nonclassical reaction kinetics in fractal media

We present an accurate description for the mean field kinetics of the reversible recombination reaction occurring in fractal media, which is in good agreement with Monte Carlo simulation results from an intermediate to a long time region. The central dynamic quantity is the Green's function for the generalized diffusion equation that is defined in a hypothetical space whose dimensionality is given by the fracton dimension d. An exact expression for the Green's function that is valid for arbitrary values of d is presented.