Sequence and Structure Signatures of Cancer Mutation Hotspots in Protein Kinases

Protein kinases are the most common protein domains implicated in cancer, where somatically acquired mutations are known to be functionally linked to a variety of cancers. Resequencing studies of protein kinase coding regions have emphasized the importance of sequence and structure determinants of cancer-causing kinase mutations in understanding of the mutation-dependent activation process. We have developed an integrated bioinformatics resource, which consolidated and mapped all currently available information on genetic modifications in protein kinase genes with sequence, structure and functional data. The integration of diverse data types provided a convenient framework for kinome-wide study of sequence-based and structure-based signatures of cancer mutations. The database-driven analysis has revealed a differential enrichment of SNPs categories in functional regions of the kinase domain, demonstrating that a significant number of cancer mutations could fall at structurally equivalent positions (mutational hotspots) within the catalytic core. We have also found that structurally conserved mutational hotspots can be shared by multiple kinase genes and are often enriched by cancer driver mutations with high oncogenic activity. Structural modeling and energetic analysis of the mutational hotspots have suggested a common molecular mechanism of kinase activation by cancer mutations, and have allowed to reconcile the experimental data. According to a proposed mechanism, structural effect of kinase mutations with a high oncogenic potential may manifest in a significant destabilization of the autoinhibited kinase form, which is likely to drive tumorigenesis at some level. Structure-based functional annotation and prediction of cancer mutation effects in protein kinases can facilitate an understanding of the mutation-dependent activation process and inform experimental studies exploring molecular pathology of tumorigenesis.     

Hierarchical Modeling of Activation Mechanisms in the ABL and EGFR Kinase Domains: Thermodynamic and Mechanistic Catalysts of Kinase Activation by Cancer Mutations

Structural and functional studies of the ABL and EGFR kinase domains have recently suggested a common mechanism of activation by cancer-causing mutations. However, dynamics and mechanistic aspects of kinase activation by cancer mutations that stimulate conformational transitions and thermodynamic stabilization of the constitutively active kinase form remain elusive. We present a large-scale computational investigation of activation mechanisms in the ABL and EGFR kinase domains by a panel of clinically important cancer mutants ABL-T315I, ABL-L387M, EGFR-T790M, and EGFR-L858R. We have also simulated the activating effect of the gatekeeper mutation on conformational dynamics and allosteric interactions in functional states of the ABL-SH2-SH3 regulatory complexes. A comprehensive analysis was conducted using a hierarchy of computational approaches that included homology modeling, molecular dynamics simulations, protein stability analysis, targeted molecular dynamics, and molecular docking. Collectively, the results of this study have revealed thermodynamic and mechanistic catalysts of kinase activation by major cancer-causing mutations in the ABL and EGFR kinase domains. By using multiple crystallographic states of ABL and EGFR, computer simulations have allowed one to map dynamics of conformational fluctuations and transitions in the normal (wild-type) and oncogenic kinase forms. A proposed multi-stage mechanistic model of activation involves a series of cooperative transitions between different conformational states, including assembly of the hydrophobic spine, the formation of the Src-like intermediate structure, and a cooperative breakage and formation of characteristic salt bridges, which signify transition to the active kinase form. We suggest that molecular mechanisms of activation by cancer mutations could mimic the activation process of the normal kinase, yet exploiting conserved structural catalysts to accelerate a conformational transition and the enhanced stabilization of the active kinase form. The results of this study reconcile current experimental data with insights from theoretical approaches, pointing to general mechanistic aspects of activating transitions in protein kinases.     

The energy landscape analysis of cancer mutations in protein kinases

