Structural and dynamic insights into the energetics of activation loop rearrangement in FGFR1 kinase

Exploring the Impact of Protein Chain Selection in Binding Energy Calculations with DFT

Calculation of binding free energies between a protein and a ligand are highly desired for computer-aided drug design. Here we approximate the binding energies of ABL1, an enzyme which is the target for drugs used in the treatment of chronic myeloid leukaemia, with minimal models and density functional theory (DFT). Starting from the crystal structures of protein-drug complexes, we estimated the binding free energies having used all available individual molecules (protein chains) within each structure, not only a single one as commonly used, in order to see if the choice of the protein chain is important in such calculations. Differences were observed between chains in the same file. Energy decomposition analysis (EDA) revealed that the most important factors for binding were exchange, repulsion and electrostatics. The desolvation term varied dramatically between the inhibitors (between 4.2 and 92.3 kcal/mol). All functionals showed similar patterns in the EDA and in discriminating between the ligands. Non-covalent interactions (NCI) analysis was used to further explain the differences between protein chains and functionals. Overall, it is shown that small minimal models of a drug binding site can be useful to infer on the suitability of an initial crystal structure for further analysis such as EDA.

Sanggenol B, a plant bioactive, as a safer alternative to tackle cancer by antagonising human FGFR

Fibroblast Growth Receptor Factor (FGFR) are a family of proteins which are, in addition to their biological role, are involved in various pathological functions, such as cancer cellular proliferation, and metastasis. Deregulation of FGFRs at various points could result in malignancy. A conformational transition of the DFG (Asp-Phe-Gly) motif can switch the enzyme from a catalytically active (DFG-in) to an inactive (DFG-out) state. There are a few FDFR inhibitors which have received approval from the FDA, but these have adverse side effects. Hence, there is a demand for safer alternatives. With this aim, Ligand and Structure based virtual screening was carried to identify suitable lead molecule. In this process, Four Featured atom-based 3D Pharmacophore with quantitative structure-activity relationship analysis (3D-QSAR) was developed. The External validation of the hypothesis was carried invoking criteria such as Area under the ROC curve. Natural plant compound databases such as the Traditional Chinese medicine, NPACT and the ZINC Natural databases were chosen for pharmacophore filtering, which was followed by virtual screening against FGFR isoforms. The compound Sanggenol B was identified as the most suitable lead molecule. Structural stability of the protein-ligand complex and interactions of the ligand (Sanggenol B & the reference compound Ponatinib) with FGFR were analysed for 1000 ns (triplicate) by means of molecular simulation and the binding free energy was calculated using MMGBSA. Sanggenol B (PubChem CID: 15233694) binds effectively at the active site with favourable energies and is proposed as a safe alternative from a natural source.Communicated by Ramaswamy H. Sarma.

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic subunit of protein kinase A (PKA-C). This long-range synergistic action is involved in substrate recognition and fidelity, and it may also regulate PKA's association with regulatory subunits and other binding partners. To date, a complete understanding of this intramolecular mechanism is still lacking. Here, we integrated NMR(Nuclear Magnetic Resonance)-restrained molecular dynamics simulations and a Markov State Model to characterize the free energy landscape and conformational transitions of PKA-C. We found that the apoenzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the αC-β4 loop. We validated the second excited state by analyzing the F100A mutant of PKA-C, assessing its structural response to ATP and substrate binding. While PKA-CF100A preserves its catalytic efficiency with Kemptide, this mutation rearranges the αC-β4 loop conformation, interrupting the coupling of the two lobes and abolishing the allosteric binding cooperativity. The highly conserved αC-β4 loop emerges as a pivotal element to control the synergistic binding of nucleotide and substrate, explaining how mutations or insertions near or within this motif affect the function and drug sensitivity in homologous kinases.

Investigating the potential of mono-chalcone compounds in targeting breast cancer receptors through network pharmacology, molecular docking, molecular dynamics simulation, antiproliferative effects, and gene expressions

The study aims to investigate various aspects of synthesized mono-chalcone compounds 5 and 8 concerning breast cancer, including network pharmacology, molecular docking, molecular dynamics (MD) simulations, antiproliferative effects, and gene expressions. Initially, the compounds underwent a network pharmacology analysis targeting breast cancer-related targets, with MalaCards, SwissTargetPrediction, and PharmMapper identifying 70 breast cancer target receptors. Subsequently, protein-protein interaction (PPI) network analysis revealed two distinct target gene clusters. Survival analysis identified seven significant target genes following Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment and Gene Ontology (GO) evaluation. Molecular docking and MD simulations were conducted on these seven target genes (AKT2, BRAF, ESR1, FGFR1, IGF1, IGF1R, and KIT), revealing that compound 8 exhibited the highest binding affinities, as well as better stability and compactness when interacting with the targeted proteins. Next, the compounds underwent cell viability assay and gene expression analysis to validate the in silico findings. Both compounds demonstrated the ability to suppress breast cancer proliferation, with compound 8 showing increased selectivity in targeting breast cancer cells while causing minimal harm to normal breast cells. The suppression of breast cancer cell proliferation was attributed to decreased expression levels of AKT2, BRAF, FGFR1, IGF1, IGF1R, KIT, and ESR1. Hence, the results provide insights into the molecular interaction responsible for the anti-breast cancer capabilities of mono-chalcone compounds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03991-y.

Molecular dynamics simulations reveal phosphorylation-induced conformational dynamics of the fibroblast growth factor receptor 1 kinase

The Fibroblast Growth Factor Receptor1 (FGFR1) kinase wields exquisite control on cell fate, proliferation, differentiation, and homeostasis. An imbalance of FGFR1 signaling leads to several pathogeneses of diseases ranging from multiple cancers to allergic and neurodegenerative disorders. In this study, we investigated the phosphorylation-induced conformational dynamics of FGFR1 in apo and ATP-bound states via all-atom molecular dynamics simulations. All simulations were performed for 2 × 2 µs. We have also investigated the energetics of the binding of ATP to FGFR1 using the molecular mechanics Poisson-Boltzmann scheme. Our study reveals that the FGFR1 kinase can reach a fully active configuration through phosphorylation and ATP binding. A 3-10 helix formation in the activation loop signifies its rearrangement leading to stability upon ATP binding. The interaction of phosphorylated tyrosine (pTyr654) with positively charged residues forms strong salt-bridge interactions, driving the compactness of the structure. The dynamic cross-correlation map reveals phosphorylation enhances correlated motions and reduces anti-correlated motions between different domains. We believe that the mechanistic understanding of large-conformational changes upon the activation of the FGFR1 kinase will aid the development of novel targeted therapeutics.Communicated by Ramaswamy H. Sarma.

