HIPPIE v2.0: enhancing meaningfulness and reliability of protein-protein interaction networks

Bioinformatics-led identification of pathophysiological hallmark genes in diabesotension via graph clustering method

Background

Diabesotension, an overlapping triad of diabetes, hypertension, and obesity, remains a diagnostic challenge due to its complex underlying molecular mechanisms. Individuals with diabesotension face twice the risk of microvascular and macrovascular complications compared to those with either condition alone. However, the complexity of diabesotension poses significant diagnostic challenges due to limited knowledge of this disease trifecta.

Methods

The protein network was constructed, and the DPClusOST algorithm was applied to determine the protein clusters with a density ranging from 0.1 to 1.0 and those relevant to the pathophysiology of diabesotension. The significance score (SScore) was computed using the p-value from Fisher's exact test to evaluate each cluster, and the clusters containing proteins associated with diabesotension were classified using receiver operating characteristic (ROC) analysis. The significant density of the cluster, as indicated by the AUC, was determined and subsequently subjected to pathway enrichment analysis using ShinyGO.

Results

At densities of 0.6 and 0.8, 14 proteins (STX3, VAMP2, STX4, SYT1, DNAJC5, HSD17B10, DLD, AIFM1, PDHA1, PDHB, DLAT, PDHX, OGDH, and STAT5A) from clusters 13 and 53 were significantly identified as potential diabesotension-related proteins. Key pathways associated with the tripartite interplay of the three pathologies were found to involve amino acid metabolism, glycolysis/gluconeogenesis, SNARE-mediated vesicle transport, insulin and salivary secretion, and the glucagon and HIF-1 signaling pathways, thus identifying novel candidates for diabesotension biomarkers and therapeutic targets.

Conclusions

This study highlights the use of graph clustering to identify potential biomarkers for the comorbid triad, which could enhance personalized future treatment strategies.

Protein Sequence Analysis landscape: A Systematic Review of Task Types, Databases, Datasets, Word Embeddings Methods, and Language Models

Protein sequence analysis examines the order of amino acids within protein sequences to unlock diverse types of a wealth of knowledge about biological processes and genetic disorders. It helps in forecasting disease susceptibility by finding unique protein signatures, or biomarkers that are linked to particular disease states. Protein Sequence analysis through wet-lab experiments is expensive, time-consuming and error prone. To facilitate large-scale proteomics sequence analysis, the biological community is striving for utilizing AI competence for transitioning from wet-lab to computer aided applications. However, Proteomics and AI are two distinct fields and development of AI-driven protein sequence analysis applications requires knowledge of both domains. To bridge the gap between both fields, various review articles have been written. However, these articles focus revolves around few individual tasks or specific applications rather than providing a comprehensive overview about wide tasks and applications. Following the need of a comprehensive literature that presents a holistic view of wide array of tasks and applications, contributions of this manuscript are manifold: It bridges the gap between Proteomics and AI fields by presenting a comprehensive array of AI-driven applications for 63 distinct protein sequence analysis tasks. It equips AI researchers by facilitating biological foundations of 63 protein sequence analysis tasks. It enhances development of AI-driven protein sequence analysis applications by providing comprehensive details of 68 protein databases. It presents a rich data landscape, encompassing 627 benchmark datasets of 63 diverse protein sequence analysis tasks. It highlights the utilization of 25 unique word embedding methods and 13 language models in AI-driven protein sequence analysis applications. It accelerates the development of AI-driven applications by facilitating current state-of-the-art performances across 63 protein sequence analysis tasks.

Explainable Machine Learning Identifies Factors for Dosage Compensation in Aneuploid Human Cancer Cells

Aneuploidy, a hallmark of cancer, leads to widespread changes in chromosome copy number, altering the abundance of hundreds or thousands of proteins. How-ever, evidence suggests that levels of proteins encoded on affected chromosomes are often buffered toward their abundances observed in diploid cells. Despite its preval-ence, the molecular mechanisms driving this protein dosage compensation remain largely unknown. It is unclear whether all proteins are buffered to a similar degree, what factors determine buffering, and whether dosage compensation varies across different cell lines or tumor types. Moreover, its potential adaptive advantage and therapeutic relevance remain unexplored. Here, we established a novel approach to quantify protein dosage buffering in a gene copy number-dependent manner, show-ing that dosage compensation is widespread but variable in cancer cell lines and in vivo tumor samples. By developing multifactorial machine learning models, we identify mean gene dependency, protein complex participation, haploinsufficiency, and mRNA decay as key predictors of buffering. We also show that dosage com-pensation can affect oncogenic potential and that higher buffering correlates with reduced proteotoxic stress and increased drug resistance. These findings highlight protein dosage compensation as a crucial regulatory mechanism and a potential therapeutic target in aneuploid cancers.

Topology-driven negative sampling enhances generalizability in protein-protein interaction prediction

MOTIVATION

Unraveling the human interactome to uncover disease-specific patterns and discover drug targets hinges on accurate protein-protein interaction (PPI) predictions. However, challenges persist in machine learning (ML) models due to a scarcity of quality hard negative samples, shortcut learning, and limited generalizability to novel proteins.

RESULTS

In this study, we introduce a novel approach for strategic sampling of protein-protein noninteractions (PPNIs) by leveraging higher-order network characteristics that capture the inherent complementarity-driven mechanisms of PPIs. Next, we introduce Unsupervised Pre-training of Node Attributes tuned for PPI (UPNA-PPI), a high throughput sequence-to-function ML pipeline, integrating unsupervised pre-training in protein representation learning with Topological PPNI (TPPNI) samples, capable of efficiently screening billions of interactions. By using our TPPNI in training the UPNA-PPI model, we improve PPI prediction generalizability and interpretability, particularly in identifying potential binding sites locations on amino acid sequences, strengthening the prioritization of screening assays and facilitating the transferability of ML predictions across protein families and homodimers. UPNA-PPI establishes the foundation for a fundamental negative sampling methodology in graph machine learning by integrating insights from network topology.

AVAILABILITY AND IMPLEMENTATION

Code and UPNA-PPI predictions are freely available at https://github.com/alxndgb/UPNA-PPI.

Elucidation of short linear motif-based interactions of the MIT and rhodanese domains of the ubiquitin-specific protease 8

Ubiquitin-specific protease 8 (USP8) is a deubiquitinating enzyme with essential functions in protein trafficking and stability. It is a multidomain protein, with an N-terminal MIT (microtubule interacting and trafficking) domain, followed by a non-catalytic rhodanese (Rhod) domain, a long intrinsically disordered region, and a C-terminal catalytic domain. The N-terminal MIT domain of USP8 is known to mediate protein-protein interactions through binding to short linear motifs. The non-catalytic Rhod domain is also involved in protein-protein interactions, however detailed insights into these interactions remain limited. In this study we explore the short linear motif-based interactions of the MIT and Rhod domains of USP8 using a combination of proteomic peptide-phage display, peptide arrays and deep mutational scanning. We show that the MIT domain can bind ligands with a general [DE][LIF]x{2,3}R[FYIL]xxL[LV] consensus motif. We uncover that the rhodanese domain of USP8 is a peptide-binding domain, and define two distinct binding motifs (Rx[LI]xGxxxPxxL and G[LV][DE][IM]WExKxxxLxE) for this domain by deep mutational scanning of two different peptide ligands. Using the motif information, we predict binding sites within known USP8 interactors and substrates and validate interactions through peptide array analysis. Our findings demonstrate that both the USP8 MIT and rhodanese domains are peptide-binding domains that can be bound by degenerate and distinct binding motifs. The detailed information on the peptide binding preference of the two N-terminal domains of USP8 provide novel insights into the molecular recognition events that underlie the function of this essential deubiquitinating enzyme.

Metabolism-associated protein network constructing and host-directed anti-influenza drug repurposing

Host-directed antivirals offer a promising strategy for addressing the challenge of viral resistance. Virus-host interactions often trigger stage-specific metabolic reprogramming in the host, and the causal links between these interactions and virus-induced metabolic changes provide valuable insights for identifying host targets. In this study, we present a workflow for repurposing host-directed antivirals using virus-induced protein networks. These networks capture the dynamic progression of viral infection by integrating host proteins directly interacting with the virus and enzymes associated with significantly altered metabolic fluxes, identified through dual-species genome-scale metabolic models. This approach reveals numerous hub nodes as potential host targets. As a case study, 50 approved drugs with potential anti-influenza virus A (IVA) activity were identified through eight stage-specific IVA-induced protein networks, each comprising 699-899 hub nodes. Lisinopril, saxagliptin, and gliclazide were further validated for anti-IVA efficacy in vitro through assays measuring the inhibition of cytopathic effects and viral titers in A549 cells infected with IVA PR8. This workflow paves the way for the rapid repurposing of host-directed antivirals.

