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.

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.

Integrative machine learning framework for the identification of cell-specific enhancers from the human genome.

Enhancers are deoxyribonucleic acid (DNA) fragments which when bound by transcription factors enhance the transcription of related genes. Due to its sporadic distribution and similar fractions, identification of enhancers from the human genome seems a daunting task. Compared to the traditional experimental approaches, computational methods with easy-to-use platforms could be efficiently applied to annotate enhancers' functions and physiological roles. In this aspect, several bioinformatics tools have been developed to identify enhancers. Despite their spectacular performances, existing methods have certain drawbacks and limitations, including fixed length of sequences being utilized for model development and cell-specificity negligence. A novel predictor would be beneficial in the context of genome-wide enhancer prediction by addressing the above-mentioned issues. In this study, we constructed new datasets for eight different cell types. Utilizing these data, we proposed an integrative machine learning (ML)-based framework called Enhancer-IF for identifying cell-specific enhancers. Enhancer-IF comprehensively explores a wide range of heterogeneous features with five commonly used ML methods (random forest, extremely randomized tree, multilayer perceptron, support vector machine and extreme gradient boosting). Specifically, these five classifiers were trained with seven encodings and obtained 35 baseline models. The output of these baseline models was integrated and again inputted to five classifiers for the construction of five meta-models. Finally, the integration of five meta-models through ensemble learning improved the model robustness. Our proposed approach showed an excellent prediction performance compared to the baseline models on both training and independent datasets in different cell types, thus highlighting the superiority of our approach in the identification of the enhancers. We assume that Enhancer-IF will be a valuable tool for screening and identifying potential enhancers from the human DNA sequences.

Computational prediction of species-specific yeast DNA replication origin via iterative feature representation.

Deoxyribonucleic acid replication is one of the most crucial tasks taking place in the cell, and it has to be precisely regulated. This process is initiated in the replication origins (ORIs), and thus it is essential to identify such sites for a deeper understanding of the cellular processes and functions related to the regulation of gene expression. Considering the important tasks performed by ORIs, several experimental and computational approaches have been developed in the prediction of such sites. However, existing computational predictors for ORIs have certain curbs, such as building only single-feature encoding models, limited systematic feature engineering efforts and failure to validate model robustness. Hence, we developed a novel species-specific yeast predictor called yORIpred that accurately identify ORIs in the yeast genomes. To develop yORIpred, we first constructed optimal 40 baseline models by exploring eight different sequence-based encodings and five different machine learning classifiers. Subsequently, the predicted probability of 40 models was considered as the novel feature vector and carried out iterative feature learning approach independently using five different classifiers. Our systematic analysis revealed that the feature representation learned by the support vector machine algorithm (yORIpred) could well discriminate the distribution characteristics between ORIs and non-ORIs when compared with the other four algorithms. Comprehensive benchmarking experiments showed that yORIpred achieved superior and stable performance when compared with the existing predictors on the same training datasets. Furthermore, independent evaluation showcased the best and accurate performance of yORIpred thus underscoring the significance of iterative feature representation. To facilitate the users in obtaining their desired results without undergoing any mathematical, statistical or computational hassles, we developed a web server for the yORIpred predictor, which is available at: http://thegleelab.org/yORIpred.

Meta-i6mA: an interspecies predictor for identifying DNA N6-methyladenine sites of plant genomes by exploiting informative features in an integrative machine-learning framework.

DNA N6-methyladenine (6mA) represents important epigenetic modifications, which are responsible for various cellular processes. The accurate identification of 6mA sites is one of the challenging tasks in genome analysis, which leads to an understanding of their biological functions. To date, several species-specific machine learning (ML)-based models have been proposed, but majority of them did not test their model to other species. Hence, their practical application to other plant species is quite limited. In this study, we explored 10 different feature encoding schemes, with the goal of capturing key characteristics around 6mA sites. We selected five feature encoding schemes based on physicochemical and position-specific information that possesses high discriminative capability. The resultant feature sets were inputted to six commonly used ML methods (random forest, support vector machine, extremely randomized tree, logistic regression, naïve Bayes and AdaBoost). The Rosaceae genome was employed to train the above classifiers, which generated 30 baseline models. To integrate their individual strength, Meta-i6mA was proposed that combined the baseline models using the meta-predictor approach. In extensive independent test, Meta-i6mA showed high Matthews correlation coefficient values of 0.918, 0.827 and 0.635 on Rosaceae, rice and Arabidopsis thaliana, respectively and outperformed the existing predictors. We anticipate that the Meta-i6mA can be applied across different plant species. Furthermore, we developed an online user-friendly web server, which is available at http://kurata14.bio.kyutech.ac.jp/Meta-i6mA/.

