Next Article in Journal
Periostin Augments Vascular Smooth Muscle Cell Calcification via β-Catenin Signaling
Next Article in Special Issue
Identification of Novel Ribonucleotide Reductase Inhibitors for Therapeutic Application in Bile Tract Cancer: An Advanced Pharmacoinformatics Study
Previous Article in Journal
The Beneficial Effects of Ultramicronized Palmitoylethanolamide in the Management of Neuropathic Pain and Associated Mood Disorders Induced by Paclitaxel in Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Potential Commercially Available Inhibitors against SARS-CoV-2 by Multi-Task Deep Learning Model

Guangdong-Hong Kong-Macao Joint Laboratory of Human-Machine Intelligence-Synergy Systems, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Biomolecules 2022, 12(8), 1156; https://doi.org/10.3390/biom12081156
Submission received: 20 July 2022 / Revised: 16 August 2022 / Accepted: 18 August 2022 / Published: 21 August 2022

Abstract

:
The outbreak of COVID-19 caused millions of deaths worldwide, and the number of total infections is still rising. It is necessary to identify some potentially effective drugs that can be used to prevent the development of severe symptoms, or even death for those infected. Fortunately, many efforts have been made and several effective drugs have been identified. The rapidly increasing amount of data is of great help for training an effective and specific deep learning model. In this study, we propose a multi-task deep learning model for the purpose of screening commercially available and effective inhibitors against SARS-CoV-2. First, we pretrained a model on several heterogenous protein–ligand interaction datasets. The model achieved competitive results on some benchmark datasets. Next, a coronavirus-specific dataset was collected and used to fine-tune the model. Then, the fine-tuned model was used to select commercially available drugs against SARS-CoV-2 protein targets. Overall, twenty compounds were listed as potential inhibitors. We further explored the model interpretability and exhibited the predicted important binding sites. Based on this prediction, molecular docking was also performed to visualize the binding modes of the selected inhibitors.

1. Introduction

As of June 2022, SARS-CoV-2 has sickened more than five hundred million and killed over six million people across the globe. SARS-CoV-2 is the seventh member of the family of coronaviruses that infect humans [1,2]. Similar to MERS-CoV and SARS-CoV, SARS-CoV-2 causes severe respiratory diseases and is capable of spreading from person to person. Unfortunately, the rapid mutation of this virus makes it difficult to develop effective vaccines, especially for Omicron. The COVID-19 pandemic will last for years. Therefore, the continuous development of effective anti-SARS-CoV-2 drugs is necessary to prevent the worsening of symptoms, or even death. As a positive-sense, single-stranded RNA beta-coronavirus, SARS-CoV-2 encodes structural, non-structural and accessory proteins. Among these, RNA-dependent RNA polymerase (RdRp), 3-chymotrypsin-like protease (3CLpro), papain-like protease (PLpro), helicase and the spike glycoprotein are supposed to be the main targets. Several compounds that targeted these viral proteins and inhibited coronavirus in vitro have been reported and moved into clinical trials [3]. For example, remdesivir is an approved HIV reverse-transcriptase inhibitor, which has broad-spectrum activities against RNA viruses such as MERS-CoV and SARS-CoV. However, it showed less effective activity in an Ebola clinical trial [4,5,6]. Previous reports showed that remdesivir inhibited SARS-CoV-2 in vitro with an EC50 of 0.77 μM [7], and was used to treat a SARS-CoV-2-infected patient in the United States [8]. However, more potential inhibitors against SARS-CoV-2 are still needed.
Computational methods, such as molecular docking, select compounds that can bind to a target protein and thus improve the success rates of drug discovery. Recently, methods based on deep learning have gained impressive performance in terms of protein–ligand binding prediction [9,10,11]. One main advantage of this algorithm is that it can extract hidden features automatically from raw data and thus significantly improve the prediction accuracy. However, deep learning models might suffer from the generalizability issue due to lack of data. Recently, pretraining models on a large-scale dataset before applying to a small dataset has emerged as a powerful paradigm for solving this issue [12].
In this study, we propose a multi-task deep learning model for selecting potential SARS-CoV-2 inhibitors. First, the model was pretrained on several heterogenous protein–ligand interactions datasets (Figure 1). The model achieved competitive results on several protein–ligand benchmark datasets. Next, a coronavirus-specific dataset was collected and used to fine-tune the model. The fine-tuned model was then used to select commercially available drugs against the SARS-CoV-2 targets. Overall, twenty drugs were listed as potential inhibitors. Furthermore, we explored the model interpretability and found that the predicted important binding sites were close to ground truth.