HHIP protein interactions in lung cells provide insight into COPD pathogenesis

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide. The primary causes of COPD are environmental, including cigarette smoking; however, genetic susceptibility also contributes to COPD risk. Genome-Wide Association Studies (GWASes) have revealed more than 80 genetic loci associated with COPD, leading to the identification of multiple COPD GWAS genes. However, the biological relationships between the identified COPD susceptibility genes are largely unknown. Genes associated with a complex disease are often in close network proximity, i.e. their protein products often interact directly with each other and/or similar proteins. In this study, we use affinity purification mass spectrometry (AP-MS) to identify protein interactions with HHIP, a well-established COPD GWAS gene which is part of the sonic hedgehog pathway, in two disease-relevant lung cell lines (IMR90 and 16HBE). To better understand the network neighborhood of HHIP, its proximity to the protein products of other COPD GWAS genes, and its functional role in COPD pathogenesis, we create HUBRIS, a protein-protein interaction network compiled from 8 publicly available databases. We identified both common and cell type-specific protein-protein interactors of HHIP. We find that our newly identified interactions shorten the network distance between HHIP and the protein products of several COPD GWAS genes, including DSP, MFAP2, TET2, and FBLN5. These new shorter paths include proteins that are encoded by genes involved in extracellular matrix and tissue organization. We found and validated interactions to proteins that provide new insights into COPD pathobiology, including CAVIN1 (IMR90) and TP53 (16HBE). The newly discovered HHIP interactions with CAVIN1 and TP53 implicate HHIP in response to oxidative stress.

Prediction and Evaluation of Coronavirus and Human Protein-Protein Interactions Integrating Five Different Computational Methods

The high lethality and infectiousness of coronaviruses, particularly SARS-Cov-2, pose a significant threat to human society. Understanding coronaviruses, especially the interactions between these viruses and humans, is crucial for mitigating the coronavirus pandemic. In this study, we conducted a comprehensive comparison and evaluation of five prevalent computational methods: interolog mapping, domain-domain interaction methodology, domain-motif interaction methodology, structure-based approaches, and machine learning techniques. These methods were assessed using unbiased datasets that include C1, C2h, C2v, and C3 test sets. Ultimately, we integrated these five methodologies into a unified model for predicting protein-protein interactions (PPIs) between coronaviruses and human proteins. Our final model demonstrates relatively better performance, particularly with the C2v and C3 test sets, which are frequently used datasets in practical applications. Based on this model, we further established a high-confidence PPI network between coronaviruses and humans, consisting of 18,012 interactions between 3843 human proteins and 129 coronavirus proteins. The reliability of our predictions was further validated through the current knowledge framework and network analysis. This study is anticipated to enhance mechanistic understanding of the coronavirus-human relationship a while facilitating the rediscovery of antiviral drug targets. The source codes and datasets are accessible at https://github.com/covhppilab/CoVHPPI.

Gene clusters linked to insulin resistance identified in a genome-wide study of the Taiwan Biobank population

This pioneering genome-wide association study examined surrogate markers for insulin resistance (IR) in 147,880 Taiwanese individuals using data from the Taiwan Biobank. The study focused on two IR surrogate markers: the triglyceride to high-density lipoprotein cholesterol (TG:HDL-C) ratio and the TyG index (the product of fasting plasma glucose and triglycerides). We identified genome-wide significance loci within four gene clusters: GCKR, MLXIPL, APOA5, and APOC1, uncovering 197 genes associated with IR. Transcriptome-wide association analysis revealed significant associations between these clusters and TyG, primarily in adipose tissue. Gene ontology analysis highlighted pathways related to Alzheimer's disease, glucose homeostasis, insulin resistance, and lipoprotein dynamics. The study identified sex-specific genes associated with TyG. Polygenic risk score analysis linked both IR markers to gout and hyperlipidemia. Our findings elucidate the complex relationships between IR surrogate markers, genetic predisposition, and disease phenotypes in the Taiwanese population, contributing valuable insights to the field of metabolic research.

In-silico proteomic analysis of the role of IL-4 and IL-10 in IVD degeneration: Protein-protein interaction networks for candidate prioritisation

Protein-protein interaction (PPI) networks provide a static map of functional protein interactions, which when combined with algorithms, can prioritize key protein candidates which experimental studies cannot capture. This study, aimed to construct knowledge-based nucleus pulposus (NP)-specific PPI networks which could be deployed to investigate complex protein interactions in human NP cells and tissues following IL-4 and IL-10 stimulation. NP-specific PPI networks were developed based on mass spectrometry (MS) and secretome datasets from human NP cells. These networks were validated using in vitro and ex vivo experimental data sets. Genes Underlying Inheritance Linked Disorders (GUILD) genome-wide network-based prioritization framework was employed for protein candidate prediction under no treatment baseline and IL-4, IL-10 and IL-1β single or combined stimulating scenarios. These secretome-based in vitro PPI networks were able to reproduce the no-treatment candidate prioritization baseline. Whereby within NP cells from discs isolated due to traumatic injury biglycan was identified whilst in degenerate samples decorin was highlighted. Furthermore, experimentally observed IL-4 pleiotropic behaviour was predicted by IL-1 receptor-like 1 prioritization. PPI network-based IL-4 and IL-10 conditions offered novel insights of potential candidates, including collagen IV and fibroblast growth factor intracellular binding protein (FIBP) as key candidates within IL-4 activation pathways, whereas urocortin 3 and neural growth factor were identified following IL-10 stimulation. Additionally, MS based PPI network propagation offered a more extensive, module-based structure networks with lower edge degree and biological variability. Overall, multiple proteomic experimental approaches are required to successfully validate in-silico prediction models to understand the complex interactions between the plethora of proteins involved in IVD degeneration.

The Potential of Peptide-Based Inhibitors in Disrupting Protein-Protein Interactions for Targeted Cancer Therapy

Protein-protein interactions (PPIs) form an intricate cellular network known as the interactome, which is essential for various cellular processes, such as gene regulation, signal transduction, and metabolic pathways. The dysregulation of this network has been closely linked to various disease states. In cancer, these aberrant PPIs, termed oncogenic PPIs (OncoPPIs), are involved in tumour formation and proliferation. Therefore, the inhibition of OncoPPIs becomes a strategy for targeted cancer therapy. Small molecule inhibitors have been the dominant strategy for PPI inhibition owing to their small size and ability to cross cell membranes. However, peptide-based inhibitors have emerged as compelling alternatives, offering distinct advantages over small molecule inhibitors. Peptides, with their larger size and flexible backbones, can effectively engage with the broad interfaces of PPIs. Their high specificity, lower toxicity, and ease of modification make them promising candidates for targeted cancer therapy. Over the past decade, significant advancements have been made in developing peptide-based inhibitors. This review discusses the critical aspects of targeting PPIs, emphasizes the significance of OncoPPIs in cancer therapy, and explores the advantages of using peptide-based inhibitors as therapeutic agents. It also highlights recent progress in peptide design aimed at overcoming the limitations of peptide therapeutics, offering a comprehensive overview of the current landscape and potential of peptide-based inhibitors in cancer treatment.

A Map of the Lipid-Metabolite-Protein Network to Aid Multi-Omics Integration

The integration of multi-omics data offers transformative potential for elucidating complex molecular mechanisms underlying biological processes and diseases. In this study, we developed a lipid-metabolite-protein network that combines a protein-protein interaction network and enzymatic and genetic interactions of proteins with metabolites and lipids to provide a unified framework for multi-omics integration. Using hyperbolic embedding, the network visualizes connections across omics layers, accessible through a user-friendly Shiny R (version 1.10.0) software package. This framework ranks molecules across omics layers based on functional proximity, enabling intuitive exploration. Application in a cardiovascular disease (CVD) case study identified lipids and metabolites associated with CVD-related proteins. The analysis confirmed known associations, like cholesterol esters and sphingomyelin, and highlighted potential novel biomarkers, such as 4-imidazoleacetate and indoleacetaldehyde. Furthermore, we used the network to analyze empagliflozin's temporal effects on lipid metabolism. Functional enrichment analysis of proteins associated with lipid signatures revealed dynamic shifts in biological processes, with early effects impacting phospholipid metabolism and long-term effects affecting sphingolipid biosynthesis. Our framework offers a versatile tool for hypothesis generation, functional analysis, and biomarker discovery. By bridging molecular layers, this approach advances our understanding of disease mechanisms and therapeutic effects, with broad applications in computational biology and precision medicine.

