key: cord-0932458-jolt5442
authors: Babaeekhou, Laleh; Ghane, Maryam; Abbas-Mohammadi, Mahdi
title: In silico targeting SARS-CoV-2 spike protein and main protease by biochemical compounds
date: 2021-09-22
journal: Biologia (Bratisl)
DOI: 10.1007/s11756-021-00881-z
sha: ced763528d3c6dc44c7b92dc671c43ef2390eda0
doc_id: 932458
cord_uid: jolt5442

Since there is no general agreement on drug treatment of SARS-CoV-2, the search for a new drug capable of treating COVID-19 is of utmost priority. This study aims to dereplicate the chemical compounds of the methanol extract of Salvia officinalis and Artemisia dracunculus, and assay the inhibitory effect of these compounds as well as the previously dereplicated components of Zingiber officinale against SARS-CoV-2 in an in-silico study. A molecular networking (MN) technique was applied to find the chemical constituents of the extracts. Docking analysis was also used to find the binding affinity of dereplicated components from S. officinalis, A. dracunculus, and Z. officinale to COV-2-SP and M(pro). 57 compounds were dereplicated from the MeOH extracts of S. officinalis and A. dracunculus which include the class of polyphenols, flavonoids, coumarins, phenylpropanoids, anthocyanins, and dihydrochalcones. Molecular docking analysis indicated a high affinity of about 27 compounds from three mentioned plants against studied targets. kaempferol 3-O-rutinoside, neodiosmin, and querciturone with docking score values of -10.575, -10.208, and − 9.904 Kcal/mol and k(i) values of 0.016606, 0.030921, and 0.051749, respectively were found to have the highest affinities against COV-2-SP. 2-phenylethyl beta-primeveroside, curcumin PE, and kaempferol 3-O-rutinoside also indicated the highest affinity against M(pro) with docking scores of -10.34, -10.126 and − 9.705 and k(i) values of 0.024726, 0.035529, and 0.072494, respectively. MN can be successfully used for the dereplication of metabolites from plant extracts. In addition, the in-silico binding energies introduced several inhibitors from Z. officinale, S. officinalis, and A. dracunculus for the treatment of SARS-CoV-2 disease. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11756-021-00881-z.

