Predicting resistance of clinical Abl mutations to targeted kinase inhibitors using alchemical free-energy calculations

The therapeutic effect of targeted kinase inhibitors can be significantly reduced by intrinsic or acquired resistance mutations that modulate the affinity of the drug for the kinase. In cancer, the majority of missense mutations are rare, making it difficult to predict their impact on inhibitor affinity. We examine the potential for alchemical free-energy calculations to predict how kinase mutations modulate inhibitor affinities to Abl, a major target in chronic myelogenous leukemia (CML). These calculations have useful accuracy in predicting resistance for eight FDA-approved kinase inhibitors across 144 clinically identified point mutations, with a root mean square error in binding free-energy changes of $$1.1_{0.9}^{1.3}$$1.10.91.3 kcal mol−1 (95% confidence interval) and correctly classifying mutations as resistant or susceptible with $$88_{82}^{93}$$888293% accuracy. This benchmark establishes the potential for physical modeling to collaboratively support the assessment and anticipation of patient mutations to affect drug potency in clinical applications.Kevin Hauser et al. accurately predict the impact of mutations in a kinase on the binding affinities of targeted kinase inhibitors using alchemical free-energy calculations. With 88% accuracy, resistance or sensitivity to therapy is computed for 144 clinically-identified point mutations in this major target in chronic myelogenous leukemia.

[1]  Xiaoyun Lu,et al.  Recent developments in the third generation inhibitors of Bcr-Abl for overriding T315I mutation. , 2011, Current Medicinal Chemistry.

[2]  P. Mahadevan,et al.  An overview , 2007, Journal of Biosciences.

[3]  Krister Wennerberg,et al.  Axitinib effectively inhibits BCR-ABL1(T315I) with a distinct binding conformation , 2015, Nature.

[4]  Jing Huang,et al.  CHARMM36 all‐atom additive protein force field: Validation based on comparison to NMR data , 2013, J. Comput. Chem..

[5]  David A. Pearlman,et al.  A Comparison of Alternative Approaches to Free Energy Calculations , 1994 .

[6]  L. O’Driscoll,et al.  Neuromedin U: a candidate biomarker and therapeutic target to predict and overcome resistance to HER-tyrosine kinase inhibitors. , 2014, Cancer research.

[7]  David L Mobley,et al.  Accurate and efficient corrections for missing dispersion interactions in molecular simulations. , 2007, The journal of physical chemistry. B.

[8]  James D. Griffin,et al.  Second generation inhibitors of BCR-ABL for the treatment of imatinib-resistant chronic myeloid leukaemia , 2007, Nature Reviews Cancer.

[9]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[10]  Matthew P. Jacobson,et al.  A Molecular Mechanics Approach to Modeling Protein-Ligand Interactions: Relative Binding Affinities in Congeneric Series , 2011, J. Chem. Inf. Model..

[11]  Tongqing Zhou,et al.  Free Energy Perturbation Calculation of Relative Binding Free Energy between Broadly Neutralizing Antibodies and the gp120 Glycoprotein of HIV-1 , 2017, Journal of molecular biology.

[12]  W L Jorgensen,et al.  Computational binding studies of human pp60c-src SH2 domain with a series of nonpeptide, phosphophenyl-containing ligands. , 2000, Bioorganic & medicinal chemistry letters.

[13]  Julien Michel,et al.  Protein-ligand binding affinity predictions by implicit solvent simulations: a tool for lead optimization? , 2006, Journal of medicinal chemistry.

[14]  A. McKenna,et al.  The genetic landscape of clinical resistance to RAF inhibition in metastatic melanoma. , 2014, Cancer discovery.

[15]  A. Czarnecka,et al.  Resistance to tyrosine kinase inhibitors in clear cell renal cell carcinoma: from the patient's bed to molecular mechanisms. , 2014, Biochimica et biophysica acta.

[16]  Berk Hess,et al.  Lennard-Jones Lattice Summation in Bilayer Simulations Has Critical Effects on Surface Tension and Lipid Properties. , 2013, Journal of chemical theory and computation.

[17]  Michael J. Bodkin,et al.  Accurate calculation of the absolute free energy of binding for drug molecules , 2015, Chemical science.

[18]  Pengyu Y. Ren,et al.  Calculation of protein–ligand binding free energy by using a polarizable potential , 2008, Proceedings of the National Academy of Sciences.

[19]  Donna Neuberg,et al.  Characterization of AMN107, a selective inhibitor of native and mutant Bcr-Abl. , 2005, Cancer cell.

[20]  T. Clackson,et al.  AP24534, a pan-BCR-ABL inhibitor for chronic myeloid leukemia, potently inhibits the T315I mutant and overcomes mutation-based resistance. , 2009, Cancer cell.

