Computer simulations of alkali-acetate solutions : Accuracy of the forcefields in difference concentrations

When proteins are solvated in electrolyte solutions that contain alkali ions, the ions interact mostlywith carboxylates on the protein surface. Correctly accounting for alkali-carboxylate interactionsis thus important for realistic simulations of proteins. Acetates are the simplest carboxylates thatare amphipathic, and experimental data for alkali acetate solutions are available and can be comparedwith observables obtained from simulations. We carried out molecular dynamics simulations of alkali acetate solutions using polarizable and non-polarizable forcefields and examined the ionacetateinteractions. In particular, activity coefficients and association constants were studied in a range of concentrations (0.03, 0.1, and 1M). In addition, quantum-mechanics (QM) based energy decomposition analysis was performed in order to estimate the contribution of polarization, electrostatics, dispersion, and QM (non-classical) effects on the cation-acetate and cation-water interactions. Simulations of Li-acetate solutions in general overestimated the binding of Li+ and acetates. In lower concentrations, the activity coefficients of alkali-acetate solutions were too high, which is suggested to be due to the simulation protocol and not the forcefields. Energy decomposition analysis suggested that improvement of the forcefield parameters to enable accurate simulations of Li-acetate solution scan be achieved but may require the use of a polarizable forcefield. Importantly, simulations with some ion parameters could not reproduce the correct ion-oxygen distances, which calls for caution in thechoice of ion parameters when protein simulations are performed in electrolyte solutions