Accurate Host-Guest Binding Free Energies Using the AMOEBA Polarizable Force Field
Moses K J Chung1,2, Ryan J Miller3, Borna Novak1
1Medical Scientist Training Program, Washington University School of Medicine, Saint Louis, Missouri 63110, United States.
Journal of Chemical Information and Modeling
|April 19, 2023
Summary
The Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field accurately predicts molecular binding affinities in simulations. This polarizable model shows excellent agreement with experimental data for host-guest systems.
Area of Science:
- Computational biophysics
- Molecular modeling
- Physical chemistry
Background:
- Accurate prediction of molecular interactions is a key challenge in computational biophysics.
- Molecular dynamics (MD) simulations are increasingly used for calculating binding affinities.
- The choice between fixed charge and polarizable force fields in MD is debated.
Purpose of the Study:
- To evaluate the performance of the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) polarizable force field.
- To compare AMOEBA against fixed charge models in host-guest binding challenges.
- To investigate the impact of ions on MD binding affinity calculations.
Main Methods:
- Utilized the AMOEBA polarizable atomic multipole force field.
- Performed molecular dynamics (MD) simulations for SAMPL7 and SAMPL8 Gibb octaacid host-guest challenges.
- Employed alchemical free energy methods and explored co-alchemical protocols for ion inclusion.
Main Results:
- AMOEBA demonstrated excellent agreement with experimental binding free energies (0.848 kcal/mol mean unsigned error).
- The polarizable force field better represents electrostatic potentials and host cavity water compared to fixed charge models.
- Higher salt concentrations significantly enhanced binding affinity due to charge screening effects.
Conclusions:
- The AMOEBA force field achieves chemical accuracy for realistic molecular systems, including host-guest interactions.
- Polarizable force fields offer significant advantages over fixed charge models for MD simulations.
- Ion concentration plays a crucial role in MD simulations of binding affinities, particularly for charged systems.
Related Concept Videos
The Equilibrium Binding Constant and Binding Strength
13.1K
The equilibrium binding constant (Kb) quantifies the strength of a protein-ligand interaction. Kb can be calculated as follows when the reaction is at equilibrium:
13.1K
Calculating Standard Free Energy Changes
21.6K
The free energy change for a reaction that occurs under the standard conditions of 1 bar pressure and at 298 K is called the standard free energy change. Since free energy is a state function, its value depends only on the conditions of the initial and final states of the system. A convenient and common approach to the calculation of free energy changes for physical and chemical reactions is by use of widely available compilations of standard state thermodynamic data. One method involves the...
21.6K
Noncovalent Attractions in Biomolecules
52.1K
Noncovalent attractions are associations within and between molecules that influence the shape and structural stability of complexes. These interactions differ from covalent bonding in that they do not involve sharing of electrons.
Four types of noncovalent interactions are hydrogen bonds, van der Waals forces, ionic bonds, and hydrophobic interactions.
Hydrogen bonding results from the electrostatic attraction of a hydrogen atom covalently bonded to a strong-electronegative atom like oxygen,...
Four types of noncovalent interactions are hydrogen bonds, van der Waals forces, ionic bonds, and hydrophobic interactions.
Hydrogen bonding results from the electrostatic attraction of a hydrogen atom covalently bonded to a strong-electronegative atom like oxygen,...
52.1K
Molecular Geometry and Dipole Moments
13.2K
The VSEPR theory can be used to determine the electron pair geometries and molecular structures as follows:
13.2K
Thermodynamic Potentials
897
Thermodynamic potentials are state functions that are extremely useful in analyzing a thermodynamic system. They have dimensions of energy. The four important thermodynamic potentials are internal energy, enthalpy, Helmholtz free energy, and Gibbs free energy. These thermodynamic potentials can be expressed using two of the following variables: pressure, volume, temperature, and entropy. These two variables are expressed as the rate of change of the thermodynamic potential with respect to other...
897
Gibbs Free Energy
33.7K
One of the challenges of using the second law of thermodynamics to determine if a process is spontaneous is that it requires measurements of the entropy change for the system and the entropy change for the surroundings. An alternative approach involving a new thermodynamic property defined in terms of system properties only was introduced in the late nineteenth century by American mathematician Josiah Willard Gibbs. This new property is called the Gibbs free energy (G) (or simply the free...
33.7K


