Diffusion‐controlled reactions. I. Molecular dynamics simulation of a noncontinuum model

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.

Dynamics theory for molecular liquids based on an interaction site model

Dynamics theories for molecular liquids based on an interaction site model have been developed over the past few decades and proved to be powerful tools to investigate various dynamical phenomena.

Diffusion-mediated geminate reactions under excluded volume interactions

In this paper, influence of crowding by inert particles on the geminate reaction kinetics is theoretically investigated. Time evolution equations for the survival probability of a geminate pair are derived from the master equation taking into account the correlation among all diffusing particles, and the results are compared with those obtained by Monte Carlo simulations. In general, excluded volume interactions by the inert particles slow down the diffusive motion of reactants. However, when the initial concentration of the inert particles is uniform and high, we show that additional influence of interference between reaction and correlated diffusion accelerates the transient decay of the survival probability in the diffusion-controlled limit. We also study the escape probability for a nonuniform initial distribution of the inert particles by taking the continuous limit in space. We show that reaction yield is increased when the reaction proceeds in the presence of a positive density gradient of the inert particles which inhibits the escape of reactants. The effect can be interpreted as a cage effect.

Diffusion-limited reactions in crowded environments

Diffusion-limited reactions are usually described within the Smoluchowski theory, which neglects interparticle interactions. We propose a simple way to incorporate excluded-volume effects building on simulations of hard sphere in the presence of a sink. For large values of the sink-to-particle size ratio R(s), the measured encounter rate is in good agreement with a simple generalization of the Smoluchowski equation at high densities. Reducing R(s), the encounter rate is substantially depressed and becomes even nonmonotonic for R(s)<<1. Concurrently with the saturation of the rate, stationary density waves set in close to the sink. A mean-field analysis helps to shed light on the subtle link between such ordering and the slowing down of the encounter dynamics. Finally, we show how an infinitesimal amount of nonreacting impurities can equally slow down dramatically the reaction.

Crowding effects on association reactions at membranes

The effect of macromolecular crowding on the binding of ligands to a receptor near membranes is studied using Brownian dynamics simulations. The receptor is modeled as a reactive patch on a hard surface and the ligands and crowding agents are modeled as spheres that interact via a steep repulsive interaction potential. When a ligand collides with the patch, it reacts with probability p(rxn). The association rate constant (k(infinity)) can be decomposed into contributions from diffusion-limited (k(D)) and reaction-limited (k(R)) rates, i.e., 1/k(infinity) = 1/k(D) + 1/k(R). The simulations show that k(D) is a nonmonotonic function of the volume fraction of crowding agents for receptors of small sizes. k(R) is always an increasing function of the volume fraction of crowding agents, and the association rate constant k(infinity) determined from both contributions has a qualitatively different dependence on the macromolecular crowding for high and low values of the reaction probability p(rxn). The simulation results are used to predict the velocity of the membrane protrusion driven by actin filament elongation. Based on the simple model where the protrusive force on the membrane is generated by the intercalation of actin monomers between the membrane and actin filament ends, we predict that crowding increases the local concentration of actin monomers near the filament ends and hence accelerates the membrane protrusion.

Rate theories for biologists

Some of the rate theories that are most useful for modeling biological processes are reviewed. By delving into some of the details and subtleties in the development of the theories, the review will hopefully help the reader gain a more than superficial perspective. Examples are presented to illustrate how rate theories can be used to generate insight at the microscopic level into biomolecular behaviors. An attempt is made to clear up a number of misconceptions in the literature regarding popular rate theories, including the appearance of Planck's constant in the transition-state theory and the Smoluchowski result as an upper limit for protein-protein and protein-DNA association rate constants. Future work in combining the implementation of rate theories through computer simulations with experimental probes of rate processes, and in modeling effects of intracellular environments so that theories can be used for generating rate constants for systems biology studies is particularly exciting.

Inertial effects in diffusion-limited reactions

Diffusion-limited reactions are commonly found in biochemical processes such as enzyme catalysis, colloid and protein aggregation and binding between different macromolecules in cells. Usually, such reactions are modeled within the Smoluchowski framework by considering purely diffusive boundary problems. However, inertial effects are not always negligible in real biological or physical media on typical observation time frames. This is all the more so for non-bulk phenomena involving physical boundaries, that introduce additional time and space constraints. In this paper, we present and test a novel numerical scheme, based on event-driven Brownian dynamics, that allows us to explore a wide range of velocity relaxation times, from the purely diffusive case to the underdamped regime. We show that our algorithm perfectly reproduces the solution of the Fokker-Planck problem with absorbing boundary conditions in all the regimes considered and is thus a good tool for studying diffusion-guided reactions in complex biological environments.

Effect of macromolecular crowding on reaction rates: a computational and theoretical study

