Dynamic 3D proteomes reveal protein functional alterations at high resolution in situ.

Biological processes are regulated by intermolecular interactions and chemical modifications that do not affect protein levels, thus escaping detection in classical proteomic screens. We demonstrate here that a global protein structural readout based on limited proteolysis-mass spectrometry (LiP-MS) detects many such functional alterations, simultaneously and in situ, in bacteria undergoing nutrient adaptation and in yeast responding to acute stress. The structural readout, visualized as structural barcodes, captured enzyme activity changes, phosphorylation, protein aggregation, and complex formation, with the resolution of individual regulated functional sites such as binding and active sites. Comparison with prior knowledge, including other 'omics data, showed that LiP-MS detects many known functional alterations within well-studied pathways. It suggested distinct metabolite-protein interactions and enabled identification of a fructose-1,6-bisphosphate-based regulatory mechanism of glucose uptake in E. coli. The structural readout dramatically increases classical proteomics coverage, generates mechanistic hypotheses, and paves the way for in situ structural systems biology.

The Global Phosphorylation Landscape of SARS-CoV-2 Infection.

The causative agent of the coronavirus disease 2019 (COVID-19) pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has infected millions and killed hundreds of thousands of people worldwide, highlighting an urgent need to develop antiviral therapies. Here we present a quantitative mass spectrometry-based phosphoproteomics survey of SARS-CoV-2 infection in Vero E6 cells, revealing dramatic rewiring of phosphorylation on host and viral proteins. SARS-CoV-2 infection promoted casein kinase II (CK2) and p38 MAPK activation, production of diverse cytokines, and shutdown of mitotic kinases, resulting in cell cycle arrest. Infection also stimulated a marked induction of CK2-containing filopodial protrusions possessing budding viral particles. Eighty-seven drugs and compounds were identified by mapping global phosphorylation profiles to dysregulated kinases and pathways. We found pharmacologic inhibition of the p38, CK2, CDK, AXL, and PIKFYVE kinases to possess antiviral efficacy, representing potential COVID-19 therapies.

MetalinksDB: a flexible and contextualizable resource of metabolite-protein interactions.

From the catalytic breakdown of nutrients to signaling, interactions between metabolites and proteins play an essential role in cellular function. An important case is cell-cell communication, where metabolites, secreted into the microenvironment, initiate signaling cascades by binding to intra- or extracellular receptors of neighboring cells. Protein-protein cell-cell communication interactions are routinely predicted from transcriptomic data. However, inferring metabolite-mediated intercellular signaling remains challenging, partially due to the limited size of intercellular prior knowledge resources focused on metabolites. Here, we leverage knowledge-graph infrastructure to integrate generalistic metabolite-protein with curated metabolite-receptor resources to create MetalinksDB. MetalinksDB is an order of magnitude larger than existing metabolite-receptor resources and can be tailored to specific biological contexts, such as diseases, pathways, or tissue/cellular locations. We demonstrate MetalinksDB's utility in identifying deregulated processes in renal cancer using multi-omics bulk data. Furthermore, we infer metabolite-driven intercellular signaling in acute kidney injury using spatial transcriptomics data. MetalinksDB is a comprehensive and customizable database of intercellular metabolite-protein interactions, accessible via a web interface (https://metalinks.omnipathdb.org/) and programmatically as a knowledge graph (https://github.com/biocypher/metalinks). We anticipate that by enabling diverse analyses tailored to specific biological contexts, MetalinksDB will facilitate the discovery of disease-relevant metabolite-mediated intercellular signaling processes.

NADH Shuttling Couples Cytosolic Reductive Carboxylation of Glutamine with Glycolysis in Cells with Mitochondrial Dysfunction.

The bioenergetics and molecular determinants of the metabolic response to mitochondrial dysfunction are incompletely understood, in part due to a lack of appropriate isogenic cellular models of primary mitochondrial defects. Here, we capitalize on a recently developed cell model with defined levels of m.8993T>G mutation heteroplasmy, mTUNE, to investigate the metabolic underpinnings of mitochondrial dysfunction. We found that impaired utilization of reduced nicotinamide adenine dinucleotide (NADH) by the mitochondrial respiratory chain leads to cytosolic reductive carboxylation of glutamine as a new mechanism for cytosol-confined NADH recycling supported by malate dehydrogenase 1 (MDH1). We also observed that increased glycolysis in cells with mitochondrial dysfunction is associated with increased cell migration in an MDH1-dependent fashion. Our results describe a novel link between glycolysis and mitochondrial dysfunction mediated by reductive carboxylation of glutamine.

