Next Article in Journal
Changes in Vaccine Hesitancy in Japan across Five Months during the COVID-19 Pandemic and Its Related Factors
Previous Article in Journal
Effectiveness of COVID-19 Vaccines against Delta (B.1.617.2) Variant: A Systematic Review and Meta-Analysis of Clinical Studies
Previous Article in Special Issue
Protection of Vaccine Preventable Diseases in a Population of HIV-Infected Children: A 3 Years Prospective Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Targeted Approaches and Drug Repurposing Reveal Possible SARS-CoV-2 Inhibitors

by
Khalid Mashay Alanazi
1,
Mohammad Abul Farah
1 and
Yan-Yan Hor
2,*
1
Department of Zoology, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
2
Department of Biotechnology, Yeungnam University, 280 Daehak-ro, Gyeongsan 38541, Gyeongbuk-do, Korea
*
Author to whom correspondence should be addressed.
Submission received: 13 November 2021 / Revised: 3 December 2021 / Accepted: 16 December 2021 / Published: 26 December 2021

Abstract

:
The COVID-19 pandemic caused by SARS-CoV-2 is unprecedented in recent memory owing to the non-stop escalation in number of infections and deaths in almost every country of the world. The lack of treatment options further worsens the scenario, thereby necessitating the exploration of already existing US FDA-approved drugs for their effectiveness against COVID-19. In the present study, we have performed virtual screening of nutraceuticals available from DrugBank against 14 SARS-CoV-2 proteins. Molecular docking identified several inhibitors, two of which, rutin and NADH, displayed strong binding affinities and inhibitory potential against SARS-CoV-2 proteins. Further normal model-based simulations were performed to gain insights into the conformational transitions in proteins induced by the drugs. The computational analysis in the present study paves the way for experimental validation and development of multi-target guided inhibitors to fight COVID-19.

1. Introduction

Starting from one patient in December 2019 at Wuhan city of China, COVID-19 has caused mayhem worldwide. As of 24 April 2020, SARS-CoV-2, the etiological agent of COVID-19 has infected 165,069,258 people causing 3,422,907 deaths globally (as of 21 May 2021) [1]. The situation is further worsened by variants of COVID-19 circulating in the global population, which have tremendously increased the transmission rate of the virus [2]. Apart from mutational, structural, and phylogenetic analyses of the SARS-CoV-2 genome, scientists have been centering on drug repurposing to develop therapeutics to combat SARS-CoV-2 contagion [3]. Various existing drugs, Remdesivir, Lopinavir/Ritonavir, Interferon beta-1a, Chloroquine/hydroxychloroquine, are under SOLIDARITY trial initiated by WHO for their inhibitory activity against different proteins of SARS-CoV-2 [4,5] nevertheless better targeted inhibitors are required for COVID-19 treatment.
SARS-CoV-2 genome is known to encode up to 14 open reading frames that translate to structural proteins, spike (S), membrane (M), envelope (E) and nucleocapsid (N); two huge non-structural proteins (NSPs) cleaving into sixteen smaller proteins along with nine accessory factors. The virus uses S protein to bind to the angiotensin-converting enzyme 2 (ACE2) receptor to enter the host cell. Both the M and E proteins are involved in forming the virus envelope and the pathogenesis of the virus, while the N protein binds to the virus’s RNA genome, creating the nucleocapsid [6]. NSPs form the replication/transcription complex that includes the papain-like proteinase (NSP3), the main proteinase (NSP5), the NSP7-NSP8 complex, the RNA-dependent RNA polymerase (NSP12), a NTPase/helicase (NSP13), an exonuclease (NSP14), an endonuclease (NSP15), and 2′O-methyltransferases (NSP16). Another group of SARS-CoV-2 is the accessory proteins, 3a, 3b, 6, 7a, 7b, 8, 9b, 9c and 10. The accessory proteins serve multitude of functions in virus replication [7]. In earlier studies, several SARS-CoV-2 proteins have been identified as having high mutational propensity, which further demands for multi-targeted inhibitors [8].
In the present study, virtual screening followed by normal model-based simulations was performed using nutraceuticals against 14 SARS-CoV-2 proteins; NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, 3a, S, E, M, 6, 7a and N that led to the identification of compounds that can be repurposed against SARS-CoV-2. The interaction analyses further exhibited the hydrogen bonding and hydrophobic network forming residues between protein and ligand. The findings of the proposed work will aid in global efforts for fight this pandemic by expediting drug development against COVID-19.

