PICKLE 3.0: Enriching the human Meta-database with the mouse protein interactome extended via mouse-human orthology

GLIMS: A two-stage gradual-learning method for cancer genes prediction using multi-omics data and co-splicing network

Identifying cancer genes is vital for cancer diagnosis and treatment. However, because of the complexity of cancer occurrence and limited cancer genes knowledge, it is hard to identify cancer genes accurately using only a few omics data, and the overall performance of existing methods is being called for further improvement. Here, we introduce a two-stage gradual-learning strategy GLIMS to predict cancer genes using integrative features from multi-omics data. Firstly, it uses a semi-supervised hierarchical graph neural network to predict the initial candidate cancer genes by integrating multi-omics data and protein-protein interaction (PPI) network. Then, it uses an unsupervised approach to further optimize the initial prediction by integrating the co-splicing network in post-transcriptional regulation, which plays an important role in cancer development. Systematic experiments on multi-omics cancer data demonstrated that GLIMS outperforms the state-of-the-art methods for the identification of cancer genes and it could be a useful tool to help advance cancer analysis.

Extracellular Matrix Interactome in Modulating Vascular Homeostasis and Remodeling

The ECM (extracellular matrix) is a major component of the vascular microenvironment that modulates vascular homeostasis. ECM proteins include collagens, elastin, noncollagen glycoproteins, and proteoglycans/glycosaminoglycans. ECM proteins form complex matrix structures, such as the basal lamina and collagen and elastin fibers, through direct interactions or lysyl oxidase-mediated cross-linking. Moreover, ECM proteins directly interact with cell surface receptors or extracellular secreted molecules, exerting matricellular and matricrine modulation, respectively. In addition, extracellular proteases degrade or cleave matrix proteins, thereby contributing to ECM turnover. These interactions constitute the ECM interactome network, which is essential for maintaining vascular homeostasis and preventing pathological vascular remodeling. The current review mainly focuses on endogenous matrix proteins in blood vessels and discusses the interaction of these matrix proteins with other ECM proteins, cell surface receptors, cytokines, complement and coagulation factors, and their potential roles in maintaining vascular homeostasis and preventing pathological remodeling.

Differential protein-protein interactions underlie signaling mediated by the TCR and a 4-1BB domain-containing CAR

Cells rely on activity-dependent protein-protein interactions to convey biological signals. For chimeric antigen receptor (CAR) T cells containing a 4-1BB costimulatory domain, receptor engagement is thought to stimulate the formation of protein complexes similar to those stimulated by T cell receptor (TCR)-mediated signaling, but the number and type of protein interaction-mediating binding domains differ between CARs and TCRs. Here, we performed coimmunoprecipitation mass spectrometry analysis of a second-generation, CD19-directed 4-1BB:ζ CAR (referred to as bbζCAR) and identified 128 proteins that increased their coassociation after target engagement. We compared activity-induced TCR and CAR signalosomes by quantitative multiplex coimmunoprecipitation and showed that bbζCAR engagement led to the activation of two modules of protein interactions, one similar to TCR signaling that was more weakly engaged by bbζCAR as compared with the TCR and one composed of TRAF signaling complexes that was not engaged by the TCR. Batch-to-batch and interindividual variations in production of the cytokine IL-2 correlated with differences in the magnitude of protein network activation. Future CAR T cell manufacturing protocols could measure, and eventually control, biological variation by monitoring these signalosome activation markers.

Protein-protein interaction network-based integration of GWAS and functional data for blood pressure regulation analysis

BACKGROUND

It is valuable to analyze the genome-wide association studies (GWAS) data for a complex disease phenotype in the context of the protein-protein interaction (PPI) network, as the related pathophysiology results from the function of interacting polyprotein pathways. The analysis may include the design and curation of a phenotype-specific GWAS meta-database incorporating genotypic and eQTL data linking to PPI and other biological datasets, and the development of systematic workflows for PPI network-based data integration toward protein and pathway prioritization. Here, we pursued this analysis for blood pressure (BP) regulation.

METHODS

The relational scheme of the implemented in Microsoft SQL Server BP-GWAS meta-database enabled the combined storage of: GWAS data and attributes mined from GWAS Catalog and the literature, Ensembl-defined SNP-transcript associations, and GTEx eQTL data. The BP-protein interactome was reconstructed from the PICKLE PPI meta-database, extending the GWAS-deduced network with the shortest paths connecting all GWAS-proteins into one component. The shortest-path intermediates were considered as BP-related. For protein prioritization, we combined a new integrated GWAS-based scoring scheme with two network-based criteria: one considering the protein role in the reconstructed by shortest-path (RbSP) interactome and one novel promoting the common neighbors of GWAS-prioritized proteins. Prioritized proteins were ranked by the number of satisfied criteria.

