key: cord-0991810-wgxomqjz
authors: Singh, Jitender; Malik, Deepti; Raina, Ashvinder
title: Computational investigation for identification of potential phytochemicals and antiviral drugs as potential inhibitors for RNA-dependent RNA polymerase of COVID-19
date: 2020-11-17
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1847688
sha: 7ce3987e9adf49a925eddca67863f14bc0852816
doc_id: 991810
cord_uid: wgxomqjz

Since the SARS/MERS epidemic, scientists across the world have been racing to identify the novel-CoVs as it has been predicted that next epidemic can very well be a result from a new mutation of CoV, for which hundred mutations have already been discovered, and the same fear has come true with world facing a raging pandemic due to COVID-19, an infectious disease caused by a newly discovered coronavirus. COVID-19 or Severe acute respiratory syndrome coronavirus2 (SARS-CoV-2), is a single stranded RNA virus. COVID −19 is highly contagious and has resulted in current global pandemic with almost no country of the world unaffected by this virus. Owing to the lack of effective therapeutics or vaccines, the best measures to control human coronaviruses remain a strong public health surveillance system coupled with rapid diagnostic testing and quarantine/social; distancing/lockdowns as and when necessary. In the present study, we have used the insilico approach for the prediction of novel drug molecules from available antiviral drugs and also from natural compounds that can be best target against RNA-dependent RNA-polymerase (Pol/RdRp) protein of SARS-CoV-2 which can be suitable drugs for the treatment of COVID-19 virus. From the current study we observed that three antiviral and three phyto-chemicals have minimum binding energy with the target protein which were further evaluated in molecular dynamics studies and could specifically bind to RdRp protein of COVID-19. Based on results we suggest that these drugs may act as best or novel inhibitor that may be used for the treatment of SARS-CoV-2. Communicated by Ramaswamy H. Sarma

An RNA-dependent RNA polymerase (RdRp) is the central catalytic subunit of the RNA-synthesizing machinery of all positive-strand RNA viruses. It is presumed that replication and transcription of the $30-kb severe acute respiratory syndrome (SARS) coronavirus (SARS-CoV) RNA genome are catalyzed by an RdRp domain in the C-terminal part of nonstructural protein 12 (nsp12), one of 16 replicase subunits which is responsible for the production of negative strand RNA (-RNA), new genome molecules and in many virus groups also subgenomic (sg) messenger RNAs (mRNAs) (Ahlquist et al., 2003; Miller & Koev, 2000) . Coronaviruses (CoVs) are a family of þ strand single stranded RNA (þssRNA) viruses that have large RNA genome and cause gamut of diseases in mammals and births to respiratory infections in humans which can be potentially lethal (Fehr & Perlman, 2015) . SARS-CoV-2 is genetically similar to bat coronaviruses, giving strong credence to the theory that bats were the primary sources for this virus, also virus found in bats have the same receptor angiotensin converting enzyme 2 (ACE2) as in humans (Benvenuto et al., 2020; Perlman, 2020; Zhou et al., 2020 ). An intermediate animal reservoir such as a pangolin is also suspected to have been the primary source of transmission of this virus to humans. SARS-CoV-2 virions are typically 50-200 nanometers in diameter . All coronaviruses including SARS-CoV-2 has four structural proteins, known as the S (spike), E (envelope), M (membrane), and N (nucleocapsid) proteins; with N protein holding the RNA genome, and the S, E, and M proteins together create the viral envelope (Hasan et al., 2020) . The spike protein with millennia of evolution is perfectly shaped to attach to the membrane of a host cell in stealth mode (Chatterjee, 2020; Lai & Cavanagh, 1997; Woo et al., 2010) . Once the CoV has attached to the host cell it gets a clear path to infiltrate the host cells even before the immune system has chance to know that an external virus has attacked the host organism. Protein modeling experiments on the spike protein of the virus soon suggested that SARS-CoV-2 has sufficient affinity to the angiotensin converting enzyme 2 (ACE2) receptors of human cells to use them as a mechanism of cell entry (Batlle et al., 2020; Chen & Hao, 2020; Walls et al., 2020) . After attaching the host cell SARS-CoV-2, the virion then releases RNA into the cell, forcing the cell to produce copies of the virus that further infect more CoV which in turn infect even more host cells.This virus is assumed to have originated in infamous wet markets of China from where it spread to humans (Andersen et al., 2020; Boopathi et al., 2020; Cyranoski, 2020) . SARS-CoV primarily infects epithelial cells within the lung. The virus is capable of entering macrophages and dendritic cells but only leads to an abortive infection (Chu et al., 2020; Conti et al., 2020) .