2. Methodology

2.1. Protein Structure Modelling and Preparation

The X-ray structures of some SARS-CoV-2 proteins, N, NSP3, NSP5, NSP9 and NSP15, are available on protein data bank (PDB) [9] corresponding to PDB ids 5MM3, 6W9C, 5RED, 6W4B and 6VWW, respectively. Homology modelling was performed for the remaining proteins using SWISS-MODEL server which identifies the structural templates followed by the alignment of query protein sequence with template structures and model generation. Prior to docking studies, the proteins were processed using Schrodinger Protein Preparation Wizard during which missing hydrogens were added, hydrogen bonds were optimized, and water molecules were deleted [10].

2.2. Nutraceuticals (Ligands) Preparation

Nutraceuticals are bioactive phytochemicals that deliver health benefits and are relatively safe for prevention and treatment against the disease. Nutraceuticals have already been used for treating several diseases, including atherosclerosis, cancer, cardiovascular diseases, diabetes, hypertension, inflammation, obesity, and others. A large number of nutraceuticals are common with FDA-approved pharmaceutical drugs. The library of 150 nutraceuticals was obtained from the DrugBank database [11]. The chemical structures of nutraceuticals were prepared with LigPrep module Schrodinger which generated diverse, accurate, and energy minimized conformations using the OPLS-2005 force field [10].

2.3. Docking Studies

A rigorous virtual screening was performed using Schrodinger Glide (grid-based ligand docking with energetics) [12]. Furthermore, docking studies for the top-scoring compounds were also done using GOLD (Genetic Optimization for Ligand Docking) [13] and AutoDock Vina [14] software.

2.4. Docking Studies Using Glide

Glide uses a range of filters and thoroughly searches the conformational, orientational, and positional space for ligand in the receptor’s binding site. The receptor was represented by creating a cubic grid of size 30 × 30 × 30 Å centered on the active site residues for each protein by Schrodinger’s Receptor Grid Generation program [12]. Further, the ligands were screened against each receptor using Glide’s high throughput virtual screening (HTVS) approach [12]. The high-ranking ligands were then subjected to stringent screening via the extra precision (XP) [12] method to eliminate false positives and acquire distinct binding modes of compounds.

2.5. Docking Studies Using GOLD

GOLD employs a genetic algorithm to examine the conformational flexibility of the ligand and partial (side-chain) flexibility of the receptor and generates the accurate binding mode of ligands [13]. The binding site for the ligands was specified as a ligand-specific pocket involving all the active site residues of protein within a 6 Å radius. The docked poses of ligand were evaluated using GoldScore fitness function, which considers protein-ligand hydrogen bonding and van der Waals energies along with ligand internal and torsional strain energy.

2.6. Docking Studies Using AutoDock Vina

AutoDock Vina uses a gradient optimization method to find low energy docked conformations of ligand [14]. A grid box of size 70 × 70 × 70 Å along X, Y and Z axes was generated considering active site residues of the protein. The lowest binding affinity ligand pose was separated and complexed with the receptor.

2.7. Interaction Analyses

LigPlus was used to compute interactions between the docked protein-ligand complexes. 2D diagrams were generated depicting hydrogen bonds and hydrophobically interacting residues.

2.8. Normal Mode Analyses

