A Combination of Conformation-Specific RAF Inhibitors Overcome Drug Resistance Brought about by RAF Overexpression

Cancer cells often adapt to targeted therapies, yet the molecular mechanisms underlying adaptive resistance remain only partially understood. Here, we explore a mechanism of RAS/RAF/MEK/ERK (MAPK) pathway reactivation through the upregulation of RAF isoform (RAFs) abundance. Using computational modeling and in vitro experiments, we show that the upregulation of RAFs changes the concentration range of paradoxical pathway activation upon treatment with conformation-specific RAF inhibitors. Additionally, our data indicate that the signaling output upon loss or downregulation of one RAF isoform can be compensated by overexpression of other RAF isoforms. We furthermore demonstrate that, while single RAF inhibitors cannot efficiently inhibit ERK reactivation caused by RAF overexpression, a combination of two structurally distinct RAF inhibitors synergizes to robustly suppress pathway reactivation.

Kinase Suppressor of RAS 1 (KSR1) Maintains the Transformed Phenotype of BRAFV600E Mutant Human Melanoma Cells

Kinase Suppressor of RAS 1 (KSR1) is a scaffolding protein for the RAS-RAF-MEK-ERK pathway, which is one of the most frequently altered pathways in human cancers. Previous results have shown that KSR1 has a critical role in mutant RAS-mediated transformation. Here, we examined the role of KSR1 in mutant BRAF transformation. We used CRISPR/Cas9 to knock out KSR1 in a BRAFV600E-transformed melanoma cell line. KSR1 loss produced a complex phenotype characterised by impaired proliferation, cell cycle defects, decreased transformation, decreased invasive migration, increased cellular senescence, and increased apoptosis. To decipher this phenotype, we used a combination of proteomic ERK substrate profiling, global protein expression profiling, and biochemical validation assays. The results suggest that KSR1 directs ERK to phosphorylate substrates that have a critical role in ensuring cell survival. The results further indicate that KSR1 loss induces the activation of p38 Mitogen-Activated Protein Kinase (MAPK) and subsequent cell cycle aberrations and senescence. In summary, KSR1 function plays a key role in oncogenic BRAF transformation.

Interactome dynamics of RAF1-BRAF kinase monomers and dimers

RAF kinases play major roles in cancer. BRAFV600E mutants drive ~6% of human cancers. Potent kinase inhibitors exist but show variable effects in different cancer types, sometimes even inducing paradoxical RAF kinase activation. Both paradoxical activation and drug resistance are frequently due to enhanced dimerization between RAF1 and BRAF, which maintains or restores the activity of the downstream MEK-ERK pathway. Here, using quantitative proteomics we mapped the interactomes of RAF1 monomers, RAF1-BRAF and RAF1-BRAFV600E dimers identifying and quantifying >1,000 proteins. In addition, we examined the effects of vemurafenib and sorafenib, two different types of clinically used RAF inhibitors. Using regression analysis to compare different conditions we found a large overlapping core interactome but also distinct condition specific differences. Given that RAF proteins have kinase independent functions such dynamic interactome changes could contribute to their functional diversification. Analysing this dataset may provide a deeper understanding of RAF signalling and mechanisms of resistance to RAF inhibitors.

Control of cell state transitions

Understanding cell state transitions and purposefully controlling them is a longstanding challenge in biology. Here we present cell state transition assessment and regulation (cSTAR), an approach for mapping cell states, modelling transitions between them and predicting targeted interventions to convert cell fate decisions. cSTAR uses omics data as input, classifies cell states, and develops a workflow that transforms the input data into mechanistic models that identify a core signalling network, which controls cell fate transitions by influencing whole-cell networks. By integrating signalling and phenotypic data, cSTAR models how cells manoeuvre in Waddington's landscape1 and make decisions about which cell fate to adopt. Notably, cSTAR devises interventions to control the movement of cells in Waddington's landscape. Testing cSTAR in a cellular model of differentiation and proliferation shows a high correlation between quantitative predictions and experimental data. Applying cSTAR to different types of perturbation and omics datasets, including single-cell data, demonstrates its flexibility and scalability and provides new biological insights. The ability of cSTAR to identify targeted perturbations that interconvert cell fates will enable designer approaches for manipulating cellular development pathways and mechanistically underpinned therapeutic interventions.

A systematic analysis of signaling reactivation and drug resistance