[21]  I. Mikkola,et al.  A novel Bcr-Abl splice isoform is associated with the L248V mutation in CML patients with acquired resistance to imatinib , 2006, Leukemia.

[22]  Julien Michel,et al.  Effects of Water Placement on Predictions of Binding Affinities for p38α MAP Kinase Inhibitors. , 2010, Journal of chemical theory and computation.

[23]  Pratyush Tiwary,et al.  How a Kinase Inhibitor Withstands Gatekeeper Residue Mutations. , 2016, Journal of the American Chemical Society.

[24]  Donavan T. Cheng,et al.  Mutational Landscape of Metastatic Cancer Revealed from Prospective Clinical Sequencing of 10,000 Patients , 2017, Nature Medicine.

[25]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[26]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[27]  Tudor I. Oprea,et al.  A comprehensive map of molecular drug targets , 2016, Nature Reviews Drug Discovery.

[28]  Sarah A Heerboth,et al.  Drug Resistance in Cancer: An Overview , 2014, Cancers.

[29]  Daniel Cappel,et al.  Relative Binding Free Energy Calculations Applied to Protein Homology Models , 2016, J. Chem. Inf. Model..

[30]  T. Clackson,et al.  Structural Mechanism of the Pan‐BCR‐ABL Inhibitor Ponatinib (AP24534): Lessons for Overcoming Kinase Inhibitor Resistance , 2011, Chemical biology & drug design.

[31]  G. Superti-Furga,et al.  Crystal structure of an SH2-kinase construct of c-Abl and effect of the SH2 domain on kinase activity. , 2015, The Biochemical journal.

[32]  David L Mobley,et al.  Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: an accurate correction scheme for electrostatic finite-size effects. , 2013, The Journal of chemical physics.

[33]  S. Boxer,et al.  Structural and Spectroscopic Analysis of the Kinase Inhibitor Bosutinib and an Isomer of Bosutinib Binding to the Abl Tyrosine Kinase Domain , 2012, PloS one.

[34]  Teresa Head-Gordon,et al.  Advanced potential energy surfaces for condensed phase simulation. , 2014, Annual review of physical chemistry.

[35]  Atli Thorarensen,et al.  Imidazotriazines: Spleen Tyrosine Kinase (Syk) Inhibitors Identified by Free‐Energy Perturbation (FEP) , 2016, ChemMedChem.

[36]  T. Clackson,et al.  Comprehensive Analysis Of The In Vitro Potency Of Ponatinib, and All Other Approved BCR-ABL Tyrosine Kinase Inhibitors (TKIs), Against a Panel Of Single and Compound BCR-ABL Mutants , 2013 .

[37]  Jennifer L. Knight,et al.  OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. , 2016, Journal of chemical theory and computation.

[38]  R. Levine,et al.  Molecular Pathways Molecular Pathways : Molecular Basis for Sensitivity and Resistance to JAK Kinase Inhibitors , 2014 .

[39]  Y. Yosaatmadja,et al.  1.85 angstrom structure of EGFR kinase domain with gefitinib , 2014 .

[40]  Dana S. Neel,et al.  Resistance is futile: overcoming resistance to targeted therapies in lung adenocarcinoma , 2017, npj Precision Oncology.

[41]  K. Shokat,et al.  Targeting the cancer kinome through polypharmacology , 2010, Nature Reviews Cancer.

[42]  T. Mikkelsen,et al.  Comprehensive mutational scanning of a kinase in vivo reveals substrate-dependent fitness landscapes , 2014, Nucleic acids research.

[43]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

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

[45]  M. Baccarani,et al.  Contribution of ABL Kinase Domain Mutations to Imatinib Resistance in Different Subsets of Philadelphia-Positive Patients: By the GIMEMA Working Party on Chronic Myeloid Leukemia , 2006, Clinical Cancer Research.

[46]  Yvonne C. Martin,et al.  Let’s not forget tautomers , 2009, J. Comput. Aided Mol. Des..

[47]  M. Deininger,et al.  Bcr-Abl kinase domain mutations, drug resistance, and the road to a cure for chronic myeloid leukemia. , 2007, Blood.

[48]  A. Iafrate,et al.  Identification of genotype-correlated sensitivity to selective kinase inhibitors by using high-throughput tumor cell line profiling , 2007, Proceedings of the National Academy of Sciences.

[49]  John Salvatier,et al.  Probabilistic programming in Python using PyMC3 , 2016, PeerJ Comput. Sci..

