Efficient Identification of Anti-SARS-CoV-2 Compounds Using Chemical Structure- and Biological Activity-Based Modeling

Improved QSAR models for PARP-1 inhibition using data balancing, interpretable machine learning, and matched molecular pair analysis

The poly (ADP-ribose) polymerase-1 (PARP-1) enzyme is an important target in the treatment of breast cancer. Currently, treatment options include the drugs Olaparib, Niraparib, Rucaparib, and Talazoparib; however, these drugs can cause severe side effects including hematological toxicity and cardiotoxicity. Although in silico models for the prediction of PARP-1 activity have been developed, the drawbacks of these models include low specificity, a narrow applicability domain, and a lack of interpretability. To address these issues, a comprehensive machine learning (ML)-based quantitative structure-activity relationship (QSAR) approach for the informed prediction of PARP-1 activity is presented. Classification models built using the Synthetic Minority Oversampling Technique (SMOTE) for data balancing gave robust and predictive models based on the K-nearest neighbor algorithm (accuracy 0.86, sensitivity 0.88, specificity 0.80). Regression models were built on structurally congeneric datasets, with the models for the phthalazinone class and fused cyclic compounds giving the best performance. In accordance with the Organization for Economic Cooperation and Development (OECD) guidelines, a mechanistic interpretation is proposed using the Shapley Additive Explanations (SHAP) to identify the important topological features to differentiate between PARP-1 actives and inactives. Moreover, an analysis of the PARP-1 dataset revealed the prevalence of activity cliffs, which possibly negatively impacts the model's predictive performance. Finally, a set of chemical transformation rules were extracted using the matched molecular pair analysis (MMPA) which provided mechanistic insights and can guide medicinal chemists in the design of novel PARP-1 inhibitors.

A low-cost machine learning framework for predicting drug-drug interactions based on fusion of multiple features and a parameter self-tuning strategy

Poly-drug therapy is now recognized as a crucial treatment, and the analysis of drug-drug interactions (DDIs) offers substantial theoretical support and guidance for its implementation. Predicting potential DDIs using intelligent algorithms is an emerging approach in pharmacological research. However, the existing supervised models and deep learning-based techniques still have several limitations. This paper proposes a novel DDI analysis and prediction framework called the Multi-View Semi-supervised Graph-based (MVSG) framework, which provides a comprehensive judgment by integrating multiple DDI features and functions without any time-consuming training process. Unlike conventional approaches, MVSG can search for the most suitable similarity (or distance) measurement among DDI data and construct graph structures for each feature. By employing a parameter self-tuning strategy, MVSG fuses multiple graphs according to the contributions of features' information. The actual anticancer drug data are extracted from the authoritative public database for evaluating the effectiveness of our framework, including 904 drugs, 7730 DDI records and 19 types of drug interactions. Validation results indicate that the prediction is more accurate when multiple features are adopted by our framework. In comparison to conventional machine learning techniques, MVSG can achieve higher performance even with less labeled data and without a training process. Finally, MVSG is employed to narrow down the search for potential valuable combinations.

On the origins of SARS-CoV-2 main protease inhibitors

In order to address the world-wide health challenge caused by the COVID-19 pandemic, the 3CL protease/SARS-CoV-2 main protease (SARS-CoV-2-Mpro) coded by its nsp5 gene became one of the biochemical targets for the design of antiviral drugs. In less than 3 years of research, 4 inhibitors of SARS-CoV-2-Mpro have actually been authorized for COVID-19 treatment (nirmatrelvir, ensitrelvir, leritrelvir and simnotrelvir) and more such as EDP-235, FB-2001 and STI-1558/Olgotrelvir or five undisclosed compounds (CDI-988, ASC11, ALG-097558, QLS1128 and H-10517) are undergoing clinical trials. This review is an attempt to picture this quite unprecedented medicinal chemistry feat and provide insights on how these cysteine protease inhibitors were discovered. Since many series of covalent SARS-CoV-2-Mpro inhibitors owe some of their origins to previous work on other proteases, we first provided a description of various inhibitors of cysteine-bearing human caspase-1 or cathepsin K, as well as inhibitors of serine proteases such as human dipeptidyl peptidase-4 or the hepatitis C protein complex NS3/4A. This is then followed by a description of the results of the approaches adopted (repurposing, structure-based and high throughput screening) to discover coronavirus main protease inhibitors.

Medicinal chemistry strategies towards the development of non-covalent SARS-CoV-2 Mpro inhibitors

