Molecular dynamics simulations show that conformational selection governs the binding preferences of imatinib for several tyrosine kinases

New ABL1 Kinase Domain Mutations in BCR::ABL1-Positive Acute Lymphoblastic Leukemia

BACKGROUND

Since the development of the first-generation Tyrosine Kinase Inhibitor (TKI), it has played a crucial role in the treatment of BCR::ABL1-positive acute lymphoblastic leukemia (ALL) and chronic myeloid leukemia (CML). However, ABL1 kinase domain (ABL1 KD) mutations confer resistance to several TKIs. These mutations have been extensively studied in chronic myeloid leukemia (CML) but less so in BCR::ABL1-positive acute lymphoblastic leukemia (ALL).

METHODS

Our study aimed to analyze the the ABL1 KD mutations in 97 consecutive newly-diagnosed adults with BCR::ABL1-positive ALL before therapy, in cytogenetic complete remission and at relapse with next generation sequencing (NGS). The relationship between ABL1 KD mutations and TKI selection was also analyzed.

RESULTS

Previously unreported ABL1 KD mutations R239G, F401V/L, R516L and K262T were the most prevalent in pre-therapy and cytogenetic remission samples, whereas T315I/P and P-loop mutations were most prevalent in relapse samples. R239G, F401V/L, R516L and K262T are related to the BCR::ABL1 structure, whereas T315I/P and P-loop mutations directly alter kinase activity. BaF3 cells transfected with ABL1 KD F401V, K262T, R239G, or R516L mutations were resistant to imatinib but strongly inhibited by olverembatinib with IC50 values of 0.73 to 1.52nM. Meanwhile, olverembatinib had advantages in increasing complete molecular response (CMR) and good prognosis.

CONCLUSION

Overall, our findings indicate the prevalence and impact of new ABL1 KD mutations in BCR::ABL1-positive ALL patients, highlighting the necessity for effective therapies targetingthese mutations.

Molecular basis for differential recognition of an allosteric inhibitor by receptor tyrosine kinases

Understanding kinase-inhibitor selectivity continues to be a major objective in kinase drug discovery. We probe the molecular basis of selectivity of an allosteric inhibitor (MSC1609119A-1) of the insulin-like growth factor-I receptor kinase (IGF1RK), which has been shown to be ineffective for the homologous insulin receptor kinase (IRK). Specifically, we investigated the structural and energetic basis of the allosteric binding of this inhibitor to each kinase by combining molecular modeling, molecular dynamics (MD) simulations, and thermodynamic calculations. We predict the inhibitor conformation in the binding pocket of IRK and highlight that the charged residues in the histidine-arginine-aspartic acid (HRD) and aspartic acid-phenylalanine-glycine (DFG) motifs and the nonpolar residues in the binding pocket govern inhibitor interactions in the allosteric pocket of each kinase. We suggest that the conformational changes in the IGF1RK residues M1054 and M1079, movement of the ⍺C-helix, and the conformational stabilization of the DFG motif favor the selectivity of the inhibitor toward IGF1RK. Our thermodynamic calculations reveal that the observed selectivity can be rationalized through differences observed in the electrostatic interaction energy of the inhibitor in each inhibitor/kinase complex and the hydrogen bonding interactions of the inhibitor with the residue V1063 in IGF1RK that are not attained with the corresponding residue V1060 in IRK. Overall, our study provides a rationale for the molecular basis of recognition of this allosteric inhibitor by IGF1RK and IRK, which is potentially useful in developing novel inhibitors with improved affinity and selectivity.

Kinematic analysis of kinases and their oncogenic mutations - Kinases and their mutation kinematic analysis

Protein kinases are crucial cellular enzymes that facilitate the transfer of phosphates from adenosine triphosphate (ATP) to their substrates, thereby regulating numerous cellular activities. Dysfunctional kinase activity often leads to oncogenic conditions. Chosen by using structural similarity to 5UG9, we selected 79 crystal structures from the PDB and based on the position of the phenylalanine side chain in the DFG motif, we classified these 79 crystal structures into 5 group clusters. Our approach applies our kinematic flexibility analysis (KFA) to explore the flexibility of kinases in various activity states and examine the impact of the activation loop on kinase structure. KFA enables the rapid decomposition of macromolecules into different flexibility regions, allowing comprehensive analysis of conformational structures. The results reveal that the activation loop of kinases acts as a "lock" that stabilizes the active conformation of kinases by rigidifying the adjacent α-helices. Furthermore, we investigate specific kinase mutations, such as the L858R mutation commonly associated with non-small cell lung cancer, which induces increased flexibility in active-state kinases. In addition, through analyzing the hydrogen bond pattern, we examine the substructure of kinases in different states. Notably, active-state kinases exhibit a higher occurrence of α-helices compared to inactive-state kinases. This study contributes to the understanding of biomolecular conformation at a level relevant to drug development.

Classifying protein kinase conformations with machine learning

Protein kinases are key actors of signaling networks and important drug targets. They cycle between active and inactive conformations, distinguished by a few elements within the catalytic domain. One is the activation loop, whose conserved DFG motif can occupy DFG-in, DFG-out, and some rarer conformations. Annotation and classification of the structural kinome are important, as different conformations can be targeted by different inhibitors and activators. Valuable resources exist; however, large-scale applications will benefit from increased automation and interpretability of structural annotation. Interpretable machine learning models are described for this purpose, based on ensembles of decision trees. To train them, a set of catalytic domain sequences and structures was collected, somewhat larger and more diverse than existing resources. The structures were clustered based on the DFG conformation and manually annotated. They were then used as training input. Two main models were constructed, which distinguished active/inactive and in/out/other DFG conformations. They considered initially 1692 structural variables, spanning the whole catalytic domain, then identified ("learned") a small subset that sufficed for accurate classification. The first model correctly labeled all but 3 of 3289 structures as active or inactive, while the second assigned the correct DFG label to all but 17 of 8826 structures. The most potent classifying variables were all related to well-known structural elements in or near the activation loop and their ranking gives insights into the conformational preferences. The models were used to automatically annotate 3850 kinase structures predicted recently with the Alphafold2 tool, showing that Alphafold2 reproduced the active/inactive but not the DFG-in proportions seen in the Protein Data Bank. We expect the models will be useful for understanding and engineering kinases.

