Human Protein Reference Database and Human Proteinpedia as resources for phosphoproteome analysis

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.

AMPK activation reduces cancer cell aggressiveness via inhibition of monoamine oxidase A (MAO-A) expression/activity

AIM

AMPK can be considered as an important target molecule for cancer for its unique ability to directly recognize cellular energy status. The main aim of this study is to explore the role of different AMPK activators in managing cancer cell aggressiveness and to understand the mechanistic details behind the process.

MAIN METHODS

First, we explored the AMPK expression pattern and its significance in different subtypes of lung cancer by accessing the TCGA data sets for LUNG, LUAD and LUSC patients and then established the correlation between AMPK expression pattern and overall survival of lung cancer patients using Kaplan-Meire plot. We further carried out several cell-based assays by employing different wet lab techniques including RT-PCR, Western Blot, proliferation, migration and invasion assays to fulfil the aim of the study.

KEY FINDINGS

SIGNIFICANCE: This study identifies the importance of AMPK activators as a repurposing agent for combating lung and colon cancer cell aggressiveness. It also suggests SRT-1720 as a potent repurposing agent for cancer treatment especially in NSCLC patients where a point mutation is present in LKB1.

Predicting gene disease associations with knowledge graph embeddings for diseases with curtailed information

Knowledge graph embeddings (KGE) are a powerful technique used in the biomedical domain to represent biological knowledge in a low dimensional space. However, a deep understanding of these methods is still missing, and, in particular, regarding their applications to prioritize genes associated with complex diseases with reduced genetic information. In this contribution, we built a knowledge graph (KG) by integrating heterogeneous biomedical data and generated KGE by implementing state-of-the-art methods, and two novel algorithms: Dlemb and BioKG2vec. Extensive testing of the embeddings with unsupervised clustering and supervised methods showed that KGE can be successfully implemented to predict genes associated with diseases and that our novel approaches outperform most existing algorithms in both scenarios. Our findings underscore the significance of data quality, preprocessing, and integration in achieving accurate predictions. Additionally, we applied KGE to predict genes linked to Intervertebral Disc Degeneration (IDD) and illustrated that functions pertinent to the disease are enriched within the prioritized gene set.

Tyr352 as a Predominant Phosphosite in the Understudied Kinase and Molecular Target, HIPK1: Implications for Cancer Therapy

Homeodomain-interacting protein kinase 1 (HIPK1) is majorly found in the nucleoplasm. HIPK1 is associated with cell proliferation, tumor necrosis factor-mediated cellular apoptosis, transcription regulation, and DNA damage response, and thought to play significant roles in health and common diseases such as cancer. Despite this, HIPK1 remains an understudied molecular target. In the present study, based on a systematic screening and mapping approach, we assembled 424 qualitative and 44 quantitative phosphoproteome datasets with 15 phosphosites in HIPK1 reported across multiple studies. These HIPK1 phosphosites were not currently attributed to any functions. Among them, Tyr352 within the kinase domain was identified as the predominant phosphosite modulated in 22 differential datasets. To analyze the functional association of HIPK1 Tyr352, we first employed a stringent criterion to derive its positively and negatively correlated protein phosphosites. Subsequently, we categorized the correlated phosphosites in known interactors, known/predicted kinases, and substrates of HIPK1, for their prioritized validation. Bioinformatics analysis identified their significant association with biological processes such as the regulation of RNA splicing, DNA-templated transcription, and cellular metabolic processes. HIPK1 Tyr352 was also identified to be upregulated in Her2+ cell lines and a subset of pancreatic and cholangiocarcinoma tissues. These data and the systems biology approach undertaken in the present study serve as a platform to explore the functional role of other phosphosites in HIPK1, and by extension, inform cancer drug discovery and oncotherapy innovation. In all, this study highlights the comprehensive phosphosite map of HIPK1 kinase and the first of its kind phosphosite-centric analysis of HIPK1 kinase based on global-level phosphoproteomics datasets derived from human cellular differential experiments across distinct experimental conditions.

Predicting Drug-Protein Interactions by Self-Adaptively Adjusting the Topological Structure of the Heterogeneous Network

Many powerful computational methods based on graph neural networks (GNNs) have been proposed to predict drug-protein interactions (DPIs). It can effectively reduce laboratory workload and the cost of drug discovery and drug repurposing. However, many clinical functions of drugs and proteins are unknown due to their unobserved indications. Therefore, it is difficult to establish a reliable drug-protein heterogeneous network that can describe the relationships between drugs and proteins based on the available information. To solve this problem, we propose a DPI prediction method that can self-adaptively adjust the topological structure of the heterogeneous networks, and name it SATS. SATS establishes a representation learning module based on graph attention network to carry out the drug-protein heterogeneous network. It can self-adaptively learn the relationships among the nodes based on their attributes and adjust the topological structure of the network according to the training loss of the model. Finally, SATS predicts the interaction propensity between drugs and proteins based on their embeddings. The experimental results show that SATS can effectively improve the topological structure of the network. The performance of SATS outperforms several state-of-the-art DPI prediction methods under various evaluation metrics. These prove that SATS is useful to deal with incomplete data and unreliable networks. The case studies on the top section of the prediction results further demonstrate that SATS is powerful for discovering novel DPIs.

Proteomic-Based Discovery of Predictive Biomarkers for Drug Therapy Response and Personalized Medicine in Chronic Immune Thrombocytopenia

Purpose

ITP is the most prevalent autoimmune blood disorder. The lack of predictive biomarkers for therapeutic response is a major challenge for physicians caring of chronic ITP patients. This study is aimed at identifying predictive biomarkers for drug therapy responses.

Methods

2D gel electrophoresis (2-DE) was performed to find differentially expressed proteins. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometer (MALDI-TOF MS) analysis was performed to identify protein spots. The Cytoscape software was employed to visualize and analyze the protein-protein interaction (PPI) network. Then, enzyme-linked immunosorbent assays (ELISA) were used to confirm the results of the proteins detected in the blood. The DAVID online software was used to explore the Gene Ontology and pathways involved in the disease.

Results

Three proteins, including APOA1, GC, and TF, were identified as hub-bottlenecks and confirmed by ELISA. Enrichment analysis results showed the importance of several biological processes and pathway, such as the PPAR signaling pathway, complement and coagulation cascades, platelet activation, vitamin digestion and absorption, fat digestion and absorption, cell adhesion molecule binding, and receptor binding. Conclusion and Clinical Relevance. Our results indicate that plasma proteins (APOA1, GC, and TF) can be suitable biomarkers for the prognosis of the response to drug therapy in ITP patients.

Neuroanatomical zones of human traumatic brain injury reveal significant differences in protein profile and protein oxidation: Implications for secondary injury events

