Accurate Host-Guest Binding Free Energies Using the AMOEBA Polarizable Force Field.

Abstract

A grand challenge of computational biophysics is accurate prediction of interactions between molecules. Molecular dynamics (MD) simulations have recently gained much interest as a tool to directly compute rigorous intermolecular binding affinities. The choice of a fixed point-charge or polarizable multipole force field used in MD is a topic of ongoing discussion. To compare alternative methods, we participated in the SAMPL7 and SAMPL8 Gibb octaacid host-guest challenges to assess the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) polarizable multipole force field. Advantages of AMOEBA over fixed charge models include improved representation of molecular electrostatic potentials and better description of water occupying the unligated host cavity. Prospective predictions for 26 host-guest systems exhibit a mean unsigned error vs experiment of 0.848 kcal/mol across all absolute binding free energies, demonstrating excellent agreement between computational and experimental results. In addition, we explore two topics related to the inclusion of ions in MD simulations: use of a neutral co-alchemical protocol and the effect of salt concentration on binding affinity. Use of the co-alchemical method minimally affects computed energies, but salt concentration significantly perturbs our binding results. Higher salt concentration strengthens binding through classical charge screening. In particular, added Na+ ions screen negatively charged carboxylate groups near the binding cavity, thereby diminishing repulsive coulomb interactions with negatively charged guests. Overall, the AMOEBA results demonstrate the accuracy available through a force field providing a detailed energetic description of the four octaacid hosts and 13 charged organic guests. Use of the AMOEBA polarizable atomic multipole force field in conjunction with an alchemical free energy protocol can achieve chemical accuracy in application to realistic molecular systems.