Development of Kinase-Centric Drugs: A Computational Perspective

Kinases are prominent drug targets in the pharmaceutical and research community due to their involvement in signal transduction, physiological responses, and upon dysregulation, in diseases such as cancer, neurological and autoimmune disorders. Several FDA-approved small-molecule drugs have been developed to combat human diseases since Gleevec was approved for the treatment of chronic myelogenous leukemia. Kinases were considered "undruggable" in the beginning. Several FDA-approved small-molecule drugs have become available in recent years. Most of these drugs target ATP-binding sites, but a few target allosteric sites. Among kinases that belong to the same family, the catalytic domain shows high structural and sequence conservation. Inhibitors of ATP-binding sites can cause off-target binding. Because members of the same family have similar sequences and structural patterns, often complex relationships between kinases and inhibitors are observed. To design and develop drugs with desired selectivity, it is essential to understand the target selectivity for kinase inhibitors. To create new inhibitors with the desired selectivity, several experimental methods have been designed to profile the kinase selectivity of small molecules. Experimental approaches are often expensive, laborious, time-consuming, and limited by the available kinases. Researchers have used computational methodologies to address these limitations in the design and development of effective therapeutics. Many computational methods have been developed over the last few decades, either to complement experimental findings or to forecast kinase inhibitor activity and selectivity. The purpose of this review is to provide insight into recent advances in theoretical/computational approaches for the design of new kinase inhibitors with the desired selectivity and optimization of existing inhibitors.

TYROSINE KINASES: COMPLEX MOLECULAR SYSTEMS CHALLENGING COMPUTATIONAL METHODOLOGIES

Classical molecular dynamics (MD) simulations based on atomic models play an increasingly important role in a wide range of applications in physics, biology, and chemistry. Nonetheless, generating genuine knowledge about biological systems using MD simulations remains challenging. Protein tyrosine kinases are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Due to the large conformational changes and long timescales involved in their function, these kinases present particularly challenging problems to modern computational and theoretical frameworks aimed at elucidating the dynamics of complex biomolecular systems. Markov state models have achieved limited success in tackling the broader conformational ensemble and biased methods are often employed to examine specific long timescale events. Recent advances in machine learning continue to push the limitations of current methodologies and provide notable improvements when integrated with the existing frameworks. A broad perspective is drawn from a critical review of recent studies.

KinaseMD: kinase mutations and drug response database

Mutations in kinases are abundant and critical to study signaling pathways and regulatory roles in human disease, especially in cancer. Somatic mutations in kinase genes can affect drug treatment, both sensitivity and resistance, to clinically used kinase inhibitors. Here, we present a newly constructed database, KinaseMD (kinase mutations and drug response), to structurally and functionally annotate kinase mutations. KinaseMD integrates 679 374 somatic mutations, 251 522 network-rewiring events, and 390 460 drug response records curated from various sources for 547 kinases. We uniquely annotate the mutations and kinase inhibitor response in four types of protein substructures (gatekeeper, A-loop, G-loop and αC-helix) that are linked to kinase inhibitor resistance in literature. In addition, we annotate functional mutations that may rewire kinase regulatory network and report four phosphorylation signals (gain, loss, up-regulation and down-regulation). Overall, KinaseMD provides the most updated information on mutations, unique annotations of drug response especially drug resistance and functional sites of kinases. KinaseMD is accessible at https://bioinfo.uth.edu/kmd/, having functions for searching, browsing and downloading data. To our knowledge, there has been no systematic annotation of these structural mutations linking to kinase inhibitor response. In summary, KinaseMD is a centralized database for kinase mutations and drug response.

A shift of dynamic equilibrium between the KIT active and inactive states causes drug resistance

Tyrosine phosphorylation, a highly regulated post-translational modification, is carried out by the enzyme tyrosine kinase (TK). TKs are important mediators in signaling cascades, facilitating diverse biological processes in response to stimuli. TKs may acquire mutations leading to malignancy and are viable targets for anti-cancer drugs. Mast/stem cell growth factor receptor KIT is a TK involved in cell differentiation, whose dysregulation leads to various types of cancer, including gastrointestinal stromal tumors, leukemia, and melanoma. KIT can be targeted by a range of inhibitors that predominantly bind to the inactive state of the enzyme. A mutation Y823D in the activation loop of KIT is known to be responsible for the loss of sensitivity to some drugs in metastatic tumors. We used all-atom molecular dynamics simulations to study the impact of Y823D on the KIT conformation and dynamics and compared it to the effect of phosphorylation of Y823. We simulated in total 6.4 μs of wild-type, mutant and phosphorylated KIT in the active- and inactive-state conformations. We found that Y823D affects the protein dynamics differently: in the active state, the mutation increases the protein stability, whereas in the inactive state it induces local destabilization, thus shifting the dynamic equilibrium towards the active state, altering the communication between distant regulatory regions. The observed dynamics of the Y823D mutant is similar to the dynamics of KIT phosphorylated at position Y823, thus we hypothesize that this mutation mimics a constitutively active kinase, which is not responsive to inhibitors that bind its inactive conformation.

Landscape of drug-resistance mutations in kinase regulatory hotspots