Traumatic brain injury (TBI) causes significant neurological deficits and long-term degenerative changes. Primary injury in TBI entails distinct neuroanatomical zones, i.e., contusion (Ct) and pericontusion (PC). Their dynamic expansion could contribute to unpredictable neurological deterioration in patients. Molecular characterization of these zones compared with away from contusion (AC) zone is invaluable for TBI management. Using proteomics-based approach, we were able to distinguish Ct, PC and AC zones in human TBI brains. Ct was associated with structural changes (blood-brain barrier (BBB) disruption, neuroinflammation, axonal injury, demyelination and ferroptosis), while PC was associated with initial events of secondary injury (glutamate excitotoxicity, glial activation, accumulation of cytoskeleton proteins, oxidative stress, endocytosis) and AC displayed mitochondrial dysfunction that could contribute to secondary injury events and trigger long-term degenerative changes. Phosphoproteome analysis in these zones revealed that certain differentially phosphorylated proteins synergistically contribute to the injury events along with the differentially expressed proteins. Non-synaptic mitochondria (ns-mito) was associated with relatively more differentially expressed proteins (DEPs) compared to synaptosomes (Syn), while the latter displayed increased protein oxidation including tryptophan (Trp) oxidation. Proteomic analysis of immunocaptured complex I (CI) from Syn revealed increased Trp oxidation in Ct > PC > AC (vs. control). Oxidized W272 in the ND1 subunit of CI, revealed local conformational changes in ND1 and the neighboring subunits, as indicated by molecular dynamics simulation (MDS). Taken together, neuroanatomical zones in TBI show distinct protein profile and protein oxidation representing different primary and secondary injury events with potential implications for TBI pathology and neurological status of the patients.

Open a new epoch of arsenic trioxide investigation: ATOdb

Arsenic trioxide (ATO) is a great discovery in the treatment of acute promyelocytic leukemia (APL), which has been used in an increasing number of malignant diseases. Systematic integrative analysis will help to precisely understand the mechanism of ATO and find new combined drugs. Therefore, we developed a one-stop comprehensive database of ATO named ATOdb by manually compiling a wealth of experimentally supported ATO-related data from 3479 articles, and integrated analysis tools. The current version of ATOdb contains 8373 associations among 2300 ATO targets, 80 conditions and 262 combined drugs. Each entry in ATOdb contains detailed information on ATO targets, therapeutic/side effects, systems, cell names, cell types, regulations, detection methods, brief descriptions, references, etc. Furthermore, ATOdb also provides data visualization and analysis results such as the drug similarities, protein-protein interactions, and miRNA-mRNA relationships. An easy-to-use web interface was deployed in ATOdb for users to easily browse, search and download the data. In conclusion, ATOdb will serve as a valuable resource for in-depth study of the mechanism of ATO, discovery of new drug combination strategies, promotion of rational drug use and individualized treatments. ATOdb is freely available at http://bio-bigdata.hrbmu.edu.cn/ATOdb/index.jsp.

Geometric graph neural networks on multi-omics data to predict cancer survival outcomes

The advance of sequencing technologies has enabled a thorough molecular characterization of the genome in human cancers. To improve patient prognosis predictions and subsequent treatment strategies, it is imperative to develop advanced computational methods to analyze large-scale, high-dimensional genomic data. However, traditional machine learning methods face a challenge in handling the high-dimensional, low-sample size problem that is shown in most genomic data sets. To address this, our group has developed geometric network analysis techniques on multi-omics data in connection with prior biological knowledge derived from protein-protein interactions (PPIs) or pathways. Geometric features obtained from the genomic network, such as Ollivier-Ricci curvature and the invariant measure of the associated Markov chain, have been shown to be predictive of survival outcomes in various cancers. In this study, we propose a novel supervised deep learning method called geometric graph neural network (GGNN) that incorporates such geometric features into deep learning for enhanced predictive power and interpretability. More specifically, we utilize a state-of-the-art graph neural network with sparse connections between the hidden layers based on known biology of the PPI network and pathway information. Geometric features along with multi-omics data are then incorporated into the corresponding layers. The proposed approach utilizes a local-global principle in such a manner that highly predictive features are selected at the front layers and fed directly to the last layer for multivariable Cox proportional-hazards regression modeling. The method was applied to multi-omics data from the CoMMpass study of multiple myeloma and ten major cancers in The Cancer Genome Atlas (TCGA). In most experiments, our method showed superior predictive performance compared to other alternative methods.

Splicing Modulators Are Involved in Human Polyglutamine Diversification via Protein Complexes Shuttling between Nucleus and Cytoplasm

Length polymorphisms of polyglutamine (polyQs) in triplet-repeat-disease-causing genes have diversified during primate evolution despite them conferring a risk of human-specific diseases. To explain the evolutionary process of this diversification, there is a need to focus on mechanisms by which rapid evolutionary changes can occur, such as alternative splicing. Proteins that can bind polyQs are known to act as splicing factors and may provide clues about the rapid evolutionary process. PolyQs are also characterized by the formation of intrinsically disordered (ID) regions, so I hypothesized that polyQs are involved in the transportation of various molecules between the nucleus and cytoplasm to regulate mechanisms characteristic of humans such as neural development. To determine target molecules for empirical research to understand the evolutionary change, I explored protein-protein interactions (PPIs) involving the relevant proteins. This study identified pathways related to polyQ binding as hub proteins scattered across various regulatory systems, including regulation via PQBP1, VCP, or CREBBP. Nine ID hub proteins with both nuclear and cytoplasmic localization were found. Functional annotations suggested that ID proteins containing polyQs are involved in regulating transcription and ubiquitination by flexibly changing PPI formation. These findings explain the relationships among splicing complex, polyQ length variations, and modifications in neural development.

Multi-omics analysis reveals distinct non-reversion mechanisms of PARPi resistance in BRCA1- versus BRCA2-deficient mammary tumors

BRCA1 and BRCA2 both function in DNA double-strand break repair by homologous recombination (HR). Due to their HR defect, BRCA1/2-deficient cancers are sensitive to poly(ADP-ribose) polymerase inhibitors (PARPis), but they eventually acquire resistance. Preclinical studies yielded several PARPi resistance mechanisms that do not involve BRCA1/2 reactivation, but their relevance in the clinic remains elusive. To investigate which BRCA1/2-independent mechanisms drive spontaneous resistance in vivo, we combine molecular profiling with functional analysis of HR of matched PARPi-naive and PARPi-resistant mouse mammary tumors harboring large intragenic deletions that prevent reactivation of BRCA1/2. We observe restoration of HR in 62% of PARPi-resistant BRCA1-deficient tumors but none in the PARPi-resistant BRCA2-deficient tumors. Moreover, we find that 53BP1 loss is the prevalent resistance mechanism in HR-proficient BRCA1-deficient tumors, whereas resistance in BRCA2-deficient tumors is mainly induced by PARG loss. Furthermore, combined multi-omics analysis identifies additional genes and pathways potentially involved in modulating PARPi response.

Graph embedding and Gaussian mixture variational autoencoder network for end-to-end analysis of single-cell RNA sequencing data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary technology to determine the precise gene expression of individual cells and identify cell heterogeneity and subpopulations. However, technical limitations of scRNA-seq lead to heterogeneous and sparse data. Here, we present autoCell, a deep-learning approach for scRNA-seq dropout imputation and feature extraction. autoCell is a variational autoencoding network that combines graph embedding and a probabilistic depth Gaussian mixture model to infer the distribution of high-dimensional, sparse scRNA-seq data. We validate autoCell on simulated datasets and biologically relevant scRNA-seq. We show that interpolation of autoCell improves the performance of existing tools in identifying cell developmental trajectories of human preimplantation embryos. We identify disease-associated astrocytes (DAAs) and reconstruct DAA-specific molecular networks and ligand-receptor interactions involved in cell-cell communications using Alzheimer's disease as a prototypical example. autoCell provides a toolbox for end-to-end analysis of scRNA-seq data, including visualization, clustering, imputation, and disease-specific gene network identification.

Systematic analysis to identify novel disease indications and plausible potential chemical leads of glutamate ionotropic receptor NMDA type subunit 1, GRIN1

