Deep learning based prediction of reversible HAT/HDAC-specific lysine acetylation

Deep learning prioritizes cancer mutations that alter protein nucleocytoplasmic shuttling to drive tumorigenesis

Genetic variants can affect protein function by driving aberrant subcellular localization. However, comprehensive analysis of how mutations promote tumor progression by influencing nuclear localization is currently lacking. Here, we systematically characterize potential shuttling-attacking mutations (SAMs) across cancers through developing the deep learning model pSAM for the ab initio decoding of the sequence determinants of nucleocytoplasmic shuttling. Leveraging cancer mutations across 11 cancer types, we find that SAMs enrich functional genetic variations and critical genes in cancer. We experimentally validate a dozen SAMs, among which R14M in PTEN, P255L in CHFR, etc. are identified to disrupt the nuclear localization signals through interfering their interactions with importins. Further studies confirm that the nucleocytoplasmic shuttling altered by SAMs in PTEN and CHFR rewire the downstream signaling and eliminate their function of tumor suppression. Thus, this study will help to understand the molecular traits of nucleocytoplasmic shuttling and their dysfunctions mediated by genetic variants.

MlyPredCSED: based on extreme point deviation compensated clustering combined with cross-scale convolutional neural networks to predict multiple lysine sites in human

In post-translational modification, covalent bonds on lysine and attached chemical groups significantly change proteins' physical and chemical properties. They shape protein structures, enhance function and stability, and are vital for physiological processes, affecting health and disease through mechanisms like gene expression, signal transduction, protein degradation, and cell metabolism. Although lysine (K) modification sites are considered among the most common types of post-translational modifications in proteins, research on K-PTMs has largely overlooked the synergistic effects between different modifications and lacked the techniques to address the problem of sample imbalance. Based on this, the Extreme Point Deviation Compensated Clustering (EPDCC) Undersampling algorithm was proposed in this study and combined with Cross-Scale Convolutional Neural Networks (CSCNNs) to develop a novel computational tool, MlyPredCSED, for simultaneously predicting multiple lysine modification sites. MlyPredCSED employs Multi-Label Position-Specific Triad Amino Acid Propensity and the physicochemical properties of amino acids to enhance the richness of sequence information. To address the challenge of sample imbalance, the innovative EPDCC Undersampling technique was introduced to adjust the majority class samples. The model's training and testing phase relies on the advanced CSCNN framework. MlyPredCSED, through cross-validation and testing, outperformed existing models, especially in complex categories with multiple modification sites. This research not only provides an efficient method for the identification of lysine modification sites but also demonstrates its value in biological research and drug development. To facilitate efficient use of MlyPredCSED by researchers, we have specifically developed an accessible free web tool: http://www.mlypredcsed.com.

In silico identification of Histone Deacetylase inhibitors using Streamlined Masked Transformer-based Pretrained features

Histone Deacetylases (HDACs) are enzymes that regulate gene expression by removing acetyl groups from histones. They are involved in various diseases, including neurodegenerative, cardiovascular, inflammatory, and metabolic disorders, as well as fibrosis in the liver, lungs, and kidneys. Successfully identifying potent HDAC inhibitors may offer a promising approach to treating these diseases. In addition to experimental techniques, researchers have introduced several in silico methods for identifying HDAC inhibitors. However, these existing computer-aided methods have shortcomings in their modeling stages, which limit their applications. In our study, we present a Streamlined Masked Transformer-based Pretrained (SMTP) encoder, which can be used to generate features for downstream tasks. The training process of the SMTP encoder was directed by masked attention-based learning, enhancing the model's generalizability in encoding molecules. The SMTP features were used to develop 11 classification models identifying 11 HDAC isoforms. We trained SMTP, a lightweight encoder, with only 1.9 million molecules, a smaller number than other known molecular encoders, yet its discriminant ability remains competitive. The results revealed that machine learning models developed using the SMTP feature set outperformed those developed using other feature sets in 8 out of 11 classification tasks. Additionally, chemical diversity analysis confirmed the encoder's effectiveness in distinguishing between two classes of molecules.

Artificial Intelligence Transforming Post-Translational Modification Research

