Abstract:Molecular dynamics simulations of biomolecules have been widely adopted in biomedical studies. As classical point-charge models continue to be used in routine biomolecular applications, there have been growing demands on developing polarizable force fields for handling more complicated biomolecular processes. Here we focus on a recently proposed polarizable Gaussian Multipole (pGM) model for biomolecular simulations. A key benefit of pGM is its screening of all short-range electrostatic interactions in a physically consistent manner, which is critical for stable charge-fitting and is needed to reproduce molecular anisotropy. Another advantage of pGM is that each atom's multipoles are represented by a single Gaussian function or its derivatives, allowing for more efficient electrostatics than other Gaussian-based models. In this study we present an efficient formulation for the pGM model defined with respect to a local frame formed with a set of covalent basis vectors. The covalent basis vectors are chosen to be along each atom's covalent bonding directions. The new local frame allows molecular flexibility during molecular simulations and facilitates an efficient formulation of analytical electrostatic forces without explicit torque computation. Subsequent numerical tests show that analytical atomic forces agree excellently with numerical finite-difference forces for the tested system. Finally, the new pGM electrostatics algorithm is interfaced with the PME implementation in Amber for molecular simulations under the periodic boundary conditions. To validate the overall pGM/PME electrostatics, we conducted an NVE simulation for a small water box of 512 water molecules. Our results show that, to achieve energy conservation in the polarizable model, it is important to ensure enough accuracy on both PME and induction iteration.

A Polarizable Valence Electron Density Based Force Field for High-Energy Interactions between Atoms and Molecules

A Simple Electron-Density Based Force Field Model for High-Energy Interactions between Atoms and Molecules

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Physically motivated, robust, ab initio force fields for CO2 and N2.

Accurate transferable polarization model derived from the monomer electron density

A fast, convenient, polarizable electrostatic model for molecular dynamics

A new polarizable force field for alkali and halide ions

Polarizable Mean-Field Model of Water for Biological Simulations with Amber and Charmm force fields

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications

Balanced polarizable Drude force field parameters for molecular anions: phosphates, sulfates, sulfamates, and oxides

Effective Inclusion of Electronic Polarization Improves the Description of Electrostatic Interactions: The prosECCo75 Biomolecular Force Field

A Polarizable Water Potential Derived from a Model Electron Density

Accounting for electronic polarization in nonpolarizable force fields

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

Development of Polarizable Water Force Fields for Phase Equilibrium Calculations

Polarizable Force Fields Based On Physical Models And Quantum Chemical Calculations

A polarizable ellipsoidal force field for halogen bonds.

Validation of Polarizable Force Field Parameters for Nucleic Acids by Inter-Molecular Interactions

Efficient Formulation of Polarizable Gaussian Multipole Electrostatics for Biomolecular Simulations