Revealing the favorable dissociation pathway of type II kinase inhibitors via enhanced sampling simulations and two-end-state calculations

How does a type II inhibitor bind to/unbind from a kinase target is still a confusing question because the small molecule occupies both the ATP pocket and the allosteric pocket of the kinase binding site. Here, by using enhanced sampling simulations (umbrella sampling, US) and two-end-state free energy calculations (MM/GSBA), we systemically studied the dissociation processes of two distinct small molecules escaping from the binding pocket of p38 MAP kinase through the allosteric channel and the ATP channel. The results show that the unbinding pathways along the allosteric channel have much lower PMF depths than those along the ATP channel, suggesting that the allosteric channel is more favorable for the dissociations of the two inhibitors and thereby supporting the general understanding that the largest channel of a target is usually the entry/exit pathway for the binding/dissociation of small molecules. Interestingly, the MM/GBSA approach yielded similar PMF profiles compared with those based on US, a much time consuming approach, indicating that for a general study, such as detecting the important transition state of a ligand binding/unbinding process, MM/GBSA may be a feasible choice.

[1]  Xuri Huang,et al.  MOLECULAR DYNAMICS AND FREE ENERGY ANALYSES OF ERK2–PYRAZOLYLPYRROLE INHIBITORS INTERACTIONS: INSIGHT INTO STRUCTURE-BASED LIGAND DESIGN , 2009 .

[2]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[3]  Tingjun Hou,et al.  Insight into Crizotinib Resistance Mechanisms Caused by Three Mutations in ALK Tyrosine Kinase using Free Energy Calculation Approaches , 2013, J. Chem. Inf. Model..

[4]  P. Kollman,et al.  Binding of a diverse set of ligands to avidin and streptavidin: an accurate quantitative prediction of their relative affinities by a combination of molecular mechanics and continuum solvent models. , 2000, Journal of medicinal chemistry.

[5]  G. Scapin.  Structural biology in drug design: selective protein kinase inhibitors. , 2002, Drug discovery today.

[6]  KumarShankar,et al.  The weighted histogram analysis method for free-energy calculations on biomolecules. I , 1992 .

[7]  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..

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

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

[10]  Tingjun Hou,et al.  P-loop Conformation Governed Crizotinib Resistance in G2032R-Mutated ROS1 Tyrosine Kinase: Clues from Free Energy Landscape , 2014, PLoS Comput. Biol..

[11]  Adrian H Elcock,et al.  Computational sampling of a cryptic drug binding site in a protein receptor: explicit solvent molecular dynamics and inhibitor docking to p38 MAP kinase. , 2006, Journal of molecular biology.

[12]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

[13]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[14]  Stewart A. Adcock,et al.  Molecular dynamics: survey of methods for simulating the activity of proteins. , 2006, Chemical reviews.

[15]  Richard H. Henchman,et al.  Standard Free Energy of Binding from a One-Dimensional Potential of Mean Force. , 2009, Journal of chemical theory and computation.

[16]  D. Barford,et al.  Mechanism of Activation of the RAF-ERK Signaling Pathway by Oncogenic Mutations of B-RAF , 2004, Cell.

[17]  N. Gray,et al.  Rational design of inhibitors that bind to inactive kinase conformations , 2006, Nature chemical biology.

[18]  Piotr Cieplak,et al.  Molecular dynamics and free energy analyses of cathepsin D-inhibitor interactions: insight into structure-based ligand design. , 2002, Journal of medicinal chemistry.

[19]  P. Kollman,et al.  Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA. , 2001, Journal of the American Chemical Society.

[20]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. , 2014, Physical chemistry chemical physics : PCCP.

[21]  P. Cohen.  Protein kinases — the major drug targets of the twenty-first century? , 2002, Nature reviews. Drug discovery.

[22]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[23]  E. Sausville,et al.  Issues and progress with protein kinase inhibitors for cancer treatment , 2003, Nature Reviews Drug Discovery.

[24]  Christophe Chipot,et al.  Standard binding free energies from computer simulations: What is the best strategy? , 2013, Journal of chemical theory and computation.

[25]  Albert C. Pan,et al.  The Dynamic Process of β2-Adrenergic Receptor Activation , 2013, Cell.

