Rigorous derivation of the long-time asymptotics for reversible binding

Diffusion-induced competitive two-site binding

The influence of diffusion on the kinetics of ligand binding to a macromolecule with two sites is considered for a simple model where, in the reaction-controlled limit, there is no cooperativity and hence the sites are independent. By applying our recently developed formalism to describe a network of coupled diffusion-influenced reactions, we show that the rate constants of chemical kinetics cannot just be renormalized. Rather a new reaction channel, which connects the two singly occupied states, must be introduced. The rate constants of this new channel depend on the committor or capture probability that a ligand that just dissociated from one site rebinds to the other. This result is rederived in an elementary way using the encounter complex model. Illustrative calculations are presented where the kinetics of the fractional saturation of one site is compared with that of a macromolecule that has only this site. If all sites are initially empty, then the second site slows down binding to the first due to competition between the sites. On the other hand, if the second site is initially occupied, the binding of the first site speeds up because of the direct diffusion-induced transitions between the two singly bound states.

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.

Dynamics of different steps of the photopyrolytic cycle of an eminent anticancer drug topotecan inside biocompatible lyotropic liquid crystalline systems

Dynamics of different steps of photopyrolytic processes of an eminent anticancer drug topotecan have been investigated inside different lyotropic liquid crystalline systems.

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.

Aging processes in reversible reaction-diffusion systems

Reversible reaction-diffusion systems display anomalous dynamics characterized by a power-law relaxation toward stationarity. In this paper we study in the aging regime the nonequilibrium dynamical properties of some model systems with reversible reactions. Starting from the exact Langevin equations describing these models, we derive expressions for two-time correlation and autoresponse functions and obtain a simple aging behavior for these quantities. The autoresponse function is thereby found to depend on the specific nature of the chosen perturbation of the system.

Self-consistent theory of reversible ligand binding to a spherical cell

Ligand binding to receptors is the initial event in many signaling processes, and a quantitative understanding of this interaction is important for modeling cell behavior. In this paper, we study the kinetics of reversible ligand binding to receptors on a spherical cell surface using a self-consistent stochastic theory. Binding, dissociation, diffusion and rebinding of ligands are incorporated into the theory in a systematic manner. We derive explicitly the time evolution of the ligand-bound receptor fraction p(t) in various regimes. Contrary to the commonly accepted view, we find that the well-known Berg-Purcell scaling for the association rate is modified as a function of time. Specifically, the effective on-rate changes non-monotonically as a function of time and equals the intrinsic rate at very early as well as late times, while being approximately equal to the Berg-Purcell value at intermediate times. The effective dissociation rate, as it appears in the binding curve or measured in a dissociation experiment, is strongly modified by rebinding events and assumes the Berg-Purcell value except at very late times, where the decay is algebraic and not exponential. In equilibrium, the ligand concentration everywhere in the solution is the same and equals its spatial mean, thus ensuring that there is no depletion in the vicinity of the cell. Implications of our results for binding experiments and numerical simulations of ligand-receptor systems are also discussed.

Corrections to the law of mass action and properties of the asymptotic t=∞ state for reversible diffusion-limited reactions

For diffusion-limited reversible A+A<==>B reactions we reexamine two fundamental concepts of classical chemical kinetics-the notion of "chemical equilibrium" and the "law of mass action." We consider a general model with distance-dependent reaction rates, such that any pair of A particles, performing standard random walks on sites of a d-dimensional lattice and being at a distance mu apart of each other at time moment t, may associate forming a B particle at the rate k+(mu). In turn, any randomly moving B particle may spontaneously dissociate at the rate k-(lambda) into a geminate pair of As "born" at a distance lambda apart of each other. Within a formally exact approach based on Gardiner's Poisson representation method we show that the asymptotic t=infinity state attained by such diffusion-limited reactions is generally not a true thermodynamic equilibrium, but rather a nonequilibrium steady state, and that the law of mass action is invalid. The classical concepts hold only in case when the ratio k+(mu)k-(mu) does not depend on mu for any mu.

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.

Power-law relaxation in a complex system: Omori law after a financial market crash

We study the relaxation dynamics of a financial market just after the occurrence of a crash by investigating the number of times the absolute value of an index return is exceeding a given threshold value. We show that the empirical observation of a power law evolution of the number of events exceeding the selected threshold (a behavior known as the Omori law in geophysics) is consistent with the simultaneous occurrence of (i) a return probability density function characterized by a power law asymptotic behavior and (ii) a power-law relaxation decay of its typical scale. Our empirical observation cannot be explained within the framework of simple and widespread stochastic volatility models.

Long-time tails in the kinetics of reversible bimolecular reactions

A new approach is developed to study the relaxation of concentrations to equilibrium in reversible bimolecular reactions. For A+B <==> C, a general analytic expression is derived for the amplitude of the power law (t(-d/2)) asymptotics for arbitrary diffusion coefficients and concentrations of the reactants. Our formalism is based on the analysis of the time correlation functions describing the equilibrium fluctuations of the concentrations. This powerful and simple procedure can be readily used to study other bimolecular reactions such as A+B <==> C+D.