Increasing evidence suggests that the reactivation of initially inhibited signaling pathways causes drug resistance. Here, we analyze how network topologies affect signaling responses to drug treatment. Network-dependent drug resistance is commonly attributed to negative and positive feedback loops. However, feedback loops by themselves cannot completely reactivate steady-state signaling. Newly synthesized negative feedback regulators can induce a transient overshoot but cannot fully restore output signaling. Complete signaling reactivation can only occur when at least two routes, an activating and inhibitory, connect an inhibited upstream protein to a downstream output. Irrespective of the network topology, drug-induced overexpression or increase in target dimerization can restore or even paradoxically increase downstream pathway activity. Kinase dimerization cooperates with inhibitor-mediated alleviation of negative feedback. Our findings inform drug development by considering network context and optimizing the design drug combinations. As an example, we predict and experimentally confirm specific combinations of RAF inhibitors that block mutant NRAS signaling.

Kholodenko et al. uncover signaling network circuitries and molecular mechanisms necessary and sufficient for complete reactivation or overshoot of steady-state signaling after kinase inhibitor treatment. The two means to revive signaling output fully are through network topology or reactivation of the kinase activity of the primary drug target. Blocking RAF dimer activity by a combination of type I½ and type II RAF inhibitors efficiently blocks mutant NRAS-driven ERK signaling.

Extensive rewiring of the EGFR network in colorectal cancer cells expressing transforming levels of KRASG13D

Protein-protein-interaction networks (PPINs) organize fundamental biological processes, but how oncogenic mutations impact these interactions and their functions at a network-level scale is poorly understood. Here, we analyze how a common oncogenic KRAS mutation (KRASG13D) affects PPIN structure and function of the Epidermal Growth Factor Receptor (EGFR) network in colorectal cancer (CRC) cells. Mapping >6000 PPIs shows that this network is extensively rewired in cells expressing transforming levels of KRASG13D (mtKRAS). The factors driving PPIN rewiring are multifactorial including changes in protein expression and phosphorylation. Mathematical modelling also suggests that the binding dynamics of low and high affinity KRAS interactors contribute to rewiring. PPIN rewiring substantially alters the composition of protein complexes, signal flow, transcriptional regulation, and cellular phenotype. These changes are validated by targeted and global experimental analysis. Importantly, genetic alterations in the most extensively rewired PPIN nodes occur frequently in CRC and are prognostic of poor patient outcomes.

PHD3 Regulates p53 Protein Stability by Hydroxylating Proline 359

Cellular p53 protein levels are regulated by a ubiquitination/de-ubiquitination cycle that can target the protein for proteasomal destruction. The ubiquitination reaction is catalyzed by a multitude of ligases, whereas the removal of ubiquitin chains is mediated by two deubiquitinating enzymes (DUBs), USP7 (HAUSP) and USP10. Here, we show that PHD3 hydroxylates p53 at proline 359, a residue that is in the p53-DUB binding domain. Hydroxylation of p53 upon proline 359 regulates its interaction with USP7 and USP10, and its inhibition decreases the association of p53 with USP7/USP10, increases p53 ubiquitination, and rapidly reduces p53 protein levels independently of mRNA expression. Our results show that p53 is a PHD3 substrate and that hydroxylation by PHD3 regulates p53 protein stability through modulation of ubiquitination.

Metabolic stress regulates ERK activity by controlling KSR-RAF heterodimerization

Altered cell metabolism is a hallmark of cancer, and targeting specific metabolic nodes is considered an attractive strategy for cancer therapy. In this study, we evaluate the effects of metabolic stressors on the deregulated ERK pathway in melanoma cells bearing activating mutations of the NRAS or BRAF oncogenes. We report that metabolic stressors promote the dimerization of KSR proteins with CRAF in NRAS-mutant cells, and with oncogenic BRAF in BRAFV600E-mutant cells, thereby enhancing ERK pathway activation. Despite this similarity, the two genomic subtypes react differently when a higher level of metabolic stress is induced. In NRAS-mutant cells, the ERK pathway is even more stimulated, while it is strongly downregulated in BRAFV600E-mutant cells. We demonstrate that this is caused by the dissociation of mutant BRAF from KSR and is mediated by activated AMPK. Both types of ERK regulation nevertheless lead to cell cycle arrest. Besides studying the effects of the metabolic stressors on ERK pathway activity, we also present data suggesting that for efficient therapies of both genomic melanoma subtypes, specific metabolic targeting is necessary.

A novel RNA sequencing data analysis method for cell line authentication

