Oncogenic Mutations Rewire Signaling Pathways by Switching Protein Recruitment to Phosphotyrosine Sites.

Tyrosine phosphorylation regulates multi-layered signaling networks with broad implications in (patho)physiology, but high-throughput methods for functional annotation of phosphotyrosine sites are lacking. To decipher phosphotyrosine signaling directly in tissue samples, we developed a mass-spectrometry-based interaction proteomics approach. We measured the in vivo EGF-dependent signaling network in lung tissue quantifying >1,000 phosphotyrosine sites. To assign function to all EGF-regulated sites, we determined their recruited protein signaling complexes in lung tissue by interaction proteomics. We demonstrated how mutations near tyrosine residues introduce molecular switches that rewire cancer signaling networks, and we revealed oncogenic properties of such a lung cancer EGFR mutant. To demonstrate the scalability of the approach, we performed >1,000 phosphopeptide pulldowns and analyzed them by rapid mass spectrometric analysis, revealing tissue-specific differences in interactors. Our approach is a general strategy for functional annotation of phosphorylation sites in tissues, enabling in-depth mechanistic insights into oncogenic rewiring of signaling networks.

Reciprocal priming between receptor tyrosine kinases at recycling endosomes orchestrates cellular signalling outputs.

Integration of signalling downstream of individual receptor tyrosine kinases (RTKs) is crucial to fine-tune cellular homeostasis during development and in pathological conditions, including breast cancer. However, how signalling integration is regulated and whether the endocytic fate of single receptors controls such signalling integration remains poorly elucidated. Combining quantitative phosphoproteomics and targeted assays, we generated a detailed picture of recycling-dependent fibroblast growth factor (FGF) signalling in breast cancer cells, with a focus on distinct FGF receptors (FGFRs). We discovered reciprocal priming between FGFRs and epidermal growth factor (EGF) receptor (EGFR) that is coordinated at recycling endosomes. FGFR recycling ligands induce EGFR phosphorylation on threonine 693. This phosphorylation event alters both FGFR and EGFR trafficking and primes FGFR-mediated proliferation but not cell invasion. In turn, FGFR signalling primes EGF-mediated outputs via EGFR threonine 693 phosphorylation. This reciprocal priming between distinct families of RTKs from recycling endosomes exemplifies a novel signalling integration hub where recycling endosomes orchestrate cellular behaviour. Therefore, targeting reciprocal priming over individual receptors may improve personalized therapies in breast and other cancers.

Dynamic lineage priming is driven via direct enhancer regulation by ERK.

Central to understanding cellular behaviour in multi-cellular organisms is the question of how a cell exits one transcriptional state to adopt and eventually become committed to another. Fibroblast growth factor-extracellular signal-regulated kinase (FGF -ERK) signalling drives differentiation of mouse embryonic stem cells (ES cells) and pre-implantation embryos towards primitive endoderm, and inhibiting ERK supports ES cell self-renewal1. Paracrine FGF-ERK signalling induces heterogeneity, whereby cells reversibly progress from pluripotency towards primitive endoderm while retaining their capacity to re-enter self-renewal2. Here we find that ERK reversibly regulates transcription in ES cells by directly affecting enhancer activity without requiring a change in transcription factor binding. ERK triggers the reversible association and disassociation of RNA polymerase II and associated co-factors from genes and enhancers with the mediator component MED24 having an essential role in ERK-dependent transcriptional regulation. Though the binding of mediator components responds directly to signalling, the persistent binding of pluripotency factors to both induced and repressed genes marks them for activation and/or reactivation in response to fluctuations in ERK activity. Among the repressed genes are several core components of the pluripotency network that act to drive their own expression and maintain the ES cell state; if their binding is lost, the ability to reactivate transcription is compromised. Thus, as long as transcription factor occupancy is maintained, so is plasticity, enabling cells to distinguish between transient and sustained signals. If ERK signalling persists, pluripotency transcription factor levels are reduced by protein turnover and irreversible gene silencing and commitment can occur.

Elucidating Fibroblast Growth Factor-Induced Kinome Dynamics Using Targeted Mass Spectrometry and Dynamic Modeling.