Prediction of protein interactions with function in protein (de-)phosphorylation

Protein-protein interactions (PPIs) form a complex network called "interactome" that regulates many functions in the cell. In recent years, there is an increasing accumulation of evidence supporting the existence of a hyperbolic geometry underlying the network representation of complex systems such as the interactome. In particular, it has been shown that the embedding of the human Protein-Interaction Network (hPIN) in hyperbolic space (H2) captures biologically relevant information. Here we explore whether this mapping contains information that would allow us to predict the function of PPIs, more specifically interactions related to post-translational modification (PTM). We used a random forest algorithm to predict PTM-related directed PPIs, concretely, protein phosphorylation and dephosphorylation, based on hyperbolic properties and centrality measures of the hPIN mapped in H2. To evaluate the efficacy of our algorithm, we predicted PTM-related PPIs of ataxin-1, a protein which is responsible for Spinocerebellar Ataxia type 1 (SCA1). Proteomics analysis in a cellular model revealed that several of the predicted PTM-PPIs were indeed dysregulated in a SCA1-related disease network. A compact cluster composed of ataxin-1, its dysregulated PTM-PPIs and their common upstream regulators may represent critical interactions for disease pathology. Thus, our algorithm may infer phosphorylation activity on proteins through directed PPIs.

Circulating causal protein networks linked to future risk of myocardial infarction

Variations in blood protein levels have been associated with a broad spectrum of complex diseases, including atherosclerotic cardiovascular disease (ACVD). These associations highlight the intricate interplay between local (e.g., cardiovascular) and systemic (non-cardiovascular) factors for the development of ACVD, emphasizing the need for a comprehensive, systems-level understanding of its etiology. To accomplish this, we developed a causal network inference framework by analyzing one of the largest serum proteomics studies to date, the Age, Gene/Environment Susceptibility-Reykjavik Study (AGES), a prospective population-based study of 7,523 serum proteins measured in 5,376 older adults. To reconstruct a causal network of serum proteins, we used cis -acting protein quantitative trait loci (pQTLs) as instrumental variables to infer causal relationships between protein pairs, while accounting for potential unobserved confounding factors. We identified 185 causal protein subnetworks (FDR = 1%, n ≥ 10 members), which collectively interacted with 5,611 target proteins, offering valuable biological insights and an overview of systemic homeostasis. Several subnetworks, many of which interact to establish a hierarchy of directional relationships, were significantly associated with future myocardial infarction and/or its long-term complications like heart failure, as well as with key cardiometabolic traits that contribute to the onset of ACVD.

Comparative proteomic landscapes elucidate human preimplantation development and failure

Understanding mammalian preimplantation development, particularly in humans, at the proteomic level remains limited. Here, we applied our comprehensive solution of ultrasensitive proteomic technology to measure the proteomic profiles of oocytes and early embryos and identified nearly 8,000 proteins in humans and over 6,300 proteins in mice. We observed distinct proteomic dynamics before and around zygotic genome activation (ZGA) between the two species. Integrative analysis with translatomic data revealed extensive divergence between translation activation and protein accumulation. Multi-omic analysis indicated that ZGA transcripts often contribute to protein accumulation in blastocysts. Using mouse embryos, we identified several transcriptional regulators critical for early development, thereby linking ZGA to the first lineage specification. Furthermore, single-embryo proteomics of poor-quality embryos from over 100 patient couples provided insights into preimplantation development failure. Our study may contribute to reshaping the framework of mammalian preimplantation development and opening avenues for addressing human infertility.

An Integrative Multitiered Computational Analysis for Better Understanding the Structure and Function of 85 Miniproteins

Background

Miniproteins, defined as polypeptides containing fewer than 50 amino acids, have recently elicited significant interest due to an emerging understanding of their diverse roles in fundamental biological processes. In addition, miniprotein dysregulation underlies human diseases and is a significant focus for biotechnology and drug development. Notably, the human genome project revealed the existence of many novel miniproteins, most of which remain uncharacterized. This study reports an approach for analyzing and scoring previously uncharacterized miniproteins by integrating knowledge from classic sequence-based bioinformatics, computational biophysics, and system biology annotations. Our results demonstrate that these approaches provide novel information on the structure-function relationship of these molecules with a particular focus on their biomedical relevance.

Methods

We identified 85 human miniproteins using a simple multi-tier approach. First, we performed a sequence-based analysis of these proteins using several algorithms to identify regions of structural and functional importance. Protein-protein interactions and gene ontology annotations were used to analyze miniprotein function. Then, we predicted miniprotein three-dimensional structures using AI-based methods and peptide modeling to determine their relative yields for these understudied polymers. Subsequently, we used several computational biophysics methods and structure-based calculations to annotate and evaluate results from both algorithms.

Results

We find several relations between predicted structure and functional properties to assign these proteins to several groups with similar properties. Sequence-based analysis leads us to identify motifs and residues that link structure-to-function for most of these proteins. We suggest novel miniprotein functions, such as thymosin beta proteins regulating the shelterin complex through TERF1 and POT1 interactions, FAM86JP and FAM66E participating in endocytic processes, and BAGE1 influencing chromatin remodeling through interaction with nuclear proteins. Further, known functions of miniproteins, such as STRIT1, STMP1, and SLN, were supported. Finally, structure-based scoring led us to build 3D models that provided complementary information to ontologies. We identify that structural propensity is not strictly dependent on polymer length. In fact, in this dataset, peptide-based algorithms may have advantages over AI-based algorithms for certain groups of miniproteins.

Conclusion

This analytic approach and resulting identification and annotation of miniproteins adds much to what is currently known about miniproteins. Our determination of novel properties of miniproteins bears significant mechanistic and biomedical relevance. We propose novel functions of miniproteins, which expands our understanding of their potential roles in cellular processes. And, we practically identify which sequence and structure-based tools provide the most information, aiding future studies of miniproteins.

SubCELL: the landscape of subcellular compartment-specific molecular interactions

The subcellular compartment-specific molecular interactions (SCSIs) are the building blocks for most molecular functions, biological processes and disease pathogeneses. Extensive experiments have therefore been conducted to accumulate the valuable information of SCSIs, but none of the available databases has been constructed to describe those data. In this study, a novel knowledge base SubCELL is thus introduced to depict the landscape of SCSIs among DNAs/RNAs/proteins. This database is UNIQUE in (a) providing, for the first time, the experimentally-identified SCSIs, (b) systematically illustrating a large number of SCSIs inferred based on well-established method and (c) collecting experimentally-determined subcellular locations for the DNAs/RNAs/proteins of diverse species. Given the essential physiological/pathological role of SCSIs, the SubCELL is highly expected to have great implications for modern molecular biological study, which can be freely accessed with no login requirement at: https://idrblab.org/subcell/.

Decoding the Therapeutic Target SVEP1: Harnessing Molecular Trait GWASs to Unravel Mechanisms of Human Disease

Although human genetics has substantial potential to illuminate novel disease pathways and facilitate drug development, identifying causal variants and deciphering their mechanisms remain challenging. We believe these challenges can be addressed, in part, by creatively repurposing the results of molecular trait genome-wide association studies (GWASs). In this review, we introduce techniques related to molecular GWASs and unconventionally apply them to understanding SVEP1, a human coronary artery disease risk locus. Our analyses highlight SVEP1's causal link to cardiometabolic disease and glaucoma, as well as the surprising discovery of SVEP1 as the first known physiologic ligand for PEAR1, a critical receptor governing platelet reactivity. We further employ these techniques to dissect the interactions between SVEP1, PEAR1, and the Ang/Tie pathway, with therapeutic implications for a constellation of diseases. This review underscores the potential of molecular GWASs to guide drug discovery and unravel the complexities of human health and disease by demonstrating an integrative approach that grounds mechanistic research in human biology.

NetREm: Network Regression Embeddings reveal cell-type transcription factor coordination for gene regulation

Motivation

Transcription factor (TF) coordination plays a key role in gene regulation via direct and/or indirect protein-protein interactions (PPIs) and co-binding to regulatory elements on DNA. Single-cell technologies facilitate gene expression measurement for individual cells and cell-type identification, yet the connection between TF-TF coordination and target gene (TG) regulation of various cell types remains unclear.

Results

