A survey on computational methods in discovering protein inhibitors of SARS-CoV-2

From docking to dynamics: Unveiling the potential non-peptide and non-covalent inhibitors of Mpro from natural products

MOTIVATION

This study aims to investigate non-covalent and non-peptide inhibitors of Mpro, a crucial protein target, by employing a comprehensive approach that integrates molecular docking, molecular dynamics simulations, and top-hits activity predictions. The focus is on elucidating the non-covalent and non-peptide binding modes of potential inhibitors with Mpro.

METHODS

We employed a semi-flexible molecular docking methodology, binding score and ADME screening, which are based on structure, to screen compounds from CMNPD and HERB in silico. These methodologies allowed us to find potential candidates depending on their binding values and interactions with the binding site of main protease. To further evaluate the stability of these interactions, we conducted molecular dynamics simulations and calculated binding energies. Ultimately, a top-hits activity prediction method was employed to prioritize compounds based on their predicted inhibitory potential.

RESULTS

Through a combination of binding energy calculations and activity predictions, we identified six potential inhibitor molecules exhibiting promising activity against Mpro. These compounds demonstrated favorable binding interactions and stability profiles, making them attractive candidates for further experimental validation and drug development efforts targeting Mpro.

Using TransR to enhance drug repurposing knowledge graph for COVID-19 and its complications

MOTIVATION

The COVID-19 pandemic has been spreading globally for four years, yet specific drugs that effectively suppress the virus remain elusive. Furthermore, the emergence of complications associated with COVID-19 presents significant challenges, making the development of therapeutics for COVID-19 and its complications an urgent task. However, traditional drug development processes are time-consuming. Drug repurposing, which involves identifying new therapeutic applications for existing drugs, presents a viable alternative.

RESULT

In this study, we construct a knowledge graph by retrieving information on genes, drugs, and diseases from databases such as DRUGBANK and GNBR. Next, we employ the TransR knowledge representation learning approach to embed entities and relationships into the knowledge graph. Subsequently, we train the knowledge graph using a graph neural network model based on TransR scoring. This trained knowledge graph is then utilized to predict drugs for the treatment of COVID-19 and its complications. Based on experimental results, we have identified 15 drugs out of the top 30 with the highest success rates associated with treating COVID-19 and its complications. Notably, out of these 15 drugs, 10 specifically aimed at treating COVID-19, such as Torcetrapib and Tocopherol, has not been previously identified in the knowledge graph. This finding highlights the potential of our model in aiding healthcare professionals in drug development and research related to this disease.

Predicting Drug-Disease Associations Through Similarity Network Fusion and Multi-View Feature Projection Representation

Predicting drug-disease associations (DDAs) through computational methods has become a prevalent trend in drug development because of their high efficiency and low cost. Existing methods usually focus on constructing heterogeneous networks by collecting multiple data resources to improve prediction ability. However, potential association possibilities of numerous unconfirmed drug-related or disease-related pairs are not sufficiently considered. In this article, we propose a novel computational model to predict new DDAs. First, a heterogeneous network is constructed, including four types of nodes (drugs, targets, cell lines, diseases) and three types of edges (associations, association scores, similarities). Second, an updating and merging-based similarity network fusion method, termed UM-SF, is presented to fuse various similarity networks with diverse weights. Finally, an intermediate layer-mediated multi-view feature projection representation method, termed IM-FP, is proposed to calculate the predicted DDA scores. This method uses multiple association scores to construct multi-view drug features, then projects them into disease space through the intermediate layer, where an intermediate layer similarity constraint is designed to learn the projection matrices. Results of comparative experiments reveal the effectiveness of our innovations. Comparisons with other state-of-the-art models by the 10-fold cross-validation experiment indicate our model's advantage on AUROC and AUPR metrics. Moreover, our proposed model successfully predicted 107 novel high-ranked DDAs.

A First Computational Frame for Recognizing Heparin-Binding Protein