Metabolic Communication by SGLT2 Inhibition.

BACKGROUND: SGLT2 (sodium-glucose cotransporter 2) inhibitors (SGLT2i) can protect the kidneys and heart, but the underlying mechanism remains poorly understood. METHODS: To gain insights on primary effects of SGLT2i that are not confounded by pathophysiologic processes or are secondary to improvement by SGLT2i, we performed an in-depth proteomics, phosphoproteomics, and metabolomics analysis by integrating signatures from multiple metabolic organs and body fluids after 1 week of SGLT2i treatment of nondiabetic as well as diabetic mice with early and uncomplicated hyperglycemia. RESULTS: Kidneys of nondiabetic mice reacted most strongly to SGLT2i in terms of proteomic reconfiguration, including evidence for less early proximal tubule glucotoxicity and a broad downregulation of the apical uptake transport machinery (including sodium, glucose, urate, purine bases, and amino acids), supported by mouse and human SGLT2 interactome studies. SGLT2i affected heart and liver signaling, but more reactive organs included the white adipose tissue, showing more lipolysis, and, particularly, the gut microbiome, with a lower relative abundance of bacteria taxa capable of fermenting phenylalanine and tryptophan to cardiovascular uremic toxins, resulting in lower plasma levels of these compounds (including p-cresol sulfate). SGLT2i was detectable in murine stool samples and its addition to human stool microbiota fermentation recapitulated some murine microbiome findings, suggesting direct inhibition of fermentation of aromatic amino acids and tryptophan. In mice lacking SGLT2 and in patients with decompensated heart failure or diabetes, the SGLT2i likewise reduced circulating p-cresol sulfate, and p-cresol impaired contractility and rhythm in human induced pluripotent stem cell-derived engineered heart tissue. CONCLUSIONS: SGLT2i reduced microbiome formation of uremic toxins such as p-cresol sulfate and thereby their body exposure and need for renal detoxification, which, combined with direct kidney effects of SGLT2i, including less proximal tubule glucotoxicity and a broad downregulation of apical transporters (including sodium, amino acid, and urate uptake), provides a metabolic foundation for kidney and cardiovascular protection.

Pan-Cancer landscape of protein activities identifies drivers of signalling dysregulation and patient survival.

Genetic alterations in cancer cells trigger oncogenic transformation, a process largely mediated by the dysregulation of kinase and transcription factor (TF) activities. While the mutational profiles of thousands of tumours have been extensively characterised, the measurements of protein activities have been technically limited until recently. We compiled public data of matched genomics and (phospho)proteomics measurements for 1,110 tumours and 77 cell lines that we used to estimate activity changes in 218 kinases and 292 TFs. Co-regulation of kinase and TF activities reflects previously known regulatory relationships and allows us to dissect genetic drivers of signalling changes in cancer. We find that loss-of-function mutations are not often associated with the dysregulation of downstream targets, suggesting frequent compensatory mechanisms. Finally, we identified the activities most differentially regulated in cancer subtypes and showed how these can be linked to differences in patient survival. Our results provide broad insights into the dysregulation of protein activities in cancer and their contribution to disease severity.

A single-sample workflow for joint metabolomic and proteomic analysis of clinical specimens.

