Premium
Three‐site and five‐site fixed‐charge water models compatible with AMOEBA force field
Author(s) -
Pan Cong,
Liu Chengwen,
Peng Junhui,
Ren Pengyu,
Huang Xuhui
Publication year - 2020
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.26151
Subject(s) - water model , force field (fiction) , multipole expansion , molecular dynamics , enthalpy of vaporization , compressibility , thermodynamics , chemistry , heat capacity , dielectric , enthalpy , chemical physics , materials science , computational chemistry , physics , quantum mechanics , optoelectronics
In a typical biomolecular simulation using Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field, the vast majority molecules in the simulation box consist of water, and these water molecules consume the most CPU power due to the explicit mutual induction effect. To improve the computational efficiency, we here develop two new nonpolarizable water models (with flexible bonds and fixed charges) that are compatible with AMOEBA solute: the 3‐site AW3C and 5‐site AW5C. To derive the force‐field parameters for AW3C and AW5C, we fit to six experimental liquid thermodynamic properties: liquid density, enthalpy of vaporization, dielectric constant, isobaric heat capacity, isothermal compressibility and thermal expansion coefficient, at a broad range of temperatures from 261.15 to 353.15 K under 1.0 atm pressure. We further validate our AW3C and AW5C water models by showing that they can well reproduce the radial distribution function g ( r ), self‐diffusion constant D , and hydration free energy from the AMOEBA03 water model and the experimental observations. Furthermore, we show that our AW3C and AW5C water models can greatly accelerate (>5 times) the bulk water as well as biomolecular simulations when compared to AMOEBA water. Specifically, we demonstrate that the applications of AW3C and AW5C water models to simulate a DNA duplex lead to a threefold acceleration, and in the meanwhile well maintain the structural properties as the fully polarizable AMOEBA water. We expect that our AW3C and AW5C water models hold great promise to be widely applied to simulate complex bio‐molecules using the AMOEBA force field.