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.

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.

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.

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.

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.

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.

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/.

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.

Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches.

Introduction: The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing. Methods: Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors. Results: Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19. Discussion: The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

Dynamic partitioning of branched-chain amino acids-derived nitrogen supports renal cancer progression.

Metabolic reprogramming is critical for tumor initiation and progression. However, the exact impact of specific metabolic changes on cancer progression is poorly understood. Here, we integrate multimodal analyses of primary and metastatic clonally-related clear cell renal cancer cells (ccRCC) grown in physiological media to identify key stage-specific metabolic vulnerabilities. We show that a VHL loss-dependent reprogramming of branched-chain amino acid catabolism sustains the de novo biosynthesis of aspartate and arginine enabling tumor cells with the flexibility of partitioning the nitrogen of the amino acids depending on their needs. Importantly, we identify the epigenetic reactivation of argininosuccinate synthase (ASS1), a urea cycle enzyme suppressed in primary ccRCC, as a crucial event for metastatic renal cancer cells to acquire the capability to generate arginine, invade in vitro and metastasize in vivo. Overall, our study uncovers a mechanism of metabolic flexibility occurring during ccRCC progression, paving the way for the development of novel stage-specific therapies.

Comparison of methods and resources for cell-cell communication inference from single-cell RNA-Seq data.

The growing availability of single-cell data, especially transcriptomics, has sparked an increased interest in the inference of cell-cell communication. Many computational tools were developed for this purpose. Each of them consists of a resource of intercellular interactions prior knowledge and a method to predict potential cell-cell communication events. Yet the impact of the choice of resource and method on the resulting predictions is largely unknown. To shed light on this, we systematically compare 16 cell-cell communication inference resources and 7 methods, plus the consensus between the methods' predictions. Among the resources, we find few unique interactions, a varying degree of overlap, and an uneven coverage of specific pathways and tissue-enriched proteins. We then examine all possible combinations of methods and resources and show that both strongly influence the predicted intercellular interactions. Finally, we assess the agreement of cell-cell communication methods with spatial colocalisation, cytokine activities, and receptor protein abundance and find that predictions are generally coherent with those data modalities. To facilitate the use of the methods and resources described in this work, we provide LIANA, a LIgand-receptor ANalysis frAmework as an open-source interface to all the resources and methods.

Reducing lipid bilayer stress by monounsaturated fatty acids protects renal proximal tubules in diabetes.

In diabetic patients, dyslipidemia frequently contributes to organ damage such as diabetic kidney disease (DKD). Dyslipidemia is associated with both excessive deposition of triacylglycerol (TAG) in lipid droplets (LDs) and lipotoxicity. Yet, it is unclear how these two effects correlate with each other in the kidney and how they are influenced by dietary patterns. By using a diabetes mouse model, we find here that high-fat diet enriched in the monounsaturated oleic acid (OA) caused more lipid storage in LDs in renal proximal tubular cells (PTCs) but less tubular damage than a corresponding butter diet with the saturated palmitic acid (PA). This effect was particularly evident in S2/S3 but not S1 segments of the proximal tubule. Combining transcriptomics, lipidomics, and functional studies, we identify endoplasmic reticulum (ER) stress as the main cause of PA-induced PTC injury. Mechanistically, ER stress is caused by elevated levels of saturated TAG precursors, reduced LD formation, and, consequently, higher membrane order in the ER. Simultaneous addition of OA rescues the cytotoxic effects by normalizing membrane order and increasing both TAG and LD formation. Our study thus emphasizes the importance of monounsaturated fatty acids for the dietary management of DKD by preventing lipid bilayer stress in the ER and promoting TAG and LD formation in PTCs.

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.

Pre-analytical processing of plasma and serum samples for combined proteome and metabolome analysis.

Metabolomic and proteomic analyses of human plasma and serum samples harbor the power to advance our understanding of disease biology. Pre-analytical factors may contribute to variability and bias in the detection of analytes, especially when multiple labs are involved, caused by sample handling, processing time, and differing operating procedures. To better understand the impact of pre-analytical factors that are relevant to implementing a unified proteomic and metabolomic approach in a clinical setting, we assessed the influence of temperature, sitting times, and centrifugation speed on the plasma and serum metabolomes and proteomes from six healthy volunteers. We used targeted metabolic profiling (497 metabolites) and data-independent acquisition (DIA) proteomics (572 proteins) on the same samples generated with well-defined pre-analytical conditions to evaluate criteria for pre-analytical SOPs for plasma and serum samples. Time and temperature showed the strongest influence on the integrity of plasma and serum proteome and metabolome. While rapid handling and low temperatures (4°C) are imperative for metabolic profiling, the analyzed proteomics data set showed variability when exposed to temperatures of 4°C for more than 2 h, highlighting the need for compromises in a combined analysis. We formalized a quality control scoring system to objectively rate sample stability and tested this score using external data sets from other pre-analytical studies. Stringent and harmonized standard operating procedures (SOPs) are required for pre-analytical sample handling when combining proteomics and metabolomics of clinical samples to yield robust and interpretable data on a longitudinal scale and across different clinics. To ensure an adequate level of practicability in a clinical routine for metabolomics and proteomics studies, we suggest keeping blood samples up to 2 h on ice (4°C) prior to snap-freezing as a compromise between stability and operability. Finally, we provide the methodology as an open-source R package allowing the systematic scoring of proteomics and metabolomics data sets to assess the stability of plasma and serum samples.

decoupleR: ensemble of computational methods to infer biological activities from omics data.