HLPpred-Fuse: improved and robust prediction of hemolytic peptide and its activity by fusing multiple feature representation.

MOTIVATION: Therapeutic peptides failing at clinical trials could be attributed to their toxicity profiles like hemolytic activity, which hamper further progress of peptides as drug candidates. The accurate prediction of hemolytic peptides (HLPs) and its activity from the given peptides is one of the challenging tasks in immunoinformatics, which is essential for drug development and basic research. Although there are a few computational methods that have been proposed for this aspect, none of them are able to identify HLPs and their activities simultaneously. RESULTS: In this study, we proposed a two-layer prediction framework, called HLPpred-Fuse, that can accurately and automatically predict both hemolytic peptides (HLPs or non-HLPs) as well as HLPs activity (high and low). More specifically, feature representation learning scheme was utilized to generate 54 probabilistic features by integrating six different machine learning classifiers and nine different sequence-based encodings. Consequently, the 54 probabilistic features were fused to provide sufficiently converged sequence information which was used as an input to extremely randomized tree for the development of two final prediction models which independently identify HLP and its activity. Performance comparisons over empirical cross-validation analysis, independent test and case study against state-of-the-art methods demonstrate that HLPpred-Fuse consistently outperformed these methods in the identification of hemolytic activity. AVAILABILITY AND IMPLEMENTATION: For the convenience of experimental scientists, a web-based tool has been established at http://thegleelab.org/HLPpred-Fuse. CONTACT: glee@ajou.ac.kr or watshara.sho@mahidol.ac.th or bala@ajou.ac.kr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Silica-coated magnetic nanoparticles activate microglia and induce neurotoxic D-serine secretion.

BACKGROUND: Nanoparticles have been studied for brain imaging, diagnosis, and drug delivery owing to their versatile properties due to their small sizes. However, there are growing concerns that nanoparticles may exert toxic effects in the brain. In this study, we assessed direct nanotoxicity on microglia, the resident macrophages of the central nervous system, and indirect toxicity on neuronal cells exerted by silica-coated magnetic nanoparticles containing rhodamine B isothiocyanate dye [MNPs@SiO2(RITC)]. METHODS: We investigated MNPs@SiO2(RITC)-induced biological changes in BV2 murine microglial cells via RNA-sequencing-based transcriptome analysis and gas chromatography-mass spectrometry-based intracellular and extracellular amino acid profiling. Morphological changes were analyzed by transmission electron microscopy. Indirect effects of MNPs@SiO2(RITC) on neuronal cells were assessed by Transwell-based coculture with MNPs@SiO2(RITC)-treated microglia. MNPs@SiO2(RITC)-induced biological changes in the mouse brain in vivo were examined by immunohistochemical analysis. RESULTS: BV2 murine microglial cells were morphologically activated and the expression of Iba1, an activation marker protein, was increased after MNPs@SiO2(RITC) treatment. Transmission electron microscopy analysis revealed lysosomal accumulation of MNPs@SiO2(RITC) and the formation of vesicle-like structures in MNPs@SiO2(RITC)-treated BV2 cells. The expression of several genes related to metabolism and inflammation were altered in 100 µg/ml MNPs@SiO2(RITC)-treated microglia when compared with that in non-treated (control) and 10 µg/ml MNPs@SiO2(RITC)-treated microglia. Combined transcriptome and amino acid profiling analyses revealed that the transport of serine family amino acids, including glycine, cysteine, and serine, was enhanced. However, only serine was increased in the growth medium of activated microglia; especially, excitotoxic D-serine secretion from primary rat microglia was the most strongly enhanced. Activated primary microglia reduced intracellular ATP levels and proteasome activity in cocultured neuronal cells, especially in primary cortical neurons, via D-serine secretion. Moreover, ubiquitinated proteins accumulated and inclusion bodies were increased in primary dopaminergic and cortical neurons cocultured with activated primary microglia. In vivo, MNPs@SiO2(RITC), D-serine, and ubiquitin aggresomes were distributed in the MNPs@SiO2(RITC)-treated mouse brain. CONCLUSIONS: MNPs@SiO2(RITC)-induced activation of microglia triggers excitotoxicity in neurons via D-serine secretion, highlighting the importance of neurotoxicity mechanisms incurred by nanoparticle-induced microglial activation.

