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).

Protein–ligand binding affinity prediction of cyclin‐dependent kinase‐2 inhibitors by dynamically averaged fragment molecular orbital‐based interaction energy

Decoding molecular mechanism of inhibitor bindings to CDK2 using molecular dynamics simulations and binding free energy calculations

Effect of Substituents in Different Positions of Aminothiazole Hinge-Binding Scaffolds on Inhibitor–CDK2 Association Probed by Interaction Entropy Method

Molecular Dynamic Simulations Give Insight into the Mechanism of Binding Between 2-Aminothiazole Inhibitors and CDK5

Novel structural features of CDK inhibition revealed by an ab initio computational method combined with dynamic simulations

Identification of natural marine compounds as potential inhibitors of CDK2 using molecular docking and molecular dynamics simulation approach

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

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

Computer-assisted Discovery and Evaluation of Potential Ribosomal Protein S6 Kinase Beta 2 Inhibitors

Free energy analysis of CDK2-inhibitor interaction

Exploring the prominent performance of CX-4945 derivatives as protein kinase CK2 inhibitors by a combined computational study.

Nonfitting Protein-Ligand Interaction Scoring Function Based on First-Principles Theoretical Chemistry Methods: Development and Application on Kinase Inhibitors.

Free Energy Profile and Kinetics of Coupled Folding and Binding of the Intrinsically Disordered Protein p53 with MDM2

Modeling Local Structural Rearrangements Using FEP/REST: Application to Relative Binding Affinity Predictions of CDK2 Inhibitors

Use of the Multilayer Fragment Molecular Orbital Method to Predict the Rank Order of Protein–Ligand Binding Affinities: A Case Study Using Tankyrase 2 Inhibitors

Potent inhibitors targeting Cyclin-dependent kinase 9 discovered via virtual high-throughput screening and absolute binding free energy calculations

Performance of protein-ligand docking with CDK4/6 inhibitors: a case study

Machine Learning Meets Physics-based Modeling: A Mass-spring System to Predict Protein-ligand Binding Affinity

Binding Mechanism of CDK5 with Roscovitine Derivatives Based on Molecular Dynamics Simulations and MM/PBSA Methods.

Design and screening of FAK, CDK 4/6 dual inhibitors by pharmacophore model, molecular docking, and molecular dynamics simulation

A Computational Analysis of the Binding Affinities of FKBP12 Inhibitors Using the MM‐PB/SA Method