RESULTS

The meta-database includes 6687 variants linked with 1167 BP-associated protein-coding genes. The GWAS-deduced PPI network includes 1065 proteins, with 672 forming a connected component. The RbSP interactome contains 1443 additional, network-deduced proteins and indicated that essentially all BP-GWAS proteins are at most second neighbors. The prioritized BP-protein set was derived from the union of the most BP-significant by any of the GWAS-based or the network-based criteria. It included 335 proteins, with ~ 2/3 deduced from the BP PPI network extension and 126 prioritized by at least two criteria. ESR1 was the only protein satisfying all three criteria, followed in the top-10 by INSR, PTN11, CDK6, CSK, NOS3, SH2B3, ATP2B1, FES and FINC, satisfying two. Pathway analysis of the RbSP interactome revealed numerous bioprocesses, which are indeed functionally supported as BP-associated, extending our understanding about BP regulation.

CONCLUSIONS

The implemented workflow could be used for other multifactorial diseases.

Deciphering ligand-receptor-mediated intercellular communication based on ensemble deep learning and the joint scoring strategy from single-cell transcriptomic data

BACKGROUND

Cell-cell communication in a tumor microenvironment is vital to tumorigenesis, tumor progression and therapy. Intercellular communication inference helps understand molecular mechanisms of tumor growth, progression and metastasis.

METHODS

Focusing on ligand-receptor co-expressions, in this study, we developed an ensemble deep learning framework, CellComNet, to decipher ligand-receptor-mediated cell-cell communication from single-cell transcriptomic data. First, credible LRIs are captured by integrating data arrangement, feature extraction, dimension reduction, and LRI classification based on an ensemble of heterogeneous Newton boosting machine and deep neural network. Next, known and identified LRIs are screened based on single-cell RNA sequencing (scRNA-seq) data in certain tissues. Finally, cell-cell communication is inferred by incorporating scRNA-seq data, the screened LRIs, a joint scoring strategy that combines expression thresholding and expression product of ligands and receptors.

RESULTS

The proposed CellComNet framework was compared with four competing protein-protein interaction prediction models (PIPR, XGBoost, DNNXGB, and OR-RCNN) and obtained the best AUCs and AUPRs on four LRI datasets, elucidating the optimal LRI classification ability. CellComNet was further applied to analyze intercellular communication in human melanoma and head and neck squamous cell carcinoma (HNSCC) tissues. The results demonstrate that cancer-associated fibroblasts highly communicate with melanoma cells and endothelial cells strong communicate with HNSCC cells.

CONCLUSIONS

The proposed CellComNet framework efficiently identified credible LRIs and significantly improved cell-cell communication inference performance. We anticipate that CellComNet can contribute to anticancer drug design and tumor-targeted therapy.

Network-Based Prediction of Side Effects of Repurposed Antihypertensive Sartans against COVID-19 via Proteome and Drug-Target Interactomes

The potential of targeting the Renin-Angiotensin-Aldosterone System (RAAS) as a treatment for the coronavirus disease 2019 (COVID-19) is currently under investigation. One way to combat this disease involves the repurposing of angiotensin receptor blockers (ARBs), which are antihypertensive drugs, because they bind to angiotensin-converting enzyme 2 (ACE2), which in turn interacts with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein. However, there has been no in silico analysis of the potential toxicity risks associated with the use of these drugs for the treatment of COVID-19. To address this, a network-based bioinformatics methodology was used to investigate the potential side effects of known Food and Drug Administration (FDA)-approved antihypertensive drugs, Sartans. This involved identifying the human proteins targeted by these drugs, their first neighbors, and any drugs that bind to them using publicly available experimentally supported data, and subsequently constructing proteomes and protein-drug interactomes. This methodology was also applied to Pfizer's Paxlovid, an antiviral drug approved by the FDA for emergency use in mild-to-moderate COVID-19 treatment. The study compares the results for both drug categories and examines the potential for off-target effects, undesirable involvement in various biological processes and diseases, possible drug interactions, and the potential reduction in drug efficiency resulting from proteoform identification.

Reconstruction of a Comprehensive Interactome and Experimental Data Analysis of FRA10AC1 May Provide Insights into Its Biological Role in Health and Disease