To study the conformational variations upon ligand binding in the protein, normal mode analyses (NMA) was performed using a rigid cluster NMA based NMSim web server [15]. NMA has been used to predict high amplitude conformational transitions and correlated atomic movements induced in the protein on the binding of a ligand, which largely occurs near the lowest energy state of the unbound protein [16,17]. The protein acquires a favorably compact state on binding a ligand, measured by radius of gyration (ROG) [18]. Hence, ROG-guided simulation implemented in NMSim was applied with default parameters. The NMSim approach is performed in three steps: a rigid cluster decomposition (RCD) is obtained in the first step (FIRST analyses) followed by calculation of normal modes using rigid cluster normal-mode analysis (RCNMA). These normal modes are then used in the simulation performed by NMSim [15]. The information regarding the conformational changes of the protein during the simulation was provided by root-mean-square-deviation (RMSD) and root-mean-square-fluctuation (RMSF) graphs.

3. Results and Discussion

3.1. Molecular Docking Analyses

The compounds that showed high binding affinity against SARS-CoV-2 proteins include Rutin (NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, ORF3a, S, E, M, ORF6, ORF7a, N); NADH (NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, ORF3a, S, E, M, ORF6, ORF7a, N); Ginsenoside Rg1 (NSP1, NSP3, NSP5, NSP9, NSP12, NSP15, E, ORF5, ORF7a and N); Ginsenoside Rb1 (NSP5, NSP12, ORF7a, NSP9, NSP15); Ginsenoside C (N, ORF6 and NSP1); Spermine (ORF6); Glutathione (NSP13); Ornithine (ORF3a), and α-tocopherol succinate (E). However, the top-scoring compounds against most of the SARS-CoV-2 proteins were Rutin (DB01698), NADH (DB00157), and Ginsenoside Rg1 (DB06750) (Table 1 and Table 2).

3.2. Simulations Analyses

The RMSDs of Cα atoms for proteins NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, ORF3a, S, E, M, ORF6, and N, in ligand-bound (closed) conformations vis-a-vis unbound (open in the absence of ligand) structures as a function of the number of conformations acquired by the proteins during the simulations, was computed after NMA (Figure 1). The graph clearly shows that large conformational transitions occurred at lower frequency (energy) modes in the case of all the proteins [19]. Among the 1500 conformations generated during each NMA run for every protein, the last 900–1500 confirmations represent the transition towards a closed structure. Upon comparison to ligand-free forms, as shown in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, it is evident that the proteins (Figure 1A, Figure 3A, Figure 4A, Figure 8A, Figure 9A, Figure 10A, Figure 12A, Figure 13A and Figure 14A) except NSP3 (Figure 2A), NSP12 (Figure 5A), NSP13 (Figure 6A), NSP15 (Figure 7A) and M (Figure 11A), did not undergo large movements to attain closed conformations.
Similar results were observed in the case of RMSF graphs of Cα atoms for open and closed structures. The magnitude of residue wise fluctuations was lower in the case of protein-ligand complexes, NSP1 (Figure 1B), NSP3 (Figure 2B), NSP5 (Figure 3B), NSP15 (Figure 7B), ORF3a (Figure 8B), S (Figure 9B) and ORF6 (Figure 12B) in contrast to the proteins alone pointing towards the stabilizing effects of ligands. The residues in NSP5, NSP15, S, and ORF6 proteins showed distinctively less degree of fluctuations in the presence of drugs. Moreover, superimposition of the open and closed conformations of proteins also revealed noticeably compact and folded proteins in the presence of drugs (Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, panel C). An entire overlap in RMSD (Figure 10A) and RMSF (Figure 10B) data points for unbound- and ligand-bound protein was observed in the case of E protein. However, superimposition of open and closed structures showed signs of protein folding and compactness in the presence of ligand (Figure 10C).

3.3. Protein-Ligand Interaction Analyses