Fibroblast growth factors (FGFs) are paracrine or endocrine signaling proteins that, activated by their ligands, elicit a wide range of health and disease-related processes, such as cell proliferation and the epithelial-to-mesenchymal transition. The detailed molecular pathway dynamics that coordinate these responses have remained to be determined. To elucidate these, we stimulated MCF-7 breast cancer cells with either FGF2, FGF3, FGF4, FGF10, or FGF19. Following activation of the receptor, we quantified the kinase activity dynamics of 44 kinases using a targeted mass spectrometry assay. Our system-wide kinase activity data, supplemented with (phospho)proteomics data, reveal ligand-dependent distinct pathway dynamics, elucidate the involvement of not earlier reported kinases such as MARK, and revise some of the pathway effects on biological outcomes. In addition, logic-based dynamic modeling of the kinome dynamics further verifies the biological goodness-of-fit of the predicted models and reveals BRAF-driven activation upon FGF2 treatment and ARAF-driven activation upon FGF4 treatment.

SubcellulaRVis: a web-based tool to simplify and visualise subcellular compartment enrichment.

Cells contain intracellular compartments, including membrane-bound organelles and the nucleus, and are surrounded by a plasma membrane. Proteins are localised to one or more of these cellular compartments; the correct localisation of proteins is crucial for their correct processing and function. Moreover, proteins and the cellular processes they partake in are regulated by relocalisation in response to various cellular stimuli. High-throughput 'omics experiments result in a list of proteins or genes of interest; one way in which their functional role can be understood is through the knowledge of their subcellular localisation, as deduced through statistical enrichment for Gene Ontology Cellular Component (GOCC) annotations or similar. We have designed a bioinformatics tool, named SubcellulaRVis, that compellingly visualises the results of GOCC enrichment for quick interpretation of the localisation of a group of proteins (rather than single proteins). We demonstrate that SubcellulaRVis precisely describes the subcellular localisation of gene lists whose locations have been previously ascertained. SubcellulaRVis can be accessed via the web (http://phenome.manchester.ac.uk/subcellular/) or as a stand-alone app (https://github.com/JoWatson2011/subcellularvis). SubcellulaRVis will be useful for experimental biologists with limited bioinformatics expertise who want to analyse data related to protein (re)localisation and location-specific modules within the intracellular protein network.

Uncovering SUMOylation dynamics during cell-cycle progression reveals FoxM1 as a key mitotic SUMO target protein.

Loss of small ubiquitin-like modification (SUMOylation) in mice causes genomic instability due to the missegregation of chromosomes. Currently, little is known about the identity of relevant SUMO target proteins that are involved in this process and about global SUMOylation dynamics during cell-cycle progression. We performed a large-scale quantitative proteomics screen to address this and identified 593 proteins to be SUMO-2 modified, including the Forkhead box transcription factor M1 (FoxM1), a key regulator of cell-cycle progression and chromosome segregation. SUMOylation of FoxM1 peaks during G2 and M phase, when FoxM1 transcriptional activity is required. We found that a SUMOylation-deficient FoxM1 mutant was less active compared to wild-type FoxM1, implying that SUMOylation of the protein enhances its transcriptional activity. Mechanistically, SUMOylation blocks the dimerization of FoxM1, thereby relieving FoxM1 autorepression. Cells deficient for FoxM1 SUMOylation showed increased levels of polyploidy. Our findings contribute to understanding the role of SUMOylation during cell-cycle progression.

Using Multilayer Heterogeneous Networks to Infer Functions of Phosphorylated Sites.

