key: cord-1009075-xjvwz1ct
authors: Pandey, Preeti; Rane, Jitendra Subhash; Chatterjee, Aroni; Kumar, Abhijeet; Khan, Rajni; Prakash, Amresh; Ray, Shashikant
title: Targeting SARS-CoV-2 spike protein of COVID-19 with naturally occurring phytochemicals: an in silico study for drug development
date: 2020-07-22
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1796811
sha: 56ec489b8504a8d0b056cb26d95909f33b26eef7
doc_id: 1009075
cord_uid: xjvwz1ct

Spike glycoprotein, a class I fusion protein harboring the surface of SARS-CoV-2 (SARS-CoV-2S), plays a seminal role in the viral infection starting from recognition of the host cell surface receptor, attachment to the fusion of the viral envelope with the host cells. Spike glycoprotein engages host Angiotensin-converting enzyme 2 (ACE2) receptors for entry into host cells, where the receptor recognition and attachment of spike glycoprotein to the ACE2 receptors is a prerequisite step and key determinant of the host cell and tissue tropism. Binding of spike glycoprotein to the ACE2 receptor triggers a cascade of structural transitions, including transition from a metastable pre-fusion to a post-fusion form, thereby allowing membrane fusion and internalization of the virus. From ancient times people have relied on naturally occurring substances like phytochemicals to fight against diseases and infection. Among these phytochemicals, flavonoids and non-flavonoids have been the active sources of different anti-microbial agents. We performed molecular docking studies using 10 potential naturally occurring compounds (flavonoids/non-flavonoids) against the SARS-CoV-2 spike protein and compared their affinity with an FDA approved repurposed drug hydroxychloroquine (HCQ). Further, our molecular dynamics (MD) simulation and energy landscape studies with fisetin, quercetin, and kamferol revealed that these molecules bind with the hACE2-S complex with low binding free energy. The study provided an indication that these molecules might have the potential to perturb the binding of hACE2-S complex. In addition, ADME analysis also suggested that these molecules consist of drug-likeness property, which may be further explored as anti-SARS-CoV-2 agents. Communicated by Ramaswamy H. Sarma

