Abstract:BackgroundThe reliable and robust estimation of ligand binding affinity continues to be a challenge in drug design. Many current methods rely on molecular mechanics (MM) calculations which do not fully explain complex molecular interactions. Full quantum mechanical (QM) computation of the electronic state of protein-ligand complexes has recently become possible by the latest advances in the development of linear-scaling QM methods such as the ab initio fragment molecular orbital (FMO) method. This approximate molecular orbital method is sufficiently fast that it can be incorporated into the development cycle during structure-based drug design for the reliable estimation of ligand binding affinity. Additionally, the FMO method can be combined with approximations for entropy and solvation to make it applicable for binding affinity prediction for a broad range of target and chemotypes.ResultsWe applied this method to examine the binding affinity for a series of published cyclin-dependent kinase 2 (CDK2) inhibitors. We calculated the binding affinity for 28 CDK2 inhibitors using the abinitio FMO method based on a number of X-ray crystal structures. The sum of the pair interaction energies (PIE) was calculated and used to explain the gas-phase enthalpic contribution to binding. The correlation of the ligand potencies to the protein-ligand interaction energies gained from FMO was examined and was seen to give a good correlation which outperformed three MM force field based scoring functions used to appoximate the free energy of binding. Although the FMO calculation allows for the enthalpic component of binding interactions to be understood at the quantum level, as it is an in vacuo single point calculation, the entropic component and solvation terms are neglected. For this reason a more accurate and predictive estimate for binding free energy was desired. Therefore, additional terms used to describe the protein-ligand interactions were then calculated to improve the correlation of the FMO derived values to experimental free energies of binding. These terms were used to account for the polar and non-polar solvation of the molecule estimated by the Poisson-Boltzmann equation and the solvent accessible surface area (SASA), respectively, as well as a correction term for ligand entropy. A quantitative structure-activity relationship (QSAR) model obtained by Partial Least Squares projection to latent structures (PLS) analysis of the ligand potencies and the calculated terms showed a strong correlation (r2 = 0.939, q2 = 0.896) for the 14 molecule test set which had a Pearson rank order correlation of 0.97. A training set of a further 14 molecules was well predicted (r2 = 0.842), and could be used to obtain meaningful estimations of the binding free energy.ConclusionsOur results show that binding energies calculated with the FMO method correlate well with published data. Analysis of the terms used to derive the FMO energies adds greater understanding to the binding interactions than can be gained by MM methods. Combining this information with additional terms and creating a scaled model to describe the data results in more accurate predictions of ligand potencies than the absolute values obtained by FMO alone.

Rigorous Free Energy Perturbation Approach to Estimating Relative Binding Affinities Between Ligands with Multiple Protonation and Tautomeric States.

Accurate Calculation of Relative Binding Free Energies Between Ligands with Different Net Charges

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

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations

Accurate calculation of absolute protein-ligand binding free energies

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

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

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

The maximal and current accuracy of rigorous protein-ligand binding free energy calculations

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Multistate Method to Efficiently Account for Tautomerism and Protonation in Alchemical Free-Energy Calculations

Accurate Estimation Of Protein-Ligand Binding Free Energies Based On Geometric Restraints

Absolute Binding Free Energy Calculation and Design of a Subnanomolar Inhibitor of Phosphodiesterase-10

Binding Free Energy Estimation for Protein-Ligand Complex Based on MM-PBSA with Various Partial Charge Models.

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

Validation of the Alchemical Transfer Method for the Estimation of Relative Binding Affinities of Molecular Series

Prediction of cyclin-dependent kinase 2 inhibitor potency using the fragment molecular orbital method

Enhancing Hit Discovery in Virtual Screening through Absolute Protein-Ligand Binding Free-Energy Calculations

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations

Non-bonded force field parameters from MBIS partitioning of the molecular electron density improve binding affinity predictions of the T4-lysozyme double mutant