key: cord-0792306-4iyk76sb
authors: Alagu Lakshmi, Selvaraj; Shafreen, Raja Mohamed Beema; Priya, Arumugam; Shunmugiah, Karutha Pandian
title: Ethnomedicines of Indian origin for combating COVID-19 infection by hampering the viral replication: using structure-based drug discovery approach
date: 2020-06-23
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1778537
sha: 9e239a682bd72830c7c0b7cd9b0c4ac5d494b7a1
doc_id: 792306
cord_uid: 4iyk76sb

In the present study, we have explored the interaction of the active components from 10 different medicinal plants of Indian origin that are commonly used for treating cold and respiratory-related disorders, through molecular docking analysis. In the current scenario, COVID-19 patients experience severe respiratory syndromes, hence it is envisaged from our study that these traditional medicines are very likely to provide a favourable effect on COVID-19 infections. The active ingredients identified from these natural products are previously reported for antiviral activities against large group of viruses. Totally 47 bioactives identified from the medicinal plants were investigated against the structural targets of SARS-CoV-2 (Mpro and spike protein) and human ACE2 receptor. The top leads were identified based on interaction energies, number of hydrogen bond and other parameters that explain their potency to inhibit SARS-CoV-2. The bioactive ligands such as Cucurbitacin E, Orientin, Bis-andrographolide, Cucurbitacin B, Isocucurbitacin B, Vitexin, Berberine, Bryonolic acid, Piperine and Magnoflorine targeted the hotspot residues of SARS-CoV-2 main protease. In fact, this protease enzyme has an essential role in mediating the viral replication and therefore compounds targeting this key enzyme are expected to block the viral replication and transcription. The top scoring conformations identified through docking analysis were further demonstrated with molecular dynamics simulation. Besides, the stability of the conformation was studied in detail by investigating the binding free energy using MM-PBSA method. Overall, the study emphasized that the proposed hit Cucurbitacin E and orientin could serve as a promising scaffold for developing anti-COVID-19 drug. Communicated by Ramaswamy H. Sarma

At the peak of COVID19 pandemic which shudder the whole world through human to human transmission of SARS-CoV-2, devouring numerous lives, medical professionals and researchers all over the globe are in pursuit for the discovery of prophylactic or anaphylactic remedy to surmount the condition (Al-Khafaji et al., 2020; Boopathi et al., 2020; Pant et al., 2020) . Thus far, no specific vaccine or antiviral therapeutic has been either discovered or anticipated as a conventional medication for this viral infection (Zumla et al., 2016) . Current treatment strategy includes supportive care which is extended with supplementation of combination of broad-spectrum antibiotics, antivirals, cortecosteroids and convalescent plasma (Chen et al., 2020; Habibzadeh & Stoneman, 2020) . SARS-CoV-2 was found to be an enveloped single-stranded RNA-type beta coronavirus which shares 79.5% sequence similarity with its antecedent class SARS-CoV (Chen et al., 2020; Wu et al., 2020) . During the previous outbreak of SARS crisis due to SARS-CoV, several research groups have conducted in vitro analysis and randomized trials on use of natural herbal medicines to control the adverse effects of infection and evidenced their use as effective (Hsu et al., 2008; Jang et al., 2009; Kim et al., 2008; Lau et al., 2005a Lau et al., , 2005b Li et al., 2005; Wen et al., 2011) . As novel coronavirus is closely related to SARS coronavirus, it can be suggested that natural herbal medicine may have potential use in existing health crisis (Aanouz et al., 2020; Elfiky, 2020b; Islam et al., 2020; Umesh et al., 2020; Wahedi et al., 2020) . In addition to this, for decades, natural herbal medicine has been expended for prevention and cure of numerous ailments concerning human wellbeing.