2. Results and Discussion

2.1. Model Performance on Benchmarks

We first trained a multi-task deep learning model on heterogenous protein–ligand interaction datasets. The model achieved competitive results on several benchmark sets. For example, the PDBbind v.2016 dataset is a widely used benchmark for evaluating protein–ligand interaction prediction methods [15]. It provides more than ten thousand protein–ligand structural complexes with binding affinity (e.g., Kd, Ki). We split the PDBbind set the same way as Pafnucy [16], a classic deep learning-based method for predicting protein–ligand affinity. Moreover, two independent test sets CASF-2013 [17] and Astex Diverse [18] were used to test the generalizability of the model.
Two metrics, RMSE (root mean square error) and Pearson′s correlation coefficient R, were used to evaluate the regression task. As shown in Figure 2, the model achieved RMSE = 1.538, R = 0.71 on the test set and performed well on both independent tests. These results are better than most classic structure-based methods. It is worth noting that our model is only sequence-based.
For the classification task, the model was evaluated on DUD-E, Human, C. elegans and KIBA datasets using 3-fold cross validation. Two traditional machine learning algorithms, SVM and random forest, were also used to compare these datasets. As shown in Figure 3, the proposed model showed excellent performance on most evaluation metrics, including AUC, accuracy, precision, recall, F1-score and specificity. Specifically, the single-task method, which indicates training on one specific label, achieved better performance than the multi-task method on more specific metrics such accuracy and precision, while the multi-task method showed better performance on recall, indicating that it covered more protein–ligand interaction data space. It should be noted that each single-task model corresponded to only one dataset (i.e., it performed well on one dataset but performed poorly on the other datasets). Therefore, the single-task model is likely to perform poorly on coronavirus-specific datasets due to the generalizability issue, even though it performed well on the benchmark dataset. In contrast, the multi-task model was trained once and achieved excellent results on all these datasets. These results suggest that the multi-task model has a better applicability by leveraging heterogenous datasets.

2.2. Screening of SARS-CoV-2 Inhibitors

Recently, pretraining a model on a large-scale dataset before applying it to a small dataset has emerged as a powerful paradigm for improving model generalizability. Inspired by this idea, we fine-tuned the pretrained model on a collected SARS-CoV-2 dataset and then selected potential SARS-CoV-2 inhibitors from a bioactive commercial library containing 10 thousand compounds. After excluding drugs with possible side effects, a total of 10 drugs with high binding affinity were identified, as listed in Table 1. Among these, abacavir (sulfate), a powerful nucleoside analog reverse-transcriptase inhibitor used to treat HIV (human immunodeficiency virus) was predicted to have high binding affinity with multiple proteins of SARS-CoV-2 including RdRp and helicase. Darunavir, a protease inhibitor used to treat HIV, was used in a clinical trial against COVID-19 (ChiCTR2000029541). It should be noted that both darunavir and darunavir (ethanolate) were not present in the training and fine-tuning sets. That is, they were “unseen” by the model. Our model also predicted that darunavir could target 3CLpro and PLpro with affinity Kd = 57.30 and 46.16 nM, respectively, while darunavir (ethanolate) binds to 3CLpro and PLpro with affinity Kd = 44.51 and 35.86 nM, respectively. These results partially prove the accuracy and generalizability of our model. In our predictions, almitrine mesylate, which is a respiratory stimulant that enhances respiration, was used in the treatment of chronic obstructive pulmonary disease. Roflumilast has anti-inflammatory effects and is used as an orally administered drug for the treatment of inflammatory conditions of the lungs such as chronic obstructive pulmonary disease. These two predicted drugs are associated with respiratory symptoms that are the main clinical symptoms of COVID-19. Interestingly, kesuting syrup and keqing capsules were used in a trial for the treatment of mild and moderate COVID-19 (ChiCTR2000029991). It is uncertain whether these drugs only help to alleviate clinical symptoms or have a direct effect on the virus. Daclatasvir is used against Hepatitis C Virus (HCV), which stops HCV viral RNA replication and protein translation by directly inhibiting HCV protein NS5A. In this study, the predicted binding affinity between daclatasvir and RdRp was 15.03 nM. Fiboflapon sodium, a high-affinity 5-lipoxygenase-activating protein inhibitor used for the treatment of asthma, was predicted to have potential affinity to PLpro with Kd = 197.63 nM.
Furthermore, the model was applied to screen lead compounds for 3CLpro inhibition from a diverse screening library containing two million compounds. This screening library, namely, the Enamine HTS Collection, encompasses versatile chemotypes developed within a couple of decades of chemical research at Enamine and its partner academic organizations. Compounds within this library frequently have unusual structures and unique properties. 3CLpro was taken as the screening target because of its wide application and large amount of corresponding data. Thus, the fine-tuned deep learning model had relatively higher accuracy for selecting 3CLpro inhibitors. For the screening, the candidate compounds were first screened by our model. Those candidates were ranked according to the predicted scores. Then, candidates with high predicted scores were further screened based on their possible side effects. For example, Panobinostat, which is a potent and orally active non-selective HDAC inhibitor and has antineoplastic activity, was predicted to inhibit RdRp at 111.19 nM. However, Panobinostat has been reported to possibly cause side effects such as severe diarrhea, headache and even severe infection. Thus, these candidates with severe side effects were excluded. Finally, a total of 10 potential lead compounds are listed in Table 2.