To address this, we introduce our innovative computational approach, Network Regression Embeddings (NetREm), to reveal cell-type TF-TF coordination activities for TG regulation. NetREm leverages network-constrained regularization, using prior knowledge of PPIs among TFs, to analyze single-cell gene expression data, uncovering cell-type coordinating TFs and identifying revolutionary TF-TG candidate regulatory network links. NetREm's performance is validated using simulation studies and benchmarked across several datasets in humans, mice, yeast. Further, we showcase NetREm's ability to prioritize valid novel human TF-TF coordination links in 9 peripheral blood mononuclear and 42 immune cell sub-types. We apply NetREm to examine cell-type networks in central and peripheral nerve systems (e.g. neuronal, glial, Schwann cells) and in Alzheimer's disease versus Controls. Top predictions are validated with experimental data from rat, mouse, and human models. Additional functional genomics data helps link genetic variants to our TF-TG regulatory and TF-TF coordination networks.

Availability and implementation

https://github.com/SaniyaKhullar/NetREm.

Deubiquitinating enzymes in breast cancer: in silico analysis of gene expression and metastatic correlation

Breast cancer, the most prevalent cancer in females, is a heterogeneous disease with various molecular subtypes, which presents challenges in diagnosis and treatment. Ubiquitination is one of the most critical post-translational protein modifications, that plays regulatory roles in numerous cellular processes including cell cycle progression, DNA replication & repair, apoptosis, transcription regulation, protein localization, trafficking and signal transduction. This modification can be reversed by deubiquitinases, or DUBs, a superfamily of cysteine proteases and metalloproteases that cleave ubiquitin-protein bonds. Dysregulation of DUBs has been associated to various diseases including cancer, making them promising targets for cancer therapy. We leveraged publicly available breast cancer datasets and employed various bioinformatics tools to identify differentially expressed DUBs in breast cancer. Our analysis identified six genes (COPS5, EIF3H, MINDY1, MINDY2, PSMD14 and USP26) with significant differential expression and survival implications. We further validated our findings experimentally and found upregulation of COPS5, EIF3H and MINDY 1 in MCF-7 and T47D breast cancer cell lines using qPCR analysis. To identify the role of these genes, EIF3H and COPS5, in disease progression, we constructed a protein-protein interaction (PPI) network with genes associated with metastasis and explored their correlation at the gene expression level in breast cancer patients. Together, this comprehensive study sheds light on DUB gene expression patterns in breast cancer with the potential to identify novel targets for therapeutic interventions.

Co-expression-wide association studies link genetically regulated interactions with complex traits

Transcriptome- and proteome-wide association studies (TWAS/PWAS) have proven successful in prioritizing genes and proteins whose genetically regulated expression modulates disease risk, but they ignore potential co-expression and interaction effects. To address this limitation, we introduce the co-expression-wide association study (COWAS) method, which can identify pairs of genes or proteins whose genetically regulated co-expression is associated with complex traits. COWAS first trains models to predict expression and co-expression conditional on genetic variation, and then tests for association between imputed co-expression and the trait of interest while also accounting for direct effects from each exposure. We applied our method to plasma proteomic concentrations from the UK Biobank, identifying dozens of interacting protein pairs associated with cholesterol levels, Alzheimer's disease, and Parkinson's disease. Notably, our results demonstrate that co-expression between proteins may affect complex traits even if neither protein is detected to influence the trait when considered on its own. We also show how COWAS can help disentangle direct and interaction effects, providing a richer picture of the molecular networks that mediate genetic effects on disease outcomes.

Emergence of power law distributions in protein-protein interaction networks through study bias

Degree distributions in protein-protein interaction (PPI) networks are believed to follow a power law (PL). However, technical and study biases affect the experimental procedures for detecting PPIs. For instance, cancer-associated proteins have received disproportional attention. Moreover, bait proteins in large-scale experiments tend to have many false-positive interaction partners. Studying the degree distributions of thousands of PPI networks of controlled provenance, we address the question if PL distributions in observed PPI networks could be explained by these biases alone. Our findings are supported by mathematical models and extensive simulations, and indicate that study bias and technical bias suffice to produce the observed PL distribution. It is, hence, problematic to derive hypotheses about the topology of the true biological interactome from the PL distributions in observed PPI networks. Our study casts doubt on the use of the PL property of biological networks as a modeling assumption or quality criterion in network biology.

LAP1 Interactome Profiling Provides New Insights into LAP1's Physiological Functions

The nuclear envelope (NE), a protective membrane bordering the nucleus, is composed of highly specialized proteins that are indispensable for normal cellular activity. Lamina-associated polypeptide 1 (LAP1) is a NE protein whose functions are just beginning to be unveiled. The fact that mutations causing LAP1 deficiency are extremely rare and pathogenic is indicative of its paramount importance to preserving human health, anticipating that LAP1 might have a multifaceted role in the cell. Mapping the LAP1 protein interactome is, thus, imperative to achieve an integrated view of its potential biological properties. To this end, we employed in silico- and mass spectrometry-based approaches to identify candidate LAP1-interacting proteins, whose functional attributes were subsequently characterized using bioinformatics tools. Our results reveal the complex and multifunctional network of protein-protein interactions associated to LAP1, evidencing a strong interconnection between LAP1 and cellular processes as diverse as chromatin and cytoskeleton organization, DNA repair, RNA processing and translation, as well as protein biogenesis and turnover, among others. Novel interactions between LAP1 and DNA repair proteins were additionally validated, strengthening the previously proposed involvement of LAP1 in the maintenance of genomic stability. Overall, this study reaffirms the biological relevance of LAP1 and the need to deepen our knowledge about this NE protein, providing new insights about its potential functional partners that will help guiding future research towards a mechanistic understanding of LAP1's functioning.

Heterogeneous network approaches to protein pathway prediction

Understanding protein-protein interactions (PPIs) and the pathways they comprise is essential for comprehending cellular functions and their links to specific phenotypes. Despite the prevalence of molecular data generated by high-throughput sequencing technologies, a significant gap remains in translating this data into functional information regarding the series of interactions that underlie phenotypic differences. In this review, we present an in-depth analysis of heterogeneous network methodologies for modeling protein pathways, highlighting the critical role of integrating multifaceted biological data. It outlines the process of constructing these networks, from data representation to machine learning-driven predictions and evaluations. The work underscores the potential of heterogeneous networks in capturing the complexity of proteomic interactions, thereby offering enhanced accuracy in pathway prediction. This approach not only deepens our understanding of cellular processes but also opens up new possibilities in disease treatment and drug discovery by leveraging the predictive power of comprehensive proteomic data analysis.

The identification of key genes and pathways in polycystic ovary syndrome by bioinformatics analysis of next-generation sequencing data

Background

Polycystic ovary syndrome (PCOS) is a reproductive endocrine disorder. The specific molecular mechanism of PCOS remains unclear. The aim of this study was to apply a bioinformatics approach to reveal related pathways or genes involved in the development of PCOS.

Methods

The next-generation sequencing (NGS) dataset GSE199225 was downloaded from the gene expression omnibus (GEO) database and NGS dataset analyzed is obtained from in vitro culture of PCOS patients’ muscle cells and muscle cells of healthy lean control women. Differentially expressed gene (DEG) analysis was performed using DESeq2. The g:Profiler was utilized to analyze the gene ontology (GO) and REACTOME pathways of the differentially expressed genes. A protein–protein interaction (PPI) network was constructed and module analysis was performed using HiPPIE and cytoscape. The miRNA-hub gene regulatory network and TF-hub gene regulatory network were constructed. The hub genes were validated by using receiver operating characteristic (ROC) curve analysis.

Results

We have identified 957 DEG in total, including 478 upregulated genes and 479 downregulated gene. GO terms and REACTOME pathways illustrated that DEG were significantly enriched in regulation of molecular function, developmental process, interferon signaling and platelet activation, signaling, and aggregation. The top 5 upregulated hub genes including HSPA5, PLK1, RIN3, DBN1, and CCDC85B and top 5 downregulated hub genes including DISC1, AR, MTUS2, LYN, and TCF4 might be associated with PCOS. The hub gens of HSPA5 and KMT2A, together with corresponding predicted miRNAs (e.g., hsa-mir-34b-5p and hsa-mir-378a-5p), and HSPA5 and TCF4 together with corresponding predicted TF (e.g., RCOR3 and TEAD4) were found to be significantly correlated with PCOS.

Conclusions

These study uses of bioinformatics analysis of NGS data to obtain hub genes and key signaling pathways related to PCOS and its associated complications. Also provides novel ideas for finding biomarkers and treatment methods for PCOS and its associated complications.