Understanding the interplay of the proteome and the metabolome helps to understand cellular regulation and response. To enable robust inferences from such multi-omics analyses, we introduced and evaluated a workflow for combined proteome and metabolome analysis starting from a single sample. Specifically, we integrated established and individually optimized protocols for metabolomic and proteomic profiling (EtOH/MTBE and autoSP3, respectively) into a unified workflow (termed MTBE-SP3), and took advantage of the fact that the protein residue of the metabolomic sample can be used as a direct input for proteome analysis. We particularly evaluated the performance of proteome analysis in MTBE-SP3, and demonstrated equivalence of proteome profiles irrespective of prior metabolite extraction. In addition, MTBE-SP3 combines the advantages of EtOH/MTBE and autoSP3 for semi-automated metabolite extraction and fully automated proteome sample preparation, respectively, thus advancing standardization and scalability for large-scale studies. We showed that MTBE-SP3 can be applied to various biological matrices (FFPE tissue, fresh-frozen tissue, plasma, serum and cells) to enable implementation in a variety of clinical settings. To demonstrate applicability, we applied MTBE-SP3 and autoSP3 to a lung adenocarcinoma cohort showing consistent proteomic alterations between tumour and non-tumour adjacent tissue independent of the method used. Integration with metabolomic data obtained from the same samples revealed mitochondrial dysfunction in tumour tissue through deregulation of OGDH, SDH family enzymes and PKM. In summary, MTBE-SP3 enables the facile and reliable parallel measurement of proteins and metabolites obtained from the same sample, benefiting from reduced sample variation and input amount. This workflow is particularly applicable for studies with limited sample availability and offers the potential to enhance the integration of metabolomic and proteomic datasets.

FUNKI: interactive functional footprint-based analysis of omics data.

MOTIVATION: Omics data are broadly used to get a snapshot of the molecular status of cells. In particular, changes in omics can be used to estimate the activity of pathways, transcription factors and kinases based on known regulated targets, that we call footprints. Then the molecular paths driving these activities can be estimated using causal reasoning on large signalling networks. RESULTS: We have developed FUNKI, a FUNctional toolKIt for footprint analysis. It provides a user-friendly interface for an easy and fast analysis of transcriptomics, phosphoproteomics and metabolomics data, either from bulk or single-cell experiments. FUNKI also features different options to visualize the results and run post-analyses, and is mirrored as a scripted version in R. AVAILABILITY AND IMPLEMENTATION: FUNKI is a free and open-source application built on R and Shiny, available at https://github.com/saezlab/ShinyFUNKI and https://saezlab.shinyapps.io/funki/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Causal integration of multi-omics data with prior knowledge to generate mechanistic hypotheses.

Multi-omics datasets can provide molecular insights beyond the sum of individual omics. Various tools have been recently developed to integrate such datasets, but there are limited strategies to systematically extract mechanistic hypotheses from them. Here, we present COSMOS (Causal Oriented Search of Multi-Omics Space), a method that integrates phosphoproteomics, transcriptomics, and metabolomics datasets. COSMOS combines extensive prior knowledge of signaling, metabolic, and gene regulatory networks with computational methods to estimate activities of transcription factors and kinases as well as network-level causal reasoning. COSMOS provides mechanistic hypotheses for experimental observations across multi-omics datasets. We applied COSMOS to a dataset comprising transcriptomics, phosphoproteomics, and metabolomics data from healthy and cancerous tissue from eleven clear cell renal cell carcinoma (ccRCC) patients. COSMOS was able to capture relevant crosstalks within and between multiple omics layers, such as known ccRCC drug targets. We expect that our freely available method will be broadly useful to extract mechanistic insights from multi-omics studies.

Increased CXCL4 expression in hematopoietic cells links inflammation and progression of bone marrow fibrosis in MPN.

Primary myelofibrosis (PMF) is a myeloproliferative neoplasm (MPN) that leads to progressive bone marrow (BM) fibrosis. Although the cellular mutations involved in the pathogenesis of PMF have been extensively investigated, the sequential events that drive stromal activation and fibrosis by hematopoietic-stromal cross-talk remain elusive. Using an unbiased approach and validation in patients with MPN, we determined that the differential spatial expression of the chemokine CXCL4/platelet factor-4 marks the progression of fibrosis. We show that the absence of hematopoietic CXCL4 ameliorates the MPN phenotype, reduces stromal cell activation and BM fibrosis, and decreases the activation of profibrotic pathways in megakaryocytes, inflammation in fibrosis-driving cells, and JAK/STAT activation in both megakaryocytes and stromal cells in 3 murine PMF models. Our data indicate that higher CXCL4 expression in MPN has profibrotic effects and is a mediator of the characteristic inflammation. Therefore, targeting CXCL4 might be a promising strategy to reduce inflammation in PMF.