The main protease (Mpro) of SARS-CoV-2 is an attractive target in anti-COVID-19 therapy for its high conservation and major role in the virus life cycle. The covalent Mpro inhibitor nirmatrelvir (in combination with ritonavir, a pharmacokinetic enhancer) and the non-covalent inhibitor ensitrelvir have shown efficacy in clinical trials and have been approved for therapeutic use. Effective antiviral drugs are needed to fight the pandemic, while non-covalent Mpro inhibitors could be promising alternatives due to their high selectivity and favorable druggability. Numerous non-covalent Mpro inhibitors with desirable properties have been developed based on available crystal structures of Mpro. In this article, we describe medicinal chemistry strategies applied for the discovery and optimization of non-covalent Mpro inhibitors, followed by a general overview and critical analysis of the available information. Prospective viewpoints and insights into current strategies for the development of non-covalent Mpro inhibitors are also discussed.

Selective, broad-spectrum antiviral photodynamic disinfection with dicationic imidazolyl chlorin photosensitizers

The COVID-19 pandemic exposes our vulnerability to viruses that acquire the ability to infect our cells. Classical disinfection methods are limited by toxicity. Existing medicines performed poorly against SARS-CoV-2 because of their specificity to targets in different organisms. We address the challenge of mitigating known and prospective viral infections with a new photosensitizer for antimicrobial photodynamic therapy (aPDT). Photodynamic inactivation is based on local oxidative stress, which is particularly damaging to enveloped viruses. We synthesized a cationic imidazolyl chlorin that reduced by > 99.999% of the percentage inhibition of amplification of SARS-CoV-2 collected from patients at 0.2 µM concentration and 4 J cm-2. Similar results were obtained in the prevention of infection of human ACE2-expressing HEK293T cells by a pseudotyped lentiviral vector exhibiting the S protein of SARS-CoV-2 at its surface. No toxicity to human epidermal keratinocytes (HaCaT) cells was found under similar conditions. aPDT with this chlorin offers fast and safe broad-spectrum photodisinfection and can be repeated with low risk of resistance.

Unravelling the Mystery of COVID-19 Pathogenesis: Spike Protein and Cu Can Synergize to Trigger ROS Production

The COVID-19 pandemic has had a devastating impact on global health, highlighting the need to understand how the SARS-CoV-2 virus damages the lungs in order to develop effective treatments. Recent research has shown that patients with COVID-19 experience severe oxidative damage to various biomolecules. We propose that the overproduction of reactive oxygen species (ROS) in SARS-CoV-2 infection involves an interaction between copper ions and the virus's spike protein. We tested two peptide fragments, Ac-ELDKYFKNH-NH2 (L1) and Ac-WSHPQFEK-NH2 (L2), derived from the spike protein of the Wuhan strain and the β variant, respectively, and found that they bind Cu(II) ions and form a three-nitrogen complexes at lung pH. Our research demonstrates that these complexes trigger the overproduction of ROS, which can break both DNA strands and transform DNA into its linear form. Using A549 cells, we demonstrated that ROS overproduction occurs in the mitochondria, not in the cytoplasm. Our findings highlight the importance of the interaction between copper ions and the virus's spike protein in the development of lung damage and may aid in the development of therapeutic procedures.

Bio-activity prediction of drug candidate compounds targeting SARS-Cov-2 using machine learning approaches

The SARS-CoV-2 3CLpro protein is one of the key therapeutic targets of interest for COVID-19 due to its critical role in viral replication, various high-quality protein crystal structures, and as a basis for computationally screening for compounds with improved inhibitory activity, bioavailability, and ADMETox properties. The ChEMBL and PubChem database contains experimental data from screening small molecules against SARS-CoV-2 3CLpro, which expands the opportunity to learn the pattern and design a computational model that can predict the potency of any drug compound against coronavirus before in-vitro and in-vivo testing. In this study, Utilizing several descriptors, we evaluated 27 machine learning classifiers. We also developed a neural network model that can correctly identify bioactive and inactive chemicals with 91% accuracy, on CheMBL data and 93% accuracy on combined data on both CheMBL and Pubchem. The F1-score for inactive and active compounds was 93% and 94%, respectively. SHAP (SHapley Additive exPlanations) on XGB classifier to find important fingerprints from the PaDEL descriptors for this task. The results indicated that the PaDEL descriptors were effective in predicting bioactivity, the proposed neural network design was efficient, and the Explanatory factor through SHAP correctly identified the important fingertips. In addition, we validated the effectiveness of our proposed model using a large dataset encompassing over 100,000 molecules. This research employed various molecular descriptors to discover the optimal one for this task. To evaluate the effectiveness of these possible medications against SARS-CoV-2, more in-vitro and in-vivo research is required.