Schizophrenia is a mental illness affecting the normal lifestyle of adults and early adolescents incurring major symptoms as jumbled speech, involvement in everyday activities eventually got reduced, patients always struggle with attention and memory, reason being both the genetic and environmental factors responsible for altered brain chemistry and structure, resulting in schizophrenia and associated orphan diseases. The network biology describes the interactions among genes/proteins encoding molecular mechanisms of biological processes, development, and diseases. Besides, all the molecular networks, protein-protein Interaction Networks have been significant in distinguishing the pathogenesis of diseases and thereby drug discovery. The present meta-analysis prioritizes novel disease indications viz. rare and orphan diseases associated with target Glutamate Ionotropic Receptor NMDA Type Subunit 1, GRIN1 using text mining knowledge-based tools. Furthermore, ZINC database was virtually screened, and binding conformation of selected compounds was performed and resulted in the identification of Narciclasine (ZINC04097652) and Alvespimycin (ZINC73138787) as potential inhibitors. Furthermore, docked complexes were subjected to MD simulation studies which suggests that the identified leads could be a better potential drug to recuperate schizophrenia.

Comprehensive molecular characterization of hypertension-related genes in cancer

Background

During cancer treatment, patients have a significantly higher risk of developing cardiovascular complications such as hypertension. In this study, we investigated the internal relationships between hypertension and different types of cancer.

Methods

First, we comprehensively characterized the involvement of 10 hypertension-related genes across 33 types of cancer. The somatic copy number alteration (CNA) and single nucleotide variant (SNV) of each gene were identified for each type of cancer. Then, the expression patterns of hypertension-related genes were analyzed across 14 types of cancer. The hypertension-related genes were aberrantly expressed in different types of cancer, and some were associated with the overall survival of patients or the cancer stage. Subsequently, the interactions between hypertension-related genes and clinically actionable genes (CAGs) were identified by analyzing the co-expressions and protein–protein interactions.

Results

We found that certain hypertension-related genes were correlated with CAGs. Next, the pathways associated with hypertension-related genes were identified. The positively correlated pathways included epithelial to mesenchymal transition, hormone androgen receptor, and receptor tyrosine kinase, and the negatively correlated pathways included apoptosis, cell cycle, and DNA damage response. Finally, the correlations between hypertension-related genes and drug sensitivity were evaluated for different drugs and different types of cancer. The hypertension-related genes were all positively or negatively correlated with the resistance of cancer to the majority of anti-cancer drugs. These results highlight the importance of hypertension-related genes in cancer.

Conclusions

This study provides an approach to characterize the relationship between hypertension-related genes and cancers in the post-genomic era.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40959-022-00136-z.

Interpretable deep learning translation of GWAS and multi-omics findings to identify pathobiology and drug repurposing in Alzheimer's disease

Translating human genetic findings (genome-wide association studies [GWAS]) to pathobiology and therapeutic discovery remains a major challenge for Alzheimer's disease (AD). We present a network topology-based deep learning framework to identify disease-associated genes (NETTAG). We leverage non-coding GWAS loci effects on quantitative trait loci, enhancers and CpG islands, promoter regions, open chromatin, and promoter flanking regions under the protein-protein interactome. Via NETTAG, we identified 156 AD-risk genes enriched in druggable targets. Combining network-based prediction and retrospective case-control observations with 10 million individuals, we identified that usage of four drugs (ibuprofen, gemfibrozil, cholecalciferol, and ceftriaxone) is associated with reduced likelihood of AD incidence. Gemfibrozil (an approved lipid regulator) is significantly associated with 43% reduced risk of AD compared with simvastatin using an active-comparator design (95% confidence interval 0.51-0.63, p < 0.0001). In summary, NETTAG offers a deep learning methodology that utilizes GWAS and multi-genomic findings to identify pathobiology and drug repurposing in AD.

Evolutionary Divergence of Phosphorylation to Regulate Interactive Protein Networks in Lower and Higher Species

The phosphorylation of proteins affects their functions in extensively documented circumstances. However, the role of phosphorylation in many interactive networks of proteins remains very elusive due to the experimental limits of exploring the transient interaction in a large complex of assembled proteins induced by stimulation. Previous studies have suggested that phosphorylation is a recent evolutionary process that differently regulates ortholog proteins in numerous lineages of living organisms to create new functions. Despite the fact that numerous phospho-proteins have been compared between species, little is known about the organization of the full phospho-proteome, the role of phosphorylation to orchestrate large interactive networks of proteins, and the intertwined phospho-landscape in these networks. In this report, we aimed to investigate the acquired role of phosphate addition in the phenomenon of protein networking in different orders of living organisms. Our data highlighted the acquired status of phosphorylation in organizing large, connected assemblages in Homo sapiens. The protein networking guided by phosphorylation turned out to be prominent in humans, chaotic in yeast, and weak in flies. Furthermore, the molecular functions of GO annotation enrichment regulated by phosphorylation were found to be drastically different between flies, yeast, and humans, suggesting an evolutionary drift specific to each species.

DGHNE: network enhancement-based method in identifying disease-causing genes through a heterogeneous biomedical network

The identification of disease-causing genes is critical for mechanistic understanding of disease etiology and clinical manipulation in disease prevention and treatment. Yet the existing approaches in tackling this question are inadequate in accuracy and efficiency, demanding computational methods with higher identification power. Here, we proposed a new method called DGHNE to identify disease-causing genes through a heterogeneous biomedical network empowered by network enhancement. First, a disease-disease association network was constructed by the cosine similarity scores between phenotype annotation vectors of diseases, and a new heterogeneous biomedical network was constructed by using disease-gene associations to connect the disease-disease network and gene-gene network. Then, the heterogeneous biomedical network was further enhanced by using network embedding based on the Gaussian random projection. Finally, network propagation was used to identify candidate genes in the enhanced network. We applied DGHNE together with five other methods into the most updated disease-gene association database termed DisGeNet. Compared with all other methods, DGHNE displayed the highest area under the receiver operating characteristic curve and the precision-recall curve, as well as the highest precision and recall, in both the global 5-fold cross-validation and predicting new disease-gene associations. We further performed DGHNE in identifying the candidate causal genes of Parkinson's disease and diabetes mellitus, and the genes connecting hyperglycemia and diabetes mellitus. In all cases, the predicted causing genes were enriched in disease-associated gene ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways, and the gene-disease associations were highly evidenced by independent experimental studies.

A systematic study of HIF1A cofactors in hypoxic cancer cells

Hypoxia inducible factor 1 alpha (HIF1A) is a transcription factor (TF) that forms highly structural and functional protein-protein interactions with other TFs to promote gene expression in hypoxic cancer cells. However, despite the importance of these TF-TF interactions, we still lack a comprehensive view of many of the TF cofactors involved and how they cooperate. In this study, we systematically studied HIF1A cofactors in eight cancer cell lines using the computational motif mining tool, SIOMICS, and discovered 201 potential HIF1A cofactors, which included 21 of the 29 known HIF1A cofactors in public databases. These 201 cofactors were statistically and biologically significant, with 19 of the top 37 cofactors in our study directly validated in the literature. The remaining 18 were novel cofactors. These discovered cofactors can be essential to HIF1A's regulatory functions and may lead to the discovery of new therapeutic targets in cancer treatment.

AnimalTFDB 4.0: a comprehensive animal transcription factor database updated with variation and expression annotations