Norepinephrine promotes triglyceride storage in macrophages via beta2-adrenergic receptor activation.

Tissue-resident macrophages are required for homeostasis, but also contribute to tissue dysfunction in pathophysiological states. The sympathetic neurotransmitter norepinephrine (NE) induces an anti-inflammatory and tissue-reparative phenotype in macrophages. As NE has a well-established role in promoting triglyceride lipolysis in adipocytes, and macrophages accumulate triglyceride droplets in various physiological and disease states, we investigated the effect of NE on primary mouse bone marrow-derived macrophage triglyceride metabolism. Surprisingly, our data show that in contrast to the canonical role of NE in stimulating lipolysis, NE acting via beta2-adrenergic receptors (B2ARs) in macrophages promotes extracellular fatty acid uptake and their storage as triglycerides and reduces free fatty acid release from triglyceride-laden macrophages. We demonstrate that these responses are mediated by a B2AR activation-dependent increase in Hilpda and Dgat1 gene expression and activity. We further show that B2AR activation favors the storage of extracellular polyunsaturated fatty acids. Finally, we present evidence that macrophages isolated from hearts after myocardial injury, for which survival critically depends on leukocyte B2ARs, have a transcriptional signature indicative of a transient triglyceride accumulation. Overall, we describe a novel and unexpected role of NE in promoting triglyceride storage in macrophages that could have potential implications in multiple diseases.

PHONEMeS: Efficient Modeling of Signaling Networks Derived from Large-Scale Mass Spectrometry Data.

Post-translational modifications of proteins play an important role in the regulation of cellular processes. The mass spectrometry analysis of proteome modifications offers huge potential for the study of how protein inhibitors affect the phosphosignaling mechanisms inside the cells. We have recently proposed PHONEMeS, a method that uses high-content shotgun phosphoproteomic data to build logical network models of signal perturbation flow. However, in its original implementation, PHONEMeS was computationally demanding and was only used to model signaling in a perturbation context. We have reformulated PHONEMeS as an Integer Linear Program (ILP) that is orders of magnitude more efficient than the original one. We have also expanded the scenarios that can be analyzed. PHONEMeS can model data upon perturbation on not only a known target but also deregulated pathways upstream and downstream of any set of deregulated kinases. Finally, PHONEMeS can now analyze data sets with multiple time points, which helps us to obtain better insight into the dynamics of the propagation of signals. We illustrate the value of the new approach on various data sets of medical relevance, where we shed light on signaling mechanisms and drug modes of action.

Converting networks to predictive logic models from perturbation signalling data with CellNOpt.

SUMMARY: The molecular changes induced by perturbations such as drugs and ligands are highly informative of the intracellular wiring. Our capacity to generate large datasets is increasing steadily. A useful way to extract mechanistic insight from the data is by integrating them with a prior knowledge network of signalling to obtain dynamic models. CellNOpt is a collection of Bioconductor R packages for building logic models from perturbation data and prior knowledge of signalling networks. We have recently developed new components and refined the existing ones to keep up with the computational demand of increasingly large datasets, including (i) an efficient integer linear programming, (ii) a probabilistic logic implementation for semi-quantitative datasets, (iii) the integration of a stochastic Boolean simulator, (iv) a tool to identify missing links, (v) systematic post-hoc analyses and (vi) an R-Shiny tool to run CellNOpt interactively. AVAILABILITY AND IMPLEMENTATION: R-package(s): https://github.com/saezlab/cellnopt. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Kinetic modelling of quantitative proteome data predicts metabolic reprogramming of liver cancer.

BACKGROUND: Metabolic alterations can serve as targets for diagnosis and cancer therapy. Due to the highly complex regulation of cellular metabolism, definite identification of metabolic pathway alterations remains challenging and requires sophisticated experimentation. METHODS: We applied a comprehensive kinetic model of the central carbon metabolism (CCM) to characterise metabolic reprogramming in murine liver cancer. RESULTS: We show that relative differences of protein abundances of metabolic enzymes obtained by mass spectrometry can be used to assess their maximal velocity values. Model simulations predicted tumour-specific alterations of various components of the CCM, a selected number of which were subsequently verified by in vitro and in vivo experiments. Furthermore, we demonstrate the ability of the kinetic model to identify metabolic pathways whose inhibition results in selective tumour cell killing. CONCLUSIONS: Our systems biology approach establishes that combining cellular experimentation with computer simulations of physiology-based metabolic models enables a comprehensive understanding of deregulated energetics in cancer. We propose that modelling proteomics data from human HCC with our approach will enable an individualised metabolic profiling of tumours and predictions of the efficacy of drug therapies targeting specific metabolic pathways.

