Overbinding and Qualitative and Quantitative Changes Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Force Fields with and without NBFIX and ECC Corrections

Overbinding of ions is a common and well-known problem in classical molecular dynamics simulations. One of its main causes is the absence of electronic polarizability in the force fields. The current approaches for minimizing overbinding typically either retain the original charges and use an ad hoc readjustment of the Lennard-Jones parameters as done in the nonbonded fix (NBFIX) approach or rescale the charges using a theoretical framework. The goal in the latter is to include shielding produced by the missing electronic polarizability as done in the electronic continuum correction (ECC) approach. NBFIX and ECC are the most common corrections, and we compare their performance to the default parameterizations provided by five different commonly used biomolecular force fields, OPLS-AA/L, CHARMM27, CHARMM36m, CHARMM22*, and AMBER99SB-ILDN. As test systems, we use poly-alpha,L-glutamic and poly-alpha,L-aspartic amino acid molecules in explicit water together with Na+ and K+ counterions. We demonstrate that the different force fields yield results that are not only quantitatively but also qualitatively different. The resulting structures of the macroions depend strongly on the model for ions. NBFIX corrections alleviate the problem of overbinding, resulting in extended peptides. The ECC corrections depend nontrivially on the original underlying model, and despite being based on a theoretical framework, they cannot always solve the problem.