The world population is facing a severe mass annihilation due to the rise of a global pandemic named Coronavirus Disease 2019 (COVID-19) (Boopathi et al., 2020; Chatterjee et al., 2020; Joshi et al., 2020; Kirchdoerfer et al., 2016) . This pandemic is caused by a novel single-stranded RNA virus belonging to the b-coronavirus genera of the coronaviridae family (Elfiky & Azzam, 2020; Enmozhi et al., 2020; Sarma et al., 2020; Sinha et al., 2020) . As this virus shares significant phylogenetic similarity and structural familiarity (about 80% nucleotide identity and 89.10% nucleotide similarity) with the severe acute respiratory syndrome coronavirus (SARS-CoV), it has been named as Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and placed in the same lineage (Subgenus Sarbecovirus) (2020; Boopathi et al., 2020; Das et al., 2020; Ou et al., 2020; . To date, no effective regime of antivirals or vaccines is available for the use of the general public to combat the effect of COVID infections, which has put the population at a more vulnerable position (Aanouz et al., 2020; Elfiky, 2020; Elmezayen et al., 2020; Muralidharan et al., 2020; Pant et al., 2020) . The transmission of the virus is taking place at a massive rate worldwide, and hence the scientists all over the world are desperately looking for effective compounds to use it as anti-CoV therapeutic agents. Natural compounds with high bioavailability and low cytotoxicity seem to be the most efficient candidates in this regard. Since ancient times, humans have resorted to the use of natural compounds, especially phytochemicals, for the treatment of different diseases and disorders (Forni et al., 2019; Islam et al., 2020) . Flavonoids are secondary metabolites produced by plants, play vital roles in plant physiology, possessing a variety of potential biological benefits such as antioxidant, anti-inflammatory, anticancer, antibacterial, antifungal and antiviral activities (Islam et al., 2020; Kumar & Pandey, 2013; Panche et al., 2016; Wang et al., 2018) . Different flavonoids, including flavones and flavonols, have been thoroughly investigated for their potential antiviral properties. Many of them showed a significant antiviral response in both in vitro and in vivo studies Nijveldt et al., 2001) . Curcumin, a component of turmeric, has been described to exhibit enhanced antiviral activity against diverse viruses such as dengue virus (serotype 2), herpes simplex virus, human immunodeficiency virus, Zika, and chikungunya viruses among others (Mounce et al., 2017; Praditya et al., 2019) . Apigenin isolated from the sweet basil (Ocimum basilicum) plant has shown potent antiviral activity against hepatitis B virus, adenoviruses, african swine fever virus, and some RNA viruses in vitro (Chiang et al., 2005; Wu et al., 2017) . Luteolin, another acclaimed flavone, has shown significant antiviral effects on both HIV-1 reactivation and inhibition of Epstein-Barr virus (EBV) reactivation in cells (Wu et al., 2016) . Besides these antiviral activities, luteolin also showed antiviral effects against Chikungunya virus, Japanese encephalitis virus (Fan et al., 2016) , severe acute respiratory syndrome coronavirus (SARS-CoV) and rhesus rotavirus. Among flavonols, the antiviral effect of quercetin has been most extensively investigated. Quercetin has been found to demonstrate an in vitro dose-dependent antiviral activity against respiratory syncytial virus, poliovirus type 1, HSV-1, and HSV-2. It was also reported that quercetin has the potential ability as a prophylactic against Ebola virus infection (Hasan et al., 2020; Kaul et al., 1985) . Kaempferol, another flavonol extracted from Ficus benjamina leaves, has shown inhibitory activity against HCMV, HSV-1, HSV-2, and influenza A virus (Lyu et al., 2005; Zakaryan et al., 2017) . Fisetin, a modified flavonol, has shown to inhibit CHIKV infection as well as the HIV-1 infection by blocking viral entry and virus-cell fusion (Zakaryan et al., 2017) . Resveratrol and pterostibenes have been found to exhibit antiviral activities against a wide range of viruses, including HIV-1 (Abba et al., 2015; Chan et al., 2017) . Resveratrol has been shown to inhibit the replication of pseudorabies virus (PVR) and Middle East Respiratory Syndrome coronavirus (MERS-CoV) and also reduced the expression of MERS-CoV nucleocapsid (N) protein (Campagna & Rivas, 2010; Docherty et al., 2005; Lin et al., 2017; Zhao et al., 2017) . Pterostilbene, a natural dimethylated analog of resveratrol found in berries and grapes, has been found to inhibit the replication of several viruses, including herpes simplex viruses (HSVs) 1 and 2, varicella-zoster virus, influenza virus, and human papillomaviruses. Considering the contagiousness of the COVID-19 and its consequences, there is an imperative need to develop an effective therapy to curtail the spread of this deadly virus and safely treat the infected individuals. In that direction, the repurposing of the FDA approved existing drugs like chloroquine (CQ) and hydroxychloroquine (HCQ) either alone or in combination with other known drugs are currently being attempted (Al-Bari, 2017; Yao et al., 2020) . Preliminary in vitro studies and clinical trials carried out by scientists on COVID-19 patients disclosed the effectiveness of HCQ, an anti-malarial drug in combination with azithromycin, a broad spectrum anti-bacterial drug in reducing the viral load . Although some of these initial studies disclosed promising results, a lot still remains to be done to analyze their compatibility, cost, accessibility, side effects, dosages, etc. Currently, scientists are indulging themselves in identifying ideal natural compounds that can target and modulate unique or novel sites like the spike glycoprotein (S) on the surface of SARS-CoV-2 (Adedeji et al., 2013; Elfiky & Azzam, 2020; Sarma et al., 2020; Walls et al., 2020) . In this in silico study, we have performed molecular docking experiments to ascertain the most potent natural compounds (flavonoids) that can bind to the functional domains of the SARS-CoV-2 spike protein (SARS-CoV-2S), a viral surface glycoprotein required for initial attachment and internalisation within host cells (Abdelli et al., 2020; Adedeji et al., 2013; Kirchdoerfer et al., 2016) . Here we found that about 10 of these compounds effectively bind to the C-terminal region of either the S1 domain or the S2 domain of SARS-CoV-2S, and their binding interaction is more stable than with that of HCQ. These natural compounds are capable of binding to either the S1 or S2 domains of the SARS-CoV-2S protein and most probably prevent it from binding to the hACE2 receptor or internalization during fusion (Hasan et al., 2020; Zhou & Simmons, 2012) . Further, our MD simulation studies reveal that phytochemicals fisetin, quercetin, and kamferol has the propensity to bind at the junction of the hACE2-S protein complex with lesser binding energy. In addition, in-depth in vitro and in vivo with these lead molecules may help to explore as novel anti-COVID-19 agents. , and Apigenin (CID: 5280443)] was downloaded from the PubChem database, and then these structures were converted in pdb format by using PyMol (DeLano, 2002) . The structure of SARS-CoV-2S protein was downloaded from the RCSB protein data bank (PDB-ID: 6VYB) (DeLano, 2002; Walls et al., 2020) . The structure of all the ligands has been provided in Table S1 (Supplementary material).