In the current COVID-19 outbreak and the previous SARS outbreak it is critical to identify an etiological agent which helps in diagnosing cases in locations where a severe CoV outbreak is occurring. The identification of cases will guide the development, of public health measures to control outbreaks. To date we have not been able to develop anti-viral therapeutics that specifically target human coronaviruses, so treatments are only supportive. The SARS and MERS outbreaks in the past and COVID-19 have stimulated research on these viruses and the ongoing research have identified a large number of suitable anti-viral targets, such as viral proteases, polymerases, and entry proteins. Significant work remains, however, to develop drugs that target these processes and have the ability to inhibit viral replication. In the present study, we have used insilico approach for the prediction of novel drug molecules from FDA approved antiviral drugs including Adefovir, Amantadine, Oseltamivir, Ribavirin and Zanamivir and seven phytochemical compounds from natural plants Curcumin, Demethoxycurcumin, Flavonoid, Isoflavonee, Terpinen-4-ol, Nimbin and Piperine that can be best targeted against RNA-dependent RNA polymerase (Pol/ RdRp) protein of severe acute respiratory syndrome coronavirus 2 which is responsible for replication and transcription of the viral RNA genome. Our study selected five antiviral drugs which have been already approved by FDA hence providing a faster access to patients suffering from COVID-19 if there is any positive correlation is found any of drug and its impact in either clearing the viral load or reducing the severity of the COVID-19 infection. Also as these drugs already have the FDA approval hence there is already sufficient evidence of safety of administering these drugs to humans. The seven natural compounds from different plants have been reported to have many properties like curcumin found to possess antiviral, antifungal, antibacterial and antimicrobial activities (Moghadamtousi et al., 2014) . Demethoxycurcumin has antiviral, antitumor, antifungal, antibacterial, antidiabetic and antimicrobial properties (Agrawal & Goel, 2016) . Flavonoid possesses antiviral effect on human viruses and protective effect on viral Infections (Kaul et al., 1985; Zakaryan et al., 2017) . Isoflavonee is shown to inhibit infectivity of enveloped or non-enveloped viruses, as well as single-stranded or doublestranded RNA or DNA viruses (Andres et al., 2009) . Terpinen-4ol is reported to have an inhibitory effect on influenza virus replication (Garozzo et al., 2011) . Nimbin is a known anti-viral with beneficial effect in viral infections (Badam et al., 1999) . Piperine have many health benefits and disease preventing properties like e.g. anti-inflammatory, antiviral, antipyretic, immune and bioavailability enhancing qualities (Mair et al., 2016) . Based on the above properties we recommend in vitro testing and clinical trials of these highly beneficial compounds for their impact on recovery of patients suffering from Covid-19, potential reduction in severity and duration of the infections. These compounds can potentially help in slowing down the replication of CoV in host providing effective therapeutic strategy against this dreaded virus.