Rutin formed six hydrogen bonds with NSP1 and NSP3, five with NSP5 and NSP9, twelve with NSP12, ten with NSP13, six with ORF3a, and nine with N protein. NADH majorly formed hydrogen bonds with NSP1 (seven), NSP3 (nine), NSP9 (twelve), NSP12 (eleven), NSP13 (thirteen), NSP15 (five), M (eight), ORF6 (five), and N protein (eight). Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, panels C and D illustrate the hydrogen bonds and hydrophobic interactions of NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, ORF3a, S, E, M, ORF6, ORF7a, and N proteins with Rutin and NADH, respectively. Figure 15 and Figure 16 illustrate the binding mode of Rutin and NADH with the 14 different SARS-CoV-2 proteins.
Several SARS-CoV-2 proteins with high mutational propensity have been identified, which further demands multi-targeted inhibitors [20]. Targeting multiple proteins with a single molecule would reduce the chances of resistance, and this approach will put forward an attractive model for drug development against COVID-19. Rutin, NADH, and Ginsenoside Rg1 showed high binding affinity against most of the SARS-CoV-2 proteins. Rutin is a flavonol glycoside found profusely in several plants and is a key component of nutritional supplements. Rutin has strong antioxidant properties and has been shown for its neuroprotective effect, anticarcinogenic and antidiabetic activities, treatment of cardiological and inflammatory diseases, and various other pharmacological activities [21,22]. NADH is a reduced coenzyme found extensively in nature and performs significant metabolic activities. NADH is being explored for its effectiveness in treating cardiovascular diseases, dementia-related Alzheimer’s and Parkinson’s disease, and chronic fatigue syndrome [23,24]. Several studies have mentioned Rutin as one of the potential inhibitors of SARS-CoV-2. Xu et al. 2020 showed the biding of Rutin with the main protease and common interaction sites, Asn142, Cys145, His164, Met165, Gln189, and Thr190 [25,26]. Another study reported interaction of Rutin with common interacting sites, Thr556, Tyr619, Lys621, Cys622, Asp623, Asn691, Asp761 of RNA-dependent RNA polymerase (NSP12); and Gly163, Arg166, Glu167, Asn267, Tyr273, of papain-like proteinase (NSP3) (Rahman et al., 2021). NADH has also been shown to be a possible inhibitor of the main protease and spike proteins of SARS-CoV-2 [27,28].
Ginsenoside Rg1 belongs to the ginsenosides class of compounds found in ginseng plants and has been used in traditional medicine for a long [29]. Ginsenoside Rg1 has been known to effect blood, cardiovascular, nervous, and immune systems performing various biological activities. The drug is also under clinical trials to treat dementia and cognitive impairments, rheumatic disorders, and stroke [30]. A large number of hydrogen bonding interactions and a robust hydrophobic network formed by these drugs with the proteins suggest their inhibitory effect on SARS-CoV-2. The candidate compounds proposed in the present study can be verified for inhibitory activity against SARS-CoV-2 and blocking viral-host interactions. This multi-targeted drug design approach will aid in global efforts to fight this pandemic by expediting drug development against COVID-19.

4. Conclusions

Studies conducted earlier in the development of SARS-CoV-2 have identified several proteins with a high mutational propensity, making multi-targeting inhibitors desirable. The present study involved the use of nutraceuticals against 14 different SARS-CoV-2 proteins, namely NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, ORF3a, S, E, M, ORF6, ORF7a, and N. Both virtual screening and model-based simulations have proved valuable in the identification of compounds that could be considered for the treatment of SARS-CoV-2. Consequently, hydrophobic networks and hydrogen bonds were observed between the protein and ligand as a consequence of the interaction analyses. The findings of this study will likely expedite the development of anti-COVID-19 drugs, which will contribute to global efforts to prevent the pandemic.

Author Contributions

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

Funding