First, the SARS-CoV-2S, and ligand were converted into pdbqt format using AutoDock tools. Polar hydrogens and gasteiger charges were added to SARS-CoV-2S and other ligand structure before docking.

The cryo-electron microscopic structure of SARS-CoV-2S was used for the molecular docking analysis. SARS-CoV-2S protein is a heterotrimer consisting of chains A, B, and C (Walls et al., 2020) . For the docking experiment, chain A of the spike protein was used. First, the SARS-CoV-2S, kamferol, and HCQ were converted into pdbqt format using AutoDock tools (Morris et al., 2009) . Polar hydrogens and gasteiger charges were added to SARS-CoV-2S, kamferol, and HCQ structure before docking. The molecular docking tool AutoDock Vina (Trott & Olson, 2010) was used to study the binding of curcumin and hydroxychloroquine on the SARS-CoV-2S. Further, blind docking of kamferol was performed to ascertain the probable binding sites. For this, the entire SARS-CoV-2S molecule was covered with a grid box of dimension 76 Å Â 92 Å Â 160 Å with grid spacing 1 Å. The SARS-CoV-2S was kept rigid while the kamferol was kept flexible. The four sets of docking were performed with exhaustiveness 100. Each set of AutoDock Vina produced 9 conformations (Trott & Olson, 2010) . Among them, 6 conformations were docked at one domain of SARS-CoV-2S, that domain is used for the direct docking. For local docking, SARS-CoV-2S was covered with the grid box of dimension 60 Å Â 54 Å Â 66 Å with grid spacing 1 Å. The exhaustiveness was kept at 100. Four sets of local docking were performed for each ligand and the site, and the site where maximum number of conformations bind was considered the binding site. The conformations with high negative binding energy are represented in the figures. The docking results were analyzed using MGL Tools 1.5.6 (Morris et al., 2009) , and the hydrophilic and hydrophobic interactions were determined using PyMol (DeLano, 2002) . Similar to kamferol, the docking of other ligand was performed on the SARS-CoV-2S protein. The docking parameters of all the ligand have been provided in Table 1 .

The biomolecular software package GROMACS-2018.1 (Kutzner et al., 2019) was used for the molecular dynamics simulations of the protein-ligand complex with the protein interactions approximated by the CHARMM36 force field and ligand parameters generated using CgenFF . For each protein-ligand complex, we place the respective complex in the center of a cubic box with 10 Å distance to the edges and fill the box with TIP3P water and counter ions (Na þ Cl À ) 0.15 M added to neutralize the system (Joung & Cheatham, 2008) . The simulation box of protein complexed with fistein, kamferol and quercontain contain water molecules, 178467, 178474 and 136669, respectively. Because of periodic boundary conditions, electrostatic interactions were evaluated using particle-Ewald summation, and a cut-off of 10 Å was used for the calculation of vdW-interactions. The resulting systems were energy-minimized by steepest descent and conjugate gradient methods, followed by a short (500 ps) equilibration in NVT ensemble and subsequent 500 ps in an NPT ensemble. Temperature and pressure were set T ¼ 310 K and 1 bar, which was controlled by a Parrinello-Danadio-Bussi thermostat and Parrinello-Rahman barostat (Bussi et al., 2007; Parrinello & Rahman, 1980) . The integration step of 2 fs was used. Each system was simulated for 100 ns, and the snapshots were saved every 10 ps for further analysis. We analyze the resulting trajectories using the tools provided by the GROMACS package (Abraham et al., 2015) . Several structural parameters that we measure and analyze include the root-meansquare deviation (RMSD), the radius of gyration (Rg), solventaccessible-surface-area (SASA), root-mean-square-fluctuation (RMSF), hydrogen bonds, hydrophobic interactions and binding free energy (Islam et al., 2020; Jain et al., 2020; Mishra et al., 2018; Sarma et al., 2020) . Hydrogen bonds were defined by a distance cut-off of 3.5 Å between the donor & acceptor atom and an angle cut-off of 30 .