2.3. Model Interpretation

To explore how the model discerns protein–ligand interaction, we conducted a method to identify the key amino acids that are critical for binding. The listed potential inhibitors were regarded as positive samples with high prediction scores. Then, we masked sub-sequences of samples to obtain a “masked” prediction score, and then calculated the importance of these masked sub-sequences. As shown in Figure 4, the critical parts for binding in protein sequences were visualized using a heatmap. In the heatmap, the brighter the region, the more important it is. For 3CLpro, the important amino acids for binding are mostly located at two main parts. It should be noted that different drugs result in different weights in these two regions. For example, roflumilast has higher weights in the first region, indicating the binding sites for roflumilast are close to the middle pocket of 3CLpro. As for abacavir (sulfate), both regions more or less affect the binding, especially the second region at the 180–200th amino acids. For PLpro, the predicted binding sites are located at the 100–120th amino acids.
Moreover, we visualized the predicted binding sites in 3D structures to show possible pockets. The 3CLpro of SARS-CoV-2 (PDB: 6lu7) was used. As mentioned above, two regions of 3CLpro contributed mostly to the binding. As shown in Figure 4A, the region in the upper part (100–200th amino acids) was predicted as the main pocket due to the high weights in most predictions. SARS-CoV-2 PLpro has a catalytic triad composed of Cys114-His275-Asp289 and a conventional zinc-binding domain of four cysteine residues: Cys192, Cys194, Cys227, and Cys229 [19]. The central “thumb” domain contributes the catalytic Cys112 to the active site. Similarly, our model predicted that the 100–120th amino acids contributed mainly to the final binding of small molecules.

2.4. Molecular Docking

Based on our prediction and known protein–ligand structural complexes, we performed the molecular docking of predicted inhibitors to 3CLpro. Four inhibitors, Z56899184, Z57728899, Z57013003 and Z1245218850, were docked to 3CLpro using ArgusLab [20]. For each compound, the best pose with the lowest binding affinity from the ArgusLab docking results was selected. As shown in Figure 5, these inhibitors all bind to amino acids that are located in the active site of 3CLpro (the predicted upper binding region in Figure 4A). Previous studies have shown that the catalytic residue Cys145 is critical for forming covalent bonds between 3CLpro and inhibitors [21,22]. As shown, compound Z56899184 also forms a covalent bond with Cys145 of 3CLpro in the docking result (with binding affinity: −8.67 kcal/mol). These results indicate that Cys145 plays a critical role in the binding process. It should be noted that the binding affinity between the protein and ligand cannot be simply considered based on docking scores (i.e., often expressed in kcal/mol).

3. Conclusions

In this work, a multi-task deep learning model was proposed to predict potential commercial drugs against SARS-CoV-2. There are two main points that should be restated. First, the model was pretrained on heterogenous protein–ligand interaction datasets before being applied to the SARS-CoV-2 inhibitor screening. This idea was inspired by the recent success of pretrained models on various fields such as natural language processing. The core logic of this method was transfer learning, which could extract features from large-scale datasets and then transfer the learned knowledge to a small set. This method has a broad application in drug discovery which is often hindered by the small data issue. Although the pretrained model had wide distribution in the protein–ligand interaction data space, the distance between existing knowledge and future application is ambiguous. Therefore, which part of the pretraining set would cause issues such as negative effect was not understood well. Second, the exploration of the “black box” within deep learning model is very important, especially in the drug discovery area. The black box model may increase the risk of subsequentially false leads, thus causing severe consequences. Our model could provide some insights into this. Overall, we are on a cheerful avenue for the application of pretrained deep learning models in drug discovery.