The growing interest in quantifying the molecular basis of protein kinase activation and allosteric regulation by cancer mutations has fueled computational studies of allosteric signaling in protein kinases. In the present study, we combined computer simulations and the energy landscape analysis of protein kinases to characterize the interplay between oncogenic mutations and locally frustrated sites as important catalysts of allostetric kinase activation. While structurally rigid kinase core constitutes a minimally frustrated hub of the catalytic domain, locally frustrated residue clusters, whose interaction networks are not energetically optimized, are prone to dynamic modulation and could enable allosteric conformational transitions. The results of this study have shown that the energy landscape effect of oncogenic mutations may be allosteric eliciting global changes in the spatial distribution of highly frustrated residues. We have found that mutation-induced allosteric signaling may involve a dynamic coupling between structurally rigid (minimally frustrated) and plastic (locally frustrated) clusters of residues. The presented study has demonstrated that activation cancer mutations may affect the thermodynamic equilibrium between kinase states by allosterically altering the distribution of locally frustrated sites and increasing the local frustration in the inactive form, while eliminating locally frustrated sites and restoring structural rigidity of the active form. The energy landsape analysis of protein kinases and the proposed role of locally frustrated sites in activation mechanisms may have useful implications for bioinformatics-based screening and detection of functional sites critical for allosteric regulation in complex biomolecular systems.     

Mutational activation of ErbB2 reveals a new protein kinase autoinhibition mechanism

Autoinhibition plays a key role in the control of protein kinase activity. ErbB2 is a unique receptor-tyrosine kinase that does not bind ligand but possesses an extracellular domain poised to engage other ErbBs. Little is known about the molecular mechanism for ErbB2 catalytic regulation. Here we show that ErbB2 kinase is strongly autoinhibited, and a loop connecting the alphaC helix and beta4 sheet within the kinase domain plays a major role in the control of kinase activity. Mutations of two Gly residues at positions 776 and 778 in this loop dramatically increase ErbB2 catalytic activity. Kinetic analysis demonstrates that mutational activation is due to approximately 10- and approximately 7-fold increases in ATP binding affinity and turnover number, respectively. Expression of the activated ErbB2 mutants in cells resulted in elevated ligand-independent ErbB2 autophosphorylation, ErbB3 phosphorylation, and stimulation of mitogen-activated protein kinase. Molecular modeling suggests that the ErbB2 kinase domain is stabilized in an inactive state via a hydrophobic interaction between the alphaC-beta4 and activation loops. Importantly, many ErbB2 human cancer mutations have been identified in the alphaC-beta4 loop, including the activating G776S mutation studied here. Our findings reveal a new kinase regulatory mechanism in which the alphaC-beta4 loop functions as an intramolecular switch that controls ErbB2 activity and suggests that loss of alphaC-beta4 loop-mediated autoinhibition is involved in oncogenic activation of ErbB2.     

In silico profiling of tyrosine kinases binding specificity and drug resistance using Monte Carlo simulations with the ensembles of protein kinase crystal structures

The molecular basis of the tyrosine kinases binding specificity and drug resistance against cancer drugs Imatinib and Dasatinib is elucidated using Monte Carlo simulations of the inhibitor-receptor binding with the ensembles of protein kinase crystal structures. In silico proteomics analysis unravels mechanisms by which structural plasticity of the tyrosine kinases is linked with the conformational preferences of Imatinib and Dasatinib in achieving effective drug binding with a distinct spectrum of the tyrosine kinome. The differences in the inhibitor sensitivities to the ABL kinase mutants are rationalized based on variations in the binding free energy profiles with the conformational states of the ABL kinase. While Imatinib binding is highly sensitive to the activation state of the enzyme, the computed binding profile of Dasatinib is remarkably tolerant to the conformational state of ABL. A comparative analysis of the inhibitor binding profiles with the clinically important ABL kinase mutants has revealed an excellent agreement with the biochemical and proteomics data. We have found that conformational adaptability of the kinase inhibitors to structurally different conformational states of the tyrosine kinases may have pharmacological relevance in acquiring a specific array of potent activities and regulating a scope of the inhibitor resistance mutations. This study outlines a useful approach for understanding and predicting the molecular basis of the inhibitor sensitivity against potential kinase targets and drug resistance.     