Discovery and Characterization of a Bicyclic Peptide (Bicycle) Binder to Thymic Stromal Lymphopoietin

Thymic stromal lymphopoietin (TSLP) is an epithelial-derived pro-inflammatory cytokine involved in the development of asthma and other atopic diseases. We used Bicycle Therapeutics' proprietary phage display platform to identify bicyclic peptides (Bicycles) with high affinity for TSLP, a target that is difficult to drug with conventional small molecules due to the extended protein-protein interactions it forms with both receptors. The hit series was shown to bind to TSLP in a hotspot, that is also used by IL-7Rα. Guided by the first X-ray crystal structure of a small peptide binding to TSLP and the identification of key metabolites, we were able to improve the proteolytic stability of this series in lung S9 fractions without sacrificing binding affinity. This resulted in the potent Bicycle 46 with nanomolar affinity to TSLP (KD = 13 nM), low plasma clearance of 6.4 mL/min/kg, and an effective half-life of 46 min after intravenous dosing to rats.

Network pharmacology, molecular docking and molecular dynamics simulation of chalcone scaffold-based compounds targeting breast cancer receptors

Compounds with a chalcone scaffold-based structure have demonstrated promising anticancer biological activity. However, the molecular interactions between chalcone scaffold-based compounds and breast cancer-associated proteins remain unclear. Through network pharmacology, molecular docking, and molecular dynamics (MD) simulation analyses, compounds with a chalcone scaffold-based structure were evaluated for their interaction with potential breast cancer targets. The compounds were retrieved from the ASINEX database, resulting in 575,302 compounds. A total of 342 compounds with chalcone scaffold-based structures were discovered. From the 342 compounds that was analysed, ten were chosen due to their adherence to Lipinski's rule, having an appropriate range of lipophilicity (LOGP), and topological polar surface area (TPSA), and absence of any toxicity. Based on target intersection, 50 target genes were found and subjected to protein-protein interaction (PPI), gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Four target genes were found to be involved in the breast cancer pathway. Consequently, molecular docking was utilised to analyse the molecular interactions between the compounds and four target protein receptors. Compound 211 exhibited the highest binding affinities for the epidermal growth factor receptor (EGFR), fibroblast growth factor receptor 1 (FGFR1), oestrogen receptor (ESR1), and cyclin dependent kinase 6 (CDK6) with values of -8.95 kcal/mol, -8.60 kcal/mol, -10.33 kcal/mol, and -9.90 kcal/mol, respectively. During MD simulation, compound 211 and its respective proteins were stable, compact, and had minimal flexibility. The findings provide foundations for future studies into the interaction underlying the anti-breast cancer potential of compounds with chalcone-based scaffold structures.Communicated by Ramaswamy H. Sarma.

A knowledge-guided pre-training framework for improving molecular representation learning

Learning effective molecular feature representation to facilitate molecular property prediction is of great significance for drug discovery. Recently, there has been a surge of interest in pre-training graph neural networks (GNNs) via self-supervised learning techniques to overcome the challenge of data scarcity in molecular property prediction. However, current self-supervised learning-based methods suffer from two main obstacles: the lack of a well-defined self-supervised learning strategy and the limited capacity of GNNs. Here, we propose Knowledge-guided Pre-training of Graph Transformer (KPGT), a self-supervised learning framework to alleviate the aforementioned issues and provide generalizable and robust molecular representations. The KPGT framework integrates a graph transformer specifically designed for molecular graphs and a knowledge-guided pre-training strategy, to fully capture both structural and semantic knowledge of molecules. Through extensive computational tests on 63 datasets, KPGT exhibits superior performance in predicting molecular properties across various domains. Moreover, the practical applicability of KPGT in drug discovery has been validated by identifying potential inhibitors of two antitumor targets: hematopoietic progenitor kinase 1 (HPK1) and fibroblast growth factor receptor 1 (FGFR1). Overall, KPGT can provide a powerful and useful tool for advancing the artificial intelligence (AI)-aided drug discovery process.

Optimization of a Protein-Targeted Medicine into an RNA-Specific Small Molecule

Protein-targeted small molecule medicines often bind RNAs and affect RNA-mediated pathways in cells. Historically, small molecule engagement and modulation of RNA have not been considered in medicine development; however, RNA should be considered both a potential on- and off-target. Kinase inhibitors have emecrged as common RNA binders with dovitinib, a classic receptor tyrosine kinase (RTK) inhibitor, inhibiting RTKs and the biogenesis of oncogenic microRNA-21 through direct engagement. In this study, we use knowledge of the molecular recognition of both protein and RNA targets by dovitinib to design molecules that specifically inhibit the RNA target but lack activity against canonical protein targets in cells. As it is now becoming apparent that RNA can be both an on- and off-target for small molecule medicines, this study lays a foundation to use design principles to maximize desired compound activity while minimizing off-target effects.

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic (C) subunit of protein kinase A (PKA). Not only this long-range synergistic action is involved in substrate recognition and fidelity, but it is likely to regulate PKA association with regulatory subunits and other binding partners. To date, a complete understanding of the molecular determinants for this intramolecular mechanism is still lacking. Here, we used an integrated NMR-restrained molecular dynamics simulations and a Markov Model to characterize the free energy landscape and conformational transitions of the catalytic subunit of protein kinase A (PKA-C). We found that the apo-enzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the tip of the αC-β4 loop. To experimentally validate the second excited state, we mutated F100 into alanine and used NMR spectroscopy to characterize the binding thermodynamics and structural response of ATP and a prototypical peptide substrate. While the activity of PKA-CF100A toward a prototypical peptide substrate is unaltered and the enzyme retains its affinity for ATP and substrate, this mutation rearranges the αC-β4 loop conformation interrupting the allosteric coupling between nucleotide and substrate. The highly conserved αC-β4 loop emerges as a pivotal element able to modulate the synergistic binding between nucleotide and substrate and may affect PKA signalosome. These results may explain how insertion mutations within this motif affect drug sensitivity in other homologous kinases.

