key: cord-0779894-erabumfn authors: Pandey, Khushboo; Lokhande, Kiran Bharat; Swamy, K Venkateswara; Nagar, Shuchi; Dake, Manjusha title: In Silico Exploration of Phytoconstituents From Phyllanthus emblica and Aegle marmelos as Potential Therapeutics Against SARS-CoV-2 RdRp date: 2021-06-24 journal: Bioinform Biol Insights DOI: 10.1177/11779322211027403 sha: 46888b2eb2a327b9aeac509725099afdf6f659b4 doc_id: 779894 cord_uid: erabumfn Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) worldwide has increased the importance of computational tools to design a drug or vaccine in reduced time with minimum risk. Earlier studies have emphasized the important role of RNA-dependent RNA polymerase (RdRp) in SARS-CoV-2 replication as a potential drug target. In our study, comprehensive computational approaches were applied to identify potential compounds targeting RdRp of SARS-CoV-2. To study the binding affinity and stability of the phytocompounds from Phyllanthus emblica and Aegel marmelos within the defined binding site of SARS-CoV-2 RdRp, they were subjected to molecular docking, 100 ns molecular dynamics (MD) simulation followed by post-simulation analysis. Furthermore, to assess the importance of features involved in the strong binding affinity, molecular field-based similarity analysis was performed. Based on comparative molecular docking and simulation studies of the selected phytocompounds with SARS-CoV-2 RdRp revealed that EBDGp possesses a stronger binding affinity (−23.32 kcal/mol) and stability than other phytocompounds and reference compound, Remdesivir (−19.36 kcal/mol). Molecular field-based similarity profiling has supported our study in the validation of the importance of the presence of hydroxyl groups in EBDGp, involved in increasing its binding affinity toward SARS-CoV-2 RdRp. Molecular docking and dynamic simulation results confirmed that EBDGp has better inhibitory potential than Remdesivir and can be an effective novel drug for SARS-CoV-2 RdRp. Furthermore, binding free energy calculations confirmed the higher stability of the SARS-CoV-2 RdRp-EBDGp complex. These results suggest that the EBDGp compound may emerge as a promising drug against SARS-CoV-2 and hence requires further experimental validation. has been referred to as SARS-CoV-2 by the World Health Organization (WHO) in the year 2019. 1 The virus belonging to the order Nidovirales and genus Betacoronavirus 2,3 has imposed a major challenge in front of human society due to global health concerns. The disease has become pandemic worldwide due to the human-to-human transmission cycle of such a pathogenic virus. The disease originated in China in 2019 has now spread across all over countries with around 170 426 245 confirmed cases of COVID-19, including 3 548 628 deaths as reported by WHO (https://covid19.who. int/) on the date June 1, 2021. Several targets currently being focused on for identifying novel inhibitors include SARS-CoV-2 spike protein, Angiotensin-converting enzyme-2 (ACE-2), human proteases including transmembrane protease, Serine 2 (TMPRSS2), furin and viral proteases like RNA-dependent RNA polymerase, papain like proteases-2 (PLpro), and main protease (M pro ). [4] [5] [6] [7] [8] [9] [10] Two types of proteins, ie, structural and non-structural proteins characterizing human coronaviruses (HCoVs). SARS-CoV-2, structural proteins include Spike (S), Nucleocapsid (N), Matrix (M), and Envelope (E), whereas non-structural proteins include Nsp1 up to Nsp16 11,12 along with the RNA-dependent RNA polymerase (RdRp) also known as Nsp12. 13, 14 The RdRp, a vital enzyme in the life cycle of RNA viruses, has been targeted in various viral infections including Hepatitis C virus (HCV), Zika virus (ZIKV), West Nile virus, and Japanese encephalitis virus ( JEV). [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] RdRp, a protein crucial for viral replication is a promising druggable antiviral target for coronavirus. 27 RdRp of SARS-CoV-2 supports the transcription and replication of a large RNA genome with approximately 30 000 nucleotides [28] [29] [30] and is the most highly conserved protein among RNA viruses. 31 RdRp activity is dependent on magnesium ions and requires the non-structural proteins Nsp7 and Nsp8 for complete activity. 32 Nucleotide and nucleoside analogs inhibiting the action of RNA polymerases would be considered promising antiviral agents. 33 Compounds showing broad-spectrum antiviral activities with their proven action against coronaviruses are under focus. Several such compounds are being evaluated through random clinical trials 34, 35 where Remdesivir (RDV) is used5 in the form of a common substrate of several viral RdRp enzymes. [36] [37] [38] [39] Remdesivir was approved by the U.S. Food and Drug Administration (US-FDA) on May 1, 2020 for emergency use to treat COVID-19. 40 Nowadays, phytochemicals are the key sources of antiviral drugs with minimal side effects and therefore are of global interest to identify and explore their potency to treat SARS-CoV-2 viral disease. 41 The extracts derived from medicinal plants have been found to inhibit replication of HSV-2, 42 HIV, 43, 44 HBV, 45, 46 and SARS-CoV-2 virus. 47 The current approach is the screening of bioactive compounds from medicinal plants, Aegle marmelos and Phyllanthus emblica that have been already reported to show antiviral activities against Human coxsackieviruses B1-B6, and herpes simplex virus, respectively. [48] [49] [50] The Aegle marmelos being an important medicinal plant with several bioactive compounds such as Seselin, Aeglein, Marmelide, and Marmelosin is chosen for antiviral study against SARS-CoV-2. 51-55 Some Phyllanthus species also exhibit inhibitory potential against HBV, HCV, HIV, and HSV [56] [57] [58] and Phyllanthus emblica against respiratory syndrome virus (PRRSV). 59 Chebulagic Acid, Pedunculagin, Azadirachtin, and Nimbolide may act as potential inhibitors of the SARS-CoV-2 ACE2 receptor and M pro . 60 Seselin showed efficacy against multiple targets of SARS-COV-2 such as spike protein and M pro . 61 However, evidence of these compounds for the treatment of COVID-19 is lacking. Therefore, these compounds were studied for their inhibitory action against SARS-CoV-2 RdRp. In the present study, in silico approach was employed to identify the therapeutic potential of a set of phytocompounds from Aegle marmelos and Phyllanthus emblica against SARS-CoV-2 RdRp. The potential of these phytocompounds (Table 1) was accessed by performing molecular docking and dynamic simulation studies for analyzing their binding affinity with SARS-CoV-2 RdRp. The molecular field-based similarity of Remdesivir with the obtained lead compound has been employed. The efficacy of these compounds has been compared with Remdesivir for experimental COVID-19 therapy. The SARS-CoV-2 RdRp structure complexed with an RNA template, and Nsp7, and Nsp8 cofactors crystallized using electron microscopy at 2.93 Å resolution has been selected as a target for the current study. This complex was downloaded from Protein Data Bank with PDB ID: 7C2K. 62 The receptor structure was prepared using the Protein Preparation Wizard tool of the Maestro program (Schrodinger Release 2020-3: Maestro) at pH 7.5. Removal of water molecules and zinc ion bound with the SARS-CoV-2 RdRp complex along with the subsequent addition of all missing hydrogen atoms and disulfide bonds was carried out. To prevent charge repulsion from the free termini and ensure conformational stability of the protein, N-termini and C-termini of the SARS-CoV-2 RdRp were capped with N-acetyl and N-methyl amide groups respectively. Finally, the prepared and minimized structure of SARS-CoV-2 RdRp was subjected to further docking and molecular dynamic simulation studies. To identify promising SARS-CoV-2 RdRp inhibitors, a set of 5 bioactive phytocompounds reported for their antiviral activity, and Remdesivir, a standard reference molecule was retrieved from the PubChem database (https://pubchem.ncbi. nlm.nih.gov/). 63 To get the energetically stable conformation of the retrieved compounds, the energy minimization was performed using the Maestro program with the OPLS3e force field. The most stable conformation having low energy for each compound was used for the docking calculation. The phytocompounds and Remdesivir with their PubChem IDs are enlisted in Table 1 . From the literature study, amino acid residues like Phe441, Asp452, Tyr455, Tyr456, Lys545, Ala547, Arg553, Arg555, Thr556, Ala558, Tyr619, Cys622, Asp623, Arg624, Thr680, Thr687, and Asp760 11,64 are important for ligand binding within the SARS-CoV2 RdRp cavity were selected as binding site residues for the docking study ( Figure 1 ). To study the mode of interactions of phytocompounds in the defined binding site of SARS-CoV-2 RdRp, molecular docking was performed using FlexX software (LeadIT 2.3.2). 65, 66 This software uses an incremental buildup algorithm guiding the flexible placement of ligand in the binding region. 67 FlexX considers ligand flexibility by changing the conformations of the ligand in the active site while making the protein rigid. 68 It is an extremely fast, robust, and highly configurable computer program for predicting protein-ligand interactions. 69 The SDF file of the compounds was uploaded in FlexX as a docking library. FlexX default docking parameters were kept with 200 conformations per iteration and maintained conformation per fragmentation. The top-ranked poses were selected for the interaction study, the FlexX rank the resulting docked poses as per the FlexX score (docking energy). The intermolecular interactions between the SARS-CoV-2 RdRp and ligands were studied using the Maestro program (Schrodinger Release 2020-3: Maestro). Binding poses of the ligands visualized in Maestro were selected for further studies based on their binding affinity and mode of interactions with SARS-CoV-2 RdRp. Subsequently, the compounds were subjected to 100 ns MD simulation to probe the binding stability. To assess the dynamic behavior and binding stability of the docked compounds in the binding cavity of the SARS-CoV-2 RdRp, 100 ns of MD simulation was carried out using the Desmond 70 module of Schrodinger software (Schrodinger Release 2020-3: Desmond). The docked complexes were solvated using an implicit solvent model in a cubic box of 10 Å spacing. The solvated systems were neutralized with counter ions, and physiological salt concentration was maintained to 0.15 M. 71 All the systems were set up at constant NPT (N[number], P[pressure], T[time]) ensemble 72 with atmospheric pressure (1.013 bar), and a constant temperature of 300k over 100 ns simulation. 73 The systems were minimized and equilibrated with the default settings of relaxation with a time step of 2 fs. The OPLS3e force field was designated to the protein-ligand complex systems. 74 The results were analyzed using MD trajectories generated during 100 ns simulation as root mean square deviation (RMSD), root mean square fluctuation (RMSF), and protein-ligand interactions. Protein-ligand interactions were recorded throughout the 100 ns MD simulation and examined various intermolecular interactions between SARS-CoV-2 RdRp and lead compounds. The protein-ligand interaction profile was normalized throughout the trajectory wherein, the binding site residues were provided that are interacted with SARS-CoV-2 RdRp by various types of interactions. The molecular mechanics generalized born and surface area (MM/GBSA approach) was used to compute the binding-free energies of the complex system using the Prime module. 75 MM/GBSA is a method to calculate binding energy, which uses energy properties of free ligand, free receptor, and receptor-ligand complex for binding affinity calculation. The more negative MM/GBSA score indicates the formation of more stable protein-ligand complexes. 76, 77 We have computed the MM/GBSA from the entire trajectory of the MD simulation run using the trajectory clustering method for each SARS-CoV-2 RdRp complex system. For this purpose, we have taken the entire 100 ns trajectory of MD simulation and extracted the coordinate file at every 10 ns interval. From these, we have calculated the ensemble-averaged binding-free energies (MM/ GBSA) value. The free energy of binding can be calculated as where ∆E elec and ∆E vdW are the electrostatic and van der Waal's contributions and G polar and G non-polar are the polar and non-polar solvation terms, respectively. 78-80 FieldTemplater, a component of Forge-Cresset software, 81 was used to perform a conformational search using the extended electron distribution (XED) force field. The technique uses "field points" as a simple and effective descriptor of the electrostatic and van der Waals maxima and minima surrounding a molecule equipped with XED charges. 82, 83 Compound with the best binding energy and mode of interactions was selected. The processing was performed using default parameters for the generation of a bioactive field template with a single common field pattern reflecting the binding requirements of the selected compound and Remdesivir. The detailed intermolecular interactions analysis is summarized in Figure 2A , the -NH group adjacent to the Phosphate group in the structure of Remdesivir is forming hydrogen bond interaction with the carboxylic oxygen atom of Asp623 amino acid of SARS-CoV-2 RdRp. The carboxylic oxygen atom of Thr556 is hydrogen-bonded with the hydroxyl group of Remdesivir. Hydrogen bonding interactions are also observed with Asp452, Cys622, and Arg624 of SARS-CoV-2 RdRp. The binding energy obtained for EBDGp with SARS-CoV-2 RdRp is −23.32 kcal/mol. The intermolecular interaction of EBDGp with SARS-CoV-2 RdRp is showing hydrogen bond interactions with the same crucial amino acids Thr556 and Asp623 as observed with Remdesivir. The carboxylic oxygen atom of amino acids Thr556 and Asp623 forms hydrogen bond interactions with the hydroxyl group of EBDGp as shown in Figure 2B with the bond distance of 1.99 and 1.60 Å, respectively. Thr556 of SARS-CoV-2 RdRp forms similar molecular interactions with both Remdesivir and EBDGp. It can be deduced from the above results that EBDGp is showing a similar mode of interactions with amino acid residues Thr556 and Asp623 as Remdesivir despite Asp623 is forming interactions with Remdesivir and EBDGp with different functional groups but forming the similar type of intermolecular interactions, ie, Hydrogen bond. In addition, Asp452, Arg555, Tyr619, Cys622, Thr687, and Asp760 are also involved in hydrogen bond interactions with EBDGp. The Marmelide shows the binding energy of −10.30 kcal/mol which is observed to be lower than that of Remdesivir and EBDGp. Interactions of Marmelide with SARS-CoV-2 RdRp are shown in Figure 3A . Hydroxyl group of Marmelide shows to form hydrogen bond interaction with an amino group of Thr556 (2.02 Å) and the secondary amino group of Arg624 with a bond distance of 1.79 Å. Six carbon aromatic ring of Seselin is forming a π-cation bond with the secondary amino group of Arg553 (3.04 Å) and the side-chain amino group of Lys621 (5.85 Å), and backbone amino group of Cys622 (1.98 Å) forms a hydrogen bond with the carbonyl oxygen atom of the Seselin ( Figure 3B ). Pedunculagin and Chebulagic acid show the lowest binding energies (−1.57 and 1.64 kcal/mol, respectively) among all the docked compounds. Interactions of Pedunculagin and Chebulagic acid with SARS-CoV-2 RdRp are shown in Figure 4A and B, respectively. Hydroxyl groups of Pedunculagin forms 2 hydrogen bonds with a carboxylic group of crucial amino acid, ie, Thr556 with a bond distance of 1.57 Å and 2.05 Å. Carbonyl oxygen of Pedunculagin forming a hydrogen bond with an amino group of Cys622 with a bond distance of 1.91 Å. Besides, Arg553 and Lys621 are also forming hydrogen To evaluate the dynamic behavior, 100 ns simulation runs for the docked compounds including Remdesivir in the defined binding pocket of SARS-CoV-2-RdRp was carried out. The information about the structural stability of the protein-ligand complex could be analyzed by RMSD. Root mean square deviation calculations were performed using changes in C-alpha atoms of SARS-CoV-2-RdRp in complex with docked phytocompounds. In each complex, it appears that stable equilibrium was reached after 5 ns. The RMSD was observed within 4.5 Å RMSD for the receptor in complex with Seselin, Marmelide, Pedunculagin, and Chebulagic acid throughout the simulation. These compounds have been shown to have lower RMSD values as compared with Remdesivir. Among these compounds, EBDGp was observed to have the lowest RMSD value below 1.5 Å where Remdesivir deviates within the RMSD range of 1.5 to 2.0 Å suggesting higher stability as compared with other phytocompounds and Remdesivir ( Figure 5B ). However, it can be seen that the receptor is least stable when in complex with Remdesivir, as shown by its highest RMSD ( Figure 5A) . The Ligand RMSD analysis shows that there is a less deviation of EBDGp (RMSD-1.25 Å) from the binding pocket of the receptor thus showing its role in the overall higher stability of SARS-CoV-2 RdRp as compared with Remdesivir (RMSD-2.25 Å). To further understand the dynamics of the backbone atoms, the RMSF values were calculated for backbone atoms at each point of the trajectories. Higher RMSF values indicate greater flexibility during the MD simulation. 84 Low RMSF values (<1.5 Å) of active site residues for all SARS-CoV-2 RdRp-complexes indicate their higher stabilities during the entire MD simulation (Figure 6 ), signifying that there are no major conformational changes seen in the binding pocket of the SARS-CoV-2 RdRp in complex with all the compounds. Also, RMSF values for the binding pocket residues (C-α atoms) were summarized in Table 3 . The RMSF of the SARS-CoV-2 RdRp binding site residues, upon binding of lead compounds, ie, EBDGp, Marmelide Seselin Pedunculagin, and Chebulagic acid is lower than 1 Å, which suggest that the SARS-CoV-2 RdRp binding pocket is more stable with minimum fluctuation during the 100 ns MD simulation. To understand the binding pocket stability, the MD trajectories captured for all systems were superimposed and analyzed using the Simulation Event Analysis tool of Desmond. Figures 7 and 8 show interacting residues of SARS-CoV-2 RdRp with the Remdesivir and lead compound during 100 ns MD simulation. Some residues make more than one specific contact with the ligand, which is represented by a darker shade of orange, according to the scale to the right of the plot. Remdesivir, Pedunculagin, EBDGp, and Chebulagic acid are showing good interactions with some active site residues. To reveal the binding stability between SARS-CoV-2 RdRp and phytocompounds like EBDGp, Marmelide, Seselin, Pedunculagin, and Chebulagic acid, hydrogen bond monitoring was done using the resulting MD trajectories from 100 ns simulation via Simulation Event Analysis module of Maestro. The plots of the hydrogen bonding profile are presented in Figure 9 . In Figure 9A , it is observed that the Remdesivir is making 8 hydrogen bonds with the active site residues of SARS-CoV-2 RdRp such as Asp623 and Asp760 throughout the 100 ns MD simulation. Interestingly, EBDGp is making the highest contacts, ie, about 10 hydrogen bonds as shown in Figure 9B till 25 ns and is observed to maintain 8-10 hydrogen bond contacts with the active site residues such as Asp452, Cys622, Asp623, and Thr680 throughout the 100 ns MD simulation, which is observed to be the highest as compared with a reference compound, Remdesivir. Whereas, other phytocompounds such as Marmelide, Seselin, Pedunculagin, and Chebulagic acid are shown to have the least hydrogen bond contacts as compared with all the other phytocompounds as can be seen in Figure 9C to F, respectively. From these results, we can conclude that the screened lead compound EBDGp forms a stable complex with SARS-CoV-2 RdRp and thus obtain the complex stability during the 100 ns MD simulation. It is observed that EBDGp is making the highest hydrogen bond contacts for a longer period with most of the active site residues as compared with Remdesivir and other phytocompounds which shows its greater potential in inhibiting SARS-CoV-2 RdRp as compared with Remdesivir. To compute binding free energies (ΔG Bind) of protein-ligand complexes, MM/GBSA calculations were performed, which gives the output in the context of VDW, hydrophobic, and solvation components. The phytocompounds and reference compound, Remdesivir, within the binding cavity of SARS-CoV-2 RdRp was subjected to ensemble-averaged Prime MM/GBSA calculations. The resulting free binding energies of each complex by taking ensemble-average MM/GBSA are summarized in Table 4 . Based on the MM/GBSA values obtained, as reported in Table 4 , the EBDGp is expected to have a strong binding affinity (−66.498 kcal/mol) as compared with other phytocompounds. However, Remdesivir shows binding free energy of −49.492 kcal/mol with SARS-CoV-2 RdRp. The binding free energy calculation signifies that EBDGp has the most favorable binding free energy (−66.498 kcal/mol) closely followed by Marmelide (−57.145 kcal/mol) as shown in Figure 10 . Binding free energy calculations of the compounds reveal that EBDGp forms a stronger and highly stable complex with the SARS-CoV-2 RdRp and all computed energies are found to be thermodynamically favorable as compared with other both phytocompounds and Remdesivir. To evaluate the importance of features involved in the strong binding affinity of EBDGp toward SARS-CoV-2 RdRp, we have performed molecular field-based similarity analysis using FieldTemplater software. It provides the necessary 3D-molecular field properties of the EBDGp in alignment with the Reference molecule, Remdesivir. FieldTemplater took Table 3 . RMSF values of the amino acids (c-α atoms) which are involved in the SaRS-coV-2 RdRp binding pocket after binding of lead compounds and Remdesivir. 2 compounds EBDGp and Remdesivir, optimally aligned their conformer fields, and yielded 89 templates ranked as per incorporated score (Structural similarity, Field similarity, and Shape similarity). The top-ranked molecular field template is presented in Figure 11 . In this study, we have explored only common fields with the aligned templates of Remdesivir and EBDGp describing electrostatic (positive and negative), Hydrophobic, and van der Waals properties. Large points indicating strong interactions as observed in field point patterns, 85 abundant in positive and negative electrostatic fields were observed in Remdesivir and EBDGp. Positive electrostatic fields are seen along with the hydroxyl group of EBDGp and amino group of Remdesivir, but interestingly the positive electrostatic fields are seen in large points along with hydroxyl groups of EBDGp. Large points of negative electrostatic fields are observed along with the carbonyl group and hydrophobic field along the benzoyl aromatic ring of both Remdesivir and EBDGp. Moreover, Van der Waals fields are also abundant along both Remdesivir and EBDGp equally. These results of molecular field-based similarity analysis show that positive electrostatic fields are largely observed along with hydroxyl groups of the EBDGp which indicates the importance of hydroxyl groups in an efficient binding with SARS-CoV-2 RdRp. This study demonstrates that the presence of the hydroxyl group can be assessed further for lead optimization and the design of a more potent lead candidate. The present study aimed to test the inhibition potency of 5 phytocompounds from Phyllanthus emblica and Aegel marmelos against SARS-CoV-2 RdRp using a computational approach. Molecular docking studies were conducted to compare binding affinities toward SARS-CoV-2 RdRp. Only one of them (EBDGp) showed higher FlexX docking energy values than other phytocompounds and a reference molecule, Remdesivir. EBDGp showed strong hydrogen bond interaction with key amino acid residues Thr556 and Asp623. The RMSD and RMSF profiles corresponding to the SARS-CoV-2 RdRp-EBDGp complex suggested that it is highly stable and experienced low conformational fluctuations. The RMSF values for the binding pocket residues of SARS-CoV-2 RdRp, upon binding of phytocompounds, was lower than 1.0 Å, thus it reveals that the binding pocket of SARS-CoV-2 RdRp was very stable during MD simulations. The protein-ligand interaction profile analysis revealed that EBDGp exhibited good interactions with the surrounded amino acid residues throughout the simulation. The pre (after docking study) and post (after MD simulation) MM/GBSA analysis of the EBDGp showed higher binding affinity than other phytocompounds and Remdesivir with SARS-CoV-2 RdRp target protein. Furthermore, molecular field-based similarity profiling has supported our study in the validation of the importance of the presence of hydroxyl groups in EBDGp, increasing its binding affinity with crucial amino acid residues, Thr556 and Asp623 of SARS-CoV-2 RdRp. This novel concept has percolated the new idea to design and develop promising drugs for effective binding resulting in inhibition of SARS-CoV-2 RdRp. Our study proves that the EBDGp can be a promising SARS-CoV-2 RdRp inhibitor by contributing to its stable and better interactions in the binding pocket of SARS-CoV-2 RdRp for a greater simulation time than that of Remdesivir. Our study for the first time reports the stabilized interactions of phytocompound EBDGp from Phyllanthus emblica with SARS-CoV-2 RdRp and confirmed the role of EBDGp as an anti-SARS-CoV-2 RdRp drug for treating the COVID-19. In the future, biological evaluation can be done to access the therapeutic potential of EBDGp against SARS-CoV-2 RdRp for proceeding to clinical trials. Cross-species transmission of the newly identified coronavirus 2019-nCoV CoV-2: an emerging coronavirus that causes a global threat Caffeine and caffeine-containing pharmaceuticals as promising inhibitors for 3-chymotrypsin-like protease of SARS-CoV-2 The proximal origin of SARS-CoV-2 TMPRSS2 and furin are both essential for proteolytic activation and spread of SARS-CoV-2 in human airway cells and provide promising drug targets The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein Inhibition of SARS-CoV-2 (previously 2019-nCoV) infection by a highly potent pan-coronavirus fusion inhibitor targeting its spike protein that harbors a high capacity to mediate membrane fusion Fusion mechanism of 2019-nCoV and fusion inhibitors targeting HR1 domain in spike protein Prediction of small molecule inhibitors targeting the severe acute respiratory syndrome coronavirus-2 RNAdependent RNA polymerase Comparative analysis based on the spike glycoproteins of SARS-CoV2 isolated from COVID 19 patients of different countries Quantitative structure-activity relationship and molecular docking revealed a potency of anti-hepatitis C virus drugs against human corona viruses A review on the cleavage priming of the spike protein on coronavirus by angiotensin-converting enzyme-2 and furin Zika viral polymerase inhibition using anti-HCV drugs both in market and under clinical trials Zika virus: novel guanosine derivatives revealed strong binding and possible inhibition of the polymerase Novel guanosine derivatives as anti-HCV NS5b polymerase: a QSAR and molecular docking study Molecular modeling comparison of the performance of NS5b polymerase inhibitor (PSI-7977) on prevalent HCV genotypes IDX-184 is a superior HCV direct acting anti-viral drug: a QSAR study Molecular dynamics simulation revealed binding of nucleotide inhibitors to ZIKV polymerase over 444 nanoseconds Molecular dynamics and docking reveal the potency of novel GTP derivatives against RNA dependent RNA polymerase of genotype 4a HCV Molecular modeling and docking revealed superiority of IDX-184 as HCV polymerase inhibitor Applications of computer-aided approaches in the development of hepatitis C antiviral agents Drug repurposing for viral infectious diseases: how far are we? Novel inhibitors of hepatitis C virus RNAdependent RNA polymerases Potent allosteric dengue virus NS5 polymerase inhibitors: mechanism of action and resistance profiling Development and validation of RdRp screen, a crystallization screen for viral RNA-dependent RNA polymerases Evaluating the 3C-like protease activity of SARS-coronavirus: recommendations for standardized assays for drug discovery The genome sequence of the SARS-associated coronavirus Mechanisms and enzymes involved in SARS coronavirus genome expression Sofosbuvir as a potential alternative to treat the SARS-CoV-2 epidemic Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors The future of antivirals: broad-spectrum inhibitors Updated approaches against SARS-CoV-2 Discovery and synthesis of a phosphoramidate prodrug of a pyrrolo[2,1-f][triazin-4-amino] adenine C-nucleoside (GS-5734) for the treatment of Ebola and emerging viruses Coronavirus susceptibility to the antiviral Remdesivir (GS-5734) is mediated by the viral polymerase and the proofreading exoribonuclease. mBio Initiation, extension, and termination of RNA synthesis by a paramyxovirus polymerase Therapeutic efficacy of the small molecule GS-5734 against Ebola virus in rhesus monkeys Mechanism of inhibition of Ebola virus RNA-dependent RNA polymerase by Remdesivir Fact Sheet for Health Care Providers: Emergency Use Authorization (EUA) of Remdesivir (GS-5734TM) Natural product-derived phytochemicals as potential agents against coronaviruses: a review Antiviral activity of Chamaecyparis lawsoniana extract: study with herpes simplex virus type 2 Anti-HIV activity against immunodeficiency virus type 1 (HIV-I) and type II (HIV-II) of compounds isolated from the stem bark of Combretum molle Herbal medicines for sexually transmitted diseases and AIDS Molecular targets of anti-HIV-1 triterpenes Inhibition of hepatitis B virus by an aqueous extract of Agrimonia eupatoria L Anti-HIV, anti-poxvirus, and anti-SARS activity of a nontoxic, acidic plant extract from the Trifollium species Secomet-V/anti-vac suggests that it contains a novel broad-spectrum antiviral In vitro antiviral activity of bael (Aegle marmelos Corr) upon human coxsackieviruses B1-B6 Enterovirus type 71 infections: a varied clinical pattern sometimes mimicking paralytic poliomyelitis In vitro anti-herpes simplex virus activity of 1, 2, 4, 6-tetra-O-galloyl-β-D-glucose from Phyllanthus emblica L Antiviral activity of Seselin from Aegle marmelos against nuclear polyhedrosis virus infection in the larvae of silkworm, Bombyx mori New alkaloids from Aegle marmelos Some alkaloids from Aegle marmelos Biological activities of crude extracts of chemical constituents of Bael, A. marmelos (L) Corr 555 Medicinal Plants: Field and Laboratory Manual (Identification With Its Phytochemical and In Vitro Studies Data) Current status of natural products from plants as anti-herpes simplex virus 1 agents Extracts and molecules from medicinal plants against herpes simplex viruses Bioactivity-guided fractionation of Phyllanthus orbicularis and identification of the principal anti HSV-2 compounds In vitro screening antiviral activity of Thai medicinal plants against porcine reproductive and respiratory syndrome virus Computational studies towards identification of lead herbal compounds of medicinal importance for development of nutraceutical against COVID-19 Inhibition of multiple SARS-CoV-2 proteins by an antiviral biomolecule, Seselin from Aegle marmelos deciphered using molecular docking analysis Structural basis for RNA replication by the SARS-CoV-2 polymerase PubChem periodic table and element pages: improving access to information on chemical elements from authoritative sources Penciclovir and Anidulafungin bind nsp12, which governs the RNA-dependent-RNA polymerase activity of SARS-CoV-2, with higher interaction energy than Remdesivir, indicating potential in the treatment of Covid-19 Prediction of binding constants of protein ligands: a fast method for the prioritization of hits obtained from de novo design or 3D database search programs Comparative study of several algorithms for flexible ligand docking Application of docking and active site analysis for enzyme linked biodegradation of textile dyes A fast flexible docking method using an incremental construction algorithm WISDOM-II: screening against multiple targets implicated in malaria using computational grid infrastructures Scalable algorithms for molecular dynamics simulations on commodity clusters Pharmacophore based virtual screening, molecular docking, molecular dynamics and MM-GBSA approach for identification of prospective SARS-CoV-2 inhibitor from natural product databases Constant pressure-constant temperature molecular dynamics: a correct constrained NPT ensemble using the molecular virial Molecular interaction studies of Deguelin and its derivatives with Cyclin D1 and Cyclin E in cancer cell signaling pathway: the computational approach Development and testing of the OPLS all atom force field on conformational energetics and properties of organic liquids A hierarchical approach to all-atom protein loop prediction MM-GBSA binding free energy of some novel isoxazole substituted 9-anilinoacridines as HER2 inhibitors targeting breast cancer Identification of potential mycolyltransferase Ag85C inhibitors of Mycobacterium tuberculosis H37Rv via virtual high throughput screening and binding free energy studies Potential therapeutic use of corticosteroids as SARS CoV-2 main protease inhibitors: a computational study Molecular docking and simulation studies on SARS-CoV-2 M pro reveals Mitoxantrone, Leucovorin, Birinapant, and Dynasore as potent drugs against COVID-19 Exploring conformational changes of PPAR-Ɣ complexed with novel kaempferol, quercetin, and resveratrol derivatives to understand binding mode assessment: a small-molecule checkmate to cancer therapy From arylamine N-acetyltransferase to folate-dependent acetyl CoA hydrolase: impact of folic acid on the activity of (HUMAN)NAT1 and its homologue (MOUSE)NAT2 Molecular field extrema as descriptors of biological activity: definition and validation Molecular field technology applied to virtual screening and finding the bioactive conformation A multi-stage virtual screening of FDA-approved drugs reveals potential inhibitors of SARS-CoV-2 main protease A selective pharmacophore model for beta(2)-adrenoceptor agonists The authors confirm contribution to the paper as follows: Study conception and design: MD; data collection: KP, KBL; analysis and interpretation of results: KBL, KP, MD; draft manuscript preparation: KP, KBL, MD; analysis of the data and critical reading and comments: KVS, SN. All authors reviewed the results and approved the final version of the manuscript. All authors approved the final version of the manuscript and the authorship list. Kiran Bharat Lokhande https://orcid.org/0000-0001-6945 -8288 Manjusha Dake https://orcid.org/0000-0001-5578-4939 All data generated or analyzed during this study are included in this article.