We have developed a novel analysis method that can interrogate the authenticity of biological samples used for generation of transcriptome profiles in public data repositories. The method uses RNA sequencing information to reveal mutations in expressed transcripts and subsequently confirms the identity of analysed cells by comparison with publicly available cell-specific mutational profiles. Cell lines constitute key model systems widely used within cancer research, but their identity needs to be confirmed in order to minimise the influence of cell contaminations and genetic drift on the analysis. Using both public and novel data, we demonstrate the use of RNA-sequencing data analysis for cell line authentication by examining the validity of COLO205, DLD1, HCT15, HCT116, HKE3, HT29 and RKO colorectal cancer cell lines. We successfully authenticate the studied cell lines and validate previous reports indicating that DLD1 and HCT15 are synonymous. We also show that the analysed HKE3 cells harbour an unexpected KRAS-G13D mutation and confirm that this cell line is a genuine KRAS dosage mutant, rather than a true isogenic derivative of HCT116 expressing only the wild type KRAS. This authentication method could be used to revisit the numerous cell line based RNA sequencing experiments available in public data repositories, analyse new experiments where whole genome sequencing is not available, as well as facilitate comparisons of data from different experiments, platforms and laboratories.

MAPK kinase signalling dynamics regulate cell fate decisions and drug resistance

The RAS/RAF/MEK/MAPK kinase pathway has been extensively studied for more than 25 years, yet we continue to be puzzled by its intricate dynamic control and plasticity. Different spatiotemporal MAPK dynamics bring about distinct cell fate decisions in normal vs cancer cells and developing organisms. Recent modelling and experimental studies provided novel insights in the versatile MAPK dynamics concerted by a plethora of feedforward/feedback regulations and crosstalk on multiple timescales. Multiple cancer types and various developmental disorders arise from persistent alterations of the MAPK dynamics caused by RAS/RAF/MEK mutations. While a key role of the MAPK pathway in multiple diseases made the development of novel RAF/MEK inhibitors a hot topic of drug development, these drugs have unexpected side-effects and resistance inevitably occurs. We review how RAF dimerization conveys drug resistance and recent breakthroughs to overcome this resistance.

Autophosphorylation on S614 inhibits the activity and the transforming potential of BRAF

The BRAF proto-oncogene serine/threonine-protein kinase, known as BRAF, belongs to the RAF kinase family. It regulates the MAPK/ERK signalling pathway affecting several cellular processes such as growth, survival, differentiation, and cellular transformation. BRAF is mutated in ~8% of all human cancers with the V600E mutation constituting ~90% of mutations. Here, we have used quantitative mass spectrometry to map and compare phosphorylation site patterns between BRAF and BRAF V600E. We identified sites that are shared as well as several quantitative differences in phosphorylation abundance. The highest difference is phosphorylation of S614 in the activation loop which is ~5fold enhanced in BRAF V600E. Mutation of S614 increases the kinase activity of both BRAF and BRAF V600E and the transforming ability of BRAF V600E. The phosphorylation of S614 is mitogen inducible and the result of autophosphorylation. These data suggest that phosphorylation at this site is inhibitory, and part of the physiological shut-down mechanism of BRAF signalling.

Substrate-Trapped Interactors of PHD3 and FIH Cluster in Distinct Signaling Pathways

Amino acid hydroxylation is a post-translational modification that regulates intra- and inter-molecular protein-protein interactions. The modifications are regulated by a family of 2-oxoglutarate- (2OG) dependent enzymes and, although the biochemistry is well understood, until now only a few substrates have been described for these enzymes. Using quantitative interaction proteomics, we screened for substrates of the proline hydroxylase PHD3 and the asparagine hydroxylase FIH, which regulate the HIF-mediated hypoxic response. We were able to identify hundreds of potential substrates. Enrichment analysis revealed that the potential substrates of both hydroxylases cluster in the same pathways but frequently modify different nodes of signaling networks. We confirm that two proteins identified in our screen, MAPK6 (Erk3) and RIPK4, are indeed hydroxylated in a FIH- or PHD3-dependent mechanism. We further determined that FIH-dependent hydroxylation regulates RIPK4-dependent Wnt signaling, and that PHD3-dependent hydroxylation of MAPK6 protects the protein from proteasomal degradation.

Phosphorylation of RAF Kinase Dimers Drives Conformational Changes that Facilitate Transactivation

RAF kinases are key players in the MAPK signaling pathway and are important targets for personalized cancer therapy. RAF dimerization is part of the physiological activation mechanism, together with phosphorylation, and is known to convey resistance to RAF inhibitors. Herein, molecular dynamics simulations are used to show that phosphorylation of a key N-terminal acidic (NtA) motif facilitates RAF dimerization by introducing several interprotomer salt bridges between the αC-helix and charged residues upstream of the NtA motif. Additionally, we show that the R-spine of RAF interacts with a conserved Trp residue in the vicinity of the NtA motif, connecting the active sites of two protomers and thereby modulating the cooperative interactions in the RAF dimer. Our findings provide a first structure-based mechanism for the auto-transactivation of RAF and could be generally applicable to other kinases, opening new pathways for overcoming dimerization-related drug resistance.