The novel strain of coronavirus (CoV) Which first emerged in Wuhan, China was identified at the end of 2019 (2019-nCoV) (Zhu et al. 2020 ) and then officially named severe acute respiratory syndrome-related coronavirus (SARS-CoV-2). Different variants of SARS-CoV-2 with different transmission and disease characteristics, and different impact on vaccine efficacy has posed one of the biggest threats to global health. Based on the WHO weekly report on Coronavirus disease 2019 (COVID-19) released on 29 June 2021, the number of confirmed cases and deaths worldwide has reached over 180 million and about 4 million respectively. Alongside tracking the newly emerged variants of the virus, one of the major priorities is the evaluation of existing vaccines for efficacy against variants (WHO organization 2021). Until 1st July 2021, 23.4 % of the world population has received at least one dose of each vaccine and in low-income countries, only 0.9 % have received at least one dose (https:// ourworldindata.org). Existing therapies like Remdesivir (Veklury) have failed for the treatment of severe forms of the disease (Goldman et al. 2020 ) and until widespread and confirmed immunity against COVID-19, prevention of further disease spread and novel therapies are needed (Voysey et al. 2020) . So, drug development should progress based on SARS-CoV-2 different molecular targets.

Coronavirus entry into host cells is initiated by binding the envelope spike glycoprotein to the cell surface receptor angiotensin-converting enzyme 2 (ACE-2) (Dimitrov 2004; Li et al. 2003) . S is a class I viral fusion protein (1,300 amino acids) that trimerizes upon its folding. It is composed of two main subunits: S1, in amino, and S2 in the carboxy-terminal. S1 includes the receptor-binding domain and S2 drives membrane fusion. In most coronaviruses, there is a cleavage site at the joining point of S1 and S2 and a proteolytic cleavage happens by host proteases at the S2 cleavage site (Bosch et al. 2003; Du et al. 2009 ). But the S1 and S2 subunits are still connected in the pre-fusion form of the S trimer. After virion attachment to the host cell receptor, a second essential cleavage by endo-lysosomal proteases occurs at the S2′ cleavage site, allowing the release of the internal fusion peptide (FP) and fusion of the spike protein envelope into the host membrane and transition of S2 into the post-fusion structure (Burkard et al. 2014; Li 2016) . So, in the viral entry process, the spike protein shows two different forms: pre-fusion or the form which is seen on mature virions, and post-fusion which is formed after membrane fusion (Shang et al. 2018; Song et al. 2018; Walls et al. 2016) . In previous studies, the potential receptor usage of the SARS-CoV-2 spike protein (COV-2-SP) is analyzed and is shown that the new strain also uses ACE2 as its receptor. This is because the sequence of COV-2-SP receptor-binding domain (RBD) that binds to ACE-2, is similar to that of severe acute respiratory syndrome coronavirus (SARS-CoV) . It is also shown that COV-2-SP RBD has improved its binding affinity to human ACE-2 by residue changes at RBD-receptor interaction spots . These findings make COV-2-SP a potential candidate that can be specifically targeted by entry blocking inhibitor drugs.

In the human coronavirus (RNA positive-stranded) replication cycle, two overlapping polyproteins that are, replicase 1a, and replicase 1ab, are encoded by the 229E replicase gene (Herold et al. 1993) . These proteins continue replication and transcription in the viral replication cycle but for the production of regulatory non-structural polypeptides from polyproteins, a 33.1-kD HCoV 229E main proteinase (M pro ) and papain-like protease (PLP) is essential (Thiel et al. 2001; Ziebuhr et al. 1995) . Because of the similarity of the cleavage site of M pro with picornavirus 3 C proteinases, it's also called 3 C-like proteinase (3CL pro ) (Anand et al. 2002) . The coronavirus M pro comprises three structural domains. Domains I (residues 8-99) and II (residues 100-183) are antiparallel β barrels representing the chymotrypsin catalytic domain. The substrate-binding site is located in a motif between these two domains. Domain III with five helices is located on the C-terminal of the enzyme (residues 200-300) and contains the proteolytic site. This latter domain is connected to domain II with a long loop (residues 184-199) (Anand et al. 2003; Sirois et al. 2007 ). There is no human protease similar to the cleavage of M pro so it is a suitable target for controlling coronaviruses (Anand et al. 2003) .

In our previous study, chemical compounds of Zingiber officinale were identified using Molecular Networking (MN) (Babaeekhou and Ghane 2020) . Considering the importance of a combination of in silico and experimental studies in drug discovery, in the present study, an MN technique based on an untargeted MS/MS analysis was used to find the chemical composition of the methanol extracts of S. officinalis and A. dracunculus. Then, the molecular docking strategy was applied to dereplicated compounds of Z. officinale, S. officinalis, and A. dracunculus to evaluate their binding energies with COV-2-SP and M pro viral targets. In continue, in this study, the structure of the SARS-CoV-2 chimeric receptor-binding domain complexed with its receptor human ACE2 (PDB: 6VW1) ); pre-fusion 2 0 1 9 -n C o V s p i k e g l y c o p r o t e i n w i t h a s i n g l e receptor-binding domain up (PDB: 6VSB) (Wrapp et al. 2020) ; and the crystal structure of COVID-19 M pro (PDB IDs 6LU7 and 6M03) (Jin et al. 2020) were selected as viral targets.

200 g of whole plants of S. officinalis and A. dracunculus were powdered and extracted using 200 mL methanol solvent (three times) via a maceration method. The extracts were then combined and evaporated to dryness under reduced pressure to obtain solvent-less extracts.

Samples were prepared in methanol at a concentration of 2 mg/mL and 5 µl were injected for each LC-MS analysis. LC-MS/MS analyses were carried out on a Waters Acquity UPLC system equipped with a Waters Xevo™ QToF mass spectrometer and an electrospray source. Mass spectrometry (MS) data were acquired simultaneously in the positive mode at a mass range of m/z 1000-1000 Da. Samples were dissolved in methanol at the concentration of 2 mg/mL and injected into a 5 μm SunFire™ C-18 column (250 × 4.6 mm). Two solvents, water + 0.1 % acetic acid (solvent A) and methanol + 0.1 % acetic acid (solvent B) were used as the mobile phase at a flow rate of 0.4 mL/min. The ESI conditions and mass spectrometry acquisition parameters were set as follows: capillary voltage, detector voltage, sampling cone, and extraction cone voltages: 3.0 kV, 2.2 kV, 25 V, and 4.0 V respectively; source and desolvation temperatures: 150 and 250 o C, respectively and cone and desolvation flow: 50 and 600 L/h, respectively. Eight most intense ions with a threshold higher than 50 were selected for data-dependent MS/MS survey scans.