Silica-coated magnetic-nanoparticle-induced cytotoxicity is reduced in microglia by glutathione and citrate identified using integrated omics.

BACKGROUND: Nanoparticles have been utilized in brain research and therapeutics, including imaging, diagnosis, and drug delivery, owing to their versatile properties compared to bulk materials. However, exposure to nanoparticles leads to their accumulation in the brain, but drug development to counteract this nanotoxicity remains challenging. To date, concerns have risen about the potential toxicity to the brain associated with nanoparticles exposure via penetration of the brain blood barrier to address this issue. METHODS: Here the effect of silica-coated-magnetic nanoparticles containing the rhodamine B isothiocyanate dye [MNPs@SiO2(RITC)] were assessed on microglia through toxicological investigation, including biological analysis and integration of transcriptomics, proteomics, and metabolomics. MNPs@SiO2(RITC)-induced biological changes, such as morphology, generation of reactive oxygen species, intracellular accumulation of MNPs@SiO2(RITC) using transmission electron microscopy, and glucose uptake efficiency, were analyzed in BV2 murine microglial cells. Each omics data was collected via RNA-sequencing-based transcriptome analysis, liquid chromatography-tandem mass spectrometry-based proteome analysis, and gas chromatography- tandem mass spectrometry-based metabolome analysis. The three omics datasets were integrated and generated as a single network using a machine learning algorithm. Nineteen compounds were screened and predicted their effects on nanotoxicity within the triple-omics network. RESULTS: Intracellular reactive oxygen species production, an inflammatory response, and morphological activation of cells were greater, but glucose uptake was lower in MNPs@SiO2(RITC)-treated BV2 microglia and primary rat microglia in a dose-dependent manner. Expression of 121 genes (from 41,214 identified genes), and levels of 45 proteins (from 5918 identified proteins) and 17 metabolites (from 47 identified metabolites) related to the above phenomena changed in MNPs@SiO2(RITC)-treated microglia. A combination of glutathione and citrate attenuated nanotoxicity induced by MNPs@SiO2(RITC) and ten other nanoparticles in vitro and in the murine brain, protecting mostly the hippocampus and thalamus. CONCLUSIONS: Combination of glutathione and citrate can be one of the candidates for nanotoxicity alleviating drug against MNPs@SiO2(RITC) induced detrimental effect, including elevation of intracellular reactive oxygen species level, activation of microglia, and reduction in glucose uptake efficiency. In addition, our findings indicate that an integrated triple omics approach provides useful and sensitive toxicological assessment for nanoparticles and screening of drug for nanotoxicity.

Decrease in membrane fluidity and traction force induced by silica-coated magnetic nanoparticles.

