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.

Dynamic regulation of RAS and RAS signaling.

RAS proteins regulate most aspects of cellular physiology. They are mutated in 30% of human cancers and 4% of developmental disorders termed Rasopathies. They cycle between active GTP-bound and inactive GDP-bound states. When active, they can interact with a wide range of effectors that control fundamental biochemical and biological processes. Emerging evidence suggests that RAS proteins are not simple on/off switches but sophisticated information processing devices that compute cell fate decisions by integrating external and internal cues. A critical component of this compute function is the dynamic regulation of RAS activation and downstream signaling that allows RAS to produce a rich and nuanced spectrum of biological outputs. We discuss recent findings how the dynamics of RAS and its downstream signaling is regulated. Starting from the structural and biochemical properties of wild-type and mutant RAS proteins and their activation cycle, we examine higher molecular assemblies, effector interactions and downstream signaling outputs, all under the aspect of dynamic regulation. We also consider how computational and mathematical modeling approaches contribute to analyze and understand the pleiotropic functions of RAS in health and disease.

B-cell and complement signature in severe hidradenitis suppurativa that does not respond to adalimumab.

BACKGROUND: Hidradenitis suppurativa (HS) is a chronic inflammatory skin disorder with significant morbidity. The pathogenesis remains incompletely understood although immune dysregulation plays an important role. It is challenging to treat and approximately 50% of patients respond clinically to adalimumab, the only licensed treatment. OBJECTIVES: To examine differences between lesional and nonlesional HS skin at baseline using bulk RNA sequencing, and to compare the transcriptome in the skin before and after 12 weeks of treatment with adalimumab. To examine transcriptomic differences between adalimumab responders and nonresponders using Hidradenitis Suppurativa Clinical Response and the International Hidradenitis Suppurativa Severity Score System (IHS4); and to compare transcriptomic differences based on disease severity (Hurley stage and IHS4). METHODS: We completed bulk RNA sequencing on lesional and nonlesional skin samples of patients before and after 12 weeks of treatment with adalimumab. RESULTS: Baseline differentially expressed genes and pathways between lesional and nonlesional skin highlighted chemokines and antimicrobial peptides produced by keratinocytes; B-cell function; T-cell-receptor, interleukin-17 and nuclear factor-κB signalling; and T-helper-cell differentiation. Transcriptomic differences were identified in lesional skin at baseline, between subsequent responders and nonresponders. Patients with severe HS who did not respond to adalimumab had enriched complement and B-cell activation pathways at baseline. In addition, logistic regression identified CCL28 in baseline lesional HS skin as a potential biomarker of treatment response. CONCLUSIONS: This highlights the potential for targeting B-cell and complement pathways in HS treatment and the potential of stratifying patients at baseline to the most suitable treatment based on the skin transcriptome. CCL28 has not previously been identified in HS skin and has potential clinical relevance due to its antimicrobial function and homing of B and T cells at epithelial surfaces. Our results provide data to inform future translational and clinical studies on therapeutics in HS.

Giving space to cell signaling.

Spatial separation of component enzymes in a pathway creates microdomains and gradients of signaling activities. Imaging techniques and computational modeling help us to understand the complex relationships between signaling and spatial information. In this issue, Neves et al. (2008) report how the propagation of spatial information is controlled by the shape of hippocampal neurons and feedback and feed-forward motifs of the beta-adrenergic receptor signaling pathway.

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.

CRABP2 - A novel biomarker for high-risk endometrial cancer.