Utilizing the MD trajectory, the binding free energies for ligands with S protein were calculated using MM-PBSA (Molecular Mechanics-Poisson-Boltzmann Surface Area). Considering the convergence issues, the snapshots at an interval of 10 ps were extracted from the last 20 ns of each trajectory and used for the MM-PBSA calculation (Batt et al., 2012; Prakash & Luthra, 2012) . The solute dielectric constant value of 2 and salt concentration of 0.15 M was used for the calculation (Kumar et al., 2019; Pandey et al., 2018) . Binding free energy of protein-ligand complexes were calculated as,

where, G complex represents the free energy of the protein-ligand complex, G receptor , the free energy of protein, G ligand as the free energy of ligand, and < > represent the ensemble average.

ADME profiling of fisetin, quercetin, and kamferol at pH 7 were determined using online software tools (Jayaram et al., 2012) . The important parameters allied with ADME properties such as Lipinski's rule of five, solubility of drug, pharmacokinetic properties, molar refractivity, and drug likeliness were deliberated (Lipinski, 2004) .

In the field of computer-aided drug designing, particularly for the identification of a lead compound (Hughes et al., 2011; Raj et al., 2019) , molecular docking is immensely employed to explore the various types of binding interaction of the prospective drug with different domains or active sites on the target molecules. Among all different types of interactions such as H-bond, p-p, amide-p interactions, etc., the binding efficacy of a ligand molecule with the active sites of a target has widely been explained by evaluating its hydrogen bonding pattern Raj et al., 2019) and the nature of residues present at the active site. The binding energy (kcal/mol) data allows us to study and compare the binding affinity of different ligands/compounds with their corresponding target receptor molecule. The lower binding energy indicates a higher affinity of the ligand for the receptor. The ligand with the highest affinity can be chosen as the potential drug for further studies. The SARS-CoV-2S is a surface glycoprotein of 2019 novel coronavirus (Ou et al., 2020; Wahedi et al., 2020) . This protein plays important roles during viral attachment, fusion, and entry into the host cells (Adedeji et al., 2013; Boopathi et al., 2020; . Biologically this protein exists in a heterotrimeric form with three separate polypeptide chains: chain A, B, and C, forming each monomer (Boopathi et al., 2020; Li et al., 2006) . For the present study, ten phytochemicals of flavonoids and non-flavonoids class with a broad range of biological activities along with one FDA-approved anti-malarial drug HCQ which exhibited efficacy against SARS-CoV-2 have been selected as ligands to investigate their binding affinity with SARS-CoV-2S chain A as the receptor target protein. The docking study revealed that out of the three different domains, namely S1-N terminal, S1-C terminal, and S2, these phytochemicals and HQC exhibited binding affinity mainly for two domains: S1-C terminal and S2 of the spike protein of SARS-CoV-2. Among all the ligands, 2-(3,4-dihydroxyphenyl)-3,7-dihydroxy-4Hchromen-4-one (fisetin) ( Table 2 ) 2-(3,4-dihydroxyphenyl)-3,5,7-trihydroxy-4H-chromen-4-one (quercetin) and (Table 2) displayed lowest and identical binding energy of À8.5 kcal/ mol as well as similar binding preferences for S2 domain of the spike protein. Despite having a similar preference for the S2 domain, the presence of an additional 5-OH group on the chromone ring of quercetin affected its hydrogen bonding interactions compared to fisetin, which was also reflected in their interaction with different residues of S2 domain. As shown in Figure 1C and Figure 1E and Table 2) .