FRA10AC1, the causative gene for the manifestation of the FRA10A fragile site, encodes a well-conserved nuclear protein characterized as a non-core spliceosomal component. Pre-mRNA splicing perturbations have been linked with neurodevelopmental diseases. FRA10AC1 variants have been, recently, causally linked with severe neuropathological and growth retardation phenotypes. To further elucidate the participation of FRA10AC1 in spliceosomal multiprotein complexes and its involvement in neurological phenotypes related to splicing, we exploited protein–protein interaction experimental data and explored network information and information deduced from transcriptomics. We confirmed the direct interaction of FRA10AC1with ESS2, a non-core spliceosomal protein, mapped their interacting domains, and documented their tissue co-localization and physical interaction at the level of intracellular protein stoichiometries. Although FRA10AC1 and SF3B2, a major core spliceosomal protein, were shown to interact under in vitro conditions, the endogenous proteins failed to co-immunoprecipitate. A reconstruction of a comprehensive, strictly binary, protein–protein interaction network of FRA10AC1 revealed dense interconnectivity with many disease-associated spliceosomal components and several non-spliceosomal regulatory proteins. The topological neighborhood of FRA10AC1 depicts an interactome associated with multiple severe monogenic and multifactorial neurodevelopmental diseases mainly referring to spliceosomopathies. Our results suggest that FRA10AC1 involvement in pre-mRNA processing might be strengthened by interconnecting splicing with transcription and mRNA export, and they propose the broader role(s) of FRA10AC1 in cell pathophysiology.

Protein network and pathway analysis in a pharmacogenetic study of cyclosporine treatment response in Greek patients with psoriasis

Although cyclosporine comprises a well-established systemic therapy for psoriasis, patients show important heterogeneity in their treatment response. The aim of our study was the pharmacogenetic analysis of 200 Greek patients with psoriasis based on the cyclosporine pathway related protein-protein interaction (PPI) network, reconstructed through the PICKLE meta-database. We genotyped 27 single nucleotide polymorphisms, mapped to 22 key protein nodes of the cyclosporine pathway, via the utilization of the iPLEX®GOLD panel of the MassARRAY® System. Single-SNP analyses showed statistically significant associations between CALM1 rs12885713 (P = 0.0108) and MALT1 rs2874116 (P = 0.0006) polymorphisms with positive response to cyclosporine therapy after correction for multiple comparisons, with the haplotype analyses further enhancing the predictive value of rs12885713 as a pharmacogenetic biomarker for cyclosporine therapy (P = 0.0173). Our findings have the potential to improve our prediction of cyclosporine efficacy and safety in psoriasis patients, as well as provide the framework for the pharmacogenetics of biological therapies in complex diseases.

MIR retrotransposons link the epigenome and the transcriptome of coding genes in acute myeloid leukemia

DNMT3A and IDH1/2 mutations combinatorically regulate the transcriptome and the epigenome in acute myeloid leukemia; yet the mechanisms of this interplay are unknown. Using a systems approach within topologically associating domains, we find that genes with significant expression-methylation correlations are enriched in signaling and metabolic pathways. The common denominator across these methylation-regulated genes is the density in MIR retrotransposons of their introns. Moreover, a discrete number of CpGs overlapping enhancers are responsible for regulating most of these genes. Established mouse models recapitulate the dependency of MIR-rich genes on the balanced expression of epigenetic modifiers, while projection of leukemic profiles onto normal hematopoiesis ones further consolidates the dependencies of methylation-regulated genes on MIRs. Collectively, MIR elements on genes and enhancers are susceptible to changes in DNA methylation activity and explain the cooperativity of proteins in this pathway in normal and malignant hematopoiesis.

How Far Are We from the Completion of the Human Protein Interactome Reconstruction?

After more than fifteen years from the first high-throughput experiments for human protein–protein interaction (PPI) detection, we are still wondering how close the completion of the genome-scale human PPI network reconstruction is, what needs to be further explored and whether the biological insights gained from the holistic investigation of the current network are valid and useful. The unique structure of PICKLE, a meta-database of the human experimentally determined direct PPI network developed by our group, presently covering ~80% of the UniProtKB/Swiss-Prot reviewed human complete proteome, enables the evaluation of the interactome expansion by comparing the successive PICKLE releases since 2013. We observe a gradual overall increase of 39%, 182%, and 67% in protein nodes, PPIs, and supporting references, respectively. Our results indicate that, in recent years, (a) the PPI addition rate has decreased, (b) the new PPIs are largely determined by high-throughput experiments and mainly concern existing protein nodes and (c), as we had predicted earlier, most of the newly added protein nodes have a low degree. These observations, combined with a largely overlapping k-core between PICKLE releases and a network density increase, imply that an almost complete picture of a structurally defined network has been reached. The comparative unsupervised application of two clustering algorithms indicated that exploring the full interactome topology can reveal the protein neighborhoods involved in closely related biological processes as transcriptional regulation, cell signaling and multiprotein complexes such as the connexon complex associated with cancers. A well-reconstructed human protein interactome is a powerful tool in network biology and medicine research forming the basis for multi-omic and dynamic analyses.

Screening and Comprehensive Analysis of Cancer-Associated tRNA-Derived Fragments

