key: cord-0847454-2bo0s1li
authors: Sharma, Aniket; Goyal, Shubham; Yadav, Arvind Kumar; Kumar, Pawan; Gupta, Lovely
title: In-silico screening of plant-derived antivirals against main protease, 3CL(pro) and endoribonuclease, NSP15 proteins of SARS-CoV-2
date: 2020-09-08
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1808077
sha: 238a9d8d98ba9a25134fe0575f978fc25e5cdd75
doc_id: 847454
cord_uid: 2bo0s1li

Novel Coronavirus or SARS-CoV-2 outbreak has developed a pandemic condition all over the world. The virus is highly infectious and spreads by human to human local transmission mode. Till date, there is no vaccination or drugs been approved for the treatment by the World Health Organisation. Henceforth, the discovery of the potential drugs is an urgent and utmost requirement for the medical fraternity. Since, the side effects of plant-derived compounds will be lower compared to synthetic/chemical drugs. The Main protease (3CL(pro) or NSP5) and endoribonuclease (NSP15) proteins are necessity for viral replication and its survival in the host cell. In the present study, in-silico approach of drug development was used to search for potential antiviral plant-derived compounds as inhibitors against SARS-CoV-2 replication proteins. Eight plant-derived compounds of which the antiviral activity was known and available, and two reported drugs against SARS-CoV-2 selected for the molecular docking analysis. The docking results suggested that bisdemethoxycurcumin, demethoxycurcumin, scutellarin, quercetin and myricetin showed least binding energy, i.e., greater than −6.5 Kcal/mol against 3CL(pro) and endoribonuclease of SARS-CoV-2. Further studies of ADME-Tox and bioavailability of drugs were also performed that exhibited efficient parameters of drug likeness. Molecular dynamics simulation calculations were performed for the most negative binding affinity of the compound to evaluate the dynamic behavior,and stability of protein-ligand complex. Our findings suggest that these compounds could be potential inhibitors of SARS‐CoV-2 main protease and endoribonuclease. However, further in-vitro and pre-clinical experiments would validate the potential inhibitors of SARS‐CoV‐2 proteins.

A mysterious breakdown of typical pneumonia by a novel coronavirus , created a pandemic situation all over the world by spreading the disease named as severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2) by the World Health Organization (WHO). This disease is more contagious than previous outbreaks of severe acute respiratory syndrome (SARS) in the year 2002 and the Middle East respiratory syndrome (MERS) in 2012 (Dhama et al., 2020; Guo et al., 2020) . By the end of the year 2019, an outbreak from COVID-19 occurred in Wuhan, Hubei Province of China.