A modified multiple-criteria decision-making approach based on a protein-protein interaction network to diagnose latent tuberculosis

BACKGROUND

DNA microarrays provide informative data for transcriptional profiling and identifying gene expression signatures to help prevent progression of latent tuberculosis infection (LTBI) to active disease. However, constructing a prognostic model for distinguishing LTBI from active tuberculosis (ATB) is very challenging due to the noisy nature of data and lack of a generally stable analysis approach.

METHODS

In the present study, we proposed an accurate predictive model with the help of data fusion at the decision level. In this regard, results of filter feature selection and wrapper feature selection techniques were combined with multiple-criteria decision-making (MCDM) methods to select 10 genes from six microarray datasets that can be the most discriminative genes for diagnosing tuberculosis cases. As the main contribution of this study, the final ranking function was constructed by combining protein-protein interaction (PPI) network with an MCDM method (called Decision-making Trial and Evaluation Laboratory or DEMATEL) to improve the feature ranking approach.

RESULTS

By applying data fusion at the decision level on the 10 introduced genes in terms of fusion of classifiers of random forests (RF) and k-nearest neighbors (KNN) regarding Yager's theory, the proposed algorithm reached a sensitivity of 0.97, specificity of 0.90, and accuracy of 0.95. Finally, with the help of cumulative clustering, the genes involved in the diagnosis of latent and activated tuberculosis have been introduced.

CONCLUSIONS

The combination of MCDM methods and PPI networks can significantly improve the diagnosis different states of tuberculosis.

CLINICAL TRIAL NUMBER

Not applicable.

Screening and identification of key biomarkers associated with endometriosis using bioinformatics and next-generation sequencing data analysis

Background

Endometriosis is a common cause of endometrial-type mucosa outside the uterine cavity with symptoms such as painful periods, chronic pelvic pain, pain with intercourse and infertility. However, the early diagnosis of endometriosis is still restricted. The purpose of this investigation is to identify and validate the key biomarkers of endometriosis.

Methods

Next-generation sequencing dataset GSE243039 was obtained from the Gene Expression Omnibus database, and differentially expressed genes (DEGs) between endometriosis and normal control samples were identified. After screening of DEGs, gene ontology (GO) and REACTOME pathway enrichment analyses were performed. Furthermore, a protein–protein interaction (PPI) network was constructed and modules were analyzed using the Human Integrated Protein–Protein Interaction rEference database and Cytoscape software, and hub genes were identified. Subsequently, a network between miRNAs and hub genes, and network between TFs and hub genes were constructed using the miRNet and NetworkAnalyst tool, and possible key miRNAs and TFs were predicted. Finally, receiver operating characteristic curve analysis was used to validate the hub genes.

Results

A total of 958 DEGs, including 479 upregulated genes and 479 downregulated genes, were screened between endometriosis and normal control samples. GO and REACTOME pathway enrichment analyses of the 958 DEGs showed that they were mainly involved in multicellular organismal process, developmental process, signaling by GPCR and muscle contraction. Further analysis of the PPI network and modules identified 10 hub genes, including vcam1, snca, prkcb, adrb2, foxq1, mdfi, actbl2, prkd1, dapk1 and actc1. Possible target miRNAs, including hsa-mir-3143 and hsa-mir-2110, and target TFs, including tcf3 (transcription factor 3) and clock (clock circadian regulator), were predicted by constructing a miRNA-hub gene regulatory network and TF-hub gene regulatory network.

Conclusions

This investigation used bioinformatics techniques to explore the potential and novel biomarkers. These biomarkers might provide new ideas and methods for the early diagnosis, treatment and monitoring of endometriosis.

Quantitative proteomic mass spectrometry of protein kinases to determine dynamic heterogeneity of the human kinome

The kinome is a dynamic system of kinases regulating signaling networks in cells and dysfunction of protein kinases contributes to many diseases. Regulation of the protein expression of kinases alters cellular responses to environmental changes and perturbations. We configured a library of 672 proteotypic peptides to quantify >300 kinases in a single LC-MS experiment using ten micrograms protein from human tissues including biopsies. This enables absolute quantitation of kinase protein abundance at attomole-femtomole expression levels, requiring no kinase enrichment and less than ten micrograms of starting protein from flash-frozen and formalin fixed paraffin embedded tissues. Breast cancer biopsies, organoids, and cell lines were analyzed using the SureQuant method, demonstrating the heterogeneity of kinase protein expression across and within breast cancer clinical subtypes. Kinome quantitation was coupled with nanoscale phosphoproteomics, providing a feasible method for novel clinical diagnosis and understanding of patient kinome responses to treatment.

Identification of PPREs and PPRE associated genes in the human genome: insights into related kinases and disease implications

Introduction

"Peroxisome Proliferator-Activated Receptors" (PPARs) belong to the class of transcription factors (TF) identified as Nuclear Receptors (NR). Upon activation by peroxisome proliferators (PPs), PPARs modulate a diverse range of genes, consequently regulating intra-cellular lipid metabolism, glucose uptake, apoptosis, and cell proliferation. Subsequent to the heterodimerization of Retinoid X Receptors (RXR) with PPARs induced by the binding of activators to PPARs, facilitates the binding of the resulting complex to Peroxisome Proliferator-Activated Receptors Response Elements (PPRE), with a consensus sequence 5'AGGTCANAGGTCA-3', and regulate the transcription of the targeted genes.

Methods

A comprehensive screening of PPRE within the whole human genome was performed using the Genome Workbench and UCSC Genome Browser to find the associated genes. Subsequently, the kinase subset was isolated from the extracted list of PPRE-related genes. Functional enrichment of the kinases was performed using FunRich, ToppGene, and ShinyGO. Network analysis and enrichment studies were then further performed using NDEx to elucidate these identified kinases' connections and significance. Additionally, the disease association of the PPRE kinases was analyzed using DisGeNET data in R studio and the COSMIC dataset.

Results

A comprehensive analysis of 1002 PPRE sequences within the human genome (T2T), yielded the identification of 660 associated genes, including 29 kinases. The engagement of these kinases in various biological pathways, such as apoptosis, platelet activation, and cytokine pathways, revealed from the functional enrichment analysis, illuminates the multifaceted role of PPAR in the regulation of cellular homeostasis and biological processes. Network analysis reveals the kinases interact with approximately 5.56% of the Human Integrated Protein-Protein Interaction rEference (HIPPIE) network. Disease association analysis using DisGeNET and COSMIC datasets revealed the significant roles of these kinases in cellular processes and disease modulation.

Discussion

This study elucidates the regulatory role of PPAR-associated genes and their association with numerous biological pathways. The involvement of the kinases with disease-related pathways highlights new potential for the development of therapeutic strategies designed for disease management and intervention.

Computational tools to predict context-specific protein complexes

Interactions between thousands of proteins define cells' protein-protein interaction (PPI) network. Some of these interactions lead to the formation of protein complexes. It is challenging to identify a protein complex in a haystack of protein-protein interactions, and it is even more difficult to predict all protein complexes of the complexome. Simulations and machine learning approaches try to crack these problems by looking at the PPI network or predicted protein structures. Clustering of PPI networks led to the first protein complex predictions, while most recently, atomistic models of protein complexes and deep-learning-based structure prediction methods have also emerged. The simulation of PPI level interactions even enables the quantitative prediction of protein complexes. These methods, the required data sources, and their potential future developments are discussed in this review.

Predicting the involvement of polyQ- and polyA in protein-protein interactions by their amino acid context

Homorepeats, specifically polyglutamine (polyQ) and polyalanine (polyA), are often implicated in protein-protein interactions (PPIs). So far, a method to predict the participation of homorepeats in protein interactions is lacking. We propose a machine learning approach to identify PPI-involved polyQ and polyA regions within the human proteome based on known interacting regions. Using the dataset of human homorepeats, we identified 157 polyQ and 745 polyA regions potentially involved in PPIs. Machine learning models, trained on amino acid context and homorepeat length, demonstrated high precision (0.90-0.98) but variable recall (0.42-0.85). Random forest outperformed other models (AUC polyQ = 0.686, AUC polyA = 0.732) using the positions surrounding the homorepeat -10 to +10. Integrating paralog information marginally improved predictions but was excluded for model simplicity. Further optimization revealed that for polyQ, using amino acid surrounding positions from -6 to +6 increased AUC to 0.715. For polyA, no improvement was found. Incorporating coiled coil overlap information enhanced polyA predictions (AUC = 0.745) but not polyQ. Finally, we applied these models to predict PPI involvement across all polyQ and polyA regions, identifying potential interactions. Case studies illustrated the method's predictive capacity, highlighting known interacting regions with high scores and elucidating potential false negatives.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here, we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced histone 3 lysine 27 acetylation (H3K27ac). Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data support clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

