Multiple diseases (acute or chronic events) occur together in a patient, which refers to the disease comorbidities, because of the multi ways associations among diseases. Due to shared genetic, molecular, environmental, and lifestyle-based risk factors, many diseases are comorbid in the same patient. Methods for integrating multiple types of omics data play an important role to identify integrative biomarkers for stratification of patients into groups with different clinical outcomes. Moreover, integrated omics and clinical information may potentially improve prediction accuracy of disease comorbidities. However, there is a lack of effective and efficient bioinformatics and statistical software for true integrative data analysis. With the availability of the wide spread huge omics, phenotype and ontology information, it is becoming more and more practical to help doctors in clinical diagnostics and comorbidity prediction by providing appropriate software tool. We developed an R software POGO to compute novel estimators of the disease comorbidity risks and patient stratification. Starting from an initial diagnosis, omics and clinical data of a patient the software identifies the association risk of disease comorbidities. The input of this software is the initial diagnosis of a patient and the output provides evidence of disease comorbidities. The functions of POGO offer flexibility for diagnostic applications to predict disease comorbidities, and can be easily integrated to high-throughput and clinical data analysis pipelines. POGO is compliant with the Bioconductor standard and it is freely available at www.cl.cam.ac.uk/~mam211/POGO/.

This study represents the first large-scale study on the chemical space of inhibitors of dipeptidyl peptidase-4 (DPP4), which is a potential therapeutic protein target for the treatment of diabetes mellitus. Herein, a large set of 2,937 compounds evaluated for their ability to inhibit DPP4 was compiled from the literature. Molecular descriptors were generated from the geometrically optimized low-energy conformers of these compounds at the semiempirical AM1 level. The origins of DPP4 inhibitory activity were elucidated from computed molecular descriptors that accounted for the unique physicochemical properties inherently present in the active and inactive sets of compounds as defined by their respective half maximal inhibitory concentration values of less than 1 μM and greater than 10 μM, respectively. Decision tree analysis revealed the importance of molecular weight, total energy of a molecule, topological polar surface area, lowest unoccupied molecular orbital, and number of hydrogen-bond donors, which correspond to molecular size, energy, surface polarity, electron acceptors, and hydrogen bond donors, respectively. The prediction model was subjected to rigorous independent testing via three external sets. Scaffold and chemical fragment analysis was also performed on these active and inactive sets of compounds to shed light on the distinguishing features of the functional moieties. Docking of representative active DPP4 inhibitors was also performed to unravel key interacting residues. The results of this study are anticipated to be useful in guiding the rational design of novel and robust DPP4 inhibitors for the treatment of diabetes.

P-glycoprotein (Pgp) is a drug transporter that plays important roles in multidrug resistance and drug pharmacokinetics. The inhibition of Pgp has become a notable strategy for combating multidrug-resistant cancers and improving therapeutic outcomes. However, the polyspecific nature of Pgp, together with inconsistent results in experimental assays, renders the determination of endpoints for Pgp-interacting compounds a great challenge. In this study, the classification of a large set of 2,477 Pgp-interacting compounds (i.e., 1341 inhibitors, 913 non-inhibitors, 197 substrates and 26 non-substrates) was performed using several machine learning methods (i.e., decision tree induction, artificial neural network modelling and support vector machine) as a function of their physicochemical properties. The models provided good predictive performance, producing MCC values in the range of 0.739-1 for internal cross-validation and 0.665-1 for external validation. The study provided simple and interpretable models for important properties that influence the activity of Pgp-interacting compounds, which are potentially beneficial for screening and rational design of Pgp inhibitors that are of clinical importance.