Post-Translational Modifications (PTMs) are covalent changes to amino acids that occur after protein synthesis, including covalent modifications on side chains and peptide backbones. Many PTMs profoundly impact cellular and molecular functions and structures, and their significance extends to evolutionary studies as well. In light of these implications, we have explored how artificial intelligence (AI) can be utilized in researching PTMs. Initially, rationales for adopting AI and its advantages in understanding the functions of PTMs are discussed. Then, various deep learning architectures and programs, including recent applications of language models, for predicting PTM sites on proteins and the regulatory functions of these PTMs are compared. Finally, our high-throughput PTM-data-generation pipeline, which formats data suitably for AI training and predictions is described. We hope this review illuminates areas where future AI models on PTMs can be improved, thereby contributing to the field of PTM bioengineering.

Current computational tools for protein lysine acylation site prediction

As a main subtype of post-translational modification (PTM), protein lysine acylations (PLAs) play crucial roles in regulating diverse functions of proteins. With recent advancements in proteomics technology, the identification of PTM is becoming a data-rich field. A large amount of experimentally verified data is urgently required to be translated into valuable biological insights. With computational approaches, PLA can be accurately detected across the whole proteome, even for organisms with small-scale datasets. Herein, a comprehensive summary of 166 in silico PLA prediction methods is presented, including a single type of PLA site and multiple types of PLA sites. This recapitulation covers important aspects that are critical for the development of a robust predictor, including data collection and preparation, sample selection, feature representation, classification algorithm design, model evaluation, and method availability. Notably, we discuss the application of protein language models and transfer learning to solve the small-sample learning issue. We also highlight the prediction methods developed for functionally relevant PLA sites and species/substrate/cell-type-specific PLA sites. In conclusion, this systematic review could potentially facilitate the development of novel PLA predictors and offer useful insights to researchers from various disciplines.

Ketogenic diet reshapes cancer metabolism through lysine β-hydroxybutyrylation

Lysine β-hydroxybutyrylation (Kbhb) is a post-translational modification induced by the ketogenic diet (KD), a diet showing therapeutic effects on multiple human diseases. Little is known how cellular processes are regulated by Kbhb. Here we show that protein Kbhb is strongly affected by the KD through a multi-omics analysis of mouse livers. Using a small training dataset with known functions, we developed a bioinformatics method for the prediction of functionally important lysine modification sites (pFunK), which revealed functionally relevant Kbhb sites on various proteins, including aldolase B (ALDOB) Lys108. KD consumption or β-hydroxybutyrate supplementation in hepatocellular carcinoma cells increases ALDOB Lys108bhb and inhibits the enzymatic activity of ALDOB. A Kbhb-mimicking mutation (p.Lys108Gln) attenuates ALDOB activity and its binding to substrate fructose-1,6-bisphosphate, inhibits mammalian target of rapamycin signalling and glycolysis, and markedly suppresses cancer cell proliferation. Our study reveals a critical role of Kbhb in regulating cancer cell metabolism and provides a generally applicable algorithm for predicting functionally important lysine modification sites.

TransPTM: a transformer-based model for non-histone acetylation site prediction

Protein acetylation is one of the extensively studied post-translational modifications (PTMs) due to its significant roles across a myriad of biological processes. Although many computational tools for acetylation site identification have been developed, there is a lack of benchmark dataset and bespoke predictors for non-histone acetylation site prediction. To address these problems, we have contributed to both dataset creation and predictor benchmark in this study. First, we construct a non-histone acetylation site benchmark dataset, namely NHAC, which includes 11 subsets according to the sequence length ranging from 11 to 61 amino acids. There are totally 886 positive samples and 4707 negative samples for each sequence length. Secondly, we propose TransPTM, a transformer-based neural network model for non-histone acetylation site predication. During the data representation phase, per-residue contextualized embeddings are extracted using ProtT5 (an existing pre-trained protein language model). This is followed by the implementation of a graph neural network framework, which consists of three TransformerConv layers for feature extraction and a multilayer perceptron module for classification. The benchmark results reflect that TransPTM has the competitive performance for non-histone acetylation site prediction over three state-of-the-art tools. It improves our comprehension on the PTM mechanism and provides a theoretical basis for developing drug targets for diseases. Moreover, the created PTM datasets fills the gap in non-histone acetylation site datasets and is beneficial to the related communities. The related source code and data utilized by TransPTM are accessible at https://www.github.com/TransPTM/TransPTM.

Emerging trends in post-translational modification: Shedding light on Glioblastoma multiforme