The MS/MS data were converted into the mzXML format using MSconvert software for spectral data processing. The mzXML data were uploaded to the Global Natural Products Social (GNPS) MN web server (http://gnps.ucsd.edu) and analyzed using the MN workflow. MS/MS spectra were filtered by choosing only the top 6 fragment ions in the +/-50Da window throughout the spectrum. Then, to create consensus spectra, a diversity of parameters such as precursor ion mass tolerance (0.2 Da), MS/MS fragment ion tolerance (0.08 Da), minimum cosine score (0.6), and minimum matched fragment ions (4 peaks) were applied. When each of the nodes appeared in each other's respective top 10 most similar nodes, edges between two nodes were kept in the network. Then, the spectra in the network were searched against GNPS' spectral libraries. Cytoscape 2.8.3 was also carried out to visualize the data as a network of nodes and edges.

Preparation of ligands and targeted enzymes, as well as the molecular docking analysis, was accomplished on the Glide of Schrodinger package 2016-2 (Vasavi et al. 2017) . The structure of COV-2-SP and M pro were taken from Protein Data Bank (PDB). They were as follows: prefusion 2019-nCoV spike glycoprotein with a single receptor-binding domain up (PDB ID: 6VSB); 2019-nCoV chimeric receptor-binding domain complexed with its receptor human ACE2 (PDB ID: 6VW1); COVID-19 main protease in the apo form (PDB ID: 6M03) and COVID-19 main protease complex with an inhibitor N3 (PDB ID: 6LU7). Nelfinavir and Lopinavir were used as interaction assessment indicators. All protein structures were prepared on protein preparation wizard in Maestro by removing the crystallized ligands, all free water molecules, and complexed proteins with targets in PDB such as human ACE2 and inhibitor N3 chains. Then energy minimization was done. The grid box of enzymes was created with a Grid generation application at particular residues of the proteins obtained from the DEPTH server. The chemical structures of ligands were drawn on ChemDraw Professional 15.0 and conformationally optimized on the LigPrep module of Maestro. Docking analysis was performed with the flexible ligand docking at an "extra precision" level. Docking scores were then used to obtain the predicted inhibition constants (K i ) via the formula: K i = exp (ΔG/RT) , where ΔG is the binding energy, R is the universal gas constant (1.98 cal mol − 1 K − 1 ) and T is the temperature (298.15 K). Finally, RMSD values were determined to compare the docked conformation of ligands with the nelfinavir conformation as a positive control.

Metabolite profiling of S. officinalis and A. dracunculus was d et er m i ne d b y U H P L C -H R M S / MS . T o pe rf o rm dereplication, the HPLC analytical conditions were initially optimized for both extracts before LC-MS and MS/MS analyses. Acquired MS/MS spectra of both extracts in the positive mode were used to generate a network for the visualization process of dereplicated constituents in an optimum manner. Analysis of the MS/MS data of S. officinalis gave rise to the identification of 637 parent ions, which were visualized as nodes in a molecular network, forming 117 clusters ranging from one to 86 connected nodes (Fig. 1) . Then, analysis of the generated network against the GNPS library resulted in annotation of thirty-seven constituents (Table 1) The most abundant compounds of this network were included in the flavonoids. Eleven glycosylated flavonoids, six un-glycosylated flavonoids, and two glycosylated anthocyanins, members of the flavonoid group of phytochemicals, were annotated. Phenyl propanoids (7 Compd.) were also identified to be the other major class of dereplicated compounds ( Fig. 2 ; Table 1 ).

The molecular networking evaluation revealed 538 parent ions for the MS/MS data of A. dracunculus, forming 61 clusters ranging from one to 43 connected nodes (Fig. 1) . Molecular networking analysis of the network of A. dracunculus against the GNPS library afforded to the tentative annotation of 20 compounds (Table 2) . Flavonoids (7 compounds), coumarins (6 compounds) phenylpropanoids (5 compounds), and dihydrochalcones (2 compounds) were the class of dereplicated compounds. The structure of dereplicated compounds has been shown in Fig. 3 .

Docking analysis was done to assess the binding energies and interaction modes of the dereplicated ligands from S. officinalis and A. dracunculus, and Z. officinale (obtained in our previous study (Babaeekhou and Ghane 2020) , against target proteins 6VW1, 6VSB (COV-2-SP), 6LU7, and 6M03 (M pro ) using Schrodinger package. The output of the docking analysis is shown in Table 3 .

To validate docking analyses, a standard ligand of each enzyme in co-crystallized complexes e.g. N3 peptide inhibitor from the M pro was removed and re-docked into the active site of the related enzyme. The same protocol such as the grid parameters and the precision level were employed in the process. It was performed to ensure the inhibitor binds exactly to the active site cleft. The superimposed analysis revealed less deviation of the re-docked complex in comparison to the actual co-crystallized complex (Table 4 ). In addition, Fig. 4 indicates the superimposed form of all docked ligands in the active site of four evaluated targets.

Based on the obtained results, the highest binding energy with COV-2-SP (6VW1) was observed for kaempferol 3-O-rutinoside (SO-36) and neodiosmin (SO-37) compounds from S. officinalis extract, followed by chicoric acid (AD-17), querciturone (AD-18) from A. dracunculus and 3-[4,5-dihydroxy-6-(hydroxymethyl)-3-[3,4,5-trihydroxyoxan-2-yl] oxyoxan-2-yl] oxy − 2-(3,4-dihydroxyphenyl)-5-hydroxy-7-methoxy chromen-4-one (ZO-13) and sissostrin The major interactions between the high potentially active ligands and the active site of proteins were observed to be the H-bond interaction of the hydroxyl groups and ion-bonds as well as π-π stacking of the aromatic rings with the enzyme (Table 5) . Regarding 6VSB which is consists of three chains and 1288 residues, a high number of H-bonds is observed between ligands and amino acids. Furthermore, binding interaction profiles of some active ligands against the M pro and COV-2-SP targets are indicated in Figs. 5 and 6.

Several studies have reported the biochemical composition of S. officinalis extract. Lima et al. (2007) evaluated the plant water and methanolic extracts for the presence of phenolic compounds using HPLC/DAD and detected 5 phenolic acids and 3 flavonoids. In a similar study, using HPLC-UV/VIS polyphenolic profile of the S. officinalis was identified for 14 compounds, and analysis of phytochemical compounds of the plant by HPLC-DAD-MSD revealed different phytochemical compounds of the plant including phenolic acids (9 compounds), flavonoids (3 compounds), phytosterols (2 compounds), saponins (6 compounds), and alkaloids (5 compounds) (Hernández- Saavedra et al. 2016 ). In the above-mentioned studies, rosmarinic acid, caffeic acid, luteolin-7-glucoside, chlorogenic acid, epicatechin, ellagic acid, quercetin were the major identified compounds. Using MN technique in the present study a wide range of components including polyphenols, alkaloids, flavonoids, terpenoids, anthraquinones, glycosides, and steroids were dereplicated and identification of some common compounds like caffeic acid, rosmarinic acid, luteolin, and quercetin in addition to 33 other compounds showed MN is an applicable and sensitive technique for chemical composition identification. The same scenario applies to the A. dracunculus and some similar compounds to the present study have been reported from previous studies (Duric et al. 2015; Mumivand et al. 2017) which confirms the accuracy of the applied technique. In the Mumivand et al. study major phenolic and flavonoid compounds of 12 Iranian A. dracunculus extracts were identified by RP-HPLC which resulted in the detection of chlorogenic, syringic, and caffeic acids and the predominant flavonoid was quercetin. It is worth mentioning that there are some identified compounds for S. officinalis and A. dracunculus which are reported for the first time in this study (Tables 1 and 2) . By targeting viral proteins in in-silico studies, we can rapidly screen plant compounds and make a shortlist of drug candidates. This helps to meet the immediate demand for an effective treatment against the 2019-nCOV infection. Therefore, in this study, we have performed a docking analysis on 104 yielded plant compounds to reach a list of potential candidates for future in vitro/vivo investigations. In this study, COV-2-SP and M pro viral potential targets were chosen for the assessment of affinity of 3 plant compounds including Z. officinale, reported in our previous study (Babaeekhou and Ghane 2020) (Table 3) .

Extensive research on the cell entry mechanisms of coronaviruses shows that different domains of COV-2-SP are participating in receptor binding and fusion. It is suggested that the receptor-binding domain plays a binary role in coronavirus entry including viral attachment to the host receptor and the fusion of the viral envelope with the cellular membrane . In this study, based on the proposed detailed structure of mouse hepatitis coronavirus (MHV) -12 -4.126 -5.269 -4.41 -5.561 spike protein and the amino acid numbering for spike protein domains Walls et al. 2016; Walls et al. 2017) , the major interactions between active ligands (high docking score) and the spike protein amino acids were assessed. 4 compounds from Z. officinale (ZO-13, ZO-31, ZO-40, ZO-41), 10 from S. officinalis and 7 from A. dracunculus (AD-13 to AD-15, AD-17 to AD-20) with high docking scores formed H-bonds, pi-pi interactions, and ion-bonds with amino acids 369-379 of 6VW1 (located on RBD of spike protein) ( Table 5 ; Fig. 5 ). So the inhibitory action of the above-mentioned compounds against the binding of COV-2-SP to its receptor can be proposed in this study. Analysis of the ligand interaction with 6VSB A, B, and C chains showed that 6 compounds from Z. officinale , 13 from S. officinalis , and 7 compounds from A. dracunculus (AD-14 to AD-20) have made H-bonds, pi-pi interactions and ion-bonds with residues located in the RBD (326-572), fusion peptide (amino acids 864-947, buried inside the pre-fusion structure) and S2ʹ cleavage site of the spike protein (Shang, Ye, et al. 2020; Walls et al. 2016; Walls et al. 2017) (Table 5 ; Fig. 5 ). So, we can assume that these plant compounds have the potential to inhibit COV-2-SP attachment, proteolysis of S2ʹsite, the transition of the spike protein to the post-fusion conformation, and the fusion of the virions. Antiviral activity of Artemisia annua against SARS-CoV is shown by Li et al. (2005) . In this study ethanol extracts of the plant could inhibit Vero E6 cells infection by Two strains of SARS-CoV (BJ001, BJ006). Also, Dihydrotanshinone, a lipophilic compound from Salvia miltiorrhiza is shown to have inhibitory activity toward the Middle East respiratory syndrome coronavirus (MERS-CoV) viral entry (Kim et al. 2018 ). There are other antiviral plant compounds that can affect spike protein. Emodin, for instance, can inhibit the binding of the protein to its receptor ACE2 and block viral penetration into cells (Ho et al. 2007; Schwarz et al. 2011) . In this present study, sissostrin (ZO-40), luteolin -5.89, -7.013/-8.216, -9.517/-9.43, -10.575/-8.943, -10.208/-9.144, -7.279/-5.975, -8.396/ -6.994, -9 .457/-9018, and − 9.256/-9.904 kcal/mol exhibited high binding affinity with both or one of studied targets for spike protein (6VW1/6VSB) and are introduced as potential viral entry inhibitors. In an in-silico study by Pandey et al. (2020) , SARS-CoV-2 spike protein with PDB-ID: 6VYB was targeted by some flavonoids and non-flavonoids compounds and similar to our results luteolin, curcumin, and quercetin could bind to S2 Domain and C-terminal of S1 domain with energies of -8.2, -7.1, and − 8.5 kcal/mol respectively. The binding energies for Kaempferol 3-O-rutinoside (SO-36) in our study were higher than the value for Kaempferol reported in Pandey et al. (2020) study (-10.575/-8.943 versus − 7.4 kcal/mol). Kiran et al. (2020) also showed high LFrank scores for luteolin from Kabasura Kudineer and quercetin from JACOM (a novel herbal formulation) and 6VSB using Cresset Flare software. Our results from iso-quercitrin (SO-31) and 6VSB with a binding energy of -9.43 kcal/mol were comparable with the result of Hiremath et al. (2021) from Phyllanthus amarus (-8.60 kcal/mol with quercitrin). In an experimental study by Yi et al. (2004) it is shown that luteolin and quercetin have inhibitory activity against SARS-CoV through interference on the fusion process and entry into host cells. Luteolin also is shown to have a high affinity to the S2 subunit of the spike protein in SARS-CoV (IC50: 10.6 µM) and it is postulating that this component can interfere with the viral cell fusion process (Wu et al. 2004 ). Docking analysis results for COVID-19 M pro structure (6LU7 and 6M03) showed more than 15 components have high binding affinity to the target structures including 3-[(2 S,3R,4 S,5 S,6R)_4,5-dihydroxy-6-(hydroxymethyl)_3-[(2 S,3R,4 S,5R)_3,4,5-rihydroxyoxan-2-yl]oxyoxan-2-yl] o x y -2 -( 3 , 4 -d i h y d r o x y p h e n y l ) _ 5 -h y d r o x y -7methoxychromen-4-one (ZO-13), [5-acetyloxy-1,7-bis (3,4-dihydroxyphenyl)heptan-3-yl] acetate (ZO-25), (2E) _5-(3-Acetoxy-6-hydroxy-5,5,8a-trimethyl-2-methylene-4 oxodeca hydro-1-naphthalenyl)_3-methyl-2-pentene-1,4diyl diacetate (ZO-31), sissostrin (ZO-40), and curcumin PE (ZO-43) from Z. officinale; SO-18, SO-19, SO-27, SO-29, SO-31, SO-35, SO-36, and SO-37 from S. officinalis, and AD-14 to AD-20 from A. dracunculus (Table 3) .

A motif between domains I and II (between amino acids 8-183) of M pro contains the substrate-binding site of the enzyme and residues 200 to 300 participate in the proteolytic activity of M pro (Anand et al. 2003) . In this study, it is shown that mentioned components with high docking scores have formed H-bonds and pi-pi interactions with the amino acids 2019; Nile et al. 2020; Tajik et al. 2017; Wu et al. 2018; Li et al. 2016) . So, the different combined usage of the reported compounds with dual interaction with COV-2-SP and M pro can be a better strategy for entry and viral replication interruption of the 2019-nCOV.