Mass spectrometry-based quantitative phosphoproteomics has become an essential approach in the study of cellular processes such as signaling. Commonly used methods to analyze phosphoproteomics datasets depend on generic, gene-centric annotations such as Gene Ontology terms, which do not account for the function of a protein in a particular phosphorylation state. Analysis of phosphoproteomics data is hampered by a lack of phosphorylated site-specific annotations. We propose a method that combines shotgun phosphoproteomics data, protein-protein interactions, and functional annotations into a heterogeneous multilayer network. Phosphorylation sites are associated to potential functions using a random walk on the heterogeneous network (RWHN) algorithm. We validated our approach against a model of the MAPK/ERK pathway and functional annotations from PhosphoSitePlus and were able to associate differentially regulated sites on the same proteins to their previously described specific functions. We further tested the algorithm on three previously published datasets and were able to reproduce their experimentally validated conclusions and to associate phosphorylation sites with known functions based on their regulatory patterns. Our approach provides a refinement of commonly used analysis methods and accurately predicts context-specific functions for sites with similar phosphorylation profiles.

Number and brightness analysis reveals that NCAM and FGF2 elicit different assembly and dynamics of FGFR1 in live cells.

Both fibroblast growth factor-2 (FGF2) and neural cell adhesion molecule (NCAM) trigger FGF receptor 1 (FGFR1) signaling; however, they induce remarkably distinct receptor trafficking and cellular responses. The molecular basis of such a dichotomy and the role of distinct types of ligand-receptor interaction remain elusive. Number of molecules and brightness (N&B) analysis revealed that FGF2 and NCAM promote different FGFR1 assembly and dynamics at the plasma membrane. NCAM stimulation elicits long-lasting cycles of short-lived FGFR1 monomers and multimers, a behavior that might reflect a rapid FGFR1 internalization and recycling. FGF2, instead, induces stable dimerization at the dose that stimulates cell proliferation. Reducing the occupancy of FGFR1 in response to low FGF2 doses causes a switch towards cyclically exposed and unstable receptor dimers, consistently with previously reported biphasic response to FGF2 and with the divergent signaling elicited by different ligand concentrations. Similar instability was observed upon altering the endocytic pathway. Thus, FGF2 and NCAM induce differential FGFR1 clustering at the cell surface, which might account for the distinct intracellular fate of the receptor and, hence, for the different signaling cascades and cellular responses.

Functional proteomics defines the molecular switch underlying FGF receptor trafficking and cellular outputs.

The stimulation of fibroblast growth factor receptors (FGFRs) with distinct FGF ligands generates specific cellular responses. However, the mechanisms underlying this paradigm have remained elusive. Here, we show that FGF-7 stimulation leads to FGFR2b degradation and, ultimately, cell proliferation, whereas FGF-10 promotes receptor recycling and cell migration. By combining mass-spectrometry-based quantitative proteomics with fluorescence microscopy and biochemical methods, we find that FGF-10 specifically induces the rapid phosphorylation of tyrosine (Y) 734 on FGFR2b, which leads to PI3K and SH3BP4 recruitment. This complex is crucial for FGFR2b recycling and responses, given that FGF-10 stimulation of either FGFR2b_Y734F mutant- or SH3BP4-depleted cells switches the receptor endocytic route to degradation, resulting in decreased breast cancer cell migration and the inhibition of epithelial branching in mouse lung explants. Altogether, these results identify an intriguing ligand-dependent mechanism for the control of receptor fate and cellular outputs that may explain the pathogenic role of deregulated FGFR2b, thus offering therapeutic opportunities.

Discovery of a Gatekeeper Residue in the C-Terminal Tail of the Extracellular Signal-Regulated Protein Kinase 5 (ERK5).

The extracellular signal-regulated protein kinase 5 (ERK5) is a non-redundant mitogen-activated protein kinase (MAPK) that exhibits a unique C-terminal extension which comprises distinct structural and functional properties. Here, we sought to elucidate the significance of phosphoacceptor sites in the C-terminal transactivation domain of ERK5. We have found that Thr732 acted as a functional gatekeeper residue controlling C-terminal-mediated nuclear translocation and transcriptional enhancement. Consistently, using a non-bias quantitative mass spectrometry approach, we demonstrated that phosphorylation at Thr732 conferred selectivity for binding interactions of ERK5 with proteins related to chromatin and RNA biology, whereas a number of metabolic regulators were associated with full-length wild type ERK5. Additionally, our proteomic analysis revealed that phosphorylation of the Ser730-Glu-Thr732-Pro motif could occur independently of dual phosphorylation at Thr218-Glu-Tyr220 in the activation loop. Collectively, our results firmly establish the significance of C-terminal phosphorylation in regulating ERK5 function. The post-translational modification of ERK5 on its C-terminal tail might be of particular relevance in cancer cells where ERK5 has be found to be hyperphosphoryated.