Identification of Potent and Selective Acetylcholinesterase/Butyrylcholinesterase Inhibitors by Virtual Screening

Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) play important roles in human neurodegenerative disorders such as Alzheimer's disease. In this study, machine learning methods were applied to develop quantitative structure-activity relationship models for the prediction of novel AChE and BChE inhibitors based on data from quantitative high-throughput screening assays. The models were used to virtually screen an in-house collection of ∼360K compounds. The optimal models achieved good performance with area under the receiver operating characteristic curve values ranging from 0.83 ± 0.03 to 0.87 ± 0.01 for the prediction of AChE/BChE inhibition activity and selectivity. Experimental validation showed that the best-performing models increased the assay hit rate by several folds. We identified 88 novel AChE and 126 novel BChE inhibitors, 25% (AChE) and 53% (BChE) of which showed potent inhibitory effects (IC50 < 5 μM). In addition, structure-activity relationship analysis of the BChE inhibitors revealed scaffolds for chemistry design and optimization. In conclusion, machine learning models were shown to efficiently identify potent and selective inhibitors against AChE and BChE and novel structural series for further design and development of potential therapeutics against neurodegenerative disorders.

Predictive Models for Human Cytochrome P450 3A7 Selective Inhibitors and Substrates

Inappropriate use of prescription drugs is potentially more harmful in fetuses/neonates than in adults. Cytochrome P450 (CYP) 3A subfamily undergoes developmental changes in expression, such as a transition from CYP3A7 to CYP3A4 shortly after birth, which provides a potential way to distinguish medication effects on fetuses/neonates and adults. The purpose of this study was to build first-in-class predictive models for both inhibitors and substrates of CYP3A7/CYP3A4 using chemical structure analysis. Three metrics were used to evaluate model performance: area under the receiver operating characteristic curve (AUC-ROC), balanced accuracy (BA), and Matthews correlation coefficient (MCC). The performance varied for each CYP3A7/CYP3A4 inhibitor/substrate model depending on the data set type, model type, rebalancing method, and specific feature set. For the active inhibitor/substrate data set, the optimal models achieved AUC-ROC values ranging from 0.77 ± 0.01 to 0.84 ± 0.01. For the selective inhibitor/substrate data set, the optimal models achieved AUC-ROC values ranging from 0.72 ± 0.02 to 0.79 ± 0.04. The predictive power of the optimal models was validated by compounds with known potencies as CYP3A7/CYP3A4 inhibitors or substrates. In addition, we identified structural features significant for CYP3A7/CYP3A4 selective or common inhibitors and substrates. In summary, the top performing models can be further applied as a tool to rapidly evaluate the safety and efficacy of new drugs separately for fetuses/neonates and adults. The significant structural features could guide the design of new therapeutic drugs as well as aid in the optimization of existing medicine for fetuses/neonates.

In Silico Identification of Anti-SARS-CoV-2 Medicinal Plants Using Cheminformatics and Machine Learning

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative pathogen of COVID-19, is spreading rapidly and has caused hundreds of millions of infections and millions of deaths worldwide. Due to the lack of specific vaccines and effective treatments for COVID-19, there is an urgent need to identify effective drugs. Traditional Chinese medicine (TCM) is a valuable resource for identifying novel anti-SARS-CoV-2 drugs based on the important contribution of TCM and its potential benefits in COVID-19 treatment. Herein, we aimed to discover novel anti-SARS-CoV-2 compounds and medicinal plants from TCM by establishing a prediction method of anti-SARS-CoV-2 activity using machine learning methods. We first constructed a benchmark dataset from anti-SARS-CoV-2 bioactivity data collected from the ChEMBL database. Then, we established random forest (RF) and support vector machine (SVM) models that both achieved satisfactory predictive performance with AUC values of 0.90. By using this method, a total of 1011 active anti-SARS-CoV-2 compounds were predicted from the TCMSP database. Among these compounds, six compounds with highly potent activity were confirmed in the anti-SARS-CoV-2 experiments. The molecular fingerprint similarity analysis revealed that only 24 of the 1011 compounds have high similarity to the FDA-approved antiviral drugs, indicating that most of the compounds were structurally novel. Based on the predicted anti-SARS-CoV-2 compounds, we identified 74 anti-SARS-CoV-2 medicinal plants through enrichment analysis. The 74 plants are widely distributed in 68 genera and 43 families, 14 of which belong to antipyretic detoxicate plants. In summary, this study provided several medicinal plants with potential anti-SARS-CoV-2 activity, which offer an attractive starting point and a broader scope to mine for potentially novel anti-SARS-CoV-2 drugs.