The molecular networking technique is an accurate tool for the differentiation of the metabolites in plant extracts. Docking analysis results of some compounds in the present study 

ARG-105

ASN-A: 388, ASP-B: 389, ASN-B: 542, GLY-C: 545 (2 Int

ZO-25 … THR-111, GLN-110

PHE-C: 543, PHE-C: 565, THR-B: 998, THR-C: 998 ZO-31 HIS-41

GLN-110, THR-111, ILE-249

GLN-A: 580, THR-B: 998, THR-C: 998 SER-371 (2 Int.), ARG-408

HIE-164, GLU-166

GLN-110 (2 Int.), THR-111, ASN-151, ILE-152

SER-371

ZO-43 HIS-41, LEU-141, ASN-142, GLU-166

MET-A: 731, LYS-A: 733, ARG-A: 815, PHE-A: 823, LEU-A: 861, ASP-A: 867, HIS-A: 1058, TYR C: 756, PHE A:970, GLN B:1002 SO-18 THR-26

ARG-015

ARG-B: 995

SO-19 HIS-41, CYS-44, SER-46, GLY-143

GLN-110

PRO-A: 1057, SO-20 ASP-C: 389 (2 Int.), THR-A: 523 (2 Int

PRO-A: 527, LYS-A: 529 (2 Int.), PHE-A: 970 (2 Int.), ARG-B: 995

CYS-379 SO-28 PHE-A: 970, PHE C: 970, ASP-A: 994, ASP-B: 994, ARG-C: 995, THR-B: 998

SO-29 HIS-41, GLU-166

ASN-C: 544, ASN-B: 542, MET-A: 731, LYA-A: 733, ARG-A: 815 (2 Int.), VAL-A: 826, LEU-A: 861, ASP-A: 867, ASP-A: 994, ARG-A: 995, THR-C: 998

LYS-378, CYS-379

PHE-A: 970, PHE C: 970, ASP-A: 994, ASP-B: 994, ARG-A: 995, ARG-C: 995, THR-B: 998

GLU-166 (2 Int.)

PHE-A: 970, PHE-C: 970, ASP-A: 994, ASP-B: 994, THR-998 B

GLN-110 triple, ASP-153

HIS-41, ASN-119, LEU-141, GLY-143, GLU-166

GLN-110, THR-111, SER-158

ASP-B: 994, ASP-C: 994 TYR-369

SER-46, LEU-141, ASN-142, GLY-143

ASP-B: 389, PRO-A: 527 (2 Int

GLU-166 (2 Int.)

ASP-C: 571 (2 Int.), LYS-C: 733, ARG-C: 815 (2 Int.), LEU-C: 861, THR-C: 866, ASP-C: 867, ASP-C: 994, ARG-A: 995, ARG-B: 995

AD-13 SER-371

AD-14 THR-25, SER-46, GLY-143, GLU-166

LYS-B: 386 (2 Int

ASN-B: 544, PRO-A: 527 (2 Int

AD-16 THR-26, ASN-142, GLY-143, HIE-164, GLU-166

GLN-110

PHE-A: 970, ASP-A: 994, ASP-B: 994, ARG-C: 995, THR-B: 998 AD-17 HIS-41

GLN-110, THR-111, ASP-153, ILE-249

ASN-B: 542, GLY-B: 545, MET-A: 731, ASP-A: 775, ARG-A: 815, PHE-A: 823, ASP-A: 867, ASP-A:994, THR-A:998, THR-B: 998, THR-C:998 TYR-369

AD-18 PHE-140, ASN-142, GLY-143

GLN-110

MET-A: 731, PHE-C:970, ASP-A: 867, ASP-994-A:994, ASP-B:994, THR-B: 998, ARG-995, HIS-A: 1058

S,3R,4 S,5 S,6R)_4,5-dihydroxy-6 -(hyy-2-( 3 ,4 -dihydroxyphenyl)

SER-144, LEU-141, CYS-145, GLU-166, and GLN-189 residues of of SARS-CoV 3CL pro and kaempferol was shown. Binding energies of luteolin (ZO-41) and 6-Methoxy luteolin (ZO-41), isoquercitrin (SO-31), quercetin (AD-13), quercetagetin (AD-14), querciturone (AD-18), kaempferol (AD-11)

Artemisia dracunculus; ACE-2, Cell surface receptor angiotensin-converting enzyme 2

Coronavirus disease 2019; M pro , Main protease; MN, Molecular Networking; RBD, Receptor-binding domain; SO

Severe Acute Respiratory Syndrome-related Coronavirus; ZO, Zingiber officinale; 2019-nCOV

C-like-protease Supplementary Information The online version contains supplementary material

Quercetagetin (AD-14) (b), kaempferol 3-O-rutinoside (SO-36

Structure of coronavirus main proteinase reveals combination of a chymotrypsin fold with an extra α-helical domain

Coronavirus main proteinase (3CLpro) structure: basis for design of anti-SARS drugs

Antimicrobial activity of ginger on cariogenic bacteria: molecular networking and molecular docking analyses

The coronavirus spike protein is a class I virus fusion protein: structural and functional characterization of the fusion core complex

Coronavirus cell entry occurs through the endo-/lysosomal pathway in a proteolysis-dependent manner

A review of the dietary flavonoid, kaempferol on human health and cancer chemoprevention

Binding interaction of quercetin-3-β-galactoside and its synthetic derivatives with SARS-CoV 3CLpro: Structure-activity relationship studies reveal salient pharmacophore features

Virus entry: molecular mechanisms and biomedical applications

The spike protein of SARS-CoV a target for vaccine and therapeutic development

Anticoagulant activity of some Artemisia dracunculus leaf extracts

Remdesivir for 5 or 10 days in patients with severe Covid-19

Phytochemical characterization and effect of Calendula officinalis, Hypericum perforatum, and Salvia officinalis infusions on obesity-associated cardiovascular risk

Nucleotide sequence of the human coronavirus 229E RNA polymerase locus

In silico docking analysis revealed the potential of phytochemicals present in Phyllanthus amarus and Andrographis paniculata, used in Ayurveda medicine in inhibiting SARS-CoV-2. 3 Biotech 11:44

Emodin blocks the SARS coronavirus spike protein and angiotensin-converting enzyme 2 interaction

Chicoric acid prevents methotrexate hepatotoxicity via attenuation of oxidative stress and inflammation and up-regulation of PPARγ and Nrf2/HO-1 signaling

Structure of M pro from SARS-CoV-2 and discovery of its inhibitors

Inhibition of SARS-CoV 3CL protease by flavonoids

Safe, highthroughput screening of natural compounds of MERS-CoV entry inhibitors using a pseudovirus expressing MERS-CoV spike protein

In silico computational screening of Kabasura Kudineer-official Siddha Formulation and JACOM against SARS-CoV-2 spike protein

Curcumin and its derivatives: their application in neuropharmacology and neuroscience in the 21st century

Structure, function, and evolution of coronavirus spike proteins

Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus

Identification of natural compounds with antiviral activities against SARS-associated coronavirus

Comparison of the antioxidant effects of quercitrin and isoquercitrin: Understanding the role of the 6 ″-OH group

Water and methanolic extracts of Salvia officinalis protect HepG2 cells from t-BHP induced oxidative damage

Evaluation of the virucidal effects of rosmarinic acid against enterovirus 71 infection via in vitro and in vivo study

Antioxidant properties and principal phenolic phytochemicals of Iranian tarragon

Flavonoid-mediated inhibition of SARS coronavirus 3C-like protease expressed in Pichia pastoris

Probing the effect of quercetin 3-glucoside from Dianthus superbus L against influenza virus infection In vitro and in silico biochemical and toxicological screening

Targeting SARS-CoV-2 spike protein of COVID-19 with naturally occurring phytochemicals: an in-silico study for drug development

Evaluation of polyphenols from Broussonetia papyrifera as coronavirus protease inhibitors

Biflavonoids from Torreya nucifera displaying SARS-CoV 3CLpro inhibition

Emodin inhibits current through SARS-associated coronavirus 3a protein

Cryoelectron microscopy structure of porcine delta-coronavirus spike protein in the prefusion state

Structure of mouse coronavirus spike protein complexed with receptor reveals mechanism for viral entry

Structural basis of receptor recognition by SARS-CoV-2

Discovery of potent anti-SARS-CoV Mpro inhibitors

Cryo-EM structure of the SARS coronavirus spike glycoprotein in complex with its host cell receptor ACE2

The potential effects of chlorogenic acid, the main phenolic components in coffee, on health: a comprehensive review of the literature

Viral replicase gene products suffice for coronavirus discontinuous transcription

Bioavailabilityenhanced Resveramax TM modulates quorum sensing and inhibits biofilm formation in Pseudomonas aeruginosa PAO1

Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: an interim analysis of four randomized controlled trials in Brazil, South Africa, and the UK

Cryo-electron microscopy structure of a coronavirus spike glycoprotein trimer

Tectonic conformational changes of a coronavirus spike glycoprotein promote membrane fusion

Specific plant terpenoids and lignoids possess potent antiviral activities against severe acute respiratory syndrome coronavirus

Weekly epidemiological update on COVID-19 -29 June

Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation

Small molecules targeting severe acute respiratory syndrome human coronavirus

Differential effects of quercetin and two of its derivatives, isorhamnetin and isorhamnetin-3-glucuronide, in inhibiting the proliferation of human breast cancer MCF-7 cells

Small molecules blocking the entry of severe acute respiratory syndrome coronavirus into host cells

Computational screening of antagonist against the SARS-CoV-2 (COVID-19) coronavirus by molecular docking

In silico screening of Chinese herbal medicines with the potential to directly inhibit 2019 novel coronavirus

A novel coronavirus from patients with pneumonia in China

Characterization of a human coronavirus (strain 229E) 3C-like proteinase activity

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Acknowledgements The authors would like to acknowledge Miss Asal Attari for editing this manuscript for the English language.

Conflict of interest The authors have no conflicts of interest to declare that are relevant to the content of this article.

This observation is in line with a recent study by Yu et al. (2020) which introduces luteolin as a potent compound against COVID-19. The binding affinity between quercetin (AD-13) from A. dracunculus and M pro (6M03) in this study was − 6.981 kcal/mol which was near to the value reported for this compound in Zhang et al. study (2020) (− 6.25 kcal/mol) ) but higher binding affinities are reported in this study for quercetin derivatives like quercetagetin (AD-14) and querciturone (AD-18) with binding energy values of -8.11 and − 9.53 kcal/mol respectively. In the mentioned study kaempferol was introduced as a candidate which can inhibit 3CLpro with a binding affinity of -6.01 kcal/mol . Higher binding affinity was observed in this study between kaempferol 3-O-rutinoside (SO-36) from S. officinalis and M pro (-9.705). Hiremath et al. (2021) have shown binding affinity of some flavonoids from Phyllanthus amarus using AutoDock Vina software and 3CLpro (6LU7) and their results showed accordance with our results from kaempferol (AD-11) and quercetin (AD-13). They have reported − 7.70 and − 7.50 Kcal/mol for mentioned compounds respectively and our obtained figures were − 6.88 and − 6.98 (Hiremath et al. 2021) .According to the reports from experimental studies, quercetin has shown inhibition activity toward 3CL pro with an IC 50 value of 73.7 µM (Nguyen et al. 2012) . Ryu et al. (2010) exhibited the protease inhibitory of luteolin, and quercetin with (IC 50 : 20.2, and 23.8 µM, respectively). Quercetin has also captured an IC 50 of 52.7 µM in Park et al. (2017) study. In the latter, IC 50 of 5.7 and 116.3 µM was reported for curcumin and kaempferol respectively and a prenylated quercetin derivative displayed the most inhibitory activity (IC 50 : 3.7). Curcumin was also reported as a 3Cl pro inhibitor in Vero E6 cells (Wen et al. 2007) .In three plants used in this study, 23 components showed high docking scores with both COV-2-SP and M pro viral targets (Table 3) , among which some flavonoids and phenylpropanoids such as curcumin PE from Z. officinale; kaempferol 3-O-rutinoside, isoquercitrin, and rosmarinic acid from S. officinalis; and chlorogenic acid, chicoric acid, querciturone and isorhamnetin 3-glucuronide from A. dracunculus have shown different therapeutics or antimicrobial properties in medicine or microbiology (Chen and Chen 2013; Hussein et al. 2020; Lee et al. 2013; Lin et al.