Abstract:The accurate prediction of binding free energy is a major challenge in structure-based drug design. Quantum mechanics (QM)-based approaches show promising potential in predicting ligand–protein binding affinity by accurately describing the behavior and structure of electrons. However, traditional QM calculations face computational limitations, hindering their practical application in drug design. Nevertheless, the fragment molecular orbital (FMO) method has gained widespread application in drug design due to its ability to reduce computational costs and achieve efficient ab initio QM calculations. Although the FMO method has demonstrated its reliability in calculating the gas phase potential energy, the binding of proteins and ligands also involves other contributing energy terms, such as solvent effects, the 'deformation energy' of a ligand's bioactive conformations, and entropy. Particularly in cases involving ionized fragments, the calculation of solvation free energy becomes particularly crucial. We conducted an evaluation of some previously reported implicit solvent methods on the same data set to assess their potential for improving the performance of the FMO method. Herein, we develop a new QM-based binding free energy calculation method called FMOScore, which enhances the performance of the FMO method. The FMOScore method incorporates linear fitting of various terms, including gas-phase potential energy, deformation energy, and solvation free energy. Compared to other widely used traditional prediction methods such as FEP+, MM/PBSA, MM/GBSA, and Autodock vina, FMOScore showed good performance in prediction accuracies. By constructing a retrospective case study, it was observed that incorporating calculations for solvation free energy and deformation energy can further enhance the precision of FMO predictions for binding affinity. Furthermore, using FMOScore-guided lead optimization against Src homology-2-containing protein tyrosine phosphatase 2 (SHP-2), we discovered a novel and potent allosteric SHP-2 inhibitor (compound 8).

Simulation on the Ligand Leaving Process of Human Heat Shock Protein

Binding Mechanism of Inhibitors to Heat Shock Protein 90 Investigated by Multiple Independent Molecular Dynamics Simulations and Prediction of Binding Free Energy

Electronic Polarization Stabilizes Tertiary Structure Prediction of HP-36.

Dynamics Revelation of Conformational Changes and Binding Modes of Heat Shock Protein 90 Induced by Inhibitor Associations

Effects of Ligand Binding on the Mechanical Stability of Protein GB1 Studied by Steered Molecular Dynamics Simulation

The effect of ligands on the thermal stability of sulfotransferases: a molecular dynamics simulation study

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Conformational dynamics modulate the catalytic activity of the molecular chaperone Hsp90

Binding Modes and Conformational Changes of FK506-binding Protein 51 Induced by Inhibitor Bindings: Insight into Molecular Mechanisms Based on Multiple Simulation Technologies

Design of A Flexible Cell-Based Assay for the Evaluation of Heat Shock Protein 70 Expression Modulators

Thermal Denaturing of Mutant Lysozyme with Both the OPLSAA and the CHARMM Force Fields

Molecular dynamics simulation of the effect of temperature on the conformation of ubiquitin protein

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors

Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Molecular Dynamics Simulation of Protein Conformation Transition in Solution

Spatially and kinetically resolved changes in the conformational dynamics of the Hsp90 chaperone machine

Conformational dynamics of free and membrane-bound human Hsp70 in model cytosolic and endo-lysosomal environments

Comparison of the Ph- and Thermally-Induced Fluctuations of a Therapeutic Antibody Fab Fragment by Molecular Dynamics Simulation

Active-State Structures of a Small Heat-Shock Protein Revealed a Molecular Switch for Chaperone Function.