The protein sequences of 4 RNA-dependent RNA polymerase (RdRp) from different variants of severe acute respiratory syndrome coronavirus (SARS-CoV) including SARS-CoV-2 (NCBI ID: YP_009725307), SARS-CoV (NCBI Reference Sequence: NP-828869.1), MERS-CoV (NCBI Reference Sequence: YP_009047223.1) and Bat-CoV (NCBI reference sequence: AID16545.1) were retrieved in FASTA format from National Centre for Biotechnology Information (http://www.ncbi.nlm. nih.gov/) database. Jalview v2 (Thompson et al., 1994) was used for the identification of the conserved region among the sequences through multiple-sequence alignment (MSA) with Clustal Omega (Hall, 1999) . BLOSUM62 distance matrix method through Jalview software was used to construct phylograms from the MSA, and to scrutinise divergence among different strains of the CoV (Waterhouse et al., 2009).

The selected target protein RdRp structure was modelled through homology methods. The template selection and validation of the target protein (RdRp) was done using BLASTP against PDB database to find the best template (https://blast. ncbi.nlm.nih.gov/Blast.cgi?PAGE=Proteins). PDB ID: 6NUR, protein name Chain A, NSP12 (Severe acute respiratory syndrome-related coronavirus) was found as a best template. After the selection of template, a 3 D protein model of target protein was generated using SWISS MODEL method (https://swissmodel.expasy.org/.) (Waterhouse et al., 2018) . Prediction of secondary structure of RdRp model was predicted through PSIPRED 4.0 (http://bioinf.cs. ucl.ac.uk/psipred/) (Jones, 1999) . This algorithm assigns secondary structural elements (a-helices, b sheets and turns).

The accuracy of predicted 3 D structure of target protein and its stereo-chemical properties were analysed using SWISS-MODEL online tool (https://swissmodel.expasy.org/) (Waterhouse et al., 2018) . SWISS MODEL structure assessment, MolProbity reports (http://molprobity.biochem.duke.edu/), QMEAN (Quantitative Model Energy Analysis), Z-score, Clash score, Ramachandran plot, Ramachandran outliers, Ramachandran favoured, Rotamer outliers, C-beta deviations, bad bonds and bad angles were used to determine stereo-chemical properties of target protein. The model structure was optimized by energy minimization through YASARA server (http://www.yasara.org/minimizationserver.htm) (Krieger et al., 2009) , which performs an energy minimization using the YASARA force field. Superimposition of RdRp model and template was done using FATCAT online web tool (http:// fatcat.godziklab.org/) (Yuzhen & Godzik, 2004) which calculates Root Mean Square Deviation (RMSD) values, p-value, differential distance matrix decomposition and Eigen value distribution.

Protein-ligand binding sites are the active sites on protein surface that perform protein functions. Thus, the identification of these binding sites is often the first step to study protein functions and structure-based drug design. Active sites (ligand binding sites) of the receptor protein were analyzed by Prank Web online tool (http://prankweb.cz/). Binding sites are the distribution of surrounding residues in the active sites and act as the catalytic residues (Jendele et al., 2019; Lukas et al., 2019 (Lam et al., 2012) . The nucleotides were used as control in the docking process.

FAFDrugs4 (http://fafdrugs4.mti.univ-paris-diderot.fr/) (David et al., 2008) and AdmetSAR (http://lmmd.ecust.edu.cn/admetsar) (Cheng et al., 2012) server was used to analyze the absorption, distribution, metabolism, and excretion properties of all ligand compounds. ADMET properties of all the ligand molecules were checked in optimal descriptors (hydrogen bonds, charge) at pH ¼ 7.4. The oral toxicity properties and LD50 was also analysed by using PROTOX web server (http://tox.charite.de/protox_II/) (Drwal et al., 2014) .

The Desmond module was used to carry out MD simulations for this study. Firstly, water model were developed and also add sodium and chloride ions to neutralize the drug and target complex by using 'system builder' tool of Desmond, Schr€ odinger. Liquid simulations optimization was minimized at OPLS3 force field. After creating frames for 100ps interval, the complexes were submitted for 50 ns simulations. We performed 50 ns molecular dynamics studies of all the ligand-target complexes. We performed molecular dynamic simulation experiment at 50 ns of all the selected 6 (3 were FDA approved and 3 were phyto-drug molecules) ligands which showed good docking score in virtual screening.

A total of 4 RNA-dependent RNA polymerase (RdRp) protein sequences from different variants of severe acute respiratory syndrome coronavirus (SARS-CoV) were retrieved in FASTA format from National Centre for Biotechnology Information database. The multiple sequence alignment (MSA) of proteins SARS-CoV 2 (NCBI ID: YP_009725307), SARS-CoV (NCBI Reference Sequence: NP_828869.1), MERS-CoV (NCBI Reference Sequence: YP_009047223.1) and Bat-CoV RdRp sequence (NCBI reference sequence: AID16545.1) was performed using Clustal omega software (https://www.ebi.ac. uk/Tools/msa/clustalo/). In the multiple sequence alignment results the black-dashed rectangle marks active site residues (successive aspartate residues D764 and D765). The active site aspartates are projecting from the beta-turn joining the b15 and b16. As implied from the MSA, the active site is highly conserved (highlighted in green colour). Also, the 5 Å . The selected conserved region from the MSA is represented in the Figure  1 . BLOSUM62 distance matrix method through Jalview software was used to construct phylograms from the MSA, with the aim to examine the divergence among the retrieved sequences (Thompson et al., 1994) . Phylogram of RNAdependent RNA polymerase is shown in Figure 2 .

Next with the aim to determine the best ligands against target protein RdRp of COVID-19 (SARS-CoV-2), the selected target protein RdRp of COVID-19 structure was modelled through homology modelling methods using the Swiss Model web server. The template selection and validation of the target protein (RdRp) was done using BLAST-P against PDB database to find the best template.The PDB ID: 6NUR, protein name Chain A, NSP12 of COVID-19 was employed as a best template. On analysis we found that 'template' protein model built in our study was the best template with 100% score (maximum score 1889 from a total score 1889), query cover 99%, E-value was 0.0 and sequence identity was 96.35% suggesting it to be best template for the further study. The MMDB ID of the selected template was 175857 with electron microscope as experimental method and the resolution of selected template was 3.1 Å (Kirchdoerfer & Ward, 2019) . In addition homology modelling was also done by SWISS-MODEL. On analysis from the Swiss Model database we found large number of models with different sequence identity, the model which exhibited a very high sequence identity (96.35%) to the template sequence of PDB ID: 6NUR of NSP12 of COVID-19 was selected. All other models with sequence identity that scored less than 96.35% were eliminated from the study. The predicted 3 D structure of the target protein is represented in Figure 3 . The structure assessment was performed by analysing the secondary structure of RdRp 3 D protein using PSIPRED 4.0 online web tool (Predict Secondary Structure) (Jones, 1999) . The results predicted large number of beta sheets, alpha helix and turns further suggesting that modelled structure were accurate ( Figure 4 ).

After the generation of 3 D structure of target protein, next we tested the quality of the model of the target protein and its stereo-chemical properties using SWISS MODEL online program. We observed that target model protein with clash score 0.24, good Global Model Quality Estimation (GMQE) score 0.83 and Qualitative Model Energy Analysis (QMEAN) score was À0.69 which guarantees the best model quality and also Z score was estimated ( Figure 5A -C). To know the active binding site(s) of 'target' protein structure, we first did the validation of the 3 D structure by Ramachandran plot and we found 97.50% residues were in the favoured region, MolProbity score was 0.69 ( Figure 5D ). In addition we found that there were no observations of any Bad bonds, Cb deviations were 11, and BAD angles observed were 62 out of 8994 angles, and thus validates the quality of the 3 D structure of the target protein (Table 1) .

The model structure was then optimized by energy minimization through YASARA Energy Minimization Server (Krieger et al., 2009) . The structure validation Z-scores and force field energies before and after the minimization were calculated through this tool. The energies of the structure before and after optimization were À351863.4 and À455799.9 kJ/mol respectively. This specifies that the protein structure was stabilized by an amount of À103,936.5 kJ/mol. Moreover, the structure validation Z-scores before and after optimization was À0.61 and À0.04 indicating good optimization of the target protein. The structure of RdRp after energy minimization is presented in Figure 6 . Superimposition of RdRp model and template was done using FATCAT online web tool. p-value is another estimator utilized to assess structural similarity in the FATCAT-flexible method. According to the FATCAT-flexible method, a P-value < 0.05 means that the two structures compared are significantly similar. These two structures shown in Figure 6 are significantly similar with p-value of 0.00e þ 00 (raw FATCAT score was 2346.50). The Root Mean Square Deviation (RMSD) is the measure of the average distance between two aligned proteins (Kufareva & Abagyan, 2012) . In general, smaller RMSD values are accompanied with protein structure pairs with greater similarity (Carugo & Pongor, 2001) . These two structures have 793 equivalent positions with an RMSD of 0.07Ð without twists. Differential distance matrix decomposition distance range from À0.5 to 0.6 Ð and Eigen value distribution range from À26.2 to 19.6 ( Figure 7A -C). These algorithms further confirmed that the target model that we generated is the best quality model.

After the validation of this predicted 3 D structure of target protein, next we were concerned in finding the ligand binding site residues in Pol/RdRp protein using PrankWeb online tool (Jendele et al., 2019) . On analysis of the 3 D structure of RdRp target protein of COVID-19, we found there were top three active sites (meta-pockets) considered as ligand binding sites in this protein. We observed aspartate active site residues in all the Metapockets (Table 2) . Since these three meta-pocket(s) were found to contain aspartate active site residues, they were then studied for further analysis. The aspartate active site residues and most of the 5 Å region surrounding it are surface accessible, though can bind to the free nucleotides (ATP, GTP, and UTP). The active site aspartates are protruding from the beta-turn joining the b15 and b16. The active site is surface accessible as we can see from the surface representation allowing the interaction with the free nucleotides passing through the nucleotide channel of the RdRp . The 3 D structure showing the position of aspartate active sites (D-760 and D-761) is represented in Figure 8 .

After finalizing the meta-pocket sites on the basis of important amino acid residues, a total of 12 compounds were selected, out of which five were FDA approved antiviral drugs: Adefovir, Amantadine, Oseltamivir, Ribavirin and Zanamivir and seven phyto-chemical compounds from natural plants were Curcumin, Demethoxycurcumin, Flavonoid, Isoflavone, Terpinen-4-ol, Nimbin and Piperine. Table 3 summarizes 7 major phyto-compounds with their pharmacological properties. All these compounds were used as ligands for docking with the selected COVID-19 target protein. The 3 D structures of all the 12 compounds were downloaded from the NCBI PubChem Database ( Figure 9A ,L). These 12 compounds were then subjected to Lipinski rule of five to evaluate drug-likeness (Lipinski, 2004) and we found that all the 12 selected compounds followed the Lipinski rule of five suggesting their role as a drug compounds (Table 4 ).

These 12 drug compounds and three control nucleotides (ATP, GTP and UTP) were then checked for their binding with target protein and docking score was evaluated using AutoDock vina PyRx software (Trott & Olson, 2010) . The target receptor file and all ligand files were prepared to AutoDock Vina for the docking with individual drug molecules. A grid map was set grid center X: 141.1104, Y: 142.3415, Z: 147.8521 and number of points X: 93, Y: 70, Z: 70, spacing was 0.3750 Angstrom. For each compound, many docking poses were obtained with the target protein.

The binding interaction(s) of all the docked 3 D and 2 D poses of receptor target protein of COVID-19 with various antiviral drugs is represented in Figure 10A -C. Similarly interaction of phytochemicals is represented in 3 D and 2 D poses in Figure 11A -C. After docking process was complete the results were further analysed for binding energy (DG), interacting residues and number of H-bonds (Table 5 ). The complex structure of antiviral and phyto drugs with target protein as a 3 D cartoon is represented in Figure 10A1 -A5 and Figure 11A1 -A7 respectively. Various hydrogen bond interactions in 3 D form shown by the antiviral drug and phto drug molecules are represented in Figure 10B1 -B5 and Figure 11B1 -B7 respectively. Similarly number of residue interactions and different bond interactions of all antiviral drug molecules and Phyto-drug molecules in 2 D form is shown in Figure 10C1 -C5 and Figure 11C1 to C7 respectively. Out of many docking poses, only those were selected which had the minimum binding energy, good hydrogen bond interaction and which interacted with aspartate active sites. From our results we found that three FDA approved . Superimposition of RdRp target protein and template (A) Representative comparative FATCAT image between RdRp protein and template. Horizontal axis represent target protein PDB and vertical axis indicate template PDB (B) Image displaying differential distance matrix decomposition distance range for both template and target protein with values ranging from À0.5-0.6 Ð (C) Eigen value distribution curve for both template and target protein ranging from À26.2-19.6. Table 2 . Table indicating amino acid positions of active sites within these meta-pockets. Blue colour residues indicate aspartate active site residues which are surface accessible. Metapockets were analysed using Prank Web online tool (http://prankweb.cz/).

Amino and binding with aspartate active residues with the target protein (RdRp target protein of COVID-19) which is shown in Table 5 . As aspartate active site residues are surface accessible residue which allows the interaction with free nucleotides passing through the nucleotide channel of RdRp protein.

We have also docked these three nucleotides (ATP, GTP and UTP) as a control drug which also showed the minimum binding energy. Graphical representation and comparison of all drugs and control molecules are shown in Figure 12 .

After the confirmation of the interaction of these drugs with target protein RdRp of COVID-19, next the aim of the study Kaul et al. (1985) , Zakaryan et al. (2017) Parsley (Petroselinum crispum)

Shown to inhibit the infectivity of enveloped or non-enveloped viruses, as well as single-stranded or double-stranded RNA or DNA viruses Andres et al. (2009) Red clover (Trifoliumpratense L.)

Have an inhibitory effect on influenza virus replication. Garozzo et al. (2011) Tea tree (Melaleucaalternifolia ) Nimbin

The antiviral and virucidal effect of methanolic extract fraction of leaves of neem Badam et al. (1999) Neem ( Demethoxycurcumin L) Flavonoid. SDF structure retrieved from PUBCHEM database was converted to PDB structure using Discovery-studio visualizer free software.

was to scrutinize the pharmacokinetics for absorption, distribution, metabolism, and excretion (ADME) properties of all these drug compounds using FAF Drugs3 (Lagorce et al., 2015) and admetSAR (Cheng et al., 2012) server. ADMET stands for Absorption, Distribution, Metabolism, Excretion and Toxicity. The prediction of the ADMET properties plays an important role in the drug design process because these properties account for the failure of about 60% of all drugs in the clinical phases. We also analysed these drugs for oral toxicity (LD50) properties by using PROTOX web server (Drwal et al., 2014) . We found that Adefovir, Amantadine, Oseltamivir, Ribavirin and Zanamivir have LD-50 values 13 mg/kg, 157 mg/kg, 260 mg/kg, 2700 mg/kg and 5000 mg/ kg respectively and the phyto-drugs including Curcumin, Demethoxycurcumin, Isoflavone, Flavonoid, Terpinen-4-ol, Nimbin and Piperine have LD-50 values 2000 mg/kg, 2000 mg/kg, 500 mg/kg, 5000 mg/kg, 1016 mg/kg, 1000 mg/ kg and 350 mg/kg. ADMET properties and toxicity calculation of drug molecules which interacted with the target protein of COVID-19 are shown in Tables 6 and 7. Our results showed that all the antiviral drug molecules and phyto-compounds we selected have promising and good ADMET properties and also showed high quality QSAR models.

We performed molecular dynamic simulation experiment at 50 ns of all the selected 6 (3 were FDA approved & 3 were phyto-drug molecules) ligand which showed good docking score in virtual screening. In our study of 50 ns MD simulation study we found that the average RMSD and RMSF value of all top drug molecules Adefovir, Ribavirin, Zanamivir, Curcumin, Demethoxycurcumin and Piperine was found to be (RMSD) 3.2 Å, 3.0 Å, 3.1 Å, 2.8 Å, 3.0 Å, 3.2 Å and (RMSF) 0.8 Å, 0.6 Å, 0.8 Å, 0.6 Å, 0.8 Å, 0.8 Å. (Figures 13-14 and Table  8 ). When comparing the RMSF value of the all the drug molecules we could find that higher fluctuations were noted and details of interactions of these ligands with the target protein are showed in Figure 15 .

Molecular docking studies explored the binding of FDA approved antiviral drugs: Adefovir, Amantadine, Oseltamivir, Ribavirin and Zanamivir and seven phyto-chemical compounds from natural plants Curcumin, Demethoxycurcumin, Flavonoid, Isoflavone, Terpinen-4-ol, Nimbin and Piperine and three control molecules (ATP, GTP and UTP) for interaction with target protein (Pol/RdRp) of SARS coronavirus 2 (COVID-19). The molecular docking study showed that selected three out of five FDA approved drug molecules (Zanamivir, Ribavirin and Adefovir) and three out of seven different phyto- (Table 5 ).

In addition, all these shortlisted drugs have also been found to interact with aspartate active site residues. Similar to our study, showed the five approved drugs i.e. Galidesivir, Remdesivir, Tenofovir, Sofosbuvir, and Ribavirin also have the ability to bind the SARS-CoV-2 RdRp, with binding energies of À7.0, À7.6, À6.9, À7.5, and À7.8 kcal/mol, Figure 12 . Graphical representation of binding energy calculated by AutoDock vina for FDA approved antiviral drugs, phyto-drug molecules and nucleotides (ATP, GTP and UTP) from minimum to maximum energy score in kcal/mol. respectively. Huang et al., 2020 also showed many herbal compounds used in traditional Chinese medicine therapy have been screened and found to have free binding affinity against RdRp protein of SARSCoV-2. In another study Pandit & Latha, 2020 reported the binding affinity of phyto-compounds derived from Silybum marianum (Silybin), Withania somnifera (Withaferin A), Tinospora cordifolia (Cordioside) and Aloe barbadensis (Catechin and Quercetin) with SARS-CoV-2 target. The binding energy of phyto-compounds was more than the widely used for hydroxychloroquine and other repurposed drugs used for treatment of COVID-19 infection. Tatar and Turhan (2020) reported the binding affinity of nafamostat, rapamycin, saracatinib, Imatinib and Camostat with binding energy À10.24, À9.28, À9.66, À9.23 and 9.07 respectively with N-Terminal domain of SARS-CoV-2 Nucleocapsid. Gupta et al. (2020) and Muralidharan et al. (2020) have also done similar studies. Ubani et al. (2020) reported best binding affinity of scopodulcic acid and 249 Dammarenolic acid with the Spike glycoprotein (6VSB) and the Mpro 250 (6FLU7) respectively. It is fascinating to find that three FDA approved antiviral drugs including Zanamivir, Ribavirin and Adefovir molecules and three Phyto-drug molecule comprising Piperin, Curcumin and Demethoxycurcumin showed minimum binding energy and hydrogen bond interactions alongwith interaction with aspartate active site with the target protein (RdRp target protein of COVID-19 (Table 5) . Interestingly, we observed that LD50 of Zanamivir and Ribavirin was 5000 mg/kg and 2700 mg/kg respectively signifying that these compounds are less toxic. Likewise Curcumin and Demethoxycurcumin compound are also found to be non-toxic with high LD50 value (LD50 > 1000 indicates that drug molecules are non-toxic while LD50 < 1000 suggests drugs are highly toxic). Furthermore, all the three FDA approved antiviral drugs and three Phyto-drug molecules fulfill all the criteria of Lipinski's rule of five and ADME/T. Similar to our study, Ahmed and Shohael, (2019) reported that the chrysophanol, aloe-emodin, rhein and emodin all fulfil the criteria of Lipinski's rule of five and ADME/T. Collectively the data and our study suggests that these selected compounds i.e. Zanamivir, Ribavirin, Adefovir and Curcumin, demehoxycurcumin, Piperine with minimum binding energy with the target protein may act as potential candidates for the design of novel drugs for the treatment of COVID-19.

In this study we studied twelve drugs like compounds five from available antiviral drugs and seven from natural chemical compounds. The three best FDA approved antiviral drugs including Zanamivir, Ribavirin and Adefovir used for other viral diseases and three phyto-drug like compounds like Curcumin, Demethoxycurcumin and Piperine could specifically bind to RdRp protein of COVID-19 showing minimum binding energy which attests to the good binding affinity of these drugs to the target protein. On the basis of molecular docking energy score of all selected three FDA approved drug molecules and three phyto-drug molecules from different plants suggest that it may act as best or novel inhibitor that may be used for the treatment of severe acute respiratory syndrome coronavirus 2. We need in-vitro proof of concept studies of these identified agents to find out a potential inhibitor of SARS-2. Hence this information may be useful for the designing of the novel drug molecules against COVID-19.

There is no conflict of interest.

No funding received as it is an insilico study.

Curcumin and its protective and therapeutic uses

Host factors in positive-strand RNA virus genome replication

In silico studies of four anthraquinones of Sennaalata l. As potential antifungal compounds

Correspondence: The proximal origin of SARS-CoV-2

Soy isoflavones and virus infections

In vitro' antiviral activity of neem (Azadirachtaindica. A. Juss) leaf extract against group B coxsackieviruses

Soluble angiotensin-converting enzyme 2: A potential approach for coronavirus infection therapy

The 2019-new coronavirus epidemic: Evidence for virus evolution

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

A normalized root-mean-square distance for comparing protein three-dimensional structures

Understanding the nature of variations in structural sequences coding for coronavirus spike, envelope, membrane and nucleocapsid proteins of SARS-CoV-2. Envelope, membrane and nucleocapsid proteins of SARS-CoV-2

The role of angiotensin converting enzyme 2 in coronaviruses/influenza viruses and cardiovascular disease

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

admetSAR: A comprehensive source and free tool for assessment of chemical ADMET properties

Comparative replication and immune activation profiles of SARS-CoV-2 and SARS-CoV in human lungs: An ex vivo study with implications for the pathogenesis of COVID-19

Induction of pro-inflammatory cytokines (IL-1 and IL-6) and lung inflammation by Coronavirus-19 (COVI-19 or SARS-CoV-2): Anti-inflammatory strategies

Mystery deepens over animal source of coronavirus

FAF-Drugs2: Free ADME/tox filtering tool to assist drug discovery and chemical biology projects

ProTox: A web server for the in silico prediction of rodent oral toxicity

Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): A molecular docking study

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

Coronaviruses: An overview of their replication and pathogenesis

Activity of Melaleuca alternifolia (tea tree) oil on Influenza virus A/PR/ 8: study on the mechanism of action

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

BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT

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

A review of therapeutic agents and Chinese herbal medicines against

PrankWeb: A web server for ligand binding site prediction and visualization

Protein secondary structure prediction based on position-specific scoring matrices

Antiviral effect of flavonoids on human viruses

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

Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8

Methods of protein structure comparison

The molecular biology of coronaviruses

FAF-Drugs3: A web server for compound property calculation and chemical library design

Genotype and subtype profiling of PSI-7977 as a nucleotide inhibitor of hepatitis C virus

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

Prank Web: A web server for ligand binding site prediction and visualization

Antiviral and anti-proliferative in vitro activities of piperamides from black pepper

Synthesis of subgenomic RNAs by positive-strand RNA viruses

A Review on antibacterial, antiviral, and antifungal activity of curcumin

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

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

Another decade, another coronavirus

Investigation of N terminal domain of SARS CoV 2 nucleocapsid protein with antiviral compounds based on molecular modeling approach

CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

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

Molecular docking analysis of some phytochemicals on two SARS-COV-2 targets

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

SWISS-MODEL: Homology modelling of protein structures and complexes

Jalview Version 2-a multiple sequence alignment editor and analysis workbench

Coronavirus genomics and bioinformatics analysis. Viruses

FATCAT: A web server for flexible structure comparison and structure similarity searching

Flavonoids: Promising natural compounds against viral infections

A pneumonia outbreak associated with a new coronavirus of probable bat origin

 III  III  III  III  III  III  III