More than 48 kinase inhibitors (KIs) have been approved by Food and Drug Administration. However, drug-resistance (DR) eventually occurs, and secondary mutations have been found in the previously targeted primary-mutated cancer cells. Cancer and drug research communities recognize the importance of the kinase domain (KD) mutations for kinasopathies. So far, a systematic investigation of kinase mutations on DR hotspots has not been done yet. In this study, we systematically investigated four types of representative mutation hotspots (gatekeeper, G-loop, αC-helix and A-loop) associated with DR in 538 human protein kinases using large-scale cancer data sets (TCGA, ICGC, COSMIC and GDSC). Our results revealed 358 kinases harboring 3318 mutations that covered 702 drug resistance hotspot residues. Among them, 197 kinases had multiple genetic variants on each residue. We further computationally assessed and validated the epidermal growth factor receptor mutations on protein structure and drug-binding efficacy. This is the first study to provide a landscape view of DR-associated mutation hotspots in kinase's secondary structures, and its knowledge will help the development of effective next-generation KIs for better precision medicine.

Cryptic pocket formation underlies allosteric modulator selectivity at muscarinic GPCRs

Allosteric modulators are highly desirable as drugs, particularly for G-protein-coupled receptor (GPCR) targets, because allosteric drugs can achieve selectivity between closely related receptors. The mechanisms by which allosteric modulators achieve selectivity remain elusive, however, particularly given recent structures that reveal similar allosteric binding sites across receptors. Here we show that positive allosteric modulators (PAMs) of the M1 muscarinic acetylcholine receptor (mAChR) achieve exquisite selectivity by occupying a dynamic pocket absent in existing crystal structures. This cryptic pocket forms far more frequently in molecular dynamics simulations of the M1 mAChR than in those of other mAChRs. These observations reconcile mutagenesis data that previously appeared contradictory. Further mutagenesis experiments validate our prediction that preventing cryptic pocket opening decreases the affinity of M1-selective PAMs. Our findings suggest opportunities for the design of subtype-specific drugs exploiting cryptic pockets that open in certain receptors but not in other receptors with nearly identical static structures.

General mathematical formula for near equilibrium relaxation kinetics of basic enzyme reactions and its applications to find conformational selection steps

A general mathematical formula of basic enzyme reactions was derived with nearly no dependence on conditions nor assumptions on relaxation kinetic processes near equilibrium in a simple single-substrate-single-product enzyme reaction. The new formula gives precise relationships between the rate constants of the elementary reaction steps and the apparent relaxation rate constant, rather than the initial velocity that is generally used to determine enzymatic parameters according to the Michaelis-Menten theory. The present formula is shown to be complementary to the Michaelis-Menten formulae in a sense that the initial velocity and the relaxation rate constant data together could determine the enzyme-substrate dissociation constant K_ES, which has been usually conditionally approximated by the Michaelis constant K_M within the framework of the Michaelis-Menten formulae. We also describe relaxation kinetics of enzyme reactions that include the conformational selection processes, in which only one enzymatic conformer among a conformational ensemble can bind with either the substrate or product. The present mathematical approaches, together with numerical computation analyses, suggested that the presence of conformational selection steps in enzymatic reactions can be experimentally detected simply by enzymatic assays with catalytic amounts of enzyme.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Advances and challenges in drug design against tuberculosis: application of in silico approaches

INTRODUCTION

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) remains the deadliest infectious disease in the world with one-third of the world's population thought to be infected. Over the years, TB mortality rate has been largely reduced; however, this progress has been threatened by the increasing appearance of multidrug-resistant Mtb. Considerable recent efforts have been undertaken to develop new generation antituberculosis drugs. Many of these attempts have relied on in silico approaches, which have emerged recently as powerful tools complementary to biochemical attempts. Areas covered: The authors review the status of pharmaceutical drug development against TB with a special emphasis on computational work. They focus on those studies that have been validated by in vitro and/or in vivo experiments, and thus, that can be considered as successful. The major goals of this review are to present target protein systems, to highlight how in silico efforts compliment experiments, and to aid future drug design endeavors. Expert opinion: Despite having access to all of the gene and protein sequences of Mtb, the search for new optimal treatments against this deadly pathogen are still ongoing. Together with the geometric growth of protein structural and sequence databases, computational methods have become a powerful technique accelerating the successful identification of new ligands.

Molecular dynamics simulations provide insights into the origin of gleevec's selectivity toward human tyrosine kinases

Protein kinases are critical drug targets against cancer. Since the discovery of Gleevec, a specific inhibitor of Abl kinase, the capability of this drug to distinguish between Abl and other tyrosine kinases, such as Src, has been intensely investigated but the origin of Gleevec's selectivity to Abl against Src is less studied. Here, we performed molecular dynamics (MD) simulations, dynamical cross-correlation matrices (DCCM), dynamical network analysis, and binding free energy calculations to explore Gleevec's selectivity based on the crystal structures of Abl, Src, and their common ancestors (ANC-AS) and the two constructed mutation systems (AS→Abl and AS→Src). MD simulations revealed that the conformation of the phosphate-binding loop (P-loop) was altered significantly in the AS→Abl system. DCCM results unraveled that mutations increased anticorrelated motions in the AS→Abl system. Community network analysis suggested that the P-loop established special contacts in the AS→Abl system that are devoid in the AS→Src system. The binding free energy calculations unveiled that the affinity of Gleevec to AS→Abl increased to near the Abl level, whereas its affinity to AS→Src decreased to near the Src level. Analysis of individual residue contributions showed that the differences were located mainly at the P-loop. This study is valuable for understanding the sensitivity of Gleevec to human tyrosine kinases. Communicated by Ramaswamy H. Sarma.

Discovery of N-(4-{[5-Fluoro-7-(2-methoxyethoxy)quinazolin-4-yl]amino}phenyl)-2-[4-(propan-2-yl)-1 H-1,2,3-triazol-1-yl]acetamide (AZD3229), a Potent Pan-KIT Mutant Inhibitor for the Treatment of Gastrointestinal Stromal Tumors