Transcription factors (TFs) are proteins that interact with specific DNA sequences to regulate gene expression and play crucial roles in all kinds of biological processes. To keep up with new data and provide a more comprehensive resource for TF research, we updated the Animal Transcription Factor Database (AnimalTFDB) to version 4.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB4/) with up-to-date data and functions. We refined the TF family rules and prediction pipeline to predict TFs in genome-wide protein sequences from Ensembl. As a result, we predicted 274 633 TF genes and 150 726 transcription cofactor genes in AnimalTFDB 4.0 in 183 animal genomes, which are 86 more species than AnimalTFDB 3.0. Besides double data volume, we also added the following new annotations and functions to the database: (i) variations (including mutations) on TF genes in various human cancers and other diseases; (ii) predicted post-translational modification sites (including phosphorylation, acetylation, methylation and ubiquitination sites) on TFs in 8 species; (iii) TF regulation in autophagy; (iv) comprehensive TF expression annotation for 38 species; (v) exact and batch search functions allow users to search AnimalTFDB flexibly. AnimalTFDB 4.0 is a useful resource for studying TF and transcription regulation, which contains comprehensive annotation and classification of TFs and transcription cofactors.

Overview of methods for characterization and visualization of a protein–protein interaction network in a multi-omics integration context

At the heart of the cellular machinery through the regulation of cellular functions, protein–protein interactions (PPIs) have a significant role. PPIs can be analyzed with network approaches. Construction of a PPI network requires prediction of the interactions. All PPIs form a network. Different biases such as lack of data, recurrence of information, and false interactions make the network unstable. Integrated strategies allow solving these different challenges. These approaches have shown encouraging results for the understanding of molecular mechanisms, drug action mechanisms, and identification of target genes. In order to give more importance to an interaction, it is evaluated by different confidence scores. These scores allow the filtration of the network and thus facilitate the representation of the network, essential steps to the identification and understanding of molecular mechanisms. In this review, we will discuss the main computational methods for predicting PPI, including ones confirming an interaction as well as the integration of PPIs into a network, and we will discuss visualization of these complex data.

MLysPRED: graph-based multi-view clustering and multi-dimensional normal distribution resampling techniques to predict multiple lysine sites

Posttranslational modification of lysine residues, K-PTM, is one of the most popular PTMs. Some lysine residues in proteins can be continuously or cascaded covalently modified, such as acetylation, crotonylation, methylation and succinylation modification. The covalent modification of lysine residues may have some special functions in basic research and drug development. Although many computational methods have been developed to predict lysine PTMs, up to now, the K-PTM prediction methods have been modeled and learned a single class of K-PTM modification. In view of this, this study aims to fill this gap by building a multi-label computational model that can be directly used to predict multiple K-PTMs in proteins. In this study, a multi-label prediction model, MLysPRED, is proposed to identify multiple lysine sites using features generated from human protein sequences. In MLysPRED, three kinds of multi-label sequence encoding algorithms (MLDBPB, MLPSDAAP, MLPSTAAP) are proposed and combined with three encoding strategies (CHHAA, DR and Kmer) to convert preprocessed lysine sequences into effective numerical features. A multidimensional normal distribution oversampling technique and graph-based multi-view clustering under-sampling algorithm were first proposed and incorporated to reduce the proportion of the original training samples, and multi-label nearest neighbor algorithm is used for classification. It is observed that MLysPRED achieved an Aiming of 92.21%, Coverage of 94.98%, Accuracy of 89.63%, Absolute-True of 81.46% and Absolute-False of 0.0682 on the independent datasets. Additionally, comparison of results with five existing predictors also indicated that MLysPRED is very promising and encouraging to predict multiple K-PTMs in proteins. For the convenience of the experimental scientists, 'MLysPRED' has been deployed as a user-friendly web-server at http://47.100.136.41:8181.

Single-cell network biology characterizes cell type gene regulation for drug repurposing and phenotype prediction in Alzheimer’s disease

Dysregulation of gene expression in Alzheimer’s disease (AD) remains elusive, especially at the cell type level. Gene regulatory network, a key molecular mechanism linking transcription factors (TFs) and regulatory elements to govern gene expression, can change across cell types in the human brain and thus serve as a model for studying gene dysregulation in AD. However, AD-induced regulatory changes across brain cell types remains uncharted. To address this, we integrated single-cell multi-omics datasets to predict the gene regulatory networks of four major cell types, excitatory and inhibitory neurons, microglia and oligodendrocytes, in control and AD brains. Importantly, we analyzed and compared the structural and topological features of networks across cell types and examined changes in AD. Our analysis shows that hub TFs are largely common across cell types and AD-related changes are relatively more prominent in some cell types (e.g., microglia). The regulatory logics of enriched network motifs (e.g., feed-forward loops) further uncover cell type-specific TF-TF cooperativities in gene regulation. The cell type networks are also highly modular and several network modules with cell-type-specific expression changes in AD pathology are enriched with AD-risk genes. The further disease-module-drug association analysis suggests cell-type candidate drugs and their potential target genes. Finally, our network-based machine learning analysis systematically prioritized cell type risk genes likely involved in AD. Our strategy is validated using an independent dataset which showed that top ranked genes can predict clinical phenotypes (e.g., cognitive impairment) of AD with reasonable accuracy. Overall, this single-cell network biology analysis provides a comprehensive map linking genes, regulatory networks, cell types and drug targets and reveals cell-type gene dysregulation in AD.

Alzheimer’s Disease (AD) is the leading cause of dementia. It affects parts of the brain that control language, behavior, and memory. The human brain is comprised of tens of billions of cells, such as neuronal cells that transmit information via electrical and chemical signals, and glial cells that maintain the brain’s immune system. Researchers have found that AD causes changes in the expression of genes within the brain cells. Gene expression is a tightly regulated process involving interconnected networks of multiple genes. Understanding how these gene networks change in AD is critical to identifying genetic biomarkers and potential drug targets. Using genomic data of post-mortem brains diagnosed with AD and healthy individuals, we identified gene networks that play a crucial role in regulating biological processes within neuronal and glial cells. We utilized these gene networks to make predictions on existing FDA approved drugs that could potentially be repurposed for AD. Furthermore, we used a machine learning strategy to identify novel genes that are more likely to be involved in AD pathology. The systems-level approach lends itself to analysis of single-cell genomics data of other human diseases.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

Protein-protein interaction networks as miners of biological discovery

Protein-protein interactions (PPIs) form the basis of a myriad of biological pathways and mechanism, such as the formation of protein-complexes or the components of signaling cascades. Here, we reviewed experimental methods for identifying PPI pairs, including yeast two-hybrid, mass spectrometry, co-localization, and co-immunoprecipitation. Furthermore, a range of computational methods leveraging biochemical properties, evolution history, protein structures and more have enabled identification of additional PPIs. Given the wealth of known PPIs, we reviewed important network methods to construct and analyze networks of PPIs. These methods aid biological discovery through identifying hub genes and dynamic changes in the network, and have been thoroughly applied in various fields of biological research. Lastly, we discussed the challenges and future direction of research utilizing the power of PPI networks. This article is protected by copyright. All rights reserved.

Decoding multilevel relationships with the human tissue-cell-molecule network

Understanding the biological functions of molecules in specific human tissues or cell types is crucial for gaining insights into human physiology and disease. To address this issue, it is essential to systematically uncover associations among multilevel elements consisting of disease phenotypes, tissues, cell types and molecules, which could pose a challenge because of their heterogeneity and incompleteness. To address this challenge, we describe a new methodological framework, called Graph Local InfoMax (GLIM), based on a human multilevel network (HMLN) that we established by introducing multiple tissues and cell types on top of molecular networks. GLIM can systematically mine the potential relationships between multilevel elements by embedding the features of the HMLN through contrastive learning. Our simulation results demonstrated that GLIM consistently outperforms other state-of-the-art algorithms in disease gene prediction. Moreover, GLIM was also successfully used to infer cell markers and rewire intercellular and molecular interactions in the context of specific tissues or diseases. As a typical case, the tissue-cell-molecule network underlying gastritis and gastric cancer was first uncovered by GLIM, providing systematic insights into the mechanism underlying the occurrence and development of gastric cancer. Overall, our constructed methodological framework has the potential to systematically uncover complex disease mechanisms and mine high-quality relationships among phenotypical, tissue, cellular and molecular elements.