Recent multi-omics studies, including proteomics, transcriptomics, genomics, and metabolomics have revealed the critical role of post-translational modifications (PTMs) in the progression and pathogenesis of Glioblastoma multiforme (GBM). Further, PTMs alter the oncogenic signaling events and offer a novel avenue in GBM therapeutics research through PTM enzymes as potential biomarkers for drug targeting. In addition, PTMs are critical regulators of chromatin architecture, gene expression, and tumor microenvironment (TME), that play a crucial function in tumorigenesis. Moreover, the implementation of artificial intelligence and machine learning algorithms enhances GBM therapeutics research through the identification of novel PTM enzymes and residues. Herein, we briefly explain the mechanism of protein modifications in GBM etiology, and in altering the biologics of GBM cells through chromatin remodeling, modulation of the TME, and signaling pathways. In addition, we highlighted the importance of PTM enzymes as therapeutic biomarkers and the role of artificial intelligence and machine learning in protein PTM prediction.

Identification and characterization of putative biomarkers and therapeutic axis in Glioblastoma multiforme microenvironment

Non-cellular secretory components, including chemokines, cytokines, and growth factors in the tumor microenvironment, are often dysregulated, impacting tumorigenesis in Glioblastoma multiforme (GBM) microenvironment, where the prognostic significance of the current treatment remains unsatisfactory. Recent studies have demonstrated the potential of post-translational modifications (PTM) and their respective enzymes, such as acetylation and ubiquitination in GBM etiology through modulating signaling events. However, the relationship between non-cellular secretory components and post-translational modifications will create a research void in GBM therapeutics. Therefore, we aim to bridge the gap between non-cellular secretory components and PTM modifications through machine learning and computational biology approaches. Herein, we highlighted the importance of BMP1, CTSB, LOX, LOXL1, PLOD1, MMP9, SERPINE1, and SERPING1 in GBM etiology. Further, we demonstrated the positive relationship between the E2 conjugating enzymes (Ube2E1, Ube2H, Ube2J2, Ube2C, Ube2J2, and Ube2S), E3 ligases (VHL and GNB2L1) and substrate (HIF1A). Additionally, we reported the novel HAT1-induced acetylation sites of Ube2S (K211) and Ube2H (K8, K52). Structural and functional characterization of Ube2S (8) and Ube2H (1) have identified their association with protein kinases. Lastly, our results found a putative therapeutic axis HAT1-Ube2S(K211)-GNB2L1-HIF1A and potential predictive biomarkers (CTSB, HAT1, Ube2H, VHL, and GNB2L1) that play a critical role in GBM pathogenesis.

qPTM: an updated database for PTM dynamics in human, mouse, rat and yeast