Similarly, lack of one hydroxyl (-OH) group at 3 0 -C of phenyl ring in case of 3,5,7-trihydroxy-2-(4-hydroxyphenyl)-4Hchromen-4-one (kamferol) ( Figure 1F , G and Table 2 ) not only lead to the reduction in binding affinity (À7.4 kcal/mol) but it also leads to binding at S1 domain rather than S2 domain which were observed in case of other three flavonols.

Like-wise the binding studies performed with hydroxychloroquine (HCQ), which has emerged as a potent therapeutic option against COVID-19 and received much attention in the last few days, revealed the lowest binding affinity (À5.6 kcal/mol) among all others (Table 2) .

It interacts with ALA 520, GLN564, PHE 565, ARG 567, and HIS 519 residues through hydrogen bond formation. It forms the hydrophobic interaction with LEU 518, LEU 517, CYS 391, LEU 546, and ALA 522. The lowest binding affinity observed with HCQ could be attributed to the presence of a conformationally labile and sterically bulkier carbon chain, which probably disrupts the interaction of the HCQ with the binding site resulting in reduced binding affinity. Overall the docking study discloses two different sets of ligands which bind at S1 and S2 domains of SARS-CoV-2S. We thought that if these molecules interact with either S1 or S2 domain strongly, it may perturb the hACE2-S protein complex interaction. So, based on the binding energy, we selected these top phytochemicals for our further in silico study.

The molecular docking analyses resulted in the selection of three phytochemicals, fistein, kamferol, and quercetin which show promising activity against hACE2-S protein complex. However, during the biological functions, proteins undergo substantial conformational changes, and aqueous environment around the protein plays a functional role in protein-ligand interactions (Luthra et al., 2009; Mongre et al., 2019; Yu and Dalby, 2018) . Thus, to better understand the conformational stability of ligand binding and the contribution of the active site in terms of binding free energy, we performed MD simulation in solvated condition at physiological To investigate the overall conformational dynamics and stability of the protein-ligand complex, we calculated the RMSD values relative to the starting structure for the backbones of all residues (Figure 2 ). The probability distribution plot of the docked structure of hACE2-S protein with fistein shows that the population of stable conformational dynamics is settled around $18-20 Å. We find a small population density between 5 and 10 Å, which indicates the initial perturbation in the structure fistein is well-adjusted at the binding pocket; thus, the dynamic equilibrium shifted to a stable population of the complex. The time evolution plot of RMSD shown in Supplementary Figure 2 shows that the complex structure of fistein achieves equilibrium at $40 ns and remains stable till the simulation ends at 100 ns. Whereas, the population of complex structure with kamferol is largely stabilized around $17-20 Å. The splitting in the RMSD distribution plot, which can be seen near 15 Å, 20 Å, and 30 Å, respectively, suggested different drifts in the trajectory during the dynamic progression. The detailed view of the time evolution plot suggested that the structure of hACE2-S protein complexed with kamferol attains equilibrium at $10 ns, which is continued up to 40 ns. A consistent increase in RMSD is observed up to 70 ns; after that, the RMSD dropped quickly and attainded stable equilibrium, which is maintained until 100 ns.

The RMSD distribution plot of complex structure with quercetin shows the higher population density >60% near 23-24 Å, which indicates the stable conformational dynamics of protein-ligand complex. The time evolution view of hACE2-S protein complexed with quercetin also provides a clear view that the complex structure attains equilibrium at $10 ns and remains stable for the remaining period of simulation. Thus, the RMSD results denote that all three ligands remain occupied at active site stabilized with potential interactions.

To evaluate the dynamic adaptability and overall compactness of the ligand-bound complex structure of hACE2-S protein, we calculated the probability distribution plot of radius of gyration (Rg) (Rocco et al., 2008; Singh et al., 2019) . Results show that the Rg values of fistein ranged between 42.0 and 43.0 Å, whereas, the higher population density of Rg values ranged between $42.0-44.0 Å and $47.0-48.0 Å for kamferol and quercetin, respectively (Figure 3) . The narrow Rg distribution ranges of highly populated conformational dynamics suggested the minimal variations in the structural compactness of protein. Thus, the overall protein-ligand complexes remain stable during the simulation, which is observed consistent with RMSD analyses.

