An MM/3D-RISM approach for ligand binding affinities

Area of Science:

• Computational Chemistry
• Biophysics
• Biomolecular Modeling

Background:

• Molecular mechanics combined with Poisson-Boltzmann or generalized Born (MM/PBSA or MM/GBSA) are popular methods for calculating binding affinities.
• These semiempirical methods have limitations in accurately capturing specific molecular interactions and solvation structures.

Purpose of the Study:

• To introduce and validate a modified computational approach using statistical-mechanical, three-dimensional molecular theory of solvation (3D-RISM-KH) coupled with molecular mechanics/dynamics.
• To assess the performance of the 3D-RISM-KH method in predicting ligand-binding affinities compared to traditional MM/PBSA and MM/GBSA approaches.

Main Methods:

• The study employed the 3D-RISM-KH theory, a first-principles statistical mechanical solvation model, integrated with molecular mechanics or molecular dynamics.
• The method was tested on the binding affinities of seven biotin analogues to avidin in aqueous solution.

Main Results:

• The 3D-RISM-KH method provides a comprehensive molecular picture of solvation, accounting for chemical specificities like hydrogen bonding and hydrophobic interactions.
• Comparison with four generalized Born and two Poisson-Boltzmann methods revealed significant differences in calculated binding energies (up to 208 kJ/mol) and mean absolute deviations in relative affinities (10-43 kJ/mol).

Conclusions:

• The 3D-RISM-KH approach offers a more accurate and detailed prediction of ligand-binding affinities compared to conventional MM/PBSA and MM/GBSA methods.
• This advanced solvation theory provides a robust framework for understanding and predicting biomolecular interactions with high chemical specificity.