Bridging the Gap between Target-Based and Cell-Based Drug Discovery with a Graph Generative Multitask Model

The development of new drugs is crucial for protecting humans from disease. In the past several decades, target-based screening has been one of the most popular methods for developing new drugs. This method efficiently screens potential inhibitors of a target protein in vitro, but it frequently fails in vivo due to insufficient activity of the selected drugs. There is a need for accurate computational methods to bridge this gap. Here, we present a novel graph multi-task deep learning model to identify compounds with both target inhibitory and cell active (MATIC) properties. On a carefully curated SARS-CoV-2 data set, the proposed MATIC model shows advantages compared with the traditional method in screening effective compounds in vivo. Following this, we investigated the interpretability of the model and discovered that the learned features for target inhibition (in vitro) or cell active (in vivo) tasks are different with molecular property correlations and atom functional attention. Based on these findings, we utilized a Monte Carlo-based reinforcement learning generative model to generate novel multiproperty compounds with both in vitro and in vivo efficacy, thus bridging the gap between target-based and cell-based drug discovery. The tool is freely accessible at https://github.com/SIAT-code/MATIC.

TranGRU: focusing on both the local and global information of molecules for molecular property prediction

Molecular property prediction is an essential but challenging task in drug discovery. The recurrent neural network (RNN) and Transformer are the mainstream methods for sequence modeling, and both have been successfully applied independently for molecular property prediction. As the local information and global information of molecules are very important for molecular properties, we aim to integrate the bi-directional gated recurrent unit (BiGRU) into the original Transformer encoder, together with self-attention to better capture local and global molecular information simultaneously. To this end, we propose the TranGRU approach, which encodes the local and global information of molecules by using the BiGRU and self-attention, respectively. Then, we use a gate mechanism to reasonably fuse the two molecular representations. In this way, we enhance the ability of the proposed model to encode both local and global molecular information. Compared to the baselines and state-of-the-art methods when treating each task as a single-task classification on Tox21, the proposed approach outperforms the baselines on 9 out of 12 tasks and state-of-the-art methods on 5 out of 12 tasks. TranGRU also obtains the best ROC-AUC scores on BBBP, FDA, LogP, and Tox21 (multitask classification) and has a comparable performance on ToxCast, BACE, and ecoli. On the whole, TranGRU achieves better performance for molecular property prediction. The source code is available in GitHub: https://github.com/Jiangjing0122/TranGRU.

Beyond the vaccines: a glance at the small molecule and peptide-based anti-COVID19 arsenal

Unprecedented efforts of the researchers have been witnessed in the recent past towards the development of vaccine platforms for the control of the COVID-19 pandemic. Albeit, vaccination stands as a practical strategy to prevent SARS-CoV-2 infection, supplementing the anti-COVID19 arsenal with therapeutic options such as small molecules/peptides and antibodies is being conceived as a prudent strategy to tackle the emerging SARS-CoV-2 variants. Noteworthy to mention that collective efforts from numerous teams have led to the generation of a voluminous library composed of chemically and mechanistically diverse small molecules as anti-COVID19 scaffolds. This review article presents an overview of medicinal chemistry campaigns and drug repurposing programs that culminated in the identification of a plethora of small molecule-based anti-COVID19 drugs mediating their antiviral effects through inhibition of proteases, S protein, RdRp, ACE2, TMPRSS2, cathepsin and other targets. In light of the evidence ascertaining the potential of small molecule drugs to approach conserved proteins required for the viral replication of all coronaviruses, accelerated FDA approvals are anticipated for small molecules for the treatment of COVID19 shortly. Though the recent attempts invested in this direction in pursuit of enrichment of the anti-COVID-19 armoury (chemical tools) are praiseworthy, some strategies need to be implemented to extract conclusive benefits of the recently reported small molecule viz. (i) detailed preclinical investigation of the generated anti-COVID19 scaffolds (ii) in-vitro profiling of the inhibitors against the emerging SARS-CoV-2 variants (iii) development of assays enabling rapid screening of the libraries of anti-COVID19 scaffold (iv) leveraging the applications of machine learning based predictive models to expedite the anti-COVID19 drug discovery campaign (v) design of antibody-drug conjugates.