Ponatinib: An update on its drug targets, therapeutic potential and safety

Leukemia is a malignancy of the hematopoietic system, and as its pathogenesis has become better understood, three generations of tyrosine kinase inhibitors (TKIs) have been developed. Ponatinib is the third-generation breakpoint cluster region (BCR) and Abelson (ABL) TKI, which has been influential in the leukemia therapy for a decade. Moreover, ponatinib is a potent multi-target kinase inhibitor that acts on various kinases, such as KIT, RET, and Src, making it a promising treatment option for triple-negative breast cancer (TNBC), lung cancer, myeloproliferative syndrome, and other diseases. The drug's significant cardiovascular toxicity poses a significant challenge to its clinical use, requiring the development of strategies to minimize its toxicity and side effects. In this article, the pharmacokinetics, targets, therapeutic potential, toxicity and production mechanism of ponatinib will be reviewed. Furthermore, we will discuss methods to reduce the drug's toxicity, providing new avenues for research to improve its safety in clinical use.

The FGF/FGFR signalling mediated anti-cancer drug resistance and therapeutic intervention

ABSTRACT Fibroblast Growth Factor (FGF) ligands and their receptors are crucial factors driving chemoresistance in several malignancies, challenging the efficacy of currently available anti-cancer drugs. The Fibroblast growth factor/receptor (FGF/FGFR) signalling malfunctions in tumor cells, resulting in a range of molecular pathways that may impact its drug effectiveness. Deregulation of cell signalling is critical since it can enhance tumor growth and metastasis. Overexpression and mutation of FGF/FGFR induce regulatory changes in the signalling pathways. Chromosomal translocation facilitating FGFR fusion production aggravates drug resistance. Apoptosis is inhibited by FGFR-activated signalling pathways, reducing multiple anti-cancer medications' destructive impacts. Angiogenesis and epithelial-mesenchymal transition (EMT) are facilitated by FGFRs-dependent signalling, which correlates with drug resistance and enhances metastasis. Further, lysosome-mediated drug sequestration is another prominent method of resistance. Inhibition of FGF/FGFR by following a plethora of therapeutic approaches such as covalent and multitarget inhibitors, ligand traps, monoclonal antibodies, recombinant FGFs, combination therapy, and targeting lysosomes and micro RNAs would be helpful. As a result, FGF/FGFR suppression treatment options are evolving nowadays. To increase positive impacts, the processes underpinning the FGF/FGFR axis' role in developing drug resistance need to be clarified, emphasizing the need for more studies to develop novel therapeutic options to address this significant problem. Communicated by Ramaswamy H. Sarma.

Online Fully Automated System for Hydrogen/Deuterium-Exchange Mass Spectrometry with Millisecond Time Resolution

Amide hydrogen/deuterium-exchange mass spectrometry (HDX-MS) is a powerful tool for analyzing the conformational dynamics of proteins in a solution. Current conventional methods have a measurement limit starting from several seconds and are solely reliant on the speed of manual pipetting or a liquid handling robot. Weakly protected regions of polypeptides, such as in short peptides, exposed loops and intrinsically disordered the protein exchange on the millisecond timescale. Typical HDX methods often cannot resolve the structural dynamics and stability in these cases. Numerous academic laboratories have demonstrated the considerable utility of acquiring HDX-MS data in the sub-second regimes. Here, we describe the development of a fully automated HDX-MS apparatus to resolve amide exchange on the millisecond timescale. Like conventional systems, this instrument boasts automated sample injection with software selection of labeling times, online flow mixing and quenching, while being fully integrated with a liquid chromatography-MS system for existing standard "bottom-up" workflows. HDX-MS's rapid exchange kinetics of several peptides demonstrate the repeatability, reproducibility, back-exchange, and mixing kinetics achieved with the system. Comparably, peptide coverage of 96.4% with 273 peptides was achieved, supporting the equivalence of the system to standard robotics. Additionally, time windows of 50 ms-300 s allowed full kinetic transitions to be observed for many amide groups; especially important are short time points (50-150 ms) for regions that are likely highly dynamic and solvent- exposed. We demonstrate that information on structural dynamics and stability can be measured for stretches of weakly stable polypeptides in small peptides and in local regions of a large enzyme, glycogen phosphorylase.

To Investigate Growth Factor Receptor Targets and Generate Cancer Targeting Inhibitors

Receptor tyrosine kinase (RTK) regulates multiple pathways, including Mitogenactivated protein kinases (MAPKs), PI3/AKT, JAK/STAT pathway, etc. which has a significant role in the progression and metastasis of tumor. As RTK activation regulates numerous essential bodily processes, including cell proliferation and division, RTK dysregulation has been identified in many types of cancers. Targeting RTK is a significant challenge in cancer due to the abnormal upregulation and downregulation of RTK receptors subfamily EGFR, FGFR, PDGFR, VEGFR, and HGFR in the progression of cancer, which is governed by multiple RTK receptor signalling pathways and impacts treatment response and disease progression. In this review, an extensive focus has been carried out on the normal and abnormal signalling pathways of EGFR, FGFR, PDGFR, VEGFR, and HGFR and their association with cancer initiation and progression. These are explored as potential therapeutic cancer targets and therefore, the inhibitors were evaluated alone and merged with additional therapies in clinical trials aimed at combating global cancer.

Electronic and structural study of T315I mutated form in DFG-out conformation of BCR-ABL inhibitors