The authors would like to extend their sincere appreciation to the Researchers Supporting Project number (RSP-2021/154), King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data available is provided in this paper.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. World Health Organization. Coronavirus Disease 2019 (COVID-19) Situation Report—87; WHO: Geneva, Switzerland, 2020. [Google Scholar]
  2. Gupta, V.; Haider, S.; Verma, M.; Singhvi, N.; Ponnusamy, K.; Malik, M.Z.; Verma, H.; Kumar, R.; Sood, U.; Hira, P.; et al. Comparative Genomics and Integrated Network Approach Unveiled Undirected Phylogeny Patterns, Co-mutational Hot Spots, Functional Cross Talk, and Regulatory Interactions in SARS-CoV-2. mSystems 2021, 6, e00030-21. [Google Scholar] [CrossRef] [PubMed]
  3. Mishra, C.B.; Pandey, P.; Sharma, R.D.; Malik, M.Z.; Mongre, R.K.; Lynn, A.M.; Prasad, R.; Jeon, R.; Prakash, A. Identifying the natural polyphenol catechin as a multi-targeted agent against SARS-CoV-2 for the plausible therapy of COVID-19: An integrated computational approach. Brief. Bioinform. 2021, 22, 1346–1360. [Google Scholar] [CrossRef] [PubMed]
  4. Dong, L.; Hu, S.; Gao, J. Discovering drugs to treat coronavirus disease 2019 (COVID-19). Drug Discov. Ther. 2020, 14, 58–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Kai Kupferschmidt, J.C. WHO launches global megatrial of the four most promising coronavirus treatments. Sci. Mag. 2020. Available online: https://www.science.org/content/article/who-launches-global-megatrial-four-most-promising-coronavirus-treatments (accessed on 12 November 2021). [CrossRef]
  6. Gordon, D.E.; Jang, G.M.; Bouhaddou, M.; Xu, J.; Obernier, K.; O’Meara, M.J.; Guo, J.Z.; Swaney, D.L.; Tummino, T.A.; Hüttenhain, R.; et al. A SARS-CoV-2-Human Protein-Protein Interaction Map Reveals Drug Targets and Potential Drug-Repurposing. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  7. Mousavizadeh, L.; Ghasemi, S. Genotype and phenotype of COVID-19: Their roles in pathogenesis. J. Microbiol. Immunol. Infect. 2020, 54, 159–163. [Google Scholar] [CrossRef]
  8. Banaganapalli, B.; Al-Rayes, N.; Awan, Z.A.; Alsulaimany, F.A.; Alamri, A.S.; Elango, R.; Malik, M.Z.; Shaik, N.A. Multilevel systems biology analysis of lung transcriptomics data identifies key miRNAs and potential miRNA target genes for SARS-CoV-2 infection. Comput. Biol. Med. 2021, 135, 104570. [Google Scholar] [CrossRef]
  9. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [Green Version]
  10. Madhavi Sastry, G.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J. Comput. -Aided Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef]
  11. Wishart, D.S.; Knox, C.; Guo, A.C.; Cheng, D.; Shrivastava, S.; Tzur, D.; Gautam, B.; Hassanali, M. DrugBank: A knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 2008, 36, D901–D906. [Google Scholar] [CrossRef]
  12. Halgren, T.A.; Murphy, R.B.; Friesner, R.A.; Beard, H.S.; Frye, L.L.; Pollard, W.T.; Banks, J.L. Glide: A new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J. Med. Chem. 2004, 47, 1750–1759. [Google Scholar] [CrossRef] [PubMed]
  13. Jones, G.; Willett, P.; Glen, R.C.; Leach, A.R.; Taylor, R. Development and validation of a genetic algorithm for flexible docking. J. Mol. Biol. 1997, 267, 727–748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
  15. Kruger, D.M.; Ahmed, A.; Gohlke, H. NMSim web server: Integrated approach for normal mode-based geometric simulations of biologically relevant conformational transitions in proteins. Nucleic Acids Res. 2012, 40, W310–W316. [Google Scholar] [CrossRef]
  16. Petrone, P.; Pande, V.S. Can conformational change be described by only a few normal modes? Biophys. J. 2006, 90, 1583–1593. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Rueda, M.; Bottegoni, G.; Abagyan, R. Consistent improvement of cross-docking results using binding site ensembles generated with elastic network normal modes. J. Chem. Inf. Model. 2009, 49, 716–725. [Google Scholar] [CrossRef]
  18. Dietzen, M.; Zotenko, E.; Hildebrandt, A.; Lengauer, T. On the applicability of elastic network normal modes in small-molecule docking. J. Chem. Inf. Model. 2012, 52, 844–856. [Google Scholar] [CrossRef]
  19. Tama, F.; Sanejouand, Y.H. Conformational change of proteins arising from normal mode calculations. Protein Eng. 2001, 14, 1–6. [Google Scholar] [CrossRef]
  20. Singh, R.K.S.; Malik, M.Z.; Singh, R.K.B. Diversity of SARS-CoV-2 isolates driven by pressure and health index. Epidemiol. Infect. 2021, 149, e38. [Google Scholar] [CrossRef]
  21. Ganeshpurkar, A.; Saluja, A.K. The Pharmacological Potential of Rutin. Saudi Pharm. J. 2017, 25, 149–164. [Google Scholar] [CrossRef] [Green Version]
  22. Al-Dhabi, N.A.; Arasu, M.V.; Park, C.H.; Park, S.U. An up-to-date review of rutin and its biological and pharmacological activities. EXCLI J 2015, 14, 59–63. [Google Scholar] [PubMed]
  23. Galsnapp, A.; Schaefer, E. Nicotinamide-Adenine Dinucleotide (NADH) in Parkinson’s Disease and Alzheimer’s Disease. Int. J. Pharm. Compd. 2000, 4, 276–279. [Google Scholar]
  24. Castro-Marrero, J.; Cordero, M.D.; Segundo, M.J.; Sáez-Francàs, N.; Calvo, N.; Román-Malo, L.; Aliste, L.; de Sevilla, T.F.; Alegre, J. Does oral coenzyme Q10 plus NADH supplementation improve fatigue and biochemical parameters in chronic fatigue syndrome? Antioxid. Redox Signal. 2015, 22, 679–685. [Google Scholar] [CrossRef] [Green Version]
  25. Rahman, F.; Tabrez, S.; Ali, R.; Alqahtani, A.S.; Ahmed, M.Z.; Rub, A. Molecular docking analysis of rutin reveals possible inhibition of SARS-CoV-2 vital proteins. J. Tradit. Complement. Med. 2021, 11, 173–179. [Google Scholar] [CrossRef]
  26. Xu, Z.; Yang, L.; Zhang, X.; Zhang, Q.; Yang, Z.; Liu, Y.; Liu, W. Discovery of Potential Flavonoid Inhibitors against COVID-19 3CL Proteinase Based on Virtual Screening Strategy. Front. Mol. Biosci. 2020, 7, 556481. [Google Scholar] [CrossRef] [PubMed]
  27. Hall, D.C., Jr.; Ji, H.F. A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease. Travel Med. Infect. Dis. 2020, 35, 101646. [Google Scholar] [CrossRef]
  28. White, M.A.; Lin, W.; Cheng, X. Discovery of COVID-19 Inhibitors Targeting the SARS-CoV2 Nsp13 Helicase. bioRxiv 2020. [Google Scholar] [CrossRef]
  29. Attele, A.S.; Wu, J.A.; Yuan, C.S. Ginseng pharmacology: Multiple constituents and multiple actions. Biochem. Pharmacol. 1999, 58, 1685–1693. [Google Scholar] [CrossRef]
  30. He, F.; Yu, C.; Liu, T.; Jia, H. Ginsenoside Rg1 as an Effective Regulator of Mesenchymal Stem Cells. Front. Pharmacol. 2019, 10, 1565. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Leader protein (NSP1). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 1. Leader protein (NSP1). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g001