In Silico Exploration of AHR-HIF Pathway Interplay: Implications for Therapeutic Targeting in ccRCC

The oxygen-sensing pathway is a crucial regulatory circuit that defines cellular conditions and is extensively exploited in cancer development. Pathogenic mutations in the von Hippel-Lindau (VHL) tumour suppressor impair its role as a master regulator of hypoxia-inducible factors (HIFs), leading to constitutive HIF activation and uncontrolled angiogenesis, increasing the risk of developing clear cell renal cell carcinoma (ccRCC). HIF hyperactivation can sequester HIF-1β, preventing the aryl hydrocarbon receptor (AHR) from correctly activating gene expression in response to endogenous and exogenous ligands such as TCDD (dioxins). In this study, we used protein-protein interaction networks and gene expression profiling to characterize the impact of VHL loss on AHR activity. Our findings reveal specific expression patterns of AHR interactors following exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and in ccRCC. We identified several AHR interactors significantly associated with poor survival rates in ccRCC patients. Notably, the upregulation of the androgen receptor (AR) and retinoblastoma-associated protein (RB1) by TCDD, coupled with their respective downregulation in ccRCC and association with poor survival rates, suggests novel therapeutic targets. The strategic activation of the AHR via selective AHR modulators (SAhRMs) could stimulate its anticancer activity, specifically targeting RB1 and AR to reduce cell cycle progression and metastasis formation in ccRCC. Our study provides comprehensive insights into the complex interplay between the AHR and HIF pathways in ccRCC pathogenesis, offering novel strategies for targeted therapeutic interventions.

Proteome-scale characterisation of motif-based interactome rewiring by disease mutations

Whole genome and exome sequencing are reporting on hundreds of thousands of missense mutations. Taking a pan-disease approach, we explored how mutations in intrinsically disordered regions (IDRs) break or generate protein interactions mediated by short linear motifs. We created a peptide-phage display library tiling ~57,000 peptides from the IDRs of the human proteome overlapping 12,301 single nucleotide variants associated with diverse phenotypes including cancer, metabolic diseases and neurological diseases. By screening 80 human proteins, we identified 366 mutation-modulated interactions, with half of the mutations diminishing binding, and half enhancing binding or creating novel interaction interfaces. The effects of the mutations were confirmed by affinity measurements. In cellular assays, the effects of motif-disruptive mutations were validated, including loss of a nuclear localisation signal in the cell division control protein CDC45 by a mutation associated with Meier-Gorlin syndrome. The study provides insights into how disease-associated mutations may perturb and rewire the motif-based interactome.

Computational analysis of pathogen-host interactome for fast and low-risk in-silico drug repurposing in emerging viral threats like Mpox

Monkeypox (Mpox), a zoonotic illness triggered by the monkeypox virus (MPXV), poses a significant threat since it may be transmitted and has no cure. This work introduces a computational method to predict Protein-Protein Interactions (PPIs) during MPXV infection. The objective is to discover prospective drug targets and repurpose current potential Food and Drug Administration (FDA) drugs for therapeutic purposes. In this work, ensemble features, comprising 2-5 node graphlet attributes and protein composition-based features are utilized for Deep Learning (DL) models to predict PPIs. The technique that is used here demonstrated an excellent prediction performance for PPI on both the Human Integrated Protein-Protein Interaction Reference (HIPPIE) and MPXV-Human PPI datasets. In addition, the human protein targets for MPXV have been identified accurately along with the detection of possible therapeutic targets. Furthermore, the validation process included conducting docking research studies on potential FDA drugs like Nicotinamide Adenine Dinucleotide and Hydrogen (NADH), Fostamatinib, Glutamic acid, Cannabidiol, Copper, and Zinc in DrugBank identified via research on drug repurposing and the Drug Consensus Score (DCS) for MPXV. This has been achieved by employing the primary crystal structures of MPXV, which are now accessible. The docking study is also supported by Molecular Dynamics (MD) simulation. The results of our study emphasize the effectiveness of using ensemble feature-based PPI prediction to understand the molecular processes involved in viral infection and to aid in the development of repurposed drugs for emerging infectious diseases such as, but not limited to, Mpox. The source code and link to data used in this work is available at: https://github.com/CMATERJU-BIOINFO/In-Silico-Drug-Repurposing-Methodology-To-Suggest-Therapies-For-Emerging-Threats-like-Mpox .

Transcriptomics and weighted protein network analyses of the LRRK2 protein interactome reveal distinct molecular signatures for sporadic and LRRK2 Parkinson's Disease

Mutations in the LRRK2 gene are the most common genetic cause of familial Parkinson's Disease (LRRK2-PD) and an important risk factor for sporadic PD (sPD). Multiple clinical trials are ongoing to evaluate the benefits associated with the therapeutical reduction of LRRK2 kinase activity. In this study, we described the changes of transcriptomic profiles (whole blood mRNA levels) of LRRK2 protein interactors in sPD and LRRK2-PD cases as compared to healthy controls with the aim of comparing the two PD conditions. We went on to model the protein-protein interaction (PPI) network centred on LRRK2, which was weighted to reflect the transcriptomic changes on expression and co-expression levels of LRRK2 protein interactors. Our results showed that LRRK2 interactors present both similar and distinct alterations in expression levels and co-expression behaviours in the sPD and LRRK2-PD cases; suggesting that, albeit being classified as the same disease based on clinical features, LRRK2-PD and sPD display significant differences from a molecular perspective. Interestingly, the similar changes across the two PD conditions result in decreased connectivity within a topological cluster of the LRRK2 PPI network associated with protein metabolism/biosynthesis and ribosomal metabolism suggesting protein homoeostasis and ribosomal dynamics might be affected in both sporadic and familial PD in comparison with controls.

Unveiling hidden connections in omics data via pyPARAGON: an integrative hybrid approach for disease network construction

Network inference or reconstruction algorithms play an integral role in successfully analyzing and identifying causal relationships between omics hits for detecting dysregulated and altered signaling components in various contexts, encompassing disease states and drug perturbations. However, accurate representation of signaling networks and identification of context-specific interactions within sparse omics datasets in complex interactomes pose significant challenges in integrative approaches. To address these challenges, we present pyPARAGON (PAgeRAnk-flux on Graphlet-guided network for multi-Omic data integratioN), a novel tool that combines network propagation with graphlets. pyPARAGON enhances accuracy and minimizes the inclusion of nonspecific interactions in signaling networks by utilizing network rather than relying on pairwise connections among proteins. Through comprehensive evaluations on benchmark signaling pathways, we demonstrate that pyPARAGON outperforms state-of-the-art approaches in node propagation and edge inference. Furthermore, pyPARAGON exhibits promising performance in discovering cancer driver networks. Notably, we demonstrate its utility in network-based stratification of patient tumors by integrating phosphoproteomic data from 105 breast cancer tumors with the interactome and demonstrating tumor-specific signaling pathways. Overall, pyPARAGON is a novel tool for analyzing and integrating multi-omic data in the context of signaling networks. pyPARAGON is available at https://github.com/netlab-ku/pyPARAGON.

TUnA: an uncertainty-aware transformer model for sequence-based protein-protein interaction prediction

Protein-protein interactions (PPIs) are important for many biological processes, but predicting them from sequence data remains challenging. Existing deep learning models often cannot generalize to proteins not present in the training set and do not provide uncertainty estimates for their predictions. To address these limitations, we present TUnA, a Transformer-based uncertainty-aware model for PPI prediction. TUnA uses ESM-2 embeddings with Transformer encoders and incorporates a Spectral-normalized Neural Gaussian Process. TUnA achieves state-of-the-art performance and, importantly, evaluates uncertainty for unseen sequences. We demonstrate that TUnA's uncertainty estimates can effectively identify the most reliable predictions, significantly reducing false positives. This capability is crucial in bridging the gap between computational predictions and experimental validation.

GPS-SUMO 2.0: an updated online service for the prediction of SUMOylation sites and SUMO-interacting motifs