Post-translational modifications (PTMs) are critical molecular mechanisms that regulate protein functions temporally and spatially in various organisms. Since most PTMs are dynamically regulated, quantifying PTM events under different states is crucial for understanding biological processes and diseases. With the rapid development of high-throughput proteomics technologies, massive quantitative PTM proteome datasets have been generated. Thus, a comprehensive one-stop data resource for surfing big data will benefit the community. Here, we updated our previous phosphorylation dynamics database qPhos to the qPTM (http://qptm.omicsbio.info). In qPTM, 11 482 553 quantification events among six types of PTMs, including phosphorylation, acetylation, glycosylation, methylation, SUMOylation and ubiquitylation in four different organisms were collected and integrated, and the matched proteome datasets were included if available. The raw mass spectrometry based false discovery rate control and the recurrences of identifications among datasets were integrated into a scoring system to assess the reliability of the PTM sites. Browse and search functions were improved to facilitate users in swiftly and accurately acquiring specific information. The results page was revised with more abundant annotations, and time-course dynamics data were visualized in trend lines. We expected the qPTM database to be a much more powerful and comprehensive data repository for the PTM research community.

Comprehensive Machine Learning Prediction of Extensive Enzymatic Reactions

New enzyme functions exist within the increasing number of unannotated protein sequences. Novel enzyme discovery is necessary to expand the pathways that can be accessed by metabolic engineering for the biosynthesis of functional compounds. Accordingly, various machine learning models have been developed to predict enzymatic reactions. However, the ability to predict unknown reactions that are not included in the training data has not been clarified. In order to cover uncertain and unknown reactions, a wider range of reaction types must be demonstrated by the models. Here, we establish 16 expanded enzymatic reaction prediction models developed using various machine learning algorithms, including deep neural network. Improvements in prediction performances over that of our previous study indicate that the updated methods are more effective for the prediction of enzymatic reactions. Overall, the deep neural network model trained with combined substrate-enzyme-product information exhibits the highest prediction accuracy with Macro F₁ scores up to 0.966 and with robust prediction of unknown enzymatic reactions that are not included in the training data. This model can predict more extensive enzymatic reactions in comparison to previously reported models. This study will facilitate the discovery of new enzymes for the production of useful substances.

PredPromoter-MF(2L): A Novel Approach of Promoter Prediction Based on Multi-source Feature Fusion and Deep Forest

Promoters short DNA sequences play vital roles in initiating gene transcription. However, it remains a challenge to identify promoters using conventional experiment techniques in a high-throughput manner. To this end, several computational predictors based on machine learning models have been developed, while their performance is unsatisfactory. In this study, we proposed a novel two-layer predictor, called PredPromoter-MF(2L), based on multi-source feature fusion and ensemble learning. PredPromoter-MF(2L) was developed based on various deep features learned by a pre-trained deep learning network model and sequence-derived features. Feature selection based on XGBoost was applied to reduce fused features dimensions, and a cascade deep forest model was trained on the selected feature subset for promoter prediction. The results both fivefold cross-validation and independent test demonstrated that PredPromoter-MF(2L) outperformed state-of-the-art methods.

Deacetylation of K481 and K484 on Penaeid Shrimp Hemocyanin Is Critical for Antibacterial Activity

Although invertebrates' innate immunity relies on several immune-like molecules, the diversity of these molecules and their immune response mechanisms are not well understood. Here, we show that Penaeus vannamei hemocyanin (PvHMC) undergoes specific deacetylation under Vibrio parahaemolyticus and LPS challenge. In vitro deacetylation of PvHMC increases its binding capacity with LPS and antibacterial activity against Gram-negative bacteria. Lysine residues K481 and K484 on the Ig-like domain of PvHMC are the main acetylation sites modulated by the acetyltransferase TIP60 and deacetylase HDAC3. Deacetylation of PvHMC on K481 and K484 allows PvHMC to form a positively charged binding pocket that interacts directly with LPS, whereas acetylation abrogates the positive charge to decrease PvHMC-LPS attraction. Besides, V. parahaemolyticus and LPS challenge increases the expression of Pvhdac3 to induce PvHMC deacetylation. This work indicates that, during bacterial infections, deacetylation of hemocyanin is crucial for binding with LPS to clear Gram-negative bacteria in crustaceans.

Mini-review: Recent advances in post-translational modification site prediction based on deep learning

Post-translational modifications (PTMs) are closely linked to numerous diseases, playing a significant role in regulating protein structures, activities, and functions. Therefore, the identification of PTMs is crucial for understanding the mechanisms of cell biology and diseases therapy. Compared to traditional machine learning methods, the deep learning approaches for PTM prediction provide accurate and rapid screening, guiding the downstream wet experiments to leverage the screened information for focused studies. In this paper, we reviewed the recent works in deep learning to identify phosphorylation, acetylation, ubiquitination, and other PTM types. In addition, we summarized PTM databases and discussed future directions with critical insights.

Accurate Models of Substrate Preferences of Post-Translational Modification Enzymes from a Combination of mRNA Display and Deep Learning

Promiscuous post-translational modification (PTM) enzymes often display nonobvious substrate preferences by acting on diverse yet well-defined sets of peptides and/or proteins. Understanding of substrate fitness landscapes for PTM enzymes is important in many areas of contemporary science, including natural product biosynthesis, molecular biology, and biotechnology. Here, we report an integrated platform for accurate profiling of substrate preferences for PTM enzymes. The platform features (i) a combination of mRNA display with next-generation sequencing as an ultrahigh throughput technique for data acquisition and (ii) deep learning for data analysis. The high accuracy (>0.99 in each of two studies) of the resulting deep learning models enables comprehensive analysis of enzymatic substrate preferences. The models can quantify fitness across sequence space, map modification sites, and identify important amino acids in the substrate. To benchmark the platform, we performed profiling of a Ser dehydratase (LazBF) and a Cys/Ser cyclodehydratase (LazDEF), two enzymes from the lactazole biosynthesis pathway. In both studies, our results point to complex enzymatic preferences, which, particularly for LazBF, cannot be reduced to a set of simple rules. The ability of the constructed models to dissect such complexity suggests that the developed platform can facilitate a wider study of PTM enzymes.

RNA modification patterns based on major RNA modifications define tumor microenvironment characteristics in glioblastoma

RNA modifications play a major role in tumorigenicity and progression, but the expression and function in glioblastoma (GBM) have not been well described. In this study, we developed a GBM score based on the differentially expressed genes (DEGs) between groups showing RNA modification patterns. We assessed the association between the GBM score and tumor microenvironment (TME) characteristics. Based on the gene expression of these regulators, we identified two clusters with distinct RNA modification patterns. Kaplan–Meier survival curves showed that patients in cluster 1 had worse survival than those in cluster 2. Kaplan–Meier and multivariate Cox regression analyses showed that GBM scores (based on DEGs between RNA modification patterns) are an independent predictive biomarker for patient prognosis. Besides, we found that samples with high scores were significantly associated with epithelial-to-mesenchymal transition and immune checkpoints, while samples with low scores were associated with cell cycle regulation. Importantly, GBM-score markedly positively correlated drug resistance, while negatively correlated with drug sensitive. The responders of anti-PD-1/PD-L1 immunotherapy tend to have a lower GBM score than non-responders. In conclusion, our comprehensive analysis of multiple RNA modifications in GBM revealed that RNA modification regulators were closely correlated with TME.

Deep Learning-Based Advances In Protein Posttranslational Modification Site and Protein Cleavage Prediction

Posttranslational modification (PTM ) is a ubiquitous phenomenon in both eukaryotes and prokaryotes which gives rise to enormous proteomic diversity. PTM mostly comes in two flavors: covalent modification to polypeptide chain and proteolytic cleavage. Understanding and characterization of PTM is a fundamental step toward understanding the underpinning of biology. Recent advances in experimental approaches, mainly mass-spectrometry-based approaches, have immensely helped in obtaining and characterizing PTMs. However, experimental approaches are not enough to understand and characterize more than 450 different types of PTMs and complementary computational approaches are becoming popular. Recently, due to the various advancements in the field of Deep Learning (DL), along with the explosion of applications of DL to various fields, the field of computational prediction of PTM has also witnessed the development of a plethora of deep learning (DL)-based approaches. In this book chapter, we first review some recent DL-based approaches in the field of PTM site prediction. In addition, we also review the recent advances in the not-so-studied PTM , that is, proteolytic cleavage predictions. We describe advances in PTM prediction by highlighting the Deep learning architecture, feature encoding, novelty of the approaches, and availability of the tools/approaches. Finally, we provide an outlook and possible future research directions for DL-based approaches for PTM prediction.

RNA adenosine modifications related to prognosis and immune infiltration in osteosarcoma

Background

RNA adenosine modifications, which are primarily mediated by “writer” enzymes (RMWs), play a key role in epigenetic regulation in various biological processes, including tumorigenesis. However, the expression and prognostic role of these genes in osteosarcoma (OS) remain unclear.

Methods

Univariate and multivariate Cox analyses were used to construct the RMW signature for OS using Target datasets. RMW expression in OS tissue was detected by qPCR analysis. Xcell and GSVA were used to determine the relationship between RMWs and immune infiltration. The DGIdb and CMap databases were used for drug prediction. In vivo and in vitro experiments showed that strophanthidin elicited antitumor activity against OS.

Results

A 3-RMW (CSTF2, ADAR and WTAP) prognostic signature in OS was constructed using the Target dataset and verified using GEO datasets and 63 independent OS tissues via qPCR analysis. High-risk OS patients had poor overall survival, and the prognostic signature was an independent prognostic factor for OS. Functional studies showed that tumour-, metabolism-, cell cycle- and immune-related pathways were related to high risk. Next, we found that RMW-derived high-risk patients exhibited increased infiltration of M2 macrophages and cDCs. Furthermore, we predicted the potential drugs for OS using the DGIdb and CMap databases. In vivo and in vitro experiments showed that strophanthidin elicited antitumor activity against OS by repressing cell growth and inducing cell cycle arrest at the G1 phase.

Conclusion

The 3-RWM-based prognostic signature established in this study is a novel gene signature associated with immune infiltration, and strophanthidin was identified as a candidate therapy for OS by repressing OS cell growth and the cell cycle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03415-6.

Comprehensive characterization of m6A methylation and its impact on prognosis, genome instability, and tumor microenvironment in hepatocellular carcinoma

BACKGROUND

N6-methyladenosine (m6A) RNA regulation was recently reported to be important in carcinogenesis and cancer development. However, the characteristics of m6A modification and its correlations with clinical features, genome instability, tumor microenvironments (TMEs), and immunotherapy responses in hepatocellular carcinoma (HCC) have not been fully explored.

METHODS

We systematically analyzed the m6A regulator-based expression patterns of 486 patients with HCC from The Cancer Genome Atlas and Gene Expression Omnibus databases, and correlated these patterns with clinical outcomes, somatic mutations, TME cell infiltration, and immunotherapy responses. The m6A score was developed by principal component analysis to evaluate m6A modifications in individual patients.

RESULTS

M6A regulators were dysregulated in HCC samples, among which 18 m6A regulators were identified as risk factors for prognosis. Three m6A regulator-based expression patterns, namely m6A clusters, were determined among HCC patients by m6A regulators with different m6A scores, somatic mutation counts, and specific TME features. Additionally, three distinct m6A regulator-associated gene-based expression patterns were also identified based on prognosis-associated genes that were differentially expressed among the three m6A clusters, showing similar properties as the m6A regulator-based expression patterns. Higher m6A scores were correlated with older age, advanced stages, lower overall survival, higher somatic mutation counts, elevated PD-L1 expression levels, and poorer responses to immune checkpoint inhibitors. The m6A score was validated as an independent and valuable prognostic factor for HCC.

CONCLUSION

M6A modification is correlated with genome instability and TME in HCC. Evaluating m6A regulator-based expression patterns and the m6A score of individual tumors may help identify candidate biomarkers for prognosis prediction and immunotherapeutic strategy selection.

A multi-network integration approach for measuring disease similarity based on ncRNA regulation and heterogeneous information

Background

Measuring similarity between complex diseases has significant implications for revealing the pathogenesis of diseases and development in the domain of biomedicine. It has been consentaneous that functional associations between disease-related genes and semantic associations can be applied to calculate disease similarity. Currently, more and more studies have demonstrated the profound involvement of non-coding RNA in the regulation of genome organization and gene expression. Thus, taking ncRNA into account can be useful in measuring disease similarities. However, existing methods ignore the regulation functions of ncRNA in biological process. In this study, we proposed a novel deep-learning method to deduce disease similarity.

Results

In this article, we proposed a novel method, ImpAESim, a framework integrating multiple networks embedding to learn compact feature representations and disease similarity calculation. We first utilize three different disease-related information networks to build up a heterogeneous network, after a network diffusion process, RWR, a compact feature learning model composed of classic Auto Encoder (AE) and improved AE model is proposed to extract constraints and low-dimensional feature representations. We finally obtain an accurate and low-dimensional feature representation of diseases, then we employed the cosine distance as the measurement of disease similarity.

Conclusion

ImpAESim focuses on extracting a low-dimensional vector representation of features based on ncRNA regulation, and gene–gene interaction network. Our method can significantly reduce the calculation bias resulted from the sparse disease associations which are derived from semantic associations.

Comparative analysis of machine learning-based approaches for identifying therapeutic peptides targeting SARS-CoV-2

Coronavirus disease 2019 (COVID-19) has impacted public health as well as societal and economic well-being. In the last two decades, various prediction algorithms and tools have been developed for predicting antiviral peptides (AVPs). The current COVID-19 pandemic has underscored the need to develop more efficient and accurate machine learning (ML)-based prediction algorithms for the rapid identification of therapeutic peptides against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Several peptide-based ML approaches, including anti-coronavirus peptides (ACVPs), IL-6 inducing epitopes and other epitopes targeting SARS-CoV-2, have been implemented in COVID-19 therapeutics. Owing to the growing interest in the COVID-19 field, it is crucial to systematically compare the existing ML algorithms based on their performances. Accordingly, we comprehensively evaluated the state-of-the-art IL-6 and AVP predictors against coronaviruses in terms of core algorithms, feature encoding schemes, performance evaluation metrics and software usability. A comprehensive performance assessment was then conducted to evaluate the robustness and scalability of the existing predictors using well-constructed independent validation datasets. Additionally, we discussed the advantages and disadvantages of the existing methods, providing useful insights into the development of novel computational tools for characterizing and identifying epitopes or ACVPs. The insights gained from this review are anticipated to provide critical guidance to the scientific community in the rapid design and development of accurate and efficient next-generation in silico tools against SARS-CoV-2.

CREB1K292 and HINFPK330 as Putative Common Therapeutic Targets in Alzheimer's and Parkinson's Disease

Integration of omics data and deciphering the mechanism of a biological regulatory network could be a promising approach to reveal the molecular mechanism involved in the progression of complex diseases, including Alzheimer's and Parkinson's. Despite having an overlapping mechanism in the etiology of Alzheimer's disease (AD) and Parkinson's disease (PD), the exact mechanism and signaling molecules behind them are still unknown. Further, the acetylation mechanism and histone deacetylase (HDAC) enzymes provide a positive direction toward studying the shared phenomenon between AD and PD pathogenesis. For instance, increased expression of HDACs causes a decrease in protein acetylation status, resulting in decreased cognitive and memory function. Herein, we employed an integrative approach to analyze the transcriptomics data that established a potential relationship between AD and PD. Data preprocessing and analysis of four publicly available microarray datasets revealed 10 HUB proteins, namely, CDC42, CD44, FGFR1, MYO5A, NUMA1, TUBB4B, ARHGEF9, USP5, INPP5D, and NUP93, that may be involved in the shared mechanism of AD and PD pathogenesis. Further, we identified the relationship between the HUB proteins and transcription factors that could be involved in the overlapping mechanism of AD and PD. CREB1 and HINFP were the crucial regulatory transcription factors that were involved in the AD and PD crosstalk. Further, lysine acetylation sites and HDAC enzyme prediction revealed the involvement of 15 and 27 potential lysine residues of CREB1 and HINFP, respectively. Our results highlighted the importance of HDAC1(K292) and HDAC6(K330) association with CREB1 and HINFP, respectively, in the AD and PD crosstalk. However, different datasets with a large number of samples and wet lab experimentation are required to validate and pinpoint the exact role of CREB1 and HINFP in the AD and PD crosstalk. It is also possible that the different datasets may or may not affect the results due to analysis parameters. In conclusion, our study potentially highlighted the crucial proteins, transcription factors, biological pathways, lysine residues, and HDAC enzymes shared between AD and PD at the molecular level. The findings can be used to study molecular studies to identify the possible relationship in the AD-PD crosstalk.

STALLION: a stacking-based ensemble learning framework for prokaryotic lysine acetylation site prediction

Protein post-translational modification (PTM) is an important regulatory mechanism that plays a key role in both normal and disease states. Acetylation on lysine residues is one of the most potent PTMs owing to its critical role in cellular metabolism and regulatory processes. Identifying protein lysine acetylation (Kace) sites is a challenging task in bioinformatics. To date, several machine learning-based methods for the in silico identification of Kace sites have been developed. Of those, a few are prokaryotic species-specific. Despite their attractive advantages and performances, these methods have certain limitations. Therefore, this study proposes a novel predictor STALLION (STacking-based Predictor for ProkAryotic Lysine AcetyLatION), containing six prokaryotic species-specific models to identify Kace sites accurately. To extract crucial patterns around Kace sites, we employed 11 different encodings representing three different characteristics. Subsequently, a systematic and rigorous feature selection approach was employed to identify the optimal feature set independently for five tree-based ensemble algorithms and built their respective baseline model for each species. Finally, the predicted values from baseline models were utilized and trained with an appropriate classifier using the stacking strategy to develop STALLION. Comparative benchmarking experiments showed that STALLION significantly outperformed existing predictor on independent tests. To expedite direct accessibility to the STALLION models, a user-friendly online predictor was implemented, which is available at: http://thegleelab.org/STALLION.

i6mA-VC: A Multi-Classifier Voting Method for the Computational Identification of DNA N6-methyladenine Sites

DNA N6-methyladenine (6 mA), as an essential component of epigenetic modification, cannot be neglected in genetic regulation mechanism. The efficient and accurate prediction of 6 mA sites is beneficial to the development of biological genetics. Biochemical experimental methods are considered to be time-consuming and laborious. Most of the established machine learning methods have a single dataset. Although some of them have achieved cross-species prediction, their results are not satisfactory. Therefore, we designed a novel statistical model called i6mA-VC to improve the accuracy for 6 mA sites. On the one hand, kmer and binary encoding are applied to extract features, and then gradient boosting decision tree (GBDT) embedded method is applied as the feature selection strategy. On the other hand, DNA sequences are represented by vectors through the feature extraction method of ring-function-hydrogen-chemical properties (RFHCP) and the feature selection strategy of ExtraTree. After fusing the two optimal features, a voting classifier based on gradient boosting decision tree (GBDT), light gradient boosting machine (LightGBM) and multilayer perceptron classifier (MLPC) is constructed for final classification and prediction. The accuracy of Rice dataset and M.musculus dataset with five-fold cross-validation are 0.888 and 0.967, respectively. The cross-species dataset is selected as independent testing dataset, and the accuracy reaches 0.848. Through rigorous experiments, it is demonstrated that the proposed predictor is convincing and applicable. The development of i6mA-VC predictor will become an effective way for the recognition of N6-methyladenine sites, and it will also be beneficial for biological geneticists to further study gene expression and DNA modification. In addition, an accessible web-server for i6mA-VC is available from http://www.zhanglab.site/ .

pCysMod: Prediction of Multiple Cysteine Modifications Based on Deep Learning Framework

Thiol groups on cysteines can undergo multiple post-translational modifications (PTMs), acting as a molecular switch to maintain redox homeostasis and regulating a series of cell signaling transductions. Identification of sophistical protein cysteine modifications is crucial for dissecting its underlying regulatory mechanism. Instead of a time-consuming and labor-intensive experimental method, various computational methods have attracted intense research interest due to their convenience and low cost. Here, we developed the first comprehensive deep learning based tool pCysMod for multiple protein cysteine modification prediction, including S-nitrosylation, S-palmitoylation, S-sulfenylation, S-sulfhydration, and S-sulfinylation. Experimentally verified cysteine sites curated from literature and sites collected by other databases and predicting tools were integrated as benchmark dataset. Several protein sequence features were extracted and united into a deep learning model, and the hyperparameters were optimized by particle swarm optimization algorithms. Cross-validations indicated our model showed excellent robustness and outperformed existing tools, which was able to achieve an average AUC of 0.793, 0.807, 0.796, 0.793, and 0.876 for S-nitrosylation, S-palmitoylation, S-sulfenylation, S-sulfhydration, and S-sulfinylation, demonstrating pCysMod was stable and suitable for protein cysteine modification prediction. Besides, we constructed a comprehensive protein cysteine modification prediction web server based on this model to benefit the researches finding the potential modification sites of their interested proteins, which could be accessed at http://pcysmod.omicsbio.info. This work will undoubtedly greatly promote the study of protein cysteine modification and contribute to clarifying the biological regulation mechanisms of cysteine modification within and among the cells.

Deep Learning in Proteomics

Proteomics, the study of all the proteins in biological systems, is becoming a data-rich science. Protein sequences and structures are comprehensively catalogued in online databases. With recent advancements in tandem mass spectrometry (MS) technology, protein expression and post-translational modifications (PTMs) can be studied in a variety of biological systems at the global scale. Sophisticated computational algorithms are needed to translate the vast amount of data into novel biological insights. Deep learning automatically extracts data representations at high levels of abstraction from data, and it thrives in data-rich scientific research domains. Here, a comprehensive overview of deep learning applications in proteomics, including retention time prediction, MS/MS spectrum prediction, de novo peptide sequencing, PTM prediction, major histocompatibility complex-peptide binding prediction, and protein structure prediction, is provided. Limitations and the future directions of deep learning in proteomics are also discussed. This review will provide readers an overview of deep learning and how it can be used to analyze proteomics data.

Pharmacological targets and molecular mechanisms of plumbagin to treat colorectal cancer: A systematic pharmacology study

Plumbagin (PL) pharmacologically plays the anti-proliferative effects in cancer cells, including effective suppression of colorectal cancer (CRC). However, the exact molecular mechanism of PL to treat CRC remains unclear. Using available SwissTargetPrediction and SuperPred databases, the anti-cancer biotargets of PL were identified, and the CRC-diseased targets were obtained through a DisGeNET database. The biological processes, and signaling pathways of PL to treat CRC were identified and visualized. Further, clinical and cell culture data were used to validate some bioinformatic findings. As shown in bioinformatics findings, 64 predictive biotargets of PL to treat CRC were collected, and 7 most important biotargets of tumor protein p53 (TP53), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), mitogen-activated protein kinase 1 (MAPK1), E1A-associated protein p300 (EP300), poly (ADP-ribose) polymerase 1 (PARP1), nuclear factor kappa p65 protein (RELA), Bcl-2 like protein 1 (BCL2L1) were identified respectively. In addition, top 20 functional biological processes, signaling pathways of PL to treat CRC were screened and prioritized. In human study, CRC samples showed elevated expressions of neoplastic MAPK1, PARP1 mRNAs and reduced EP300 mRNA level. In cell culture study, PL-treated CRC cells resulted in down-regulated MAPK1, PARP1 mRNA expressions and up-regulation of EP300 mRNA level, characterized with suppressed cell proliferation. Taken together, the therapeutic biotargets and molecular mechanisms of PL to treat CRC were screened and identified by using a systematic pharmacology analysis, and some bioinformatic findings were validated in clinical and cell line experiments. Potentially, these hub biotargets may be the biomarkers for CRC detection and treatment.