Heparin-binding protein (HBP) is a cationic antibacterial protein derived from multinuclear neutrophils and an important biomarker of infectious diseases. The correct identification of HBP is of great significance to the study of infectious diseases. This work provides the first HBP recognition framework based on machine learning to accurately identify HBP. By using four sequence descriptors, HBP and non-HBP samples were represented by discrete numbers. By inputting these features into a support vector machine (SVM) and random forest (RF) algorithm and comparing the prediction performances of these methods on training data and independent test data, it is found that the SVM-based classifier has the greatest potential to identify HBP. The model could produce an auROC of 0.981 ± 0.028 on training data using 10-fold cross-validation and an overall accuracy of 95.0% on independent test data. As the first model for HBP recognition, it will provide some help for infectious diseases and stimulate further research in related fields.

Optimization of potential non-covalent inhibitors for the SARS-CoV-2 main protease inspected by a descriptor of the subpocket occupancy

The main protease is regarded as an essential drug target for treating Coronavirus Disease 2019. In the present study, 13 marketed drugs were investigated to explore the possible binding mechanism, utilizing molecular docking, molecular dynamics simulation, and MM-PB(GB)SA binding energy calculations. Our results suggest that fusidic acid, polydatin, SEN-1269, AZD6482, and UNC-2327 have high binding affinities of more than 23 kcal mol^-1. A descriptor was defined for the energetic occupancy of the subpocket, and it was found that S4 had a low occupancy of less than 10% on average. The molecular optimization of ADZ6482 via reinforcement learning algorithms was carried out to screen out three lead compounds, in which slight structural changes give more considerable binding energies and an occupancy of the S4 subpocket of up to 43%. The energetic occupancy could be a useful descriptor for evaluating the local binding affinity for drug design.

Analysis and modeling of myopia-related factors based on questionnaire survey

With the rapid development of science and technology, the trend of low age myopia is becoming increasingly significant. The latest national survey done by the Chinese government found that more than 80% of Chinese teenagers suffer from myopia. Adolescent myopia is closely related to living environment, heredity, and living habits. Quantifying the relationship between myopia and living environment, heredity, and living habits is conductive to the prevention and intervention of adolescent myopia. In this study, we investigated the relationships between four main factors (environment, habits, parental vision, and demographic) and myopia status by analyzing the questionnaire data. Data were collected from Chengdu, China in 2021, including 2808 myopia samples and 5693 non-myopia samples, with a total of 22 features. Then, these 22 features were inputted into three machine learning algorithms to discriminate the two classes of samples. Results show that the computational model could produce an AUC of 0.768. To pick out the most important features which play important roles in classification, we used incremental feature selection strategy to screen the 22 features. As a result, we found that the 4 most influential features with XGBoost could achieve a competitive AUC of 0.764. To further investigate the risk and protective factors affecting adolescent myopia, we used OR values derived from MLE-LR to analyze the relationship between 22 features and adolescent myopia. Results showed that the age variable was the most significant risk factor for myopia, followed by the myopia status of parents. The most protective factor for eyesight is the measure taken by the children, followed by the distance between books and eyes when reading. These discoveries can guide the prevention and control of myopia in children and adolescents.

TCM2COVID: A resource of anti-COVID-19 traditional Chinese medicine with effects and mechanisms

In China, traditional Chinese medicine (TCM) has been widely used for coronavirus infectious disease 2019 (COVID-19) prevention, treatment, and recovery and has played a part in the battle against the disease. A variety of TCM treatments have been recommended for different stages of COVID-19. But, to the best of our knowledge, a comprehensive database for storing and organizing anti-COVID TCM treatments is still lacking. Herein, we developed TCM2COVID, a manually curated resource of anti-COVID TCM formulas, natural products (NPs), and herbs. The current version of TCM2COVID (1) documents over 280 TCM formulas (including over 300 herbs) with detailed clinical evidence and therapeutic mechanism information; (2) records over 80 NPs with detailed potential therapeutic mechanisms; and (3) launches a useful web server for querying, analyzing and visualizing documented formulas similar to those supplied by the user (formula similarity analysis). In summary, TCM2COVD provides a user-friendly and practical platform for documenting, querying, and browsing anti-COVID TCM treatments, and will help in the development and elucidation of the mechanisms of action of new anti-COVID TCM therapies to support the fight against the COVID-19 epidemic. TCM2COVID is freely available at http://zhangy-lab.cn/tcm2covid/.

