The Effect of Conformational Flexibility on Binding Free Energy Estimation between Kinases and Their Inhibitors

Accurate prediction of binding affinities of drug candidates to their targets remains challenging because of protein flexibility in solution. Conformational flexibility of the ATP-binding site in the CDK2 and ERK2 kinases was identified using molecular dynamics simulations. The binding free energy (ΔG) of twenty-four ATP-competitive inhibitors toward these kinases was assessed using an alchemical free energy perturbation method, MP-CAFEE. However, large calculation errors of 2-3 kcal/mol were observed using this method, where the free energy simulation starts from a single equilibrated conformation. Here, we developed a new ΔG computation method, where the starting structure was set to multiconformations to cover flexibility. The calculation accuracy was successfully improved, especially for larger molecular size compounds, leading to reliable prediction of a broader range of drug candidates. The present study demonstrates that conformational flexibility of interactions between a compound and the glycine-rich loop in the kinases is a key factor in ΔG estimation.

[1]  Andrea Cavalli,et al.  Steered Molecular Dynamics Simulations for Studying Protein-Ligand Interaction in Cyclin-Dependent Kinase 5 , 2014, J. Chem. Inf. Model..

[2]  R. Zwanzig.  High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[3]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[4]  K. Ohno,et al.  Docking study and binding free energy calculation of poly (ADP-ribose) polymerase inhibitors , 2011, Journal of molecular modeling.

[5]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[6]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[7]  Cheng Zhang,et al.  Structure and Function of an Irreversible Agonist-β2 Adrenoceptor complex , 2010, Nature.

[8]  B. Roux,et al.  Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials. , 2006, Biophysical journal.

[9]  Benoît Roux,et al.  Computational Study of Gleevec and G6G Reveals Molecular Determinants of Kinase Inhibitor Selectivity , 2014, Journal of the American Chemical Society.

[10]  V. Pande,et al.  Activation pathway of Src kinase reveals intermediate states as novel targets for drug design , 2014, Nature Communications.

[11]  Hagai Meirovitch,et al.  New method for calculating the absolute free energy of binding: the effect of a mobile loop on the avidin/biotin complex. , 2011, The journal of physical chemistry. B.

[12]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[13]  B. Roux,et al.  Computations of standard binding free energies with molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[14]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[15]  Eric T. Kim,et al.  How does a drug molecule find its target binding site? , 2011, Journal of the American Chemical Society.

[16]  Shantenu Jha,et al.  Computing Clinically Relevant Binding Free Energies of HIV-1 Protease Inhibitors , 2014, Journal of chemical theory and computation.

[17]  Joseph A. Bank,et al.  Supporting Online Material Materials and Methods Figs. S1 to S10 Table S1 References Movies S1 to S3 Atomic-level Characterization of the Structural Dynamics of Proteins , 2022 .

[18]  Benoît Roux,et al.  Free Energy Perturbation Hamiltonian Replica-Exchange Molecular Dynamics (FEP/H-REMD) for Absolute Ligand Binding Free Energy Calculations. , 2010, Journal of chemical theory and computation.

[19]  A. Cavalli,et al.  Role of Molecular Dynamics and Related Methods in Drug Discovery. , 2016, Journal of medicinal chemistry.

[20]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[21]  Haruki Nakamura,et al.  Application of MDGRAPE‐3, a special purpose board for molecular dynamics simulations, to periodic biomolecular systems , 2009, J. Comput. Chem..

[22]  David L Mobley,et al.  The Confine-and-Release Method: Obtaining Correct Binding Free Energies in the Presence of Protein Conformational Change. , 2007, Journal of chemical theory and computation.

[23]  Junichi Goto,et al.  ASEDock-Docking Based on Alpha Spheres and Excluded Volumes , 2008, J. Chem. Inf. Model..

[24]  K. Dill,et al.  Predicting absolute ligand binding free energies to a simple model site. , 2007, Journal of molecular biology.

[25]  K. Schulten,et al.  Calculating potentials of mean force from steered molecular dynamics simulations. , 2004, The Journal of chemical physics.

[26]  Benoît Roux,et al.  Calculation of Standard Binding Free Energies:  Aromatic Molecules in the T4 Lysozyme L99A Mutant. , 2006, Journal of chemical theory and computation.

[27]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[28]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[29]  Susan S. Taylor,et al.  Surface comparison of active and inactive protein kinases identifies a conserved activation mechanism , 2006, Proceedings of the National Academy of Sciences.

[30]  R. Dror,et al.  A conserved protonation-dependent switch controls drug binding in the Abl kinase , 2009, Proceedings of the National Academy of Sciences.

[31]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[32]  Masahiko Nakatsui,et al.  Constructing a Foundational Platform Driven by Japan’s K Supercomputer for Next‐Generation Drug Design , 2014, Molecular informatics.

[33]  Bo Wang,et al.  Molecular Recognition in a Diverse Set of Protein-Ligand Interactions Studied with Molecular Dynamics Simulations and End-Point Free Energy Calculations , 2013, J. Chem. Inf. Model..

[34]  Hiroshi Yamashita,et al.  Prediction of the binding affinity of compounds with diverse scaffolds by MP-CAFEE. , 2013, Biophysical Chemistry.

[35]  Todd J. A. Ewing,et al.  DOCK 4.0: Search strategies for automated molecular docking of flexible molecule databases , 2001, J. Comput. Aided Mol. Des..

[36]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[37]  Nobuyuki Matubayasi,et al.  Evaluation of protein-protein docking model structures using all-atom molecular dynamics simulations combined with the solution theory in the energy representation. , 2012, The Journal of chemical physics.

[38]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[39]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[40]  Richard D. Smith,et al.  CSAR Benchmark Exercise of 2010: Selection of the Protein–Ligand Complexes , 2011, J. Chem. Inf. Model..

[41]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998 .

[42]  R. Altman,et al.  Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways , 2013, Nature chemistry.

[43]  A. Iafrate,et al.  Two Novel ALK Mutations Mediate Acquired Resistance to the Next-Generation ALK Inhibitor Alectinib , 2014, Clinical Cancer Research.

[44]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[45]  Joseph A Morrone,et al.  Advances in free-energy-based simulations of protein folding and ligand binding. , 2016, Current opinion in structural biology.

[46]  Tetsuya J. Kobayashi,et al.  Perturbation analyses of intermolecular interactions. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[47]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[48]  Haruki Nakamura,et al.  Protein-Ligand Binding Free Energy Calculation by the Smooth Reaction Path Generation (SRPG) Method , 2009, J. Chem. Inf. Model..

[49]  Holger Gohlke,et al.  MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[50]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[51]  Michael R. Shirts,et al.  Equilibrium free energies from nonequilibrium measurements using maximum-likelihood methods. , 2003, Physical review letters.

[52]  P. Kollman,et al.  Computational Alanine Scanning To Probe Protein−Protein Interactions: A Novel Approach To Evaluate Binding Free Energies , 1999 .

[53]  Haruki Nakamura,et al.  Computational Study of Drug Binding Affinity to Influenza A Neuraminidase Using Smooth Reaction Path Generation (SRPG) Method , 2015, J. Chem. Inf. Model..

[54]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[55]  H. Fujitani,et al.  Massively parallel computation of absolute binding free energy with well-equilibrated states. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[56]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[57]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .