key: cord-0963162-te7uh3sb
authors: Hiremath, Shridhar; Kumar, H. D. Vinay; Nandan, M.; Mantesh, M.; Shankarappa, K. S.; Venkataravanappa, V.; Basha, C. R. Jahir; Reddy, C. N. Lakshminarayana
title: 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
date: 2021-01-11
journal: 3 Biotech
DOI: 10.1007/s13205-020-02578-7
sha: d2b8865e3efa922d7f3e9ee870aab57f5ce34119
doc_id: 963162
cord_uid: te7uh3sb

The Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has resulted in outbreak of global pandemic, fatal pneumonia in human referred as Coronavirus Disease-2019 (Covid-19). Ayurveda, the age old practice of treating human ailments in India, can be considered against SARS-CoV-2. Attempt was made to provide preliminary evidences for interaction of 35 phytochemicals from two plants (Phyllanthus amarus and Andrographis paniculata used in Ayurveda) with SARS-CoV-2 proteins (open & closed state S protein, 3CLpro, PLpro and RdRp) through in silico docking analysis. The nucleotide analogue remdesivir, being used in treatment of SARS-CoV-2, was used as a positive control. The results revealed that 18 phytochemicals from P. amarus and 14 phytochemicals from A. paniculata shown binding energy affinity/dock score < − 6.0 kcal/mol, which is considered as minimum threshold for any compound to be used for drug development. Phytochemicals used for docking studies in the current study from P. amarus and A. paniculata showed binding affinity up to − 9.10 kcal/mol and − 10.60 kcal/mol, respectively. There was no significant difference in the binding affinities of these compounds with closed and open state S protein. Further, flavonoids (astragalin, kaempferol, quercetin, quercetin-3-O-glucoside and quercetin) and tannins (corilagin, furosin and geraniin) present in P. amarus have shown more binding affinity (up to − 10.60 kcal/mol) than remdesivir (up to − 9.50 kcal/mol). The pharmacokinetic predictions suggest that compounds from the two plants species studied in the current study are found to be non-carcinogenic, water soluble and biologically safe. The phytochemicals present in the extracts of P. amarus and A. paniculata might have synergistic effect with action on multiple target sites of SARS-CoV-2. The information generated here might serve as preliminary evidence for anti SARS-CoV-2 activity of phytochemicals present from P. amarus and A. paniculata and the potential of Ayurveda medicine in combating the virus. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-020-02578-7.