BACKGROUND: Nanoparticles are being increasingly used in biomedical applications owing to their unique physical and chemical properties and small size. However, their biophysical assessment and evaluation of side-effects remain challenging. We addressed this issue by investigating the effects of silica-coated magnetic nanoparticles containing rhodamine B isothiocyanate [MNPs@SiO2(RITC)] on biophysical aspects, such as membrane fluidity and traction force of human embryonic kidney 293 (HEK293) cells. We further extended our understanding on the biophysical effects of nanoparticles on cells using a combination of metabolic profiling and transcriptomic network analysis. RESULTS: Overdose (1.0 µg/µL) treatment with MNPs@SiO2(RITC) induced lipid peroxidation and decreased membrane fluidity in HEK293 cells. In addition, HEK293 cells were morphologically shrunk, and their aspect ratio was significantly decreased. We found that each traction force (measured in micropillar) was increased, thereby increasing the total traction force in MNPs@SiO2(RITC)-treated HEK293 cells. Due to the reduction in membrane fluidity and elevation of traction force, the velocity of cell movement was also significantly decreased. Moreover, intracellular level of adenosine triphosphate (ATP) was also decreased in a dose-dependent manner upon treatment with MNPs@SiO2(RITC). To understand these biophysical changes in cells, we analysed the transcriptome and metabolic profiles and generated a metabotranscriptomics network, which revealed relationships among peroxidation of lipids, focal adhesion, cell movement, and related genes and metabolites. Furthermore, in silico prediction of the network showed increment in the peroxidation of lipids and suppression of focal adhesion and cell movement. CONCLUSION: Taken together, our results demonstrated that overdose of MNPs@SiO2(RITC) impairs cellular movement, followed by changes in the biophysical properties of cells, thus highlighting the need for biophysical assessment of nanoparticle-induced side-effects.

Machine intelligence in peptide therapeutics: A next-generation tool for rapid disease screening.

Discovery and development of biopeptides are time-consuming, laborious, and dependent on various factors. Data-driven computational methods, especially machine learning (ML) approach, can rapidly and efficiently predict the utility of therapeutic peptides. ML methods offer an array of tools that can accelerate and enhance decision making and discovery for well-defined queries with ample and sophisticated data quality. Various ML approaches, such as support vector machines, random forest, extremely randomized tree, and more recently deep learning methods, are useful in peptide-based drug discovery. These approaches leverage the peptide data sets, created via high-throughput sequencing and computational methods, and enable the prediction of functional peptides with increased levels of accuracy. The use of ML approaches in the development of peptide-based therapeutics is relatively recent; however, these techniques are already revolutionizing protein research by unraveling their novel therapeutic peptide functions. In this review, we discuss several ML-based state-of-the-art peptide-prediction tools and compare these methods in terms of their algorithms, feature encodings, prediction scores, evaluation methodologies, and software utilities. We also assessed the prediction performance of these methods using well-constructed independent data sets. In addition, we discuss the common pitfalls and challenges of using ML approaches for peptide therapeutics. Overall, we show that using ML models in peptide research can streamline the development of targeted peptide therapies.

mAHTPred: a sequence-based meta-predictor for improving the prediction of anti-hypertensive peptides using effective feature representation.

MOTIVATION: Cardiovascular disease is the primary cause of death globally accounting for approximately 17.7 million deaths per year. One of the stakes linked with cardiovascular diseases and other complications is hypertension. Naturally derived bioactive peptides with antihypertensive activities serve as promising alternatives to pharmaceutical drugs. So far, there is no comprehensive analysis, assessment of diverse features and implementation of various machine-learning (ML) algorithms applied for antihypertensive peptide (AHTP) model construction. RESULTS: In this study, we utilized six different ML algorithms, namely, Adaboost, extremely randomized tree (ERT), gradient boosting (GB), k-nearest neighbor, random forest (RF) and support vector machine (SVM) using 51 feature descriptors derived from eight different feature encodings for the prediction of AHTPs. While ERT-based trained models performed consistently better than other algorithms regardless of various feature descriptors, we treated them as baseline predictors, whose predicted probability of AHTPs was further used as input features separately for four different ML-algorithms (ERT, GB, RF and SVM) and developed their corresponding meta-predictors using a two-step feature selection protocol. Subsequently, the integration of four meta-predictors through an ensemble learning approach improved the balanced prediction performance and model robustness on the independent dataset. Upon comparison with existing methods, mAHTPred showed superior performance with an overall improvement of approximately 6-7% in both benchmarking and independent datasets. AVAILABILITY AND IMPLEMENTATION: The user-friendly online prediction tool, mAHTPred is freely accessible at http://thegleelab.org/mAHTPred. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Silica-coated magnetic nanoparticles induce glucose metabolic dysfunction in vitro via the generation of reactive oxygen species.