In this work, the four main drugs for the treatment of chronic myeloid leukemia were analyzed, being imatinib, dasatinib, nilotinib and ponatinib followed by four derivative molecules of nilotinib and ponatinib. For these derivative molecules, the fluorine atoms were replaced by hydrogen and chlorine atoms in order to shade light to the structural effects on this set of inhibitors. Electronic studies were performed at density functional theory level with the B3LYP functional and 6-311+G(d,p) basis set. The frontier molecular orbitals, gap HOMO-LUMO, and NBO were analyzed and compared to docking studies for mutant T315I tyrosine kinase protein structure code 3IK3, in the DFG-out conformation. Structural similarities were pointed out, such as the presence of groups common to all inhibitors and modifications raised up on new generations of imatinib-based inhibitors. One of them is the trifluoromethyl group present in nilotinib and later included in ponatinib, in addition to the 1-methylpiperazin-1-ium group that is present in imatinib and ponatinib. The frontier molecular orbitals of imatinib and ponatinib are contributing to the same amino acid residues, and the ineffectiveness of imatinib against the T315I mutation was discussed.Communicated by Ramaswamy H. Sarma.

Affinities and Kinetics Detection of Protein-Small Molecule Interactions with a Monolayer MoS2 -Based Optical Imaging Platform

Quantifying the binding kinetics and affinities of protein-small molecule interactions is critical for biomarker validation, drug discovery, and deep understanding of various biological processes at the molecular-scale. Novel approaches are demanded as most common label-free techniques are mass-sensitive, which are not suitable for the detection of small molecule interactions. Here, an optical imaging platform is developed to measure the binding kinetics of both protein-small molecules and protein-ions based on monolayer MoS₂ , an ultra-thin 2D material whose optical absorption is extremely sensitive to charge. A model is established to calibrate the optical response due to the charged analyte binding and it is applied to quantify the interactions between abl1 kinase and different small-molecule inhibitors. Such a presented method is capable of distinguishing different inhibitors binding to a wild or mutated kinase, which provides guidance for drug evaluation and drug mechanism exploration. The binding kinetics of calcium ions to calmodulin is also measured, further broadening the application field of the method. In addition, the imaging capability allows mapping the local binding kinetics of the molecular interactions with a high resolution, which reveals visible spatial variability and offers a promising tool for studying heterogeneous local interfacial interactions.

Design, synthesis and biological evaluation of a series of novel pyrrolo[2,3-d]pyrimidin/pyrazolo[3,4-d]pyrimidin-4-amine derivatives as FGFRs-dominant multi-target receptor tyrosine kinase inhibitors for the treatment of gastric cancer

Gastric cancer is the second most lethal cancer across the world. With the progress in therapeutic approaches, the 5-year survival rate of early gastric cancer can reach > 95%. However, the prognosis and survival time of advanced gastric cancer is still somber. Therefore, more effective targeted therapies for gastric cancer treatment are urgently needed. FGFR, VEGFR and other receptor tyrosine kinases have recently been suggested as potential targets for gastric cancer treatment. We herein report the discovery of pyrrolo[2,3-d]pyrimidin/pyrazolo[3,4-d]pyrimidin-4-amine derivatives as a new class of FGFRs-dominant multi-target receptor tyrosine kinase inhibitors. SAR assessment identified the most active compounds 8f and 8k, which showed excellent inhibitory activity against a variety of receptor tyrosine kinases. Moreover, 8f and 8k displayed excellent potency in the SNU-16 gastric cancer cell line. Furthermore, 8f and 8k could inhibit FGFR1 phosphorylation and downstream signaling pathways as well as induce cell apoptosis. In vivo, 8f and 8k suppress tumor growth in the SNU-16 xenograft model without inducing obvious toxicity. These findings raise the possibility that compounds 8f and 8k might serve as potential agents for the treatment of gastric cancer.

Reduced Cardiotoxicity of Ponatinib-Loaded PLGA-PEG-PLGA Nanoparticles in Zebrafish Xenograft Model

Tyrosine kinase inhibitors (TKIs) are the new generation of anti-cancer drugs with high potential against cancer cells’ proliferation and growth. However, TKIs are associated with severe cardiotoxicity, limiting their clinical value. One TKI that has been developed recently but not explored much is Ponatinib. The use of nanoparticles (NPs) as a better therapeutic agent to deliver anti-cancer drugs and reduce their cardiotoxicity has been recently considered. In this study, with the aim to reduce Ponatinib cardiotoxicity, Poly(D,L-lactide-co-glycolide)-b-poly(ethyleneoxide)-b-poly(D,L-lactide-co-glycolide) (PLGA-PEG-PLGA) triblock copolymer was used to synthesize Ponatinib in loaded PLGA-PEG-PLGA NPs for chronic myeloid leukemia (CML) treatment. In addition to physicochemical NPs characterization (NPs shape, size, size distribution, surface charge, dissolution rate, drug content, and efficacy of encapsulation) the efficacy and safety of these drug-delivery systems were assessed in vivo using zebrafish. Zebrafish are a powerful animal model for investigating the cardiotoxicity associated with anti-cancer drugs such as TKIs, to determine the optimum concentration of smart NPs with the least side effects, and to generate a xenograft model of several cancer types. Therefore, the cardiotoxicity of unloaded and drug-loaded PLGA-PEG-PLGA NPs was studied using the zebrafish model by measuring the survival rate and cardiac function parameters, and therapeutic concentration for in vivo efficacy studies was optimized in an in vivo setting. Further, the efficacy of drug-loaded PLGA-PEG-PLGA NPs was tested on the zebrafish cancer xenograft model, in which human myelogenous leukemia cell line K562 was transplanted into zebrafish embryos. Our results demonstrated that the Ponatinib-loaded PLGA-PEG-PLGA NPs at a concentration of 0.001 mg/mL are non-toxic/non-cardio-toxic in the studied zebrafish xenograft model.

From Fragment to Lead: De Novo Design and Development toward a Selective FGFR2 Inhibitor

Fibroblast growth factor receptors (FGFRs) are implicated in a range of cancers with several pan-kinase and selective-FGFR inhibitors currently being evaluated in clinical trials. Pan-FGFR inhibitors often cause toxic side effects and few examples of subtype-selective inhibitors exist. Herein, we describe a structure-guided approach toward the development of a selective FGFR2 inhibitor. De novo design was carried out on an existing fragment series to yield compounds predicted to improve potency against the FGFRs. Subsequent iterative rounds of synthesis and biological evaluation led to an inhibitor with nanomolar potency that exhibited moderate selectivity for FGFR2 over FGFR1/3. Subtle changes to the lead inhibitor resulted in a complete loss of selectivity for FGFR2. X-ray crystallographic studies revealed inhibitor-specific morphological differences in the P-loop which were posited to be fundamental to the selectivity of these compounds. Additional docking studies have predicted an FGFR2-selective H-bond which could be utilized to design more selective FGFR2 inhibitors.

