Toward the correction of effective electrostatic forces in explicit-solvent molecular dynamics simulations: restraints on solvent-generated electrostatic potential and solvent polarization

Abstract

Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations.Electronic supplementary materialThe online version of this article (doi:10.1007/s00214-014-1600-8) contains supplementary material, which is available to authorized users.

Compatibility of GROMOS-Derived Atomic Parameters for Lipopolysaccharide Membranes with the SPC/E Water Model and Alternative Long-Range Electrostatic Treatments Using Single Nonbonded Cutoff and Atom-Based Charge Schemes

“…Nevertheless, in simulation conditions where periodic solute copies show unphysical interactions, the RF correction allows for a user-defined tuning of the lengthscale of allowed interactions. 81 Previously, a systematic comparison of MD simulations of phosphatidylcholine bilayers using an extension of the GROMOS 53A6 82 and the SPC model have not found any significant difference between structural and dynamical properties derived from simulations using the RF or PME truncation. 52 This is consistent with our present findings for OM simulations, despite the difference in charge between the LPS (−8 e) and 1,2-dipalmitoyl-sn-giycero-3-phosphocholine (DPPC, zwitterionic) molecules.…”

Correcting electrostatic artifacts due to net‐charge changes in the calculation of ligand binding free energies

“…As a consequence, dominant errors arising from finite‐size effects might become negligible. Also, a hybrid method between postsimulation corrections and instantaneous correction schemes was described . Here, corrections were imposed on the forces on‐the‐fly during the MD simulation.…”

The effect of different cutoff schemes in molecular simulations of proteins

“…The effects on system properties have been described repeatedly. [6][7][8][9][10] Another way to treat electrostatic and van der Waals interactions are cutoff schemes, in which interactions are only computed up to a fixed atomic or molecular distance. Since a straight truncation leads to major artifacts, [11,12] a reaction-field contribution combined with shifting or switching functions are used to ensure that the energy approaches zero at the cutoff distance.…”