Structure-based virtual screening of bioactive compounds from Indonesian medical plants against severe acute respiratory syndrome coronavirus-2.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a virus that causes the infectious disease coronavirus disease-2019. Currently, there is no effective drug for the prevention and treatment of this virus. This study aimed to identify secondary metabolites that potentially inhibit the key proteins of SARS-CoV-2. This was an in silico molecular docking study of several secondary metabolites of Indonesian herbal plant compounds and other metabolites with antiviral testing history. Virtual screening using AutoDock Vina of 216 Lipinski rule-compliant plant metabolites was performed on 3C-like protease (3CLpro), RNA-dependent RNA polymerase (RdRp), and spike glycoprotein. Ligand preparation was performed using JChem and Schrödinger's software, and virtual protein elucidation was performed using AutoDockTools version 1.5.6. Virtual screening identified several RdRp, spike, and 3CLpro inhibitors. Justicidin D had binding affinities of -8.7, -8.1, and -7.6 kcal mol-1 on RdRp, 3CLpro, and spike, respectively. 10-methoxycamptothecin had binding affinities of -8.5 and -8.2 kcal mol-1 on RdRp and spike, respectively. Inoxanthone had binding affinities of -8.3 and -8.1 kcal mol-1 on RdRp and spike, respectively, while binding affinities of caribine were -9.0 and -7.5 mol-1 on 3CLpro and spike, respectively. Secondary metabolites of compounds from several plants were identified as potential agents for SARS-CoV-2 therapy. Copyright:

INTRODUCTION

Coronavirus Disease-2019 (COVID-19) is a disease similar to pneumonia that caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Since the first case on December 12, 2019, the number of cases and deaths due to COVID-19 has continued to increase. As of December 27, 2020, there have been 79.2 million cases, and over 1.7 million deaths have been reported.[1]

SARS-CoV-2 is a positive-strand RNA virus. Patients infected with this virus experience symptoms such as fever, dry cough, and shortness of breath. However, many patients are carriers and show no symptoms. This is one of the reasons for the widespread transmission of this virus.[2] The treatment of COVID-19 patients currently uses drugs that have previously received approval from drug regulatory agencies, such as hydroxychloroquine, chloroquine, and remdesivir.[3]

The SARS-CoV-2 genome encodes several proteins that are important in the virulence process. These proteins include SARS-CoV-2 3C-like protease (3CLpro) and RNA-dependent RNA polymerase (RdRp), which are useful for viral replication and transcription, as well as spike glycoproteins that mediate virus-host binding and viral fusion.[2] Inhibition of these proteins in SARS-CoV-2 could lead to a potential therapy for COVID-19.

This study aimed to identify secondary metabolites that potentially inhibit the key proteins of SARS-CoV-2 by virtual screening of plant metabolites against the 3CLpro, RdRp proteins, and spike glycoprotein. These metabolites were derived from several species of Indonesian traditional medicine as well as other metabolites that had an antiviral testing history.

MATERIALS AND METHODS

Ligand and protein preparation

The ligands used were metabolites of Cymbopogon citratus, Kaempferia galangal, Curcuma longa, Curcuma xanthorrhiza, and Zingiber officinale. Plant metabolites that had a history of antiviral testing were also used as ligands. All compounds were obtained from http://www.knapsackfamily.com. Screening using SwissADME (http://www.swissadme.ch/index.php) according to the five Lipinski rules was performed on 240 metabolites. The ligand structure was formed using JChem software and optimized using the LigPrep module from Schrödinger's software (Schrödinger LLC, New York, NY, USA). Atomic protonation was adjusted to pH 7.0 with Epik software (Schrödinger LLC), and geometry optimization used OPLS_2005 Force Field software (Schrödinger LLC).