Summary: Many methods allow us to extract biological activities from omics data using information from prior knowledge resources, reducing the dimensionality for increased statistical power and better interpretability. Here, we present decoupleR, a Bioconductor and Python package containing computational methods to extract these activities within a unified framework. decoupleR allows us to flexibly run any method with a given resource, including methods that leverage mode of regulation and weights of interactions, which are not present in other frameworks. Moreover, it leverages OmniPath, a meta-resource comprising over 100 databases of prior knowledge. Using decoupleR, we evaluated the performance of methods on transcriptomic and phospho-proteomic perturbation experiments. Our findings suggest that simple linear models and the consensus score across top methods perform better than other methods at predicting perturbed regulators. Availability and implementation: decoupleR's open-source code is available in Bioconductor (https://www.bioconductor.org/packages/release/bioc/html/decoupleR.html) for R and in GitHub (https://github.com/saezlab/decoupler-py) for Python. The code to reproduce the results is in GitHub (https://github.com/saezlab/decoupleR_manuscript) and the data in Zenodo (https://zenodo.org/record/5645208). Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Systems approach reveals distinct and shared signaling networks of the four PGE2 receptors in T cells.

Prostaglandin E2 (PGE2) promotes an immunosuppressive microenvironment in cancer, partly by signaling through four receptors (EP1, EP2, EP3, and EP4) on T cells. Here, we comprehensively characterized PGE2 signaling networks in helper, cytotoxic, and regulatory T cells using a phosphoproteomics and phosphoflow cytometry approach. We identified ~1500 PGE2-regulated phosphosites and several important EP1­4 signaling nodes, including PKC, CK2, PKA, PI3K, and Src. T cell subtypes exhibited distinct signaling pathways, with the strongest signaling in EP2-stimulated CD8+ cells. EP2 and EP4, both of which signal through Gαs, induced similar signaling outputs, but with distinct kinetics and intensity. Functional predictions from the observed phosphosite changes revealed PGE2 regulation of key cellular and immunological processes. Last, network modeling suggested signal integration between the receptors and a substantial contribution from G protein­independent signaling. This study offers a comprehensive view of the different PGE2-regulated phosphoproteomes in T cell subsets, providing a valuable resource for further research on this physiologically and pathophysiologically important signaling system.

SREBP1-induced fatty acid synthesis depletes macrophages antioxidant defences to promote their alternative activation.

Macrophages exhibit a spectrum of activation states ranging from classical to alternative activation1. Alternatively, activated macrophages are involved in diverse pathophysiological processes such as confining tissue parasites2, improving insulin sensitivity3 or promoting an immune-tolerant microenvironment that facilitates tumour growth and metastasis4. Recently, the metabolic regulation of macrophage function has come into focus as both the classical and alternative activation programmes require specific regulated metabolic reprogramming5. While most of the studies regarding immunometabolism have focussed on the catabolic pathways activated to provide energy, little is known about the anabolic pathways mediating macrophage alternative activation. In this study, we show that the anabolic transcription factor sterol regulatory element binding protein 1 (SREBP1) is activated in response to the canonical T helper 2 cell cytokine interleukin-4 to trigger the de novo lipogenesis (DNL) programme, as a necessary step for macrophage alternative activation. Mechanistically, DNL consumes NADPH, partitioning it away from cellular antioxidant defences and raising reactive oxygen species levels. Reactive oxygen species serves as a second messenger, signalling sufficient DNL, and promoting macrophage alternative activation. The pathophysiological relevance of this mechanism is validated by showing that SREBP1/DNL is essential for macrophage alternative activation in vivo in a helminth infection model.

Macrophage beta2-adrenergic receptor is dispensable for the adipose tissue inflammation and function.

OBJECTIVE: Neuroimmune interactions between the sympathetic nervous system (SNS) and macrophages are required for the homeostasis of multiple tissues, including the adipose tissue. It has been proposed that the SNS maintains adipose tissue macrophages (ATMs) in an anti-inflammatory state via direct norepinephrine (NE) signaling to macrophages. This study aimed to investigate the physiological importance of this paradigm by utilizing a mouse model in which the adrenergic signaling from the SNS to macrophages, but not to other adipose tissue cells, was disrupted. METHODS: We generated a macrophage-specific B2AR knockout mouse (Adrb2ΔLyz2) by crossing Adrb2fl/fl and Lyz2Cre/+ mice. We have previously shown that macrophages isolated from Adrb2ΔLyz2 animals do not respond to NE stimulation in vitro. Herein we performed a metabolic phenotyping of Adrb2ΔLyz2 mice on either chow or high-fat diet (HFD). We also assessed the adipose tissue function of Adrb2ΔLyz2 animals during fasting and cold exposure. Finally, we transplanted Adrb2ΔLyz2 bone marrow to low-density lipoprotein receptor (LDLR) knockout mice and investigated the development of atherosclerosis during Western diet feeding. RESULTS: We demonstrated that SNS-associated ATMs have a transcriptional profile indicative of activated beta-2 adrenergic receptor (B2AR), the main adrenergic receptor isoform in myeloid cells. However, Adrb2ΔLyz2 mice have unaltered energy balance on a chow or HFD. Furthermore, Adrb2ΔLyz2 mice show similar levels of adipose tissue inflammation and function during feeding, fasting, or cold exposure, and develop insulin resistance during HFD at the same rate as controls. Finally, macrophage-specific B2AR deletion does not affect the development of atherosclerosis on an LDL receptor-null genetic background. CONCLUSIONS: Overall, our data suggest that the SNS does not directly modulate the phenotype of adipose tissue macrophages in either lean mice or mouse models of cardiometabolic disease. Instead, sympathetic nerve activity exerts an indirect effect on adipose tissue macrophages through the modulation of adipocyte function.

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.

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.