While the treatment of gastrointestinal stromal tumors (GISTs) has been revolutionized by the application of targeted tyrosine kinase inhibitors capable of inhibiting KIT-driven proliferation, diverse mutations to this kinase drive resistance to established therapies. Here we describe the identification of potent pan-KIT mutant kinase inhibitors that can be dosed without being limited by the tolerability issues seen with multitargeted agents. This effort focused on identification and optimization of an existing kinase scaffold through the use of structure-based design. Starting from a series of previously reported phenoxyquinazoline and quinoline based inhibitors of the tyrosine kinase PDGFRα, potency against a diverse panel of mutant KIT driven Ba/F3 cell lines was optimized, with a particular focus on reducing activity against a KDR driven cell model in order to limit the potential for hypertension commonly seen in second and third line GIST therapies. AZD3229 demonstrates potent single digit nM growth inhibition across a broad cell panel, with good margin to KDR-driven effects. Selectivity over KDR can be rationalized predominantly by the interaction of water molecules with the protein and ligand in the active site, and its kinome selectivity is similar to the best of the approved GIST agents. This compound demonstrates excellent cross-species pharmacokinetics, shows strong pharmacodynamic inhibition of target, and is active in several in vivo models of GIST.

Using molecular simulation to explore the nanoscale dynamics of the plant kinome

Eukaryotic protein kinases (PKs) are a large family of proteins critical for cellular response to external signals, acting as molecular switches. PKs propagate biochemical signals by catalyzing phosphorylation of other proteins, including other PKs, which can undergo conformational changes upon phosphorylation and catalyze further phosphorylations. Although PKs have been studied thoroughly across the domains of life, the structures of these proteins are sparsely understood in numerous groups of organisms, including plants. In addition to efforts towards determining crystal structures of PKs, research on human PKs has incorporated molecular dynamics (MD) simulations to study the conformational dynamics underlying the switching of PK function. This approach of experimental structural biology coupled with computational biophysics has led to improved understanding of how PKs become catalytically active and why mutations cause pathological PK behavior, at spatial and temporal resolutions inaccessible to current experimental methods alone. In this review, we argue for the value of applying MD simulation to plant PKs. We review the basics of MD simulation methodology, the successes achieved through MD simulation in animal PKs, and current work on plant PKs using MD simulation. We conclude with a discussion of the future of MD simulations and plant PKs, arguing for the importance of molecular simulation in the future of plant PK research.

A change in structural integrity of c-Kit mutant D816V causes constitutive signaling

Several signaling pathways, ligands, and genes that regulate proliferative and self-renewal properties of the Hematopoietic Stem Cells (HSCs) have been studied meticulously. One of the signaling pathways that play a crucial role in the process of hematopoiesis is the Stem Cell Factor (SCF) mediated c-Kit pathway. The c-Kit is a Receptor Tyrosine Kinase (RTK), which is expressed in the cells including HSCs. It undergoes dimerization upon binding with its cognate ligand SCF. As a result, phosphorylation of the Juxtamembrane (JM) domain of c-Kit takes place at Tyr568 and Tyr570 residues. These phosphorylated residues become the docking sites for protein tyrosine phosphatases (PTPs) namely SHP-1 and SHP-2, which in turn cause dephosphorylation and negative regulation of the downstream signaling responsible for the cell proliferation. Interestingly, it has been reported that the mutation of c-Kit at D816V makes it independent of SCF stimulation and SHP-1/SHP-2 inhibition, thereby, causing its constitutive activation. The present study was commenced to elucidate the structural behavior of this mutation in the JM and A-loop region of c-Kit using Molecular Dynamics (MD) simulations of the wild-type and mutant c-Kit in unphosphorylated and phosphorylated states. The energy difference computed between the wild type and mutant (D816V) c-Kit, and protein-protein docking and complex analysis revealed the impact of this single residue mutation on the integrity dynamics of c-Kit that makes it independent of SHP-1/SHP-2 negative regulation.

Dynamic Equilibrium of the Aurora A Kinase Activation Loop Revealed by Single-Molecule Spectroscopy

The conformation of the activation loop (T-loop) of protein kinases underlies enzymatic activity and influences the binding of small-molecule inhibitors. By using single-molecule fluorescence spectroscopy, we have determined that phosphorylated Aurora A kinase is in dynamic equilibrium between a DFG-in-like active T-loop conformation and a DFG-out-like inactive conformation, and have measured the rate constants of interconversion. Addition of the Aurora A activating protein TPX2 shifts the equilibrium towards an active T-loop conformation whereas addition of the inhibitors MLN8054 and CD532 favors an inactive T-loop. We show that Aurora A binds TPX2 and MLN8054 simultaneously and provide a new model for kinase conformational behavior. Our approach will enable conformation-specific effects to be integrated into inhibitor discovery across the kinome, and we outline some immediate consequences for structure-based drug discovery.

Structural and Dynamic Insights of the Interaction between Tritrpticin and Micelles: An NMR Study

A large number of antimicrobial peptides (AMPs) acts with high selectivity and specificity through interactions with membrane lipid components. These peptides undergo complex conformational changes in solution; upon binding to an interface, one major conformation is stabilized. Here we describe a study of the interaction between tritrpticin (TRP3), a cathelicidin AMP, and micelles of different chemical composition. The peptide's structure and dynamics were examined using one-dimensional and two-dimensional NMR. Our data showed that the interaction occurred by conformational selection and the peptide acquired similar structures in all systems studied, despite differences in detergent headgroup charge or dipole orientation. Fluorescence and paramagnetic relaxation enhancement experiments showed that the peptide is located in the interface region and is slightly more deeply inserted in 1-myristoyl-2-hydroxy-sn-glycero-3-phospho-1'-rac-glycerol (LMPG, anionic) than in 1-lauroyl-2-hydroxy-sn-glycero-3-phosphocholine (LLPC, zwitterionic) micelles. Moreover, the tilt angle of an assumed helical portion of the peptide is similar in both systems. In previous work we proposed that TRP3 acts by a toroidal pore mechanism. In view of the high hydrophobic core exposure, hydration, and curvature presented by micelles, the conformation of TRP3 in these systems could be related to the peptide's conformation in the toroidal pore.

Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases

Protein kinases are enzymes that catalyze the covalent transfer of the γ-phosphate of an adenosine triphosphate (ATP) molecule onto a tyrosine, serine, threonine, or histidine residue in the substrate and thus send a chemical signal to networks of downstream proteins. They are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Unregulated protein kinase activity is often associated with a wide range of diseases, therefore making protein kinases major therapeutic targets. A prototypical system of central interest to understand the regulation of kinase activity is provided by tyrosine kinase c-Src, which belongs to the family of Src-related non-receptor tyrosine kinases (SFKs). Although the broad picture of autoinhibition via the regulatory domains and via the phosphorylation of the C-terminal tail is well characterized from a structural point of view, a detailed mechanistic understanding at the atomic-level is lacking. Advanced computational methods based on all-atom molecular dynamics (MD) simulations are employed to advance our understanding of tyrosine kinase activation. The computational studies suggest that the isolated kinase domain (KD) is energetically most favorable in the inactive conformation when the activation loop (A-loop) of the KD is not phosphorylated. The KD makes transient visits to a catalytically competent active-like conformation. The process of bimolecular trans-autophosphorylation of the A-loop eventually locks the KD in the active state. Activating point mutations may act by slightly increasing the population of the active-like conformation, enhancing the availability of the A-loop to be phosphorylated. The Src-homology 2 (SH2) and Src-homology 3 (SH3) regulatory domains, depending upon their configuration, either promote the inactive or the active state of the kinase domain. In addition to the roles played by the SH3, SH2, and KD, the Src-homology 4-Unique domain (SH4-U) region also serves as a key moderator of substrate specificity and kinase function. Thus, a fundamental understanding of the conformational propensity of the SH4-U region and how this affects the association to the membrane surface are likely to lead to the discovery of new intermediate states and alternate strategies for inhibition of kinase activity for drug discovery. The existence of a multitude of KD conformations poses a great challenge aimed at the design of specific inhibitors. One promising computational strategy to explore the conformational flexibility of the KD is to construct Markov state models from aggregated MD data.

Insight on Mutation-Induced Resistance from Molecular Dynamics Simulations of the Native and Mutated CSF-1R and KIT

The receptors tyrosine kinases (RTKs) for the colony stimulating factor-1, CSF-1R, and for the stem cell factor, SCFR or KIT, are important mediators of signal transduction. The abnormal function of these receptors, promoted by gain-of-function mutations, leads to their constitutive activation, associated with cancer or other proliferative diseases. A secondary effect of the mutations is the alteration of receptors' sensitivity to tyrosine kinase inhibitors, compromising effectiveness of these molecules in clinical treatment. In particular, the mutation V560G in KIT increases its sensitivity to Imatinib, while the D816V in KIT, and D802V in CSF-1R, triggers resistance to the drug. We analyzed the Imatinib binding affinity to the native and mutated KIT (mutations V560G, S628N and D816V) and CSF-1R (mutation D802V) by using molecular dynamics simulations and energy calculations of Imatinib•target complexes. Further, we evaluated the sensitivity of the studied KIT receptors to Imatinib by measuring the inhibition of KIT phosphorylation. Our study showed that (i) the binding free energy of Imatinib to the targets is highly correlated with their experimentally measured sensitivity; (ii) the electrostatic interactions are a decisive factor affecting the binding energy; (iii) the most deleterious impact to the Imatinib sensitivity is promoted by D802V (CSF-1R) and D816V (KIT) mutations; (iv) the role of the juxtamembrane region, JMR, in the imatinib binding is accessory. These findings contribute to a better description of the mutation-induced effects alternating the targets sensitivity to Imatinib.

Structural propensities of kinase family proteins from a Potts model of residue co-variation

Understanding the conformational propensities of proteins is key to solving many problems in structural biology and biophysics. The co-variation of pairs of mutations contained in multiple sequence alignments of protein families can be used to build a Potts Hamiltonian model of the sequence patterns which accurately predicts structural contacts. This observation paves the way to develop deeper connections between evolutionary fitness landscapes of entire protein families and the corresponding free energy landscapes which determine the conformational propensities of individual proteins. Using statistical energies determined from the Potts model and an alignment of 2896 PDB structures, we predict the propensity for particular kinase family proteins to assume a "DFG-out" conformation implicated in the susceptibility of some kinases to type-II inhibitors, and validate the predictions by comparison with the observed structural propensities of the corresponding proteins and experimental binding affinity data. We decompose the statistical energies to investigate which interactions contribute the most to the conformational preference for particular sequences and the corresponding proteins. We find that interactions involving the activation loop and the C-helix and HRD motif are primarily responsible for stabilizing the DFG-in state. This work illustrates how structural free energy landscapes and fitness landscapes of proteins can be used in an integrated way, and in the context of kinase family proteins, can potentially impact therapeutic design strategies.

Conformation-Selective Analogues of Dasatinib Reveal Insight into Kinase Inhibitor Binding and Selectivity

In the kinase field, there are many widely held tenets about conformation-selective inhibitors that have yet to be validated using controlled experiments. We have designed, synthesized, and characterized a series of kinase inhibitor analogues of dasatinib, an FDA-approved kinase inhibitor that binds the active conformation. This inhibitor series includes two Type II inhibitors that bind the DFG-out inactive conformation and two inhibitors that bind the αC-helix-out inactive conformation. Using this series of compounds, we analyze the impact that conformation-selective inhibitors have on target binding and kinome-wide selectivity.

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.

Advanced molecular dynamics simulation methods for kinase drug discovery

