Scalable improvement of SPME multipolar electrostatics in anisotropic polarizable molecular mechanics using a general short‐range penetration correction up to quadrupoles |
| |
Authors: | Christophe Narth Louis Lagardère Étienne Polack Nohad Gresh Qiantao Wang David R Bell Joshua A Rackers Jay W Ponder Pengyu Y Ren Jean‐Philip Piquemal |
| |
Institution: | 1. UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, Paris, France;2. UPMC Univ. Paris 06, Institut du Calcul et de la Simulation, Paris, France;3. UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, Paris, France;4. UPMC Univ. Paris 06, UMR 7598, Laboratoire Jacques‐Louis Lions, Paris, France;5. Chemistry and Biology Nucleo(s)tides and immunology for Therapy (CBNIT), UMR 8601 CNRS, UFR Biomédicale, Paris, France;6. Department of Biomedical Engineering, The University of Texas at Austin, Texas;7. Computational and Molecular Biophysics Program, Washington University, St. Louis, Missouri;8. Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri |
| |
Abstract: | We propose a general coupling of the Smooth Particle Mesh Ewald SPME approach for distributed multipoles to a short‐range charge penetration correction modifying the charge‐charge, charge‐dipole and charge‐quadrupole energies. Such an approach significantly improves electrostatics when compared to ab initio values and has been calibrated on Symmetry‐Adapted Perturbation Theory reference data. Various neutral molecular dimers have been tested and results on the complexes of mono‐ and divalent cations with a water ligand are also provided. Transferability of the correction is adressed in the context of the implementation of the AMOEBA and SIBFA polarizable force fields in the TINKER‐HP software. As the choices of the multipolar distribution are discussed, conclusions are drawn for the future penetration‐corrected polarizable force fields highlighting the mandatory need of non‐spurious procedures for the obtention of well balanced and physically meaningful distributed moments. Finally, scalability and parallelism of the short‐range corrected SPME approach are addressed, demonstrating that the damping function is computationally affordable and accurate for molecular dynamics simulations of complex bio‐ or bioinorganic systems in periodic boundary conditions. © 2016 Wiley Periodicals, Inc. |
| |
Keywords: | polarizable force fields electrostatics energy decomposition analysis |
|
|