An outbreak, caused by an RNA virus, SARS-CoV-2 named COVID-19 has become pandemic with a magnitude which is daunting to all public health institutions in the absence of specific antiviral treatment. Surface glycoprotein and nucleocapsid phosphoprotein are two important proteins of this virus facilitating its entry into host cell and genome replication. Small interfering RNA (siRNA) is a prospective tool of the RNA interference (RNAi) pathway for the control of human viral infections by suppressing viral gene expression through hybridization and neutralization of target complementary mRNA. So, in this study, the power of RNA interference technology was harnessed to develop siRNA molecules against specific target genes namely, nucleocapsid phosphoprotein gene and surface glycoprotein gene. Conserved sequence from 139 SARS-CoV-2 strains from around the globe was collected to construct 78 siRNA that can inactivate nucleocapsid phosphoprotein and surface glycoprotein genes. Finally, based on GC content, free energy of folding, free energy of binding, melting temperature, efficacy prediction and molecular docking analysis, 8 siRNA molecules were selected which are proposed to exert the best action. These predicted siRNAs should effectively silence the genes of SARS-CoV-2 during siRNA mediated treatment assisting in the response against SARS-CoV-2.

To identify key gene expression pathways altered with infection of the novel coronavirus SARS-CoV-2, we performed the largest comparative genomic and transcriptomic analysis to date. We compared the novel pandemic coronavirus SARS-CoV-2 with SARS-CoV and MERS-CoV, as well as influenza A strains H1N1, H3N2 and H5N1. Phylogenetic analysis confirms that SARS-CoV-2 is closely related to SARS-CoV at the level of the viral genome. RNAseq analyses demonstrate that human lung epithelial cell responses to SARS-CoV-2 infection are distinct. Extensive Gene Expression Omnibus literature screening and drug predictive analyses show that SARS-CoV-2 infection response pathways are closely related to those of SARS-CoV and respiratory syncytial virus infections. We validated SARS-CoV-2 infection response genes as disease-associated using Kaplan-Meier survival estimates in lung disease patient data. We also analysed COVID-19 patient peripheral blood samples, which identified signalling pathway concordance between the primary lung cell and blood cell infection responses.

Understanding of quorum-sensing peptides (QSPs) in their functional mechanism plays an essential role in finding new opportunities to combat bacterial infections by designing drugs. With the avalanche of the newly available peptide sequences in the post-genomic age, it is highly desirable to develop a computational model for efficient, rapid and high-throughput QSP identification purely based on the peptide sequence information alone. Although, few methods have been developed for predicting QSPs, their prediction accuracy and interpretability still requires further improvements. Thus, in this work, we proposed an accurate sequence-based predictor (called iQSP) and a set of interpretable rules (called IR-QSP) for predicting and analyzing QSPs. In iQSP, we utilized a powerful support vector machine (SVM) cooperating with 18 informative features from physicochemical properties (PCPs). Rigorous independent validation test showed that iQSP achieved maximum accuracy and MCC of 93.00% and 0.86, respectively. Furthermore, a set of interpretable rules IR-QSP was extracted by using random forest model and the 18 informative PCPs. Finally, for the convenience of experimental scientists, the iQSP web server was established and made freely available online. It is anticipated that iQSP will become a useful tool or at least as a complementary existing method for predicting and analyzing QSPs.

Cancer imposes a global health burden as it represents one of the leading causes of morbidity and mortality while also giving rise to significant economic burden owing to the associated expenditures for its monitoring and treatment. In spite of advancements in cancer therapy, the low success rate and recurrence of tumor has necessitated the ongoing search for new therapeutic agents. Aside from drugs based on small molecules and protein-based biopharmaceuticals, there has been an intense effort geared towards the development of peptide-based therapeutics owing to its favorable and intrinsic properties of being relatively small, highly selective, potent, safe and low in production costs. In spite of these advantages, there are several inherent weaknesses that are in need of attention in the design and development of therapeutic peptides. An abundance of data on bioactive and therapeutic peptides have been accumulated over the years and the burgeoning area of artificial intelligence has set the stage for the lucrative utilization of machine learning to make sense of these large and high-dimensional data. This review summarizes the current state-of-the-art on the application of machine learning for studying the bioactivity of anticancer peptides along with future outlook of the field. Data and R codes used in the analysis herein are available on GitHub at https://github.com/Shoombuatong2527/anticancer-peptides-review.