Interest in the application of molecular dynamics (MD) simulations has increased in the field of protein kinase (PK) drug discovery. PKs belong to an important drug target class because they are directly involved in a number of diseases, including cancer. MD methods simulate dynamic biological and chemical events at an atomic level. This information can be combined with other in silico and experimental methods to efficiently target selected receptors. In this review, we present common and advanced methods of MD simulations and we focus on the recent applications of MD-based methodologies that provided significant insights into the elucidation of biological mechanisms involving PKs and into the discovery of novel kinase inhibitors.

Towards a Molecular Understanding of the Link between Imatinib Resistance and Kinase Conformational Dynamics

Due to its inhibition of the Abl kinase domain in the BCR-ABL fusion protein, imatinib is strikingly effective in the initial stage of chronic myeloid leukemia with more than 90% of the patients showing complete remission. However, as in the case of most targeted anti-cancer therapies, the emergence of drug resistance is a serious concern. Several drug-resistant mutations affecting the catalytic domain of Abl and other tyrosine kinases are now known. But, despite their importance and the adverse effect that they have on the prognosis of the cancer patients harboring them, the molecular mechanism of these mutations is still debated. Here by using long molecular dynamics simulations and large-scale free energy calculations complemented by in vitro mutagenesis and microcalorimetry experiments, we model the effect of several widespread drug-resistant mutations of Abl. By comparing the conformational free energy landscape of the mutants with those of the wild-type tyrosine kinases we clarify their mode of action. It involves significant and complex changes in the inactive-to-active dynamics and entropy/enthalpy balance of two functional elements: the activation-loop and the conserved DFG motif. What is more the T315I gatekeeper mutant has a significant impact on the binding mechanism itself and on the binding kinetics.

Molecular Dynamics: New Frontier in Personalized Medicine

The field of drug discovery has witnessed infinite development over the last decade with the demand for discovery of novel efficient lead compounds. Although the development of novel compounds in this field has seen large failure, a breakthrough in this area might be the establishment of personalized medicine. The trend of personalized medicine has shown stupendous growth being a hot topic after the successful completion of Human Genome Project and 1000 genomes pilot project. Genomic variant such as SNPs play a vital role with respect to inter individual's disease susceptibility and drug response. Hence, identification of such genetic variants has to be performed before administration of a drug. This process requires high-end techniques to understand the complexity of the molecules which might bring an insight to understand the compounds at their molecular level. To sustenance this, field of bioinformatics plays a crucial role in revealing the molecular mechanism of the mutation and thereby designing a drug for an individual in fast and affordable manner. High-end computational methods, such as molecular dynamics (MD) simulation has proved to be a constitutive approach to detecting the minor changes associated with an SNP for better understanding of the structural and functional relationship. The parameters used in molecular dynamic simulation elucidate different properties of a macromolecule, such as protein stability and flexibility. MD along with docking analysis can reveal the synergetic effect of an SNP in protein-ligand interaction and provides a foundation for designing a particular drug molecule for an individual. This compelling application of computational power and the advent of other technologies have paved a promising way toward personalized medicine. In this in-depth review, we tried to highlight the different wings of MD toward personalized medicine.

Molecular Determinants Underlying Binding Specificities of the ABL Kinase Inhibitors: Combining Alanine Scanning of Binding Hot Spots with Network Analysis of Residue Interactions and Coevolution

Quantifying binding specificity and drug resistance of protein kinase inhibitors is of fundamental importance and remains highly challenging due to complex interplay of structural and thermodynamic factors. In this work, molecular simulations and computational alanine scanning are combined with the network-based approaches to characterize molecular determinants underlying binding specificities of the ABL kinase inhibitors. The proposed theoretical framework unveiled a relationship between ligand binding and inhibitor-mediated changes in the residue interaction networks. By using topological parameters, we have described the organization of the residue interaction networks and networks of coevolving residues in the ABL kinase structures. This analysis has shown that functionally critical regulatory residues can simultaneously embody strong coevolutionary signal and high network centrality with a propensity to be energetic hot spots for drug binding. We have found that selective (Nilotinib) and promiscuous (Bosutinib, Dasatinib) kinase inhibitors can use their energetic hot spots to differentially modulate stability of the residue interaction networks, thus inhibiting or promoting conformational equilibrium between inactive and active states. According to our results, Nilotinib binding may induce a significant network-bridging effect and enhance centrality of the hot spot residues that stabilize structural environment favored by the specific kinase form. In contrast, Bosutinib and Dasatinib can incur modest changes in the residue interaction network in which ligand binding is primarily coupled only with the identity of the gate-keeper residue. These factors may promote structural adaptability of the active kinase states in binding with these promiscuous inhibitors. Our results have related ligand-induced changes in the residue interaction networks with drug resistance effects, showing that network robustness may be compromised by targeted mutations of key mediating residues. This study has outlined mechanisms by which inhibitor binding could modulate resilience and efficiency of allosteric interactions in the kinase structures, while preserving structural topology required for catalytic activity and regulation.

Evolution and intelligent design in drug development

Sophisticated protein kinase networks, empowering complexity in higher organisms, are also drivers of devastating diseases such as cancer. Accordingly, these enzymes have become major drug targets of the twenty-first century. However, the holy grail of designing specific kinase inhibitors aimed at specific cancers has not been found. Can new approaches in cancer drug design help win the battle with this multi-faced and quickly evolving enemy? In this perspective we discuss new strategies and ideas that were born out of a recent breakthrough in understanding the molecular basis underlying the clinical success of the cancer drug Gleevec. An "old" method, stopped-flow kinetics, combined with old enzymes, the ancestors dating back up to about billion years, provides an unexpected outlook for future intelligent design of drugs.

Binding energies of tyrosine kinase inhibitors: Error assessment of computational methods for imatinib and nilotinib binding