To monitor the structural stability of the ligand-bound complex of hACE2-S protein, we also computed the probability distribution plot of the solvent-accessible surface area (SASA), shown in Figure 4 . This plot shows that the distribution of highly populated structure of all three complex restricted around $64,000 Å 2 , which provide elegant evidence of stable conformational dynamics of the protein-ligand complex in an aqueous environment. Thus, the results collectively obtained from probability distribution plots of RMSD, Rg, and SASA suggested that the complexation of ligands at the active site enhances the conformational stability hACE2-S protein.

In addition, to evaluate the molecular interaction, we also examine the time evolution plot of hydrogen bond (H-bond) formed between the protein and ligands ( Figure 5 ). H-bond acts as a key interaction in determining the binding affinity, selectivity, and stabilization effect of ligands with protein.

From Figure 5 , we find that the initial trajectory of fistein shows the appearance of five H-bonds with hACE2-S protein, out of which three H-bonds were observed consistently during the simulation. The H-bond plot of kamferol shows the occurrence of three-four H-bonds with protein; however, two H-bond interactions were found more stable, which can be seen retained throughout the simulation. Whereas, the H-bond plot of quercetin shows a maximum of six H-bond interactions with hACE2-S protein, during the initial 0-15 ns. During the progression of simulation, three H-bonds were observed up to 40 ns, and only two remained long-lived, which were maintained during the entire simulation time. 

To probe the conformational dynamics of protein-ligand complexes during the simulation, we applied principal component analysis (PCA). The PCA is an efficient tool which quantitatively defines the collective motion and measures the movement directions (David & Jacobs, 2014; Kumar et al., 2020; Sarma et al., 2020) . The principal components, PC1 and PC2, were obtained from diagonalization of the covariance matrix of the Ca atomic fluctuation and were used as reaction coordinates to calculate the Gibbs free energy landscape (FEL). FEL plot ( Figure 6 ) displayed that the complex of hACE2-S protein with three compounds, fisetin, kamferol, and quercetin, spans different subspace of structures during the period of 100 ns simulation. It is observed that the collective motion of fistein is confined to a small conformational space. The population of complex structures occupying the small energy basins represents the initial conformation of the optimization phase, the ligand well settled at the binding pocket which is subsequentially visited is broad energy basins, having the population of stable protein-ligand complex. The structure of hACE2-S protein complexed with kamferol navigated a broad conformational space, the collective motion of protein in the three different energy wells represents the different conformational stability of protein-ligand complex. The small energy basin represents the initial equilibration phase; then, the complex structure shifted to another energy basin, which is relatability more stable. Finally, the structure transverse through the small energy barrier to a broader energy basin, having population of a more stable complex. It may be the reason we observed two drifts at $20 ns and $75 ns in RMSD trajectory, representing the two different but stable conformational states between $25-60 ns and $80-100 ns. Quercetin experiences wide conformational space, visiting through the three energy basins. The structure crosses smoothly from one subspace to another subspace through the small energy barrier; thus, we observed that stable conformation of protein-ligand during the initial 0-40 ns, due to excursion protein jump out from the stable energy basin, navigated the broad conformational, representing the rugged energy basin. However, during the dynamics progression, protein acquired equilibrium with the dropdown in RMSD and entered to deep energy basin, suggesting the population of stable protein-ligand structures. Collectively, the structural dynamics analyses and FEL suggested the stable conformation dynamics of the docked complex of hACE2-S protein complex with ligands, fistein, kamferol, and quercetin.

In order to understand the contribution of active site residues which influence the binding free energy and structural stability of ligands during the dynamic progression, we applied MM-PBSA Pandey et al., 2018) . MM-PBSA determines the strength of binding by calculating the absolute binding free energies (DG bind ) which is sum of the electrostatic, Van der Waals interaction, polar and nonpolar solvation free energies (Supplementary Figure 3) . In addition, the contribution of each residue to the binding free energy was also calculated for the aforesaid energy terms. The estimated energy components for binding of fistein, kamferol, and quercetin with hACE2-S protein are summarized in Table 2 . The compound quercetin shows the lowest binding free energy DG bind ¼ À22.17 ± 3.04 kcal mol À1 which indicates the highest binding affinity for hACE2-S protein.