Cancer remains one of the major causes of death worldwide. Angiogenesis is crucial for the pathogenesis of various human diseases, especially solid tumors. The discovery of anti-angiogenic peptides is a promising therapeutic route for cancer treatment. Thus, reliably identifying anti-angiogenic peptides is extremely important for understanding their biophysical and biochemical properties that serve as the basis for the discovery of new anti-cancer drugs. This study aims to develop an efficient and interpretable computational model called TargetAntiAngio for predicting and characterizing anti-angiogenic peptides. TargetAntiAngio was developed using the random forest classifier in conjunction with various classes of peptide features. It was observed via an independent validation test that TargetAntiAngio can identify anti-angiogenic peptides with an average accuracy of 77.50% on an objective benchmark dataset. Comparisons demonstrated that TargetAntiAngio is superior to other existing methods. In addition, results revealed the following important characteristics of anti-angiogenic peptides: (i) disulfide bond forming Cys residues play an important role for inhibiting blood vessel proliferation; (ii) Cys located at the C-terminal domain can decrease endothelial formatting activity and suppress tumor growth; and (iii) Cyclic disulfide-rich peptides contribute to the inhibition of angiogenesis and cell migration, selectivity and stability. Finally, for the convenience of experimental scientists, the TargetAntiAngio web server was established and made freely available online.

Welding involves exposure to fumes, gases and radiant energy that can be hazardous to human health. Welding fumes (WFs) comprise a complex mixture of metallic oxides, silicates and fluorides that may result in different health effects. Inhalation of WFs in large quantities over a long periods may pose a risk of developing neurodegenerative diseases (NDGDs), but the nature of this risk is poorly understood. To address this we performed transcriptomic analysis to identify links between WF exposure and NDGDs.

Coronavirus Disease 2019 (COVID-19), although most commonly demonstrates respiratory symptoms, but there is a growing set of evidence reporting its correlation with the digestive tract and faeces. Interestingly, recent studies have shown the association of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection with gastrointestinal symptoms in infected patients but any sign of respiratory issues. Moreover, some studies have also shown that the presence of live SARS-CoV-2 virus in the faeces of patients with COVID-19. Therefore, the pathophysiology of digestive symptoms associated with COVID-19 has raised a critical need for comprehensive investigative efforts. To address this issue we have developed a bioinformatics pipeline involving a system biological framework to identify the effects of SARS-CoV-2 messenger RNA expression on deciphering its association with digestive symptoms in COVID-19 positive patients. Using two RNA-seq datasets derived from COVID-19 positive patients with celiac (CEL), Crohn's (CRO) and ulcerative colitis (ULC) as digestive disorders, we have found a significant overlap between the sets of differentially expressed genes from SARS-CoV-2 exposed tissue and digestive tract disordered tissues, reporting 7, 22 and 13 such overlapping genes, respectively. Moreover, gene set enrichment analysis, comprehensive analyses of protein-protein interaction network, gene regulatory network, protein-chemical agent interaction network revealed some critical association between SARS-CoV-2 infection and the presence of digestive disorders. The infectome, diseasome and comorbidity analyses also discover the influences of the identified signature genes in other risk factors of SARS-CoV-2 infection to human health. We hope the findings from this pathogenetic analysis may reveal important insights in deciphering the complex interplay between COVID-19 and digestive disorders and underpins its significance in therapeutic development strategy to combat against COVID-19 pandemic.

COVID-19 has emerged as global health threats. Chronic kidney disease (CKD) patients are immune-compromised and may have a high risk of infection by the SARS-CoV-2. We aimed to detect common transcriptomic signatures and pathways between COVID-19 and CKD by systems biology analysis. We analyzed transcriptomic data obtained from peripheral blood mononuclear cells (PBMC) infected with SARS-CoV-2 and PBMC of CKD patients. We identified 49 differentially expressed genes (DEGs) which were common between COVID-19 and CKD. The gene ontology and pathways analysis showed the DEGs were associated with "platelet degranulation", "regulation of wound healing", "platelet activation", "focal adhesion", "regulation of actin cytoskeleton" and "PI3K-Akt signalling pathway". The protein-protein interaction (PPI) network encoded by the common DEGs showed ten hub proteins (EPHB2, PRKAR2B, CAV1, ARHGEF12, HSP90B1, ITGA2B, BCL2L1, E2F1, TUBB1, and C3). Besides, we identified significant transcription factors and microRNAs that may regulate the common DEGs. We investigated protein-drug interaction analysis and identified potential drugs namely, aspirin, estradiol, rapamycin, and nebivolol. The identified common gene signature and pathways between COVID-19 and CKD may be therapeutic targets in COVID-19 patients with CKD comorbidity.