Nanoparticles are a useful material in biomedicine given their unique properties and biocompatibility; however, there is increasing concern regarding the potential toxicity of nanoparticles with respect to cell metabolism. Some evidence suggests that nanoparticles can disrupt glucose and energy homeostasis. In this study, we investigated the metabolomic, transcriptomic, and integrated effects of silica-coated magnetic nanoparticles containing rhodamine B isothiocyanate dye [MNPs@SiO2(RITC)] on glucose metabolism in human embryonic kidney 293 (HEK293) cells. Using gas chromatography-tandem mass spectrometry, we analysed the metabolite profiles of 14 organic acids (OAs), 20 amino acids (AAs), and 13 fatty acids (FAs) after treatment with 0.1 or 1.0 µg/µl MNPs@SiO2(RITC) for 12 h. The metabolic changes were highly related to reactive oxygen species (ROS) generation and glucose metabolism. Additionally, effects on the combined metabolome and transcriptome or "metabotranscriptomic network" indicated a relationship between ROS generation and glucose metabolic dysfunction. In the experimental validation, MNPs@SiO2(RITC) treatment significantly decreased the amount of glucose in cells and was associated with a reduction in glucose uptake efficiency. Decreased glucose uptake efficiency was also related to ROS generation and impaired glucose metabolism in the metabotranscriptomic network. Our results suggest that exposure to high concentrations of MNPs@SiO2(RITC) produces maladaptive alterations in glucose metabolism and specifically glucose uptake as well as related metabolomic and transcriptomic disturbances via increased ROS generation. These findings further indicate that an integrated metabotranscriptomics approach provides useful and sensitive toxicological assessment for nanoparticles.

A Molecular Dynamics Approach to Explore the Intramolecular Signal Transduction of PPAR-α.

Dynamics and functions of the peroxisome proliferator-activated receptor (PPAR)-α are modulated by the types of ligands that bind to the orthosteric sites. While several X-ray crystal structures of PPAR-α have been determined in their agonist-bound forms, detailed structural information in their apo and antagonist-bound states are still lacking. To address these limitations, we apply unbiased molecular dynamics simulations to three different PPAR-α systems to determine their modulatory mechanisms. Herein, we performed hydrogen bond and essential dynamics analyses to identify the important residues involved in polar interactions and conformational structural variations, respectively. Furthermore, betweenness centrality network analysis was carried out to identify key residues for intramolecular signaling. The differences observed in the intramolecular signal flow between apo, agonist- and antagonist-bound forms of PPAR-α will be useful for calculating maps of information flow and identifying key residues crucial for signal transductions. The predictions derived from our analysis will be of great help to medicinal chemists in the design of effective PPAR-α modulators and additionally in understanding their regulation and signal transductions.

Evolution of In Silico Strategies for Protein-Protein Interaction Drug Discovery.

The advent of advanced molecular modeling software, big data analytics, and high-speed processing units has led to the exponential evolution of modern drug discovery and better insights into complex biological processes and disease networks. This has progressively steered current research interests to understanding protein-protein interaction (PPI) systems that are related to a number of relevant diseases, such as cancer, neurological illnesses, metabolic disorders, etc. However, targeting PPIs are challenging due to their "undruggable" binding interfaces. In this review, we focus on the current obstacles that impede PPI drug discovery, and how recent discoveries and advances in in silico approaches can alleviate these barriers to expedite the search for potential leads, as shown in several exemplary studies. We will also discuss about currently available information on PPI compounds and systems, along with their usefulness in molecular modeling. Finally, we conclude by presenting the limits of in silico application in drug discovery and offer a perspective in the field of computer-aided PPI drug discovery.

Harnessing the Therapeutic Potential of Capsaicin and Its Analogues in Pain and Other Diseases.