Antifungal activity of an artificial peptide aptamer SNP-D4 against Fusarium oxysporum

Fusarium oxysporum f. sp. cubense (FOC4) is a pathogen of banana fusarium wilt, which is a serious problem that has plagued the tropical banana industry for many years. The pathogenic mechanism is complex and unclear, so the prevention and control in agricultural production applications is ineffective. SNP-D4, an artificial peptide aptamer, was identified and specifically inhibited FOC4. To evaluate the efficacy of SNP-D4, FoC4 spores were treated with purified SNP-D4 to calculate the germination and fungicide rates. Damage of FOC4 spores was observed by staining with propidium iodide (PI). Eight proteins of FOC4 were identified to have high affinity for SNP-D4 by a pull-down method combined with Q-Exactive mass spectrometry. Of these eight proteins, A0A5C6SPC6, the aldehyde dehydrogenase of FOC4, was selected as an example to scrutinize the interaction sites with SNP-D4. Molecular docking revealed that Thr66 on the peptide loop of SNP-D4 bound with Tyr437 near the catalytic center of A0A5C6SPC6. Subsequently 42 spore proteins which exhibited associations with the eight proteins were retrieved for protein-protein interaction analysis, demonstrating that SNP-D4 interfered with pathways including 'translation', 'folding, sorting and degradation', 'transcription', 'signal transduction' and 'cell growth and death', eventually causing the inhibition of growth of FOC4. This study not only investigated the possible pathogenic mechanism of FOC4, but also provided a potential antifungal agent SNP-D4 for use in the control of banana wilt disease.

Computational analyses of mechanism of action (MoA): data, methods and integration

The elucidation of a compound's Mechanism of Action (MoA) is a challenging task in the drug discovery process, but it is important in order to rationalise phenotypic findings and to anticipate potential side-effects. Bioinformatic approaches, advances in machine learning techniques and the increasing deposition of high-throughput data in public databases have significantly contributed to recent advances in the field, but it is not straightforward to decide which data and methods are most suitable to use in a given case. In this review, we focus on these methods and data and their applications in generating MoA hypotheses for subsequent experimental validation. We discuss compound-specific data such as -omics, cell morphology and bioactivity data, as well as commonly used supplementary prior knowledge such as network and pathway data, and provide information on databases where this data can be accessed. In terms of methodologies, we discuss both well-established methods (connectivity mapping, pathway enrichment) as well as more developing methods (neural networks and multi-omics integration). Finally, we review case studies where the MoA of a compound was successfully suggested from computational analysis by incorporating multiple data modalities and/or methodologies. Our aim for this review is to provide researchers with insights into the benefits and drawbacks of both the data and methods in terms of level of understanding, biases and interpretation - and to highlight future avenues of investigation which we foresee will improve the field of MoA elucidation, including greater public access to -omics data and methodologies which are capable of data integration.

Biomedical data, computational methods and tools for evaluating disease-disease associations

In recent decades, exploring potential relationships between diseases has been an active research field. With the rapid accumulation of disease-related biomedical data, a lot of computational methods and tools/platforms have been developed to reveal intrinsic relationship between diseases, which can provide useful insights to the study of complex diseases, e.g. understanding molecular mechanisms of diseases and discovering new treatment of diseases. Human complex diseases involve both external phenotypic abnormalities and complex internal molecular mechanisms in organisms. Computational methods with different types of biomedical data from phenotype to genotype can evaluate disease-disease associations at different levels, providing a comprehensive perspective for understanding diseases. In this review, available biomedical data and databases for evaluating disease-disease associations are first summarized. Then, existing computational methods for disease-disease associations are reviewed and classified into five groups in terms of the usages of biomedical data, including disease semantic-based, phenotype-based, function-based, representation learning-based and text mining-based methods. Further, we summarize software tools/platforms for computation and analysis of disease-disease associations. Finally, we give a discussion and summary on the research of disease-disease associations. This review provides a systematic overview for current disease association research, which could promote the development and applications of computational methods and tools/platforms for disease-disease associations.

Shared Molecular Mechanisms between Atherosclerosis and Periodontitis by Analyzing the Transcriptomic Alterations of Peripheral Blood Monocytes

OBJECTIVE

This study investigated the nature of shared transcriptomic alterations in PBMs from periodontitis and atherosclerosis to unravel molecular mechanisms underpinning their association.

METHODS

Gene expression data from PBMs from patients with periodontitis and those with atherosclerosis were each downloaded from the GEO database. Differentially expressed genes (DEGs) in periodontitis and atherosclerosis were identified through differential gene expression analysis. The disease-related known genes related to periodontitis and atherosclerosis each were downloaded from the DisGeNET database. A Venn diagram was constructed to identify crosstalk genes from four categories: DEGs expressed in periodontitis, periodontitis-related known genes, DEGs expressed in atherosclerosis, and atherosclerosis-related known genes. A weighted gene coexpression network analysis (WGCNA) was performed to identify significant coexpression modules, and then, coexpressed gene interaction networks belonging to each significant module were constructed to identify the core crosstalk genes.

RESULTS

Functional enrichment analysis of significant modules obtained by WGCNA analysis showed that several pathways might play the critical crosstalk role in linking both diseases, including bacterial invasion of epithelial cells, platelet activation, and Mitogen-Activated Protein Kinases (MAPK) signaling. By constructing the gene interaction network of significant modules, the core crosstalk genes in each module were identified and included: for GSE23746 dataset, RASGRP2 in the blue module and VAMP7 and SNX3 in the green module, as well as HMGB1 and SUMO1 in the turquoise module were identified; for GSE61490 dataset, SEC61G, PSMB2, SELPLG, and FIBP in the turquoise module were identified.

CONCLUSION

Exploration of available transcriptomic datasets revealed core crosstalk genes (RASGRP2, VAMP7, SNX3, HMGB1, SUMO1, SEC61G, PSMB2, SELPLG, and FIBP) and significant pathways (bacterial invasion of epithelial cells, platelet activation, and MAPK signaling) as top candidate molecular linkage mechanisms between atherosclerosis and periodontitis.

HFIP: an integrated multi-omics data and knowledge platform for the precision medicine of heart failure

As the terminal clinical phenotype of almost all types of cardiovascular diseases, heart failure (HF) is a complex and heterogeneous syndrome leading to considerable morbidity and mortality. Existing HF-related omics studies mainly focus on case/control comparisons, small cohorts of special subtypes, etc., and a large amount of multi-omics data and knowledge have been generated. However, it is difficult for researchers to obtain biological and clinical insights from these scattered data and knowledge. In this paper, we built the Heart Failure Integrated Platform (HFIP) for data exploration, fusion analysis and visualization by collecting and curating existing multi-omics data and knowledge from various public sources and also provided an auto-updating mechanism for future integration. The developed HFIP contained 253 datasets (7842 samples), multiple analysis flow, and 14 independent tools. In addition, based on the integration of existing databases and literature, a knowledge base for HF was constructed with a scoring system for evaluating the relationship between molecular signals and HF. The knowledge base includes 1956 genes and annotation information. The literature mining module was developed to assist the researcher to overview the hotspots and contexts in basic and clinical research. HFIP can be used as a data-driven and knowledge-guided platform for the basic and clinical research of HF.