Cylindromatosis Tumor Suppressor Protein (CYLD) Deubiquitinase is Necessary for Proper Ubiquitination and Degradation of the Epidermal Growth Factor Receptor.

Cylindromatosis tumor suppressor protein (CYLD) is a deubiquitinase, best known as an essential negative regulator of the NFkB pathway. Previous studies have suggested an involvement of CYLD in epidermal growth factor (EGF)-dependent signal transduction as well, as it was found enriched within the tyrosine-phosphorylated complexes in cells stimulated with the growth factor. EGF receptor (EGFR) signaling participates in central cellular processes and its tight regulation, partly through ubiquitination cascades, is decisive for a balanced cellular homeostasis. Here, using a combination of mass spectrometry-based quantitative proteomic approaches with biochemical and immunofluorescence strategies, we demonstrate the involvement of CYLD in the regulation of the ubiquitination events triggered by EGF. Our data show that CYLD regulates the magnitude of ubiquitination of several major effectors of the EGFR pathway by assisting the recruitment of the ubiquitin ligase Cbl-b to the activated EGFR complex. Notably, CYLD facilitates the interaction of EGFR with Cbl-b through its Tyr15 phosphorylation in response to EGF, which leads to fine-tuning of the receptor's ubiquitination and subsequent degradation. This represents a previously uncharacterized strategy exerted by this deubiquitinase and tumors suppressor for the negative regulation of a tumorigenic signaling pathway.

Immature truncated O-glycophenotype of cancer directly induces oncogenic features.

Aberrant expression of immature truncated O-glycans is a characteristic feature observed on virtually all epithelial cancer cells, and a very high frequency is observed in early epithelial premalignant lesions that precede the development of adenocarcinomas. Expression of the truncated O-glycan structures Tn and sialyl-Tn is strongly associated with poor prognosis and overall low survival. The genetic and biosynthetic mechanisms leading to accumulation of truncated O-glycans are not fully understood and include mutation or dysregulation of glycosyltransferases involved in elongation of O-glycans, as well as relocation of glycosyltransferases controlling initiation of O-glycosylation from Golgi to endoplasmic reticulum. Truncated O-glycans have been proposed to play functional roles for cancer-cell invasiveness, but our understanding of the biological functions of aberrant glycosylation in cancer is still highly limited. Here, we used exome sequencing of most glycosyltransferases in a large series of primary and metastatic pancreatic cancers to rule out somatic mutations as a cause of expression of truncated O-glycans. Instead, we found hypermethylation of core 1 ß3-Gal-T-specific molecular chaperone, a key chaperone for O-glycan elongation, as the most prevalent cause. We next used gene editing to produce isogenic cell systems with and without homogenous truncated O-glycans that enabled, to our knowledge, the first polyomic and side-by-side evaluation of the cancer O-glycophenotype in an organotypic tissue model and in xenografts. The results strongly suggest that truncation of O-glycans directly induces oncogenic features of cell growth and invasion. The study provides support for targeting cancer-specific truncated O-glycans with immunotherapeutic measures.

Systems-wide analysis of BCR signalosomes and downstream phosphorylation and ubiquitylation.

B-cell receptor (BCR) signaling is essential for the development and function of B cells; however, the spectrum of proteins involved in BCR signaling is not fully known. Here we used quantitative mass spectrometry-based proteomics to monitor the dynamics of BCR signaling complexes (signalosomes) and to investigate the dynamics of downstream phosphorylation and ubiquitylation signaling. We identify most of the previously known components of BCR signaling, as well as many proteins that have not yet been implicated in this system. BCR activation leads to rapid tyrosine phosphorylation and ubiquitylation of the receptor-proximal signaling components, many of which are co-regulated by both the modifications. We illustrate the power of multilayered proteomic analyses for discovering novel BCR signaling components by demonstrating that BCR-induced phosphorylation of RAB7A at S72 prevents its association with effector proteins and with endo-lysosomal compartments. In addition, we show that BCL10 is modified by LUBAC-mediated linear ubiquitylation, and demonstrate an important function of LUBAC in BCR-induced NF-κB signaling. Our results offer a global and integrated view of BCR signaling, and the provided datasets can serve as a valuable resource for further understanding BCR signaling networks.