IBPred: a sequence-based predictor for identifying ion binding protein in phage

Ion binding proteins (IBPs) can selectively and non-covalently interact with ions. IBPs in phages also play an important role in biological processes. Therefore, accurate identification of IBPs is necessary for understanding their biological functions and molecular mechanisms that involve binding to ions. Since molecular biology experimental methods are still labor-intensive and cost-ineffective in identifying IBPs, it is helpful to develop computational methods to identify IBPs quickly and efficiently. In this work, a random forest (RF)-based model was constructed to quickly identify IBPs. Based on the protein sequence information and residues’ physicochemical properties, the dipeptide composition combined with the physicochemical correlation between two residues were proposed for the extraction of features. A feature selection technique called analysis of variance (ANOVA) was used to exclude redundant information. By comparing with other classified methods, we demonstrated that our method could identify IBPs accurately. Based on the model, a Python package named IBPred was built with the source code which can be accessed at https://github.com/ShishiYuan/IBPred.

Design of SARS-CoV-2 Main Protease Inhibitors Using Artificial Intelligence and Molecular Dynamic Simulations

Drug design is a time-consuming and cumbersome process due to the vast search space of drug-like molecules and the difficulty of investigating atomic and electronic interactions. The present paper proposes a computational drug design workflow that combines artificial intelligence (AI) methods, i.e., an evolutionary algorithm and artificial neural network model, and molecular dynamics (MD) simulations to design and evaluate potential drug candidates. For the purpose of illustration, the proposed workflow was applied to design drug candidates against the main protease of severe acute respiratory syndrome coronavirus 2. From the ∼140,000 molecules designed using AI methods, MD analysis identified two molecules as potential drug candidates.

Heterogeneity Analysis of Bladder Cancer Based on DNA Methylation Molecular Profiling

Bladder cancer is a highly complex and heterogeneous malignancy. Tumor heterogeneity is a barrier to effective diagnosis and treatment of bladder cancer. Human carcinogenesis is closely related to abnormal gene expression, and DNA methylation is an important regulatory factor of gene expression. Therefore, it is of great significance for bladder cancer research to characterize tumor heterogeneity by integrating genetic and epigenetic characteristics. This study explored specific molecular subtypes based on DNA methylation status and identified subtype-specific characteristics using patient samples from the TCGA database with DNA methylation and gene expression were measured simultaneously. The results were validated using an independent cohort from GEO database. Four DNA methylation molecular subtypes of bladder cancer were obtained with different prognostic states. In addition, subtype-specific DNA methylation markers were identified using an information entropy-based algorithm to represent the unique molecular characteristics of the subtype and verified in the test set. The results of this study can provide an important reference for clinicians to make treatment decisions.

Computational study on the affinity of potential drugs to SARS-CoV-2 main protease

Herein, we report a computational investigation of the binding affinity of dexamethasone, betamethasone, chloroquine and hydroxychloroquine to SARS-CoV-2 main protease using molecular and quantum mechanics as well as molecular docking methodologies. We aim to provide information on the anti-COVID-19 mechanism of the abovementioned potential drugs against SARS-CoV-2 coronavirus. Hence, the 6w63 structure of the SARS-CoV-2 main protease was selected as potential target site for the docking analysis. The study includes an initial conformational analysis of dexamethasone, betamethasone, chloroquine and hydroxychloroquine. For the most stable conformers, a spectroscopic analysis has been carried out. In addition, global and local reactivity indexes have been calculated to predict the chemical reactivity of these molecules. The molecular docking results indicate that dexamethasone and betamethasone have a higher affinity than chloroquine and hydroxychloroquine for their theoretical 6w63 target. Additionally, dexamethasone and betamethasone show a hydrogen bond with the His41 residue of the 6w63 protein, while the interaction between chloroquine and hydroxychloroquine with this amino acid is weak. Thus, we confirm the importance of His41 amino acid as a target to inhibit the SARS-CoV-2 Mpro activity.