Figure 2. Papain-like protease (NSP3). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 2. Papain-like protease (NSP3). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g002
Figure 3. Main protease (NSP5). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 3. Main protease (NSP5). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g003
Figure 4. RNA binding protein (NSP9). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 4. RNA binding protein (NSP9). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g004
Figure 5. RNA-dependent RNA polymerase (NSP12). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 5. RNA-dependent RNA polymerase (NSP12). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g005
Figure 6. Helicase (NSP13). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 6. Helicase (NSP13). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g006
Figure 7. Endoribonuclease (NSP15). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 7. Endoribonuclease (NSP15). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g007
Figure 8. ORF 3a. (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 8. ORF 3a. (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g008
Figure 9. Spike (S). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 9. Spike (S). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g009
Figure 10. Envelope (E). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 10. Envelope (E). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g010
Figure 11. Membrane (M). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 11. Membrane (M). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g011
Figure 12. ORF6. (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 12. ORF6. (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g012
Figure 13. ORF7a. (A) Molecular interactions of protein with Rutin (B) Molecular interactions of protein with NADH.
Figure 13. ORF7a. (A) Molecular interactions of protein with Rutin (B) Molecular interactions of protein with NADH.
Vaccines 10 00024 g013
Figure 14. Nucleocapsid (N). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Figure 14. Nucleocapsid (N). (A) RMSD plot of unbound and ligand-bound protein (B) RMSF plot of unbound and ligand-bound protein (C) Superimposition of unbound and ligand-bound protein (D) Molecular interactions of protein with Rutin (E) Molecular interactions of protein with NADH.
Vaccines 10 00024 g014
Figure 15. The binding mode of rutin (pink) and NADH (yellow) with the 14 different SARS-CoV-2 proteins. (A) NSP1. (B) NSP3. (C) NSP5. (D) NSP9. (E) NSP12. (F) NSP13. (G) NSP15. (H) ORF3a. (I) S. (J) E. (K) M. (L). ORF6. (M) ORF7a. (N) N Proteins.
Figure 15. The binding mode of rutin (pink) and NADH (yellow) with the 14 different SARS-CoV-2 proteins. (A) NSP1. (B) NSP3. (C) NSP5. (D) NSP9. (E) NSP12. (F) NSP13. (G) NSP15. (H) ORF3a. (I) S. (J) E. (K) M. (L). ORF6. (M) ORF7a. (N) N Proteins.
Vaccines 10 00024 g015
Figure 16. The binding of rutin (pink) and NADH (yellow) in the pockets of 14 different SARS-CoV-2 proteins. (A) NSP1. (B) NSP3. (C) NSP5. (D) NSP9. (E) NSP12. (F) NSP13. (G) NSP15. (H) ORF3a. (I) S. (J) E. (K) M. (L). ORF6. (M) ORF7a. (N) N Proteins.
Figure 16. The binding of rutin (pink) and NADH (yellow) in the pockets of 14 different SARS-CoV-2 proteins. (A) NSP1. (B) NSP3. (C) NSP5. (D) NSP9. (E) NSP12. (F) NSP13. (G) NSP15. (H) ORF3a. (I) S. (J) E. (K) M. (L). ORF6. (M) ORF7a. (N) N Proteins.
Vaccines 10 00024 g016
Table 1. Activity of drugs against SARS-CoV-2 proteins.
Table 1. Activity of drugs against SARS-CoV-2 proteins.
DrugDescriptionActivity against ORFs
RutinAn existing USFDA-approved drug used for strengthening weakened capillaries. Additionally, it has a powerful antioxidant with potential biological effect in reducing post-thrombotic syndrome, veins insufficiency or endothelial dysfunction.NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, ORF3a, Spike, Envelope, Membrane, ORF6, ORF7a, Nucleocapsid
NADHNADH plays essential metabolic roles and has been used to combat chronic fatigue syndrome. It is also being explored to be used against dementia and improving mental health.NSP1, NSP3, NSP5, NSP9, NSP12, NSP13, NSP15, ORF3a, Spike, Envelope, Membrane, ORF6, ORF7a, Nucleocapsid
Ginsenoside Rg1, protopanaxatriolGinsenoside is a major component of the root and stem of ginseng plant. It possesses a broad spectrum of pharmacological properties such as neuroprotection, anti-inflammation, anti-aging, anti-fatigue and memory-enhancing properties.NSP5, Nucleocapsid, NSP1, Envelope, NSP12, ORF5, ORF7a, NSP3, NSP9, NSP15
Table 2. Virtual screening results of nutraceuticals against SARS-CoV-2 proteins.
Table 2. Virtual screening results of nutraceuticals against SARS-CoV-2 proteins.
ProteinDrug NameGlide Score (kcal/mol)GOLD ScoreAutoDock Score (kcal/mol)
1Spike (S)NADH−11.3180.76−8.4
Rutin−9.9495.52−6.7
2Main protease (NSP5)Ginsenoside Rb1−9.37132.14−19.9
Rutin−9.5896.35−7.4
NADH−8.6585.20−8.5
Ginsenoside Rg1−8.16110.03−10.7
3Nucleocapsid (N)Rutin−9.3487.13−5.6
Ginsenoside C−8.82122.35−12.0
NADH−8.1370.55−8.6
Ginsenoside Rg1−5.17110.79−10.7
4ORF6Ginsenoside C−7.51109.10−13.5
Spermine−7.2143.27−6.8
Rutin−6.5987.57−5.7
NADH−4.9860.11−8.6
5Leader protein (NSP1)Ginsenoside C−9.30107.85−11.6
NADH−7.3865.42−7.7
Rutin−7.0980.41−4.7
Ginsenoside Rg1−5.8185.13−9.4
6Envelope (E)Ginsenoside Rg1−8.3094.13−11.5
α−tocopherol succinate−8.1139.97−16.7
Rutin−3.8694.13−11.5
NADH−3.0868.11−9.5
7RNA−dependent RNA polymerase (NSP12)Ginsenoside Rb1−11.00123.74−19.6
Rutin−10.98113.81−7.0
Ginsenoside Rg1−10.14143.41−12.2
NADH−9.0077.89−9.5
8ORF 3aRutin−11.4793.67−7.5
Ornithine−9.1039.31−6.3
NADH−3.3491.53−10.0
9Membrane (M)Ginsenoside Rg1−9.3581.77−9.4
Rutin−7.6277.90−5.5
NADH−5.4058.68−8.0
10ORF 7aGinsenoside Rb1−8.70166.98−16.8
Ginsenoside Rg1−7.00135.96−10.7
NADH−6.8378.90−7.6
Rutin−4.96126.00−5.0
11Papain-like protease (NSP3)Ginsenoside Rg1−10.37103.78−10.2
Rutin−8.2291.30−6.7
NADH−6.6380.28−8.0
12Helicase (NSP13)Glutathione−10.3160.00−6.4
NADH−9.5487.21−10.5
Rutin−7.8585.62−7.1
13RNA binding protein (Orf1ab, nsp9)Ginsenoside Rb1−7.92102.08−20.1
Ginsenoside Rg1−6.8491.46−9.6
Rutin−6.5670.25−5.1
NADH−5.6655.60−7.9
14Endoribonuclease (NSP15)Ginsenoside Rb1−13.50110.76−24.4
Ginsenoside Rg1−11.00111.16−12.5
NADH−9.9578.94−9.2
Rutin−9.60100.05−7.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alanazi, K.M.; Farah, M.A.; Hor, Y.-Y. Multi-Targeted Approaches and Drug Repurposing Reveal Possible SARS-CoV-2 Inhibitors. Vaccines 2022, 10, 24. https://0-doi-org.brum.beds.ac.uk/10.3390/vaccines10010024

AMA Style

Alanazi KM, Farah MA, Hor Y-Y. Multi-Targeted Approaches and Drug Repurposing Reveal Possible SARS-CoV-2 Inhibitors. Vaccines. 2022; 10(1):24. https://0-doi-org.brum.beds.ac.uk/10.3390/vaccines10010024

Chicago/Turabian Style

Alanazi, Khalid Mashay, Mohammad Abul Farah, and Yan-Yan Hor. 2022. "Multi-Targeted Approaches and Drug Repurposing Reveal Possible SARS-CoV-2 Inhibitors" Vaccines 10, no. 1: 24. https://0-doi-org.brum.beds.ac.uk/10.3390/vaccines10010024

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