DNA repair glycosylase hNEIL1 triages damaged bases via competing interaction modes

DNA glycosylases must distinguish the sparse damaged sites from the vast expanse of normal DNA bases. However, our understanding of the nature of nucleobase interrogation is still limited. Here, we show that hNEIL1 (human endonuclease VIII-like 1) captures base lesions via two competing states of interaction: an activated state that commits catalysis and base excision repair, and a quarantine state that temporarily separates and protects the flipped base via auto-inhibition. The relative dominance of the two states depends on key residues of hNEIL1 and chemical properties (e.g. aromaticity and hydrophilicity) of flipped bases. Such a DNA repair mechanism allows hNEIL1 to recognize a broad spectrum of DNA damage while keeps potential gratuitous repair in check. We further reveal the molecular basis of hNEIL1 activity regulation mediated by post-transcriptional modifications and provide an example of how exquisite structural dynamics serves for orchestrated enzyme functions.

hNEIL1 (human endonuclease VIII-like 1) is a broadly specific DNA glycosylase for base excision repair. Here, the authors show that hNEIL1 can assume activated or triage conformations: the structural basis for the mechanism that enables broad specificity and reduces futile repair of normal bases.

Classical MD and metadynamics simulations on back-pocket binders of CDK2 and VEGFR2: a guidepost to design novel small-molecule dual inhibitors

Cyclin-Dependent Kinase 2 (CDK2) and Vascular-Endothelial Growth Factor Receptor 2 (VEGFR2) are promising targets for the design of novel inhibitors in anticancer therapeutics. In a recent work, our group designed a set of potential dual inhibitors predicted to occupy an allosteric back pocket near the active site of both enzymes, but their dynamic and unbinding behavior was unclear. Here, we used molecular dynamics (MD) and metadynamics (meta-D) simulations to study two of these virtual candidates (herein called IQ2 and IQ3). Their binding mode was predicted to be similar to that observed in LQ5 and BAX, well-known back-pocket binders of CDK2 and VEGFR2, respectively, including H-bonding with critical residues such as Leu83/Cys113 and Asp145/Asp190 (but excepting H-bonding with Glu51/Glu111) in CDK2/VEGFR2, correspondingly. Likewise, while LQ5 and BAX unbound through the allosteric channel as expected for type-IIA inhibitors, IQ2 and IQ3 unbound via the ATP channel (except for CDK2-IQ2) as expected for type-I½A inhibitors. Interestingly, a C-C single/double bond difference between IQ2/IQ3, respectively, resulted associated with differences in the AS/T loop flexibility observed for CDK2. These insights will help developing scaffold modifications during an optimization stage, serving as a starting point to develop dual kinase inhibitors in challenging biological targets with a promising anticancer potential.Communicated by Ramaswamy H. Sarma.

Mapping Ligand-Shape Space for Protein-Ligand Systems: Distinguishing Key-in-Lock and Hand-in-Glove Proteins

Many of the recently developed methods to study the shape of molecules permit one conformation of one molecule to be compared to another conformation of the same or a different molecule: a relative shape. Other methods provide an absolute description of the shape of a conformation that does not rely on comparisons or overlays. Any absolute description of shape can be used to generate a self-organizing map (shape map) that places all molecular shapes relative to one another; in the studies reported here, the shape fingerprint and ultrafast shape recognition methods are employed to create such maps. In the shape maps, molecules that are near one another have similar shapes, and the maps for the 102 targets in the DUD-E set have been generated. By examining the distribution of actives in comparison with their physical-property-matched decoys, we show that the proteins of key-in-lock type (relatively rigid receptor and ligand) can be distinguished from those that are more of a hand-in-glove type (more flexible receptor and ligand). These are linked to known differences in protein flexibility and binding-site size.

The application of the MM/GBSA method in the binding pose prediction of FGFR inhibitors

The success of a structure-based drug is highly dependent on a known binding pose of the protein-ligand system. However, this is not always available. In this study, we set out to explore the applicability of the popular and easy-to-use MD-based MM/GBSA method to determine the binding poses of known FGFR inhibitors. It was found that MM/GBSA combined with 100 ns of MD simulation significantly improved the success rate of docking methods from 30-40% to 70%. This work demonstrates a way that the MM/GBSA method can be more accurate than it is in ligand ranking, filling a gap in structure-based drug discovery when the binding pose is unknown.

The structure of human GCN2 reveals a parallel, back-to-back kinase dimer with a plastic DFG activation loop motif

When activated by amino acid starvation, the stress sensing protein kinase GCN2 phosphorylates the eukaryotic initiation factor 2 alpha, inhibiting translation to conserve energy and facilitate cell survival. Amino acid starvation, particularly of tryptophan and arginine, affects immune tolerance by suppressing differentiation and proliferation of T-cells via activation of GCN2 kinase. In addition, the GCN2 pathway mediates cancer survival directly within the context of metabolic stress. Here, we report the first crystal structures of the human GCN2 kinase domain (KD) in complex with two inhibitors of different size, shape, and chemical scaffold. Three novel activation loop conformations representative of different activation states of the kinase are described. In addition, a novel dimerization organization for GCN2 is observed. This arrangement is consistent with the hypothesis that the GCN2 KD forms an antiparallel inactive dimer until uncharged tRNA binds to it and triggers conformational changes that shift the equilibrium to the active parallel dimer.

Utilizing Stimulated Raman Scattering Microscopy To Study Intracellular Distribution of Label-Free Ponatinib in Live Cells

Stimulated Raman scattering (SRS) microscopy represents a powerful method for imaging label-free drug distribution with high resolution. SRS was applied to image label-free ponatinib with high sensitivity and specificity in live human chronic myeloid leukemia (CML) cell lines. This was achieved at biologically relevant, nanomolar concentrations, allowing determination of ponatinib uptake and sequestration into lysosomes during the development of acquired drug resistance and an improved understanding of target engagement.