The lowest binding affinity was observed with kamferol which shows DG bind ¼ À15.07 ± 2.42 kcal mol À1 . Whereas, the binding of fistein with hACE2-S protein shows DG bind ¼ À21.11 ± 3.49 kcal mol À1 . This suggests that the binding affinity of quercetin and fistein to hACE2-S protein remains almost with equal potency, which is consistent with molecular docking results, showing similar docked score DG ¼ À8.50 kcal mol À1 . Thus, MM-PBSA results strengthen the molecular docking and conformational dynamics observation, which suggested the stable binding of fistein, kamferol, and quercetin with hACE2-S protein. These results may contribute to the understanding of structure-based designing in the development of COVID-19 therapy.

It is generally used to Gupta et al., 2020; Jayaram et al., 2012; Lipinski, 2004 ) evaluate potential interactions between the drug and other target non-drug molecules. It evaluates the propensity of a compound with a certain pharmacological or biological activity to be used as a potential drug. It acts as a filter to screen potential therapeutic agents/drugs just at the initiation of the program, thereby minimizing the labour and cost of exercises involving clinical drug development and to a large extent preventing late-stage clinical failures. The rule mainly determines the molecular properties of a compound, which are its prime requisition to be a potential drug, like absorption, distribution, metabolism, and excretion (ADME). Lipinski's rule states that for any compound to be selected as a potential drug it should have (a) Molecular mass < 500 Dalton (b) high lipophilicity (expressed as LogP 5) (c) less than 5 hydrogen bond donors (d) Less than 10 hydrogen bond acceptors (e) molar refractivity between 40 and 130. If a compound of interest possesses more than two of the aforementioned criteria, then the compound is likely to be a potential candidate for drug development. All three selected phytochemicals used in this study were found to pass all the five criteria mentioned above in Lipinski's rule (Table 3 ). Thus we suggest that all of these phytochemicals have the potential ability to work effectively as novel drugs.

In vitro and in vivo experiments with our selected group of potential anti-viral compounds will strengthen our perspectives that natural products based therapeutic interventions is the need of the hour. Herbal medicines formulated from the phytochemicals we predicted, will be extremely essential for the identification of anti-COVID-19 inhibitors. 

Molecular docking and MD simulation is the best platform to screen potential drug candidates. In this study, we described the binding ability of naturally occurring phytochemicals with the SARS-CoV2-S protein. Further, using MD simulation, we showed that kamferol, quercetin, and fisetin bind to the hACE2-S-protein complex, near the interface of hACE2 and S-protein binding. Our results were supported by other analyses like the free energy estimation of binding using MM-PBSA and free energy landscape. Further, we found that these compounds satisfy all criteria of Lipinski's rule of five. Therefore, we believe that these compounds could be potential anti-COVID drugs; hence the researcher in this field should verify this using in vitro experimental techniques.

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

Moroccan Medicinal plants as inhibitors against SARS-CoV-2 main protease: Computational investigations

Antiviral Activity of Resveratrol against Human and Animal Viruses

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

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers

Novel inhibitors of severe acute respiratory syndrome coronavirus entry that act by three distinct mechanisms

Targeting endosomal acidification by chloroquine analogs as a promising strategy for the treatment of emerging viral diseases

Structural basis of inhibition of Mycobacterium tuberculosis DprE1 by benzothiazinone inhibitors

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

Canonical sampling through velocity rescaling

Antiviral activity of resveratrol

Potent Inhibition of HIV-1 replication in resting CD4 T cells by resveratrol and pterostilbene

The 2019 novel coronavirus disease (COVID-19) pandemic: A review of the current evidence

Regulation of protein-ligand binding affinity by hydrogen bond pairing

Antiviral activities of extracts and selected pure constituents of Ocimum basilicum

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

Principal component analysis: A method for determining the essential dynamics of proteins

The PyMOL Molecular Graphics System

Effect of resveratrol on herpes simplex virus vaginal infection in the mouse

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

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

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

Antiviral activity of luteolin against Japanese encephalitis virus