The WHO estimates that, 80% of population from underdeveloped realm rely mostly on the traditional medicines. In addition to this, WHO listed around 21,000 plants which are known for their therapeutic potential around the world among which 2500 varieties are found around India including 150 species with large scale commercial usage (Shukla et al., 2019) . In India, Siddha is one among the three traditional medical systems practiced and the other two being Ayurveda and Unani. In Siddha medicine Nilavembu Kudineer Chooranam (NKC) has been recommended for the preclusion and control of all the types of viral infectious diseases. NKC is an herbal formulation comprising nine varieties of plant materials viz. Nilavembu (Andrographis paniculata), Vetiver (Vetiveria zizanioides), Vilamiccam ver (Plectranthus vettiveroides), Santanam (Santalum album), Pei pudal (Trichosanthes dioica), Korai kizhangu (Cyperus rotandus), Chukku (Zingeber officinale), Milagu (Piper nigrum) and Parpatakam (Mollugo cerviana) (Kavinilavan et al., 2017) . NKC has been reported to have therapeutic applications such as antiviral, antipyretic, antibacterial, antiulcer, antioxidant, analgesic etc . On the other hand, it has also been stated to act as an immunostimulant and immunomodulant which enhances the host immune system and modulates the defence response that helps to protect the complications of infections (Kamalarajan et al., 2019; Nakkeeran et al., 2016; Ramanathan et al., 2019) . Previously, NKC has been shown to alleviate the symptoms of dengue and chikungunya (Anbarasu et al., 2011; Mattummal et al., 2018) . In fact, Government of Tamil Nadu have taken initiatives in supplying decoction of NKC to patients visiting the primary health care centres and also promoted the drink for prevention and control of morbidity level during outbreak of viral fever (Ramanathan et al., 2019) . Owing to their medicinal properties and specifically antiviral potential, the present study aims at investigating the therapeutic efficacy of active components of NKC and natural constituents from other known medicinal plants (Tinospora cordifolia and Anethumsowa) which are also used as combinations to treat respiratory diseases. Phytomolecules from the components of NKC and the aforementioned medicinal plants against the two major viral targets spike glycoprotein and protease and a host target human angiotensin-converting enzyme 2 (ACE2) receptor, which is the entry point for viral encounter were examined with the prospects of identifying potential drug candidate(s) against COVID19 infection. Thus, ethnomedicines identified from the study through in silico approaches can combat the global health crisis instigated by SARS-CoV-2 (Abdelli et al., 2020; Elfiky & Azzam, 2020; Elmezayen et al., 2020; Enayatkhani et al., 2020; Enmozhi et al., 2020 , Elfiky, 2020a Gupta et al., 2020; Sarma et al., 2020; Sinha et al., 2020) . Unlike modern medicines, the antique herbal remedies do not possess side effects, indeed they establish natural resistance to alleviate the viral infection by boosting the immune system. Thus, the main objective of the study is to identify a potent phytocompound as inhibitor against SARS-CoV-2 drug targets from vintage medicinal plants through computational investigations.