Alzheimer's disease (AD) is the commonest progressive neurodegenerative condition in humans, and is currently incurable. A wide spectrum of comorbidities, including other neurodegenerative diseases, are frequently associated with AD. How AD interacts with those comorbidities can be examined by analysing gene expression patterns in affected tissues using bioinformatics tools. We surveyed public data repositories for available gene expression data on tissue from AD subjects and from people affected by neurodegenerative diseases that are often found as comorbidities with AD. We then utilized large set of gene expression data, cell-related data and other public resources through an analytical process to identify functional disease links. This process incorporated gene set enrichment analysis and utilized semantic similarity to give proximity measures. We identified genes with abnormal expressions that were common to AD and its comorbidities, as well as shared gene ontology terms and molecular pathways. Our methodological pipeline was implemented in the R platform as an open-source package and available at the following link: https://github.com/unchowdhury/AD_comorbidity. The pipeline was thus able to identify factors and pathways that may constitute functional links between AD and these common comorbidities by which they affect each others development and progression. This pipeline can also be useful to identify key pathological factors and therapeutic targets for other diseases and disease interactions.

Schizophrenia (SCZ) is a psychiatric disorder characterized by both positive symptoms (i.e., psychosis) and negative symptoms (such as apathy, anhedonia, and poverty of speech). Epidemiological data show a high likelihood of early onset of type 2 diabetes mellitus (T2DM) in SCZ patients. However, the molecular processes that could explain the epidemiological association between SCZ and T2DM have not yet been characterized. Therefore, in the present study, we aimed to identify underlying common molecular pathogenetic processes and pathways between SCZ and T2DM. To this aim, we analyzed peripheral blood mononuclear cell (PBMC) transcriptomic data from SCZ and T2DM patients, and we detected 28 differentially expressed genes (DEGs) commonly modulated between SCZ and T2DM. Inflammatory-associated processes and membrane trafficking pathways as common biological processes were found to be in common between SCZ and T2DM. Analysis of the putative transcription factors involved in the regulation of the DEGs revealed that STAT1 (Signal Transducer and Activator of Transcription 1), RELA (v-rel reticuloendotheliosis viral oncogene homolog A (avian)), NFKB1 (Nuclear Factor Kappa B Subunit 1), and ERG (ETS-related gene) are involved in the expression of common DEGs in SCZ and T2DM. In conclusion, we provide core molecular signatures and pathways that are shared between SCZ and T2DM, which may contribute to the epidemiological association between them.

N7-methylguanosine (m7G) is an essential, ubiquitous, and positively charged modification at the 5' cap of eukaryotic mRNA, modulating its export, translation, and splicing processes. Although several machine learning (ML)-based computational predictors for m7G have been developed, all utilized specific computational framework. This study is the first instance we explored four different computational frameworks and identified the best approach. Based on that we developed a novel predictor, THRONE (A three-layer ensemble predictor for identifying human RNA N7-methylguanosine sites) to accurately identify m7G sites from the human genome. THRONE employs a wide range of sequence-based features inputted to several ML classifiers and combines these models through ensemble learning. The three-step ensemble learning is as follows: 54 baseline models were constructed in the first layer and the predicted probability of m7G was considered as a new feature vector for the sequential step. Subsequently, six meta-models were created using the new feature vector and their predicted probability was yet again considered as novel features. Finally, random forest was deemed as the best super classifier learner for the final prediction using a systematic approach incorporated with novel features. Interestingly, THRONE outperformed other existing methods in the prediction of m7G sites on both cross-validation analysis and independent evaluation. The proposed method is publicly accessible at: http://thegleelab.org/THRONE/ and expects to help the scientific community identify the putative m7G sites and formulate a novel testable biological hypothesis.

