key: cord-0824019-27l14n9s
authors: Das, Sanjib Kumar; Mahanta, Saurov; Tanti, Bhaben; Tag, Hui; Hui, Pallabi Kalita
title: Identification of phytocompounds from Houttuynia cordata Thunb. as potential inhibitors for SARS-CoV-2 replication proteins through GC–MS/LC–MS characterization, molecular docking and molecular dynamics simulation
date: 2021-05-07
journal: Mol Divers
DOI: 10.1007/s11030-021-10226-2
sha: 046b3ad0ca5ae02dd3c0bd6cce39b5db48108426
doc_id: 824019
cord_uid: 27l14n9s

The COVID-19 pandemic caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a massive viral disease outbreak of international concerns. The present study is mainly intended to identify the bioactive phytocompounds from traditional antiviral herb Houttuynia cordata Thunb. as potential inhibitors for three main replication proteins of SARS-CoV-2, namely Main protease (Mpro), Papain-Like protease (PLpro) and ADP ribose phosphatase (ADRP) which control the replication process. A total of 177 phytocompounds were characterized from H. cordata using GC–MS/LC–MS and they were docked against three SARS-CoV-2 proteins (receptors), namely Mpro, PLpro and ADRP using Epic, LigPrep and Glide module of Schrödinger suite 2020-3. During docking studies, phytocompounds (ligand) 6-Hydroxyondansetron (A104) have demonstrated strong binding affinity toward receptors Mpro (PDB ID 6LU7) and PLpro (PDB ID 7JRN) with G-score of − 7.274 and − 5.672, respectively, while Quercitrin (A166) also showed strong binding affinity toward ADRP (PDB ID 6W02) with G-score -6.788. Molecular Dynamics Simulation (MDS) performed using Desmond module of Schrödinger suite 2020–3 has demonstrated better stability in the ligand–receptor complexes A104-6LU7 and A166-6W02 within 100 ns than the A104-7JRN complex. The ADME-Tox study performed using SwissADMEserver for pharmacokinetics of the selected phytocompounds 6-Hydroxyondansetron (A104) and Quercitrin (A166) demonstrated that 6-Hydroxyondansetron passes all the required drug discovery rules which can potentially inhibit Mpro and PLpro of SARS-CoV-2 without causing toxicity while Quercitrin demonstrated less drug-like properties but also demonstrated as potential inhibitor for ADRP. Present findings confer opportunities for 6-Hydroxyondansetron and Quercitrin to be developed as new therapeutic drug against COVID-19. [Image: see text]

Sudden outbreak of COVID-19 has become life threatening for millions of human population across the globe. World Health Organization (WHO) declared COVID-19 as the 6th public health emergency of international concern on 30th January 2020. WHO reported 760.2 million active cases of COVID-19 followed by 16.6 millions of confirmed deaths globally till the month of December, 2020 [1] . The COVID-19 is fundamentally caused by the infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the human host and transmits to another human being through physical contact or droplets. Although antiviral drug and vaccine development are being undertaken for the prevention of the disease, presently, there are no effective therapeutic drug molecules available to date against COVID-19. SARS-CoV-2 is an enveloped virus that belongs to family of Coronaviridae containing single positive-stranded RNA and replicate inside the cytoplasm of human cell [2, 3] . SARS-CoV-2 genome encodes several proteins that are involved in the replication of viral genome. Inhibition of the expression of these receptor proteins may lead to the possibilities of anti-SARS-CoV-2 drug development for effective treatment of COVID-19. Three proteins of SARS-CoV-2, namely Main protease (Mpro), Papain-Like protease (PLpro) and ADP ribose phosphatase (ADRP), are reported to be mainly responsible for SARS-CoV-2 replication process [4] . The Main protease (Mpro) is present in SARS-CoV-2 and other coronaviruses (CoVs) which cleave two replicase polyprotein responsible for viral replication and maturation [5] . The papain-like protease (PLpro) is another proteolytic enzyme which is essential for processing of SARS-CoV-2 polyproteins to generate a functional replicase complex and initiate viral spread. PLpro is also involved in cleaving protein during post-translational modifications on host proteins as an eluding mechanism against host antiviral immune responses. PLpro cleaves ubiquitin and ISG15 which are known regulators of host innate immune pathways [6] . On the other hand, ADP ribose phosphatase (ADRP) is known enzymes of SARS-CoV-2 replication that convert ADP-ribose 1″-monophosphate (Appr-1″-p) to ADP-ribose (Appr), which regulate replication of the virus [7] . Therefore, inhibition of Mpro, PLpro and ADRP activities through bioactive phytochemical drug molecules would prevent replication mechanism of SARS-CoV-2.