Database URL: http://heartfailure.medical-bigdata.com

CPLM 4.0: an updated database with rich annotations for protein lysine modifications

Here, we reported the compendium of protein lysine modifications (CPLM 4.0, http://cplm.biocuckoo.cn/), a data resource for various post-translational modifications (PTMs) specifically occurred at the side-chain amino group of lysine residues in proteins. From the literature and public databases, we collected 450 378 protein lysine modification (PLM) events, and combined them with the existing data of our previously developed protein lysine modification database (PLMD 3.0). In total, CPLM 4.0 contained 592 606 experimentally identified modification events on 463 156 unique lysine residues of 105 673 proteins for up to 29 types of PLMs across 219 species. Furthermore, we carefully annotated the data using the knowledge from 102 additional resources that covered 13 aspects, including variation and mutation, disease-associated information, protein-protein interaction, protein functional annotation, DNA & RNA element, protein structure, chemical-target relation, mRNA expression, protein expression/proteomics, subcellular localization, biological pathway annotation, functional domain annotation, and physicochemical property. Compared to PLMD 3.0 and other existing resources, CPLM 4.0 achieved a >2-fold increase in collection of PLM events, with a data volume of ∼45GB. We anticipate that CPLM 4.0 can serve as a more useful database for further study of PLMs.

Network-based analysis of virulence factors for uncovering Aeromonas veronii pathogenesis

BACKGROUND

Aeromonas veronii is a bacterial pathogen in aquaculture, which produces virulence factors to enable it colonize and evade host immune defense. Given that experimental verification of virulence factors is time-consuming and laborious, few virulence factors have been characterized. Moreover, most studies have only focused on single virulence factors, resulting in biased interpretation of the pathogenesis of A. veronii.

RESULTS

In this study, a PPI network at genome-wide scale for A. veronii was first constructed followed by prediction and mapping of virulence factors on the network. When topological characteristics were analyzed, the virulence factors had higher degree and betweenness centrality than other proteins in the network. In particular, the virulence factors tended to interact with each other and were enriched in two network modules. One of the modules mainly consisted of histidine kinases, response regulators, diguanylate cyclases and phosphodiesterases, which play important roles in two-component regulatory systems and the synthesis and degradation of cyclic-diGMP. Construction of the interspecies PPI network between A. veronii and its host Oreochromis niloticus revealed that the virulence factors interacted with homologous proteins in the host. Finally, the structures and interacting sites of the virulence factors during interaction with host proteins were predicted.

CONCLUSIONS

The findings here indicate that the virulence factors probably regulate the virulence of A. veronii by involving in signal transduction pathway and manipulate host biological processes by mimicking and binding competitively to host proteins. Our results give more insight into the pathogenesis of A. veronii and provides important information for designing targeted antibacterial drugs.

Measuring pathway database coverage of the phosphoproteome

Protein phosphorylation is one of the best known post-translational mechanisms playing a key role in the regulation of cellular processes. Over 100,000 distinct phosphorylation sites have been discovered through constant improvement of mass spectrometry based phosphoproteomics in the last decade. However, data saturation is occurring and the bottleneck of assigning biologically relevant functionality to phosphosites needs to be addressed. There has been finite success in using data-driven approaches to reveal phosphosite functionality due to a range of limitations. The alternate, more suitable approach is making use of prior knowledge from literature-derived databases. Here, we analysed seven widely used databases to shed light on their suitability to provide functional insights into phosphoproteomics data. We first determined the global coverage of each database at both the protein and phosphosite level. We also determined how consistent each database was in its phosphorylation annotations compared to a global standard. Finally, we looked in detail at the coverage of each database over six experimental datasets. Our analysis highlights the relative strengths and weaknesses of each database, providing a guide in how each can be best used to identify biological mechanisms in phosphoproteomic data.

Elucidation of the Mechanisms and Molecular Targets of Qishen Yiqi Formula for the Treatment of Pulmonary Arterial Hypertension using a Bioinformatics/Network Topology-based Strategy

BACKGROUND AND OBJECTIVE

Qishen Yiqi formula (QSYQ) is used to treat cardiovascular disease in the clinical practice of traditional Chinese medicine. However, few studies have explored whether QSYQ affects pulmonary arterial hypertension (PAH), and the mechanisms of action and molecular targets of QSYQ for the treatment of PAH are unclear. A bioinformatics/network topology-based strategy was used to identify the bioactive ingredients, putative targets, and molecular mechanisms of QSYQ in PAH.

METHODS

A network pharmacology-based strategy was employed by integrating active component gathering, target prediction, PAH gene collection, network topology, and gene enrichment analysis to systematically explore the multicomponent synergistic mechanisms.

RESULTS

In total, 107 bioactive ingredients of QSYQ and 228 ingredient targets were identified. Moreover, 234 PAH-related differentially expressed genes with a |fold change| >2 and an adjusted P value < 0.005 were identified between the PAH patient and control groups, and 266 therapeutic targets were identified. The pathway enrichment analysis indicated that 85 pathways, including the PI3K-Akt, MAPK, and HIF-1 signaling pathways, were significantly enriched. TP53 was the core target gene, and 7 other top genes (MAPK1, RELA, NFKB1, CDKN1A, AKT1, MYC, and MDM2) were the key genes in the gene-pathway network based on the effects of QSYQ on PAH.

CONCLUSION

An integrative investigation based on network pharmacology may elucidate the multicomponent synergistic mechanisms of QSYQ in PAH and lay a foundation for further animal experiments, human clinical trials and rational clinical applications of QSYQ.

Identification of the Key Role of NF-κB Signaling Pathway in the Treatment of Osteoarthritis With Bushen Zhuangjin Decoction, a Verification Based on Network Pharmacology Approach

This study aimed to identify whether the NF-κB signaling pathway plays a key role in the treatment of osteoarthritis (OA) with Bushen Zhuangjin Decoction (BZD) based on a typical network pharmacology approach (NPA). Four sequential experiments were performed: 1) conventional high-performance liquid chromatography (HPLC), 2) preliminary observation of the therapeutic effects of BZD, 3) NPA using three OA-related gene expression profiles, and 4) verification of the key pathway identified by NPA. Only one HPLC-verified compound (paeoniflorin) was identified from the candidate compounds discovered by NPA. The genes verified in the preliminary observation were also identified by NPA. NPA identified a key role for the NF-κB signaling pathway in the treatment of OA with BZD, which was confirmed by conventional western blot analysis. This study identified and verified NF-κB signaling pathway as the most important inflammatory signaling pathway involved in the mechanisms of BZD for treating OA by comparing the NPA results with conventional methods. Our findings also indicate that NPA is a powerful tool for exploring the molecular targets of complex herbal formulations, such as BZD.

Exploring the Imbalance of Periodontitis Immune System From the Cellular to Molecular Level