The effect of macromolecular crowding on the rates of association reactions are investigated using theory and computer simulations. Reactants and crowding agents are both hard spheres, and when two reactants collide they form product with a reaction probability, p(rxn). A value of p(rxn) < 1 crudely mimics the fact that proteins must be oriented properly for an association reaction to occur. The simulations show that the dependence of the reaction rate on the volume fraction of crowding agents varies with the reaction probability. For reaction probabilities close to unity where most of encounters between reactants lead to a reaction, the reaction rate always decreases as the volume fraction of crowding agents is increased due to the reduced diffusion coefficient of reactants. On the other hand, for very small reaction probabilities where, in most of encounters, the reaction does not occur, the reaction rate increases with the volume fraction of crowding agents--in this case, due to the increase probability of a recollision. The Smoluchowski theory refined with the radiation boundary condition and the radial distribution function at contact is in quantitative agreement with simulations for the reaction rate constant and allows the quantitative analysis of both effects separately.

Toward realistic modeling of dynamic processes in cell signaling: quantification of macromolecular crowding effects

One of the major factors distinguishing molecular processes in vivo from biochemical experiments in vitro is the effect of the environment produced by macromolecular crowding in the cell. To achieve a realistic modeling of processes in the living cell based on biochemical data, it becomes necessary, therefore, to consider such effects. We describe a protocol based on Brownian dynamics simulation to characterize and quantify the effect of various forms of crowding on diffusion and bimolecular association in a simple model of interacting hard spheres. We show that by combining the elastic collision method for hard spheres and the mean field approach for hydrodynamic interaction (HI), our simulations capture the correct dynamics of a monodisperse system. The contributions from excluded volume effect and HI to the crowding effect are thus quantified. The dependence of the results on size distribution of each component in the system is illustrated, and the approach is applied as well to the crowding effect on electrostatic-driven association in both neutral and charged environments; values for effective diffusion constants and association rates are obtained for the specific conditions. The results from our simulation approach can be used to improve the modeling of cell signaling processes without additional computational burdens.

Substrate concentration dependence of the diffusion-controlled steady-state rate constant

The Smoluchowski approach to diffusion-controlled reactions is generalized to interacting substrate particles by including the osmotic pressure and hydrodynamic interactions of the nonideal particles in the Smoluchoswki equation within a local-density approximation. By solving the strictly linearized equation for the time-independent case with absorbing boundary conditions, we present an analytic expression for the diffusion-limited steady-state rate constant for small substrate concentrations in terms of an effective second virial coefficient B2*. Comparisons to Brownian dynamics simulations excluding hydrodynamic interactions show excellent agreement up to bulk number densities of B2*rho0 < approximately = 0.4 for hard sphere and repulsive Yukawa-like interactions between the substrates. Our study provides an alternative way to determine the second virial coefficient of interacting macromolecules experimentally by measuring their steady-state rate constant in diffusion-controlled reactions at low densities.

An efficient molecular dynamics simulation method for calculating the diffusion-influenced reaction rates

We present a molecular dynamics (MD) simulation method for calculating the diffusion-influenced reaction rates in the limit of low reactant concentrations. To calculate the reaction rate coefficient, we use MD trajectories of a nonreactive equilibrium system that are initiated with a pair of reactant molecules in reactive configuration. Hence reaction systems involving complicated reactant molecules with geometrically restricted reactivities can be treated with comparable efficiency as the simple hard-sphere reaction system. Compared to the similar MD method proposed by Van Beijeren, Dong, and Bocquet [J. Chem. Phys. 114, 6265 (2001)], the present method has a couple of advantages. First, reactions involving more general sink functions can be treated. Second, more accurate results can be obtained when the reaction probability upon collision is less than unity. As an application, we investigate the effects of nondiffusive dynamics and hydrodynamic interaction of reactants on the reaction rate.

Application of Fokker-Planck-Kramers equation treatment for short-time dynamics of diffusion-controlled reaction in supercritical Lennard-Jones fluids over a wide density range

The validity of a Fokker-Planck-Kramers equation (FPKE) treatment of the rate of diffusion-controlled reaction at short times [K. Ibuki and M. Ueno, J. Chem. Phys. 119, 7054 (2003)] is tested in a supercritical Lennard-Jones fluid over a wide density range by comparing it with the Langevin dynamics and molecular dynamics simulations and other theories. The density n range studied is 0.323n(c)< or =n< or =2.58n(c) and the temperature 1.52T(c), where n(c) and T(c) are the critical density and temperature, respectively. For the rate of bimolecular reactions, the transition between the collision-limited and diffusion-limited regimes is expected to take place in this density range. The simulations show that the rate constant decays with time extensively at high densities, and that the magnitude of decay decreases gradually with decreasing density. The decay profiles of the rate constants obtained by the simulations are reproduced reasonably well by the FPKE treatment in the whole density range studied if a continuous velocity distribution is used in solving the FPKE approximately. If a discontinuous velocity distribution is used instead of the continuous one, the FPKE treatment leads to a rate constant much larger than the simulation results at medium and low densities. The rate constants calculated from the Smoluchowski-Collins-Kimball (SCK) theory based on the diffusion equation are somewhat smaller than the simulation results in medium and low densities when the intrinsic rate constant is chosen to adjust the steady state rate constant in the low density limit to that derived by the kinetic collision theory. The discrepancy is relatively small, so that the SCK theory provides a useful guideline for a qualitative discussion of the density effect on the rate constant.