The study focuses predominantly on bioactive ligands from ten traditional medicinal plants that are primarily used in the treatment of respiratory illness. The active components were selected and screened against SARS-CoV-2 main protease, spike (S) protein and host cell ACE2 receptor. The list of active leads used in the present study is shown in Table 1 . Ligand 3D structure was obtained in .sdf format from PubChem database and converted into .pdb format through online program (https://cactus.nci.nih.gov/translate/). PubChem ID of bioactives are provided in Supplementary Table 1 .

Crystal Structures of SARS-CoV-2 main protease (PDB ID: 5R82.pdb), spike protein (PDB ID: 6VYB.pdb) and human entry receptor ACE2 (PDB ID:1R42.pdb) were retrieved from RCSB PDB database (http://www.rcsb.org/pdb). In order to minimize the energy, PDB structures of the target proteins were refined through ModRefiner algorithm (Xu & Zhang, 2011) . Ligands were optimized by adding hydrogen atoms and energy minimization was performed with MMFF94 Force Field in Autodock v4.2 (Norgan et al., 2011) . Grid box was defined on binding pocket of target protein and the grid points were extended in all directions to encompass the binding region. Molecular docking was performed through Autodock v4.2 by tethering ligand to target protein and the binding affinity was determined through kcal mol À1 . The pose with highest binding energy was extracted and analysed in Biovia Discovery Studio 4.5 software (Systemes, 2015) .

The top two ligands showing interaction with Mpro (revealed through docking analysis) were further subjected to 20 ns MD simulation studies. The energy and structural properties of the ligand-receptor complexes were determined through Gromacs 5.1.4 simulation package using CHARMM force field. The protein-ligand complex was inserted into a cubic box of 10 Å dimension and solvated using SPCE water model and the system was neutralized by adding counter ions (0.15 M NaCl). An initial minimization of 1000-step was carried out by applying steepest descent algorithm followed by 2000-step using conjugate gradient minimization. The system was gradually heated to 300 K for 100 ps and then equilibrated using NVT ensemble for 100 ps and retained at constant temperature. Sequentially, the system was equilibrated using NPT ensemble for 1000 ps. MD production run for 20 ns was performed with a time step of 2 ps and the trajectories were recorded at every step (Beema Shafreen et al., 2014) . The structural behaviour of Mpro protein in the presence and absence of the drug molecules were calculated through Root mean square deviation (RMSD) and Root mean square fluctuation (RMSF). Overall methodology is represented as a flowchart in Figure 1 (Borgio et al., 2020) .

The binding free energy between the protein and drug molecules were calculated using molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) method in Gromacs (g_mmpbsa) software. The following equation was used to calculate the binding free energy.

Where DG Bind is the binding free energy; DE MM is the interaction energies such as DE vdw (van der Waals interaction energy) and DE elec (electrostatic interaction energy) between the protein and drug complex; DG sol is the solvation energy comprising DG GB (polar components) and DG SA (non-polar components); TDS is the entropy change at temperature T.

Crystal structures of main protease, spike protein and human entry receptor (ACE2) were used to define the molecular interactions (Meng et al., 2011) with the scrutinized bioactives. The top ten scoring bioactive ligands are represented in Table 2 with binding energy, dock score, van der Waals and H-bond energy. Interestingly, among the bioactives screened, the compounds Orientin, Vitexin, Berberine, Bryonolic acid and Magnoflorine were found to interact with all the drug targets of SARS-CoV-2 including the human entry receptor ACE2 (Figure 2 ). Although the aforesaid compounds had interaction with all the drug targets of SARS-CoV-2, interaction with active site residues of main protease Table 1 . List of natural products that are commonly used to treat fever, cold and respiratory related ailments in traditional South Indian medicine system. The active ingredients extracted from these natural products are used for docking with SARS-CoV-2 drug targets.

Tamil name Scientific name Accepted Scientific name Active ingredients was found to be much stronger than the other two drug targets (spike protein and ACE2) used in the study (Hasan et al., 2020) . The highest binding affinity between the ligandsreceptor complex was confirmed through binding energy and other possible non-covalent interactions. Based on the scoring values it is confirmed that these compounds can significantly inhibit SARS-CoV-2 viral replication by binding strongly with the main protease.

Among 47 bioactives used in the study, Orientin, a C-glycosyl compound, showed the highest binding affinity with main protease with a dock score of À134.14 and binding energy of À90.2 kcal mol À1 (Figure 2-IB) . It formed seven conventional hydrogen bonds with five residues (Thr26, His41, Gly143, Cys145 and Glu 166). These hydrogen bond interactions indeed enhance the binding affinity of the complex by stabilizing the ligand at the active site of the target and thereby help in increasing the biological activity of the complex. Notably, among the active site residues the catalytic dyad of main protease (His41 and Cys145) was mainly involved in hydrogen bond interactions. These interactions are expected to prevent the virus from processing polyprotein which consequently hinder the viral replication. The other interactions that stabilizes the Orientin and main protease complex includes pi-pi stacked interaction with His163, pi-pi T-shaped interaction with His41, pi-alkyl interaction with Cys145, pi-sulfur interaction with Met49. In addition, the residue Asn142 formed a pi-donor hydrogen bond with the ligand molecule. The conformation had the van der Waals and H-bond energy of À63.65 and À26.58 kcal mol À1 , respectively. All these crucial interactions facilitate in augmenting the efficacy of the drug molecules by stabilizing the binding structures. Orientin docked with spike protein (6VYB) revealed the docking and binding energy score of À110.05 and À72.3 kcal mol À1 , respectively. Besides, it formed six hydrogen bond interactions with four amino acids (Tyr28, Asn30, Phe59 and Asn61) and one carbon-hydrogen bond with Thr29. H-bond and van der Waals interaction energies were found to be À13.23 and À59.02 kcal mol À1 , respectively. Similarly, Orientin docked with ACE2 receptor revealed the docking and binding energy score of À101.17 and À70.6 kcal mol À1 , respectively. The binding affinity was associated with five hydrogen bond interactions formed by three residues (Glu22, Lys26 and Asn90), one carbon hydrogen bond interaction with Glu22 and pi-alkyl interaction with Lys94. The structure possesses van der Waals and H-bond interactions energy of À58.84 and À11.76 kcal mol À1 , respectively.

Vitexin, an apigenin flavone glycoside, was found to have the binding energy of À85.7 kcal mol À1 and the dock score of À98.03 during the interaction with main protease. The strong affinity was explicable from six hydrogen bonds formed with six residues of main protease (Phe140, Ser144, Cys145, Asn142, Leu167 and Gln189). Vitexin-Mpro complex was stabilized by pi-anion (Cys145), pi-sulfur (Glu166), pi-lone pair (Asn142), pipi T-shaped (His41) and pi-alkyl (Met165) interactions. The van der Waals and H-bond energy was found to be À72.97 and À12.71 kcal mol À1 , respectively. Vitexin docked with spike protein reported À98.34 and À67.7 kcal mol À1 as docking and binding energy score, respectively. It formed four hydrogen bonds with four residues (Tyr28, Asn30, Phe59 and Asn61) and one carbon hydrogen bond with Thr29. The conformation had the van der Waals and H-bond energy of À57.78 and À9.93, respectively. ACE2 receptor with Vitexin had a docking and binding energy score of À85.24 and À68.1 kcal mol À1 , respectively. Conformational stability of the complex (Vitexin-ACE2) is maintained by hydrogen bond (two hydrogen bonds from Glu22 and Asn90), van der Waals (Lys26, Leu29, Gln89 and Gln96) and pi-alkyl interactions (Val93 and Lys 94). H-bond and van der Waals energy was found to be À57.61 and À10.5 kcal mol À1 , respectively.

Main protease and Berberine (an organic heteropentacyclic compound) interaction showed a binding energy of À83.2 kcal mol À1 and a dock score of À103.89. The interaction between ligand with the receptor binding site of Mpro was stabilized by pi-pi T-shaped interaction with His163, Pi-alkyl interaction with Cys145 & His41 and alkyl interaction with His163. The conformation showed van der Waals and H-bond energy of À74.84 and À8.34 kcal mol À1 , respectively. Berberine with spike protein revealed the docking score, binding and interaction energies such as H-bond and van der Waals as À82.48, À69.7, À4.89 and À64.84 kcal mol À1 , respectively. The other stabilizing interactions includes van der Waals (Asn30), amide pi-stacked (Phe59), pi-pi stacked (Thr29) and pi-alkyl (Phe59) interactions. Berberine with ACE2 receptor had a docking score of À95.33, binding energy of À71.5 kcal mol À1 and H-bond and van der Waals interaction energies of À58.759 and À12.728, respectively. It formed two conventional hydrogen bond interactions with Lys26 and Thr27. Further stabilization occurred with pication (Glu22), pi-anion (Glu22), pi-sulfur (Val93) and alkyl (Val93) interactions. Bryonolic acid formed a strong interaction with the receptor binding sites of main protease with a binding energy score of À81.0 kcal mol À1 and dock score of À91.63. Bryonolic acid formed two hydrogen bond interactions with Arg188 and Gln189. Further conformational interactions such as carbon hydrogen bond interaction was observed with Ser139 and His163, alkyl interaction with Tyr118, Cys145, His163 and Met49 and pi-alkyl interaction with His41, His163, His172 and Met165 which stabilized the binding of the bryonolic with Mpro. The van der Waals and H-bond interaction energy was À74.25 and À6.71 kcal mol À1 , respectively. Bryonolic acid docked with spike protein reported À80.57, À63.9, À60.37 and À3.5 kcal mol À1 as docking, binding, Hbond and van der Waals energies, respectively. Key interactions between Bryonolic acid and spike protein includes one conventional hydrogen bond with Tyr28 and pi-alkyl interaction with Phe58. Interaction with ACE2 receptor was found to have a docking score of À80.99, binding energy score of À65.6 kcal mol À1 and van der Waals and H-bond interaction score of À54.032 and À11.596 kcal mol À1 , respectively. Energetic contributions to the ligand include conventional hydrogen bonds (four bonds with Glu22, Thr27, Gln89 and Asn 90) and pi-alkyl (Leu29 and Lys94) interactions.

Magnoflorine docked with main protease showed a binding energy score of À77 kcal mol À1 and a dock score of À92.41. It formed six hydrogen bonds with four amino acids (Leu141, Asn142, Ser144 and Cys145). Besides, it also forms carbon hydrogen bond with Phe140, His163 and Glu166 and pi-donor hydrogen bond with Cys145. The extra stabilizing interaction associated with ligand and main protease complex includes pi-cation interaction with His41, pi-alkyl interaction with Leu27 and Cys145 and pi-sigma interaction with His41. The van der Waals and H-bond interaction energy was À64.02 and À12.98 kcal mol À1 , respectively. Magnoflorine with spike protein has a dock score of À90.88, binding energy score of À65.5 kcal mol À1 and van der Waals and Hbond interactions score of À57.44 and À8.02 kcal mol À1 , respectively. It formed conventional hydrogen bonds with Asn61 and carbon hydrogen and pi-donor hydrogen bond with Tyr28. Other significant interactions comprise pi-sigma (Tyr28) and pi-pi stacked (Tyr28) interactions. Magnoflorine with ACE2 receptors exhibited a docking score of À95.46, binding energy score of À64.4 kcal mol À1 , H-bond and van der Waals interaction scores of À15.52 and À48.90 kcal mol À1 , respectively. It formed two hydrogen bond interactions with Thr27 and Leu29 and carbon hydrogen bond with Val93 and Glu22. Other energetic contribution includes pication (Lys26 and Glu22) and pi-alkyl (Leu29 and Val23) interactions. In addition, it formed electrostatic interaction with Lys26 and Glu22, thereby stabilizing the complex formed.

Besides, the compound Cucurbitacin E, Bis-andrographolide, Cucurbitacin B, Isocucurbitacin B and Piperine had strong binding affinity to main protease of SARS-CoV-2. Cucurbitacin E was the top scorer among the main protease inhibitors. It had the highest binding energy score of À91.0 kcal mol À1 with a dock score of À113.4. The H-bond and van der Waals interactions energy was found to be À19.48 and À71.53 kcal mol À1 , respectively. It formed two conventional hydrogen bond interactions with Asn142 and Cys145 and one carbon hydrogen bond with Ser46. Besides, it had alkyl and pi-alkyl interaction with Met49, His41, Cys44 and His163. The compounds Cucurbitacin B and Isocucurbitacin B had similar kind of pi-alkyl and alkyl interactions as Cucurbitacin E. But in case of conventional hydrogen bond interactions both lacked Cys145 interactions. In addition, Isocucurbitacin B lacked carbon hydrogen bond interaction with Ser46. Notably, Cucurbitacin B formed a specific conventional hydrogen bond with Thr25.

Bis-andrographolide had the docking and binding energy score of À100.3 and À90 kcal mol À1 , respectively. It formed conventional hydrogen bond with Glu166 and alkyl interaction with Met49. The conformation had the van der Waals and H-bond energy score of À88.60 and À1.38 kcal mol À1 , respectively. Piperine interacted with main protease with a dock score of À90.95 and binding energy score of À78.10 kcal mol À1 . It formed one conventional hydrogen bond with hotspot residue His41 and other stabilizing interactions include pi-sulfur, pi-sigma, pi-pi T-shaped and alkyl interactions.

The compounds Palmatine, Piperine, Stigmasterol, Andrographolide and b-sitosterol and the components 2-monolinolenin, Andrographolide, Isocucurbitacin B, Cucurbitacin E and B had strong binding affinity to spike and ACE2 receptor, respectively. The schematic representation summarizing the overall scoring functions of top ten bioactive molecules is shown in Figure 3 .

Among the ligands screened through docking studies, Cucurbitacin E and Orientin revealed significant interaction with most critical residues of Mpro with high binding affinities. Therefore, these two compounds were selected as the best inhibitors and used for further analysis. MD simulation of 20 ns was performed for the two selected drug molecules to get deep insights into the stability of the Mpro-drug complex. The overall stability was determined through RMSD and RMSF analysis. The energy profile in the presence and absence of ligands were computed during 20 ns MD simulation. Wherein, the potential energy profile of the system was found to be stable throughout the simulation and the lowest energy profile in the presence of ligands indicates that the system is considerably stable than in the absence of ligand (Figure 4(A) ). Similarly, RMSD analysis during each time step unveiled the stability of the system. Mpro in the absence of drug molecule maintained the steady state throughout the simulation with an average RMSD value of 3.71 Å, indicating the convergence of the system during 20 ns simulation. However, Mpro in complex with drug molecules exhibited little variations during the simulation run (Figure 4(A) ). Orientin presented an average RMSD of 4.6 Å whereas Cucurbitacin E presented an average RMSD of 4.9 Å. Although, the system started with increasing RMSD profile, it favored to equilibrium state after 2 ns and successfully maintained the stable state throughout the remaining simulation run (Figure 4(B) ). RMSF provides the visibility for most flexible and mobile regions in the protein during the simulation. In fact, these structural fluctuations reveal the binding and interaction of the drug molecules with the protein motif. RMSF analysis of Mpro in the presence of ligand unveiled the fluctuating residues at 1-12, 128-146 and 185-192 amino acid regions corresponding to the loop region (Figure 4(C) ) of Mpro. The average RMSF value for Mpro was found to be 1.53 Å revealing that there was little fluctuation in the loop region during MD run. However, Mpro in complex with Cucurbitacin E and Orientin was found to have high fluctuation in the druggable pocket with an average RMSF value of 2.10 Å and 1.79 Å, respectively. The Mpro-Cucurbitacin E complex was found to have significant fluctuation with the catalytic residues (Cys145) of Mpro in comparison to Orientin. Since, Cys145 resides in the loop region, the fluctuations are favorable and advantageous in accommodating the drug molecule at the receptor site.

MM-PBSA is the most reliable method for predicting the binding interaction between protein-drug complexes. The binding free energy for Cucurbitacin E and Orientin was À152.96 and À123.62 kcal/mol respectively. The binding energies obtained from MD trajectories further confirms the accommodation of inhibitors into the druggable pocket of Mpro. Moreover, the binding pocket with polar microenvironment reinforces the strong hydrogen and hydrophobic interactions with the drug molecule and seizes them in the active site of Mpro. Overall, the results were suggestive that Mpro maintains a stable binding state with the inhibitors (Cucurbitacin E and Orientin) for plausible interactions.

Siddha, one of the oldest traditional systems of medicine in India, has been shown to be effective in treating respiratoryrelated infections (Ram et al., 2009) . Several natural bioactives isolated from medicinal plants of Siddha system are reported for antiviral property against respiratory viruses. It is thus likely that these compounds may exert an antiviral activity against SARS-CoV-2 (Adeoye et al., 2020; Oliveira et al., 2020) . Hence, in the present study anti-COVID-19 drug candidates were identified from Siddha medicinal plant, with special attention to NKC combination plants, due to its antiviral, immunostimulant and immunomodulatory activities (Kavinilavan et al., 2017) . The potent anti-COVID-19 drug candidates were screened from 47 bioactives, among them 5 (Orientin, Vitexin, Berberine, Bryonolic acid and Magnoflorine) bioactives had potential interactions with all drug targets (mentioned in the present study). All these drug candidates exhibited strong affinity towards the target site of the receptor molecule with high binding energy. The top scoring receptor-ligand complex was stabilized by non-covalent interactions such as hydrogen-bonding, van der Waals and electrostatics interactions (Chen & Kurgan, 2009 ). These interactions are indeed prerequisites for biological functions and successful drug developments. Moreover, the top two scoring ligands Cucurbitacin E and Orientin were further substantiated through 20 ns MD simulation studies.

Mpro, one of the best characterized drug targets of SARS-CoV-2, is essential for polyprotein processing (translated from viral RNA) that enables viral replication and gene expressions. The catalytic dyad His41 and Cys145 is essential for substrate hydrolysis. Thus, compounds inhibiting the active site of this protease are expected to prevent replication and transcription of viral components (Das et al., 2020; Gyebi et al., 2020; Joshi et al., 2020; Khan et al., 2020a; 2020b; Kumar et al., 2020; Muralidharan et al., 2020; Zhang et al., 2020a; Zhang et al., 2020b) . Other drug targets of SARS-CoV-2 include spike glycoprotein and host cell receptor ACE2. In general, SARS-CoV-2 makes an entry into host cell through ACE2 receptors and it is mediated by receptor binding domain of homo trimeric spike glycoprotein. The crystal structure of spike protein in complex with human ACE2 receptor revealed that residue Ser19, Gln24,Thr27, Phe28, Asp30, Lys31, His34, Glu35,Glu37, Asp38, Tyr41, Gln42, Leu45, Leu79, Met82, Tyr83,Asn330,Lys353, Gly354, Asp355, Arg357 and Arg393 of human ACE2 interacted with residues Ala475/Gly476, Ala474/G476/Asn487, Phe56/ Tyr473/Ala475/Tyr489, Tyr489, Lys417/Leu455/Phe456, Leu455/Phe456/Glu484/Tyr489/Phe490/Gln493, Tyr453/ Leu455/Gln493, Gln493, Tyr505, Tyr449/Gly496/Gln498, Gln498/Thr500/Asn501, Gly446/Tyr449/Gln498, Gln498/ Thr500, Phe486, Phe486, Phe486/Asn487/Tyr489, Thr500, Gly496/Asn501/Gly502/Tyr505, Tyr502/Tyr505, Thr500,Gly502, Thr500, Tyr505, respectively, of spike glycoprotein . Therefore, bioactives targeting these hotspot residues are likely to become a potent inhibitor of SARS-CoV-2. However, spike protein and ACE2 receptor showed no significant inhibition. In case of ACE2 drug target, only Thr27 (hotspot residue) of ACE2 had interactions with bioactives (Veeramachaneni et al., 2020) .

Orientin possesses wide range of pharmacological effects such as antioxidant, antiaging, antibacterial, anti-inflammation, free radical scavenging, cardioprotective, radioprotective, neuroprotective, antidepressant, anticancer and thrombocytopenia activities (Lam et al., 2016; Sharma et al., 2016; Yadav et al., 2018) . It has also been documented for antiviral activity against parainfluenza type 3 virus (PIV-3) (Li et al., 2002) . Like Orientin, the compound Vitexin was also reported for multiple drug targets including antioxidant, anticancer, anti-inflammatory, antihyperalgesic and neuroprotective effects (He et al., 2016) . It has also been reported for antiviral activity against several viruses such as herpes simplex virus type 1 (HSV-1), Hepatitis A Virus (HAV-H10), PIV-3 and rotavirus (Fahmy et al., 2020; Knipping et al., 2012; Li et al., 2002) . Bryonolic acid, a pentacyclic triterpenoids, exhibits biological activities such as antiallergic, antitumor, antioxidant and anti-inflammatory (Gatbonton-Schwager et al., 2012; Khallouki et al., 2018; Tanaka et al., 1991) . An aporphine alkaloid, Magnoflorine, has been reported for antidiabetic and anticancer property (Patel & Mishra, 2012) . In addition, Magnoflorine has been previously reported for their antiviral effects against HSV-1 and poliovirus type-1 (Mohamed et al., 2010) . Berberine, an organic hetero pentacyclic compound has been reported for numerous therapeutic values such as antioxidant, anti-inflammatory, antidiarrhea, antilipemic, antibiofilm, immunomodulatory, cardioprotective, nephroprotective, hepatoprotective and glucose metabolism (Neag et al., 2018; Zhang & Shen, 1989) . Most importantly, it has been reported for antiviral activity against numerous viruses such as influenza virus , enterovirus 71 (Wang et al., 2017) , human cytomegalovirus (Hayashi et al., 2007) , hepatitis C virus (Hung et al., 2019) , HSV-1 (Dkhil & Al-Quraishy, 2014) and chikungunya virus (Varghese et al., 2016) . Besides, it has been proven to be an antibacterial agent against multiple drug resistance bacteria such as methicillin-resistant Staphylococcus aureus and Streptococcus agalactiae (Chu et al., 2016; Peng et al., 2015; Tan et al., 2019; Yu et al., 2005) and antifungal agent against Candida spp. (Xie et al., 2020) .

The compounds Piperine, Cucurbitacin B, Cucurbitacin E, Isocucurbitacin B and Bis-andrographolide have also been reported for numerous pharmacological effects (Chao & Lin, 2010; Derosa et al., 2016; Kaushik et al., 2015) . In addition, Piperine has been reported for antiviral activity against Hepatitis B (Jiang et al., 2013) . Cucurbitacin B showed antiviral activity against HSV-1 (Hassan et al., 2017) and Bovine viral diarrhea virus (BVDV) (Alsayari et al., 2012) . Cucurbitacin E was reported against BVDV (Alsayari et al., 2012) and Bisandrographolide against HIV virus (Reddy et al., 2005) .

The multiple pharmacological effects exhibited by these components strongly suggest that these phytochemicals would certainly find their way into the arsenal of antiviral drugs. Moreover, the active components Vitexin, Berberine, Bryonolic acid, Magnoflorine, Cucurbitacin, Piperine, Palmatine, Stigmasterol, Andrographolide, Sitosterol and 2monolinolenin were FDA approved drugs. Hence, these phytochemicals can be further taken for preclinical studies to treat COVID-19 infection (Lobo-Galo et al., 2020) . The docking studies expounded with MD simulation and binding affinity analyses strongly suggest that the compound Cucurbitacin E and Orientin are promising drug candidate to combat COVID-19 infections. However, in vitro and in vivo experimental evidences are further required to substantiate the postulations observed from the present study. In conclusion, consumption of NKC will helps in developing an immune response, elimination of toxins from the system, soothes the body and improves the overall health condition.

The present study explored the important bioactive constituents of Siddha medicine, especially NKC through molecular docking and simulation analysis for the prevention and cure of COVID19 infection. The compounds Cucurbitacin E, Orientin, Bis-andrographolide, Cucurbitacin B, Isocucurbitacin B, Vitexin, Berberine, Bryonolic acid, Piperine and Magnoflorine were identified as potential lead molecules that have been shown to possess the ability to interact with the components that block the viral replication in SARS-CoV-2. Moreover, the immune-enhancing properties of these compounds without any adverse side effects could provide natural immune power to resist COVID-19 infections. However, preclinical studies would provide deep insight into mechanism of actions and also drugability of these bioactives for the treatment of COVID19 infection.

The authors declare that they have no conflicts of interest.

The authors sincerely acknowledge the computational and bioinformatics analysis provided by the Bioinformatics Infrastructure Facility 

The data that support the findings of this study are openly available in GenBank database. 

Moroccan Medicinal plants as inhibitors of COVID-19: Computational investigations

In silico study the inhibition of Angiotensin converting enzyme 2 receptor of COVID-19 by Ammoides verticillata components harvested from western Algeria

Repurposing of chloroquine and some clinically approved antiviral drugs as effective therapeutics to prevent cellular entry and replication of coronavirus

Using Integrated Computational Approaches to Identify Safe and Rapid Treatment for SARS-CoV-2

Antibovine viral diarrhea virus activity of cucurbitacins as new potential antiviral agents

Antipyretic, anti-inflammatory and analgesic properties of nilavembu kudineer choornam: A classical preparation used in the treatment of chikungunya fever

In silico and in vitro studies of cinnamaldehyde and their derivatives against LuxS in Streptococcus pyogenes: Effects on biofilm and virulence genes

Novel 2019 coronavirus structure, mechanism of action, antiviral drug promises and rule out against its treatment

State-of-the-art tools unveil potent drug targets amongst clinically approved drugs to inhibit helicase in SARS-CoV-2

Isolation and identification of bioactive compounds in Andrographis paniculata (Chuanxinlian)

Investigation of atomic level patterns in protein-small ligand interactions

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Role of berberine in the treatment of methicillin-resistant Staphylococcus aureus infections

An investigation into the identification of potential inhibitors of SARS-CoV-2 main protease using molecular docking study

Piperine and its role in chronic diseases

Evaluation of antiviral activity of berberine against herpes simplex viruses

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: An in silico perspective

Natural products may interfere with SARS-CoV-2 attachment to the host cell

Novel guanosine derivatives against MERS CoV polymerase: An in silico perspective

Drug repurposing for coronavirus (COVID-19): in silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes

Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: An in silico study

Andrographolide as a potential inhibitor of SARS-CoV-2 main protease: An in silico approach

Breaking down the barriers to a natural antiviral agent: Antiviral activity and molecular docking of erythrina speciosa extract, fractions, and the major compound

Bryonolic acid transcriptional control of anti-inflammatory and antioxidant genes in macrophages in vitro and in vivo

In-silico approaches to detect inhibitors of the human severe acute respiratory syndrome coronavirus envelope protein ion channel

Potential Inhibitors of Coronavirus 3-Chymotrypsin-Like Protease (3CLpro): An in silico screening of Alkaloids and Terpenoids from African medicinal plants

The novel coronavirus: A Bird's eye view

A review on the cleavage priming of the spike protein on coronavirus by angiotensin-converting enzyme-2 and furin

Antiviral activity of berberine and related compounds against human cytomegalovirus

A review on the pharmacological effects of vitexin and isovitexin

An evaluation of the additive effect of natural herbal medicine on SARS or SARS-like infectious diseases in 2003: A randomized, double-blind, and controlled pilot study

Berberine inhibits hepatitis C virus entry by targeting the viral E2 glycoprotein

A molecular modeling approach to identify effective antiviral phytochemicals against the main protease of SARS-CoV-2

Antiviral activity of ethanolic extract of Nilavembu Kudineer against dengue and chikungunya virus through in vitro evaluation

An overview of the herbal remedies for Severe Acute Respiratory Syndrome (SARS) in WHO official report

Anti-HBV active constituents from Piper longum

Discovery of potential multi-target-directed ligands by targeting host-specific SARS-CoV-2 structurally conserved main protease

Phytochemical investigation of nilavembu kudineer chooranam ethyl acetate extract and its ability to reduce intracellular antioxidant levels in THP-I cells

Cucurbitacins-An insight into medicinal leads from nature

Exploration of immunomodulatory effect of nilavembu kudineer chooranam against newcastle disease virus in backyard chicken

Bryonolic acid blocks cancer cell clonogenicity and invasiveness through the inhibition of fatty acid: Cholesteryl ester formation

Targeting SARS-Cov-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like Proteinase and 2'-O-RiboseMethyltransferase

Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach

In vitro inhibition of coronavirus replications by the traditionally used medicinal herbal extracts, Cimicifuga rhizoma, Meliae cortex, Coptidis rhizoma, and Phellodendron cortex

An evaluation of the inhibitory effects against rotavirus infection of edible plant extracts

Understanding the binding affinity of noscapines with protease of SARS-CoV-2 for COVID-19 using MD simulations at different temperatures

A review on medicinal properties of orientin

The use of an herbal formula by hospital care workers during the severe acute respiratory syndrome epidemic in Hong Kong to prevent severe acute respiratory syndrome transmission, relieve influenza-related symptoms, and improve quality of life: A prospective cohort study

Using herbal medicine as a means of prevention experience during the SARS crisis

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

FDA-approved thiol-reacting drugs that potentially bind into the SARS-CoV-2 main protease, essential for viral replication

Bioactive molecules in Siddha Polyherbal Nilavembu Kudineer alleviating symptoms of Dengue/Chikugunya

Molecular docking: A powerful approach for structure-based drug discovery

Cytotoxic and antiviral activities of aporphine alkaloids of Magnolia grandiflora L

Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 Protease against COVID-19

FTIR analysis on Nilavembu Kudineer Churanam and acetominaphen

Berberine: Botanical occurrence, traditional uses, extraction methods, and relevance in cardiovascular, metabolic, hepatic, and renal disorders

Multilevel parallelization of AutoDock 4.2

Repurposing approved drugs as inhibitors of SARS-CoV-2 S-protein from molecular modeling and virtual screening

Peptide-like and small-molecule inhibitors against Covid-19

Magnoflorine from Tinospora cordifolia stem inhibits a-glucosidase and is antiglycemic in rats

Antibacterial activity and mechanism of berberine against Streptococcus agalactiae

Medicinal plants from Siddha system of medicine useful for treating respiratory diseases

Formulation and evaluation of Nilavembu kudineer capsules

A new bis-andrographolide ether from Andrographis paniculata nees and evaluation of anti-HIV activity

In-silico homology assisted identification of inhibitor of RNA binding against 2019-nCoV N-protein (N terminal domain)

Structure-activity relationship studies on holy basil (Ocimum sanctum L.) based flavonoid orientin and its analogue for cytotoxic activity in liver cancer cell line HepG2

Indian herbal plants used as antipyretic: A review

An in-silico evaluation of different Saikosaponins for their potency against SARS-CoV-2 using NSP15 and fusion spike glycoprotein as targets

BIOVIA, discovery studio modeling environment

Antimicrobial characteristics of Berberine against prosthetic joint infection-related Staphylococcus aureus of different multi-locus sequence types

Anti-allergic effect of bryonolic acid from Luffa cylindrica cell suspension cultures

Identification of new anti-nCoV drug chemical compounds from Indian spices exploiting SARS-CoV-2 main protease as target

The antiviral alkaloid berberine reduces chikungunya virus-induced mitogenactivated protein kinase signaling

Structural and Simulation analysis of hot spot residues interactions of SARS-CoV 2 with human ACE2 receptor

Stilbene-based natural compounds as promising drug candidates against COVID-19

Berberine inhibits enterovirus 71 replication by downregulating the MEK/ERK signaling pathway and autophagy

Structural and functional basis of SARS-CoV-2 entry by using human ACE2

Traditional Chinese medicine herbal extracts of Cibotium barometz, Gentiana scabra, Dioscorea batatas, Cassia tora, and Taxillus chinensis inhibit SARS-CoV replication

Genome composition and divergence of the novel coronavirus (2019-nCoV) originating in China

In vivo and in vitro antiviral effects of berberine on influenza virus

In vitro antifungal effects of berberine against Candida spp

Improving the physical realism and structural accuracy of protein models by a two-step atomic-level energy minimization

Ocimum flavone Orientin as a countermeasure for thrombocytopenia

Antimicrobial activity of berberine alone and in combination with ampicillin or oxacillin against methicillin-resistant Staphylococcus aureus

a-Ketoamides as broad-spectrum inhibitors of coronavirus and enterovirus replication: Structure-based design, synthesis, and activity assessment

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved a-ketoamide inhibitors

Antidiarrheal and anti-inflammatory effects of berberine

Coronaviruses-drug discovery and therapeutic options