[26]  Gabriele Costantino,et al.  Unbinding Pathways of VEGFR2 Inhibitors Revealed by Steered Molecular Dynamics , 2014, J. Chem. Inf. Model..

[27]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[28]  Tingjun Hou,et al.  Molecular dynamics and free energy studies on the wild-type and double mutant HIV-1 protease complexed with amprenavir and two amprenavir-related inhibitors: mechanism for binding and drug resistance. , 2007, Journal of medicinal chemistry.

[29]  P. Seeburg,et al.  Structural mechanism for STI-571 inhibition of abelson tyrosine kinase. , 2000, Science.

[30]  Tingjun Hou,et al.  Molecular Principle of Topotecan Resistance by Topoisomerase I Mutations through Molecular Modeling Approaches , 2013, J. Chem. Inf. Model..

[31]  Eric Darve,et al.  Adaptive biasing force method for scalar and vector free energy calculations. , 2008, The Journal of chemical physics.

[32]  L. Tong,et al.  Inhibition of p38 MAP kinase by utilizing a novel allosteric binding site , 2002, Nature Structural Biology.

[33]  Jörg Weiser,et al.  Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO) , 1999, J. Comput. Chem..

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

[35]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[36]  Jiye Shi,et al.  Exploring transition pathway and free-energy profile of large-scale protein conformational change by combining normal mode analysis and umbrella sampling molecular dynamics. , 2014, The journal of physical chemistry. B.

[37]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[38]  T. Hunter,et al.  The Protein Kinase Complement of the Human Genome , 2002, Science.

[39]  Normand Mousseau,et al.  Use of Umbrella Sampling to Calculate the Entrance/Exit Pathway for Z-Pro-Prolinal Inhibitor in Prolyl Oligopeptidase. , 2011, Journal of chemical theory and computation.

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

[41]  G. Narahari Sastry,et al.  Sequence, Structure, and Active Site Analyses of p38 MAP Kinase: Exploiting DFG-out Conformation as a Strategy to Design New Type II Leads , 2011, J. Chem. Inf. Model..

[42]  Ying Yang,et al.  Molecular Dynamics Simulation and Free Energy Calculation Studies of the Binding Mechanism of Allosteric Inhibitors with p38α MAP Kinase , 2011, J. Chem. Inf. Model..

[43]  Federico Filomia,et al.  Insights into MAPK p38alpha DFG flip mechanism by accelerated molecular dynamics. , 2010, Bioorganic & medicinal chemistry.

[44]  G. Torrie,et al.  Monte Carlo study of a phase‐separating liquid mixture by umbrella sampling , 1977 .

[45]  Hui-Yong Sun,et al.  Finding chemical drugs for genetic diseases. , 2014, Drug discovery today.

[46]  J. Kästner.  Umbrella sampling , 2011 .

[47]  Lin-Li Li,et al.  Steered Molecular Dynamics Simulations Reveal the Likelier Dissociation Pathway of Imatinib from Its Targeting Kinases c-Kit and Abl , 2009, PloS one.

[48]  J. Kuriyan,et al.  The Conformational Plasticity of Protein Kinases , 2002, Cell.

[49]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[50]  Jing Chen,et al.  Free energy landscape for the binding process of Huperzine A to acetylcholinesterase , 2013, Proceedings of the National Academy of Sciences.

[51]  Albert C. Pan,et al.  Activation mechanism of the β2-adrenergic receptor , 2011, Proceedings of the National Academy of Sciences.

[52]  Huanxiang Liu,et al.  Molecular modeling study on the resistance mechanism of HCV NS3/4A serine protease mutants R155K, A156V and D168A to TMC435. , 2012, Antiviral research.

[53]  J. Bajorath,et al.  Quo vadis, virtual screening? A comprehensive survey of prospective applications. , 2010, Journal of medicinal chemistry.

[54]  Christophe Chipot,et al.  Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables. , 2010, Journal of chemical theory and computation.

[55]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

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

[57]  Peter A. Kollman,et al.  Calculating Structures and Free Energies of Complex Molecules: Combining Molecular Mechanics and Continuum Models , 2001 .

[58]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[59]  Tingjun Hou,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. , 2014, Physical chemistry chemical physics : PCCP.

[60]  Zhong-Yi Wang,et al.  Structural and Energetic Analyses of SNPs in Drug Targets and Implications for Drug Therapy , 2013, J. Chem. Inf. Model..