OBJECTIVE: Investigate the clinical and functional implications of elevated CRABP2 expression in endometrial cancer (EC) patients. METHODS: Patients were stratified into high and low CRABP2 expression groups using a decision tree classifier. Univariate and multivariate statistical analyses determined the prognostic and clinicopathological consequences of increased CRABP2 expression. A CRABP2 gene signature was generated using differential expression analysis, and analyzed using network-based approaches. The findings were validated in The Clinical Proteomic Tumor Analysis Consortium (CPTAC), a newly generated cohort of 120 endometrial tissues, and The Cancer Dependency Map (DepMap). RESULTS: 60 (11%) patients in TCGA had high CRABP2 expression, whilst 468 (89%) had low expression. High expression was associated with serous EC, reduced overall survival, advanced stage and grade. Downstream retinoic acid receptors (RARG and RARA) were correlated with CRABP2 expression and were associated with worse prognosis in serous EC. The CRABP2 gene signature was enriched for Polycomb target gene sets, and was regulated by ELP3 and BMP7. BMP7 expression was increased in the CRABP2-high group, was associated with worse prognosis, and CRISPR-Cas9 screens revealed correlations in its cell-fitness score with CRABP2 following gene knockout. The opposite was true for ELP3, suggesting opposing effects from both master regulators. CONCLUSIONS: CRABP2 expression is associated with poor prognosis and advanced EC. The expression of RARA and RARG correlates with CRABP2 and are associated with worse prognosis in advanced histological subtypes. Polycomb target gene sets and two master regulators, ELP3 and BMP7, were identified as functionally relevant mechanisms driving aberrant CRABP2 expression.

Signalling ballet in space and time.

Although we have amassed extensive catalogues of signalling network components, our understanding of the spatiotemporal control of emergent network structures has lagged behind. Dynamic behaviour is starting to be explored throughout the genome, but analysis of spatial behaviours is still confined to individual proteins. The challenge is to reveal how cells integrate temporal and spatial information to determine specific biological functions. Key findings are the discovery of molecular signalling machines such as Ras nanoclusters, spatial activity gradients and flexible network circuitries that involve transcriptional feedback. They reveal design principles of spatiotemporal organization that are crucial for network function and cell fate decisions.

Accurate prediction of kinase-substrate networks using knowledge graphs.

Phosphorylation of specific substrates by protein kinases is a key control mechanism for vital cell-fate decisions and other cellular processes. However, discovering specific kinase-substrate relationships is time-consuming and often rather serendipitous. Computational predictions alleviate these challenges, but the current approaches suffer from limitations like restricted kinome coverage and inaccuracy. They also typically utilise only local features without reflecting broader interaction context. To address these limitations, we have developed an alternative predictive model. It uses statistical relational learning on top of phosphorylation networks interpreted as knowledge graphs, a simple yet robust model for representing networked knowledge. Compared to a representative selection of six existing systems, our model has the highest kinome coverage and produces biologically valid high-confidence predictions not possible with the other tools. Specifically, we have experimentally validated predictions of previously unknown phosphorylations by the LATS1, AKT1, PKA and MST2 kinases in human. Thus, our tool is useful for focusing phosphoproteomic experiments, and facilitates the discovery of new phosphorylation reactions. Our model can be accessed publicly via an easy-to-use web interface (LinkPhinder).

Loss of RAF kinase inhibitor protein is involved in myelomonocytic differentiation and aggravates RAS-driven myeloid leukemogenesis.

RAS-signaling mutations induce the myelomonocytic differentiation and proliferation of hematopoietic stem and progenitor cells. Moreover, they are important players in the development of myeloid neoplasias. RAF kinase inhibitor protein (RKIP) is a negative regulator of RAS-signaling. As RKIP loss has recently been described in RAS-mutated myelomonocytic acute myeloid leukemia, we now aimed to analyze its role in myelomonocytic differentiation and RAS-driven leukemogenesis. Therefore, we initially analyzed RKIP expression during human and murine hematopoietic differentiation and observed that it is high in hematopoietic stem and progenitor cells and lymphoid cells but decreases in cells belonging to the myeloid lineage. By employing short hairpin RNA knockdown experiments in CD34+ umbilical cord blood cells and the undifferentiated acute myeloid leukemia cell line HL-60, we show that RKIP loss is indeed functionally involved in myelomonocytic lineage commitment and drives the myelomonocytic differentiation of hematopoietic stem and progenitor cells. These results could be confirmed in vivo, where Rkip deletion induced a myelomonocytic differentiation bias in mice by amplifying the effects of granulocyte macrophage-colony-stimulating factor. We further show that RKIP is of relevance for RAS-driven myelomonocytic leukemogenesis by demonstrating that Rkip deletion aggravates the development of a myeloproliferative disease in NrasG12D -mutated mice. Mechanistically, we demonstrate that RKIP loss increases the activity of the RAS-MAPK/ERK signaling module. Finally, we prove the clinical relevance of these findings by showing that RKIP loss is a frequent event in chronic myelomonocytic leukemia, and that it co-occurs with RAS-signaling mutations. Taken together, these data establish RKIP as novel player in RAS-driven myeloid leukemogenesis.

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.