Periodontitis is a common chronic inflammatory disease of periodontal tissue, mostly concentrated in people over 30 years old. Statistics show that compared with foreign countries, the prevalence of periodontitis in China is as high as 40%, and the prevalence of periodontal disease is more than 90%, which must arouse our great attention. Diagnosis and treatment of periodontitis currently rely mainly on clinical criteria, and the exploration of the etiologic criteria is relatively lacking. We, therefore, have explored the pathogenesis of periodontitis from the perspective of immune imbalance. By predicting the fraction of 22 immune cells in periodontitis tissues and comparing them with normal tissues, we found that multiple immune cell infiltration in periodontitis tissues was inhibited and this feature can clearly distinguish periodontitis from normal tissues. Further, protein interaction network (PPI) and transcription regulation network have been constructed based on differentially expressed genes (DEGs) to explore the interaction function modules and regulation pathways. Three functional modules have been revealed and top TFs such as EGR1 and ETS1 have been shown to regulate the expression of periodontitis-related immune genes that play an important role in the formation of the immunosuppressive microenvironment. The classifier was also used to verify the reliability of periodontitis features obtained at the cellular and molecular levels. In conclusion, we have revealed the immune microenvironment and molecular characteristics of periodontitis, which will help to better understand the mechanism of periodontitis and its application in clinical diagnosis and treatment.

Proteome analysis using machine learning approaches and its applications to diseases

With the tremendous developments in the fields of biological and medical technologies, huge amounts of data are generated in the form of genomic data, images in medical databases or as data on protein sequences, and so on. Analyzing this data through different tools sheds light on the particulars of the disease and our body's reactions to it, thus, aiding our understanding of the human health. Most useful of these tools is artificial intelligence and deep learning (DL). The artificially created neural networks in DL algorithms help extract viable data from the datasets, and further, to recognize patters in these complex datasets. Therefore, as a part of machine learning, DL helps us face all the various challenges that come forth during protein prediction, protein identification and their quantification. Proteomics is the study of such proteins, their structures, features, properties and so on. As a form of data science, Proteomics has helped us progress excellently in the field of genomics technologies. One of the major techniques used in proteomics studies is mass spectrometry (MS). However, MS is efficient with analysis of large datasets only with the added help of informatics approaches for data analysis and interpretation; these mainly include machine learning and deep learning algorithms. In this chapter, we will discuss in detail the applications of deep learning and various algorithms of machine learning in proteomics.

Structural similarity-based prediction of host factors associated with SARS-CoV-2 infection and pathogenesis

The current pandemic resulted from SARS-CoV-2 still remains as the major public health concern globally. The precise mechanism of viral pathogenesis is not fully understood, which remains a major hurdle for medical intervention. Here we generated an interactome profile of protein-protein interactions based on host and viral protein structural similarities information. Further computational biological study combined with Gene enrichment analysis predicted key enriched pathways associated with viral pathogenesis. The results show that axon guidance, membrane trafficking, vesicle-mediated transport, apoptosis, clathrin-mediated endocytosis, Vpu mediated degradation of CD4 T cell, and interferon-gamma signaling are key events associated in SARS-CoV-2 life cycle. Further, degree centrality analysis reveals that IRF1/9/7, TP53, and CASP3, UBA52, and UBC are vital proteins for IFN-γ-mediated signaling, apoptosis, and proteasomal degradation of CD4, respectively. We crafted chronological events of the virus life cycle. The SARS-CoV-2 enters through clathrin-mediated endocytosis, and the genome is trafficked to the early endosomes in a RAB5-dependent manner. It is predicted to replicate in a double-membrane vesicle (DMV) composed of the endoplasmic reticulum, autophagosome, and ERAD machinery. The SARS-CoV-2 down-regulates host translational machinery by interacting with protein kinase R, PKR-like endoplasmic reticulum kinase, and heme-regulated inhibitor and can phosphorylate eIF2a. The virion assembly occurs in the ER-Golgi intermediate compartment (ERGIC) organized by the spike and matrix protein. Collectively, we have established a spatial link between viral entry, RNA synthesis, assembly, pathogenesis, and their associated diverse host factors, those could pave the way for therapeutic intervention.

Multivariate genome wide association and network analysis of subcortical imaging phenotypes in Alzheimer's disease

BACKGROUND

Genome-wide association studies (GWAS) have identified many individual genes associated with brain imaging quantitative traits (QTs) in Alzheimer's disease (AD). However single marker level association discovery may not be able to address the underlying biological interactions with disease mechanism.

RESULTS

In this paper, we used the MGAS (Multivariate Gene-based Association test by extended Simes procedure) tool to perform multivariate GWAS on eight AD-relevant subcortical imaging measures. We conducted multiple iPINBPA (integrative Protein-Interaction-Network-Based Pathway Analysis) network analyses on MGAS findings using protein-protein interaction (PPI) data, and identified five Consensus Modules (CMs) from the PPI network. Functional annotation and network analysis were performed on the identified CMs. The MGAS yielded significant hits within APOE, TOMM40 and APOC1 genes, which were known AD risk factors, as well as a few new genes such as LAMA1, XYLB, HSD17B7P2, and NPEPL1. The identified five CMs were enriched by biological processes related to disorders such as Alzheimer's disease, Legionellosis, Pertussis, and Serotonergic synapse.

CONCLUSIONS

The statistical power of coupling MGAS with iPINBPA was higher than traditional GWAS method, and yielded new findings that were missed by GWAS. This study provides novel insights into the molecular mechanism of Alzheimer's Disease and will be of value to novel gene discovery and functional genomic studies.

Association Between DCC Polymorphisms and Susceptibility to Autism Spectrum Disorder

Autism spectrum disorder (ASD) represents a group of childhood-onset lifelong neuro-developmental disorders. However, the association between single nucleotide polymorphisms (SNPs) in the deleted in colorectal carcinoma (DCC) gene and ASD susceptibility remains unclear. We investigated the association between ASD susceptibility and seven SNPs in DCC on the basis of a case-control study (231 ASD cases and 242 controls) in Chinese Han. We found that there was no association between ASD susceptibility and the seven SNPs in DCC; however, T-A haplotype (rs2229082-rs2270954), T-A-T-C haplotype (rs2229082-rs2270954-rs2292043-rs2292044), C-G-T-C-T haplotype (rs934345-rs17753970-rs2229082-rs2270954-rs2292043), C-G-T-C-T-G haplotype (rs934345-rs17753970-rs2229082-rs2270954-rs2292043-rs2292044), and G-G-T-C-C-C-C haplotype (rs934345-rs17753970-rs2229082-rs2270954-rs2292043-rs2292044-rs16956878) were associated with ASD susceptibility. Our results indicate that the haplotypes formed on the basis of the seven SNPs in DCC may be implicated in ASD.

Is the Proteome of Bronchoalveolar Lavage Extracellular Vesicles a Marker of Advanced Lung Cancer?

Simple Summary

Bronchoalveolar lavage is routinely collected during bronchoscopy for cytology analysis in the diagnostic of lung cancer. Due to low sensitivity of this method, early-stage cancers are undetected, lowering the treatment success. In this study, we analyzed extracellular vesicles isolated from bronchoalveolar lavage of lung cancer suspects by mass spectrometry-based proteomics. The protein composition of bronchoalveolar lavage extracellular vesicles of late-stage cancer showed a higher proteome complexity associated with mortality within the two year follow-up period. We identified a potential therapeutic target DNMT3B complex which was significantly expressed in bronchoalveolar lavage extracellular vesicles as well as in tumor tissue. Bronchoalveolar lavage extracellular vesicles proteome analysis of immune markers indicates the presence of markers of innate immune and fibroblast cells.

Abstract