Recent findings and technological advances in phosphoproteomics for cells and tissues.

Site-specific phosphorylation is a fast and reversible covalent post-translational modification that is tightly regulated in cells. The cellular machinery of enzymes that write, erase and read these modifications (kinases, phosphatases and phospho-binding proteins) is frequently deregulated in different diseases, including cancer. Large-scale studies of phosphoproteins - termed phosphoproteomics - strongly rely on the use of high-performance mass spectrometric instrumentation. This powerful technology has been applied to study a great number of phosphorylation-based phenotypes. Nevertheless, many technical and biological challenges have to be overcome to identify biologically relevant phosphorylation sites in cells and tissues. This review describes different technological strategies to identify and quantify phosphorylation sites with high accuracy, without significant loss of analysis speed and reproducibility in tissues and cells. Moreover, computational tools for analysis, integration and biological interpretation of phosphorylation events are discussed.

Off-line high-pH reversed-phase fractionation for in-depth phosphoproteomics.

Protein phosphorylation is an important post-translational modification (PTM) involved in embryonic development, adult homeostasis, and disease. Over the past decade, several advances have been made in liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based technologies to identify thousands of phosphorylation sites. However, in-depth phosphoproteomics often require off-line enrichment and fractionation techniques. In this study, we provide a detailed analysis of the physicochemical characteristics of phosphopeptides, which have been fractionated by off-line high-pH chromatography (HpH) before subsequent titanium dioxide (TiO2) enrichment and LC-MS/MS analysis. Our results demonstrate that HpH is superior to standard strong-cation exchange (SCX) fractionation in the total number of phosphopeptides detected when analyzing the same number of fractions by identical LC-MS/MS gradients. From 14 HpH fractions, we routinely identified over 30,000 unique phosphopeptide variants, which is more than twice the number of that obtained from SCX fractionation. HpH chromatography displayed an exceptional ability to fractionate singly phosphorylated peptides, with minor benefits for doubly phosphorylated peptides over that with SCX. Further optimizations in the pooling and concatenation strategy increased the total number of multiphosphorylated peptides detected after HpH fractionation. In conclusion, we provide a basic framework and resource for performing in-depth phosphoproteome studies utilizing off-line basic reversed-phased fractionation. Raw data is available at ProteomeXchange (PXD001404).

TIMP-1 increases expression and phosphorylation of proteins associated with drug resistance in breast cancer cells.

Tissue inhibitor of metalloproteinase 1 (TIMP-1) is a protein with a potential biological role in drug resistance. To elucidate the unknown molecular mechanisms underlying the association between high TIMP-1 levels and increased chemotherapy resistance, we employed SILAC-based quantitative mass spectrometry to analyze global proteome and phosphoproteome differences of MCF-7 breast cancer cells expressing high or low levels of TIMP-1. In TIMP-1 high expressing cells, 312 proteins and 452 phosphorylation sites were up-regulated. Among these were the cancer drug targets topoisomerase 1, 2A, and 2B, which may explain the resistance phenotype to topoisomerase inhibitors that was observed in cells with high TIMP-1 levels. Pathway analysis showed an enrichment of proteins from functional categories such as apoptosis, cell cycle, DNA repair, transcription factors, drug targets and proteins associated with drug resistance or sensitivity, and drug transportation. The NetworKIN algorithm predicted the protein kinases CK2a, CDK1, PLK1, and ATM as likely candidates involved in the hyperphosphorylation of the topoisomerases. Up-regulation of protein and/or phosphorylation levels of topoisomerases in TIMP-1 high expressing cells may be part of the mechanisms by which TIMP-1 confers resistance to treatment with the widely used topoisomerase inhibitors in breast and colorectal cancer.

