Molecular Environment-Specific Atomic Charges Improve Binding Affinity Predictions of SAMPL5 Host-Guest Systems
Duván González1, Luis Macaya1, Esteban Vöhringer-Martinez1
1Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, 4070386 Concepción, Chile.
Journal of Chemical Information and Modeling
|August 31, 2021
Summary
This study introduces a new method for calculating atomic charges in host-guest systems using quantum mechanics/molecular mechanics. This approach improves the prediction of binding affinity, especially for charged molecules.
Area of Science:
- Computational Chemistry
- Molecular Modeling
- Biophysics
Background:
- Host-guest systems are crucial for benchmarking computational methods in predicting binding free energies.
- Accurate force fields are essential for reliable binding affinity predictions in alchemical free energy calculations.
- Current methods for deriving atomic charges can limit the accuracy of these predictions.
Purpose of the Study:
- To develop a novel methodology for deriving atomic charges in host-guest systems.
- To improve the accuracy of absolute binding free energy predictions.
- To assess the impact of quantum mechanics/molecular mechanics-derived charges on host-guest binding affinity.
Main Methods:
- Utilized quantum mechanics/molecular mechanics (QM/MM) calculations with minimal basis iterative stockholder (MBIS) partitioning.
- Developed an interface between OpenMM and ORCA software for D-MBIS charge derivation.
- Calculated the energetic cost of guest polarization in both bound and unbound states.
Main Results:
- D-MBIS charges, when used with the general Amber force field, enhanced binding affinity predictions for six guests in octa acid hosts.
- The new method outperformed the AM1-BCC charge set, particularly after correcting for polarization energy differences.
- The polarization cost correction significantly impacted the binding affinity predictions for anionic guests.
Conclusions:
- The D-MBIS charge derivation methodology offers improved accuracy for host-guest binding free energy calculations.
- Accounting for the energetic cost of guest polarization is critical for accurate predictions, especially for charged species.
- This work provides a more robust computational approach for studying host-guest interactions and molecular recognition.
More Related Videos
Related Concept Videos
The Equilibrium Binding Constant and Binding Strength
14.3K
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:
14.3K
Noncovalent Attractions in Biomolecules
59.6K
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,...
59.6K
Conserved Binding Sites
4.7K
Many proteins’ biological role depends on their interactions with their ligands, small molecules that bind to specific locations on the protein known as ligand-binding sites. Ligand-binding sites are often conserved among homologous proteins as these sites are critical for protein function.
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
Binding sites are often located in large pockets, and if their location on a protein’s surface is unknown, it can be predicted using various approaches. The energetic method computationally...
4.7K
Ligand Binding Sites
14.2K
Proteins are dynamic macromolecules that carry out a wide variety of essential processes; however, the activities of most proteins depend on their interactions with other molecules or ions, known as ligands.
Protein-ligand interactions are quite specific; even though numerous potential ligands surround a cellular protein at any given time, only a particular ligand can bind to that protein. Moreover, a ligand binds only to a dedicated area on the surface of the protein, known as the...
Protein-ligand interactions are quite specific; even though numerous potential ligands surround a cellular protein at any given time, only a particular ligand can bind to that protein. Moreover, a ligand binds only to a dedicated area on the surface of the protein, known as the...
14.2K
Crystal Field Theory - Octahedral Complexes
28.5K
Crystal Field Theory
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
28.5K
Molecular Geometry and Dipole Moments
15.8K
The VSEPR theory can be used to determine the electron pair geometries and molecular structures as follows:
15.8K