Deciphering adaptor specificity in GFL-dependent RET-mediated proliferation and neurite outgrowth

Glial cell derived neurotrophic factor (GDNF)-dependent receptor tyrosine kinase RET activity is required for proper development of the nervous system and genitourinary tract. Loss-of-function mutations in RET are associated with enteric nervous system abnormalities (Hirschsprung disease) and renal deficits (Potter's syndrome), whereas activating mutations lead to hereditary cancer syndromes (multiple endocrine neoplasia type 2A and type 2B). RET activation is crucial for the proper regulation of a variety of cellular processes including cell migration, proliferation and neurite outgrowth. By analyzing a series of RET mutants we found that Y1062 was critical for stimulating GDNF-mediated proliferation as well as proliferation stimulated by GDNF-independent oncogenic forms of RET. Studies using small interfering RNA driven by lentivirus to knock-down expression of particular adaptor proteins that interact with RET phospho-Y1062, demonstrated that only Src-homology 2 and growth factor receptor binding protein 2 were necessary for RET-mediated proliferation by wild type and oncogenic forms of RET. Interestingly, we discovered that Y1062 was also required for GDNF-stimulated neurite outgrowth. However, small interfering RNAs to either Src-homology 2 or growth factor receptor binding protein 2 or a panel of other adaptor proteins known to interact with RET Y1062 were incapable of blocking GDNF-stimulated neurite formation, indicating that differential use of intracellular adaptors is responsible for regulating alternative RET-stimulated cellular events such as proliferation versus a differentiation response like neurite outgrowth.     

Identification of Novel Small Molecule Inhibitors of Oncogenic RET Kinase

Oncogenic mutation of the RET receptor tyrosine kinase is observed in several human malignancies. Here, we describe three novel type II RET tyrosine kinase inhibitors (TKI), ALW-II-41-27, XMD15-44 and HG-6-63-01, that inhibit the cellular activity of oncogenic RET mutants at two digit nanomolar concentration. These three compounds shared a 3-trifluoromethyl-4-methylpiperazinephenyl pharmacophore that stabilizes the ‘DFG-out’ inactive conformation of RET activation loop. They blocked RET-mediated signaling and proliferation with an IC50 in the nM range in fibroblasts transformed by the RET/C634R and RET/M918T oncogenes. They also inhibited autophosphorylation of several additional oncogenic RET-derived point mutants and chimeric oncogenes. At a concentration of 10 nM, ALW-II-41-27, XMD15-44 and HG-6-63-01 inhibited RET kinase and signaling in human thyroid cancer cell lines carrying oncogenic RET alleles; they also inhibited proliferation of cancer, but not non-tumoral Nthy-ori-3-1, thyroid cells, with an IC50 in the nM range. The three compounds were capable of inhibiting the ‘gatekeeper’ V804M mutant which confers substantial resistance to established RET inhibitors. In conclusion, we have identified a type II TKI scaffold, shared by ALW-II-41-27, XMD15-44 and HG-6-63-01, that may be used as novel lead for the development of novel agents for the treatment of cancers harboring oncogenic activation of RET.     

Inducible dimerization of RET reveals a specific AKT deregulation in oncogenic signaling

Dominant-activating mutations in the RET (rearranged during transfection) proto-oncogene, a receptor tyrosine kinase, are causally associated with the development of multiple endocrine neoplasia type 2A (MEN2A) syndrome. Such oncogenic RET mutations induce its ligand-independent constitutive activation, but whether it spreads identical signaling to ligand-induced signaling is uncertain. To address this question, we designed a cellular model in which RET can be activated either by its natural ligand, or alternatively, by controlled dimerization of the protein that mimics MEN2A dimerization. We have shown that controlled dimerization leaves proximal RET signaling intact but impacts substantially on the tuning of the distal AKT kinase activation (delayed and sustained). In marked contrast, distal activation of ERK remained unaffected. We further demonstrated that specific temporal adjustment of ligand-induced AKT activation is dependent upon a lipid-based cholesterol-sensitive environment, and this control step is bypassed by MEN2A RET mutants. Therefore, these studies revealed that MEN2A mutations propagate previously unappreciated subtle differences in signaling pathways and unravel a role for lipid rafts in the temporal regulation of AKT activation.     

Docking protein FRS2 links the protein tyrosine kinase RET and its oncogenic forms with the mitogen-activated protein kinase signaling cascade

The receptor tyrosine kinase RET functions as the signal transducing receptor for the GDNF (for "glial cell-derived neurotrophic factors") family of ligands. Mutations in the RET gene were implicated in Hirschsprung disease (HSCR), multiple endocrine neoplasia type 2 (MEN 2), and thyroid carcinomas. In this report we demonstrate that the docking protein FRS2 is tyrosine phosphorylated by ligand-stimulated and by constitutively activated oncogenic forms of RET. Complex formation between RET and FRS2 is mediated by binding of the phosphotyrosine-binding domain of FRS2 to pY1062, a residue in RET that also functions as a binding site for Shc. However, overexpression of FRS2 but not Shc potentiates mitogen-activated protein (MAP) kinase activation by RET oncoproteins. We demonstrate that oncogenic RET-PTC proteins are associated with FRS2 constitutively, leading to tyrosine phosphorylation of FRS2, MAP kinase stimulation, and cell proliferation. However, loss-of-function HSCR-associated RET mutants exhibit impaired FRS2 binding and reduced MAP kinase activation. These experiments demonstrate that FRS2 couples both ligand-regulated and oncogenic forms of RET, with the MAP kinase signaling cascade as part of the response of RET under normal biological conditions and pathological conditions, such as MEN 2 and papillary thyroid carcinomas.     

Expression, purification, and inhibition of human RET tyrosine kinase

Tyrosine kinases are emerging as frequent targets of primary oncogenic events and therefore represent an optimal focus of therapeutical intervention. Genetic alterations that cause dysregulated activation of the RET tyrosine kinase are responsible for a significant fraction of thyroid carcinomas. In an effort towards therapeutic RET inactivation, we have developed a method for expression and purification of recombinant RET catalytic domain for structural purposes and for use in the screening of potential inhibitors of RET kinase activity. His-tagged RET kinase domain was purified from Sf9 insect cell lysate using a two-step chromatographic protocol and characterised. Purified recombinant RET phosphorylated itself and exogenous substrates at physiological pH. A specific peptide substrate, derived from RET activation loop, was identified and experimentally validated. These reagents were used to develop a rapid ELISA-based kinase assay for screening potential inhibitors. Novel RET inhibitors were identified using this assay.     

A Src-like inactive conformation in the abl tyrosine kinase domain

The improper activation of the Abl tyrosine kinase results in chronic myeloid leukemia (CML). The recognition of an inactive conformation of Abl, in which a catalytically important Asp-Phe-Gly (DFG) motif is flipped by approximately 180 degrees with respect to the active conformation, underlies the specificity of the cancer drug imatinib, which is used to treat CML. The DFG motif is not flipped in crystal structures of inactive forms of the closely related Src kinases, and imatinib does not inhibit c-Src. We present a structure of the kinase domain of Abl, determined in complex with an ATP-peptide conjugate, in which the protein adopts an inactive conformation that resembles closely that of the Src kinases. An interesting aspect of the Src-like inactive structure, suggested by molecular dynamics simulations and additional crystal structures, is the presence of features that might facilitate the flip of the DFG motif by providing room for the phenylalanine to move and by coordinating the aspartate side chain as it leaves the active site. One class of mutations in BCR-Abl that confers resistance to imatinib appears more likely to destabilize the inactive Src-like conformation than the active or imatinib-bound conformations. Our results suggest that interconversion between distinctly different inactive conformations is a characteristic feature of the Abl kinase domain.     

Genetic factors predisposing to the development of papillary thyroid cancer

According to classic theory of neogenesis, cancer arises from well-differentiated cell that in response to variety of factors de-differentiates, becomes able to proliferate without control and/or loses its ability to undergo apoptosis. According to another theory, cancers (at least cancers of some organs) originate from stem cells, which "by definition" are poorly differentiated and able to proliferate indefinitely. Therefore a lower number of abnormal events is necessary for these cells to escape proliferation-controlling mechanisms. With regard to papillary thyroid cancers it is still thought that it arises from well-differentiated thyreocyte. One of the characteristic features of cancer cell is chromosomal instability. Lowest number of such abnormalities is observed in well-differentiated thyroid cancers (including papillary cancer), intermediate - in poorly-differentiated cancers, while highest - in anaplastic cancers. Microarray analysis shows that despite of clinical heterogeneity, gene expression profiles of papillary cancers are very similar. Genetic anomalies predisposing to the development of papillary cancer most commonly regard proteins that possess kinase activity. Kinases phosphorylate other proteins, and play an extremely important role in signal transduction from outside the cell as well as inside the cell. Constitutive activation of some kinases may lead to the excessive and/or permanent activation of some transduction pathways specific for mitogens or growth factors. This results in excessive proliferation. The best known protein of such type which function is altered in papillary thyroid cancers is RET - a membrane-located growth factor-receptor with kinase activity. RET gene undergoes different rearrangements in this type of cancer. There are approximately 10 RET rearrangements known, with RET/PTC3 and RET/PTC1 being most common. In this anomaly kinase domain-encoding 3' end of RET gene is aberrantly bound to 5' end of another gene. Fusion protein synthesized on such hybrid template is not present in the cell membrane but in the cytoplasm, where it permanently activates transduction pathway specific for RET. NTRK1 gene encoding a member of family of neuronal growth factor receptors containing thyrosine kinase domain is also rearranged in papillary cancers. However, genes fused to its kinase domain-encoding sequence are different from the ones fused to RET. MET, a gene encoding another membrane protein with thyrosine kinase activity, which acts as a growth factor-receptor, is overexpressed in 70%-90% of papillary thyroid cancers. BRAF gene encoding another yet kinase transducing signals from RAS and RAF to the cell is mutated at position 1796 (T/A, amino acid substitution V599E) in 38-69% of papillary cancers. The presence of this activatory mutation is associated with higher degree of clinical advancement of the disease. In addition, in majority of papillary cancers tested, mutations of the genes encoding nuclear triiodothyronine receptors were found. Transgenic mice with both TRB allele replaced with dominant-negative TRB mutants develop aggressive thyroid cancers. Progression from papillary to anaplastic cancer is most possibly caused by the occurrence of additional anomalies within P53, RAS, NM23,b-catenin gene and other genes.     

MNK,EIF4E and targeting translation for therapy

Deregulation of protein translation is a common event in cancer and occurs frequently as a result of mutational activation of the AKT signaling pathway. We had previously reported the in vivo oncogenic activity of the translation initiation factor eIF4E, which acts downstream AKT and mTOR. We now identified an absolute requirement for Ser209 phosphorylation by the MNK1/2 kinases for eIF4E’s oncogenic action. MNK1/2 kinases are dispensable for normal development in mammals. This potential difference between normal and cancer cells may provide a therapeutic avenue for targeting translational requirements in cancer.     

Differences in signaling pathways and expression level of the phosphoinositide phosphatase SHIP1 between two oncogenic mutants of the receptor tyrosine kinase KIT

Oncogenic mutations of the receptor tyrosine kinase KIT are encountered in myeloid leukemia and various solid tumors, including gastrointestinal stromal tumors. We previously identified the human oncogenic germ line mutant KIT(K642E), a substitution in the tyrosine kinase 1 domain (TK1D) in a familial form of gastrointestinal stromal tumors. The effects of oncogenic KIT mutants on cell signaling and regulation are complex. Cellular models are valuable basic tools to tailor novel strategies on specific cellular and molecular bases for tumors expressing KIT oncogenic mutants. Murine KIT(WT) and the murine homologues of human KIT oncogenic mutants, further referred to as KIT(K641E) and KIT(del559), a point deletion in the juxtamembrane domain (JMD), were stably expressed in IL-3-dependent Ba/F3 cells. Major differences in the constitutively activation of Akt/PKB, MAP kinases and STATs pathways were observed between KIT(K641E) and KIT(del559), whereas KIT ligand elicited responses in both mutants. Noteworthy, the protein level of the phosphoinositide phosphatase SHIP1, but not SHIP2 and PTEN, was reduced in KIT(K641E) only while inhibition of KIT phosphorylation reversibly raised SHIP1 level in both JMD and TK1D oncogenic mutants, unraveling the control of SHIP protein level by KIT phosphorylation.     

Ras-mutant Cancer Cells Display B-Raf Binding to Ras That Activates Extracellular Signal-regulated Kinase and Is Inhibited by Protein Kinase A Phosphorylation

The small G protein Ras regulates proliferation through activation of the mitogen-activated protein (MAP) kinase (ERK) cascade. The first step of Ras-dependent activation of ERK signaling is Ras binding to members of the Raf family of MAP kinase kinase kinases, C-Raf and B-Raf. Recently, it has been reported that in melanoma cells harboring oncogenic Ras mutations, B-Raf does not bind to Ras and does not contribute to basal ERK activation. For other types of Ras-mutant tumors, the relative contributions of C-Raf and B-Raf are not known. We examined non-melanoma cancer cell lines containing oncogenic Ras mutations and express both C-Raf and B-Raf isoforms, including the lung cancer cell line H1299 cells. Both B-Raf and C-Raf were constitutively bound to oncogenic Ras and contributed to Ras-dependent ERK activation. Ras binding to B-Raf and C-Raf were both subject to inhibition by the cAMP-dependent protein kinase PKA. cAMP inhibited the growth of H1299 cells and Ras-dependent ERK activation via PKA. PKA inhibited the binding of Ras to both C-Raf and B-Raf through phosphorylations of C-Raf at Ser-259 and B-Raf at Ser-365, respectively. These studies demonstrate that in non-melanocytic Ras-mutant cancer cells, Ras signaling to B-Raf is a significant contributor to ERK activation and that the B-Raf pathway, like that of C-Raf, is a target for inhibition by PKA. We suggest that cAMP and hormones coupled to cAMP may prove useful in dampening the effects of oncogenic Ras in non-melanocytic cancer cells through PKA-dependent actions on B-Raf as well as C-Raf.     

Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases

The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level.     

Conformational dynamics of the EGFR kinase domain reveals structural features involved in activation

The epidermal growth factor receptor (EGFR) has been the focus of intensive studies because of its importance in cancer research. Thus, a broader understanding of the molecular mechanism of activation of the EGFR kinase will have profound significance for the development of novel therapeutics. Numerous crystal structures of EGFR kinase, including the structure of the activating‐kinase dimer, have provided snapshots of the specific pathway. Herein, we performed unrestrained‐, as well as targeted‐molecular dynamics simulations based on these data, to gain further insight into the conformational changes responsible for activation. Comparison of the monomer‐ versus activating‐EGFR‐dimer simulations indicates that the dimerization is stabilizing structural elements associated with the activated state and predicts new salt‐bridge interactions involving activation‐loop residues that may also be associated with that state. Targeted molecular dynamics simulations of the inactive‐to‐active EGFR transition, as well as the reverse pathway, confirm the formation of conserved structural features of functional importance for the activity or stabilization of either conformation. Interestingly, simulations of the L834R mutant, which is associated with cancer, suggest that the structural basis of the activation induced by that mutation might be the ability of the mutated R834 residue to consecutively form salt bridges with neighboring acidic residues and cause destabilization of a hydrophobic cluster in the inactive state. Proteins 2009. © 2008 Wiley‐Liss, Inc.