Abstract
Accurate force field parameters for ions are essential for meaningful simulation studies of proteins and nucleic acids. Currently accepted models of ions, especially for divalent ions, do not necessarily reproduce the right physiological behavior of Ca(2+) and Mg(2+) ions. Saxena and Sept (J. Chem. Theor. Comput. 2013, 9, 3538-3542) described a model, called the multisite-ion model, where instead of treating the ions as an isolated sphere, the charge was split into multiple sites with partial charge. This model provided accurate inner shell coordination of the ion with biomolecules and predicted better free energies for proteins and nucleic acids. Here, we expand and refine the multisite model to describe the behavior of divalent ions in concentrated MgCl2 and CaCl2 electrolyte solutions, eliminating the unusual ion-ion pairing and clustering of ions which occurred in the original model. We calibrate and improve the parameters of the multisite model by matching the osmotic pressure of concentrated solutions of MgCl2 to the experimental values and then use these parameters to test the behavior of CaCl2 solutions. We find that the concentrated solutions of both divalent ions exhibit the experimentally observed behavior with correct osmotic pressure, the presence of solvent separated ion pairs instead of direct ion pairs, and no aggregation of ions. The improved multisite model for (Mg(2+) and Ca(2+)) can be used in classical simulations of biomolecules at physiologically relevant salt concentrations.