The binding energies of imatinib and nilotinib to tyrosine kinase have been determined by quantum mechanical (QM) computations, and compared with literature binding energy studies using molecular mechanics (MM). The potential errors in the computational methods include these critical factors. Errors in X-ray structures such as structural distortions and steric clashes give unrealistically high van der Waals energies, and erroneous binding energies.MM optimization gives a very different configuration to the QM optimization for nilotinib, whereas the imatinib ion gives similar configurations. Solvation energies are a major component of the overall binding energy. The QM based solvent model (PCM/SMD) gives different values from those used in the implicit PBSA solvent MM models. A major error in inhibitor—kinase binding lies in the non-polar solvation terms. Solvent transfer free energies and the required empirical solvent accessible surface area factors for nilotinib and imatinib ion to give the transfer free energies have been reverse calculated. These values differ from those used in the MM PBSA studies.An intertwined desolvation—conformational binding selectivity process is a balance of thermodynamic desolvation and intramolecular conformational kinetic control.The configurational entropies (TΔS) are minor error sources.

Kinase dynamics. Using ancient protein kinases to unravel a modern cancer drug's mechanism

Macromolecular function is rooted in energy landscapes, where sequence determines not a single structure but an ensemble of conformations. Hence, evolution modifies a protein's function by altering its energy landscape. Here, we recreate the evolutionary pathway between two modern human oncogenes, Src and Abl, by reconstructing their common ancestors. Our evolutionary reconstruction combined with x-ray structures of the common ancestor and pre-steady-state kinetics reveals a detailed atomistic mechanism for selectivity of the successful cancer drug Gleevec. Gleevec affinity is gained during the evolutionary trajectory toward Abl and lost toward Src, primarily by shifting an induced-fit equilibrium that is also disrupted in the clinical T315I resistance mutation. This work reveals the mechanism of Gleevec specificity while offering insights into how energy landscapes evolve.

Targeting conformational plasticity of protein kinases

The quest for ever more selective kinase inhibitors as potential future drugs has yielded a large repertoire of chemical probes that are selective for specific kinase conformations. These probes have been useful tools to obtain structural snapshots of kinase conformational plasticity. Similarly, kinetic and thermodynamic inhibitor binding experiments provide glimpses at the time scales and energetics of conformational interconversions. These experimental insights are complemented by computational predictions of conformational energy landscapes and simulations of conformational transitions and of the process of inhibitors binding to the protein kinase domain. A picture emerges in which highly selective inhibitors capitalize on the dynamic nature of kinases.

Computational study of the "DFG-flip" conformational transition in c-Abl and c-Src tyrosine kinases

Protein tyrosine kinases are crucial to cellular signaling pathways regulating cell growth, proliferation, metabolism, differentiation, and migration. To maintain normal regulation of cellular signal transductions, the activities of tyrosine kinases are also highly regulated. The conformation of a three-residue motif Asp-Phe-Gly (DFG) near the N-terminus of the long “activation” loop covering the catalytic site is known to have a critical impact on the activity of c-Abl and c-Src tyrosine kinases. A conformational transition of the DFG motif can switch the enzyme from an active (DFG-in) to an inactive (DFG-out) state. In the present study, the string method with swarms-of-trajectories was used to computationally determine the reaction pathway connecting the two end-states, and umbrella sampling calculations were carried out to characterize the thermodynamic factors affecting the conformations of the DFG motif in c-Abl and c-Src kinases. According to the calculated free energy landscapes, the DFG-out conformation is clearly more favorable in the case of c-Abl than that of c-Src. The calculations also show that the protonation state of the aspartate residue in the DFG motif strongly affects the in/out conformational transition in c-Abl, although it has a much smaller impact in the case of c-Src due to local structural differences.

Small-world networks of residue interactions in the Abl kinase complexes with cancer drugs: topology of allosteric communication pathways can determine drug resistance effects

Computational modelling of efficiency and robustness of the residue interaction networks and allosteric pathways in kinase structures can characterize protein kinase sensitivity to drug binding and drug resistance effects.

Conformational analysis of the DFG-out kinase motif and biochemical profiling of structurally validated type II inhibitors

Structural coverage of the human kinome has been steadily increasing over time. The structures provide valuable insights into the molecular basis of kinase function and also provide a foundation for understanding the mechanisms of kinase inhibitors. There are a large number of kinase structures in the PDB for which the Asp and Phe of the DFG motif on the activation loop swap positions, resulting in the formation of a new allosteric pocket. We refer to these structures as “classical DFG-out” conformations in order to distinguish them from conformations that have also been referred to as DFG-out in the literature but that do not have a fully formed allosteric pocket. We have completed a structural analysis of almost 200 small molecule inhibitors bound to classical DFG-out conformations; we find that they are recognized by both type I and type II inhibitors. In contrast, we find that nonclassical DFG-out conformations strongly select against type II inhibitors because these structures have not formed a large enough allosteric pocket to accommodate this type of binding mode. In the course of this study we discovered that the number of structurally validated type II inhibitors that can be found in the PDB and that are also represented in publicly available biochemical profiling studies of kinase inhibitors is very small. We have obtained new profiling results for several additional structurally validated type II inhibitors identified through our conformational analysis. Although the available profiling data for type II inhibitors is still much smaller than for type I inhibitors, a comparison of the two data sets supports the conclusion that type II inhibitors are more selective than type I. We comment on the possible contribution of the DFG-in to DFG-out conformational reorganization to the selectivity.

Computational study of Gleevec and G6G reveals molecular determinants of kinase inhibitor selectivity