MaBoSS 2.0: an environment for stochastic Boolean modeling.

MOTIVATION: Modeling of signaling pathways is an important step towards the understanding and the treatment of diseases such as cancers, HIV or auto-immune diseases. MaBoSS is a software that allows to simulate populations of cells and to model stochastically the intracellular mechanisms that are deregulated in diseases. MaBoSS provides an output of a Boolean model in the form of time-dependent probabilities, for all biological entities (genes, proteins, phenotypes, etc.) of the model. RESULTS: We present a new version of MaBoSS (2.0), including an updated version of the core software and an environment. With this environment, the needs for modeling signaling pathways are facilitated, including model construction, visualization, simulations of mutations, drug treatments and sensitivity analyses. It offers a framework for automated production of theoretical predictions. AVAILABILITY AND IMPLEMENTATION: MaBoSS software can be found at https://maboss.curie.fr , including tutorials on existing models and examples of models. CONTACT: gautier.stoll@upmc.fr or laurence.calzone@curie.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

LIANA+ provides an all-in-one framework for cell-cell communication inference.

The growing availability of single-cell and spatially resolved transcriptomics has led to the development of many approaches to infer cell-cell communication, each capturing only a partial view of the complex landscape of intercellular signalling. Here we present LIANA+, a scalable framework built around a rich knowledge base to decode coordinated inter- and intracellular signalling events from single- and multi-condition datasets in both single-cell and spatially resolved data. By extending and unifying established methodologies, LIANA+ provides a comprehensive set of synergistic components to study cell-cell communication via diverse molecular mediators, including those measured in multi-omics data. LIANA+ is accessible at https://github.com/saezlab/liana-py with extensive vignettes ( https://liana-py.readthedocs.io/ ) and provides an all-in-one solution to intercellular communication inference.

Multiomic profiling of breast cancer cells uncovers stress MAPK-associated sensitivity to AKT degradation.

More than 50% of human tumors display hyperactivation of the serine/threonine kinase AKT. Despite evidence of clinical efficacy, the therapeutic window of the current generation of AKT inhibitors could be improved. Here, we report the development of a second-generation AKT degrader, INY-05-040, which outperformed catalytic AKT inhibition with respect to cellular suppression of AKT-dependent phenotypes in breast cancer cell lines. A growth inhibition screen with 288 cancer cell lines confirmed that INY-05-040 had a substantially higher potency than our first-generation AKT degrader (INY-03-041), with both compounds outperforming catalytic AKT inhibition by GDC-0068. Using multiomic profiling and causal network integration in breast cancer cells, we demonstrated that the enhanced efficacy of INY-05-040 was associated with sustained suppression of AKT signaling, which was followed by induction of the stress mitogen-activated protein kinase (MAPK) c-Jun N-terminal kinase (JNK). Further integration of growth inhibition assays with publicly available transcriptomic, proteomic, and reverse phase protein array (RPPA) measurements established low basal JNK signaling as a biomarker for breast cancer sensitivity to AKT degradation. Together, our study presents a framework for mapping the network-wide signaling effects of therapeutically relevant compounds and identifies INY-05-040 as a potent pharmacological suppressor of AKT signaling.

Collaborative effect of Csnk1a1 haploinsufficiency and mutant p53 in Myc induction can promote leukemic transformation.

