Most cited article - PubMed ID 22630587
Accurate description of aqueous carbonate ions: an effective polarization model verified by neutron scattering
prosECCo75 is an optimized force field effectively incorporating electronic polarization via charge scaling. It aims to enhance the accuracy of nominally nonpolarizable molecular dynamics simulations for interactions in biologically relevant systems involving water, ions, proteins, lipids, and saccharides. Recognizing the inherent limitations of nonpolarizable force fields in precisely modeling electrostatic interactions essential for various biological processes, we mitigate these shortcomings by accounting for electronic polarizability in a physically rigorous mean-field way that does not add to computational costs. With this scaling of (both integer and partial) charges within the CHARMM36 framework, prosECCo75 addresses overbinding artifacts. This improves agreement with experimental ion binding data across a broad spectrum of systems─lipid membranes, proteins (including peptides and amino acids), and saccharides─without compromising their biomolecular structures. prosECCo75 thus emerges as a computationally efficient tool providing enhanced accuracy and broader applicability in simulating the complex interplay of interactions between ions and biomolecules, pivotal for improving our understanding of many biological processes.
The inclusion of electronic polarization is of crucial importance in molecular simulations of systems containing charged moieties. When neglected, as often done in force field simulations, charge-charge interactions in solution may become severely overestimated, leading to unrealistically strong bindings of ions to biomolecules. The electronic continuum correction introduces electronic polarization in a mean-field way via scaling of charges by the reciprocal of the square root of the high-frequency dielectric constant of the solvent environment. Here, we use ab initio molecular dynamics simulations to quantify the effect of electronic polarization on pairs of like-charged ions in a model nonaqueous environment where electronic polarization is the only dielectric response. Our findings confirm the conceptual validity of this approach, underlining its applicability to complex aqueous biomolecular systems. Simultaneously, the results presented here justify the potential employment of weaker charge scaling factors in force field development.
- Publication type
- Journal Article MeSH
We employed density functional theory-based ab initio molecular dynamics simulations to examine the hydration structure of several common alkali and alkali earth metal cations. We found that the commonly used atom pairwise dispersion correction scheme D3, which assigns dispersion coefficients based on the neutral form of the atom rather than its actual oxidation state, leads to inaccuracies in the hydration structures of these cations. We evaluated this effect for lithium, sodium, potassium, and calcium and found that the inaccuracies are particularly pronounced for sodium and potassium compared to the experiment. To remedy this issue, we propose disabling the D3 correction specifically for all cation-including pairs, which leads to a much better agreement with experimental data.
- Publication type
- Journal Article MeSH
Models of the hydrogenoxalate (bioxalate, charge -1) and oxalate (charge -2) anions were developed for classical molecular dynamics (CMD) simulations and parametrized against ab initio molecular dynamics (AIMD) data from our previous study (Kroutil et al. (2016) J Mol Model 22:210). The interactions of the anions with water were described using charges scaled according to the electronic continuum correction approach with rescaling of nonbonded parameters (ECCR), and those descriptions of anion interactions were found to agree well with relevant AIMD and experimental results. The models with full RESP charges showed excessively strong electrostatic interactions between the solute and water molecules, leading to an overstructured solvation shell around the anions and thus to a diffusion coefficient that was much too low. The effect of charge scaling was more evident for the oxalate dianion than for the hydrogenoxalate anion. Our work provides CMD models for ions of oxalic acid and extends previous studies that showed the importance of ECCR for modeling divalent ions and ions of organic compounds. Graphical abstract The radial distribution function between a water oxygen (Ow) and an oxygen of the oxalate dianion (Ox) significantly improved when scaled charges were applied to the anion.
- Keywords
- AIMD, Ab initio molecular dynamics, CMD, Classical molecular dynamics, ECCR, Electronic continuum correction, Hydrogenoxalate, Oxalate, Oxalic acid anions,
- Publication type
- Journal Article MeSH
