Molecular Dynamics in Physiological Solutions: Force Fields, Alkali Metal Ions, and Ionic Strength. 2010

Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
German Research School for Simulation Sciences, FZ-Juelich/RWTH Aachen University, Aachen, Germany, Pacific Northwest National Laboratory, 902 Battelle Boulevard, Richland, Washington 99352, Rush University Medical Center, 1653 W. Congress Parkway, Chicago, Illinois 60612, and SISSA, CNR-INFN-DEMOCRITOS, and Italian Institue of Technology (IIT), SISSA Unit, Trieste, Italy.

The monovalent ions Na(+) and K(+) and Cl(-) are present in any living organism. The fundamental thermodynamic properties of solutions containing such ions is given as the excess (electro-)chemical potential differences of single ions at finite ionic strength. This quantity is key for many biological processes, including ion permeation in membrane ion channels and DNA-protein interaction. It is given by a chemical contribution, related to the ion activity, and an electric contribution, related to the Galvani potential of the water/air interface. Here we investigate molecular dynamics based predictions of these quantities by using a variety of ion/water force fields commonly used in biological simulation, namely the AMBER (the newly developed), CHARMM, OPLS, Dang95 with TIP3P, and SPC/E water. Comparison with experiment is made with the corresponding values for salts, for which data are available. The calculations based on the newly developed AMBER force field with TIP3P water agrees well with experiment for both KCl and NaCl electrolytes in water solutions, as previously reported. The simulations based on the CHARMM-TIP3P and Dang95-SPC/E force fields agree well for the KCl and NaCl solutions, respectively. The other models are not as accurate. Single cations excess (electro-)chemical potential differences turn out to be similar for all the force fields considered here. In the case of KCl, the calculated electric contribution is consistent with higher level calculations. Instead, such agreement is not found with NaCl. Finally, we found that the calculated activities for single Cl(-) ions turn out to depend clearly on the type of counterion used, with all the force fields investigated. The implications of these findings for biomolecular systems are discussed.

UI MeSH Term Description Entries

Related Publications

Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
May 2013, Journal of chemical theory and computation,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
July 2009, The journal of physical chemistry. B,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
July 2019, Chemical reviews,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
December 2006, Biophysical chemistry,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
January 2016, PloS one,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
January 2019, Chemical science,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
July 1995, Environmental science & technology,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
July 2013, The journal of physical chemistry. B,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
November 2019, ACS omega,
Chao Zhang, and Simone Raugei, and Bob Eisenberg, and Paolo Carloni
August 2010, The Journal of chemical physics,
Copied contents to your clipboard!