Gleevec is a potent inhibitor of Abl tyrosine kinase but not of the highly homologous c-Src kinase. Because the ligand binds to an inactive form of the protein in which an Asp-Phe-Gly structural motif along the activation loop adopts a so-called DFG-out conformation, it was suggested that binding specificity was controlled by a "conformational selection" mechanism. In this context, the binding affinity displayed by the kinase inhibitor G6G poses an intriguing challenge. Although it possesses a chemical core very similar to that of Gleevec, G6G is a potent inhibitor of both Abl and c-Src kinases. Both inhibitors bind to the DFG-out conformation of the kinases, which seems to be in contradiction with the conformational selection mechanism. To address this issue and display the hidden thermodynamic contributions affecting the binding selectivity, molecular dynamics free energy simulations with explicit solvent molecules were carried out. Relative to Gleevec, G6G forms highly favorable van der Waals dispersive interactions upon binding to the kinases via its triazine functional group, which is considerably larger than the corresponding pyridine moiety in Gleevec. Upon binding of G6G to c-Src, these interactions offset the unfavorable free energy cost of the DFG-out conformation. When binding to Abl, however, G6G experiences an unfavorable free energy penalty due to steric clashes with the phosphate-binding loop, yielding an overall binding affinity that is similar to that of Gleevec. Such steric clashes are absent when G6G binds to c-Src, due to the extended conformation of the phosphate-binding loop.

Energetic dissection of Gleevec's selectivity toward human tyrosine kinases

Protein kinases are obvious drug targets against cancer due to their central role in cellular regulation. Since the discovery of Gleevec, a potent and specific inhibitor of Abl kinase, as a highly successful cancer therapeutic, the ability of this drug to distinguish between Abl and other tyrosine kinases like Src has been intensely investigated, but without much success. Using NMR and fast kinetics, we establish a novel model that solves this longstanding question of two tyrosine kinases adopting almost identical structures when bound to Gleevec, yet having vastly different affinities. In contrast to all other proposed models we show that the origin of Abl’s high affinity lies predominantly in a conformational change after binding. An energy landscape that provides tight affinity via an induced-fit and binding plasticity via conformational selection mechanism is likely to be general for many inhibitors.

BCR-ABL residues interacting with ponatinib are critical to preserve the tumorigenic potential of the oncoprotein

Patients with chronic myeloid leukemia in whom tyrosine kinase inhibitors (TKIs) fail often present mutations in the BCR-ABL catalytic domain. We noticed a lack of substitutions involving 4 amino acids (E286, M318, I360, and D381) that form hydrogen bonds with ponatinib. We therefore introduced mutations in each of these residues, either preserving or altering their physicochemical properties. We found that E286, M318, I360, and D381 are dispensable for ABL and BCR-ABL protein stability but are critical for preserving catalytic activity. Indeed, only a "conservative" I360T substitution retained kinase proficiency and transforming potential. Molecular dynamics simulations of BCR-ABL(I360T) revealed differences in both helix αC dynamics and protein-correlated motions, consistent with a modified ATP-binding pocket. Nevertheless, this mutant remained sensitive to ponatinib, imatinib, and dasatinib. These results suggest that changes in the 4 BCR-ABL residues described here would be selected against by a lack of kinase activity or by maintained responsiveness to TKIs. Notably, amino acids equivalent to those identified in BCR-ABL are conserved in 51% of human tyrosine kinases. Hence, these residues may represent an appealing target for the design of pharmacological compounds that would inhibit additional oncogenic tyrosine kinases while avoiding the emergence of resistance due to point mutations.

Monte Carlo simulations of proteins at constant pH with generalized Born solvent, flexible sidechains, and an effective dielectric boundary

Titratable residues determine the acid/base behavior of proteins, strongly influencing their function; in addition, proton binding is a valuable reporter on electrostatic interactions. We describe a method for pK(a) calculations, using constant-pH Monte Carlo (MC) simulations to explore the space of sidechain conformations and protonation states, with an efficient and accurate generalized Born model (GB) for the solvent effects. To overcome the many-body dependency of the GB model, we use a "Native Environment" approximation, whose accuracy is shown to be good. It allows the precalculation and storage of interactions between all sidechain pairs, a strategy borrowed from computational protein design, which makes the MC simulations themselves very fast. The method is tested for 12 proteins and 167 titratable sidechains. It gives an rms error of 1.1 pH units, similar to the trivial "Null" model. The only adjustable parameter is the protein dielectric constant. The best accuracy is achieved for values between 4 and 8, a range that is physically plausible for a protein interior. For sidechains with large pKa shifts, ≥2, the rms error is 1.6, compared to 2.5 with the Null model and 1.5 with the empirical PROPKA method.

Computational analysis of the binding specificity of Gleevec to Abl, c-Kit, Lck, and c-Src tyrosine kinases

Gleevec, a well-known cancer therapeutic agent, is an effective inhibitor of several tyrosine kinases, including Abl and c-Kit, but displays less potency to inhibit closely homologous tyrosine kinases, such as Lck and c-Src. Because many structural features of the binding site are highly conserved in these homologous kinases, the molecular determinants responsible for the binding specificity of Gleevec remain poorly understood. To address this issue, free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent was used to compute the binding affinity of Gleevec to Abl, c-Kit, Lck, and c-Src. The results of the FEP/MD calculations are in good agreement with experiments, enabling a detailed and quantitative dissection of the absolute binding free energy in terms of various thermodynamic contributions affecting the binding specificity of Gleevec to the kinases. Dominant binding free energy contributions arises from the van der Waals dispersive interaction, compensating about two-thirds of the unfavorable free energy penalty associated with the loss of translational, rotational, and conformational freedom of the ligand upon binding. In contrast, the contributions from electrostatic and repulsive interactions nearly cancel out due to solvent effects. Furthermore, the calculations show the importance of the conformation of the kinase activation loop. Among the kinases examined, Abl provides the most favorable binding environment for Gleevec via optimal protein-ligand interactions and a small free energy cost for loss of the translational, rotational, and conformational freedom upon ligand binding. The FEP/MD calculations additionally reveal that Lck and c-Src provide similar nonbinding interactions with the bound-Gleevec, but the former pays less entropic penalty for the ligand losing its translational, rotational, and conformational motions to bind, examining the empirically observed differential binding affinities of Gleevec between the two Src-family kinases.