Capsaicin is the most predominant and naturally occurring alkamide found in Capsicum fruits. Since its discovery in the 19th century, the therapeutic roles of capsaicin have been well characterized. The potential applications of capsaicin range from food flavorings to therapeutics. Indeed, capsaicin and few of its analogues have featured in clinical research covered by more than a thousand patents. Previous records suggest pleiotropic pharmacological activities of capsaicin such as an analgesic, anti-obesity, anti-pruritic, anti-inflammatory, anti-apoptotic, anti-cancer, anti-oxidant, and neuro-protective functions. Moreover, emerging data indicate its clinical significance in treating vascular-related diseases, metabolic syndrome, and gastro-protective effects. The dearth of potent drugs for management of such disorders necessitates the urge for further research into the pharmacological aspects of capsaicin. This review summarizes the historical background, source, structure and analogues of capsaicin, and capsaicin-triggered TRPV1 signaling and desensitization processes. In particular, we will focus on the therapeutic roles of capsaicin and its analogues in both normal and pathophysiological conditions.

Unveiling local and global conformational changes and allosteric communications in SOD1 systems using molecular dynamics simulation and network analyses.

BACKGROUND: Amyotrophic lateral sclerosis (ALS) is a serious neurodegenerative disorder affecting nerve cells in the brain and spinal cord that is caused by mutations in the superoxide dismutase 1 (SOD1) enzyme. ALS-related mutations cause misfolding, dimerisation instability, and increased formation of aggregates. The underlying allosteric mechanisms, however, remain obscure as far as details of their fundamental atomistic structure are concerned. Hence, this gap in knowledge limits the development of novel SOD1 inhibitors and the understanding of how disease-associated mutations in distal sites affect enzyme activity. METHODS: We combined microsecond-scale based unbiased molecular dynamics (MD) simulation with network analysis to elucidate the local and global conformational changes and allosteric communications in SOD1 Apo (unmetallated form), Holo, Apo_CallA (mutant and unmetallated form), and Holo_CallA (mutant form) systems. To identify hotspot residues involved in SOD1 signalling and allosteric communications, we performed network centrality, community network, and path analyses. RESULTS: Structural analyses showed that unmetallated SOD1 systems and cysteine mutations displayed large structural variations in the catalytic sites, affecting structural stability. Inter- and intra H-bond analyses identified several important residues crucial for maintaining interfacial stability, structural stability, and enzyme catalysis. Dynamic motion analysis demonstrated more balanced atomic displacement and highly correlated motions in the Holo system. The rationale for structural disparity observed in the disulfide bond formation and R143 configuration in Apo and Holo systems were elucidated using distance and dihedral probability distribution analyses. CONCLUSION: Our study highlights the efficiency of combining extensive MD simulations with network analyses to unravel the features of protein allostery.

SEP-AlgPro: An efficient allergen prediction tool utilizing traditional machine learning and deep learning techniques with protein language model features.

Allergy is a hypersensitive condition in which individuals develop objective symptoms when exposed to harmless substances at a dose that would cause no harm to a "normal" person. Most current computational methods for allergen identification rely on homology or conventional machine learning using limited set of feature descriptors or validation on specific datasets, making them inefficient and inaccurate. Here, we propose SEP-AlgPro for the accurate identification of allergen protein from sequence information. We analyzed 10 conventional protein-based features and 14 different features derived from protein language models to gauge their effectiveness in differentiating allergens from non-allergens using 15 different classifiers. However, the final optimized model employs top 10 feature descriptors with top seven machine learning classifiers. Results show that the features derived from protein language models exhibit superior discriminative capabilities compared to traditional feature sets. This enabled us to select the most discriminatory baseline models, whose predicted outputs were aggregated and used as input to a deep neural network for the final allergen prediction. Extensive case studies showed that SEP-AlgPro outperforms state-of-the-art predictors in accurately identifying allergens. A user-friendly web server was developed and made freely available at https://balalab-skku.org/SEP-AlgPro/, making it a powerful tool for identifying potential allergens.

ADP-Fuse: A novel two-layer machine learning predictor to identify antidiabetic peptides and diabetes types using multiview information.