A combined computational and experimental strategy identifies mutations conferring resistance to drugs targeting the BCR-ABL fusion protein

Drug resistance is of increasing concern, especially during the treatments of infectious diseases and cancer. To accelerate the drug discovery process in combating issues of drug resistance, here we developed a computational and experimental strategy to predict drug resistance mutations. Using BCR-ABL as a case study, we successfully recaptured the clinically observed mutations that confer resistance imatinib, nilotinib, dasatinib, bosutinib, and ponatinib. We then experimentally tested the predicted mutants in vitro. We found that although all mutants showed weakened binding strength as expected, the binding constants alone were not a good indicator of drug resistance. Instead, the half-maximal inhibitory concentration (IC50) was shown to be a good indicator of the incidence of the predicted mutations, together with change in catalytic efficacy. Our suggested strategy for predicting drug-resistance mutations includes the computational prediction and in vitro selection of mutants with increased IC50 values beyond the drug safety window.

A Dynamic Switch in Inactive p38γ Leads to an Excited State on the Pathway to an Active Kinase

The inactive state of mitogen-activated protein kinases (MAPKs) adopts an open conformation while the active state exists in a compact form stabilized by phosphorylation. In the active state, eukaryotic kinases undergo breathing motions related to substrate binding and product release that have not previously been detected in the inactive state. However, docking interactions of partner proteins with inactive MAPK kinases exhibit allostery in binding of activating kinases. Interactions at a site distant from the activation loop are coupled to the configuration of the activation loop, suggesting that the inactive state may also undergo concerted dynamics. X-ray crystallographic studies of nonphosphorylated, inactive p38γ reveal differences in domain orientations and active site structure in the two molecules in the asymmetric unit. One molecule resembles an inactive kinase with an open active site. The second molecule has a rotation of the N-lobe that leads to partial compaction of the active site, resulting in a conformation that is intermediate between the inactive open state and the fully closed state of the activated kinase. Although the compact state of apo p38γ displays several of the features of the activated enzyme, it remains catalytically inert. In solution, the kinase fluctuates on a millisecond time scale between the open ground state and a weakly populated excited state that is similar in structure to the compact state observed in the crystal. The nuclear magnetic resonance and crystal structure data imply that interconversion between the open and compact states involves a molecular switch associated with the DFG loop.

Peptide-based NTA(Ni)-nanodiscs for studying membrane enhanced FGFR1 kinase activities

Tyrosine autophosphorylation plays a crucial regulatory role in the kinase activities of fibroblast growth factor receptors (FGFRs), and in the recruitment and activation of downstream intracellular signaling pathways. Biophysical and biochemical investigations of FGFR kinase domains in membrane environments offer key insights into phosphorylation mechanisms. Hence, we constructed nickel chelating nanodiscs based on a 22-residue peptide. The spontaneous anchoring of N-terminal His₆-tagged FGFR1c kinase domain (FGFR1K) onto peptide nanodiscs grants FGFR1K orientations occurring on native plasma membranes. Following membrane incorporation, the autophosphorylation of FGFR1K, as exemplified by Y653 and Y654 in the A-loop and the total tyrosine phosphorylation, increase significantly. This in vitro reconstitution system may be applicable to studies of other membrane associated phenomena.

Application of a Substrate-Mediated Selection with c-Src Tyrosine Kinase to a DNA-Encoded Chemical Library

As aberrant activity of protein kinases is observed in many disease states, these enzymes are common targets for therapeutics and detection of activity levels. The development of non-natural protein kinase substrates offers an approach to protein substrate competitive inhibitors, a class of kinase inhibitors with promise for improved specificity. Also, kinase activity detection approaches would benefit from substrates with improved activity and specificity. Here, we apply a substrate-mediated selection to a peptidomimetic DNA-encoded chemical library for enrichment of molecules that can be phosphorylated by the protein tyrosine kinase, c-Src. Several substrates were identified and characterized for activity. A lead compound (SrcDEL10) showed both the ability to serve as a substrate and to promote ATP hydrolysis by the kinase. In inhibition assays, compounds displayed IC_50's ranging from of 8-100 µM. NMR analysis of SrcDEL10 bound to the c-Src:ATP complex was conducted to characterize the binding mode. An ester derivative of the lead compound demonstrated cellular activity with inhibition of Src-dependent signaling in cell culture. Together, the results show the potential for substrate-mediated selections of DNA-encoded libraries to discover molecules with functions other than simple protein binding and offer a new discovery method for development of synthetic tyrosine kinase substrates.

Fibroblast Growth Factor Receptors (FGFRs): Structures and Small Molecule Inhibitors

Fibroblast growth factor receptors (FGFRs) are a family of receptor tyrosine kinases expressed on the cell membrane that play crucial roles in both developmental and adult cells. Dysregulation of FGFRs has been implicated in a wide variety of cancers, such as urothelial carcinoma, hepatocellular carcinoma, ovarian cancer and lung adenocarcinoma. Due to their functional importance, FGFRs have been considered as promising drug targets for the therapy of various cancers. Multiple small molecule inhibitors targeting this family of kinases have been developed, and some of them are in clinical trials. Furthermore, the pan-FGFR inhibitor erdafitinib (JNJ-42756493) has recently been approved by the U.S. Food and Drug Administration (FDA) for the treatment of metastatic or unresectable urothelial carcinoma (mUC). This review summarizes the structure of FGFR, especially its kinase domain, and the development of small molecule FGFR inhibitors.

Detection of Drug Binding to a Target Protein Using EVV 2DIR Spectroscopy

We demonstrate that electron-vibration-vibration two-dimensional infrared spectroscopy (EVV 2DIR) can be used to detect the binding of a drug to a target protein-active site. The EVV 2DIR spectrum of the FGFR1 kinase target protein is found to have ∼200 detectable cross-peaks in the spectral region 1250-1750 cm^-1/2600-3400 cm^-1, with additional 63 peaks caused by the addition of a drug, SU5402. Of these 63 new peaks, it is shown that only six are due to protein-drug interactions, with the other 57 being due to vibrational coupling within the drug itself. Quantum mechanical calculations employing density functional theory are used to support assignment of the six binding-dependent peaks, with one being assigned to a known interaction between the drug and a backbone carbonyl group which forms part of the binding site. None of the 57 intramolecular coupling peaks associated with the drug molecule change substantially in either intensity or frequency when the drug binds to the target protein. This strongly suggests that the structure of the drug in the target binding site is essentially identical to that when it is not bound.

