Kinetics of reversible diffusion influenced reactions: The self-consistent relaxation time approximation

Mobile trap algorithm for zinc detection using protein sensors

We present a mobile trap algorithm to sense zinc ions using protein-based sensors such as carbonic anhydrase (CA). Zinc is an essential biometal required for mammalian cellular functions although its intracellular concentration is reported to be very low. Protein-based sensors like CA molecules are employed to sense rare species like zinc ions. In this study, the zinc ions are mobile targets, which are sought by the mobile traps in the form of sensors. Particle motions are modeled using random walk along with the first passage technique for efficient simulations. The association reaction between sensors and ions is incorporated using a probability (p1) upon an ion-sensor collision. The dissociation reaction of an ion-bound CA molecule is modeled using a second, independent probability (p2). The results of the algorithm are verified against the traditional simulation techniques (e.g., Gillespie's algorithm). This study demonstrates that individual sensor molecules can be characterized using the probability pair (p1,p2), which, in turn, is linked to the system level chemical kinetic constants, kon and koff. Further investigations of CA-Zn reaction using the mobile trap algorithm show that when the diffusivity of zinc ions approaches that of sensor molecules, the reaction data obtained using the static trap assumption differ from the reaction data obtained using the mobile trap formulation. This study also reveals similar behavior when the sensor molecule has higher dissociation constant. In both the cases, the reaction data obtained using the static trap formulation reach equilibrium at a higher number of complex molecules (ion-bound sensor molecules) compared to the reaction data from the mobile trap formulation. With practical limitations on the number sensors that can be inserted/expressed in a cell and stochastic nature of the intracellular ionic concentrations, fluorescence from the number of complex sensor molecules at equilibrium will be the measure of the intracellular ion concentration. For reliable detection of zinc ions, it is desirable that the sensors must not bind all the zinc ions tightly, but should rather bind and unbind. Thus for a given fluorescence and with association-dissociation reactions between ions and sensors, the static trap approach will underestimate the number of zinc ions present in the system.

When does the Michaelis-Menten equation hold for fluctuating enzymes?

Enzymes are dynamic entities: both their conformation and catalytic activity fluctuate over time. When such fluctuations are relatively fast, it is not surprising that the classical Michaelis-Menten (MM) relationship between the steady-state enzymatic velocity and the substrate concentration still holds. However, recent single-molecule experiments have shown that this is the case even for an enzyme whose catalytic activity fluctuates on the 10(-4)-10 s range. The purpose of this paper is to examine various scenarios in which slowly fluctuating enzymes would still obey the MM relationship. Specifically, we consider (1) the quasi-static condition (e.g., the conformational fluctuation of the enzyme-substrate complex is much slower than binding, catalysis, and the conformational fluctuations of the free enzyme), (2) the quasi-equilibrium condition (when the substrate dissociation is much faster than catalysis, irrespective of the time scales or amplitudes of conformational fluctuations), and (3) the conformational-equilibrium condition (when the dissociation and catalytic rates depend on the conformational coordinate in the same way). For each of these scenarios, the physical meaning of the apparent Michaelis constant and catalytic rate constant is provided. Finally, as an example, the theoretical analysis of a recent single-molecule enzyme assay is considered in light of the perspectives presented in this paper.

Time-dependent rate coefficients for diffusion-influenced reactions with centrosymmetric potentials

Simple closed-form expressions are presented for the time-dependent rate coefficients of diffusion-influenced reactions in the presence of spherically symmetric potentials. For diffusion-controlled contact reactions, our expression reproduces the first two terms in both the short- and long-time expansions of the rate coefficient. At intermediate times, agreement with numerical results for the Debye-Hückel potential is found to be within a few percent for a wide range of parameters. For diffusion-influenced contact reactions (described by the radiation boundary condition), the agreement is even better. When the reactivity depends on the distance between the reactants (e.g., exponentially), our analytic result is less accurate, because it reproduces the two terms in the long-time expansion only to the linear order of the reciprocal of the diffusion coefficient. Our results should prove useful in the analysis of experimental data for diffusion-influenced reactions with centrosymmetric interaction potentials.

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.

Controlled bimolecular collisions allow sub-diffraction limited microscopy of lipid vesicles

The concentration and vesicle size-controlled collisions of single molecules with target biological assemblies allow sub-diffraction limited optical images to be obtained that are not subject to the usual photobleaching problems with single molecule experiments. For example, single molecules of the probe Nile Red in aqueous solution emit a burst of fluorescence when they collide with a 50 nm hydrophobic vesicle situated on the surface in the laser focus. The bimolecular kinetics of the bursts is defined by their on- and off-time distribution functions which depend on the concentration and diffusion of the probe and the vesicle size. The mean burst frequency changes much more sharply than does the fluorescence intensity when a vesicle is raster scanned through the laser focus. This sharpness allows the spatial resolution of two objects to be improved and separations less than the diffraction limited resolution of the conventional optical microscope to be measured. The principle of this method of trajectory time distribution optical microscopy (TTDOM) could be used in a far field optical microscopic system with a resolution of several nanometers.

Molecular-dynamics simulations for nonclassical kinetics of diffusion-controlled bimolecular reactions

Molecular-dynamics simulations are presented for the diffusion-controlled bimolecular reaction A+B<==>C in two and three dimensions. The reactants and solvent molecules are modeled as spheres interacting via continuous potential-energy functions. The interaction potential between two reactants contains a deep well that results in a reaction. When the solvent concentration is low and the reactant dynamics is essentially ballistic, the system reaches equilibrium rapidly, and the reaction follows classical kinetics with exponential decay to the equilibrium. When the solvent concentration is high the particles enter the normal diffusion regime quickly and nonclassical behavior is observed, i.e., the reactant concentrations approach equilibrium as t(-d/2) where d is the dimensionality of space. When the reaction well depth is large, however, the reaction becomes irreversible within the simulation time. In this case the reactant concentrations decay as t(-d/4). Interestingly this behavior is also observed at intermediate times for reversible reactions.

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.

Coupling of protein relaxation to ligand binding and migration in myoglobin

Protein relaxation, ligand binding, and ligand migration into a hydrophobic cavity in myoglobin are unified by a bounded diffusion model which produces an accurate fit to complex ligand rebinding data over eight decades in time and a 160 K temperature range, in qualitative agreement with time-resolved x-ray crystallography. Protein relaxation operates in a cyclic manner to move the ligand away from the binding site.

Enhancement of association rates by nonspecific binding to DNA and cell membranes

A comprehensive analytic theory is developed for the kinetics of reversible association with specific sites on DNA and receptors on cell membranes in the presence of nonspecific binding to the target surfaces. Nonspecific binding is treated as a short-range attractive potential, which is more fundamental and realistic than the surface sliding model. The presence of a surface potential around the target enhances the rate of specific association and for reversible reactions leads to deviations from single exponential relaxation.

Exact solution of the excited-state geminate A*+B⇄C*+D reaction with two different lifetimes and quenching

The authors obtain, in the Laplace transform space, the exact analytic solution for the Green function and survival probabilities for the excited-state diffusion-influenced reversible geminate reaction, A*+B <==> C*+D, with two different lifetimes and in the presence of an added quenching process. This extends a previous investigation by Popov and Agmon [J. Chem. Phys. 117, 5770 (2002)] of the ground-state reaction without quenching. The long-time asymptotic behavior of the survival probabilities is obtained in the time domain. It is found to be different from the equal-lifetime case. This paper also provides a useful short-time approximation for the kinetics.

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.