Accurate biomolecular simulations account for electronic polarization

International audienceIn this perspective, we discuss where and how accounting for electronic many-body polarization affects the accuracy of classical molecular dynamics simulations of biomolecules.While the effects of electronic polarization are highly pronounced for molecules with an opposite total charge, they are also non-negligible for interactions with overall neutral molecules. For instance, neglecting these effects in important biomolecules like amino acids and phospholipids affects the structure of proteins and membranes having a large impact on interpreting experimental data as well as building coarse grained models. With the combined advances in theory, algorithms and computational power it is currently realistic to perform simulations with explicit polarizable dipoles on systems with relevant sizes and complexity. Alternatively, the effects of electronic polarization can also be included at zero additional computational cost compared to standard fixed-charge force fields using the electronic continuum correction, as was recently demonstrated for several classes of biomolecules

A Physics-Based Intermolecular Potential for Biomolecular Simulation

The grand challenge of biophysics is to use the fundamental laws of physics to predict how biological molecules will move and interact. The atomistic HIPPO (Hydrogen-like Intermolecular Polarizable Potential) force field is meant to address this challenge. It does so by breaking down the intermolecular potential energy function of biomolecular interactions into physically meaningful components (electrostatics, polarization, dispersion, and exchangerepulsion) and using this function to drive molecular dynamics simulations. This force field is able to achieve accuracy within 1 kcal/mol for each component when compared with ab initio Symmetry Adapted Perturbation Theory calculations. HIPPO is capable of this accuracy because it introduces a model electron density on every atom in the molecular system. Since the model is built on first-principles physics, it is transferable from small model systems to bulk phase. In the first test case, the HIPPO force field for water was able to reproduce the experimental density, heat of vaporization and dielectric constant to within 1%. Importantly, HIPPO has been shown to be only 10% more computationally expensive than the widely-used AMOEBA force field, meaning that more accurate simulations of larger biological molecules are well within reach

Development and analysis of Tinker-OpenMM as a GPU-based free energy perturbation engine

The utilization of computational technologies for the lead optimization process is one of the biggest challenges in the computational chemistry field. In this dissertation, I describe the addition of GPU-based absolute and relative free energy calculation methods using polarizable force field AMOEBA to Tinker-OpenMM. I then proceed to test the capabilities of this platform by studying the binding free energy and binding structures of derivatives of the MELK inhibitor IN17. Also, I present the implementation of virial-based pressure control to the Tinker-OpenMM platform that is needed for performing isobaric simulations.Cellular and Molecular Biolog