Abstract CN05-03: Personalized cancer diagnostics and therapeutics based on the computational modeling of signal transduction networks

Signalling pathways exert finely tuned control over cell fate decisions that ultimately determine the behaviour of cancer cells. It could therefore be expected that the dynamics of key pathway activation may contain prognostically relevant information over and above that which is contained in the static nature of traditional biomarkers. To investigate this hypothesis we have characterised the network architecture regulating JNK stress signalling in neuroblastoma cells and applied an experimentally calibrated and validated computational model of this network to extract prognostic information from neuroblastoma patient-specific simulations of JNK activation. Survival analysis based upon the dynamics of these simulations revealed that an inability to initiate switch-like JNK signalling in silico was significantly associated with poor overall survival for both MYCN amplified and non-amplified neuroblastoma patients. Furthermore, our analysis demonstrated that in order to extract prognostic information from a signalling pathway, deciphering the extant network structure is a vital consideration in model development. We also show that network based analysis can lead to the discovery of new therapeutic targets.

Citation Format: Walter Kolch, Dirk Fey, Melinda Halasz, Nora Rauch, Amaya Garcia Munoz, Ruth Pilkington, Boris N. Kholodenko, David R. Croucher, Sean P. Kennedy, Jordan F. Hastings, Frank Westermann, Daniel Dreidax, Matthias Fischer, David Duffy, Aleksandar Krstic, Thomas Schwarzl. Personalized cancer diagnostics and therapeutics based on the computational modeling of signal transduction networks. [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2015 Nov 5-9; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2015;14(12 Suppl 2):Abstract nr CN05-03.

Nonreactive spreading at high temperature: molten metals and oxides on molybdenum

The spontaneous spreading of small liquid metal (Cu, Ag, Au) and oxide drops on Mo substrates has been studied using a drop transfer setup combined with high-speed video. Under the experimental conditions used in this work, spreading occurs in the absence of interfacial reactions or ridging. The analysis of the spreading data indicates that dissipation at the triple junction (that can be described in terms of a triple-line friction) is playing a dominant role in the movement of the liquid front. This is due, in part, to the much stronger atomic interactions in high-temperature systems when compared to organic liquids. As a result of this analysis, a comprehensive view of spreading emerges in which the strength of the atomic interactions (solid-liquid, liquid-liquid) determines the relative roles of viscous impedance and dissipation at the triple junction in spreading kinetics.

Characterization of cyclic AMP-dependent phosphorylation of neuronal membrane proteins

We examined the patterns of cyclic AMP-dependent protein phosphorylation in membranes prepared from rat cortical synaptosomes following gel electrophoresis and autoradiography. We determined the optimum pH (6.2), time (20 s), Mg2+ concentration (10 mM) and cyclic AMP concentration (5 microM) for the reaction. We also found that the detergents Triton X-100 and gramicidin S enhanced cyclic AMP-dependent protein phosphorylation. Inhibitors of the Na+, K+ ATPase (ouabain, NaF, vanadate) enhanced protein phosphorylation. This effect occurred in the presence but not in the absence of detergent. The addition of purified bovine brain cyclic AMP-dependent protein kinase catalytic subunit enhanced membrane protein phosphorylation. The addition of homogeneous neural (bovine brain) and non-neural (bovine skeletal muscle) cyclic AMP-dependent protein kinase type II regulatory subunit partially inhibited protein phosphorylation. Both neural and non-neural regulatory subunits behaved similarly. In addition to cyclic AMP-dependent phosphorylation, the alpha-subunit of pyruvate dehydrogenase (Mr = 41,000) is phosphorylated in a cyclic AMP-independent fashion. We also examined the phosphorylation pattern of membranes prepared from rat heart and found that the number of acceptor substrates was much less than that from the nervous system.

The influence of portosystemic shunting on zinc and vitamin A metabolism in liver cirrhosis

To elucidate the possible role of portosystemic shunting on zinc and vitamin A deficiency which has been described in patients with cirrhosis of the liver, a study on 37 hospitalized patients with liver cirrhosis was performed. Patients with surgical portosystemic shunt were found to have a significantly lower levels of zinc, vitamin A and retinol-binding protein (RBP) than controls and patients with cirrhosis without shunt. Patients with portal hypertension--considered to have spontaneous shunting--also has lower levels than those without this symptom. A significant correlation between zinc and vitamin A and RBP levels, respectively, was found. Also an increased renal zinc output was demonstrated. An influence of portosystemic shunting on zinc deficiency and subsequent vitamin A deficiency by decreased RBP release is concluded. The importance of these metabolic disorders for clinical symptoms is discussed.