ABSTRACT: It is still not fully understood how genetic haploinsufficiency in del(5q) myelodysplastic syndrome (MDS) contributes to malignant transformation of hematopoietic stem cells. We asked how compound haploinsufficiency for Csnk1a1 and Egr1 in the common deleted region on chromosome 5 affects hematopoietic stem cells. Additionally, Trp53 was disrupted as the most frequently comutated gene in del(5q) MDS using CRISPR/Cas9 editing in hematopoietic progenitors of wild-type (WT), Csnk1a1-/+, Egr1-/+, Csnk1a1/Egr1-/+ mice. A transplantable acute leukemia only developed in the Csnk1a1-/+Trp53-edited recipient. Isolated blasts were indefinitely cultured ex vivo and gave rise to leukemia after transplantation, providing a tool to study disease mechanisms or perform drug screenings. In a small-scale drug screening, the collaborative effect of Csnk1a1 haploinsufficiency and Trp53 sensitized blasts to the CSNK1 inhibitor A51 relative to WT or Csnk1a1 haploinsufficient cells. In vivo, A51 treatment significantly reduced blast counts in Csnk1a1 haploinsufficient/Trp53 acute leukemias and restored hematopoiesis in the bone marrow. Transcriptomics on blasts and their normal counterparts showed that the derived leukemia was driven by MAPK and Myc upregulation downstream of Csnk1a1 haploinsufficiency cooperating with a downregulated p53 axis. A collaborative effect of Csnk1a1 haploinsufficiency and p53 loss on MAPK and Myc upregulation was confirmed on the protein level. Downregulation of Myc protein expression correlated with efficient elimination of blasts in A51 treatment. The "Myc signature" closely resembled the transcriptional profile of patients with del(5q) MDS with TP53 mutation.

Phosphoproteomic analysis of metformin signaling in colorectal cancer cells elucidates mechanism of action and potential therapeutic opportunities.

BACKGROUND: The biguanide drug metformin is a safe and widely prescribed drug for type 2 diabetes. Interestingly, hundreds of clinical trials have been set to evaluate the potential role of metformin in the prevention and treatment of cancer including colorectal cancer (CRC). However, the "metformin signaling" remains controversial. AIMS AND METHODS: To interrogate cell signaling induced by metformin in CRC and explore the druggability of the metformin-rewired phosphorylation network, we performed integrative analysis of phosphoproteomics, bioinformatics, and cell proliferation assays on a panel of 12 molecularly heterogeneous CRC cell lines. Using the high-resolute data-independent analysis mass spectrometry (DIA-MS), we monitored a total of 10,142 proteins and 56,080 phosphosites (P-sites) in CRC cells upon a short- and a long-term metformin treatment. RESULTS AND CONCLUSIONS: We found that metformin tended to primarily remodel cell signaling in the long-term and only minimally regulated the total proteome expression levels. Strikingly, the phosphorylation signaling response to metformin was highly heterogeneous in the CRC panel, based on a network analysis inferring kinase/phosphatase activities and cell signaling reconstruction. A "MetScore" was determined to assign the metformin relevance of each P-site, revealing new and robust phosphorylation nodes and pathways in metformin signaling. Finally, we leveraged the metformin P-site signature to identify pharmacodynamic interactions and confirmed a number of candidate metformin-interacting drugs, including navitoclax, a BCL-2/BCL-xL inhibitor. Together, we provide a comprehensive phosphoproteomic resource to explore the metformin-induced cell signaling for potential cancer therapeutics. This resource can be accessed at https://yslproteomics.shinyapps.io/Metformin/.

Mapping cardiac remodeling in chronic kidney disease.

Patients with advanced chronic kidney disease (CKD) mostly die from sudden cardiac death and recurrent heart failure. The mechanisms of cardiac remodeling are largely unclear. To dissect molecular and cellular mechanisms of cardiac remodeling in CKD in an unbiased fashion, we performed left ventricular single-nuclear RNA sequencing in two mouse models of CKD. Our data showed a hypertrophic response trajectory of cardiomyocytes with stress signaling and metabolic changes driven by soluble uremia-related factors. We mapped fibroblast to myofibroblast differentiation in this process and identified notable changes in the cardiac vasculature, suggesting inflammation and dysfunction. An integrated analysis of cardiac cellular responses to uremic toxins pointed toward endothelin-1 and methylglyoxal being involved in capillary dysfunction and TNFα driving cardiomyocyte hypertrophy in CKD, which was validated in vitro and in vivo. TNFα inhibition in vivo ameliorated the cardiac phenotype in CKD. Thus, interventional approaches directed against uremic toxins, such as TNFα, hold promise to ameliorate cardiac remodeling in CKD.