4. Methods

4.1. Data

A few benchmark protein–ligand interaction datasets were used for model pretraining, including DUD-E [23], Human, C. elegans [24], KIBA [25] and PDBbind [15]. For classification, DUD-E, Human, C. elegans and kinase inhibitor bioactivity (KIBA) were used. A total of 91,220 samples from the DUD-E dataset were used in this study, in which the negative-to-positive ratio was set to 3:1 to avoid an imbalance of data. Similarly, the negative-to-positive ratios of Human and C. elegans sets were also set to 3:1. For the KIBA, protein–ligand interactions with values bigger than 12.1 were regarded as positive samples. For regression, the PDBbind v.2016 and kinase dataset Davis were used. Specifically, 13,196 protein–ligand complexes from the PDBbind set were split the same way as Pafnucy. Similarly, the kinase dataset Davis consisting of a total of 9125 interactions with the corresponding binding affinity (Kd) was used here.
For drug screening, several important SARS-CoV-2 proteins were used as targets, such as RNA-dependent RNA polymerase (RdRp), 3-chymotrypsin-like protease, papain-like protease and helicase. The sequences of these proteins were extracted from NCBI (NC_045512.2). The SARS-CoV/SARS-CoV-2 specific dataset was collected from various papers [21,26,27,28,29] and public datasets such as GHDDI and PubChem with filtration. After removing duplicates, a collection of 10 thousand approved and bioactive compounds was used for screening. Similarly, a diverse Enamine HTS collection containing 2 million compounds was also used.

4.2. Model

Basically, the proposed model consists of three parts: protein feature extraction by word2vec and transformer, drug feature extraction by node2vec and interaction prediction by multi-task modules. The last dense layer was activated as output. The loss functions were defined as binary cross entropy (BCE) and mean squared error (MSE) for classification and regression, respectively. Figure 1 displays the basic architecture of the model, and more details are described below.
Preprocessing. Given the drug SMILES and protein sequence, we first constructed a graph and words, respectively, and then obtained the dense vectors by the embedding method. Let L i d x = x 1 ,   x 2 , ,   x L 1 ,   x L be a word sequence of each protein, where idx is the sequence index, L is the word sequence length and x i R d 1 is the d1-dimensional embedding of the ith word that we split by overlapping 3-gram on the original sequence. For drugs, let G = ( V ,   E ) be a given drug SMILES and f : V R d 2 be the mapping function from V i d x = { v 1 ,   v 2 ,   ,   v l 1 ,   v l } to feature representation. Where i d x is the sequence index, l is the drug length and v i is the d2-dimensional embedding of the i th node on the original sequence [30].
Pretraining. With the building of a large-scale word-segmented corpus of protein sequence, we then used each current word as an input to a log-linear classifier with a continuous projection layer, and predicted words within a certain range before and after the current word, proposed as a skip-gram model [13], to obtain the final d 1 -dimensional embedding words. Meanwhile, to identify the specific drug network nodes, we considered node2vec [14] as the representation learning method to design a flexible neighborhood sampling strategy to explore a random walk, such as { s ,   v 3 ,   v 5 ,   v 9 ,   v 8 ,   v 1 ,   v 6 ,   v 4 }   of start node s = v 2 and walking length l = 8; thus, we can consider the skip-gram architecture mentioned above as a feature learning method.
Multi-task module. This module consists of two parts: shared layers and task-specific layers. The shared layers are designed to learn a joint representation for all tasks. Task-specific layers are used to learn the weights of specific blocks based on the joint representation. In this study, two related tasks are defined: classification and regression. The model was fine-tuned by a virus-specific dataset to acquire robust results for coronavirus.
The model weights were initialized using a Glorot uniform initializer, which keeps the variance of output and input by correcting the variance of uniform. Then, the Adam optimize was used to train the model. The initial learning rate and batch size was set to 10-4 and 256. During the training process, early stopping was used to monitor the training and validation loss in order to select the final model without overfitting.

4.3. Biological Interpretation