Structure, activation and dysregulation of fibroblast growth factor receptor kinases: perspectives for clinical targeting

The receptor tyrosine kinase family of fibroblast growth factor receptors (FGFRs) play crucial roles in embryonic development, metabolism, tissue homeostasis and wound repair via stimulation of intracellular signalling cascades. As a consequence of FGFRs' influence on cell growth, proliferation and differentiation, FGFR signalling is frequently dysregulated in a host of human cancers, variously by means of overexpression, somatic point mutations and gene fusion events. Dysregulation of FGFRs is also the underlying cause of many developmental dysplasias such as hypochondroplasia and achondroplasia. Accordingly, FGFRs are attractive pharmaceutical targets, and multiple clinical trials are in progress for the treatment of various FGFR aberrations. To effectively target dysregulated receptors, a structural and mechanistic understanding of FGFR activation and regulation is required. Here, we review some of the key research findings from the last couple of decades and summarise the strategies being explored for therapeutic intervention.

NMR backbone assignments of the tyrosine kinase domain of human fibroblast growth factor receptor 3 in apo state and in complex with inhibitor PD173074

Fibroblast growth factors receptors (FGFR) are transmembrane protein tyrosine kinases involved in many cellular process, including growth, differentiation and angiogenesis. Dysregulation of FGFR enzymatic activity is associated with developmental disorders and cancers; therefore FGFRs have become attractive targets for drug discovery, with a number of agents in late-stage clinical trials. Here, we present the backbone resonance assignments of FGFR3 tyrosine kinase domain in the ligand-free form and in complex with the canonical FGFR kinase inhibitor PD173074. Analysis of chemical shift changes upon inhibitor binding highlights a characteristic pattern of allosteric network perturbations that is of relevance for future drug discovery activities aimed at development of conformationally-selective FGFR inhibitors.

Altered conformational landscape and dimerization dependency underpins the activation of EGFR by αC-β4 loop insertion mutations

Mutational activation of epidermal growth factor receptor (EGFR) in human cancers involves both point mutations and complex mutations (insertions and deletions). In particular, short in-frame insertion mutations within a conserved αC-β4 loop in the EGFR kinase domain are frequently observed in tumor samples and patients harboring these mutations are insensitive to first-generation EGFR inhibitors. Despite the prevalence and clinical relevance of insertion mutations, the mechanisms by which these mutations regulate EGFR activity and contribute to drug sensitivity are poorly understood. Using cell-based mutation screening, we find that the precise location, length, and sequence of the inserted segment are critical for ligand-independent EGFR activation and downstream signaling. We identify three insertion mutations (N771_P772insN, D770_N771insG, and D770>GY) that activate EGFR in a unique way by relying more on the "acceptor" interface for kinase activation. Our drug inhibition studies indicate that these activating insertion mutations respond more favorably to osimertinib, a recently Food and Drug Administration-approved EGFR inhibitor for T790M-positive patients with lung cancer. Molecular dynamics simulations and umbrella sampling of WT and mutant EGFR suggest a model in which activating insertion mutations increase catalytic activity by relieving key autoinhibitory interactions associated with αC-helix movement and by lowering the transition free energy ([Formula: see text]) between active and inactive states. Our studies also identify a transition state sampled by activating insertion mutations that can be exploited in the design of mutant-selective EGFR inhibitors.

Disease Variants of FGFR3 Reveal Molecular Basis for the Recognition and Additional Roles for Cdc37 in Hsp90 Chaperone System

Receptor tyrosine kinase FGFR3 is involved in many signaling networks and is frequently mutated in developmental disorders and cancer. The Hsp90/Cdc37 chaperone system is essential for function of normal and neoplastic cells. Here we uncover the mechanistic inter-relationships between these proteins by combining approaches including NMR, HDX-MS, and SAXS. We show that several disease-linked mutations convert FGFR3 to a stronger client, where the determinant underpinning client strength involves an allosteric network through the N-lobe and at the lobe interface. We determine the architecture of the client kinase/Cdc37 complex and demonstrate, together with site-specific information, that binding of Cdc37 to unrelated kinases induces a common, extensive conformational remodeling of the kinase N-lobe, beyond localized changes and interactions within the binary complex. As further shown for FGFR3, this processing by Cdc37 deactivates the kinase and presents it, in a specific orientation established in the complex, for direct recognition by Hsp90.

Protein conformational flexibility modulates kinetics and thermodynamics of drug binding

Structure-based drug design has often been restricted by the rather static picture of protein-ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery.

An Alignment-Independent 3D-QSAR Study of FGFR2 Tyrosine Kinase Inhibitors

Purpose: Receptor tyrosine kinase (RTK) inhibitors are widely used pharmaceuticals in cancer therapy. Fibroblast growth factor receptors (FGFRs) are members of RTK superfamily which are highly expressed on the surface of carcinoma associate fibroblasts (CAFs). The involvement of FGFRs in different types of cancer makes them promising target in cancer therapy and hence, the identification of novel FGFR inhibitors is of great interest. In the current study we aimed to develop an alignment independent three dimensional quantitative structure-activity relationship (3D-QSAR) model for a set of 26 FGFR2 kinase inhibitors allowing the prediction of activity and identification of important structural features for these inhibitors.

Methods: Pentacle software was used to calculate grid independent descriptors (GRIND) for the active conformers generated by docking followed by the selection of significant variables using fractional factorial design (FFD). The partial least squares (PLS) model generated based on the remaining descriptors was assessed by internal and external validation methods.

Results: Six variables were identified as the most important probes-interacting descriptors with high impact on the biological activity of the compounds. Internal and external validations were lead to good statistical parameters (r² values of 0.93 and 0.665, respectively).

Conclusion: The results showed that the model has good predictive power and may be used for designing novel FGFR2 inhibitors.

Genetics of Pheochromocytomas and Paragangliomas: An Overview on the Recently Implicated Genes MERTK, MET, Fibroblast Growth Factor Receptor 1, and H3F3A