Diabetes mellitus has become a major public health concern associated with high mortality and reduced life expectancy and can cause blindness, heart attacks, kidney failure, lower limb amputations, and strokes. A new generation of antidiabetic peptides (ADPs) that act on ß-cells or T-cells to regulate insulin production is being developed to alleviate the effects of diabetes. However, the lack of effective peptide-mining tools has hampered the discovery of these promising drugs. Hence, novel computational tools need to be developed urgently. In this study, we present ADP-Fuse, a novel two-layer prediction framework capable of accurately identifying ADPs or non-ADPs and categorizing them into type 1 and type 2 ADPs. First, we comprehensively evaluated 22 peptide sequence-derived features coupled with eight notable machine learning algorithms. Subsequently, the most suitable feature descriptors and classifiers for both layers were identified. The output of these single-feature models, embedded with multiview information, was trained with an appropriate classifier to provide the final prediction. Comprehensive cross-validation and independent tests substantiate that ADP-Fuse surpasses single-feature models and the feature fusion approach for the prediction of ADPs and their types. In addition, the SHapley Additive exPlanation method was used to elucidate the contributions of individual features to the prediction of ADPs and their types. Finally, a user-friendly web server for ADP-Fuse was developed and made publicly accessible (https://balalab-skku.org/ADP-Fuse), enabling the swift screening and identification of novel ADPs and their types. This framework is expected to contribute significantly to antidiabetic peptide identification.

Therapeutic applications of nucleic acids and their analogues in Toll-like receptor signaling.

Toll-like receptors (TLRs) belong to a family of innate immune receptors that detect and clear invading microbial pathogens. Specifically intracellular TLRs such as TLR3, TLR7, TLR8 and TLR9 recognize nucleic acids such as double-stranded RNA, single-stranded RNA and CpG DNA respectively derived from microbial components. Upon infection, nucleic acid sensing TLRs signal within endosomal compartment triggering the induction of essential proinflammatory cytokines and type I interferons to initiate innate immune responses thereby leading to a critical role in the development of adaptive immune responses. Thus, stimulation of TLRs by nucleic acids is a promising area of research for the development of novel therapeutic strategies against pathogenic infection, allergies, malignant neoplasms and autoimmunity. This review summarizes the therapeutic applications of nucleic acids or nucleic acid analogues through the modulation of TLR signaling pathways.

Amyotrophic lateral sclerosis disease-related mutations disrupt the dimerization of superoxide dismutase 1 - A comparative molecular dynamics simulation study.

More than 150 genes are involved in amyotrophic lateral sclerosis (ALS), with superoxide dismutase 1 (SOD1) being one of the most studied. Mutations in SOD1 gene, which encodes the enzyme SOD1 is the second most prevalent and studied cause of familial ALS. SOD1 is a ubiquitous, homodimeric metalloenzyme that forms a critical component of the cellular defense against reactive oxygen species. Several mutations in the SOD1 enzyme cause misfolding, dimerization instability, and increased aggregate formation in ALS. However, there is a lack of information on the dimerization of SOD1 monomers and the mechanistic underpinnings on how the pathogenic mutations disrupt the dimerization mechanism. Here, we presented microsecond-scale molecular dynamics (MD) simulations to unravel how interface-based mutations compromise SOD1 dimerization and provide mechanistic understanding into the corresponding process using WT and three interface-based mutant systems (A4V, T54R, and I113T). Structural stability analysis showed that the mutant systems displayed disparate variations in the catalytic sites which may directly alter the stability and activity of the SOD1 enzyme. Based on the dynamic network analysis and principal component analysis, it has been identified that the mutations weakened the correlated motions along the dimer interface and altered the protein conformational behavior, thus weakening the stability of dimer formation. Moreover, the simulation results identified crucial residues such as G51, D52, G114, I151, and Q153 in establishing the dimerization interaction network, which were weakened or absent in the presence of interfacial mutants. Surface potential analysis on mutant systems also displayed changes in the dimerization potential, thus showing the unfavorable dimer formation. Furthermore, network analysis identified the hotspot residues necessary for SOD1 signal transduction which were surprisingly found in the catalytic sites rather than the anticipated dimerization interface.