A non-parametric method, “occlusion”, used in our previous study [11], was applied to explore which parts of the input protein sequences were critical to the task. The intuitive explanation of this method is: (i) mask a subsequence (e.g., 10–20th amino acids) of a full-length protein sequence; (ii) predict a result of this modified sequence; (iii) calculate the difference between this result and the ground truth; (iv) repeat (i) to (iii) and find the most important amino acids. Specifically, si from test samples (i = 0, 1, 2,..., n−1, here, n is sample size of test set) was expressed as tuple (protein input_i, compound input_i). While maintaining compound input_i as unchanged, we systematically masked the protein input_i in si to track the changes of the output. Then, the importance of each sub-sequence in the sequence to the prediction could be calculated.

Author Contributions

Conceptualization, F.H. and P.Y.; methodology, F.H.; software, F.H. and J.J.; validation, F.H., J.J. and P.Y.; formal analysis, F.H.; investigation, F.H.; resources, F.H.; data curation, F.H.; writing—original draft preparation, F.H.; writing—review and editing, F.H. and P.Y.; visualization, F.H. and J.J.; supervision, P.Y.; project administration, F.H. and P.Y.; funding acquisition, P.Y.. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Strategic Priority Research Program of Chinese Academy of Sciences (NO. XDB 38040200), the Shenzhen Science and Technology Innovation Committee (JCYJ20180703145002040) and Shenzhen Science and Technology Program (Grant No. SGDX20201103095603009).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The source code and data are available at https://github.com/AI-bio/multi-task-for-cov, accessed on 17 August 2022.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef] [PubMed]
  2. Gorbalenya, A.E.; Baker, S.C.; Baric, R.S.; de Groot, R.J.; Drosten, C.; Gulyaeva, A.A.; Haagmans, B.L.; Lauber, C.; Leontovich, A.M.; Neuman, B.W.; et al. Severe Acute Respiratory Syndrome-Related Coronavirus: The Species and Its Viruses-a Statement of the Coronavirus Study Group. bioRxiv 2020, 1, 1–15. [Google Scholar] [CrossRef] [Green Version]
  3. Li, G.; De Clercq, E. Therapeutic Options for the 2019 Novel Coronavirus (2019-NCoV). Nat. Rev. Drug Discov. 2020, 19, 149–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Brown, A.J.; Won, J.J.; Graham, R.L.; Dinnon, K.H.; Sims, A.C.; Feng, J.Y.; Cihlar, T.; Denison, M.R.; Baric, R.S.; Sheahan, T.P. Broad Spectrum Antiviral Remdesivir Inhibits Human Endemic and Zoonotic Deltacoronaviruses with a Highly Divergent RNA Dependent RNA Polymerase. Antivir. Res. 2019, 169, 104541. [Google Scholar] [CrossRef]
  5. Sheahan, T.P.; Sims, A.C.; Leist, S.R.; Schäfer, A.; Won, J.; Brown, A.J.; Montgomery, S.A.; Hogg, A.; Babusis, D.; Clarke, M.O.; et al. Comparative Therapeutic Efficacy of Remdesivir and Combination Lopinavir, Ritonavir, and Interferon Beta against MERS-CoV. Nat. Commun. 2020, 11, 222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Tchesnokov, E.; Feng, J.; Porter, D.; Götte, M. Mechanism of Inhibition of Ebola Virus RNA-Dependent RNA Polymerase by Remdesivir. Viruses 2019, 11, 326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Wang, M.; Cao, R.; Zhang, L.; Yang, X.; Liu, J.; Xu, M.; Shi, Z.; Hu, Z.; Zhong, W.; Xiao, G. Remdesivir and Chloroquine Effectively Inhibit the Recently Emerged Novel Coronavirus (2019-NCoV) in Vitro. Cell Res. 2020, 30, 269–271. [Google Scholar] [CrossRef]
  8. Holshue, M.L.; DeBolt, C.; Lindquist, S.; Lofy, K.H.; Wiesman, J.; Bruce, H.; Spitters, C.; Ericson, K.; Wilkerson, S.; Tural, A.; et al. First Case of 2019 Novel Coronavirus in the United States. N. Engl. J. Med. 2020, 382, 929–936. [Google Scholar] [CrossRef]
  9. Öztürk, H.; Özgür, A.; Ozkirimli, E. DeepDTA: Deep Drug–Target Binding Affinity Prediction. Bioinformatics 2018, 34, i821–i829. [Google Scholar] [CrossRef] [Green Version]
  10. Tsubaki, M.; Tomii, K.; Sese, J. Compound-Protein Interaction Prediction with End-to-End Learning of Neural Networks for Graphs and Sequences. Bioinformatics 2018, 35, 309–318. [Google Scholar] [CrossRef]
  11. Hu, F.; Jiang, J.; Wang, D.; Zhu, M.; Yin, P. Multi-PLI: Interpretable Multi-task Deep Learning Model for Unifying Protein–Ligand Interaction Datasets. J. Cheminform. 2021, 13, 30. [Google Scholar] [CrossRef] [PubMed]
  12. Devlin, J.; Chang, M.-W.; Lee, K.; Toutanova, K. BERT: Pre-Training of Deep Bidirectional Transformers for Language Understanding. arXiv 2018, arXiv:1810.04805. [Google Scholar]
  13. Mikolov, T.; Chen, K.; Corrado, G.; Dean, J. Efficient Estimation of Word Representations in Vector Space. arXiv 2013, arXiv:1301.3781. [Google Scholar]
  14. Grover, A.; Leskovec, J. Node2vec. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; ACM: New York, NY, USA, 2016; pp. 855–864. [Google Scholar]
  15. Wang, R.; Fang, X.; Lu, Y.; Wang, S. The PDBbind Database: Collection of Binding Affinities for Protein−Ligand Complexes with Known Three-Dimensional Structures. J. Med. Chem. 2004, 47, 2977–2980. [Google Scholar] [CrossRef]
  16. Stepniewska-Dziubinska, M.M.; Zielenkiewicz, P.; Siedlecki, P. Development and Evaluation of a Deep Learning Model for Protein–Ligand Binding Affinity Prediction. Bioinformatics 2018, 34, 3666–3674. [Google Scholar] [CrossRef] [Green Version]
  17. Li, Y.; Han, L.; Liu, Z.; Wang, R. Comparative Assessment of Scoring Functions on an Updated Benchmark: 2. Evaluation Methods and General Results. J. Chem. Inf. Model. 2014, 54, 1717–1736. [Google Scholar] [CrossRef]
  18. Hartshorn, M.J.; Verdonk, M.L.; Chessari, G.; Brewerton, S.C.; Mooij, W.T.M.; Mortenson, P.N.; Murray, C.W. Diverse, High-Quality Test Set for the Validation of Protein-Ligand Docking Performance. J. Med. Chem. 2007, 50, 726–741. [Google Scholar] [CrossRef]
  19. Gao, X.; Qin, B.; Chen, P.; Zhu, K.; Hou, P.; Wojdyla, J.A.; Wang, M.; Cui, S. Crystal Structure of SARS-CoV-2 Papain-like Protease. Acta Pharm. Sin. B 2021, 11, 237–245. [Google Scholar] [CrossRef]
  20. Bitencourt-Ferreira, G.; de Azevedo, W.F. Molecular Docking Simulations with ArgusLab. In Docking Screens for Drug Discovery; Humana: New York, NY, USA, 2019; pp. 203–220. [Google Scholar]
  21. Zhu, W.; Xu, M.; Chen, C.Z.; Guo, H.; Shen, M.; Hu, X.; Shinn, P.; Klumpp-Thomas, C.; Michael, S.G.; Zheng, W. Identification of SARS-CoV-2 3CL Protease Inhibitors by a Quantitative High-Throughput Screening. ACS Pharmacol. Transl. Sci. 2020, 3, 1008–1016. [Google Scholar] [CrossRef]
  22. Hu, F.; Wang, L.; Hu, Y.; Wang, D.; Wang, W.; Jiang, J.; Li, N.; Yin, P. A Novel Framework Integrating AI Model and Enzymological Experiments Promotes Identification of SARS-CoV-2 3CL Protease Inhibitors and Activity-Based Probe. Brief. Bioinform. 2021, 22, bbab301. [Google Scholar] [CrossRef]
  23. Mysinger, M.M.; Carchia, M.; Irwin, J.J.; Shoichet, B.K. Directory of Useful Decoys, Enhanced (DUD-E): Better Ligands and Decoys for Better Benchmarking. J. Med. Chem. 2012, 55, 6582–6594. [Google Scholar] [CrossRef] [PubMed]
  24. Liu, H.; Sun, J.; Guan, J.; Zheng, J.; Zhou, S. Improving Compound-Protein Interaction Prediction by Building up Highly Credible Negative Samples. Bioinformatics 2015, 31, i221–i229. [Google Scholar] [CrossRef] [PubMed]
  25. Tang, J.; Szwajda, A.; Shakyawar, S.; Xu, T.; Hintsanen, P.; Wennerberg, K.; Aittokallio, T. Making Sense of Large-Scale Kinase Inhibitor Bioactivity Data Sets: A Comparative and Integrative Analysis. J. Chem. Inf. Model. 2014, 54, 735–743. [Google Scholar] [CrossRef] [PubMed]
  26. Jin, Z.; Du, X.; Xu, Y.; Deng, Y.; Liu, M.; Zhao, Y.; Zhang, B.; Li, X.; Zhang, L.; Peng, C.; et al. Structure of Mpro from SARS-CoV-2 and Discovery of Its Inhibitors. Nature 2020, 582, 289–293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Fu, L.; Ye, F.; Feng, Y.; Yu, F.; Wang, Q.; Wu, Y.; Zhao, C.; Sun, H.; Huang, B.; Niu, P.; et al. Both Boceprevir and GC376 Efficaciously Inhibit SARS-CoV-2 by Targeting Its Main Protease. Nat. Commun. 2020, 11, 4417. [Google Scholar] [CrossRef]
  28. Ma, C.; Sacco, M.D.; Hurst, B.; Townsend, J.A.; Hu, Y.; Szeto, T.; Zhang, X.; Tarbet, B.; Marty, M.T.; Chen, Y.; et al. Boceprevir, GC-376, and Calpain Inhibitors II, XII Inhibit SARS-CoV-2 Viral Replication by Targeting the Viral Main Protease. Cell Res. 2020, 30, 678–692. [Google Scholar] [CrossRef]
  29. Riva, L.; Yuan, S.; Yin, X.; Martin-Sancho, L.; Matsunaga, N.; Pache, L.; Burgstaller-Muehlbacher, S.; De Jesus, P.D.; Teriete, P.; Hull, M.V.; et al. Discovery of SARS-CoV-2 Antiviral Drugs through Large-Scale Compound Repurposing. Nature 2020, 586, 113–119. [Google Scholar] [CrossRef]
  30. Dong, Q.-W.; Wang, X.-L.; Lin, L. Application of Latent Semantic Analysis to Protein Remote Homology Detection. Bioinformatics 2006, 22, 285–290. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of the proposed model. The model involves two parts: protein/ligand feature extraction and their interaction prediction. First, the protein sequence is processed in turn by word2vec [13], multi-heads residual layer and Bi-GRU (bidirectional gated recurrent unit) modules. Ligand smiles is processed by node2vec [14]. Then, their representations are fed into a shared layer and task-specific layer and connected to the labels.