Functional selectivity of Receptor Tyrosine Kinases regulates distinct cellular outputs.

Functional selectivity refers to the activation of differential signalling and cellular outputs downstream of the same membrane-bound receptor when activated by two or more different ligands. Functional selectivity has been described and extensively studied for G-protein Coupled Receptors (GPCRs), leading to specific therapeutic options for dysregulated GPCRs functions. However, studies regarding the functional selectivity of Receptor Tyrosine Kinases (RTKs) remain sparse. Here, we will summarize recent data about RTK functional selectivity focusing on how the nature and the amount of RTK ligands and the crosstalk of RTKs with other membrane proteins regulate the specificity of RTK signalling. In addition, we will discuss how structural changes in RTKs upon ligand binding affects selective signalling pathways. Much remains to be known about the integration of different signals affecting RTK signalling specificity to orchestrate long-term cellular outcomes. Recent advancements in omics, specifically quantitative phosphoproteomics, and in systems biology methods to study, model and integrate different types of large-scale omics data have increased our ability to compare several signals affecting RTK functional selectivity in a global, system-wide fashion. We will discuss how such methods facilitate the exploration of important signalling hubs and enable data-driven predictions aiming at improving the efficacy of therapeutics for diseases like cancer, where redundant RTK signalling pathways often compromise treatment efficacy.

Oncogenic ERRB2 signals through the AP-1 transcription factor to control mesenchymal-like properties of oesophageal adenocarcinoma.

Oesophageal adenocarcinoma (OAC) is a deadly disease with poor survival statistics and few targeted therapies available. One of the most common molecular aberrations in OAC is amplification or activation of the gene encoding the receptor tyrosine kinase ERBB2, and ERBB2 is targeted in the clinic for this subset of patients. However, the downstream consequences of these ERBB2 activating events are not well understood. Here we used a combination of phosphoproteomics, open chromatin profiling and transcriptome analysis on cell line models and patient-derived datasets to interrogate the molecular pathways operating downstream from ERBB2. Integrated analysis of these data sets converge on a model where dysregulated ERBB2 signalling is mediated at the transcriptional level by the transcription factor AP-1. AP-1 in turn controls cell behaviour by acting on cohorts of genes that regulate cell migration and adhesion, features often associated with EMT. Our study therefore provides a valuable resource for the cancer cell signalling community and reveals novel molecular determinants underlying the dysregulated behaviour of OAC cells.

Fibroblast growth factor receptor signalling dysregulation and targeting in breast cancer.

Fibroblast Growth Factor Receptor (FGFR) signalling plays a critical role in breast embryonal development, tissue homeostasis, tumorigenesis and metastasis. FGFR, its numerous FGF ligands and signalling partners are often dysregulated in breast cancer progression and are one of the causes of resistance to treatment in breast cancer. Furthermore, FGFR signalling on epithelial cells is affected by signals from the breast microenvironment, therefore increasing the possibility of breast developmental abnormalities or cancer progression. Increasing our understanding of the multi-layered roles of the complex family of FGFRs, their ligands FGFs and their regulatory partners may offer novel treatment strategies for breast cancer patients, as a single agent or rational co-target, which will be explored in depth in this review.

Data integration and mechanistic modelling for breast cancer biology: Current state and future directions.

Breast cancer is one of the most common cancers threatening women worldwide. A limited number of available treatment options, frequent recurrence, and drug resistance exacerbate the prognosis of breast cancer patients. Thus, there is an urgent need for methods to investigate novel treatment options, while taking into account the vast molecular heterogeneity of breast cancer. Recent advances in molecular profiling technologies, including genomics, epigenomics, transcriptomics, proteomics and metabolomics data, enable approaching breast cancer biology at multiple levels of omics interaction networks. Systems biology approaches, including computational inference of 'big data' and mechanistic modelling of specific pathways, are emerging to identify potential novel combinations of breast cancer subtype signatures and more diverse targeted therapies.