Targeting MAPK Signaling in Cancer: Mechanisms of Drug Resistance and Sensitivity.

Mitogen-activated protein kinase (MAPK) pathways represent ubiquitous signal transduction pathways that regulate all aspects of life and are frequently altered in disease. Here, we focus on the role of MAPK pathways in modulating drug sensitivity and resistance in cancer. We briefly discuss new findings in the extracellular signaling-regulated kinase (ERK) pathway, but mainly focus on the mechanisms how stress activated MAPK pathways, such as p38 MAPK and the Jun N-terminal kinases (JNK), impact the response of cancer cells to chemotherapies and targeted therapies. In this context, we also discuss the role of metabolic and epigenetic aberrations and new therapeutic opportunities arising from these changes.

Targeting promiscuous heterodimerization overcomes innate resistance to ERBB2 dimerization inhibitors in breast cancer.

BACKGROUND: The oncogenic receptor tyrosine kinase (RTK) ERBB2 is known to dimerize with other EGFR family members, particularly ERBB3, through which it potently activates PI3K signalling. Antibody-mediated inhibition of this ERBB2/ERBB3/PI3K axis has been a cornerstone of treatment for ERBB2-amplified breast cancer patients for two decades. However, the lack of response and the rapid onset of relapse in many patients now question the assumption that the ERBB2/ERBB3 heterodimer is the sole relevant effector target of these therapies. METHODS: Through a systematic protein-protein interaction screen, we have identified and validated alternative RTKs that interact with ERBB2. Using quantitative readouts of signalling pathway activation and cell proliferation, we have examined their influence upon the mechanism of trastuzumab- and pertuzumab-mediated inhibition of cell growth in ERBB2-amplified breast cancer cell lines and a patient-derived xenograft model. RESULTS: We now demonstrate that inactivation of ERBB3/PI3K by these therapeutic antibodies is insufficient to inhibit the growth of ERBB2-amplified breast cancer cells. Instead, we show extensive promiscuity between ERBB2 and an array of RTKs from outside of the EGFR family. Paradoxically, pertuzumab also acts as an artificial ligand to promote ERBB2 activation and ERK signalling, through allosteric activation by a subset of these non-canonical RTKs. However, this unexpected activation mechanism also increases the sensitivity of the receptor network to the ERBB2 kinase inhibitor lapatinib, which in combination with pertuzumab, displays a synergistic effect in single-agent resistant cell lines and PDX models. CONCLUSIONS: The interaction of ERBB2 with a number of non-canonical RTKs activates a compensatory signalling response following treatment with pertuzumab, although a counter-intuitive combination of ERBB2 antibody therapy and a kinase inhibitor can overcome this innate therapeutic resistance.

Transcriptional and metabolic rewiring of colorectal cancer cells expressing the oncogenic KRASG13D mutation.

BACKGROUND: Activating mutations in KRAS frequently occur in colorectal cancer (CRC) patients, leading to resistance to EGFR-targeted therapies. METHODS: To better understand the cellular reprogramming which occurs in mutant KRAS cells, we have undertaken a systems-level analysis of four CRC cell lines which express either wild type (wt) KRAS or the oncogenic KRASG13D allele (mtKRAS). RESULTS: RNAseq revealed that genes involved in ribosome biogenesis, mRNA translation and metabolism were significantly upregulated in mtKRAS cells. Consistent with the transcriptional data, protein synthesis and cell proliferation were significantly higher in the mtKRAS cells. Targeted metabolomics analysis also confirmed the metabolic reprogramming in mtKRAS cells. Interestingly, mtKRAS cells were highly transcriptionally responsive to EGFR activation by TGFα stimulation, which was associated with an unexpected downregulation of genes involved in a range of anabolic processes. While TGFα treatment strongly activated protein synthesis in wtKRAS cells, protein synthesis was not activated above basal levels in the TGFα-treated mtKRAS cells. This was likely due to the defective activation of the mTORC1 and other pathways by TGFα in mtKRAS cells, which was associated with impaired activation of PKB signalling and a transient induction of AMPK signalling. CONCLUSIONS: We have found that mtKRAS cells are substantially rewired at the transcriptional, translational and metabolic levels and that this rewiring may reveal new vulnerabilities in oncogenic KRAS CRC cells that could be exploited in future.

