phuEGO: A Network-Based Method to Reconstruct Active Signaling Pathways From Phosphoproteomics Datasets.

Signaling networks are critical for virtually all cell functions. Our current knowledge of cell signaling has been summarized in signaling pathway databases, which, while useful, are highly biased toward well-studied processes, and do not capture context specific network wiring or pathway cross-talk. Mass spectrometry-based phosphoproteomics data can provide a more unbiased view of active cell signaling processes in a given context, however, it suffers from low signal-to-noise ratio and poor reproducibility across experiments. While progress in methods to extract active signaling signatures from such data has been made, there are still limitations with respect to balancing bias and interpretability. Here we present phuEGO, which combines up-to-three-layer network propagation with ego network decomposition to provide small networks comprising active functional signaling modules. PhuEGO boosts the signal-to-noise ratio from global phosphoproteomics datasets, enriches the resulting networks for functional phosphosites and allows the improved comparison and integration across datasets. We applied phuEGO to five phosphoproteomics data sets from cell lines collected upon infection with SARS CoV2. PhuEGO was better able to identify common active functions across datasets and to point to a subnetwork enriched for known COVID-19 targets. Overall, phuEGO provides a flexible tool to the community for the improved functional interpretation of global phosphoproteomics datasets.

Integrated systems biology approach identifies gene targets for endothelial dysfunction.

Endothelial dysfunction (ED) is critical in the development and progression of cardiovascular (CV) disorders, yet effective therapeutic targets for ED remain elusive due to limited understanding of its underlying molecular mechanisms. To address this gap, we employed a systems biology approach to identify potential targets for ED. Our study combined multi omics data integration, with siRNA screening, high content imaging and network analysis to prioritise key ED genes and identify a pro- and anti-ED network. We found 26 genes that, upon silencing, exacerbated the ED phenotypes tested, and network propagation identified a pro-ED network enriched in functions associated with inflammatory responses. Conversely, 31 genes ameliorated ED phenotypes, pointing to potential ED targets, and the respective anti-ED network was enriched in hypoxia, angiogenesis and cancer-related processes. An independent screen with 17 drugs found general agreement with the trends from our siRNA screen and further highlighted DUSP1, IL6 and CCL2 as potential candidates for targeting ED. Overall, our results demonstrate the potential of integrated system biology approaches in discovering disease-specific candidate drug targets for endothelial dysfunction.

Large-scale phage-based screening reveals extensive pan-viral mimicry of host short linear motifs.

Viruses mimic host short linear motifs (SLiMs) to hijack and deregulate cellular functions. Studies of motif-mediated interactions therefore provide insight into virus-host dependencies, and reveal targets for therapeutic intervention. Here, we describe the pan-viral discovery of 1712 SLiM-based virus-host interactions using a phage peptidome tiling the intrinsically disordered protein regions of 229 RNA viruses. We find mimicry of host SLiMs to be a ubiquitous viral strategy, reveal novel host proteins hijacked by viruses, and identify cellular pathways frequently deregulated by viral motif mimicry. Using structural and biophysical analyses, we show that viral mimicry-based interactions have similar binding strength and bound conformations as endogenous interactions. Finally, we establish polyadenylate-binding protein 1 as a potential target for broad-spectrum antiviral agent development. Our platform enables rapid discovery of mechanisms of viral interference and the identification of potential therapeutic targets which can aid in combating future epidemics and pandemics.

Identification of phenotype-specific networks from paired gene expression-cell shape imaging data.

The morphology of breast cancer cells is often used as an indicator of tumor severity and prognosis. Additionally, morphology can be used to identify more fine-grained, molecular developments within a cancer cell, such as transcriptomic changes and signaling pathway activity. Delineating the interface between morphology and signaling is important to understand the mechanical cues that a cell processes in order to undergo epithelial-to-mesenchymal transition and consequently metastasize. However, the exact regulatory systems that define these changes remain poorly characterized. In this study, we used a network-systems approach to integrate imaging data and RNA-seq expression data. Our workflow allowed the discovery of unbiased and context-specific gene expression signatures and cell signaling subnetworks relevant to the regulation of cell shape, rather than focusing on the identification of previously known, but not always representative, pathways. By constructing a cell-shape signaling network from shape-correlated gene expression modules and their upstream regulators, we found central roles for developmental pathways such as WNT and Notch, as well as evidence for the fine control of NF-kB signaling by numerous kinase and transcriptional regulators. Further analysis of our network implicates a gene expression module enriched in the RAP1 signaling pathway as a mediator between the sensing of mechanical stimuli and regulation of NF-kB activity, with specific relevance to cell shape in breast cancer.

Dynamic regulation of hypoxia-inducible factor-1α activity is essential for normal B cell development.

B lymphocyte development and selection are central to adaptive immunity and self-tolerance. These processes require B cell receptor (BCR) signaling and occur in bone marrow, an environment with variable hypoxia, but whether hypoxia-inducible factor (HIF) is involved is unknown. We show that HIF activity is high in human and murine bone marrow pro-B and pre-B cells and decreases at the immature B cell stage. This stage-specific HIF suppression is required for normal B cell development because genetic activation of HIF-1α in murine B cells led to reduced repertoire diversity, decreased BCR editing and developmental arrest of immature B cells, resulting in reduced peripheral B cell numbers. HIF-1α activation lowered surface BCR, CD19 and B cell-activating factor receptor and increased expression of proapoptotic BIM. BIM deletion rescued the developmental block. Administration of a HIF activator in clinical use markedly reduced bone marrow and transitional B cells, which has therapeutic implications. Together, our work demonstrates that dynamic regulation of HIF-1α is essential for normal B cell development.