Corona virus disease (COVID-19) has been emerged as pandemic infectious disease. The recent epidemiological data suggest that the smokers are more vulnerable to infection with COVID-19; however, the influence of smoking (SMK) on the COVID-19 infected patients and the mortality is not known yet. In this study, we aimed to discern the influence of SMK on COVID-19 infected patients utilizing the transcriptomics data of COVID-19 infected lung epithelial cells and transcriptomics data smoking matched with controls from lung epithelial cells. The bioinformatics based analysis revealed the molecular insights into the level of transcriptional changes and pathways which are important to identify the impact of smoking on COVID-19 infection and prevalence. We compared differentially expressed genes (DEGs) between COVID-19 and SMK and 59 DEGs were identified as consistently dysregulated at transcriptomics levels. The correlation network analyses were constructed for these common genes using WGCNA R package to see the relationship among these genes. Integration of DEGs with network analysis (protein-protein interaction) showed the presence of 9 hub proteins as key so called "candidate hub proteins" overlapped between COVID-19 patients and SMK. The Gene Ontology and pathways analysis demonstrated the enrichment of inflammatory pathway such as IL-17 signaling pathway, Interleukin-6 signaling, TNF signaling pathway and MAPK1/MAPK3 signaling pathways that might be the therapeutic targets in COVID-19 for smoking persons. The identified genes, pathways, hubs genes, and their regulators might be considered for establishment of key genes and drug targets for SMK and COVID-19.

Anticancer peptides (ACPs) have emerged as a new class of therapeutic agent for cancer treatment due to their lower toxicity as well as greater efficacy, selectivity and specificity when compared to conventional small molecule drugs. However, the experimental identification of ACPs still remains a time-consuming and expensive endeavor. Therefore, it is desirable to develop and improve upon existing computational models for predicting and characterizing ACPs. In this study, we present a bioinformatics tool called the ACPred, which is an interpretable tool for the prediction and characterization of the anticancer activities of peptides. ACPred was developed by utilizing powerful machine learning models (support vector machine and random forest) and various classes of peptide features. It was observed by a jackknife cross-validation test that ACPred can achieve an overall accuracy of 95.61% in identifying ACPs. In addition, analysis revealed the following distinguishing characteristics that ACPs possess: (i) hydrophobic residue enhances the cationic properties of α-helical ACPs resulting in better cell penetration; (ii) the amphipathic nature of the α-helical structure plays a crucial role in its mechanism of cytotoxicity; and (iii) the formation of disulfide bridges on β-sheets is vital for structural maintenance which correlates with its ability to kill cancer cells. Finally, for the convenience of experimental scientists, the ACPred web server was established and made freely available online.

This work focuses on the computational modelling of osteomyelitis, a bone pathology caused by bacteria infection (mostly Staphylococcus aureus). The infection alters the RANK/RANKL/OPG signalling dynamics that regulates osteoblasts and osteoclasts behaviour in bone remodelling, i.e. the resorption and mineralization activity. The infection rapidly leads to severe bone loss, necrosis of the affected portion, and it may even spread to other parts of the body. On the other hand, osteoporosis is not a bacterial infection but similarly is a defective bone pathology arising due to imbalances in the RANK/RANKL/OPG molecular pathway, and due to the progressive weakening of bone structure.

Currently, monomeric fluorescent proteins (FP) are ideal markers for protein tagging. The prediction of oligomeric states is helpful for enhancing live biomedical imaging. Computational prediction of FP oligomeric states can accelerate the effort of protein engineering efforts of creating monomeric FPs. To the best of our knowledge, this study represents the first computational model for predicting and analyzing FP oligomerization directly from the amino acid sequence.

Diagnosing epileptic seizures using electroencephalogram (EEG) in combination with deep learning computational methods has received much attention in recent years. However, to date, deep learning techniques in seizure detection have not been effectively harnessed due to sub-optimal classifier design and improper representation of the time-domain signal.

Accurate prediction of the disease severity of patients with COVID-19 would greatly improve care delivery and resource allocation and thereby reduce mortality risks, especially in less developed countries. Many patient-related factors, such as pre-existing comorbidities, affect disease severity and can be used to aid this prediction.