Acellular bronchoalveolar lavage (BAL) proteomics can partially separate lung cancer from non-lung cancer patients based on principal component analysis and multivariate analysis. Furthermore, the variance in the proteomics data sets is correlated mainly with lung cancer status and, to a lesser extent, smoking status and gender. Despite these advances BAL small and large extracellular vehicles (EVs) proteomes reveal aberrant protein expression in paracrine signaling mechanisms in cancer initiation and progression. We consequently present a case-control study of 24 bronchoalveolar lavage extracellular vesicle samples which were analyzed by state-of-the-art liquid chromatography-mass spectrometry (LC-MS). We obtained evidence that BAL EVs proteome complexity correlated with lung cancer stage 4 and mortality within two years´ follow-up (p value = 0.006). The potential therapeutic target DNMT3B complex is significantly up-regulated in tumor tissue and BAL EVs. The computational analysis of the immune and fibroblast cell markers in EVs suggests that patients who deceased within the follow-up period display higher marker expression indicative of innate immune and fibroblast cells (four out of five cases). This study provides insights into the proteome content of BAL EVs and their correlation to clinical outcomes.

In silico and in vitro studies reveal complement system drives coagulation cascade in SARS-CoV-2 pathogenesis

The emergence and continued spread of SARS-CoV-2 have resulted in a public health emergency across the globe. The lack of knowledge on the precise mechanism of viral pathogenesis is impeding medical intervention. In this study, we have taken both in silico and in vitro experimental approaches to unravel the mechanism of viral pathogenesis associated with complement and coagulation pathways. Based on the structural similarities of viral and host proteins, we initially generated a protein-protein interactome profile. Further computational analysis combined with Gene Ontology (GO) analysis and KEGG pathway analysis predicted key annotated pathways associated with viral pathogenesis. These include MAPK signaling, complement, and coagulation cascades, endocytosis, PD-L1 expression, PD-1 checkpoint pathway in cancer and C-type lectin receptor signaling pathways. Degree centrality analysis pinned down to MAPK1, MAPK3, AKT1, and SRC are crucial drivers of signaling pathways and often overlap with the associated pathways. Most strikingly, the complement and coagulation cascade and platelet activation pathways are interconnected, presumably directing thrombotic activity observed in severe or critical cases of COVID-19. This is complemented by in vitro studies of Huh7 cell infection and analysis of the transcriptome and proteomic profile of gene candidates during viral infection. The most known candidates associated with complement and coagulation cascade signaling by KEGG pathway analysis showed significant up-regulated fold change during viral infection. Collectively both in silico and in vitro studies suggest complement and coagulation cascade signaling are a mechanism for intravascular coagulation, thrombotic changes, and associated complications in severe COVID-19 patients.

Evolution of Sequence-based Bioinformatics Tools for Protein-protein Interaction Prediction

Protein-protein interactions (PPIs) are the physical connections between two or more proteins via electrostatic forces or hydrophobic effects. Identification of the PPIs is pivotal, which contributes to many biological processes including protein function, disease incidence, and therapy design. The experimental identification of PPIs via high-throughput technology is time-consuming and expensive. Bioinformatics approaches are expected to solve such restrictions. In this review, our main goal is to provide an inclusive view of the existing sequence-based computational prediction of PPIs. Initially, we briefly introduce the currently available PPI databases and then review the state-of-the-art bioinformatics approaches, working principles, and their performances. Finally, we discuss the caveats and future perspective of the next generation algorithms for the prediction of PPIs.

Advancing from protein interactomes and gene co-expression networks towards multi-omics-based composite networks: approaches for predicting and extracting biological knowledge

Prediction of biological interaction networks from single-omics data has been extensively implemented to understand various aspects of biological systems. However, more recently, there is a growing interest in integrating multi-omics datasets for the prediction of interactomes that provide a global view of biological systems with higher descriptive capability, as compared to single omics. In this review, we have discussed various computational approaches implemented to infer and analyze two of the most important and well studied interactomes: protein-protein interaction networks and gene co-expression networks. We have explicitly focused on recent methods and pipelines implemented to infer and extract biologically important information from these interactomes, starting from utilizing single-omics data and then progressing towards multi-omics data. Accordingly, recent examples and case studies are also briefly discussed. Overall, this review will provide a proper understanding of the latest developments in protein and gene network modelling and will also help in extracting practical knowledge from them.

GPS-PBS: A Deep Learning Framework to Predict Phosphorylation Sites that Specifically Interact with Phosphoprotein-Binding Domains

Protein phosphorylation is essential for regulating cellular activities by modifying substrates at specific residues, which frequently interact with proteins containing phosphoprotein-binding domains (PPBDs) to propagate the phosphorylation signaling into downstream pathways. Although massive phosphorylation sites (p-sites) have been reported, most of their interacting PPBDs are unknown. Here, we collected 4458 known PPBD-specific binding p-sites (PBSs), considerably improved our previously developed group-based prediction system (GPS) algorithm, and implemented a deep learning plus transfer learning strategy for model training. Then, we developed a new online service named GPS-PBS, which can hierarchically predict PBSs of 122 single PPBD clusters belonging to two groups and 16 families. By comparison, GPS-PBS achieved a highly competitive accuracy against other existing tools. Using GPS-PBS, we predicted 371,018 mammalian p-sites that potentially interact with at least one PPBD, and revealed that various PPBD-containing proteins (PPCPs) and protein kinases (PKs) can simultaneously regulate the same p-sites to orchestrate important pathways, such as the PI3K-Akt signaling pathway. Taken together, we anticipate GPS-PBS can be a great help for further dissecting phosphorylation signaling networks.

dbPTM in 2019: exploring disease association and cross-talk of post-translational modifications

The dbPTM (http://dbPTM.mbc.nctu.edu.tw/) has been maintained for over 10 years with the aim to provide functional and structural analyses for post-translational modifications (PTMs). In this update, dbPTM not only integrates more experimentally validated PTMs from available databases and through manual curation of literature but also provides PTM-disease associations based on non-synonymous single nucleotide polymorphisms (nsSNPs). The high-throughput deep sequencing technology has led to a surge in the data generated through analysis of association between SNPs and diseases, both in terms of growth amount and scope. This update thus integrated disease-associated nsSNPs from dbSNP based on genome-wide association studies. The PTM substrate sites located at a specified distance in terms of the amino acids encoded from nsSNPs were deemed to have an association with the involved diseases. In recent years, increasing evidence for crosstalk between PTMs has been reported. Although mass spectrometry-based proteomics has substantially improved our knowledge about substrate site specificity of single PTMs, the fact that the crosstalk of combinatorial PTMs may act in concert with the regulation of protein function and activity is neglected. Because of the relatively limited information about concurrent frequency and functional relevance of PTM crosstalk, in this update, the PTM sites neighboring other PTM sites in a specified window length were subjected to motif discovery and functional enrichment analysis. This update highlights the current challenges in PTM crosstalk investigation and breaks the bottleneck of how proteomics may contribute to understanding PTM codes, revealing the next level of data complexity and proteomic limitation in prospective PTM research.

Mass-spectrometry-based draft of the Arabidopsis proteome

Plants are essential for life and are extremely diverse organisms with unique molecular capabilities¹. Here we present a quantitative atlas of the transcriptomes, proteomes and phosphoproteomes of 30 tissues of the model plant Arabidopsis thaliana. Our analysis provides initial answers to how many genes exist as proteins (more than 18,000), where they are expressed, in which approximate quantities (a dynamic range of more than six orders of magnitude) and to what extent they are phosphorylated (over 43,000 sites). We present examples of how the data may be used, such as to discover proteins that are translated from short open-reading frames, to uncover sequence motifs that are involved in the regulation of protein production, and to identify tissue-specific protein complexes or phosphorylation-mediated signalling events. Interactive access to this resource for the plant community is provided by the ProteomicsDB and ATHENA databases, which include powerful bioinformatics tools to explore and characterize Arabidopsis proteins, their modifications and interactions.