Genomic studies conducted by different centers have uncovered various new genes mutated in pheochromocytomas and paragangliomas (PPGLs) at germline, mosaic, and/or somatic levels, greatly expanding our knowledge of the genetic events occurring in these tumors. The current review focuses on very new findings and discusses the previously not recognized role of MERTK, MET, fibroblast growth factor receptor 1, and H3F3A genes in syndromic and nonsyndromic PPGLs. These 4 new genes were selected because although their association with PPGLs is very recent, mounting evidence was generated that rapidly consolidated the prominence of these genes in the molecular pathogenesis of PPGLs.

Kinetics for Drug Discovery: an industry-driven effort to target drug residence time

A considerable number of approved drugs show non-equilibrium binding characteristics, emphasizing the potential role of drug residence times for in vivo efficacy. Therefore, a detailed understanding of the kinetics of association and dissociation of a target-ligand complex might provide crucial insight into the molecular mechanism-of-action of a compound. This deeper understanding will help to improve decision making in drug discovery, thus leading to a better selection of interesting compounds to be profiled further. In this review, we highlight the contributions of the Kinetics for Drug Discovery (K4DD) Consortium, which targets major open questions related to binding kinetics in an industry-driven public-private partnership.

Conformational transition of FGFR kinase activation revealed by site-specific unnatural amino acid reporter and single molecule FRET

Protein kinases share significant structural similarity; however, structural features alone are insufficient to explain their diverse functions. Thus, bridging the gap between static structure and function requires a more detailed understanding of their dynamic properties. For example, kinase activation may occur via a switch-like mechanism or by shifting a dynamic equilibrium between inactive and active states. Here, we utilize a combination of FRET and molecular dynamics (MD) simulations to probe the activation mechanism of the kinase domain of Fibroblast Growth Factor Receptor (FGFR). Using genetically-encoded, site-specific incorporation of unnatural amino acids in regions essential for activation, followed by specific labeling with fluorescent moieties, we generated a novel class of FRET-based reporter to monitor conformational differences corresponding to states sampled by non phosphorylated/inactive and phosphorylated/active forms of the kinase. Single molecule FRET analysis in vitro, combined with MD simulations, shows that for FGFR kinase, there are populations of inactive and active states separated by a high free energy barrier resulting in switch-like activation. Compared to recent studies, these findings support diversity in features of kinases that impact on their activation mechanisms. The properties of these FRET-based constructs will also allow further studies of kinase dynamics as well as applications in vivo.

An overview of the binding models of FGFR tyrosine kinases in complex with small molecule inhibitors

The fibroblast growth factor receptor (FGFR) family receptor tyrosine kinase (RTK) includes four structurally related members, termed as FGFR1, FGFR2, FGFR3, and FGFR4. Given its intimate role in the progression of several solid tumors, excessive FGFR signaling provides an opportunity for anticancer therapy. Along with extensive pharmacological studies validating the therapeutic potential of targeting the FGFRs for cancer treatment, co-crystal structures of FGFRs/inhibitors are continuously coming up to study the mechanism of actions and explore new inhibitors. Herein, we review the reported co-crystals of FGFRs in complex with the corresponding inhibitors, main focusing our attention on the binding models and the pharmacological activities of the inhibitors.

Biophysical characterization of drug-resistant mutants of fibroblast growth factor receptor 1

Over-expression and aberrant activation of tyrosine kinases occur frequently in human cancers. Various tyrosine kinase inhibitors (TKIs) are under clinical use, but acquisition of resistance to these drugs is a major problem. Here, we studied the interaction between two drug-resistant mutants of fibroblast growth factor receptor 1 (FGFR1), N546K and V561M, and four ATP-competitive inhibitors, ponatinib, dovitinib, PD173074 and BGJ-398. Among these protein-drug systems, the only marked reduction in affinity was that of PD173074 for the V561M mutant. We also examined the interaction of these FGFR1 variants to AMP-PNP, a nonhydrolyzable analogue of ATP, and showed that N546K showed increased affinity for the ATP analogue as compared with the wild type. These findings will help to clarify the mechanism of drug resistance in mutant tyrosine kinases.

Landscape of activating cancer mutations in FGFR kinases and their differential responses to inhibitors in clinical use

Frequent genetic alterations discovered in FGFRs and evidence implicating some as drivers in diverse tumors has been accompanied by rapid progress in targeting FGFRs for anticancer treatments. Wider assessment of the impact of genetic changes on the activation state and drug responses is needed to better link the genomic data and treatment options. We here apply a direct comparative and comprehensive analysis of FGFR3 kinase domain variants representing the diversity of point-mutations reported in this domain. We reinforce the importance of N540K and K650E and establish that not all highly activating mutations (for example R669G) occur at high-frequency and conversely, that some "hotspots" may not be linked to activation. Further structural characterization consolidates a mechanistic view of FGFR kinase activation and extends insights into drug binding. Importantly, using several inhibitors of particular clinical interest (AZD4547, BGJ-398, TKI258, JNJ42756493 and AP24534), we find that some activating mutations (including different replacements of the same residue) result in distinct changes in their efficacy. Considering that there is no approved inhibitor for anticancer treatments based on FGFR-targeting, this information will be immediately translatable to ongoing clinical trials.

Conformational Selection and Induced Fit Mechanisms in the Binding of an Anticancer Drug to the c-Src Kinase

Understanding the conformational changes associated with the binding of small ligands to their biological targets is a fascinating and meaningful question in chemistry, biology and drug discovery. One of the most studied and important is the so-called "DFG-flip" of tyrosine kinases. The conserved three amino-acid DFG motif undergoes an "in to out" movement resulting in a particular inactive conformation to which "type II" kinase inhibitors, such as the anti-cancer drug Imatinib, bind. Despite many studies, the details of this prototypical conformational change are still debated. Here we combine various NMR experiments and surface plasmon resonance with enhanced sampling molecular dynamics simulations to shed light into the conformational dynamics associated with the binding of Imatinib to the proto-oncogene c-Src. We find that both conformational selection and induced fit play a role in the binding mechanism, reconciling opposing views held in the literature. Moreover, an external binding pose and local unfolding (cracking) of the aG helix are observed.