Rate constants for ion pair formation and dissociation in water

A molecular dynamics investigation of actinyl–ligand speciation in aqueous solution

Actinyl ions (AnO₂ⁿ⁺), the form in which actinides are commonly found in aqueous solution, are important species in the nuclear fuel cycle.

Nonexponential kinetics of ion pair dissociation in electrofreezing water

Temporally- or spatially-heterogeneous environments can participate in many kinetic processes, from chemical reactions and self-assembly to the forced dissociation of biomolecules. Here, we simulate the molecular dynamics of a model ion pair forced to dissociate in an explicit, aqueous solution. Triggering dissociation with an external electric field causes the surrounding water to electrofreeze and the ion pair population to decay nonexponentially. To further probe the role of the aqueous environment in the kinetics, we also simulate dissociation events under a purely mechanical force on the ion pair. In this case, regardless of whether the surrounding water is a liquid or already electrofrozen, the ion pair population decays exponentially with a well-defined rate constant that is specific to the medium and applied force. These simulation data, and the rate parameters we extract, suggest the disordered kinetics in an electrofreezing medium are a result of the comparable time scales of two concurrent processes, electrofreezing and dissociation.

Multiscale Estimation of Binding Kinetics Using Brownian Dynamics, Molecular Dynamics and Milestoning

The kinetic rate constants of binding were estimated for four biochemically relevant molecular systems by a method that uses milestoning theory to combine Brownian dynamics simulations with more detailed molecular dynamics simulations. The rate constants found using this method agreed well with experimentally and theoretically obtained values. We predicted the association rate of a small charged molecule toward both a charged and an uncharged spherical receptor and verified the estimated value with Smoluchowski theory. We also calculated the kon rate constant for superoxide dismutase with its natural substrate, O2-, in a validation of a previous experiment using similar methods but with a number of important improvements. We also calculated the kon for a new system: the N-terminal domain of Troponin C with its natural substrate Ca2+. The kon calculated for the latter two systems closely resemble experimentally obtained values. This novel multiscale approach is computationally cheaper and more parallelizable when compared to other methods of similar accuracy. We anticipate that this methodology will be useful for predicting kinetic rate constants and for understanding the process of binding between a small molecule and a protein receptor.

Reaction coordinates, one-dimensional Smoluchowski equations, and a test for dynamical self-consistency

We propose a method for identifying accurate reaction coordinates among a set of trial coordinates. The method applies to special cases where motion along the reaction coordinate follows a one-dimensional Smoluchowski equation. In these cases the reaction coordinate can predict its own short-time dynamical evolution, i.e., the dynamics projected from multiple dimensions onto the reaction coordinate depend only on the reaction coordinate itself. To test whether this property holds, we project an ensemble of short trajectory swarms onto trial coordinates and compare projections of individual swarms to projections of the ensemble of swarms. The comparison, quantified by the Kullback-Leibler divergence, is numerically performed for each isosurface of each trial coordinate. The ensemble of short dynamical trajectories is generated only once by sampling along an initial order parameter. The initial order parameter should separate the reactants and products with a free energy barrier, and distributions on isosurfaces of the initial parameter should be unimodal. The method is illustrated for three model free energy landscapes with anisotropic diffusion. Where exact coordinates can be obtained from Kramers-Langer-Berezhkovskii-Szabo theory, results from the new method agree with the exact results. We also examine characteristics of systems where the proposed method fails. We show how dynamical self-consistency is related (through the Chapman-Kolmogorov equation) to the earlier isocommittor criterion, which is based on longer paths.

Computing accurate potentials of mean force in electrolyte solutions with the generalized gradient-augmented harmonic Fourier beads method

We establish the accuracy of the novel generalized gradient-augmented harmonic Fourier beads (ggaHFB) method in computing free-energy profiles or potentials of mean force (PMFs) through comparison with two independent conventional techniques. In particular, we employ umbrella sampling with one dimensional weighted histogram analysis method (WHAM) and free molecular dynamics simulation of radial distribution functions to compute the PMF for the Na(+)-Cl(-) ion-pair separation to 16 A in 1.0M NaCl solution in water. The corresponding ggaHFB free-energy profile in six dimensional Cartesian space is in excellent agreement with the conventional benchmarks. We then explore changes in the PMF in response to lowering the NaCl concentration to physiological 0.3 and 0.1M, and dilute 0.0M concentrations. Finally, to expand the scope of the ggaHFB method, we formally develop the free-energy gradient approximation in arbitrary nonlinear coordinates. This formal development underscores the importance of the logarithmic Jacobian correction to reconstruct true PMFs from umbrella sampling simulations with either WHAM or ggaHFB techniques when nonlinear coordinate restraints are used with Cartesian propagators. The ability to employ nonlinear coordinates and high accuracy of the computed free-energy profiles further advocate the use of the ggaHFB method in studies of rare events in complex systems.

Water exchange dynamics of manganese(II), cobalt(II), and nickel(II) ions in aqueous solution

The first row transition metal ions Mn(2+), Co(2+), and Ni(2+) have been studied by classical umbrella sampling molecular dynamics simulations. The water exchange mechanisms, estimates of reaction rates, as well as structural changes during the activation process are discussed. Mn(2+) was found to react via an I(A) mechanism, whereas Co(2+) and Ni(2+) both proceed via I(D). Reaction rate constants are generally higher than those obtained by experiment but the simply constructed metal(II) ion-water potential reproduces the relative order quite well.