Beneficial Role of Phytochemicals on Oxidative Stress and Age-Related Diseases

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

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

CHARMM36m: An improved force field for folded and intrinsically disordered proteins

Principles of early drug discovery

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

Comparative structural and functional analysis of STL and SLL, chitin-binding lectins from Solanum spp

Sanjeevini: A freely accessible web-server for target directed lead molecule discovery

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

Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations

Antiviral effect of flavonoids on human viruses

Targeting SARS-CoV-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like proteinase and 2 0 -O-ribose methyltransferase

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

Pre-fusion structure of a human coronavirus spike protein

Structurebased virtual screening, molecular dynamics simulation and MM-PBSA toward identifying the inhibitors for two-component regulatory system protein NarL of Mycobacterium tuberculosis

Chemistry and biological activities of flavonoids: An overview

Structurebased virtual screening, molecular dynamics simulation and MM-PBSA toward identifying the inhibitors for two component regulatory system protein NarL of Mycobacterium Tuberculosis

More bang for your buck: Improved use of GPU nodes for GROMACS 2018

Conformational states of the severe acute respiratory syndrome coronavirus spike protein ectodomain

Effective inhibition of MERS-CoV infection by resveratrol

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

In silico study of naphtha [1, 2-d] thiazol-2-amine with adenosine A 2A receptor and its role in antagonism of haloperidolinduced motor impairments in mice

Antiherpetic activities of flavonoids against herpes simplex virus type 1 (HSV-1) and type 2 (HSV-2) in vitro

Discovery of novel Methylsulfonyl phenyl derivatives as potent human Cyclooxygenase-2 inhibitors with effective anticonvulsant action: Design, synthesis, in-silico, in-vitro and in-vivo evaluation

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

Novel Carbazole-Piperazine Hybrid Small Molecule Induces Apoptosis by Targeting BCL-2 and Inhibits Tumor Progression in Lung Adenocarcinoma in Vitro and Xenograft Mice Model

Curcumin inhibits Zika and chikungunya virus infection by inhibiting cell binding

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

Flavonoids: A review of probable mechanisms of action and potential applications

Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV

Flavonoids: An overview

Comparison of molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics-three-dimensional reference interaction site model (MM-3D-RISM) method to calculate the binding free energy of protein-ligand complexes: Effect of metal ion and advance statistical test

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

Crystal structure and pair potentials: A Molecular-dynamics study

Anti-infective properties of the golden spice curcumin

Insilico study of the A(2A)R-D (2)R kinetics and interfacial contact surface for face for heteromerization

Identification of lead molecules against potential drug target protein MAPK4 from L. donovani: An in-silico approach using docking, molecular dynamics and binding free energy calculation

BCG vaccination strategy implemented to reduce the impact of COVID-19: Hype or Hope?

Combating COVID-19 with mesenchymal stem cell therapy

Characterization of the protein unfolding processes induced by urea and temperature

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

Delineating the conformational dynamics of intermediate structures on the unfolding pathway of beta-lactoglobulin in aqueous urea and dimethyl sulfoxide

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

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

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

Structure, Function, and Antigenicity of the SARS-CoV

Bioactive flavonoids in medicinal plants: Structure, activity and biological fate

End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design

Inhibition of Epstein-Barr virus reactivation by the flavonoid apigenin

Luteolin inhibits Epstein-Barr virus lytic reactivation by repressing the promoter activities of immediate-early genes

In vitro antiviral activity and projection of optimized dosing design of hydroxychloroquine for the treatment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)

Exploiting correlated molecular-dynamics networks to counteract enzyme activity-stability trade-off

Flavonoids: Promising natural compounds against viral infections

Antiviral properties of resveratrol against pseudorabies virus are associated with the inhibition of IjB kinase activation

Development of novel entry inhibitors targeting emerging viruses

We thank the Indian Institute of Technology Delhi for providing online facilities of Supercomputing facility for Bioinformatics and Computational Biology. AK and SR thanks Mahatma Gandhi Central University Motihari, Bihar.

The authors report no conflicts of interest.

http://orcid.org/0000-0002-4821-1654 Shashikant Ray http://orcid.org/0000-0001-8838-398X