tRNA-derived fragments (tRFs) constitute a novel class of small non-coding RNA cleaved from tRNAs. In recent years, researches have shown the regulatory roles of a few tRFs in cancers, illuminating a new direction for tRF-centric cancer researches. Nonetheless, more specific screening of tRFs related to oncogenesis pathways, cancer progression stages and cancer prognosis is continuously demanded to reveal the landscape of the cancer-associated tRFs. In this work, by combining the clinical information recorded in The Cancer Genome Atlas (TCGA) and the tRF expression profiles curated by MINTbase v2.0, we systematically screened 1,516 cancer-associated tRFs (ca-tRFs) across seven cancer types. The ca-tRF set collectively combined the differentially expressed tRFs between cancer samples and control samples, the tRFs significantly correlated with tumor stage and the tRFs significantly correlated with patient survival. By incorporating our previous tRF-target dataset, we found the ca-tRFs tend to target cancer-associated genes and onco-pathways like ATF6-mediated unfolded protein response, angiogenesis, cell cycle process regulation, focal adhesion, PI3K-Akt signaling pathway, cellular senescence and FoxO signaling pathway across multiple cancer types. And cell composition analysis implies that the expressions of ca-tRFs are more likely to be correlated with T-cell infiltration. We also found the ca-tRF expression pattern is informative to prognosis, suggesting plausible tRF-based cancer subtypes. Together, our systematic analysis demonstrates the potentially extensive involvements of tRFs in cancers, and provides a reasonable list of cancer-associated tRFs for further investigations.

Mouse Genome Informatics (MGI): latest news from MGD and GXD

The Mouse Genome Informatics (MGI) database system combines multiple expertly curated community data resources into a shared knowledge management ecosystem united by common metadata annotation standards. MGI's mission is to facilitate the use of the mouse as an experimental model for understanding the genetic and genomic basis of human health and disease. MGI is the authoritative source for mouse gene, allele, and strain nomenclature and is the primary source of mouse phenotype annotations, functional annotations, developmental gene expression information, and annotations of mouse models with human diseases. MGI maintains mouse anatomy and phenotype ontologies and contributes to the development of the Gene Ontology and Disease Ontology and uses these ontologies as standard terminologies for annotation. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are MGI's two major knowledgebases. Here, we highlight some of the recent changes and enhancements to MGD and GXD that have been implemented in response to changing needs of the biomedical research community and to improve the efficiency of expert curation. MGI can be accessed freely at http://www.informatics.jax.org .

The Dynamic Network of RNP RNase P Subunits

Ribonuclease P (RNase P) is an important ribonucleoprotein (RNP), responsible for the maturation of the 5′ end of precursor tRNAs (pre-tRNAs). In all organisms, the cleavage activity of a single phosphodiester bond adjacent to the first nucleotide of the acceptor stem is indispensable for cell viability and lies within an essential catalytic RNA subunit. Although RNase P is a ribozyme, its kinetic efficiency in vivo, as well as its structural variability and complexity throughout evolution, requires the presence of one protein subunit in bacteria to several protein partners in archaea and eukaryotes. Moreover, the existence of protein-only RNase P (PRORP) enzymes in several organisms and organelles suggests a more complex evolutionary timeline than previously thought. Recent detailed structures of bacterial, archaeal, human and mitochondrial RNase P complexes suggest that, although apparently dissimilar enzymes, they all recognize pre-tRNAs through conserved interactions. Interestingly, individual protein subunits of the human nuclear and mitochondrial holoenzymes have additional functions and contribute to a dynamic network of elaborate interactions and cellular processes. Herein, we summarize the role of each RNase P subunit with a focus on the human nuclear RNP and its putative role in flawless gene expression in light of recent structural studies.

mPPI: a database extension to visualize structural interactome in a one-to-many manner

Protein-protein interaction (PPI) databases with structural information are useful to investigate biological functions at both systematic and atomic levels. However, most existing PPI databases only curate binary interactome. From the perspective of the display and function of PPI, as well as the structural binding interface, the related database and resources are summarized. We developed a database extension, named mPPI, for PPI structural visualization. Comparing with the existing structural interactomes that curate resolved PPI conformation in pairs, mPPI can visualize target protein and its multiple interactors simultaneously, which facilitates multi-target drug discovery and structure prediction of protein macro-complexes. By employing a protein-protein docking algorithm, mPPI largely extends the coverage of structural interactome from experimentally resolved complexes. mPPI is designed to be a customizable and convenient plugin for PPI databases. It possesses wide potential applications for various PPI databases, and it has been used for a neurodegenerative disease-related PPI database as demonstration. Scripts and implementation guidelines of mPPI are documented at the database tool website. Database URL  http://bis.zju.edu.cn/mppi/.