Mutant K-Ras activation of the proapoptotic MST2 pathway is antagonized by wild-type K-Ras.

K-Ras mutations are frequent in colorectal cancer (CRC), albeit K-Ras is the only Ras isoform that can elicit apoptosis. Here, we show that mutant K-Ras directly binds to the tumor suppressor RASSF1A to activate the apoptotic MST2-LATS1 pathway. In this pathway LATS1 binds to and sequesters the ubiquitin ligase Mdm2 causing stabilization of the tumor suppressor p53 and apoptosis. However, mutant Ras also stimulates autocrine activation of the EGF receptor (EGFR) which counteracts mutant K-Ras-induced apoptosis. Interestingly, this protection requires the wild-type K-Ras allele, which inhibits the MST2 pathway in part via AKT activation. Confirming the pathophysiological relevance of the molecular findings, we find a negative correlation between K-Ras mutation and MST2 expression in human CRC patients and CRC mouse models. The small number of tumors with co-expression of mutant K-Ras and MST2 has elevated apoptosis rates. Thus, in CRC, mutant K-Ras transformation is supported by the wild-type allele.

Vascular Endothelial Growth Factor (VEGF) Promotes Assembly of the p130Cas Interactome to Drive Endothelial Chemotactic Signaling and Angiogenesis.

p130Cas is a polyvalent adapter protein essential for cardiovascular development, and with a key role in cell movement. In order to identify the pathways by which p130Cas exerts its biological functions in endothelial cells we mapped the p130Cas interactome and its dynamic changes in response to VEGF using high-resolution mass spectrometry and reconstruction of protein interaction (PPI) networks with the aid of multiple PPI databases. VEGF enriched the p130Cas interactome in proteins involved in actin cytoskeletal dynamics and cell movement, including actin-binding proteins, small GTPases and regulators or binders of GTPases. Detailed studies showed that p130Cas association of the GTPase-binding scaffold protein, IQGAP1, plays a key role in VEGF chemotactic signaling, endothelial polarization, VEGF-induced cell migration, and endothelial tube formation. These findings indicate a cardinal role for assembly of the p130Cas interactome in mediating the cell migratory response to VEGF in angiogenesis, and provide a basis for further studies of p130Cas in cell movement.

MST2-RASSF protein-protein interactions through SARAH domains.

The detailed, atomistic-level understanding of molecular signaling along the tumor-suppressive Hippo signaling pathway that controls tissue homeostasis by balancing cell proliferation and death through apoptosis is a promising avenue for the discovery of novel anticancer drug targets. The activation of kinases such as Mammalian STE20-Like Protein Kinases 1 and 2 (MST1 and MST2)-modulated through both homo- and heterodimerization (e.g. interactions with Ras association domain family, RASSF, enzymes)-is a key upstream event in this pathway and remains poorly understood. On the other hand, RASSFs (such as RASSF1A or RASSF5) act as important apoptosis activators and tumor suppressors, although their exact regulatory roles are also unclear. We present recent molecular studies of signaling along the Ras-RASSF-MST pathway, which controls growth and apoptosis in eukaryotic cells, including a variety of modern molecular modeling and simulation techniques. Using recently available structural information, we discuss the complex regulatory scenario according to which RASSFs perform dual signaling functions, either preventing or promoting MST2 activation, and thus control cell apoptosis. Here, we focus on recent studies highlighting the special role being played by the specific interactions between the helical Salvador/RASSF/Hippo (SARAH) domains of MST2 and RASSF1a or RASSF5 enzymes. These studies are crucial for integrating atomistic-level mechanistic information about the structures and conformational dynamics of interacting proteins, with information available on their system-level functions in cellular signaling.

SARAH Domain-Mediated MST2-RASSF Dimeric Interactions.