The outbreak of highly transmissible fatal pneumonia in human referred as Coronavirus Disease-2019 (Covid-19) has spread globally (215 countries and territories) with over 16.81 M confirmed cases and over 662,095 deaths worldwide as on 31 July 2020 (https ://www.world omete rs.info/coron aviru s), whereas, in India, 15, 83,792 people with 34,968 deaths (https ://www.world omete rs.info/coron aviru s). According to the International Committee on Taxonomy of Viruses (ICTV), the novel zoonotic pathogenic virus associated with COVID-19 is Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the virus species is severe acute respiratory syndrome-related coronavirus (Gorbalenya et al. 2020 ). The break out of SARS-CoV-2 has been first reported during December 2019 in Wuhan, Hubei Province, China, with a fatality rate of 4 percent Wu et al. 2020; Zhu et al. 2020) . Sore throat, dry cough, fever, chills, repeated shaking with chills, muscle pain, headache, loss of taste or smell and respiratory 1 3 44 Page 2 of 18 distress are major symptoms induced by SARS-CoV-2 (Cheng et al. 2007 ). In most severe cases, the clinical presentation of SARS-CoV-2 infection includes pneumonia that progresses to multi-organ failure and even death (Zaki et al. 2012; Machhi et al. 2020) . Covid-19 has become global pandemic and created the wide spread fear of death everywhere in the world resulting lockdown in human society (Mandal 2020) . The increasing pandemic potential of SARS-CoV-2 is due to its high human-to-human transmissible efficiency which makes it difficult to contain (Wu et al. 2020; Zhao et al. 2020) .

The SARS-CoV-2 is a β-coronavirus, belongs to family Coronaviridae of order Nidovirales ). It has enveloped positive sense single-stranded RNA genome with 33 kilo base (Kb) encoding for 14 open reading frames (ORFs) (Fehr and Perlman 2015; Guo et al. 2020) . The genome of SARS-CoV-2 lacks the hemagglutinin-esterase gene. However, it comprises two flanking untranslated regions (UTRs) at 5′ end and 3′ with 265 and 358 nucleotides, respectively. The first ORF comprises approximately 67 percent of the genome and encodes 16 non-structural proteins (nsps), which were processed by proteolytic cleavage, whereas the remaining 13 ORFs encode for structural proteins and protein accessory (Boopathi et al. 2020; Chan et al. 2020) . The nsps mainly include helicase, RNA-dependent RNA polymerase (RdRp), 3-Chymotrypsin-like protease (3CLpro), Papain-like protease (PLpro) and others likely to be involved in the transcription and replication of SARS-CoV-2 (Simmons et al. 2005; Boopathi et al. 2020) . Four major structural proteins are spike surface glycoprotein (S protein), membrane protein (M), nucleocapsid protein (N) and envelope (E). M protein possess N-terminal glycosylated ectodomain at the N-terminal end that consists of three transmembrane domains (TM) and a long C-terminal CT domain (Boopathi et al. 2020; Chan et al. 2020) . The M and E proteins are required for virus morphogenesis, assembly and budding, whereas S protein is a fusion viral protein comprising two subunits S1 and S2. The S1 subunit consists of receptor-binding domain (RBD), N-terminal domain (NTD) and signal peptide. The S2 subunit comprises two heptad repeat regions known as HR-N and HR-C, which form the coiled-coil structures surrounded by the protein ectodomain (Walls et al. 2020) . A furin cleavage site (PRRARS'V) is present in S protein at the interface between S1 and S2 subunits that will be processed during the biogenesis (Coutard et al. 2020) . SARS-CoV-2 S protein facilitates the entry of viral particles into the host cell after binding with the host angiotensin-converting enzyme 2 (ACE-2) (Tortorici and Veesler 2019) and also determines the host range and cellular tropism (Simmons et al. 2005) . Entry of SARS-CoV-2 into host cells is mediated by the cleavage of viral spike protein by host proteases, viz. serine proteases (TMPRSS2) or cysteine proteases (Cathepsin B/L) (Boopathi et al. 2020; Yan et al. 2020) . The ACE2 and TMPRSS2 proteins are prime targets of the coronavirus (Matsuyama et al. 2010; Glowacka et al. 2011; Zhao et al. 2020) .

The highly conserved spike protein (S protein) of SARS-CoV-2 enters human cells by targeting the host ACE2. Hence, S protein is one of the drug target to neutralize the virus (Walker and Burton 2018; Karakus et al. 2020) . Proteases (3CLpro and PLpro) translated from viral RNA are essential for polyprotein processing of SARS-CoV-2 and involved in replication and transcription of SARS-CoV-2 can serve as potent therapeutic targets (Prussia et al. 2011; Sheahan et al. 2020) . The SARS-CoV-2 proteases act on polyproteins (PP1a and PP1ab) to release the nsps which are encoded by open reading frame (ORF1a/b) of virus (Sheahan et al. 2020) . Inhibiting virus entry into the cell and its replication combined with modulating the immune system could be a potential target for drug therapy (Sanders et al. 2020) . Currently, there are no available vaccines or specific medicines for the treatment of COVID-19.

Several countries and pharma giants in the world are working to come out with either vaccines or drugs which will be helpful in containing this global pandemic (Begum et al. 2020) . In India, attempts were made to address this problem through Ayurveda and Siddha practices as well apart from development of drugs and vaccines, which were still widely used among the Indian population to combat diseases caused by viral infections such as jaundice, malaria, cholera, etc. (Gomathi et al. 2020) . Ayurveda is the age old practice of treating human ailments in India and provided remedies against many medical conditions. Medicinal herbs from India are promising in treating various illnesses (Gomathi et al. 2020) . Further, identification of phytochemicals from medicinal plants can be used as drugs for viral infection targeting viral receptors (Chang and Woo 2003; Keyaerts et al. 2004 ), viral integration (Kim et al. 2010) , reverse transcription, viral replication and viral protein translation (Mansouri et al. 2009 ).

Extracts from Andrographis paniculata Wall (Kalmegh) and Phyllanthus amarus Schum. & Thonn were used in Ayurvedic medicinal systems and proven to have numerous medicinal properties (Niranjan et al. 2010; Patel et al. 2011) .

A. paniculata also known as 'King of bitters' belonging to Acanthaceae family (Pholphana et al. 2013 ) is widely used to treat fever, diarrhoea and common cold. It consists of a wide range of pharmacologically active substances such as andrographolides and its derivatives that have been showed anticancer, antimalarial, hepatoprotective, antioxidant activities (Niranjan et al. 2010 ) and reported to have potent antiviral activity against diverse group of viruses belonging to different families (Gupta et al. 2017) . It also been implicated in suppression of increased NOD-like receptor protein 3 (NLRP3), caspase-1 and interleukin-1β molecules which are extensively involved in the pathogenesis of SARS-COV and likely SARS-CoV-2 as well (Liu et al. 2020a, b) . Phyllanthus amarus Schum. & Thonn. (Family: Euphorbiaceae) is one of the medicinal herb widely distributed in tropical and subtropical countries from Africa to Asia, South America and the West Indies. It is a rich source of secondary metabolites such as alkaloids, flavonoids, hydrolysable tannins, lignins, polyphenols, triterpenes, sterols and volatile oils (Patel et al. 2011) . It exhibits an array of ethanopharmacological activities such as anti-inflammatory, hepatoprotective, nephroprotective, anti-amnesia, anticancer, diuretic, antioxidant, antiviral, anti-bacterial, anti-hyperglycemic, anti-hypercholesterolemia and so on (Calixto et al. 1998; Patel et al. 2011) .

Several Ayurvedic practitioners in India claiming Ayurveda medicine can cure Covid-19, which are not supported by enough scientific evidences. Nevertheless, the role of herbal extracts in addressing the human ailments cannot be ignored totally. Lack of proper experimental evidences might have made regulatory authorities to enforce not to sale the Ayurveda medicines as a curative agent of Covid-19. There are several Ayurveda tablets available in the Indian market claiming to boost immune system in the backdrop of Covid-19 outbreak. Here we discussed on the potential Ayurveda medicine prepared using A. paniculata and P. amarus against SARS-CoV-2 through the results of in silico analysis of interaction between SARS-CoV-2 proteins and secondary metabolites of these plants available in the databases. This provides some preliminary scientific evidence to harness the potential of Ayurveda medicines in combating SARS-CoV-2. However, this needs to be to be validated with experimental evidences required by the regulatory authorities such as large-scale clinical trials.

The list of phytochemicals present in P. amarus and A. paniculata was obtained from previously published scientific literature (Chao and Lin 2010; Gupta and Vaghela 2019; Joseph and Raj 2011; Nisar et al. 2018; Pereira et al. 2016; Tan et al. 2016) . A total of 35 compounds were selected for the molecular docking studies (21 from P. amarus and 14 from A. paniculata) (Table 1) . Remdesivir, an antagonist of SARS-CoV-2 RdRp, was chosen as a positive control (Yu et al. 2020) .

The three-dimensional (3D) conformers of the selected ligands were retrieved from PubChem (https ://pubch em.ncbi.nlm.nih.gov) database in SDF and MOL formats. Similarly, the 3D structures of proteins determined using a biophysical technique such as X-ray crystallography or NMR spectroscopy are openly available in protein data bank were retrieved (www.rcsb.org). Four different viral targets were considered for docking studies, viz. spike protein (S protein), 3-chymotrypsin-like protease (3CLpro), papain-like protease (PLpro) and RNA-dependent RNA polymerase (RdRp). The crystal structure of S protein and 3CLpro of SARS-CoV-2 available in the protein databank were retrieved for the current study. The PDB ID for S protein and 3CLpro are PDB ID: 6VSB, chains A, B, C, 6VXX (closed), 6VYB (open) (Wrapp et al. 2020 ) and PDB ID: 6LU7, chain A , respectively. Crystal protein structures of PLpro and RdRp of SARS-CoV-2 were not available in the protein databank at the time of analysis. Hence, the crystal structure developed for the both proteins using the X-ray structure of PLpro of human coronavirus papain-like proteases ( 

The retrieved 3D structures of the ligands were optimized using Avogadro version 1.2.0 (Hanwell et al. 2012 ) with force filed type MMFF94 and converted to PDB format using Open Babel version 3.1.1. (O'Boyle et al. 2011) . For further simplification of analysis, the optimized ligands with the lowest energy were loaded to AutoDock-MGL Tools (Morris et al. 2009 ), the Gasteiger charges were added and standard processes were used to obtain the PDBQT files. The downloaded 3D structures were checked for the presence of any improper bonds, side-chain anomalies and missing hydrogens using PyMOL software (Delano 2009 ). All the water molecules, complex molecules, ions and ligands of the proteins were removed in Biovia Discovery Studio 2020 (Biovia 2020) . The optimized PDB structures were uploaded to AutoDock-MGL Tool (Morris et al. 2009 ), succinctly, polar hydrogens were added and standard processes were used to obtain the PDBQT files.

The active sites of the proteins were determined using the Computed Atlas for Surface Topography of Proteins (CASTp) (http://sts.bioe.uic.edu/castp /index .html?201l) and Biovia Discovery Studio 2020. For closed and open conformational state of spike protein, the active sites were chosen based on the receptor-binding domain, which interacts with the ACE2 receptors on the host cell membrane of humans (Gur et al. 2020) . In case of 3CLpro protein, the binding site was selected based on the location of co-crystallized inhibitor N3, which was developed using computer-aided drug design and has shown specific inhibition to SARS and MERS with good binding affinity in the binding pockets of crystal structure of 3CLpro ). The optimized Similarly, 2D structure was visualized in Biovia Discovery studio 2020. The 3D visualization represents the exact location or the binding pocket of the target protein, whereas the 2D structure visualization represents the different bonds formed between the amino acid residues of the viral target protein and the ligand. The workflow of in silico molecular docking analysis is represented in Fig. 1 .

The drug-likeness analysis was carried out based on Lipinski's rule of five (Lipinski 2004) . Absorption, distribution, metabolism and excretion (ADME) are pharmacokinetic properties of a ligand that deals primarily with its absorption, distribution, metabolism and excretion in the human body. The major parameter, Lipinski's rule of five was generated using SWISSADME (Daina et al. 2017 ) for evaluating the pharmacological significance. The ligands in PDB format were converted to canonical SMILES format and uploaded to SWISSADME prediction (http://www.swiss adme.ch). Properties such as solubility (logS), human intestinal absorption (HIA) and carcinogen were obtained from admetSAR (http://lmmd.ecust .edu.cn/admet sar2).

Impact of Covid-19 is having catastrophic effect on human health and life. The scientific communities from the nuke and corner of the world have scrupulously engaged in developing the drug or vaccine against SARS-CoV-2, causal agent of Covid-19. The current study discusses on the potential of phytochemicals present in two medicinal plants used in the Ayurveda, an age old medical practice in India. Plants/herbal extracts part of Ayurveda medicines are known to boost immunity in humans, which is crucial for fighting against microbial enemies The docking score of respective phytochemicals higher than the positive control (remdesivir) has been denoted in bold *The phytochemicals showing better binding affinity with all the viral proteins in A. paniculata were used for molecular docking simulation with four SARS-CoV-2 receptor targets.

Prior to the docking analysis, the ligands were optimized by minimizing the energy with force field type MMFF94, this helps in removing clashes among atoms and develops a stable starting pose of the ligands for binding interaction (Klebe et al. 2006) . The docking coupled with a scoring function can be utilized to evaluate a large number of potential drugs in silico. In the process of molecular docking, generally, the binding affinity lesser than the upper threshold (− 6 kcal/mol) is regarded as a cut off value for determining good binding affinity between ligand and protein (Shityakov et al. 2014) . Of the 35 compounds screened from two plant species (Table 1) , three compounds (cinnamic acid, phenylalanine and farnesyl farnesol) did not qualify as per the threshold binding affinity to neutralize SARS-CoV-2 proteins.

In the recent clinical trials, remdesivir, a nucleotide analogue of adenosine triphosphate has been reported to be a potential prodrug against SARS-CoV-2 infection (Gao et al. 2020) . It was originally developed for treating viral infection caused by Ebola and Marburg virus but was also found effective against SARS and MERS (Elfiky 2020). Hence, remdesivir has been used as a positive control in the current study. During virus replication, the RdRp mistakenly incorporates remdesivir into its genome, which results in stalling of replication . In our study, remdesivir has shown good interaction with RdRp protein (Table 1) with the binding affinity of − 7.6 kcal/mol. Similarly, it has shown a good binding affinity with 3CLpro (− 7.90 kcal/ mol), S protein (− 8.00 kcal/mol) and PLpro (− 9.50 kcal/ mol) (Table 1 ) and might be having some role in inhibiting SARS-CoV-2. Many recent computational investigations also showed inhibition activity of remdesivir with various viral receptors of SARS-CoV-2 (Elfiky 2020; Hall and Feng 2020; Murugan et al. 2020) . This shows that the prodrug might be effective against all the four viral receptors.

The 3D and 2D visualization of top three phytochemicals along with positive control remdesivir based on binding affinity with respective viral receptor protein is represented in Figs. 2, 3 , respectively. Hydrogen bond energy is the major contributor in dock score. Hydrogen bond formation between ligand structures and viral receptors is responsible for inhibition of the target protein, and it reflects the firmness of bonding between the protein and ligand (Chen et al. 2016) . Selected phytochemicals from P. amarus and A. paniculata used in the current docking studies against SARS-CoV-2 proteins shown very good dock score above the threshold cut-off of − 6 kcal/mol. Ligand structures and necessary hydrogen bond formation between phytochemical with respective viral receptor is illustrated in Table 2 and supplementary Fig. 1.   Fig. 1 Workflow of molecular docking analysis of phytochemicals present in P. amarus and A. paniculata with four proteins (S protein, 3CLpro, PLpro and RdRp) targets of SARS-CoV-2. The phytochemicals of the two plant species were obtained from PubChem (https :// pubch em.ncbi.nlm.nih.gov) database in SDF and MOL formats. The 3D structures of proteins available in protein data bank were retrieved (www.rcsb.org). The detailed protocol is given in the material and methods

S protein of the SARS-CoV-2 latches onto human angiotensin-converting enzyme 2 (ACE2) receptor and undergoes structural changes that allow the viral membrane to fuse with the cell membrane of the host (Walls et al. 2020) . Despite the high similarity between sequences and structures of SARS and SARS CoV-2, the drugs successfully used to cure SARS are ineffective in case of this novel virus (Fani et al. 2020) . This indicates that more number of potential molecules should be screened to find an effective molecule.

The molecular docking analysis of S protein of the SARS-CoV-2 with 21 selected phytochemicals from P. amarus showed that the majority of the compounds bind to S protein with favourable binding energy ranging from − 6.10 kcal/mol (for phyllanthin) to − 10.60 kcal/mol (for geraniin) (Table 1) . When the binding energy of phytochemicals was investigated, nine compounds, viz. geraniin (− 10.60 kcal/mol), furosin (− 10.20 kcal/mol), astragalin (− 9.80 kcal/mol), quercetin-3-O-glucoside (− 9.80 kcal/ mol), quercitrin (− 9.20 kcal/mol), corilagin (− 8.80 kcal/ mol), hinokinin (− 8.40 kcal/mol), kaempferol (− 8.30 kcal/ mol) and quercetin (− 8.10 kcal/mol), exhibited highest binding affinity. Interestingly, the phytochemicals belonging to the flavonoids and tannins of P. amarus have shown a very high affinity with the target protein. Molecular dynamic simulation studies also reported that quercetin and kaempferol bind with the ACE2-S complex with low binding free energy (Pandey et al. 2020) . The antiviral bioactivities of many flavonoids from medicinal plants were reported by several studies (Thayil et al. 2016; Zakaryan et al. 2017; Jo et al. 2020 ). Among all compounds, geraniin exhibited (− 10.60 kcal/mol) highest docking score with S protein. THR302, ILE312, GLN314 and ARG765 amino acid residues were involved in forming four hydrogen bonds in the binding pocket of S protein (Table 2, Fig. 3) .

In binding interaction between phytochemical compounds (14) from A. paniculata with S protein of the SARS-CoV-2, the docking score ranged from − 7.30 (for andrograpanin) to − 9.10 kcal/mol (for isoandrographolide) (Table 1) . Of the 14 compounds, five compounds: isoandrographolide (− 9.10 kcal/mol), 19-O-acetyl-14-deoxy-11, 12-didehydroandrographolide (− 8.40 kcal/mol), neoandrographolide (− 8.20 kcal/mol) and 5-hydroxy-7, 8, 2′, 5′-tetramethoxy flavone (− 8.20 kcal/mol) shown better binding affinity as compared to remdesivir (− 8.00 kcal/mol). Isoandrographolide has properly positioned into the binding pocket constructed by hydrogen bonding with GLU773 and ARG1019 amino acids of S protein exhibited highest (− 9.10 kcal/mol) binding affinity (Table 2, Fig. 2) . Totally, 13 compounds from P. amarus and A. paniculata have shown better inhibitory activity against S protein.

The S protein undergoes conformational changes to produce closed and open state based on down and up positions, respectively, in the receptor-binding domain. Hence, docking analysis was carried out in both the conformational states with 21 selected phytochemicals from P. amarus and 14 phytochemicals from A. paniculata. There was no significant difference in the binding affinities of the compounds with both closed and open state S protein. The range of binding affinities of 21 selected phytochemicals from P. amarus with closed state S protein of SARS-CoV-2 was − 6.60 (for melatonin) to − 10.10 kcal/mol (for geraniin). Seven compounds, geraniin (− 10.10 kcal/mol), furosin (− 9.60 kcal/ mol), corilagin (− 9.60 kcal/mol), quercetin (− 8.80 kcal/ mol), quercetin-3-O-glucoside (− 8.60 kcal/mol), quercitrin (− 8.60 kcal/mol) and hinokinin (− 8.50 kcal/mol) have shown better binding affinity than remdesivir (− 8.40 kcal/ mol). Among these, geraniin with maximum binding affinity of − 10.10 kcal/mol has formed two hydrogen bonds with CYS291, SER750 amino acid residues of closed state S protein.

In case of open state S protein, the dock score of 21 phytochemicals from P. amarus ranged between − 6.10 kcal/ mol (for niranthin) and − 10.10 kcal/mol (for geraniin). Six compounds, geraniin (− 10.10 kcal/mol), corilagin (− 9.50 kcal/mol), furosin (− 9.40 kcal/mol), quercetin (− 8.80 kcal/mol), quercetin-3-O-glucoside (− 8.60 kcal/ mol) and quercitrin (− 8.60 kcal/mol) have shown favourable binding affinity than remdesivir (− 8.50 kcal/mol). Geraniin with maximum dock score of − 10.10 kcal/mol form five hydrogen bonds with amino acid residues, GLU988, ARG995, TYR756, THR998 of closed state S protein.

All the 15 phytochemicals of A. paniculata docked with both closed and open state S protein falls between − 6.00 kcal/mol and − 8.4 kcal/mol for closed state and − 8.50 kcal/mol (remdesivir). However, these phytochemicals might also have neutralizing effects on the activity of S protein since their dock score is more than − 6.00 kcal/mol, which is a cut-off value for considering the effective binding affinity. The binding affinities ranged from − 7.10 kcal/ mol (for andrograpanin) to − 8.20 kcal/mol (for14-deoxy-14,15-didehydroandrographolide and 14-deoxyandrographolide) for closed state of S protein. For the open state S protein, dock score ranged between − 6.60 kcal/mol (for andrograpanin) and − 8.40 kcal/mol (for neoandrographolide 

The SARS-CoV-2 protease enzyme, 3CLpro, involved in C-terminal cleavage of virus polyprotein is a potential target for drug discovery to SARS-CoV-2, since it controls replication (Qmar et al. 2020). Hence, this potential target was evaluated with the selected 35 compounds from P. amarus and A. paniculata. Before docking simulation, the target site for viral receptor protein was selected based on the location of co-crystallized inhibitor N3 as this target site has been considered to be the most commending site for ligand interaction (Murugan et al. 2020) .

The binding affinity of compounds to 3CLpro ranged from − 6.10 kcal/mol (for andrograpanin) to − 9.30 kcal/mol (for geraniin) (Table 1) . Among 21 phytochemicals derived from P. amarus, six compounds, viz. geraniin (− 9.30 kcal/ mol), furosin (− 8.70 kcal/mol), corilagin (− 8.70 kcal/ mol), quercitrin (− 8.60 kcal/mol), astragalin (− 8.40 kcal/ mol) and quercetin-3-O-glucoside (− 8.20 kcal/mol), have shown better binding affinity as compared to remdesivir (− 7.90 kcal/mol) ( Table 1 ). The compounds belong to tannins (geraniin, furosin and corilagin) and flavonoids (quercitrin, astragalin and quercetin-3-O-glucoside) are found to have good inhibitory action on 3CLpro. Several reports have been shown antiviral effect of flavonoids (Gescher et al. 2011; Chen et al. 2012) . Recently, Adem et al. (2020) evaluated 80 flavonoids compounds present in different medicinal plants against 3CLpro of SARS-CoV-2. Among them, hesperidin, rutin, diosmin, apiin, diacetylcurcumin, (E)-1-(2-Hydroxy-4-methoxyphenyl)-3-[3-[(E)-3-(2-hydroxy-4methoxyphenyl)-3-oxoprop-1-enyl] phenyl] prop-2-en-1-one and beta, beta'-(4-Methoxy-1,3-phenylene) bis(2′-hydroxy-4′,6′-dimethoxyacrylophenone have shown effective binding affinity with SARS-CoV-2 3Clpro. The flavonoids from the two plant species used in the present study have higher binding affinity to protease similar results were reported by Adem et al. (2020) . Geraniin, with the highest binding affinity of − 9.30 kcal/mol, has formed five hydrogen bonds with the amino acid residues, LEU141, GLY143, SER144, CYS145 and THR190 of 3CLpro (Table 2, Fig. 3) .

The docking score for compounds present in A. paniculata ranged from − 6.10 kcal/mol (for andrograpanin) to − 7.80 kcal/mol (for neoandrographolide) (Table 1) which is on par with the remdesivir showing the score of − 7.90 kcal/ mol. During replication of the viral genome, the 3CLpro cleaves the C-terminal of polypeptide protein at 11 distinct sites to generate various nsps (Qmar et al. 2020) . Hence, the phytochemicals present in the extracts of P. amarus might hamper the activity of the protease enzyme there by inhibiting the process of virus multiplication.

The SARS-CoV-2 another protease enzyme, PLpro known to be involved in cleaving of N-terminal of viral polyprotein which is also an important antiviral drug target. It also plays a crucial role in its replication along with 3CLpro (Baez-Santos et al. 2014) . The 3D structure of PLpro used in the docking simulation is the x-ray structure of SARS-CoV, which has the highest similarity with the PLpro of SARS-CoV-2.

The binding energy of phytochemicals present in P. amarus ranged from -6.30 kcal/mol (for cinnamic acid) to − 10.30 kcal/mol (for quercetin-3-O-glucoside) ( Table 1) . Four compounds; quercetin-3-O-glucoside (− 10.30 kcal/ mol), hinokinin (− 9.80 kcal/mol), astragalin (− 9.70 kcal/ mol) and kaempferol (− 9.60 kcal/mol) from P. amarus have exhibited better binding energies compared to remdesivir (− 9.50 kcal/mol). Apart from this, the binding affinities of furosin (− 9.20 kcal/mol) and geraniin (− 9.20 kcal/mol) are reasonably well with that of the positive control, remdesivir. Quercetin-3-O-glucoside has the maximum dock score of − 10.30 kcal/mol and shown its binding to HIS74, ARG83, TYR155, ASN157 and HIS176 amino acid residues of viral protein at the binding pocket of PLpro with eight hydrogen bonds (Table 2, Fig. 3) .

The binding affinity of compounds from A. paniculata with PLpro ranged from -6.80 kcal/mol (for andrograpanin) to − 9.00 kcal/mol (for 5-Hydroxy-7, 8-dimethoxyflavone). Even though binding affinities of all the 15 compounds falls between -6.00 kcal/mol and less than the dock score of positive control remdesivir (− 9.5 kcal/mol), they also might have inhibitory effects on the activity of PLpro. The results suggest that, among the promising phytochemicals from P. amarus and A. paniculata docked with PLpro, flavonoids (quercetin-3-O-glucoside and astragalin) and phenylpropanoid (hinokinin) have exhibited better binding affinity. The SARS-CoV-2 PLpro is an essential for virus replication (Baez-Santos et al. 2014 ) and any compound having good binding affinity with inhibitory effects will be of immense potential in containing this virus. The flavonoids and phenylpropanoids showing good affinity binding abilities might serve as lead molecules for the discovery of novel drugs to SARS-CoV-2.

The RdRp plays a major role in the replication and transcription of SARS-CoV-2 genome (Aftab et al. 2020 ). The RdRp used in the current study is the solved structure of SARS-CoV having 97.08 percent sequence identity with the RdRp of SARS-CoV-2 (Kirchdoeffer and Ward 2020). Among 21 compounds from P. amarus screened against RdRp, epibubbialine has shown lowest dock score of − 6.30 kcal/mol (below − 6.0 kcal/mol is considered as having good binding affinity) and corilagin has shown highest dock score of − 9.30 kcal/mol (Table 1) . Seven compounds: corilagin (− 9.30 kcal/mol), furosin (− 9.10 kcal/mol), geraniin (− 9.10 kcal/mol), quercetin (− 8.60 kcal/mol), hinokinin (− 8.20 kcal/mol), quercetin (− 8.10 kcal/mol) and quercetin-3-O-glucoside (− 7.70 kcal/mol) from P. amarus exhibited better dock score compared to the RdRp antagonist remdesivir (− 7.60 kcal/mol). Corilagin with high binding affinity of − 9.30 kcal/mol has interacted with active site residues, TYR129, HIS133, LEU207, ASN705 and GLN724 of RdRp of virus through five hydrogen bonds (Table 2, Fig. 3 ).

Five compounds, neoandrographolide (-8.50 kcal/ mol), isoandrog rapholide (− 8.30 kcal/mol), 14-deoxy-14,15-didehydroandrographolide (− 7.90 kcal/ mol), 14-Deoxyandrographolide (− 7.80 kcal/mol) and 5-hydroxy-7, 8, 2′,5′-tetramethoxy flavone (− 7.80 kcal/ mol) from A. paniculata, have shown better binding energy affinity compared to remdesivir (− 7.60 kcal/mol) (Table 1) . Among these, neoandrographolide exhibited higher dock score of − 8.50 kcal/mol and interacted with ALA125, ASP126, ARG132 and HIS133 amino acid residues of RdRp through four hydrogen bonds (Table 2, Fig. 3 ).

In general, the flavonoids (astragalin, kaempferol, quercetin, quercetin-3-O-glucoside and quercetin), a phenylpropanoid (hinkonin) and tannins (corilagin, furosin and geraniin) from P. amarus and a bioactive compound (isoandrographolide), diterpenes ,12 didehydroandrographolide, neoandrographolide) and a flavonoid (5-hydroxy-7,8,2′,5′tetramethoxy flavone) from A. paniculata have shown good binding affinities with all the four SARS-CoV-2 target proteins. The natural compounds obtained from the medicinal plants have been used in case of two previous coronavirus outbreaks of MERS-CoV and SARS-CoV (Li et al. 2005; Lin et al. 2005) . Similarly, strong interaction of tannins (polyphenols) from green tea, epigallocatechin gallate, epicatechingallate and gallocatechin-3-gallate with viral receptor (3CLpro) of SARS-Cov-2 was shown (Ghosh et al. 2020) . The flavonoids and tannins interaction with proteins of SARS-CoV-2 in vitro may provide more insights into exploiting these compounds as potential SARS-CoV-2 inhibiting compounds and lead to novel drug development. Further, among all the phytochemicals, the geraniin has shown consistent interaction and very good binding affinity with all the viral protein receptors. The geraniin, extracted from P. amarus, has shown to inhibit HIV-1 replication in HeLa CD4 + cells with 50% effective concentration (EC50) values ranging from 0.9 to 7.6 μg/ml (Notka et al. 2004 ). This suggests that geraniin with multiple target binding abilities in docking analysis with high dock score than remdesivir can be a potential candidate to harness as anti SARS-CoV-2 drug. However, this needs in vitro evaluation to authenticate the prediction studies.

The pharmacokinetic evaluation and ADME (absorption, distribution, metabolism and excretion) properties predictions of pharmacological compounds in advance are important methods to save huge monetary expenses in drug discovery (Zhu et al. 2011) . Apart from that, it mainly helps in reducing the potential risk during clinical development as it provides very crucial information to determine whether a compound is suitable to proceed to clinical stage. SWISS-ADME (http://www.swiss adme.ch) and admetSAR (http:// lmmd.ecust .edu.cn/admet sar2) web servers were used to obtain the ADME and other pharmacokinetic properties of the selected compounds. Among the 35 compounds from P. amarus and A. paniculata evaluated against SARS-CoV-2 by in silico docking analysis, 15 compounds (9 compounds from P. amarus and 6 compounds from A. paniculata) were selected for ADME and other pharmacokinetic properties predictions. All these 15 compounds selected showed better binding ability than the remdesivir. The ADME properties of the selected 15 compounds have been represented in supplementary Fig. 2 . The ADME and pharmacokinetic properties of selected phytochemicals are given in Table 3 . Lipinski rule of five (Lipinski 2004 ) is a thumb rule to evaluate the drug likeliness or to determine the pharmacological activities that would make it orally active drug in humans. According to the rule, the compound should abide following criteria, i.e., (1) molecular mass less than 500 Daltons, (2) high lipophilicity (expressed as logP less than 5), (3) less than five hydrogen bond donors and (4) less than ten hydrogen bond acceptors (Petit et al. 2012) . Of 15 compounds along with remdesivir screened, eight compounds (kaempferol, quercetin, hinokinin, isoandrographolide, 14-deoxy-14, 15-didehydroandrographolide, 14-Deoxyandrographolide, 19-O-Acetyl-14-deoxy-11, 12-didehydroandrographolide and neoandrographolide) did not violate any of the above rule. However, four compounds, viz. astragalin (HBA > 10, HBD > 5), quercetin-3-O-glucoside (HBA > 10, HBD > 5), quercitrin (HBA > 10, HBD > 5) and 5-hydroxy-7, 8, 2′, 5′-tetramethoxy flavone (MW > 500, HBA > 10) along with positive control remdesivir, have violated two rules. Three compounds belong to tannins, i.e., corilagin (MW > 500, HBA > 10, HBD > 5), furosin (MW > 500, HBA > 10, HBD > 5) and geraniin (MW > 500, HBA > 10, HBD > 5), have violated three rules ( Table 3 ). The violation of two or more criteria predicts a molecule to be a non-orally available drug (Lipinski 2004 ). However, this rule applies only for oral drugs, and the compounds violating the Lipinski's rule of five can be administered through other means after determining its side effects (Benet et al. 2016) .

The logS value is directly related to the solubility of the drug in water which ideally ranged between − 6.5 and 0.5 (Nisha et al. 2016) . All the screened compounds lie within this range. The solubility of the molecules in water significantly facilitates drug developmental activities (Diana et al. 2017) . The selected phytochemicals have passed the carcinogenic test and all exhibit non-carcinogenic nature (Table 3) . Among 15 phytochemicals, seven phytochemicals (kaempferol, quercetin, hinokinin, isoandrographolide, 14-deoxy-14,15-didehydroandrographolide, 14-Deoxyandrographolide and 19-O-acetyl-14-deoxy-11, 12-didehydroandrographolide) have shown high gastrointestinal or human intestinal absorption. Remaining eight phytochemicals (astragalin, quercetin-3-O-glucoside, quercetin, corilagin, furosin, geraniin, neoandropholidae and 5-hydroxy-7, 8, 2′, 5′-tetramethoxy flavone) have exhibited low gastrointestinal or human intestinal absorption. The positive control remdesivir also has low human intestinal absorption value. As all the compounds have exhibited non-carcinogenic nature and good pharmacokinetic properties these compounds individually or in combination may be effective against the SARS-CoV-2. Further, inflammation of lungs resulting in breathing problem is a major symptom in the critical patients suffering from Covid-19 (Li and Ma 2020) . In addition to all the prediction studies, the plant extracts of P. amarus have shown anti-inflammatory properties in rat Kupffer cells (Kiemer et al. 2003) . The pharmacokinetic predictions suggest that compounds from the two plants species studied in the current study are found to be biologically safe and can be utilized as therapeutics against SARS-CoV-19. However, more studies need to be conducted to decipher the exact Table 3 Absorption, distribution, metabolism and excretion (ADME) and other pharmacokinetic properties of selected phytochemicals from Phyllanthus amarus mechanism involved in the interaction between these phytochemicals and the viral receptor proteins.

To summarize, all these current studies were done by in silico predictions of individual phytochemical originated from P. amarus and A. paniculata. The information generated here may provide some insights into exploring for potential drugs against SARS-CoV-2 from plants, which can result in discovery of novel drugs. Further, the Ayurveda medications using herbal extracts cannot be undermined in addressing human ailments, because these extracts contain several compounds, which might have synergistic effect with action on multi-target sites of pathogen. However, the biggest problem in Ayurveda medicine is fixing the dose against the target pathogens. Optimization of this issue supported by scientific evidences for getting required clearances which will have far fetching impact in the harnessing Indian ancient knowledge of medication "Ayurveda" in treating the human ailments including Covid-19 caused by SARS-CoV-2, which threatened the mankind with global pandemic outbreak.

Identification of Potent COVID-19 Main protease (Mpro) inhibitors from natural polyphenols: an in silico strategy unveils a hope against Corona

Analysis of SARS-CoV-2 RNA-dependent RNA polymerase as a potential therapeutic drug target using a computational approach

X-ray structural and biological evaluation of a series of potent and highly selective inhibitors of human coronavirus papain-like proteases

Challenges and prospects of COVID-19 vaccine development based on the progress made in SARS and MERS vaccine development

BDDCS, the rule of 5 and drugability

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

A review of the plants of the genus Phyllanthus: their chemistry, pharmacology, and therapeutic potential

Genomic characterization of the 2019 novel humanpathogenic coronavirus isolated from a patient with a typical pneumonia after visiting Wuhan

Korean medicinal plants inhibiting to Human Immunodeficiency Virus type 1 (HIV-1) fusion

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

Epigallocatechin-3-gallate inhibits the replication cycle of hepatitis C virus

Regulation of protein-ligand binding affinity by hydrogen bond pairing

Severe acute respiratory syndrome coronavirus as an agent of emerging and re-emerging infection

The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade

SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

The PyMOL Molecular Graphics System

Rehabilitating the addicted brain with transcranial magnetic stimulation

Anti-HCV, nucleotide inhibitors, repurposing against COVID-19

Comparison of the COVID-2019 (SARS-CoV-2) pathogenesis with SARS-CoV and MERS-CoV infections

Coronaviruses: an overview of their replication and pathogenesis

Computational protein-ligand docking and virtual drug screening with the AutoDock suite

Structure of the RNA-dependent RNA polymerase from COVID-19 virus

Oligomeric proanthocyanidins from Rumex acetosa L. inhibit the attachment of herpes simplex virus type-1

Evaluation of green tea polyphenols as novel corona virus (SARS CoV-2) main protease (Mpro) inhibitors-an in silico docking and molecular dynamics simulation study

Drug studies on Rett syndrome: from bench to bedside

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak-an update on the status

Recent advances in pharmacological and phytochemistry studies on Phyllanthus amarus

Broad-spectrum antiviral properties of andrographolide

Conformational transition of SARS-CoV-2 spike glycoprotein between its closed and open states

A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease

Avogadro: an advanced semantic chemical editor, visualization, and analysis platform

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

Inhibition of SARS-CoV 3CL protease by flavonoids

An overview: pharmacognostic properties of Phyllanthus amarus Linn

Breaking the convention: sialoglycan variants, coreceptors, and alternative receptors for influenza a virus entry

In vitro inhibition of severe acute respiratory syndrome coronavirus by chloroquine

Phyllanthus amarus has anti-inflammatory potential by inhibition of iNOS, COX-2, and cytokines via the NF-κB pathway

Inhibition of influenza virus replication by plant-derived isoquercetin

Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors

Virtual ligand screening: strategies, perspectives and limitations

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

Acute respiratory failure in COVID-19: is it "typical" ARDS?

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

Anti-SARS coronavirus 3C-like protease effects of Isatis indigotica root and plant-derived phenolic compounds

Lead-and drug-like compounds: the rule-of-five revolution

Research and development on therapeutic agents and vaccines for COVID-19 and related human coronavirus diseases

Composition and divergence of coronavirus spike proteins and host ACE2 receptors predict potential intermediate hosts of SARS-CoV-2

The natural history, pathobiology, and clinical manifestations of SARS-CoV-2 infections

The global emergence of severe acute respiratory syndrome coronavirus 2 in human

Suppression of human T-cell leukemia virus I gene expression by pokeweed antiviral protein

Efficient activation of the severe acute respiratory syndrome coronavirus spike protein by the transmembrane protease TMPRSS2

AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility

Computational investigation on Andrographis paniculata phytochemicals to evaluate their potency against SARS-CoV-2 in comparison to known antiviral compounds in drug trials

Biological activities of kalmegh (Andrographis paniculata Nees) and its active principles: a review

Chemical components and biological activities of the genus Phyllanthus: a review of the recent literature

Molecular docking and in silico ADMET study reveals acylguanidine 7a as a potential inhibitor of β-secretase

Concerted inhibitory activities of Phyllanthus amarus on HIV replication in vitro and ex vivo

Open Babel: an open chemical toolbox

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

Phyllanthus amarus: ethnomedicinal uses, phytochemistry and pharmacology: a review

Extraction of lignans from Phyllanthus amarus Schum. & Thonn using pressurized liquids and low pressure methods

Softening the rule of five-where to draw the line

Changes in the contents of four active diterpenoids at different growth stages in Andrographis paniculata

Systematic approaches towards the development of host-directed antiviral therapeutics

Structural basis of SARS-CoV-2 3CLpro and anti-COVID-19 drug discovery from medicinal plants

Pharmacologic treatments for coronavirus disease 2019 (COVID-19): a review

Comparative therapeutic efficacy of remdesivir and combination lopinavir, ritonavir, and interferon beta against MERS-CoV

Ionization states, cellular toxicity and molecular modeling studies of midazolam complexed with trimethyl-β-cyclodextrin

Inhibitors of cathepsin L prevent severe acute respiratory syndrome coronavirus entry

Chemical constituents of Andrographis paniculata

Pa-9: a flavonoid extracted from Plectranthus amboinicus inhibits HIV-1 protease

Structural insights into coronavirus entry

AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

Passive immunotherapy of viral infections, 'super-antibodies' enter the fray

Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein

Review of the 2019 novel coronavirus (SARS-CoV-2) based on current evidence

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

A new coronavirus associated with human respiratory disease in China

Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2

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

Flavonoids: promising natural compounds against viral infections

Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia

Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCov

Recent developments of in silico predictions of oral bioavailability

A novel coronavirus from patients with pneumonia in China