In silico method employed for the screening of potential drug molecules is proven rapid and cost effective while compared with the trial and error methods using experimental studies. Molecular docking is an in silico approach with ability to screen and identify potential drug molecules from large and massive bioactive compound libraries. Currently, several molecular docking studies were carried out against SARS-CoV-2 receptors with some selected bioactive compounds (Table 1 ) and have been able to identify few potential compounds drug molecules effective against targeted receptors of COVID-19 [7, [9] [10] [11] [12] [13] [14] . However, most of these docking studies relied on structure-based drug design (SBDD), as well as similarity searching and quantitative structure-activity relationship (QSAR) modeling [8] . Some recent cases of molecular docking studies of different molecular targets for SARS-CoV-2 are presented in Table 1 .

The traditional healer of different ethnic communities used medicinal plants for treatment and prevention of various diseases including viral infections. Isolation of these bioactive phytocompounds from traditional medicinal plants may lead to the development of antiviral drug for effective treatment of COVID-19 [15] . Houttuynia cordata Thunb. is a herb belonging to the family Saururaceae, distributed in China, India, Bangladesh, Cambodia, East Himalaya, Hainan, Japan, Korea, Myanmar, Nepal, Taiwan, Thailand, Tibet (China), Vietnam, West Himalaya [16] . The whole plant is used in folk medicine for the treatment of cough, pneumonia, bronchitis, dysentery, dropsy, leukorrhea, uteritis, eczema, herpes simplex, acne, chronic sinusitis, stomach ulcer, infection, control wrinkle, chapped skin, septic, febrifuge, heatstroke, malaria, lung disorder, tonsillitis, skin ulcer, diarrhea, dysentery arthritis, appendicitis, and snake bite in Japan, Korea, Indonesia and Myanmar [17] [18] [19] . In Indian sub-continent, H. cordata is mainly used for the treatment and prevention of pneumonia, stomach disorder, sinusitis and heart disorders by the indigenous tribes of North East Hesperidin [14] India [20] [21] [22] . The major bioactive phytocompounds including flavonoid, alkaloid, as well as essential oil, have been reported from various parts of H. cordata which possesses antimicrobial, hepatoprotective, anti-inflammatory, antioxidative, anticancer, anti-diabetic, anti-obesity and antiviral properties [17] . Houttuynia cordata was one of the most important ingredient of the herbal formulation used for the treatment of Severe Acute Respiratory Syndrome (SARS) outbreak of Southern China in 2003 [16] . Several scientific studies also reveal that the extract of H. cordata has the potential to inhibit replication of the various viral strains including SARS coronavirus (SARS-CoV), Chikungunya, Herpes simplex viruses, dengue virus serotype 2 (DEN-2), Influenza neuraminidase, pseudorabies herpes virus (PrV), Human noroviruses (HuNoVs), murine coronavirus and dengue virus infection and also possess innate immune modulation activities [16, 18, [23] [24] [25] [26] [27] [28] [29] . Therefore, considering the therapeutic importance of H. cordata with strong ethnopharmacological use background, the present study is mainly intended to characterize the bioactive phytocompounds by using gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) studies, and to perform molecular docking studies of bioactive phytocompounds against the three main proteins (target receptors) of SARS-CoV-2, namely Mpro, PLpro and ADRP, responsible for replication process of the virus, and to identify the bioactive phytocompounds (ligands) potential to inhibit SARS-CoV-2 replication proteins by determining the binding affinities of the ligands against selected receptors. We performed molecular docking study between the active phytochemicals of H. cordata using Epic, LigPrep and Glide module of Schrödinger suite 2020-3 followed by Molecular Dynamics simulation by Desmond module of Schrödinger suite 2020-3. We also performed Adsorption, Digestion, Metabolism, Excretion and Toxicity (ADME-Tox) study of the shortlisted phytocompounds using SwissADMEserver to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of the selected bioactive phytocompounds.

The herb H. cordata was collected on 13 November 2019 from Ziro valley of Lower Subansiri District of Arunachal Pradesh, India (26°55′-28°21′ N and 92°40′-94°21′ E). The plant materials were identified and authenticated at ASSAM Herbarium, Botanical Survey of India, Eastern Regional Center, Shillong, Meghalaya, India. The accepted names and global distribution ranges were verified in http:// www. plant softh eworl donli ne. org (POWO) hosted by Royal Botanic Garden, Kew, UK. The voucher specimen No. SD/ HAU-0110 dated 13.11.2019 was deposited in the Herbarium of Arunachal University (HAU), Department of Botany, Rajiv Gandhi University, Rono Hills, Doimukh, Arunachal Pradesh, India for future references. The fresh herb harvested was washed in running tap water to remove microflora, dirt and soil particles and then shade dried at normal room temperature (20.0-27.0 °C) with proper ventilation. The minimum moisture content of the dried samples was maintained at 15% using Material Moisture Meters (Model Testo 606-1, Germany). The powdered samples were passed through 40-mesh sieve and stored for further use.

The crude extract was prepared by cold maceration of powder materials with methanol (solvent ratio 1:5). 50 g of powder sample was soaked in 250 ml of methanol and placed on magnetic stirrer for 48 h [30] . After 48 h, filtrate was 

The GC-MS characterization of crude extract of H. cordata was performed by injecting 2 µl of sample in a GC-MS-2010 Shimadzu instrument operating in EI mode at 70 eV with Restek-5MS column (30 × 0.25 mm film thickness 0.25 μm). The oven temperature was programmed as follows: kept at 60 °C for 2 min, then increased to 210 °C, at 3 °C min −1 and from 210-280 °C at the rate of 8 °C/min and held for 8 min at 280 °C. The injector temperature was 260 °C with normal injection mode. The flow rate of carrier helium gas was 1.21 ml/min. For the plant samples, the temperature program was set according to Jiang et al. [31] . Mass spectra were acquired using full scan monitoring mode with a mass scan range of 40-650 m/z. The chromatogram and mass spectra were evaluated using the Xcalibur™ software embedded in the GC-MS/LC-MS system. Interpretation on mass spectrum GC-MS was conducted using the database of National Institute Standard and Technology (NIST), Wiley having more than 62,000 patterns and Dr. Duke's phytochemical ethnobotanical databases (Duke 2018; phyto chem. nal. usda. gov). The spectrum of the unknown component was compared with the spectrum of the known components stored in the NIST library. The name, molecular weight and structure of the compounds of the test materials were ascertained.

The LC-MS characterization of the H. cordata was performed by injecting 10 µl of sample in a LC-MS Thermo Scientific Plus with Dionex Ultimate 3000 HPLC with column Hypersil Gold (C18) Diameter 150 × 2.1 mm, Particle Size 1.9 μ at room temperature. The mobile phases used were A (acetonitrile) and B (0.2% aqueous acetic acid, v/v). The run time was set for 20 min followed by flow rate of 0.6 mL/min. DAD detector was set at 280 nm for acquiring chromatograms. The mass spectrometer used was a Triple quadruple Mass Spectrometer (Thermo Scientific) equipped with ion sources ESI with Mass range for full scans m/z 50-6000. All raw data acquired were processed by METLIN database [32] .

The crystal structure of SARS-CoV-2 main protease in complex with an inhibitor N3 having PDB ID 6LU7 [33] (resolution 2.16 Å), wild type SARS-CoV-2 papain-like protease (PLPro) with inhibitor GRL0617 having PDB ID 7JRN (resolution 2.48 Å) and ADP ribose phosphatase of NSP3 from SARS-CoV-2 in the complex with ADP ribose having PDB ID 6W02 (resolution 1.5 Å) were downloaded from RCSB Protein Data Bank [34] and same were processed through protein preparation wizard of Epic module [35] of Schrödinger suite 2020-3. The proteins were prepared by removing similar binding sites, unnecessary water molecules and also refining bond orders. Missing chain atoms were added by using the prime module of Schrödinger suite 2020-3. Energy minimization of the proteins was performed using optimized potentials for liquid simulations-3 (OPLS3e) molecular force field with root-mean-square difference (RMSD) of crystallographic heavy atoms kept at 0.3 Å. The ligand binding information of N3, GRL0617 and ligand binding pocket of ADP ribose were used for prediction of the active site of the proteins.

The ligand library from the GC-MS and LC-MS analysis of H. cordata, having traditional medicinal linkage with antiviral treatment, was processed through the LigPrep module of Schrödinger suite 2020-3. During preparation of the ligands for docking, Lipinski Drug Discovery rule [36] was kept as filtering criteria to screen out unfavorable molecules from the dataset. 2D structures were converted to 3D structures, and optimized for their geometry, desalted and their chirality was corrected. The ionization and tautomeric states were generated between pH 6.8 and 7.2 by using Epik module of Schrödinger suite 2020-3. The ligands were minimized using OPLSe-3 force field in Schrödinger suite 2020-3 until a RMSD of 2.0 Å was achieved and the optimized ligands were used for docking analysis.

The filtered phytocompounds were prepared by LigPrep module of Schrödinger suite before proceeding for docking, and the docking was performed against each of the three selected receptor proteins, Mpro (PDB ID 6LU7), PLpro (PDB ID 7JRN) and ADRP (PDB ID 6W02) using Glide module of Schrödinger suite 2020-3. During the docking process, the receptor grid maps were prepared from the active site information available from the previously bound ligands with the receptor in crystal structure available from RCSB PDB. The binding modes with best docking score and best Glide gscores (G-Score) were selected. These scores perceive positive lipophilic, hydrogen bonding and metal-ligand associations and punish steric conflicts. The G-scoring capacity is mainly dependent on docking parameters like lipophilic perseverance in which the ligands are covered in the lipophilic pocket. The electrostatic forces and hydrogen bonding with ligands are other parameters to increase the binding affinity.

For Molecular Dynamics (MD) Simulation, the Desmond module through Maestro of Schrödinger suite 2020-3 was used in the Linux platform supported by 6 GB NVIDIA GeForce GTX graphics card. To study the stability of the docked complex with the ligand, an MD Simulation study of 100 nano second (ns) was performed [37] . The complex in the explicit solvent system with the OPLS3e force field was studied using the Desmond module of Schrödinger suite 2020-3. The molecular system was solvated with crystallographic water (TIP3P) molecules under orthorhombic periodic boundary conditions for a 10 Å buffer region. The overlapping water molecules are deleted, and the system was neutralized by adding Na + as counter ions. An ensemble (NPT) of Nose-Hoover thermostat [38] and barostat was applied to maintain the constant temperature (300 K) and pressure (1 bar) of the systems, respectively. A hybrid energy minimization algorithm with 1000 steps of steepest descent followed by conjugate gradient algorithms was utilized.

The structure and dynamic properties of the protein-ligand complexes were analyzed as the backbone RMSDs during the simulation period of 100 ns. The RMSD was measured as the average distance between the backbone atoms of the protein-ligand structures and it was derived from the following equation:

where N represents the total number of atoms considered in the calculation and ẟ represents the distance between the N pairs of equivalent atoms.

The root-mean-square fluctuations (RMSF) were assessed and plotted to equate the flexibility of each residue in the ligand-protein complexes. The RMSF of the protein-ligand complex denoted the minimized fluctuation for all the complexes.

Adsorption, Digestion, Metabolism, Excretion and Toxicity (ADME-Tox) study of the shortlisted compounds was performed using SwissADMEserver (http:// www. swiss adme. ch) [39] . 

GC-MS analysis of H. cordata has identified and quantified 80 phytocompounds and the result is summarized in Table 2 . The chromatograms of GC-MS result are given as Fig. 1 .

The non-volatile compounds were identified by LC-MS. Ninety-seven phytocompounds were identified by LC-MS. The result of LC-MS is summarized in Table 3 and chromatogram is presented as Fig. 2 .

The molecular docking of ligands (bioactive phytocompounds from H. cordata) was carried out against three different receptors of SARS-CoV-2, namely Mpro (PDB ID 6LU7), PLpro (PDB ID 7JRN) and ADRP (PDB ID 6W02), respectively, by using an advanced molecular docking program, Glide of Schrödinger suite 2020-3 to determine the binding affinities of the ligands against selected receptors. The binding affinities of the ligands toward receptors were determined based on Docking Score and Glide gscore. Lower docking score and Glide gscore indicate better affinity. The results of the molecular docking study revealed that the bioactive phytocompound coded with A107 from H. cordata showed highest binding affinity toward Mpro (6LU7) with G-Score − 7.929 ( Table 4 ). The bioactive phytocompounds A104, A120, A99 and A127 from H. cordata also showed significant binding affinities against Mpro (6LU7) with G-score − 7.274, − 7.024, − 6.883 and− 6.793, respectively (Table 4) . Similarly, bioactive phytocompounds A105 showed highest binding affinity against PLpro (7JRN) with G-Score − 5.843, and bioactive phytocompounds A140, A104, A100 and A106 showed potential binding affinities against PLpro (7JRN) with G-score − 5.779, − 5.672, − 5.602 and − 5.227, respectively (Table 4 ). ADRP (6W02) is another important protein that is required for SARS-CoV-2 replication against which the bioactive phytocompounds A166 showed highest binding affinity with G-Score (− 6.788). The phytocompounds A165, A163, A167, A164 showed binding affinity against ADRP (6W02) with G − Score − 6.099, − 5.703, − 5.703 and − 5.206, respectively (Table 4 ). From the present study, it is suggested that the bioactive phytocompound A104 is capable of binding both Mpro (6LU7) and PLpro (7JRN) which has demonstrated multi-targets ability and have been identified as a viable drug molecule potential to inhibit replication proteins Mpro (6LU7) and PLpro (7JRN) of SARS-CoV-2. Although A107 and A105/A140 have the highest predicted affinity toward Mpro (6LU7) and PLPro (7JRN), the phytocompound A104 was selected for further analysis because of its ability to bind both the receptors. Binding affinities with formation of H-bonds and other non-bonding interactions of the complexes are presented in Fig. 3a, b. On the other hand, the best binding ligand A166 was selected for further analysis against the target protein ADRP and the interactions of A166 with ADRP (6W02) are presented in Fig. 3c .

The protein-ligand complexes with the shortlisted ligands A104 and A166 were performed by MD simulation study to examine the stability of the compounds to the binding site of 6LU7, 7JRN and 6W02 according to their predicted affinity as mentioned in Table 4 . MD simulation deals with studying the behavior of protein and ligand for a particular time. The whole simulation was subjected for 100 ns in the production phase for the ligand complexes. The structure and dynamic properties of the protein-ligand complexes were analyzed as the backbone RMSDs during the simulation period of 100 ns. The RMSD plot for 6LU7 and 6LU7 bounded with A104 presented in Fig. 4a reveals that during the initial 60 ns of MD simulation, the ligand-bound protein showed significant stability, however, it was slightly unstable after 60 ns of simulation. Similarly, RMSD plot of 7JRN and 7JRN bonded with A104 is presented in Fig. 4b , where the ligand-bound complex established a stabilization between 50 and 70 ns of simulation. However, the RMSD plot 6W02 bonded with A166 presented in Fig. 4c showed instability from 0 to 19 ns of simulation, and after 20 ns, it became stabilized up to 30 ns. Although a slight deviation was observed between 30 and 60 ns, the complexes stabilized after 60 ns. Ligand RMSD, the radius of Gyration (rGyr), Molecular Surface Area (MolSA), Solvent Accessible Surface Area (SASA), Polar Surface Area (PSA) of ligands with respect to the reference conformation were also studied during the MD simulation which is presented in Fig. 5a-c. The root-mean-square fluctuations (RMSF) were assessed and plotted to equate the flexibility of each residue in the ligand-protein complexes. The RMSF of the protein-ligand complex denoted the minimized fluctuation for all the complexes. The RMSF did not deviate much during the simulation period of 100 ns, and the average RMSF values were kept constant for all the complexes and the same is presented in Fig. 6a -c. The hydrogen bond interaction between the bioactive compounds and the active site residues of the proteins was also analyzed during the entire duration of simulation. For Bold words signify main compound's potential to inhibit SARS-CoV-2 replication proteins 

SwissADME server was used to study the ADME-Tox properties of A104 (6-Hydroxyondansetron) and A166 (Quercitrin), and the results are presented in Table 5 :

Bioinformatics tools and techniques are widely used in drug discovery sector. In silico-based drug design is an approach which plays an important role in modern drug discovery. It reduces the cost and takes less time to search out lead molecules. This structural-based drug designing study provides the probable drug compound for a particular disease which can be further confirmed by in vivo, as well as in vitro experimental studies [40, 41] . The outbreak of COVID-19 caused by SARS-CoV-2 has thrown millions of human lives out of gear within one year timescale. The virus has not only disrupted the global healthcare networks but it also severely affected the global economy [42] . Therefore, there is an urgent need of antiviral drugs from bioactive phytocompounds for the inhibition of SARS-CoV-2 replication proteins. The SARS-CoV-2 is a single-stranded RNA virus which requires mainly three proteins, namely Main protease (Mpro), Papain-like protease (PLpro) and ADP ribose phosphatase (ADRP), for their replication. The inhibition of the expression of these proteins may lead to the development of potent antiviral effective against treatment of COVID-19. Bioactive phytocompounds from traditional medicinal plants are reported to possess several therapeutic properties. In the present study, we explored the possibilities and demonstrated that at least two bioactive phytocompounds, namely 6-Hydroxyondansetron (Fig. 9a) and Quercitrin (Fig. 9b) , quantified from Houttuynia cordata Thunb. possess potential anti-SARS-CoV-2 activities and demonstrated potential to inhibit Main protease (Mpro), Papain-like protease (PLpro) and ADP ribose phosphatase (ADRP). Previous study reported that the extract of H. cordata inactivates SARS-CoV 3CL Pro and inhibits the viral replication [9] . It also demonstrated antiviral activity against Chikungunya viral strains [18] , Herpes simplex viruses-1 (HSV-1) and Herpes simplex viruses-2 (HSV-2) [21] , and dengue virus serotype 2 (DEN-2) [24] . (6W02), respectively, revealed that the bioactive phytocompounds (ligand) coded with A104 (6-Hydroxyondansetron) have shown more binding affinities toward Mpro (6LU7) and PLpro (7JRN) with G-score − 7.274 and − 5.672, respectively, while bioactive phytocompound A166 (Quercitrin) has demonstrated the best binding affinity toward ADRP (6W02) with G-score − 6.788. Therefore, the molecular dynamics simulation study of both the bioactive phytocompounds A104 and A166 was performed. From the MD simulation analysis, it is evident that the complexes A104-6LU7 and A166-6W02 have shown better stability than the A104-7JRN complex within 100 ns duration. The interaction analysis during the entire duration of MDs, as presented in Fig. 7a-c and Fig. 8a -c, and its subsequent analysis also support good binding efficiency of the selected ligands to the receptor proteins. MDs study ultimately identified and confirmed A104 and A166 as potential inhibitors for Mpro, PLpro and ADRP which can prevent replication machinery of SARS-CoV-2. ADME-Tox study also reveals that the bioactive phytocompound A104 (6-Hydroxyondansetron) passes all the required drug discovery rules and possesses multi-target capabilities which can become a potential inhibitor for the Mpro (6LU7) and PLpro (7JRN) protein of SARS-CoV-2 without causing toxicity with minimum 

The present studies have characterized 177 phytocompounds from H. cordata through GC-MS/LC-MS and they were docked against three receptor proteins, namely Main protease (Mpro), Papain-like protease (PLpro) and ADP ribose as another promising inhibitor as it has shown best binding affinity toward protein ADRP (PDB ID 6W02) of SARS-CoV-2. Among the two-identified bioactive phytocompounds (ligands), 6-Hydroxyondansetron (A104) has demonstrated more promising properties as drug candidate since it has shown more binding affinity toward multiple receptors. This is very significant as chances of development of drug resistance are much lower in case of multi-target capabilities as inhibitor. The Molecular Dynamic Simulation analysis confirmed that the complex A104-6LU7 and A166-6W02 have shown better stability than the A104-7JRN complex within 100 ns duration. ADME-Tox study further demonstrated that out of two bioactive phytocompounds evaluated, 6-Hydroxyondansetron (A104) passes all the required drug discovery rules, and thus confirmed as viable drug candidate which can be used as potential inhibitor for the two SARS-CoV-2 replication proteins (Mpro and PLpro). However, Quercitrin (A166) has demonstrated relatively less druglike properties but can also be used as potential inhibitor for ADRP. 

World Health Organization (2020) Weekly operational Update on COVID-19

Al Zarzour R (2020) COVID-19 and therapy with essential oils having antiviral, antiinflammatory, and immunomodulatory properties

From SARS to MERS, thrusting coronaviruses intothe spotlight

Proteolytic processing, deubiquitinase and interferon antagonist activities of Middle East respiratory syndrome coronavirus papain-like protease

Structure-based screening of novel lichen compounds against SARS Coronavirus main protease (Mpro) as potentials inhibitors of COVID-19

Papain-like protease regulates SARS-CoV-2 viral spread and innate immunity

Essential oils as antiviral agents, potential of essential oils to treat SARS-CoV-2 infection: an in-silico investigation

Consensus virtual screening of dark chemical matter and food chemicals uncover potential inhibitors of SARS-CoV-2 main protease

Identification of potent COVID-19 main protease (Mpro) inhibitors from natural polyphenols: an in-silico strategy unveils a hope against CORONA

Recognition of natural products as potential inhibitors of COVID-19 main protease (Mpro): In-silico evidences

In silico studies reveal potential antiviral activity of phytochemicals from medicinal plants for the treatment of COVID-19 infection

Structure-based drug designing for potential antiviral activity of selected natural products from Ayurveda against SARS-CoV-2 spike glycoprotein and its cellular receptor

Virtual screening of natural metabolites and antiviral drugs with potential inhibitory activity against 3CL-PRO and PL-PRO

Molecular docking study of potential phytochemicals and their effects on the complex of SARS-CoV2 spike protein and human ACE2

Antiviral activity of plants and their isolated bioactive compounds: an update

Immunomodulatory and anti-SARS activities of Houttuynia cordata

Houttuynia cordata Thunb: a review of phytochemistry and pharmacology and quality control

Antiviral, antifungal and antibacterial activities of the Chinese medicinal plants, Houttuynia cordata, Lobelia chinensis and Selaginella uncinata

Botanical remedies of the former Dutch East Indies (Indonesia): part I: Eumycetes, pteridophyta, gymnospermae, angiospermae (monocotyledones only)

First report on the ethnopharmacological uses of medicinal plants by Monpa tribe from the Zemithang Region of Arunachal Pradesh, Eastern Himalayas

Ethnomedicinal plant knowledge among the Adi Tribe of Yingkiong and Mariyang Valley

Ethnomedicinal botany of the Apatani in the Eastern Himalayan region of India

Houttuynia cordata targets the beginning stage of herpes simplex virus infection

The inhibitory actions of Houttuynia cordata aqueous extract on dengue virus and dengueinfected cells

Flavonoids from Houttuynia cordata attenuate H1N1-induced acute lung injury in mice via inhibition of influenza virus and Toll-like receptor signalling

The effect of Houttuynia cordata injection on pseudorabies herpes virus (PrV) infection in vitro

Antiviral effects of Houttuyniacordata polysaccharide extract on Murine Norovirus-1 (MNV-1): a human norovirus surrogate

Evaluation of antiviral activities of Houttuynia cordata Thunb. extract, quercetin, quercetrin and cinanserin on murine coronavirus and dengue virus infection

Houttuynia cordata modulates oral innate immune mediators: potential role of herbal plant on oral health

Amphiphile behavior in mixed solvent media I: self-aggregation and ion association of sodium dodecylsulfate in 1, 4-dioxane-water and methanol-water media

Characterization of gingerol-related compounds in ginger rhizome (Zingiber officinale Rosc.) by high-performance liquid chromatography/electrospray ionization mass spectrometry

Evaluation of changes induced in rice metabolome by Cd and Cu exposure using LC-MS with XCMS and MCR-ALS data analysis strategies

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

Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments

Going further than Lipinski's rule in drug design

The natural way forward: Molecular dynamics simulation analysis of phytochemicals from Indian medicinal plants as potential inhibitors of SARS-CoV-2 targets

Reversible multiple time scale molecular dynamics

How computational chemistry and drug delivery techniques can support the development of new anticancer drugs

In silico tools to study molecular targets of neglected diseases: inhibition of TcSir 2 rp, an epigenetic enzyme of Trypanosomacruzi

Informatics for chemistry, biology and biomedical sciences

Modelling and forecasting of COVID-19 spread using wavelet-coupled random vector functional link networks

Conflict of interest Authors declare no conflict of interest.