RASSF enzymes act as key apoptosis activators and tumor suppressors, being downregulated in many human cancers, although their exact regulatory roles remain unknown. A key downstream event in the RASSF pathway is the regulation of MST kinases, which are main effectors of RASSF-induced apoptosis. The regulation of MST1/2 includes both homo- and heterodimerization, mediated by helical SARAH domains, though the underlying molecular interaction mechanism is unclear. Here, we study the interactions between RASSF1A, RASSF5, and MST2 SARAH domains by using both atomistic molecular simulation techniques and experiments. We construct and study models of MST2 homodimers and MST2-RASSF SARAH heterodimers, and we identify the factors that control their high molecular stability. In addition, we also analyze both computationally and experimentally the interactions of MST2 SARAH domains with a series of synthetic peptides particularly designed to bind to it, and hope that our approach can be used to address some of the challenging problems in designing new anti-cancer drugs.

Comparison of different statistical approaches for urinary peptide biomarker detection in the context of coronary artery disease.

BACKGROUND: When combined with a clinical outcome variable, the size, complexity and nature of mass-spectrometry proteomics data impose great statistical challenges in the discovery of potential disease-associated biomarkers. The purpose of this study was thus to evaluate the effectiveness of different statistical methods applied for urinary proteomic biomarker discovery and different methods of classifier modelling in respect of the diagnosis of coronary artery disease in 197 study subjects and the prognostication of acute coronary syndromes in 368 study subjects. RESULTS: Computing the discovery sub-cohorts comprising [Formula: see text] of the study subjects based on the Wilcoxon rank sum test, t-score, cat-score, binary discriminant analysis and random forests provided largely different numbers (ranging from 2 to 398) of potential peptide biomarkers. Moreover, these biomarker patterns showed very little overlap limited to fragments of type I and III collagens as the common denominator. However, these differences in biomarker patterns did mostly not translate into significant differently performing diagnostic or prognostic classifiers modelled by support vector machine, diagonal discriminant analysis, linear discriminant analysis, binary discriminant analysis and random forest. This was even true when different biomarker patterns were combined into master-patterns. CONCLUSION: In conclusion, our study revealed a very considerable dependence of peptide biomarker discovery on statistical computing of urinary peptide profiles while the observed diagnostic and/or prognostic reliability of classifiers was widely independent of the modelling approach. This may however be due to the limited statistical power in classifier testing. Nonetheless, our study showed that urinary proteome analysis has the potential to provide valuable biomarkers for coronary artery disease mirroring especially alterations in the extracellular matrix. It further showed that for a comprehensive discovery of biomarkers and thus of pathological information, the results of different statistical methods may best be combined into a master pattern that then can be used for classifier modelling.

Splicing factor hnRNP A2 activates the Ras-MAPK-ERK pathway by controlling A-Raf splicing in hepatocellular carcinoma development.

In recent years, it has become clear that splicing factors play a direct role in cancer development. We showed previously that splicing factors SRSF1, SRSF6, and hnRNP A2/B1 are up-regulated in several cancers and can act as oncogenes when up-regulated. Here we examined the role of splicing factors hnRNP A1/A1b and hnRNP A2/B1 in hepatocellular carcinoma (HCC). We show that the splicing factors hnRNP A1 and hnRNP A2 are up-regulated in HCC tumors derived from inflammation-induced liver cancer mouse model. Overexpression of hnRNP A1 or hnRNP A2, but not the splicing isoform hnRNP B1, induced tumor formation of immortalized liver progenitor cells, while knockdown of these proteins inhibited anchorage-independent growth and tumor growth of human liver cancer cell lines. In addition, we found that cells overexpressing hnRNP A2 showed constitutive activation of the Ras-MAPK-ERK pathway. In contrast, knockdown of hnRNP A2 inhibited the Ras-MAPK-ERK pathway and prevented ERK1/2 activation by EGF. Moreover, we found that hnRNP A2 regulates the splicing of A-Raf, reducing the production of a short dominant-negative isoform of A-Raf and elevating the full-length A-Raf transcript. Taken together, our data suggest that hnRNP A2 up-regulation in HCC induces an alternative splicing switch that down-regulates a dominant-negative isoform of A-Raf, leading to activation of the Raf-MEK-ERK pathway and cellular transformation.