It was reported to infect over 3.3 million people and killed more than 0.23 million people in around 215 countries as of May 03, 2020 (https://www.who.int/emergencies/diseases/ novel-coronavirus-2019). According to WHO situation report-71, India reported 1071 positive cases of SARS-CoV-2 and after one month, WHO situation report-99 reported 29435 positive cases of coronavirus affected patients. India's death rate of COVID-19 is still lower than the global death rate of around 3.4 percent, according to data and statistics website worldometer (https://www.worldometers.info/coronavirus/). SARS-CoV-2 is considered to be a seventh known human coronavirus (hCoV) from the same family after 229E, NL63, OC43, HKU1, MERS-CoV and SARS-CoV that belongs to the b-coronavirus family (coronaviridae). It contains a singlestranded positive-sense (þ) RNA encapsulated within the lipid membrane. It has been reported that the genome of SARS-CoV-2 was 96.2%, 79%, and 50% nucleotide identical with bat-CoV RaTG13, SARS-CoV, and MERS genome, respectively (Dhama et al., 2020; Guo et al., 2020) . Therefore, the bat was suspected to be a natural host of virus origin that infected human via pangolins (Manisjavanica) as an intermediate host.

The genome of SARS-CoV-2 is divided into two parts: open reading frames (ORFs) that encode accessory and structural proteins; and non-structural proteins (NSP). Structural protein involves spike (S) glycoprotein, small envelope (E) protein, matrix (M) protein, and nucleocapsid (N) protein (Boopathi et al., 2020) . Sixteen NSP's are known to encode various enzymes required for replication and transcription process such as main protease (NSP5 and NSP3), RNA dependent-RNA polymerase (NSP12), a helicase/triphosphatase (NSP13), an exoribonuclease (NSP14), an endonuclease (NSP15) . Studies have proved that protein NSP15 of SARS-CoV-2 has 88% and $50% sequence identity, 95% and 65% similarity with SARS-CoV and MERS, respectively . NSP15 protein, an endoribonuclease, one of the crucial NSP for viral replication and transcription. It is highly conserved in all known CoVs (Hackbart et al., 2020) . It is a nidoviral RNA uridylate-specific (NendoU) endoribonuclease containing 346 amino acid and made up of 2 domains NSP-15A1 and NSP-15B. It is an interferon (IFN) antagonist and inhibits interferon-b production via an endoribonuclease activity-independent mechanism. (Liu et al., 2019) There are two proteases encoded in the SARS-CoV-2 RNA strand: papain-like protease (PL P or NSP3) and 3-Chymotrypsin-like protease (3CL pro or NSP5) co-translationally cleave the two polypeptides into mature NSPs (Lim et al., 2000) . The 3CL pro , also known as M pro , has a dominant role in the post-translational processing of the replicase gene. Protease M pro has been reported to be a conserved protein, in favour of new drug target in (þ) strand coronaviruses (Chandel et al., 2020) . It has two homodimer NSP5A and NSP5B. It belongs to the cysteine protease class of protein and responsible for the cleavage of viral peptides into functional units for replication and packaging within the host cells . It is produced by ORF pp1ab only. Targeting M pro enzyme will inhibit the viral maturation and enhance the host innate immune response against COVID-19 (Elmezayen et al., 2020) . Other studies suggested that the sequence similarity of NSP5 of SARS-CoV-2 reported nearly 95% with SARS-CoV . In the current study, two protein targets (3CL pro and endoribonuclease) of SARS-CoV-2 selected based on conserved sequence and their importance in virus replication according to the literature.

For the treatment of COVID-19, doctors are prescribing antimalarial, anti-HIV, anti-influenza drugs and their combinations, but these drugs are not permanent cure of coronavirus infection (Joshi et al., 2020) . Drugs such as Hydroxychloroquine, Chloroquine, Nelfinavir, Rhein, Withaferin A, Withanolide D, Enoxacin, Mercaptopurine and Aloe-emodin are tested against various target proteins of viruses (Muralidharan et al., 2020) . Conversely, various phytochemicals reported in the literature to possess potential antiviral activity that could be an alternative to restrict the coronavirus replication (Elfiky, 2020) . Natural compounds consist of high chemical diversity, lower cost of production than that of biotechnological products or compounds synthesised from combinatorial chemistry and also possess milder or inexistent side effects than chemical drugs (Kitazato et al., 2007) . Therefore, in the present study attempt has been made to screen and evaluate the plant-derived compounds that have therapeutic applications against different genetically and functionally diverse viruses using molecular docking study.

Retrieval of the target structure and their preparation

The protein structure of 3-chymotrypsin-like cysteine protease (3CL pro ) enzyme (PDB ID: 6LU7) with 2.1A˚resolution and endoribonuclease, non-structural protein 15 (NSP15) protein (PDB ID: 6VWW) with 2.2A˚resolution of SARS-CoV-2 were downloaded from RCSB protein databank (http://www.rcsb. org). Initially, different tasks were performed for protein preparation such as charge assignment, solvation parameters and fragmental volumes using SPDVB-4.10 version (Morris et al., 2009) . The protein molecules were further optimized using AutoDock4 Tool for the molecular docking (Morris et al., 2009 ).

The secondary structure of proteins was speculated using the Self-Optimized Prediction Method with Alignment (SOPMA) server at https://npsa-prabi.ibcp.fr/cgi-bin/npsa_ automat.pl?page=/NPSA/npsa_sopma.html (Geourjon & Del eage, 1995) . SOPMA predicts four conformational states including Helix, Beta sheets and bridges, Turns and Coils for each antigen.

The natural compounds for the ligand dataset were selected from the previous literature (Table 2) . Eight selected natural compounds and two drugs were obtained from the PubChem compound database (https://www.ncbi.nlm.nih. gov/pmc/articles/PMC4702940/) in InChI (Inch key) format. The two-dimensional (2D) ligand structures were generated using ACD/Chemsketch and save as mol file format. The three-dimensional (3D) structure of ligands were converted using OpenBabel tool (O'Boyle et al., 2011) .

The docking analysis was performed by molecular docking program AutoDock4 Tool (Morris et al., 2009 ) to explore the poses of potential inhibitors within the targets active site.

For protein-ligand interactions, the Lamarckian genetic algorithm was utilized to perform the docking with the preset parameters. The total number of poses was set to 50. Poses were further clustered using all atom RMSD cutoff of 0.3 Å to remove redundancy and on average 20 cluster representatives were kept. All other parameters for docking and scoring were used as default sets. The protein structure was kept rigid in all steps. All the docking poses and interaction analysis has been performed using the UCSF Chimera (Pettersen et al., 2004) and Discovery Studio Visualizer programs (http:// accelrys.com/products/collaborative-science/biovia-discoverystudio/visualization-download.php). The amino acids residues that displayed interactions with the ligand were documented in Table 3 and Figure S1 and S2 (Supplementary material).

The qualitative measurement of absorption, distribution, metabolism, excretion, and toxicity profile of selected compounds were predicted by using SwissADME program (http:// www.swissadme.ch/index.php) (Hoda et al., 2019) . The major ADME-Tox properties parameters include Molecular weight, H-bond acceptor, H-bond donors, Predicted octanol/water partition coefficient (MLogP), Total Polar Surface Area (TPSA), Lipinski (drug likeness), GI absorption, BBB permeant, Bioavailability score.

To analyze the conformational dynamics of the unbound and bound state of 3CL pro and endoribonuclease enzyme, molecular dynamics simulation has been performed using the Amber16 package (Gajula et al., 2016) . After the preparation of the docked ligand molecules using the antechamber (Wang et al., 2006) , tleap module is used to prepare the cubic simulation box of size 10 Å considering the protein at the center, the box was filled with TIP3P water molecules (Jorgensen et al., 1983) for solvation under periodic boundary condition. Some water molecules were replaced by the counterions to neutralize the systems. The Generalized Amber Force Field (GAFF) for ligand (Wang et al., 2004) and FF14SB FF for protein (Maier et al., 2015) was used during the simulation setup. For all electrostatic interactions calculation, Particle Mesh Ewald (PME) method (Essmann et al., 1995) was used and 2 fs time step is used with the SHAKE algorithm (Ryckaert et al., 1977) to constrain the hydrogen bonds.

Prepared simulation box was initially minimized for 500,000 steps (400,000SD þ100,000CG) to relax the system. The minimized system was gradually heated from the 0 to 300 K over the period of 6000ps. A 1 ns simulation run was further used to equilibrate the system under NVT ensemble. The final MD production run of 50 ns was performed with an integration time step of 2.0 fs. Simulation conformations were collected after every 50ps and all trajectories were analyzed for the RMSD, RMSF, SASA and hydrogen bond dynamics using the CPPTRAJ module (Roe & Cheatham, 2013 ) of the Amber software.

To perform the target protein analysis, FASTA sequences and 3 D structure of 3CL pro and endoribonuclease were retrieved from PDB databank (PDB ID: 6LU7 and 6VWW) ( Figure 1 ). The secondary structure prediction indicated that protein 3CL pro consisted 29.08% a-helix (h), 27.12% extended strand (e), 11.44% beta-turn (t), 32.35% random coil (c) elements and NSP15 endoribonuclease comprised of 24.26% a-helix (h), 26.95% extended strand (e), 9.70% beta-turn (t), 39.08% random coil (c). There are two graphs shown in the result-1) to visualize the prediction; 2) score curves for all predicted states ( Figure 2 ). It also shows the parameters such as window width, number of states etc. that are used for the prediction.

Molecular docking was performed to evaluate the interaction of plant-derived compounds, having reported pharmacological activities against two non-structural proteins of SARS-CoV-2. Non-structural proteins, 3CL pro and endoribonuclease were targeted and the docking pocket was formed at the catalytic sites of the respective proteins as mentioned in Table 1 . Also, the binding affinity of interacted compounds at the active site were assessed. Known anti-viral drugs such as remdesivir and ribavirin were docked against 3CL pro as well as endoribonuclease proteins. Plant-derived active compounds, reported with their pharmacological properties were listed in Table 2 .

Interaction of plant-derived compounds with active site of protein 3CL pro Bisdemethoxycurcumin and Demethoxycurcumin the analogues of curcumin, exhibited the docked score of À7.3Kcal/ mol and À7.02Kcal/mol, and formed five hydrogen bonds with amino acid residues Lys-290, Gly-248, His-235, Thr-341, and Glu-340 and one hydrogen bond with His-163 amino acid residue, respectively ( Figure 3 ). Scutellarin represented the binding affinity of À7.13 Kcal/mol was found to form two hydrogen bonds with amino acid residues Gly-143 and His-163, followed by flavonoid quercetin with binding energy of À6.58Kcal/mol and hydrogen bonds formed with amino acids Asn-142, Gly-143, Ser-144, Leu-141, Cys-145, His-41, and Thr-26. Among all, four plant-derived compounds with binding affinity more than À6.5Kcal/mol were myricetin (-6.15Kcal/mol), curcumin (-6.04Kcal/mol), bergapten (-5.98Kcal/mol) and isoflavone (-5.69Kcal/mol) with amino acids present in the docking pocket of 3CL pro shown in Table 3 . Furthermore, the anti-viral drugs remedesivir and ribavirin showed þ196 Kcal/mol and À5.43 Kcal/mol binding affinity, respectively.

Compounds represented a similar kind of affinity with endoribonuclease as observed with 3CL pro . Based on binding affinities, three compounds were in common with both the protein targeted sites. Demethoxycurcumin with binding affinity of À7.51Kcal/mol and formed five hydrogen bonds with Lys-290, Gly-248, His-235, Thr-341, and Glu-340 amino acid residue. Another compound scutellarin, with binding affinity of À6.97Kcal/mol with Lys-290, His-250, Gly-248, and Glu-340 amino acids formed four hydrogen bonds within docking pocket followed by bisdemethoxycurcumin, binding affinity of À6.56Kcal/mol and one hydrogen bond formation with Lys-290. The flavonoid group of myricetin and quercetin exhibited negligible difference in their binding energies À6.52 and À6.49 Kcal/mol, but hydrogen bond formed with same amino acids residues in the pocket-Lys-290, His-235, His-250, and Glu-340. Bergapten (-5.89Kcal/mol) and isoflavone (-5.69Kcal/mol) also showed the minimum binding energy and formed hydrogen bonds. The anti-viral drugs remedesivir docked with À7.7Kcal/mol with one hydrogen bond, whereas ribavirin binding energy (-5.84Kcal/mol) was found close to that of 3CL pro target site (Table 3; Figure 4 ).

All the plant-derived compounds studied were analysed for their physicochemical properties and ADME to screen bioavailability via oral route. The eight compounds listed in Table 4 showed the five key criteria to qualify a compound to possess drug-like properties. The results explained that the plant-derived compounds studied would be more potent to be used as bioactive compounds as compared to the existing anti-viral drug remedesivir. All the compounds have molecular weight less than 500 g/mol, lipophilicity (MlogP) value less than 5, hydrogen bond donor and acceptor are less than 5 and 10, respectively which corresponds to no violations of Lipinski's rule of five. Thus, it suggested that all of these studied compounds have the potential ability to work effectively as novel inhibitors.

The potential inhibitors of enzymes involved in replication and transcription of SARS-CoV-2 studied were analysed for their oral bioavailability as well as GI tract absorption parameters. The potent compounds extracted from various natural sources were listed in the Table S1 (Supplementary material). Natural sources such as herbs, oils have many active compounds: Curcumin, Quercetin, Bergapten etc. which possess to have potential as anti-SARS-CoV-2 activity. The bioavailability of all the plant-based compounds is equal to 0.55 whereas remdesivir has lower bioavailability of 0.17. The absorption of the compounds by GI tract is high except for myricetin while the drugs have lower absorption rate. 

Molecular dynamics simulation has been preformed for 50 ns simulation time for four systems, two from 3CL pro (receptor only and with selected dock ligand complex) and two from endoribonuclease (receptor and docked ligand complex). For each case, receptor only dynamics is used as a reference to assess the stability dynamics of the ligand bound complex. RMSD ( Figure 5 ), RMSF ( Figure 6 ) and SASA (Figure 7) is calculated to assess the time dependent evolution of the systems and stability dynamics. At the same time, Hydrogen-bond analysis has been done to understand the engagement of the docked ligand with the corresponding receptor (Figure 8 ).

The rapidly spreading coronavirus disease SARS-CoV-2, also known as COVID-19, is declared as a pandemic outbreak by the WHO. It has affected the public as well as social health of more than 200 countries around the globe (https://www. who.int/emergencies/diseases/novel-coronavirus-2019). Hence, it's the need of an hour that has demanded the scientific community from all over the world to collaborate and make an extraordinary effort for the rapid detection of an effective anti-COVID-19 drug.

The NSPs, endoribonuclease (NSP15) and main protease enzyme encoded by 3CL pro have been sequenced and their X-ray crystallography 3-D structures are available in PDB databank (Jin et al., 2020; Kim et al., 2020) . Each monomeric unit of endoribonuclease protein is comprised of $345 amino acids which are folded into three domains: N-terminal, middle and NendoU C-terminal catalytic domain. The NendoU C-terminal catalytic domain reported to possess two b-sheets in antiparallel direction (Sinha et al., 2020) . This protein is required by the virus as it takes part in various viral mechanisms including the formation of the replicase and transcriptase complex (RTC) (Deng et al., 2017) . Also, in innate immune response, therefore NSP15 considered to be essential for biological progression of coronavirus. Study has reported that the main active sites of Nsp15 protein, conserved among SARS-CoV-2, SARS-CoV and MERS-CoV are located in a shallow groove between the two b-sheets that comprise six amino acids (His235, His250, Lys290, Thr341, Tyr343, and Ser294) (Sinha et al., 2020) . The main protease, 3CL pro protein is a homodimer with Cys-His catalytic dyad on the active site which is responsible for the protease activity and hence inhibiting the catalytic site present in enzyme would lead to the inhibition of activity of the protein (Islam et al., 2020; Lobo-Galo et al., 2020) . The secondary structures of 3CL pro , has three domins: Domains I and II which have an antiparallel b-barrel structure and Domain III contains five a-helices arranged into a largely antiparallel globular cluster, and is connected to domain II by means of a long loop region. These features are similar to those of other M pro reported previously (Ren et al., 2013; Wang et al., 2016; Xue et al., 2008) .

At present, there are no approved medications from different drug administrations around the world for the treatment of novel coronavirus (SARS-CoV-2). Hence, the primary focus has been on clinical management which includes the prevention of infection, control measures and intensive supportive care. Recently, anti-malarial drugs like chloroquine, hydroxychloroquine have given positive results in the cell culture systems but has to proceed for the clinical trials to be established as an effective drug (Enmozhi et al., 2020) . Several other anti-viral drugs such as remdesivir, ribavirin used in hepatitis C infection and Ebola, and many other vaccines are under rapid study to be used against coronavirus infection.

An anti-viral activity of curcumin and its analogues demethoxycurcumin and bisdemethoxycurcumin was reported against several different viruses including hepatitis viruses, influenza viruses and emerging arboviruses like the Zika virus (ZIKV) (Mounce et al., 2017) or chikungunya virus (CHIKV) (von Rhein et al., 2016) . These analogues have described role as a metabolite, an anti-neoplastic agent and an anti-inflammatory agent in PubChem databank. They possess much greater chemical stability at physiological pH. Scutellarin is a glycosyloxyflavone, a glucosiduronic acid, a trihydroxyflavone and a monosaccharide derivative. (Cheng et al., 2013) observed that it could inactivate porcine reproductive and respiratory syndrome virus (PRRSV) and restrict the early stage of PRRSV replication under in vitro conditions. Similarly, Zhang et al.(2005) Figure 5. RMSD plot of the 3CL pro system (above) and endoribonuclease system (below). For both cases, protein-ligand complex is showing better stability compared to the receptor alone. Figure 6 . RMSF analysis plot for 3CL pro system (above) and endoribonuclease system (below). As expected, except the binding site residues, nearly similar fluctuations are recorded for the both unbound and bound system. Binding site residues are less fluctuating in ligand complex compared to receptor only dynamics. et al., 2015) . Quercetin is a polyphenol, known as flavonoids, its various pharmacological activities described in the literature such as anti-inflammatory, anti-oxidant, anti-enzymatic activities. (Polansky & Lori, 2020) concluded that quercetin and its derivatives were reported to inhibit the SARS and MERS virus replication via protease inhibition (Cheng et al., 2013) . Myricetin is structurally related to several well-known phenolic compounds namely quercetin, morin, kaempferol and fisetin. Myricetin demonstrated antiviral activity as a strong inhibitor of reverse transcriptase against Rauscher murine leukemia virus (RLV) and human immunodeficiency virus (HIV). Some studies reported that bergapten exhibited potent anti-HIV activity via inhibition of HIV-1 replication with therapeutic index value greater than 5.0 (Shikishima et al., 2001) . According to literature, it has been reported that 5,7,3 0 ,4 0 -tetrahydroxy-2 0 -(3,3-dimethylallyl) isoflavone used as an anti-leishmanial agent (Salem & Werbovetz, 2006) , and it is also described in traditional Chinese medicine records (Zhou et al., 2011) . The preliminary screening of potential molecules having anti-viral activity via computational approach would help in providing the fast selection of promising compounds towards the development of therapeutics against SARS-CoV-2. To elucidate the binding affinity of plant-derived compounds with known properties of anti-viral, were docked against two NSPs of COVID-19 using AutoDock4 software (Table 2; Figure 1 and 2). The binding scores from the docking study revealed the affinity of a specific ligand and strength by which a compound interacts with and binds to the pocket of a target protein. Docking results were successful with least binding scores and presence of significant hydrogen bonding with demethoxycurcumin, bisdemethoxycurcumin and a flavone scutellarin. Compounds having binding affinity of À6.5 Kcal/mol or less are considered to be a good inhibitor of enzymatic activities (Shah et al., 2020) . In the present study, five plant-derived compounds demonstrated binding affinity more negative to À6.5Kcal/mol with 3CL pro , whereas with NSP15 endoribonuclease, five different plant-derived compounds exhibited the same results. These compounds could be potent drugs for inhibition of two essential viral proteins and thus inhibiting the viral replication in the infected cells. Additionally, the in-silico experiment represented the inability of remdesivir to bind with 3CL pro , as binding affinity of remdesivir was þ196 Kcal/mol, which showed instability. In contrast, Remdesivir contributes À6.58 kcal/mol to the binding free energy at the electrostatic S1 subsite of protease protein .

Further, the physicochemical properties of these compounds were studied by ADME-Tox following the Lipinski's rule of five for evaluating drug likeness. Lipinski's rule determines the molecular properties which are important for a drug's pharmacokinetics in the human body such as absorption, distribution, metabolism, and excretion (ADME). Lipinski's rule of five criteria for an ideal drug are (i) a Figure 7 . SASA plot for the 3CL pro system (above) and endoribonuclease system (below). For both systems, SASA are nearly consistent over time and reflects that after ligand binding receptor proteins do not undergo any large conformational change. molecular mass less than 500 Daltons (ii) no more than 5 hydrogen bond donors, (iii) no more than 10 hydrogen bond acceptors, (iv) an octanol-water partition coefficient log P not greater than 5. Three or more than 3 violations of Lipinski's rule do not fit into the criteria of drug likeliness via oral route administration (Benet et al., 2016) . ADME-Tox studies of docked eight plant-derived compounds showed that all the compounds had virtual hits and were successful at passing through these ADME-Tox test filters, stating that they could be potential drugs which can be easily administered orally (Table 4) .

The other crucial experiments like bioavailability and GI absorptions are of detecting the low bioavailability is a major concern for drug design, as it may lead to failure of a novel drug in the clinical trials, in spite of the fact that the compounds have high efficacy in in-vitro and/or in-vivo tests (Gyebi et al., 2020) . Therefore, oral bioavailability of a compound is essential to be taken into account when predicting the compound as a drug candidate. The compounds in the present study exhibited bioavailability of 0.55 with low as well as high GI absorption. The absorption into GI was found to be at a high rate in all compounds except that of Myricetin.

Protein flexibility plays an essential role in the inhibitor designing process (Teague, 2003) and as molecular docking has been performed using the single structure, molecular dynamic simulation can provide an in-depth understanding of the protein-ligand interaction information. 3CL pro and endoribonuclease proteins, along with selected compounds, were used for MDS analysis, which suggests that the ligandbound receptor adopt stability over the receptor only dynamics (Figure 5 and 7) . H-bond based interaction analysis suggests that both docked compounds with respective proteins utilize the available binding interactions and maintain these interactions throughout the simulation time (Figure 8) . Apart from the receptor-based interaction, docked ligand also formed on average two H-bonds with the nearly water molecules which additionally provide the stability to the docked ligand (Supplementary material Figure S3 ).

In the current scenario, there is an urgent requirement of the potential treatment of SARS-CoV-2 infection. Compounds extracted from natural sources might be a source in the development of new potent drugs which can improve the novel coronavirus treatment. Various plant-derived compounds have been reported for their anti-viral activity, may found to be effective against SARS-CoV-2. The bioinformatics Figure 8 . Hydrogen bond dynamics between Bisdemethoxycurcumin-3CL pro complex (above) and Demethoxycurcumin-endoribonuclease complex (below). This analysis shows that selected ligands engage with receptors with nearly two hydrogen bonds and do not deviate much from the initial structure. approach could be a very useful and efficient tool to identify potent inhibitors against the specific targets of novel coronavirus. The docking study concluded that plant-derived compounds (Demethoxycurcumin, Bisdemethoxycurcumin, Scutellarin, Quercetin and Myricetin) exhibited potential drug candidates against two enzymes, main protease (3CL pro ) and endoribonuclease (NSP15) proteins of SARS-CoV-2. Also, molecular dynamics and simulations of two best docked protein-ligand structures revealed the dynamics information of their stability in the biological system. Instability of these two enzyme activities would lead to inhibition of viral replication, which in response would restrict the infection in the host cell.

Effect of genistein on replication of bovine herpesvirus type 1

Curcumin and natural derivatives inhibit Ebola viral proteins: An In silico approach

BDDCS, the Rule of 5 and drugability

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

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

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster

In silico identification of potent COVID-19 main protease inhibitors from FDA approved antiviral compounds and active phytochemicals through molecular docking: A drug repurposing approach

In vitro screening for compounds derived from traditional Chinese medicines with antiviral activities against porcine reproductive and respiratory syndrome virus

Coronavirus nonstructural protein 15 mediates evasion of dsRNA sensors and limits apoptosis in macrophages

Coronavirus Disease 2019-COVID-19

Target Identification and Polypharmacology of Nutraceuticals

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

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

A smooth particle mesh Ewald method

Protocol for molecular dynamics simulations of proteins

SOPMA: Significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments

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

Potential inhibitors of coronavirus 3-chymotrypsin-like protease (3CLpro): An in silico screening of Alkaloids and Terpenoids from African medicinal plants

Coronavirus endoribonuclease targets viral polyuridine sequences to evade activating host sensors

Inhibition of aspergillus fumigatus biofilm and cytotoxicity study of natural compound cis-9-hexadecenal

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

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

Comparison of simple potential functions for simulating liquid water

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

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

Crystal structure of Nsp15 endoribonuclease NendoU from SARS-CoV-2

Crystal structure of Nsp15 endoribonuclease NendoU from SARS-CoV-2

Viral infectious disease and natural products with antiviral activity

Curcumin inhibits herpes simplex virus immediate-early gene expression by a mechanism independent of p300/CBP histone acetyltransferase activity

Identification of a novel cleavage activity of the first papain-like proteinase domain encoded by open reading frame 1a of the coronavirus avian infectious bronchitis virus and characterization of the cleavage products

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

Porcine deltacoronavirus nsp15 antagonizes interferon-b production independently of its endoribonuclease activity

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

ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

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

Open Babel: An open chemical toolbox

UCSF Chimera-a visualization system for exploratory research and analysis

Coronavirus disease 2019 (COVID-19): first indication of efficacy of Gene-Eden-VIR/Novirin in SARS-CoV-2 infection

Antiviral effect of compounds derived from Angelica archangelica L. on Herpes simplex virus-1 and Coxsackievirus B3 infections. Food and Chemical Toxicology

The newly emerged SARS-Like coronavirus HCoV-EMC also has an "Achilles' heel": current effective inhibitor targeting a 3C-like protease

PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data

Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes

Isoflavonoids and other compounds from Psorothamnus arborescens with antiprotozoal activities

In silico studies on therapeutic agents for COVID-19: Drug repurposing approach

Chemical constituents of Prangos tschiniganica; structure elucidation and absolute configuration of coumarin and furanocoumarin derivatives with anti-HIV activity

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

Antiviral activities of natural compounds derived from traditional chinese medicines against porcine circovirus type 2 (PCV2)

Implications of protein flexibility for drug discovery

Synergistic action of quercetin and murine alpha/beta interferon in the treatment of Mengo virus infection in mice

Curcumin and Boswellia serrata gum resin extract inhibit chikungunya and vesicular stomatitis virus infections in vitro

Structure of Main Protease from Human Coronavirus NL63: Insights for Wide Spectrum Anti-Coronavirus Drug Design

Automatic atom type and bond type perception in molecular mechanical calculations

Development and testing of a general amber force field

Structures of two coronavirus main proteases: Implications for substrate binding and antiviral drug design

Identification of myricetin and scutellarein as novel chemical inhibitors of the SARS coronavirus helicase, nsP13

Encyclopedia of traditional Chinese medicines -molecular structures, pharmacological activities, natural sources and applications

The authors would like to thank Amity University Uttar Pradesh, India for providing facilities to carry out this work. The authors also would like to thank the researchers who generated and shared the sequencing and protein data of main protease (3CL pro ) and endoribonuclease (NSP15) proteins of SARS-CoV-2 in Protein Data Bank website (https://www.rcsb. org/). The authors would like to acknowledge Ms. Sapna Ratnakar, Ph.D. Scholar, School of Life Sciences, JNU for providing her support to carry out MDS. P.K. would like to thank DBT Project for supporting the research work and also thank Dr. Debabsisa Mohanty for providing the NII computational facility to carry out the simulation analysis.

The authors declare no potential conflicts of interest to disclose.

All authors conceived, conducted and analyzed the results of experiment(s) equally. L.G. performed docking studies, manuscript editing; A.S. and S.G. drafted manuscript, performed ADMET studies; A.K.Y. manuscript editing, result analysis and P.K. performed molecular dynamics and simulations.

http://orcid.org/0000-0002-2610-3142