Protein structures targeted for docking were obtained from the Protein Data Bank (https://www.rcsb.org/). The proteins used were 3CLpro (PDB ID: 6M2N chain A), RdRp (PDB ID: 6M71 chain A), and spike glycoprotein (PDB ID: 6VXX chain B). Protein preparation was performed using AutoDockTools version 1.5.6. Protein validation for 3CLpro was performed by calculating the root mean square deviation (RMSD) of the protein with its co-crystalline ligand using PyMOL software version 2.3.4 (Schrödinger LLC).

Virtual screening and molecular docking

Virtual screening was performed using an Intel® Celeron® 2955U at 1.40 GHZ processor, 2.00 GB of RAM 64-bit operating system. Targeted docking was performed on the 3CLpro protein. The grid box was 40 Å × 40 Å × 40 Å centered at −32,981, −65,436, and 41.404. The validity of the docking protocol was assessed based on the RMSD. For the RdRp protein and spike glycoprotein, blind docking and targeted docking were performed. Blind docking was useful for finding receptor pockets that bind positive control compounds (remdesivir, hydroxychloroquine, and chloroquine). The chains targeted for blind docking were chains that had binding pockets with high (top five) drug score, as determined by DoGSiteScorer (https://proteinsplus.zbh.uni-hamburg.de/#dogsite), which were chain A of RdRp and chain B of spike. Blind docking was performed by making the grid box sufficient to cover all sides of the protein, and the docking was performed three times in a row. Remdesivir had the lowest binding affinity at both the receptors and was consistent in one binding pocket, so the remdesivir position after docking was used as a grid box for RdRp and spike glycoprotein. The grid box for RdRp was 40 Å × 40 Å × 40 Å centered at 132,754, 112,571, and 103,907, while spike glycoproteins was 40 Å × 40 Å × 40 Å centered at 175,288, 177,632, and 235,662. Molecular docking was performed using AutoDock Vina, and the ligand–receptor interaction was visualized using Discovery Studio Visualizer v17.2.016349 software (Dassault Systèmes, San Diego, California, USA).

RESULTS

Two hundred and sixteen compounds that complied with Lipinski's rules were subjected to virtual screening against 3CLpro, RdRp, and the spike glycoprotein [Figure 1]. The compounds that potentially become inhibitors of the proteins were compounds that had same or lower than each benchmark's binding affinities [Table 1]. Virtual screening showed that four compounds had the potential to become RdRp inhibitors [Figure 2 and Table 2]. On the other hand, there were 79 compounds which had the potential to become spike inhibitors. Of these 79 compounds, 11 were selected with the lowest binding affinities [Figure 3 and Table 3]. While for the 3CLpro receptor, they were 11 compounds that potentially inhibit the protein [Figure 4 and Table 4].

Figure 1

Cartoon representation of spike glycoprotein (a), 3C-like protease (b), and RNA-dependent RNA polymerase (c) of Severe acute respiratory syndrome coronavirus-2, with their binding sites and selected grid box as a targeted docking. The binding site were determined with DoGSiteScorer

Table 1

Binding affinity (kcal/mol) of anti-severe acute respiratory syndrome-coronavirus-2 comparator on the targeted protein RNA-dependent RNA polymerase, spike glycoprotein, and 3C-like protease

Compounds	Receptor
	RdRp	Spike	3CLpro
5,6,7-trihydroxy-2-phenyl-4H-chromen- 4-one	ND	ND	−7.8
Remdesivir	−8.2	−6.6	−6.7
Hydroxychloroquine	−6.7	−5.7	−6.5
Chloroquine	−5.8	−5.9	−6.0

ND: Not determined, 3CLpro: 3C-like protease, RdRp: RNA-dependent RNA polymerase

Figure 2

Molecular interaction of selected ligands with RNA-dependent RNA polymerase of severe acute respiratory syndrome coronavirus-2. Remdesivir (a), hydroxychloroquine (b), chloroquine (c), justicidin D (d), 10-methoxycamptothecin (e), inoxanthone (f), 3-O-Caffeoylquinic acid (g)

Table 2

Binding affinity of anti-severe acute respiratory syndrome-coronavirus-2 RNA-dependent RNA polymerase candidates

C_ID	Compounds	Plant	Binding affinity (kcal/mol)
C00000720	Justicidin D	J. procumbens[5]	−8.7
C00029372	10-Methoxycamptothecin	C. acuminata[6]	−8.5
C00035841	Inoxanthone	C. inophyllum[7]	−8.3
C00002724	3-O-Caffeoylquinic acid	Schefflera heptaphylla. L.[4]	−8.2

J. procumbens: Justicia procumbens, C. acuminata: Camptotheca acuminata, C. inophyllum: Calophyllum inophyllum

Figure 3

Molecular interaction of selected ligands with spike glycoprotein of severe acute respiratory syndrome coronavirus-2. Remdesivir (a), hydroxychloroquine (b), chloroquine (c), (−)-beta-Sitosterol (d), 10-methoxycamptothecin (e), inoxanthone (f), emetine (g), alpha-cedrene (h), zonarene (i), justicidin D (j), bis-demethoxycurcumin (k), curcumenol (l), beta-caryophyllene epoxide (m), caribine (n)

Table 3

Binding affinity of anti-Severe acute respiratory syndrome-coronavirus-2 spike glycoprotein candidates

C_ID	Compounds	Plant	Binding affinity (kcal/mol)
C00003672	(−)-beta-sitosterol	Z. officinale, T. capitatus[8,9]	−8.4
C00029372	10-Methoxycamptothecin	C. acuminata[6]	−8.2
C00035841	Inoxanthone	C. inophyllum[7]	−8.1
C00001849	Emetine	C. ipecacuanha[10]	−7.7
C00003111	Alpha-cedrene	Z. officinale[11]	−7.6
C00032554	Zonarene	Z. officinale[12]	−7.6
C00000720	Justicidin D	J. procumbens[4]	−7.6
C00032769	Bis-demethoxycurcumin	C. longa[13]	−7.5
C00020351	Curcumenol	C. longa[14]	−7.5
C00012483	Beta-caryophyllene epoxide	C. zanthorrhiza[15]	−7.5
C00001566	Caribine	H. arenicola[16]	−7.5

Z. officinale: Zingiber officinale, T. capitatus: Thymus capitatus, C. acuminata: Camptotheca acuminata, C. inophyllum: Calophyllum inophyllum, J. procumbens: Justicia procumbens, C. longa: Curcuma longa, C. zanthorrhiza: Curcuma zanthorrhiza, H. arenicola: Hymenocallis arenicola, C. ipecacuanha: Carapichea ipecacuanha

Figure 4

Molecular interaction of selected ligands with 3C-like protease of severe acute respiratory syndrome coronavirus-2. 5,6,7-trihydroxy-2-phenyl-4H-chromen-4-one (a), remdesivir (b), hydroxychloroquine (c), chloroquine (d), caribine (e), cryptopleurine (f), justicidin D (g), diphyllin (h), quercetin 3-O-methyl ether (i), quercetin (j), dihydroquercetin (k), lycorine (l), luteolin (m), moralbanone (n), (2R,3R)-fustin (o)

Table 4

Binding affinity of anti-severe acute respiratory syndrome-coronavirus-2 3C-like protease candidates

C_ID	Compounds	Plant	Binding affinity (kcal mol−1)
C00001566	Caribine	H. arenicola[16]	−9.0
C00002331	Cryptopleurine	T. indica[17]	−8.2
C00000720	Justicidin D	J. procumbens[7]	−8.1
C00002600	Diphyllin	J. gendarussa[18]	−8.0
C00004632	Quercetin 3-O-methyl ether	Widespread in medicinal plants	−7.9
C00004631	Quercetin	Widespread in medicinal plants	−7.9
C00000677	Dihydroquercetin	Widespread in medicinal plants	−7.9
C00001576	Lycorine	L. radiata[19]	−7.9
C00000674	Luteolin	Widespread in medicinal plants	−7.9
C00013411	Moralbanone	Morus alba L.[20]	−7.8
C00000963	(2R,3R)-fustin	R. verniciflua[21]	−7.8

H. arenicola: Hymenocallis arenicola, T. indica: Tylophora indica, Justicia procumbens: J. procumbens, J. gendarussa: Justicia gendarussa, L. radiata: Lycoris radiata, R. verniciflua: Rhus verniciflua

DISCUSSION

RNA-dependent RNA polymerase

As with other positive-strand RNA viruses, genome replication and transcription of SARS-CoV-2 are dependent on RdRp.[2] In this study, the potential SARS-CoV-2 RdRp inhibitors were 3-O-caffeoylquinic acid, justicidin D, 10-methoxycamptothecin, and inoxanthone. 3-O-caffeoylquinic acid is a caffeoylquinic acid derivative with a history of antiviral activity testing against respiratory syncytial virus; however, its antiviral activity is not as potent as the other two derivatives of caffeoylquinic acid.[4] Justicidin D is a metabolite isolated from Justicia procumbens var. leucantha and is known to have antiviral activity in vesicular stomatitis virus.[5] However, there have been no recent reports of these compounds. Moreover, methoxycamptothecin is a camptothecin derivative found in Camptotheca acuminata and has been reported to have antiviral activity.[6] In addition, inoxanthone is a xanthone derivative found in Calophyllum inophyllum and has known cytotoxicity. However, reports of antiviral activity of inoxanthone compound have not been found.[78]

Spike glycoprotein

Spike (S) glycoprotein is the protein that mediated the entry of SARS-CoV-2 into the host.[3] Based on the present study, some of the compounds that potentially inhibit the spike protein of SARS-CoV-2 were emetine and carbine. Emetine is an alkaloid found in ipecacuanha (ipecac) root and recently was reported to show activity against SARS-CoV-2 in vitro.[10] Carbine is a metabolite isolated from Hymenocallis arenicola, but no recent reports on this compound have been found.[16]

Three metabolites of Z. officinale are also potential spike protein inhibitors: (−)-beta-sitosterol, alpha-cedrene, and zonarene.[1112] (−)-beta-sitosterol is also found in Thymus capitatus and known to have activity against herpes simplex virus type 2 (HSV-2).[9] Alpha-cedrene is a compound that improves and prevents liver steatosis in vivo.[22] On the other hand, zonarene is also found in essential oils and is known to have antifungal, mosquito larvicidal, and biting deterrent activities.[23] Based on the present study, Z. officinale may be useful as an herbal treatment for SARS-CoV-2. Nonetheless, further research for this is needed

Furthermore, two metabolites of C. longa, bis-demethoxycurcumin and curcumenol, are also potential inhibitors of spike glycoproteins. These compounds are known to have anticancer activity. However, activity against SARS-CoV-2 has been more widely reported for its analog, curcumin.[1314] Moreover, one of the metabolites of C. xanthorrhiza, beta-caryophyllene epoxide, is also a potential inhibitor of the spike glycoprotein. This compound was previously reported to have activity against HSV type 1 in vitro.[1524]

10-methoxycamptothecin, inoxanthone, and justicidin D were in the 11 compounds that had the lowest binding affinities for the spike. These compounds also had good binding affinity to the RdRp, as has been previously stated. Therefore, broad antiviral activity of these three compounds was possible in SARS-CoV-2.

3C-like protease

3CLpro is a protease that cuts polyproteins. Based on the present study, quercetin, quercetin 3-O-methyl ether and dihydroquercetin, has the potential to be 3CLpro SARS-CoV-2 inhibitors. These compounds are known to have various activities, including antiviral against SARS-CoV.[25] Lycorine and luteolin also have the potential to inhibit this protein. Lycorine is an alkaloid derived from Lycoris radiata, whereas luteolin is a metabolite that is widespread in medicinal plants.[19] Lycorine has antiviral activity against SARS-CoV-2 in vitro, with an EC50 of 0.31 μM.[1926] On the other hand, luteolin has antiviral activity against SARS-CoV.[27]

In this study, cryptopleurine, diphyllin, fustin, and moralbanone potentially inhibit 3CLpro. Cryptopleurine is an alkaloid isolated from Tylophora indica and has been shown to have antiviral activity.[17] Diphyllin is a metabolite of Justicia gendarussa, a ATPase blocker that has activity against viruses.[1828] Fustin is a metabolite of Rhus verniciflua and was reported to have antiviral activity in vitro.[2129] Furthermore, moralbanone, a compound isolated from Morus alba L, also has the potential to inhibit 3CLpro.[20] However, no recent reports of moralbanone have been found.

Furthermore, virtual screening suggested that caribine and justicidin D had the potential to inhibit the 3CLpro receptor. As previously explained, this compound also may inhibit both the RdRp and spike receptors, so these two compounds may have broad-spectrum antiviral activity by targeting those three receptors.

CONCLUSION

Virtual screenings targeting the RdRp, spike glycoprotein, and 3CLpro receptors showed that several plant metabolites had the same or lower binding affinities of compounds currently used as SARS-CoV-2 therapeutic agents. Several compounds such as emetine, lycorine, luteolin and quercetin are widely known to have antiviral activity. On the other hand, other compounds such as caribine and justicidin D that are also predicted to have broad antiviral SARS-CoV-2 activity have not much pharmacological activities reported. These metabolites could potentially be used in the management of COVID-19. However, further studies are needed.

Financial support and sponsorship

Nil.

Conflicts of interest

There are no conflicts of interest.