Small ubiquitin-like modifiers (SUMOs) are tiny but important protein regulators involved in orchestrating a broad spectrum of biological processes, either by covalently modifying protein substrates or by noncovalently interacting with other proteins. Here, we report an updated server, GPS-SUMO 2.0, for the prediction of SUMOylation sites and SUMO-interacting motifs (SIMs). For predictor training, we adopted three machine learning algorithms, penalized logistic regression (PLR), a deep neural network (DNN), and a transformer, and used 52 404 nonredundant SUMOylation sites in 8262 proteins and 163 SIMs in 102 proteins. To further increase the accuracy of predicting SUMOylation sites, a pretraining model was first constructed using 145 545 protein lysine modification sites, followed by transfer learning to fine-tune the model. GPS-SUMO 2.0 exhibited greater accuracy in predicting SUMOylation sites than did other existing tools. For users, one or multiple protein sequences or identifiers can be input, and the prediction results are shown in a tabular list. In addition to the basic statistics, we integrated knowledge from 35 public resources to annotate SUMOylation sites or SIMs. The GPS-SUMO 2.0 server is freely available at https://sumo.biocuckoo.cn/. We believe that GPS-SUMO 2.0 can serve as a useful tool for further analysis of SUMOylation and SUMO interactions.

Review: Cancer and neurodevelopmental disorders: multi-scale reasoning and computational guide

The connection and causality between cancer and neurodevelopmental disorders have been puzzling. How can the same cellular pathways, proteins, and mutations lead to pathologies with vastly different clinical presentations? And why do individuals with neurodevelopmental disorders, such as autism and schizophrenia, face higher chances of cancer emerging throughout their lifetime? Our broad review emphasizes the multi-scale aspect of this type of reasoning. As these examples demonstrate, rather than focusing on a specific organ system or disease, we aim at the new understanding that can be gained. Within this framework, our review calls attention to computational strategies which can be powerful in discovering connections, causalities, predicting clinical outcomes, and are vital for drug discovery. Thus, rather than centering on the clinical features, we draw on the rapidly increasing data on the molecular level, including mutations, isoforms, three-dimensional structures, and expression levels of the respective disease-associated genes. Their integrated analysis, together with chromatin states, can delineate how, despite being connected, neurodevelopmental disorders and cancer differ, and how the same mutations can lead to different clinical symptoms. Here, we seek to uncover the emerging connection between cancer, including pediatric tumors, and neurodevelopmental disorders, and the tantalizing questions that this connection raises.

TSG-6 Is Involved in Fibrous Structural Remodeling after the Injection of Adipose-derived Stem Cells

Background

Although aesthetic treatments can rejuvenate the skin, they often cause specific forms of tissue damage. Unlike wounding, which typically results in fibrotic scar tissue, damage from aesthetic treatments induces a distinct histological rejuvenation. The mechanisms that drive this rejuvenation are not yet fully understood. Here, we were interested in cellular responses following aesthetic treatments injecting adipose-derived stem cells (ASCs) subcutaneously. Through investigation with an ex vivo experimental model, a key gene was identified that orchestrates fibrous structural changes and tissue remodeling.

Methods

Using fresh human subcutaneous adipose tissue co-cultured with ASCs, the changes in the fibrous architecture of the tissue were sequentially mapped. The key regulatory genes involved in remodeling were identified using gene expression and computational analyses.

Results

We identified the regulatory elements that are crucial for tissue remodeling. Among those, we found that tumor necrosis factor-stimulated gene-6 (TSG-6) is a paracrine mediator essential for the collagen activity. It not only alleviates tissue inflammation but also promotes collagen replacement ex vivo. This is primarily achieved by inhibiting the formation of neutrophil extracellular traps, which are known to promote fibrosis.

Conclusions

TSG-6 is a key factor modulating tissue inflammation. As our results demonstrate, after ASCs treatment, this factor directs skin healing away from fibrosis by reducing neutrophil extracellular trap formation in subcutaneous adipose tissue and promotes fibrous rejuvenation.

Role of Protein Phosphatases in Tumor Angiogenesis: Assessing PP1, PP2A, PP2B and PTPs Activity

Tumor angiogenesis, the formation of new blood vessels to support tumor growth and metastasis, is a complex process regulated by a multitude of signaling pathways. Dysregulation of signaling pathways involving protein kinases has been extensively studied, but the role of protein phosphatases in angiogenesis within the tumor microenvironment remains less explored. However, among angiogenic pathways, protein phosphatases play critical roles in modulating signaling cascades. This review provides a comprehensive overview of the involvement of protein phosphatases in tumor angiogenesis, highlighting their diverse functions and mechanisms of action. Protein phosphatases are key regulators of cellular signaling pathways by catalyzing the dephosphorylation of proteins, thereby modulating their activity and function. This review aims to assess the activity of the protein tyrosine phosphatases and serine/threonine phosphatases. These phosphatases exert their effects on angiogenic signaling pathways through various mechanisms, including direct dephosphorylation of angiogenic receptors and downstream signaling molecules. Moreover, protein phosphatases also crosstalk with other signaling pathways involved in angiogenesis, further emphasizing their significance in regulating tumor vascularization, including endothelial cell survival, sprouting, and vessel maturation. In conclusion, this review underscores the pivotal role of protein phosphatases in tumor angiogenesis and accentuate their potential as therapeutic targets for anti-angiogenic therapy in cancer.

Genopyc: a Python library for investigating the functional effects of genomic variants associated to complex diseases

MOTIVATION

Integrative Biomedicl Informatics, Research Program on Biomedical Informatics (IBI - GRIB), Hospital Del Mar Medical Research Institute (IMIM), Department of Experimental and Health Sciences, Universitat Pompeu Fabra (UPF) C/ del Dr. Aiguader 88 Barcelona 08003 Spain.Understanding the genetic basis of complex diseases is one of the main challenges in modern genomics. However, current tools often lack the versatility to efficiently analyze the intricate relationships between genetic variations and disease outcomes. To address this, we introduce Genopyc, a novel Python library designed for comprehensive investigation of how the variants associated to complex diseases affects downstream pathways. Genopyc offers an extensive suite of functions for heterogeneous data mining and visualization, enabling researchers to delve into and integrate biological information from large-scale genomic datasets.

RESULTS

In this work, we present the Genopyc library through application to real-world genome wide association studies variants. Using Genopyc to investigate the functional consequences of variants associated to intervertebral disc degeneration enabled a deeper understanding of the potential dysregulated pathways involved in the disease, which can be explored and visualized by exploiting the functionalities featured in the package. Genopyc emerges as a powerful asset for researchers, facilitating the investigation of complex diseases paving the way for more targeted therapeutic interventions.

AVAILABILITY AND IMPLEMENTATION

Genopyc is available on pip https://pypi.org/project/genopyc/.The source code of Genopyc is available at https://github.com/freh-g/genopyc. A tutorial notebook is available at https://github.com/freh-g/genopyc/blob/main/tutorials/Genopyc_tutorial_notebook.ipynb. Finally, a detailed documentation is available at: https://genopyc.readthedocs.io/en/latest/.

Polyglutamine disease proteins: Commonalities and differences in interaction profiles and pathological effects

Currently, nine polyglutamine (polyQ) expansion diseases are known. They include spinocerebellar ataxias (SCA1, 2, 3, 6, 7, 17), spinal and bulbar muscular atrophy (SBMA), dentatorubral-pallidoluysian atrophy (DRPLA), and Huntington's disease (HD). At the root of these neurodegenerative diseases are trinucleotide repeat mutations in coding regions of different genes, which lead to the production of proteins with elongated polyQ tracts. While the causative proteins differ in structure and molecular mass, the expanded polyQ domains drive pathogenesis in all these diseases. PolyQ tracts mediate the association of proteins leading to the formation of protein complexes involved in gene expression regulation, RNA processing, membrane trafficking, and signal transduction. In this review, we discuss commonalities and differences among the nine polyQ proteins focusing on their structure and function as well as the pathological features of the respective diseases. We present insights from AlphaFold-predicted structural models and discuss the biological roles of polyQ-containing proteins. Lastly, we explore reported protein-protein interaction networks to highlight shared protein interactions and their potential relevance in disease development.

Deciphering the similarities and disparities of molecular mechanisms behind respiratory epithelium response to HCoV-229E and SARS-CoV-2 and drug repurposing, a systems biology approach

PURPOSE

Identifying the molecular mechanisms behind SARS-CoV-2 disparities and similarities will help find new treatments. The present study determines networks' shared and non-shared (specific) crucial elements in response to HCoV-229E and SARS-CoV-2 viruses to recommend candidate medications.

METHODS

