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

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

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.

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.

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.

Precise prediction of phase-separation key residues by machine learning

Understanding intracellular phase separation is crucial for deciphering transcriptional control, cell fate transitions, and disease mechanisms. However, the key residues, which impact phase separation the most for protein phase separation function have remained elusive. We develop PSPHunter, which can precisely predict these key residues based on machine learning scheme. In vivo and in vitro validations demonstrate that truncating just 6 key residues in GATA3 disrupts phase separation, enhancing tumor cell migration and inhibiting growth. Glycine and its motifs are enriched in spacer and key residues, as revealed by our comprehensive analysis. PSPHunter identifies nearly 80% of disease-associated phase-separating proteins, with frequent mutated pathological residues like glycine and proline often residing in these key residues. PSPHunter thus emerges as a crucial tool to uncover key residues, facilitating insights into phase separation mechanisms governing transcriptional control, cell fate transitions, and disease development.

TCMNPAS: a comprehensive analysis platform integrating network formulaology and network pharmacology for exploring traditional Chinese medicine

The application of network formulaology and network pharmacology has significantly advanced the scientific understanding of traditional Chinese medicine (TCM) treatment mechanisms in disease. The field of herbal biology is experiencing a surge in data generation. However, researchers are encountering challenges due to the fragmented nature of the data and the reliance on programming tools for data analysis. We have developed TCMNPAS, a comprehensive analysis platform that integrates network formularology and network pharmacology. This platform is designed to investigate in-depth the compatibility characteristics of TCM formulas and their potential molecular mechanisms. TCMNPAS incorporates multiple resources and offers a range of functions designed for automated analysis implementation, including prescription mining, molecular docking, network pharmacology analysis, and visualization. These functions enable researchers to analyze and obtain core herbs and core formulas from herbal prescription data through prescription mining. Additionally, TCMNPAS facilitates virtual screening of active compounds in TCM and its formulas through batch molecular docking, allowing for the rapid construction and analysis of networks associated with "herb-compound-target-pathway" and disease targets. Built upon the integrated analysis concept of network formulaology and network pharmacology, TCMNPAS enables quick point-and-click completion of network-based association analysis, spanning from core formula mining from clinical data to the exploration of therapeutic targets for disease treatment. TCMNPAS serves as a powerful platform for uncovering the combinatorial rules and mechanism of TCM formulas holistically. We distribute TCMNPAS within an open-source R package at GitHub ( https://github.com/yangpluszhu/tcmnpas ), and the project is freely available at http://54.223.75.62:3838/ .

Structured Tandem Repeats in Protein Interactions

Tandem repeats (TRs) in protein sequences are consecutive, highly similar sequence motifs. Some types of TRs fold into structural units that pack together in ensembles, forming either an (open) elongated domain or a (closed) propeller, where the last unit of the ensemble packs against the first one. Here, we examine TR proteins (TRPs) to see how their sequence, structure, and evolutionary properties favor them for a function as mediators of protein interactions. Our observations suggest that TRPs bind other proteins using large, structured surfaces like globular domains; in particular, open-structured TR ensembles are favored by flexible termini and the possibility to tightly coil against their targets. While, intuitively, open ensembles of TRs seem prone to evolve due to their potential to accommodate insertions and deletions of units, these evolutionary events are unexpectedly rare, suggesting that they are advantageous for the emergence of the ancestral sequence but are early fixed. We hypothesize that their flexibility makes it easier for further proteins to adapt to interact with them, which would explain their large number of protein interactions. We provide insight into the properties of open TR ensembles, which make them scaffolds for alternative protein complexes to organize genes, RNA and proteins.

Proteome-wide assessment of human interactome as a source of capturing domain-motif and domain-domain interactions

Protein-protein interactions (PPIs) play a crucial role in various biological processes by establishing domain-motif (DMI) and domain-domain interactions (DDIs). While the existence of real DMIs/DDIs is generally assumed, it is rarely tested; therefore, this study extensively compared high-throughput methods and public PPI repositories as sources for DMI and DDI prediction based on the assumption that the human interactome provides sufficient data for the reliable identification of DMIs and DDIs. Different datasets from leading high-throughput methods (Yeast two-hybrid [Y2H], Affinity Purification coupled Mass Spectrometry [AP-MS], and Co-fractionation-coupled Mass Spectrometry) were assessed for their ability to capture DMIs and DDIs using known DMI/DDI information. High-throughput methods were not notably worse than PPI databases and, in some cases, appeared better. In conclusion, all PPI datasets demonstrated significant enrichment in DMIs and DDIs (p-value <0.001), establishing Y2H and AP-MS as reliable methods for predicting these interactions. This study provides valuable insights for biologists in selecting appropriate methods for predicting DMIs, ultimately aiding in SLiM discovery.

Illuminating function of the understudied druggable kinome

The human kinome, with more than 500 proteins, is crucial for cell signaling and disease. Yet, about one-third of kinases lack in-depth study. The Data and Resource Generating Center for Understudied Kinases has developed multiple resources to address this challenge including creation of a heavy amino acid peptide library for parallel reaction monitoring and quantitation of protein kinase expression, use of understudied kinases tagged with a miniTurbo-biotin ligase to determine interaction networks by proximity-dependent protein biotinylation, NanoBRET probe development for screening chemical tool target specificity in live cells, characterization of small molecule chemical tools inhibiting understudied kinases, and computational tools for defining kinome architecture. These resources are available through the Dark Kinase Knowledgebase, supporting further research into these understudied protein kinases.

Parkinson's disease and schizophrenia interactomes contain temporally distinct gene clusters underlying comorbid mechanisms and unique disease processes

Genome-wide association studies suggest significant overlaps in Parkinson's disease (PD) and schizophrenia (SZ) risks, but the underlying mechanisms remain elusive. The protein-protein interaction network ('interactome') plays a crucial role in PD and SZ and can incorporate their spatiotemporal specificities. Therefore, to study the linked biology of PD and SZ, we compiled PD- and SZ-associated genes from the DisGeNET database, and constructed their interactomes using BioGRID and HPRD. We examined the interactomes using clustering and enrichment analyses, in conjunction with the transcriptomic data of 26 brain regions spanning foetal stages to adulthood available in the BrainSpan Atlas. PD and SZ interactomes formed four gene clusters with distinct temporal identities (Disease Gene Networks or 'DGNs'1-4). DGN1 had unique SZ interactome genes highly expressed across developmental stages, corresponding to a neurodevelopmental SZ subtype. DGN2, containing unique SZ interactome genes expressed from early infancy to adulthood, correlated with an inflammation-driven SZ subtype and adult SZ risk. DGN3 contained unique PD interactome genes expressed in late infancy, early and late childhood, and adulthood, and involved in mitochondrial pathways. DGN4, containing prenatally-expressed genes common to both the interactomes, involved in stem cell pluripotency and overlapping with the interactome of 22q11 deletion syndrome (comorbid psychosis and Parkinsonism), potentially regulates neurodevelopmental mechanisms in PD-SZ comorbidity. Our findings suggest that disrupted neurodevelopment (regulated by DGN4) could expose risk windows in PD and SZ, later elevating disease risk through inflammation (DGN2). Alternatively, variant clustering in DGNs may produce disease subtypes, e.g., PD-SZ comorbidity with DGN4, and early/late-onset SZ with DGN1/DGN2.

Active natural compounds perturb the melanoma risk-gene network

Cutaneous melanoma is an aggressive type of skin cancer with a complex genetic landscape caused by the malignant transformation of melanocytes. This study aimed at providing an in silico network model based on the systematic profiling of the melanoma-associated genes considering germline mutations, somatic mutations, and genome-wide association study signals accounting for a total of 232 unique melanoma risk genes. A protein-protein interaction network was constructed using the melanoma risk genes as seeds and evaluated to describe the functional landscape in which the melanoma genes operate within the cellular milieu. Not only were the majority of the melanoma risk genes able to interact with each other at the protein level within the core of the network, but this showed significant enrichment for genes whose expression is altered in human melanoma specimens. Functional annotation showed the melanoma risk network to be significantly associated with processes related to DNA metabolism and telomeres, DNA damage and repair, cellular ageing, and response to radiation. We further explored whether the melanoma risk network could be used as an in silico tool to predict the efficacy of anti-melanoma phytochemicals, that are considered active molecules with potentially less systemic toxicity than classical cytotoxic drugs. A significant portion of the melanoma risk network showed differential expression when SK-MEL-28 human melanoma cells were exposed to the phytochemicals harmine and berberine chloride. This reinforced our hypothesis that the network modeling approach not only provides an alternative way to identify molecular pathways relevant to disease but it may also represent an alternative screening approach to prioritize potentially active compounds.

Shared Proteins and Pathways of Cardiovascular and Cognitive Diseases: Relation to Vascular Cognitive Impairment

One of the primary goals of systems medicine is the detection of putative proteins and pathways involved in disease progression and pathological phenotypes. Vascular cognitive impairment (VCI) is a heterogeneous condition manifesting as cognitive impairment resulting from vascular factors. The precise mechanisms underlying this relationship remain unclear, which poses challenges for experimental research. Here, we applied computational approaches like systems biology to unveil and select relevant proteins and pathways related to VCI by studying the crosstalk between cardiovascular and cognitive diseases. In addition, we specifically included signals related to oxidative stress, a common etiologic factor tightly linked to aging, a major determinant of VCI. Our results show that pathways associated with oxidative stress are quite relevant, as most of the prioritized vascular cognitive genes and proteins were enriched in these pathways. Our analysis provided a short list of proteins that could be contributing to VCI: DOLK, TSC1, ATP1A1, MAPK14, YWHAZ, CREB3, HSPB1, PRDX6, and LMNA. Moreover, our experimental results suggest a high implication of glycative stress, generating oxidative processes and post-translational protein modifications through advanced glycation end-products (AGEs). We propose that these products interact with their specific receptors (RAGE) and Notch signaling to contribute to the etiology of VCI.

Identification of potential therapeutic targets for COVID-19 through a structural-based similarity approach between SARS-CoV-2 and its human host proteins

Background: The COVID-19 pandemic caused by SARS-CoV-2 has led to millions of deaths worldwide, and vaccination efficacy has been decreasing with each lineage, necessitating the need for alternative antiviral therapies. Predicting host-virus protein-protein interactions (HV-PPIs) is essential for identifying potential host-targeting drug targets against SARS-CoV-2 infection. Objective: This study aims to identify therapeutic target proteins in humans that could act as virus-host-targeting drug targets against SARS-CoV-2 and study their interaction against antiviral inhibitors. Methods: A structure-based similarity approach was used to predict human proteins similar to SARS-CoV-2 ("hCoV-2"), followed by identifying PPIs between hCoV-2 and its target human proteins. Overlapping genes were identified between the protein-coding genes of the target and COVID-19-infected patient's mRNA expression data. Pathway and Gene Ontology (GO) term analyses, the construction of PPI networks, and the detection of hub gene modules were performed. Structure-based virtual screening with antiviral compounds was performed to identify potential hits against target gene-encoded protein. Results: This study predicted 19,051 unique target human proteins that interact with hCoV-2, and compared to the microarray dataset, 1,120 target and infected group differentially expressed genes (TIG-DEGs) were identified. The significant pathway and GO enrichment analyses revealed the involvement of these genes in several biological processes and molecular functions. PPI network analysis identified a significant hub gene with maximum neighboring partners. Virtual screening analysis identified three potential antiviral compounds against the target gene-encoded protein. Conclusion: This study provides potential targets for host-targeting drug development against SARS-CoV-2 infection, and further experimental validation of the target protein is required for pharmaceutical intervention.

Lifelong exercise training promotes the remodelling of the immune system and prostate signalome in a rat model of prostate carcinogenesis

This work aimed to understand how lifelong exercise training promotes the remodelling of the immune system and prostate signalome in a rat model of PCa. Fifty-five male Wistar rats were divided into four groups: control sedentary, control exercised, induced PCa sedentary and induced PCa exercised. Exercised animals were trained in a treadmill for 53 weeks. Pca induction consisted on the sequential administration of flutamide, N-methyl-N-nitrosourea and testosterone propionate implants. Serum concentrations of C-reactive protein (CRP) and tumor necrosis factor (TNF)-like weak inducer of apoptosis (TWEAK) were not different among groups. Peripheral levels of γδ T cells were higher in Pca exercised group than in the PCa sedentary group (p < 0.05). Exercise training also induced Oestrogen Receptor (ESR1) upregulation and Mitogen-activated Protein Kinase 13 (MAPK13) downregulation, changed the content of the phosphorylated (at Ser-104) form of this receptor (coded by the gene ESR1) and seemed to increase Erα phosphorylation and activity in exercised PCa rats when compared with sedentary PCa rats. Our data highlight the exercise-induced remodelling of peripheral lymphocyte subpopulations and lymphocyte infiltration in prostate tissue. Moreover, exercise training promotes the remodelling prostate signalome in this rat model of prostate carcinogenesis.

Selection of collagen IV fragments forming the outer sphere of the native protein: Assessment of biological activity for regenerative medicine

The aim of this research was to select the fragments that make up the outer layer of the collagen IV (COL4A6) protein and to assess their potential usefulness for regenerative medicine. It was expected that because protein-protein interactions take place via contact between external domains, the set of peptides forming the outer sphere of collagen IV will determine its interaction with other proteins. Cellulose-immobilized protein fragment libraries treated with polyclonal anti-collagen IV antibodies were used to select the peptides forming the outer sphere of collagen IV. In the first test, 33 peptides that strongly interacted with the polyclonal anti-collagen IV antibodies were selected from a library of non-overlapping fragments of collagen IV. The selected fragments of collagen IV (cleaved from the cellulose matrix) were tested for their cytotoxicity, their effects on cell viability and proliferation, and their impact on the formation of reactive oxygen species (ROS). The studies used RAW 264.7 mouse macrophage cells and Hs 680.Tr human fibroblasts. PrestoBlue, ToxiLight™, and ToxiLight 100% Lysis Control assays were conducted. The viability of fibroblasts cultured with the addition of increasing concentrations of the peptide mix did not show statistically significant differences from the control. Fragments 161-170, 221-230, 721-730, 1331-1340, 1521-1530, and 1661-1670 of COL4A6 were examined for cytotoxicity against BJ normal human foreskin fibroblasts. None of the collagen fragments were found to be cytotoxic. Further research is underway on the potential uses of collagen IV fragments in regenerative medicine.

AGF-PPIS: A protein-protein interaction site predictor based on an attention mechanism and graph convolutional networks

Protein-protein interactions play an important role in various biological processes. Interaction among proteins has a wide range of applications. Therefore, the correct identification of protein-protein interactions sites is crucial. In this paper, we propose a novel predictor for protein-protein interactions sites, AGF-PPIS, where we utilize a multi-head self-attention mechanism (introducing a graph structure), graph convolutional network, and feed-forward neural network. We use the Euclidean distance between each protein residue to generate the corresponding protein graph as the input of AGF-PPIS. On the independent test dataset Test_60, AGF-PPIS achieves superior performance over comparative methods in terms of seven different evaluation metrics (ACC, precision, recall, F1-score, MCC, AUROC, AUPRC), which fully demonstrates the validity and superiority of the proposed AGF-PPIS model. The source codes and the steps for usage of AGF-PPIS are available at https://github.com/fxh1001/AGF-PPIS.

miRNAs in Heart Development and Disease

Cardiovascular diseases (CVD) are a group of disorders that affect the heart and blood vessels. They include conditions such as myocardial infarction, coronary artery disease, heart failure, arrhythmia, and congenital heart defects. CVDs are the leading cause of death worldwide. Therefore, new medical interventions that aim to prevent, treat, or manage CVDs are of prime importance. MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression at the posttranscriptional level and play important roles in various biological processes, including cardiac development, function, and disease. Moreover, miRNAs can also act as biomarkers and therapeutic targets. In order to identify and characterize miRNAs and their target genes, scientists take advantage of computational tools such as bioinformatic algorithms, which can also assist in analyzing miRNA expression profiles, functions, and interactions in different cardiac conditions. Indeed, the combination of miRNA research and bioinformatic algorithms has opened new avenues for understanding and treating CVDs. In this review, we summarize the current knowledge on the roles of miRNAs in cardiac development and CVDs, discuss the challenges and opportunities, and provide some examples of recent bioinformatics for miRNA research in cardiovascular biology and medicine.

Emerging roles of the Protein Phosphatase 1 (PP1) in the context of viral infections

Protein Phosphatase 1 (PP1) is a major serine/threonine phosphatase in eukaryotes, participating in several cellular processes and metabolic pathways. Due to their low substrate specificity, PP1's catalytic subunits do not exist as free entities but instead bind to Regulatory Interactors of Protein Phosphatase One (RIPPO), which regulate PP1's substrate specificity and subcellular localization. Most RIPPOs bind to PP1 through combinations of short linear motifs (4-12 residues), forming highly specific PP1 holoenzymes. These PP1-binding motifs may, hence, represent attractive targets for the development of specific drugs that interfere with a subset of PP1 holoenzymes. Several viruses exploit the host cell protein (de)phosphorylation machinery to ensure efficient virus particle formation and propagation. While the role of many host cell kinases in viral life cycles has been extensively studied, the targeting of phosphatases by viral proteins has been studied in less detail. Here, we compile and review what is known concerning the role of PP1 in the context of viral infections and discuss how it may constitute a putative host-based target for the development of novel antiviral strategies.

Cracking the black box of deep sequence-based protein-protein interaction prediction

Identifying protein-protein interactions (PPIs) is crucial for deciphering biological pathways. Numerous prediction methods have been developed as cheap alternatives to biological experiments, reporting surprisingly high accuracy estimates. We systematically investigated how much reproducible deep learning models depend on data leakage, sequence similarities and node degree information, and compared them with basic machine learning models. We found that overlaps between training and test sets resulting from random splitting lead to strongly overestimated performances. In this setting, models learn solely from sequence similarities and node degrees. When data leakage is avoided by minimizing sequence similarities between training and test set, performances become random. Moreover, baseline models directly leveraging sequence similarity and network topology show good performances at a fraction of the computational cost. Thus, we advocate that any improvements should be reported relative to baseline methods in the future. Our findings suggest that predicting PPIs remains an unsolved task for proteins showing little sequence similarity to previously studied proteins, highlighting that further experimental research into the 'dark' protein interactome and better computational methods are needed.

Hubs and Bottlenecks in Protein-Protein Interaction Networks

Protein-protein interaction networks (PPINs) represent the physical interactions among proteins in a cell. These interactions are critical in all cellular processes, including signal transduction, metabolic regulation, and gene expression. In PPINs, centrality measures are widely used to identify the most critical nodes. The two most commonly used centrality measures in networks are degree and betweenness centralities. Degree centrality is the number of connections a node has in the network, and betweenness centrality is the measure of the extent to which a node lies on the shortest paths between pairs of other nodes in the network. In PPINs, proteins with high degree and betweenness centrality are referred to as hubs and bottlenecks respectively. Hubs and bottlenecks are topologically and functionally essential proteins that play crucial roles in maintaining the network's structure and function. This article comprehensively reviews essential literature on hubs and bottlenecks, including their properties and functions.

Unagi: Deep Generative Model for Deciphering Cellular Dynamics and In-Silico Drug Discovery in Complex Diseases

Human diseases are characterized by intricate cellular dynamics. Single-cell sequencing provides critical insights, yet a persistent gap remains in computational tools for detailed disease progression analysis and targeted in-silico drug interventions. Here, we introduce UNAGI, a deep generative neural network tailored to analyze time-series single-cell transcriptomic data. This tool captures the complex cellular dynamics underlying disease progression, enhancing drug perturbation modeling and discovery. When applied to a dataset from patients with Idiopathic Pulmonary Fibrosis (IPF), UNAGI learns disease-informed cell embeddings that sharpen our understanding of disease progression, leading to the identification of potential therapeutic drug candidates. Validation via proteomics reveals the accuracy of UNAGI's cellular dynamics analyses, and the use of the Fibrotic Cocktail treated human Precision-cut Lung Slices confirms UNAGI's predictions that Nifedipine, an antihypertensive drug, may have antifibrotic effects on human tissues. UNAGI's versatility extends to other diseases, including a COVID dataset, demonstrating adaptability and confirming its broader applicability in decoding complex cellular dynamics beyond IPF, amplifying its utility in the quest for therapeutic solutions across diverse pathological landscapes.

Intranuclear inclusions of polyQ-expanded ATXN1 sequester RNA molecules

Spinocerebellar ataxia type 1 (SCA1) is an autosomal dominant neurodegenerative disease caused by a trinucleotide (CAG) repeat expansion in the ATXN1 gene. It is characterized by the presence of polyglutamine (polyQ) intranuclear inclusion bodies (IIBs) within affected neurons. In order to investigate the impact of polyQ IIBs in SCA1 pathogenesis, we generated a novel protein aggregation model by inducible overexpression of the mutant ATXN1(Q82) isoform in human neuroblastoma SH-SY5Y cells. Moreover, we developed a simple and reproducible protocol for the efficient isolation of insoluble IIBs. Biophysical characterization showed that polyQ IIBs are enriched in RNA molecules which were further identified by next-generation sequencing. Finally, a protein interaction network analysis indicated that sequestration of essential RNA transcripts within ATXN1(Q82) IIBs may affect the ribosome resulting in error-prone protein synthesis and global proteome instability. These findings provide novel insights into the molecular pathogenesis of SCA1, highlighting the role of polyQ IIBs and their impact on critical cellular processes.

COVID-19: A novel holistic systems biology approach to predict its molecular mechanisms (in vitro) and repurpose drugs

PURPOSE

COVID-19 strangely kills some youth with no history of physical weakness, and in addition to the lungs, it may even directly harm other organs. Its complex mechanism has led to the loss of any significantly effective drug, and some patients with severe forms still die daily. Common methods for identifying disease mechanisms and drug design are often time-consuming or reductionist. Here, we use a novel holistic systems biology approach to predict its molecular mechanisms (in vitro), significant molecular relations with SARS, and repurpose drugs.

METHODS

We have utilized its relative phylogenic similarity to SARS. Using the available omics data for SARS and the fewer data for COVID-19 to decode the mechanisms and their significant relations, We applied the Cytoscape analyzer, MCODE, STRING, and DAVID tools to predict the topographically crucial molecules, clusters, protein interaction mappings, and functional analysis. We also applied a novel approach to identify the significant relations between the two infections using the Fischer exact test for MCODE clusters. We then constructed and analyzed a drug-gene network using PharmGKB and DrugBank (retrieved using the dgidb).

RESULTS

Some of the shared identified crucial molecules, BPs and pathways included Kaposi sarcoma-associated herpesvirus infection, Influenza A, and NOD-like receptor signaling pathways. Besides, our identified crucial molecules specific to host response against SARS-CoV-2 included FGA, BMP4, PRPF40A, and IFI16.

CONCLUSION

We also introduced seven new repurposed candidate drugs based on the drug-gene network analysis for the identified crucial molecules. Therefore, we suggest that our newly recommended repurposed drugs be further investigated in Vitro and in Vivo against COVID-19.

Artificial intelligence for neurodegenerative experimental models

INTRODUCTION

Experimental models are essential tools in neurodegenerative disease research. However, the translation of insights and drugs discovered in model systems has proven immensely challenging, marred by high failure rates in human clinical trials.

METHODS

Here we review the application of artificial intelligence (AI) and machine learning (ML) in experimental medicine for dementia research.

RESULTS

Considering the specific challenges of reproducibility and translation between other species or model systems and human biology in preclinical dementia research, we highlight best practices and resources that can be leveraged to quantify and evaluate translatability. We then evaluate how AI and ML approaches could be applied to enhance both cross-model reproducibility and translation to human biology, while sustaining biological interpretability.

DISCUSSION

AI and ML approaches in experimental medicine remain in their infancy. However, they have great potential to strengthen preclinical research and translation if based upon adequate, robust, and reproducible experimental data.

HIGHLIGHTS

There are increasing applications of AI in experimental medicine. We identified issues in reproducibility, cross-species translation, and data curation in the field. Our review highlights data resources and AI approaches as solutions. Multi-omics analysis with AI offers exciting future possibilities in drug discovery.

Structural coverage of the human interactome

Complex biological processes in cells are embedded in the interactome, representing the complete set of protein-protein interactions. Mapping and analyzing the protein structures are essential to fully comprehending these processes' molecular details. Therefore, knowing the structural coverage of the interactome is important to show the current limitations. Structural modeling of protein-protein interactions requires accurate protein structures. In this study, we mapped all experimental structures to the reference human proteome. Later, we found the enrichment in structural coverage when complementary methods such as homology modeling and deep learning (AlphaFold) were included. We then collected the interactions from the literature and databases to form the reference human interactome, resulting in 117 897 non-redundant interactions. When we analyzed the structural coverage of the interactome, we found that the number of experimentally determined protein complex structures is scarce, corresponding to 3.95% of all binary interactions. We also analyzed known and modeled structures to potentially construct the structural interactome with a docking method. Our analysis showed that 12.97% of the interactions from HuRI and 73.62% and 32.94% from the filtered versions of STRING and HIPPIE could potentially be modeled with high structural coverage or accuracy, respectively. Overall, this paper provides an overview of the current state of structural coverage of the human proteome and interactome.

Neurodevelopmental disorders and cancer networks share pathways, but differ in mechanisms, signaling strength, and outcome

Epidemiological studies suggest that individuals with neurodevelopmental disorders (NDDs) are more prone to develop certain types of cancer. Notably, however, the case statistics can be impacted by late discovery of cancer in individuals afflicted with NDDs, such as intellectual disorders, autism, and schizophrenia, which may bias the numbers. As to NDD-associated mutations, in most cases, they are germline while cancer mutations are sporadic, emerging during life. However, somatic mosaicism can spur NDDs, and cancer-related mutations can be germline. NDDs and cancer share proteins, pathways, and mutations. Here we ask (i) exactly which features they share, and (ii) how, despite their commonalities, they differ in clinical outcomes. To tackle these questions, we employed a statistical framework followed by network analysis. Our thorough exploration of the mutations, reconstructed disease-specific networks, pathways, and transcriptome levels and profiles of autism spectrum disorder (ASD) and cancers, point to signaling strength as the key factor: strong signaling promotes cell proliferation in cancer, and weaker (moderate) signaling impacts differentiation in ASD. Thus, we suggest that signaling strength, not activating mutations, can decide clinical outcome.

An androgen receptor regulated gene score is associated with epithelial to mesenchymal transition features in triple negative breast cancers

BACKGROUND

Androgen receptor (AR) is considered a marker of better prognosis in hormone receptor positive breast cancers (BC), however, its role in triple negative breast cancer (TNBC) is controversial. This may be attributed to intrinsic molecular differences or scoring methods for AR positivity. We derived AR regulated gene score and examined its utility in BC subtypes.

METHODS

AR regulated genes were derived by applying a bioinformatic pipeline on publicly available microarray data sets of AR+ BC cell lines and gene score was calculated as average expression of six AR regulated genes. Tumors were divided into AR high and low based on gene score and associations with clinical parameters, circulating androgens, survival and epithelial to mesenchymal transition (EMT) markers were examined, further evaluated in invitro models and public datasets.

RESULTS

53% (133/249) tumors were classified as AR gene score high and were associated with significantly better clinical parameters, disease-free survival (86.13 vs 72.69 months, log rank p = 0.032) when compared to AR low tumors. 36% of TNBC (N = 66) were AR gene score high with higher expression of EMT markers (p = 0.024) and had high intratumoral levels of 5α-reductase, enzyme involved in intracrine androgen metabolism. In MDA-MB-453 treated with dihydrotestosterone, SLUG expression increased, E-cadherin decreased with increase in migration and these changes were reversed with bicalutamide. Similar results were obtained in public datasets.

CONCLUSION

Deciphering the role of AR in BC is difficult based on AR protein levels alone. Our results support the context dependent function of AR in driving better prognosis in ER positive tumors and EMT features in TNBC tumors.

Low-dose Bisphenol A exposure alters the functionality and cellular environment in a human cardiomyocyte model

Early embryonic development represents a sensitive time-window during which the foetus might be vulnerable to the exposure of environmental contaminants, potentially leading to heart diseases also later in life. Bisphenol A (BPA), a synthetic chemical widely used in plastics manufacturing, has been associated with heart developmental defects, even in low concentrations. This study aims to investigate the effects of environmentally relevant doses of BPA on developing cardiomyocytes using a human induced pluripotent stem cell (hiPSC)-derived model. Firstly, a 2D in vitro differentiation system to obtain cardiomyocytes from hiPSCs (hiPSC-CMs) have been established and characterised to provide a suitable model for the early stages of cardiac development. Then, the effects of a repeated BPA exposure, starting from the undifferentiated stage throughout the differentiation process, were evaluated. The chemical significantly decreased the beat rate of hiPSC-CMs, extending the contraction and relaxation time in a dose-dependent manner. Quantitative proteomics analysis revealed a high abundance of basement membrane (BM) components (e.g., COL4A1, COL4A2, LAMC1, NID2) and a significant increase in TNNC1 and SERBP1 proteins in hiPSC-CMs treated with BPA. Network analysis of proteomics data supported altered extracellular matrix remodelling and provided a disease-gene association with well-known pathological conditions of the heart. Furthermore, upon hypoxia-reoxygenation challenge, hiPSC-CMs treated with BPA showed higher rate of apoptotic events. Taken together, our results revealed that a long-term treatment, even with low doses of BPA, interferes with hiPSC-CMs functionality and alters the surrounding cellular environment, providing new insights about diseases that might arise upon the toxin exposure. Our study contributes to the current understanding of BPA effects on developing human foetal cardiomyocytes, in correlation with human clinical observations and animal studies, and it provides a suitable model for New Approach Methodologies (NAMs) for environmental chemical hazard and risk assessment.

DensePPI: A Novel Image-Based Deep Learning Method for Prediction of Protein-Protein Interactions

Protein-protein interactions (PPI) are crucial for understanding the behaviour of living organisms and identifying disease associations. This paper proposes DensePPI, a novel deep convolution strategy applied to the 2D image map generated from the interacting protein pairs for PPI prediction. A colour encoding scheme has been introduced to embed the bigram interaction possibilities of Amino Acids into RGB colour space to enhance the learning and prediction task. The DensePPI model is trained on 5.5 million sub-images of size 128×128 generated from nearly 36,000 interacting and 36,000 non-interacting benchmark protein pairs. The performance is evaluated on independent datasets from five different organisms; Caenorhabditis elegans, Escherichia coli, Helicobacter Pylori, Homo sapiens and Mus Musculus. The proposed model achieves an average prediction accuracy score of 99.95% on these datasets, considering inter-species and intra-species interactions. The performance of DensePPI is compared with the state-of-the-art methods and outperforms those approaches in different evaluation metrics. Improved performance of DensePPI indicates the efficiency of the image-based encoding strategy of sequence information with the deep learning architecture in PPI prediction. The enhanced performance on diverse test sets shows that the DensePPI is significant for intra-species interaction prediction and cross-species interactions. The dataset, supplementary file, and the developed models are available at https://github.com/Aanzil/DensePPI for academic use only.

Construction of HBV-HCC prognostic model and immune characteristics based on potential genes mining through protein interaction networks

OBJECTIVE

To search for human protein-coding genes related to hepatocellular carcinoma (HCC) in the context of hepatitis B virus (HBV) infection, and perform prognosis risk assessment.

METHODS

Genes related to HBV-HCC were selected through literature screening and protein-protein interaction (PPI) network database analysis. Prognosis potential genes (PPGs) were identified using Cox regression analysis. Patients were divided into high-risk and low-risk groups based on PPGs, and risk scores were calculated. Kaplan-Meier plots were used to analyze overall survival rates, and the results were predicted based on clinicopathological variables. Association analysis was also conducted with immune infiltration, immune therapy, and drug sensitivity. Experimental verification of the expression of PPGs was done in patient liver cancer tissue and normal liver tissue adjacent to tumors.

RESULTS

The use of a prognosis potential genes risk assessment model can reliably predict the prognosis risk of patients, demonstrating strong predictive ability. Kaplan-Meier analysis showed that the overall survival rate of the low-risk group was significantly higher than that of the high-risk group. There were significant differences between the two subgroups in terms of immune infiltration and IC50 association analysis. Experimental verification revealed that CYP2C19, FLNC, and HNRNPC were highly expressed in liver cancer tissue, while UBE3A was expressed at a lower level.

CONCLUSION

PPGs can be used to predict the prognosis risk of HBV-HCC patients and play an important role in the diagnosis and treatment of liver cancer. They also reveal their potential role in the tumor immune microenvironment, clinical-pathological characteristics, and prognosis.

Plasma proteomic associations with genetics and health in the UK Biobank

The Pharma Proteomics Project is a precompetitive biopharmaceutical consortium characterizing the plasma proteomic profiles of 54,219 UK Biobank participants. Here we provide a detailed summary of this initiative, including technical and biological validations, insights into proteomic disease signatures, and prediction modelling for various demographic and health indicators. We present comprehensive protein quantitative trait locus (pQTL) mapping of 2,923 proteins that identifies 14,287 primary genetic associations, of which 81% are previously undescribed, alongside ancestry-specific pQTL mapping in non-European individuals. The study provides an updated characterization of the genetic architecture of the plasma proteome, contextualized with projected pQTL discovery rates as sample sizes and proteomic assay coverages increase over time. We offer extensive insights into trans pQTLs across multiple biological domains, highlight genetic influences on ligand-receptor interactions and pathway perturbations across a diverse collection of cytokines and complement networks, and illustrate long-range epistatic effects of ABO blood group and FUT2 secretor status on proteins with gastrointestinal tissue-enriched expression. We demonstrate the utility of these data for drug discovery by extending the genetic proxied effects of protein targets, such as PCSK9, on additional endpoints, and disentangle specific genes and proteins perturbed at loci associated with COVID-19 susceptibility. This public-private partnership provides the scientific community with an open-access proteomics resource of considerable breadth and depth to help to elucidate the biological mechanisms underlying proteo-genomic discoveries and accelerate the development of biomarkers, predictive models and therapeutics1.

Protein network analysis links the NSL complex to Parkinson's disease via mitochondrial and nuclear biology

Whilst the majority of Parkinson's Disease (PD) cases are sporadic, much of our understanding of the pathophysiological basis of the disease can be traced back to the study of rare, monogenic forms of PD. In the past decade, the availability of genome-wide association studies (GWAS) has facilitated a shift in focus, toward identifying common risk variants conferring increased risk of developing PD across the population. A recent mitophagy screening assay of GWAS candidates has functionally implicated the non-specific lethal (NSL) complex in the regulation of PINK1-mitophagy. Here, a bioinformatics approach has been taken to investigate the proteome of the NSL complex, to unpick its relevance to PD pathogenesis. The NSL interactome has been built, using 3 online tools: PINOT, HIPPIE and MIST, to mine curated, literature-derived protein-protein interaction (PPI) data. We built (i) the 'mitochondrial' NSL interactome exploring its relevance to PD genetics and (ii) the PD-oriented NSL interactome to uncover biological pathways underpinning the NSL/PD association. In this study, we find the mitochondrial NSL interactome to be significantly enriched for the protein products of PD-associated genes, including the Mendelian PD genes LRRK2 and VPS35. In addition, we find nuclear processes to be amongst those most significantly enriched within the PD-associated NSL interactome. These findings strengthen the role of the NSL complex in sporadic and familial PD, mediated by both its mitochondrial and nuclear functions.

LEVELNET to visualize, explore, and compare protein-protein interaction networks

Physical interactions between proteins are central to all biological processes. Yet, the current knowledge of who interacts with whom in the cell and in what manner relies on partial, noisy, and highly heterogeneous data. Thus, there is a need for methods comprehensively describing and organizing such data. LEVELNET is a versatile and interactive tool for visualizing, exploring, and comparing protein-protein interaction (PPI) networks inferred from different types of evidence. LEVELNET helps to break down the complexity of PPI networks by representing them as multi-layered graphs and by facilitating the direct comparison of their subnetworks toward biological interpretation. It focuses primarily on the protein chains whose 3D structures are available in the Protein Data Bank. We showcase some potential applications, such as investigating the structural evidence supporting PPIs associated to specific biological processes, assessing the co-localization of interaction partners, comparing the PPI networks obtained through computational experiments versus homology transfer, and creating PPI benchmarks with desired properties.

AutoCore: A network-based definition of the core module of human autoimmunity and autoinflammation

Although research on rare autoimmune and autoinflammatory diseases has enabled definition of nonredundant regulators of homeostasis in human immunity, because of the single gene-single disease nature of many of these diseases, contributing factors were mostly unveiled in sequential and noncoordinated individual studies. We used a network-based approach for integrating a set of 186 inborn errors of immunity with predominant autoimmunity/autoinflammation into a comprehensive map of human immune dysregulation, which we termed "AutoCore." The AutoCore is located centrally within the interactome of all protein-protein interactions, connecting and pinpointing multidisease markers for a range of common, polygenic autoimmune/autoinflammatory diseases. The AutoCore can be subdivided into 19 endotypes that correspond to molecularly and phenotypically cohesive disease subgroups, providing a molecular mechanism-based disease classification and rationale toward systematic targeting for therapeutic purposes. Our study provides a proof of concept for using network-based methods to systematically investigate the molecular relationships between individual rare diseases and address a range of conceptual, diagnostic, and therapeutic challenges.

Proteomics and disease network associations evaluation of environmentally relevant Bisphenol A concentrations in a human 3D neural stem cell model

Bisphenol A (BPA) exposure is associated with a plethora of neurodevelopmental abnormalities and brain disorders. Previous studies have demonstrated BPA-induced perturbations to critical neural stem cell (NSC) characteristics, such as proliferation and differentiation, although the underlying molecular mechanisms remain under debate. The present study evaluated the effects of a repeated-dose exposure of environmentally relevant BPA concentrations during the in vitro 3D neural induction of human induced pluripotent stem cells (hiPSCs), emulating a chronic exposure scenario. Firstly, we demonstrated that our model is suitable for NSC differentiation during the early stages of embryonic brain development. Our morphological image analysis showed that BPA exposure at 0.01, 0.1 and 1 µM decreased the average spheroid size by day 21 (D21) of the neural induction, while no effect on cell viability was detected. No alteration to the rate of the neural induction was observed based on the expression of key neural lineage and neuroectodermal transcripts. Quantitative proteomics at D21 revealed several differentially abundant proteins across all BPA-treated groups with important functions in NSC proliferation and maintenance (e.g., FABP7, GPC4, GAP43, Wnt-8B, TPPP3). Additionally, a network analysis demonstrated alterations to the glycolytic pathway, potentially implicating BPA-induced changes to glycolytic signalling in NSC proliferation impairments, as well as the pathophysiology of brain disorders including intellectual disability, autism spectrum disorders, and amyotrophic lateral sclerosis (ALS). This study enhances the current understanding of BPA-related NSC aberrations based mostly on acute, often high dose exposures of rodent in vivo and in vitro models and human GWAS data in a novel human 3D cell-based model with real-life scenario relevant prolonged and low-level exposures, offering further mechanistic insights into the ramifications of BPA exposure on the developing human brain and consequently, later life neurological disorders.

Transcription factors organize into functional groups on the linear genome and in 3D chromatin

Transcription factors (TFs) and their binding sites have evolved to interact cooperatively or competitively with each other. Here we examine in detail, across multiple cell lines, such cooperation or competition among TFs both in sequential and spatial proximity (using chromatin conformation capture assays), considering in vivo binding data as well as TF binding motifs in DNA. We ascertain significantly co-occurring ("attractive") or avoiding ("repulsive") TF pairs using robust randomized models that retain the essential characteristics of the experimental data. Across human cell lines TFs organize into two groups, with intra-group attraction and inter-group repulsion. This is true for both sequential and spatial proximity, and for both in vivo binding and sequence motifs. Attractive TF pairs exhibit significantly more physical interactions suggesting an underlying mechanism. The two TF groups differ significantly in their genomic and network properties, as well in their function-while one group regulates housekeeping function, the other potentially regulates lineage-specific functions, that are disrupted in cancer. Weaker binding sites tend to occur in spatially interacting regions of the genome. Our results suggest that a complex pattern of spatial cooperativity of TFs and chromatin has evolved with the genome to support housekeeping and lineage-specific functions.

Uncharacterized Proteins CxORFx: Subinteractome Analysis and Prognostic Significance in Cancers

Functions of about 10% of all the proteins and their associations with diseases are poorly annotated or not annotated at all. Among these proteins, there is a group of uncharacterized chromosome-specific open-reading frame genes (CxORFx) from the 'Tdark' category. The aim of the work was to reveal associations of CxORFx gene expression and ORF proteins' subinteractomes with cancer-driven cellular processes and molecular pathways. We performed systems biology and bioinformatic analysis of 219 differentially expressed CxORFx genes in cancers, an estimation of prognostic significance of novel transcriptomic signatures and analysis of subinteractome composition using several web servers (GEPIA2, KMplotter, ROC-plotter, TIMER, cBioPortal, DepMap, EnrichR, PepPSy, cProSite, WebGestalt, CancerGeneNet, PathwAX II and FunCoup). The subinteractome of each ORF protein was revealed using ten different data sources on physical protein-protein interactions (PPIs) to obtain representative datasets for the exploration of possible cellular functions of ORF proteins through a spectrum of neighboring annotated protein partners. A total of 42 out of 219 presumably cancer-associated ORF proteins and 30 cancer-dependent binary PPIs were found. Additionally, a bibliometric analysis of 204 publications allowed us to retrieve biomedical terms related to ORF genes. In spite of recent progress in functional studies of ORF genes, the current investigations aim at finding out the prognostic value of CxORFx expression patterns in cancers. The results obtained expand the understanding of the possible functions of the poorly annotated CxORFx in the cancer context.

Network medicine framework identified drug-repurposing opportunities of pharmaco-active compounds of Angelica acutiloba (Siebold & Zucc.) Kitag. for skin aging

Increasing incidence of skin aging has highlighted the importance of identifying effective drugs with repurposed opportunities for skin aging. We aimed to identify pharmaco-active compounds with drug-repurposing opportunities for skin aging from Angelica acutiloba (Siebold & Zucc.) Kitag. (AAK). The proximity of network medicine framework (NMF) firstly identified 8 key AAK compounds with repurposed opportunities for skin aging, which may exert by regulating 29 differentially expressed genes (DGEs) of skin aging, including 13 up-regulated targets and 16 down-regulated targets. Connectivity MAP (cMAP) analysis revealed 8 key compounds were involved in regulating the process of cell proliferation and apoptosis, mitochondrial energy metabolism and oxidative stress of skin aging. Molecular docking analysis showed that 8 key compounds had a high docked ability with AR, BCHE, HPGD and PI3, which were identified as specific biomarker for the diagnosis of skin aging. Finally, the mechanisms of these key compounds were predicted to be involved in inhibiting autophagy pathway and activating Phospholipase D signaling pathway. In conclusion, this study firstly elucidated the drug-repurposing opportunities of AAK compounds for skin aging, providing a theoretical reference for identifying repurposing drugs from Chinese medicine and new insights for our future research.

On the identification of potential novel therapeutic targets for spinocerebellar ataxia type 1 (SCA1) neurodegenerative disease using EvoPPI3

EvoPPI (http://evoppi.i3s.up.pt), a meta-database for protein-protein interactions (PPI), has been upgraded (EvoPPI3) to accept new types of data, namely, PPI from patients, cell lines, and animal models, as well as data from gene modifier experiments, for nine neurodegenerative polyglutamine (polyQ) diseases caused by an abnormal expansion of the polyQ tract. The integration of the different types of data allows users to easily compare them, as here shown for Ataxin-1, the polyQ protein involved in spinocerebellar ataxia type 1 (SCA1) disease. Using all available datasets and the data here obtained for Drosophila melanogaster wt and exp Ataxin-1 mutants (also available at EvoPPI3), we show that, in humans, the Ataxin-1 network is much larger than previously thought (380 interactors), with at least 909 interactors. The functional profiling of the newly identified interactors is similar to the ones already reported in the main PPI databases. 16 out of 909 interactors are putative novel SCA1 therapeutic targets, and all but one are already being studied in the context of this disease. The 16 proteins are mainly involved in binding and catalytic activity (mainly kinase activity), functional features already thought to be important in the SCA1 disease.

Structure Learning for Gene Regulatory Networks

Inference of biological network structures is often performed on high-dimensional data, yet is hindered by the limited sample size of high throughput "omics" data typically available. To overcome this challenge, often referred to as the "small n, large p problem," we exploit known organizing principles of biological networks that are sparse, modular, and likely share a large portion of their underlying architecture. We present SHINE-Structure Learning for Hierarchical Networks-a framework for defining data-driven structural constraints and incorporating a shared learning paradigm for efficiently learning multiple Markov networks from high-dimensional data at large p/n ratios not previously feasible. We evaluated SHINE on Pan-Cancer data comprising 23 tumor types, and found that learned tumor-specific networks exhibit expected graph properties of real biological networks, recapture previously validated interactions, and recapitulate findings in literature. Application of SHINE to the analysis of subtype-specific breast cancer networks identified key genes and biological processes for tumor maintenance and survival as well as potential therapeutic targets for modulating known breast cancer disease genes.