Prediction of Signed Protein Kinase Regulatory Circuits.

Complex networks of regulatory relationships between protein kinases comprise a major component of intracellular signaling. Although many kinase-kinase regulatory relationships have been described in detail, these tend to be limited to well-studied kinases whereas the majority of possible relationships remains unexplored. Here, we implement a data-driven, supervised machine learning method to predict human kinase-kinase regulatory relationships and whether they have activating or inhibiting effects. We incorporate high-throughput data, kinase specificity profiles, and structural information to produce our predictions. The results successfully recapitulate previously annotated regulatory relationships and can reconstruct known signaling pathways from the ground up. The full network of predictions is relatively sparse, with the vast majority of relationships assigned low probabilities. However, it nevertheless suggests denser modes of inter-kinase regulation than normally considered in intracellular signaling research. A record of this paper's transparent peer review process is included in the Supplemental Information.

Systems analysis of RhoGEF and RhoGAP regulatory proteins reveals spatially organized RAC1 signalling from integrin adhesions.

Rho GTPases are central regulators of the cytoskeleton and, in humans, are controlled by 145 multidomain guanine nucleotide exchange factors (RhoGEFs) and GTPase-activating proteins (RhoGAPs). How Rho signalling patterns are established in dynamic cell spaces to control cellular morphogenesis is unclear. Through a family-wide characterization of substrate specificities, interactomes and localization, we reveal at the systems level how RhoGEFs and RhoGAPs contextualize and spatiotemporally control Rho signalling. These proteins are widely autoinhibited to allow local regulation, form complexes to jointly coordinate their networks and provide positional information for signalling. RhoGAPs are more promiscuous than RhoGEFs to confine Rho activity gradients. Our resource enabled us to uncover a multi-RhoGEF complex downstream of G-protein-coupled receptors controlling CDC42-RHOA crosstalk. Moreover, we show that integrin adhesions spatially segregate GEFs and GAPs to shape RAC1 activity zones in response to mechanical cues. This mechanism controls the protrusion and contraction dynamics fundamental to cell motility. Our systems analysis of Rho regulators is key to revealing emergent organization principles of Rho signalling.

Proteomics and phosphoproteomics in precision medicine: applications and challenges.

Recent advances in proteomics allow the accurate measurement of abundances for thousands of proteins and phosphoproteins from multiple samples in parallel. Therefore, for the first time, we have the opportunity to measure the proteomic profiles of thousands of patient samples or disease model cell lines in a systematic way, to identify the precise underlying molecular mechanism and discover personalized biomarkers, networks and treatments. Here, we review examples of successful use of proteomics and phosphoproteomics data sets in as well as their integration other omics data sets with the aim of precision medicine. We will discuss the bioinformatics challenges posed by the generation, analysis and integration of such large data sets and present potential reasons why proteomics profiling and biomarkers are not currently widely used in the clinical setting. We will finally discuss ways to contribute to the better use of proteomics data in precision medicine and the clinical setting.

ATtRACT-a database of RNA-binding proteins and associated motifs.

RNA-binding proteins (RBPs) play a crucial role in key cellular processes, including RNA transport, splicing, polyadenylation and stability. Understanding the interaction between RBPs and RNA is key to improve our knowledge of RNA processing, localization and regulation in a global manner. Despite advances in recent years, a unified non-redundant resource that includes information on experimentally validated motifs, RBPs and integrated tools to exploit this information is lacking. Here, we developed a database named ATtRACT (available athttp://attract.cnic.es) that compiles information on 370 RBPs and 1583 RBP consensus binding motifs, 192 of which are not present in any other database. To populate ATtRACT we (i) extracted and hand-curated experimentally validated data from CISBP-RNA, SpliceAid-F, RBPDB databases, (ii) integrated and updated the unavailable ASD database and (iii) extracted information from Protein-RNA complexes present in Protein Data Bank database through computational analyses. ATtRACT provides also efficient algorithms to search a specific motif and scan one or more RNA sequences at a time. It also allows discoveringde novomotifs enriched in a set of related sequences and compare them with the motifs included in the database.Database URL:http:// attract. cnic. es.

A knowledge base for Vitis vinifera functional analysis.

BACKGROUND: Vitis vinifera (Grapevine) is the most important fruit species in the modern world. Wine and table grapes sales contribute significantly to the economy of major wine producing countries. The most relevant goals in wine production concern quality and safety. In order to significantly improve the achievement of these objectives and to gain biological knowledge about cultivars, a genomic approach is the most reliable strategy. The recent grapevine genome sequencing offers the opportunity to study the potential roles of genes and microRNAs in fruit maturation and other physiological and pathological processes. Although several systems allowing the analysis of plant genomes have been reported, none of them has been designed specifically for the functional analysis of grapevine genomes of cultivars under environmental stress in connection with microRNA data. DESCRIPTION: Here we introduce a novel knowledge base, called BIOWINE, designed for the functional analysis of Vitis vinifera genomes of cultivars present in Sicily. The system allows the analysis of RNA-seq experiments of two different cultivars, namely Nero d'Avola and Nerello Mascalese. Samples were taken under different climatic conditions of phenological phases, diseases, and geographic locations. The BIOWINE web interface is equipped with data analysis modules for grapevine genomes. In particular users may analyze the current genome assembly together with the RNA-seq data through a customized version of GBrowse. The web interface allows users to perform gene set enrichment by exploiting third-party databases. CONCLUSIONS: BIOWINE is a knowledge base implementing a set of bioinformatics tools for the analysis of grapevine genomes. The system aims to increase our understanding of the grapevine varieties and species of Sicilian products focusing on adaptability to different climatic conditions, phenological phases, diseases, and geographic locations.