[50]  Jeremy R. Greenwood,et al.  Epik: a software program for pKa prediction and protonation state generation for drug-like molecules , 2007, J. Comput. Aided Mol. Des..

[51]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[52]  S. Laufer,et al.  Fighting cancer drug resistance: Opportunities and challenges for mutation-specific EGFR inhibitors. , 2015, Drug resistance updates : reviews and commentaries in antimicrobial and anticancer chemotherapy.

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

[54]  Webster K. Cavenee,et al.  Precision cancer therapy is impacted by oncogene-dependent epigenome remodeling , 2017, npj Precision Oncology.

[55]  Jan H. Jensen,et al.  Calculating pH and salt dependence of protein-protein binding. , 2008, Current pharmaceutical biotechnology.

[56]  A. Drilon,et al.  Structural insight into selectivity and resistance profiles of ROS1 tyrosine kinase inhibitors , 2015, Proceedings of the National Academy of Sciences.

[57]  Sarat Chandarlapaty,et al.  AKT inhibition relieves feedback suppression of receptor tyrosine kinase expression and activity. , 2011, Cancer cell.

[58]  Jennifer L. Knight,et al.  Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. , 2015, Journal of the American Chemical Society.

[59]  Meining Wang,et al.  Alectinib: a novel second generation anaplastic lymphoma kinase (ALK) inhibitor for overcoming clinically-acquired resistance , 2015, Acta pharmaceutica Sinica. B.

[60]  K. Sauer,et al.  Mechanisms of drug resistance in kinases , 2011, Expert opinion on investigational drugs.

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

[62]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[63]  Lihua Huang,et al.  Mechanisms of resistance to EGFR tyrosine kinase inhibitors , 2015, Acta pharmaceutica Sinica. B.

[64]  W. Sherman,et al.  Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. , 2010, Journal of chemical theory and computation.

[65]  A. Roitberg,et al.  Long-Time-Step Molecular Dynamics through Hydrogen Mass Repartitioning. , 2015, Journal of chemical theory and computation.

[66]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[67]  Robert Abel,et al.  Accelerating drug discovery through tight integration of expert molecular design and predictive scoring. , 2017, Current opinion in structural biology.

[68]  Ravi Radhakrishnan,et al.  Erlotinib binds both inactive and active conformations of the EGFR tyrosine kinase domain , 2012, The Biochemical journal.

[69]  B. Taylor,et al.  Implementing Genome-Driven Oncology , 2017, Cell.

[70]  Mindy I. Davis,et al.  Comprehensive analysis of kinase inhibitor selectivity , 2011, Nature Biotechnology.

[71]  D. Raleigh,et al.  Experimental and Computational Analysis of Protein Stabilization by Gly-to-d-Ala Substitution: A Convolution of Native State and Unfolded State Effects. , 2016, Journal of the American Chemical Society.

[72]  G. Superti-Furga,et al.  Structural Basis for the Autoinhibition of c-Abl Tyrosine Kinase , 2003, Cell.

[73]  A. Redig,et al.  Basket trials and the evolution of clinical trial design in an era of genomic medicine. , 2015, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[74]  J. Kuriyan,et al.  Multiple BCR-ABL kinase domain mutations confer polyclonal resistance to the tyrosine kinase inhibitor imatinib (STI571) in chronic phase and blast crisis chronic myeloid leukemia. , 2002, Cancer cell.

[75]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

[76]  Melissa Coates Ford,et al.  Examining the Feasibility of Using Free Energy Perturbation (FEP+) in Predicting Protein Stability , 2017, J. Chem. Inf. Model..

[77]  P. Johnston,et al.  Cancer drug resistance: an evolving paradigm , 2013, Nature Reviews Cancer.

[78]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[79]  B. Druker,et al.  Combined Abl Inhibitor Therapy for Minimizing Drug Resistance in Chronic Myeloid Leukemia: Src/Abl Inhibitors Are Compatible with Imatinib , 2005, Clinical Cancer Research.

[80]  Lingle Wang,et al.  On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities , 2012, Proceedings of the National Academy of Sciences.

[81]  Woody Sherman,et al.  Predicting the Effect of Amino Acid Single-Point Mutations on Protein Stability-Large-Scale Validation of MD-Based Relative Free Energy Calculations. , 2017, Journal of molecular biology.

[82]  Carey-Ann D. Burnham,et al.  Evaluation of Machine Learning and Rules-Based Approaches for Predicting Antimicrobial Resistance Profiles in Gram-negative Bacilli from Whole Genome Sequence Data , 2016, Front. Microbiol..

[83]  Emil Alexov,et al.  Protonation and pK changes in protein–ligand binding , 2013, Quarterly Reviews of Biophysics.