We retrieved the omics data on respiratory cells infected with HCoV-229E and SARS-CoV-2, constructed PPIN and GRN, and detected clusters and motifs. Using a drug-gene interaction network, we determined the similarities and disparities of mechanisms behind their host response and drug-repurposed.

RESULTS

CXCL1, KLHL21, SMAD3, HIF1A, and STAT1 were the shared DEGs between both viruses' protein-protein interaction network (PPIN) and gene regulatory network (GRN). The NPM1 was a specific critical node for HCoV-229E and was a Hub-Bottleneck shared between PPI and GRN in HCoV-229E. The HLA-F, ADCY5, TRIM14, RPF1, and FGA were the seed proteins in subnetworks of the SARS-CoV-2 PPI network, and HSPA1A and RPL26 proteins were the seed in subnetworks of the PPI network of HCOV-229E. TRIM14, STAT2, and HLA-F played the same role for SARS-CoV-2. Top enriched KEGG pathways included cell cycle and proteasome in HCoV-229E and RIG-I-like receptor, Chemokine, Cytokine-cytokine, NOD-like receptor, and TNF signaling pathways in SARS-CoV-2. We suggest some candidate medications for COVID-19 patient lungs, including Noscapine, Isoetharine mesylate, Cycloserine, Ethamsylate, Cetylpyridinium, Tretinoin, Ixazomib, Vorinostat, Venetoclax, Vorinostat, Ixazomib, Venetoclax, and epoetin alfa for further in-vitro and in-vivo investigations.

CONCLUSION

We suggested CXCL1, KLHL21, SMAD3, HIF1A, and STAT1, ADCY5, TRIM14, RPF1, and FGA, STAT2, and HLA-F as critical genes and Cetylpyridinium, Cycloserine, Noscapine, Ethamsylate, Epoetin alfa, Isoetharine mesylate, Ribavirin, and Tretinoin drugs to study further their importance in treating COVID-19 lung complications.

Decoding common genetic alterations between Barrett's esophagus and esophageal adenocarcinoma: A bioinformatics analysis

Background

Esophageal adenocarcinoma (EAC) is a common cancer with a poor prognosis in advanced stages. Therefore, early EAC diagnosis and treatment have gained attention in recent decades. It has been found that various pathological changes, particularly Barrett's Esophagus (BE), can occur in the esophageal tissue before the development of EAC. In this study, we aimed to identify the molecular contributor in BE to EAC progression by detecting the essential regulatory genes that are differentially expressed in both BE and EAC.

Materials and methods

We conducted a comprehensive bioinformatics analysis to detect BE and EAC-associated genes. The common differentially expressed genes (DEGs) and common single nucleotide polymorphisms (SNPs) were detected using the GEO and DisGeNET databases, respectively. Then, hub genes and the top modules within the protein-protein interaction network were identified. Moreover, the co-expression network of the top module by the HIPPIE database was constructed. Additionally, the gene regulatory network was constructed based on miRNAs and circRNAs. Lastly, we inspected the DGIdb database for possible interacted drugs.

Results

Our microarray dataset analysis identified 92 common DEGs between BE and EAC with significant enrichment in skin and epidermis development genes. The study also identified 22 common SNPs between BE and EAC. The top module of PPI network analysis included SCEL, KRT6A, SPRR1A, SPRR1B, SPRR3, PPL, SPRR2B, EVPL, and CSTA. We constructed a ceRNA network involving three specific mRNAs, 23 miRNAs, and 101 selected circRNAs. According to the results from the DGIdb database, TD101 was found to interact with the KRT6A gene.

Conclusion

The present study provides novel potential candidate genes that may be involved in the molecular association between Esophageal adenocarcinoma and Barrett's Esophagus, resulting in developing the diagnostic tools and therapeutic targets to prevent progression of BE to EAC.

Mapping protein-protein interactions by mass spectrometry

Protein-protein interactions (PPIs) are essential for numerous biological activities, including signal transduction, transcription control, and metabolism. They play a pivotal role in the organization and function of the proteome, and their perturbation is associated with various diseases, such as cancer, neurodegeneration, and infectious diseases. Recent advances in mass spectrometry (MS)-based protein interactomics have significantly expanded our understanding of the PPIs in cells, with techniques that continue to improve in terms of sensitivity, and specificity providing new opportunities for the study of PPIs in diverse biological systems. These techniques differ depending on the type of interaction being studied, with each approach having its set of advantages, disadvantages, and applicability. This review highlights recent advances in enrichment methodologies for interactomes before MS analysis and compares their unique features and specifications. It emphasizes prospects for further improvement and their potential applications in advancing our knowledge of PPIs in various biological contexts.

DNA-Binding Proteins and Passenger Proteins in Plasma DNA-Protein Complexes: Imprint of Parental Cells or Key Mediators of Carcinogenesis Processes?

Knowledge of the composition of proteins that interact with plasma DNA will provide a better understanding of the homeostasis of circulating nucleic acids and the various modes of interaction with target cells, which may be useful in the development of gene targeted therapy approaches. The goal of the present study is to shed light on the composition and architecture of histone-containing nucleoprotein complexes (NPCs) from the blood plasma of healthy females (HFs) and breast cancer patients (BCPs) and to explore the relationship of proteins with crucial steps of tumor progression: epithelial-mesenchymal transition (EMT), cell proliferation, invasion, cell migration, stimulation of angiogenesis, and immune response. MALDI-TOF mass spectrometric analysis of NPCs isolated from blood samples using affine chromatography was performed. Bioinformatics analysis showed that the shares of DNA-binding proteins in the compositions of NPCs in normal and cancer patients are comparable and amount to 40% and 33%, respectively; in total, we identified 38 types of DNA-binding motifs. Functional enrichment analysis using FunRich 3.13 showed that, in BCP blood, the share of DNA-binding proteins involved in nucleic acid metabolism increased, while the proportion of proteins involved in intercellular communication and signal transduction decreased. The representation of NPC passenger proteins in breast cancer also changes: the proportion of proteins involved in transport increases and the share of proteins involved in energy biological pathways decreases. Moreover, in the HF blood, proteins involved in the processes of apoptosis were more represented in the composition of NPCs and in the BCP blood-in the processes of active secretion. For the first time, bioinformatics approaches were used to visualize the architecture of circulating NPCs in the blood and to show that breast cancer has an increased representation of passenger proteins involved in EMT, cell proliferation, invasion, cell migration, and immune response. Using breast cancer protein data from the Human Protein Atlas (HPA) and DEPC, we found that 86% of NPC proteins in the blood of BCPs were not previously annotated in these databases. The obtained data may indirectly indicate directed protein sorting in NPCs, which, along with extracellular vesicles, can not only be diagnostically significant molecules for liquid biopsy, but can also carry out the directed transfer of genetic material from donor cells to recipient cells.

Wnt, glucocorticoid and cellular prion protein cooperate to drive a mesenchymal phenotype with poor prognosis in colon cancer

BACKGROUND

The mesenchymal subtype of colorectal cancer (CRC), associated with poor prognosis, is characterized by abundant expression of the cellular prion protein PrPC, which represents a candidate therapeutic target. How PrPC is induced in CRC remains elusive. This study aims to elucidate the signaling pathways governing PrPC expression and to shed light on the gene regulatory networks linked to PrPC.

METHODS

We performed in silico analyses on diverse datasets of in vitro, ex vivo and in vivo models of mouse CRC and patient cohorts. We mined ChIPseq studies and performed promoter analysis. CRC cell lines were manipulated through genetic and pharmacological approaches. We created mice combining conditional inactivation of Apc in intestinal epithelial cells and overexpression of the human prion protein gene PRNP. Bio-informatic analyses were carried out in two randomized control trials totalizing over 3000 CRC patients.

RESULTS

In silico analyses combined with cell-based assays identified the Wnt-β-catenin and glucocorticoid pathways as upstream regulators of PRNP expression, with subtle differences between mouse and human. We uncover multiple feedback loops between PrPC and these two pathways, which translate into an aggravation of CRC pathogenesis in mouse. In stage III CRC patients, the signature defined by PRNP-CTNNB1-NR3C1, encoding PrPC, β-catenin and the glucocorticoid receptor respectively, is overrepresented in the poor-prognosis, mesenchymal subtype and associates with reduced time to recurrence.

CONCLUSIONS

An unleashed PrPC-dependent vicious circle is pathognomonic of poor prognosis, mesenchymal CRC. Patients from this aggressive subtype of CRC may benefit from therapies targeting the PRNP-CTNNB1-NR3C1 axis.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced H3K27ac. Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data supports clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.