Figure 1. Schematic of the proposed model. The model involves two parts: protein/ligand feature extraction and their interaction prediction. First, the protein sequence is processed in turn by word2vec [13], multi-heads residual layer and Bi-GRU (bidirectional gated recurrent unit) modules. Ligand smiles is processed by node2vec [14]. Then, their representations are fed into a shared layer and task-specific layer and connected to the labels.
Biomolecules 12 01156 g001
Figure 2. Model performance on the PDBbind v.2016 set and two independent sets. (A) PDBbind validation set. (B) PDBbind v.2016 core set. (C) CASF-2013. (D) Astex Diverse set. RMSE: root mean square error. MAE: mean absolute error. R: Pearson′s correlation coefficient R. SD: stand deviation.
Figure 2. Model performance on the PDBbind v.2016 set and two independent sets. (A) PDBbind validation set. (B) PDBbind v.2016 core set. (C) CASF-2013. (D) Astex Diverse set. RMSE: root mean square error. MAE: mean absolute error. R: Pearson′s correlation coefficient R. SD: stand deviation.
Biomolecules 12 01156 g002
Figure 3. Model performance on (A) DUD-E. (B) Human and C. elegans. (C)KIBA datasets.
Figure 3. Model performance on (A) DUD-E. (B) Human and C. elegans. (C)KIBA datasets.
Biomolecules 12 01156 g003
Figure 4. The predicted binding sites. (A) 3CLpro (PDB ID: 6lu7); (B) PLpro (PDB ID: 6wx4). The abscissa axis is the length of protein sequences. Red indicates the predicted important sites and black dotted circles indicate possible pockets.
Figure 4. The predicted binding sites. (A) 3CLpro (PDB ID: 6lu7); (B) PLpro (PDB ID: 6wx4). The abscissa axis is the length of protein sequences. Red indicates the predicted important sites and black dotted circles indicate possible pockets.
Biomolecules 12 01156 g004
Figure 5. Molecular docking of predicted inhibitors to SARS-CoV-2 3CLpro. (A) Z56899184, −8.67 kcal/mol; (B) Z57728899, −7.86 kcal/mol; (C) Z57013003, −7.45 kcal/mol; (D) Z1245218850, −8.41 kcal/mol. 3CLpro is represented as translucent gray cartoon. The involved amino acids are represented in cyan.
Figure 5. Molecular docking of predicted inhibitors to SARS-CoV-2 3CLpro. (A) Z56899184, −8.67 kcal/mol; (B) Z57728899, −7.86 kcal/mol; (C) Z57013003, −7.45 kcal/mol; (D) Z1245218850, −8.41 kcal/mol. 3CLpro is represented as translucent gray cartoon. The involved amino acids are represented in cyan.
Biomolecules 12 01156 g005
Table 1. Potential inhibitors for SARS-CoV-2.
Table 1. Potential inhibitors for SARS-CoV-2.
DrugCASTargetPredicted Affinity (nM)
Abacavir (sulfate)188062-50-2RdRp
helicase
3.03
3.06
Darunavir206361-99-13CLpro
PLpro
57.30
46.16
Darunavir (ethanolate)635728-49-33CLpro
PLpro
44.51
35.86
Itraconazole84625-61-6PLpro
RdRp
127.98
16.90
Almitrine mesylate29608-49-93CLpro29.31
Daclatasvir1009119-64-5RdRp15.03
Daclatasvir (dihydrochloride)1009119-65-6RdRp19.87
Metoprolol tartrate56392-17-7PLpro153.23
Fiboflapon sodium1196070-26-4PLpro197.63
Roflumilast162401-32-33CLpro248.89
Table 2. The predicted lead compounds for SARS-CoV-2 3CLpro.
Table 2. The predicted lead compounds for SARS-CoV-2 3CLpro.
IDSMILESProbability
Z56899184O=[N+]([O-])c1cccc(SSc2cccc([N+](=O)[O-])c2)c10.969
Z229622170N#Cc1cccc(CN2C(=O)C(=O)c3cccc(Br)c32)c10.947
Z57728899COc1ccc(CN2C(=O)C(=O)c3cc(Br)ccc32)cc10.918
Z90667629[Na+].[O-][n+]1ccccc1[S-]0.863
Z1238998507N#Cc1cc([N+](=O)[O-])ccc1Oc1cncc(Cl)c10.782
Z56833036O=C1c2ccccc2C(=O)c2c1cccc2S(=O)(=O)N1CCOCC10.724
Z57013003COc1ccc(N2CCN(C(=O)c3cc(=O)[nH]c4ccccc34)CC2)cc10.711
Z1245218850N#CCc1cccc(C(=O)Oc2cncc(Cl)c2)c10.697
Z56785091c1nc(SSc2nc[nH]n2)n[nH]10.682
Z1776036493Cc1cc(Br)cc2c1N(CCBr)C(=O)C2=O0.676
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hu, F.; Jiang, J.; Yin, P. Prediction of Potential Commercially Available Inhibitors against SARS-CoV-2 by Multi-Task Deep Learning Model. Biomolecules 2022, 12, 1156. https://doi.org/10.3390/biom12081156

AMA Style

Hu F, Jiang J, Yin P. Prediction of Potential Commercially Available Inhibitors against SARS-CoV-2 by Multi-Task Deep Learning Model. Biomolecules. 2022; 12(8):1156. https://doi.org/10.3390/biom12081156

Chicago/Turabian Style

Hu, Fan, Jiaxin Jiang, and Peng Yin. 2022. "Prediction of Potential Commercially Available Inhibitors against SARS-CoV-2 by Multi-Task Deep Learning Model" Biomolecules 12, no. 8: 1156. https://doi.org/10.3390/biom12081156

APA Style

Hu, F., Jiang, J., & Yin, P. (2022). Prediction of Potential Commercially Available Inhibitors against SARS-CoV-2 by Multi-Task Deep Learning Model. Biomolecules, 12(8), 1156. https://doi.org/10.3390/biom12081156

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop