Parameters of monovalent ions in the AMBER-99 forcefield: Assessment of inaccuracies and proposed improvements

The monovalent ion parameters used by the AMBER-99 forcefield are shown to exhibit physically inaccurate behavior in molecular dynamics simulations of strong 1:1 electrolytes. These errors arise from an ad hoc adaptation of Aqvist's cation parameters. The result is the rapid formation of large, unphysical clusters at concentrations that are well below solubility limits. The observed unphysical behavior poses a serious challenge for simulating ions around highly charged polymers such as nucleic acids. In this communication, we explain the source of this unphysical behavior. To facilitate the continued use of the popular AMBER parameters, we prescribe a simple fix whereby Aqvist's cations and anions are used in conjunction with the AMBER forcefield for nucleic acids. A preliminary test of this strategy suggests that the proposed fix is reasonable and is likely to be generalizable for simulating diffuse and specific ion binding to nucleic acids.