About cookies on this site Our websites require some cookies to function properly (required). In addition, other cookies may be used with your consent to analyze site usage, improve the user experience and for advertising. For more information, please review your options. By visiting our website, you agree to our processing of information as described in IBM’sprivacy statement. To provide a smooth navigation, your cookie preferences will be shared across the IBM web domains listed here.
Publication
Journal of Physical Chemistry B
Paper
Simulation of Ca2+ and Mg2+ solvation using polarizable atomic multipole potential
Abstract
The alkaline earth metals calcium and magnesium are critically involved in many biomolecular processes. To understand the hydration thermodynamics of these ions, we have performed molecular dynamics simulations using a polarizable potential. Particle-mesh Ewald for point multipoles has been applied to the calculation of electrostatic interactions. The parameters in this model have been determined from an ab initio quantum mechanical calculation of dimer interactions between ions and water. Two methods for ion solvation free energy calculation, free energy perturbation, and the Bennett acceptance ratio have been compared. Both predict results consistent with other theoretical estimations while the Bennett approach leads to a much smaller statistical error. Based on the Born theory and the ion-oxygen radial distribution functions, we estimate the effective size of the ions in solution, concluding that K+ > Na+ ≈ Ca2+ > Mg 2+. There appears to be much stronger perturbation in water structure, dynamics, and dipole moment around the divalent cations than the monovalent K+ and Na+. The average water coordination numbers for Ca2+ and Mg2+ are 7.3 and 6, respectively. The lifetime of water molecules in the first solvation shell of Mg2+ is on the order of hundreds of picoseconds, in contrast to only few picoseconds for Ca2+, K+or Na+. © 2006 American Chemical Society.