key: cord-0950975-dmruvypg
authors: Vivek-Ananth, R. P.; Krishnaswamy, Sankaran; Samal, Areejit
title: Potential phytochemical inhibitors of SARS-CoV-2 helicase Nsp13: a molecular docking and dynamic simulation study
date: 2021-06-12
journal: Mol Divers
DOI: 10.1007/s11030-021-10251-1
sha: b4d569f6ca0b55973104ef06577a394b832ea6f1
doc_id: 950975
cord_uid: dmruvypg

The SARS-CoV-2 helicase Nsp13 is a promising target for developing anti-COVID drugs. In the present study, we have identified potential natural product inhibitors of SARS-CoV-2 Nsp13 targeting the ATP-binding site using molecular docking and molecular dynamics (MD) simulations. MD simulation of the prepared crystal structure of SARS-CoV-2 Nsp13 was performed to generate an ensemble of structures of helicase Nsp13 capturing the conformational diversity of the ATP-binding site. A natural product library of more than 14,000 phytochemicals from Indian medicinal plants was used to perform virtual screening against the ensemble of Nsp13 structures. Subsequently, a two-stage filter, first based on protein–ligand docking binding energy value and second based on protein residues in the ligand-binding site and non-covalent interactions between the protein residues and the ligand in the best-docked pose, was used to identify 368 phytochemicals as potential inhibitors of SARS-CoV-2 helicase Nsp13. MD simulations of the top inhibitors complexed with protein were performed to confirm stable binding, and to compute MM-PBSA based binding energy. From among the 368 potential phytochemical inhibitors, the top identified potential inhibitors of SARS-CoV-2 helicase Nsp13 namely, Picrasidine M, (+)-Epiexcelsin, Isorhoeadine, Euphorbetin and Picrasidine N, can be taken up initially for experimental studies. [Image: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11030-021-10251-1.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an enveloped positive-sense RNA virus belonging to the family Coronaviridae that is the causative pathogen of the ongoing coronavirus disease 2019 (COVID- 19) pandemic [1] . In humans, the clinical features of this new disease ranges from asymptomatic state to acute respiratory distress syndrome [2] . Several small molecules with strong in vitro antiviral activity against SARS-CoV-2 have been identified [3, 4] and are under different stage of clinical trials. Remdesivir, one such small molecule, was granted emergency use authorization by US Food and Drug Administration (FDA) for the treatment of COVID-19 [5] . There is an increasing need to develop anti-COVID drugs, especially, those targeting the key proteins of SARS-CoV-2.

The SARS-CoV-2 genome has been reported to have 14 open reading frames which encode for the 27 viral proteins [6] . The 5'-terminus of the genome has two genes orf1ab and orf1a which code for two polyproteins pp1ab and pp1a, respectively. pp1ab and pp1a are cleaved by two proteases namely, papain-like protease (PLpro or Nsp3) and 3C-like cysteine protease (3CLpro or Nsp5), into 15 non-structural proteins, Nsp1 to Nsp10 and Nsp12 to Nsp16 [6] . In contrast, the 3'-terminus of the genome codes for the 4 structural proteins, namely, spike (S), envelope (E), matrix (M) and nucleocapsid (N), and 8 accessory proteins (3a, 3b, p6, 7a, 7b, 8b, 9b and orf14) [6] . The 15 non-structural proteins in SARS-CoV-2 comprise the viral replication and transcription complex which is essential for the coronavirus life cycle [7] . Of the 15 non-structural proteins, papain-like protease (PLpro or Nsp3), 3C-like cysteine protease (3CLpro or Nsp5), RNA-dependent RNA polymerase (RdRp or Nsp12) and helicase (Nsp13) are among the highly studied targets for identifying anti-COVID drugs [8] .

SARS-CoV-2 helicase Nsp13 has both ATPase and helicase activity, as it unwinds the RNA helices in an ATPdependent manner [9] . Notably, due to its high sequence conservation across the coronavirus family, Nsp13 is considered an attractive target for the development of antiviral drugs [10, 11] . Also, it was shown that SARS-CoV-2 helicase Nsp13 can hydrolyze all types of NTPs including ATP to unwind the RNA helices [9] . Therefore, the known ATP-binding site of the helicase Nsp13 is a promising target for effective inhibition. In this direction, the recently deposited crystal structure of SARS-CoV-2 helicase Nsp13 (PDB 6ZSL) has made development of anti-COVID drugs via targeting of Nsp13 more viable. Notably, the Nsp13 of SARS-CoV and SARS-CoV-2 share 99.8% sequence identity [11, 12] . Moreover, similar to the SARS-CoV Nsp13 structure [10] , the SARS-CoV-2 Nsp13 adopts a triangular pyramid shape with five domains namely, the RecAlike domains 1A and 2A, the 2B domain, the zinc-binding domain (ZBD) and the stalk domain ( Fig. 1) . Thus, potential drugs against SARS-CoV-2 Nsp13 may have broad activity against β-coronaviruses.

Given the urgent need to identify anti-COVID drugs, several computational studies in the past year have undertaken virtual screening of small molecule libraries against key SARS-CoV-2 non-structural proteins especially 3CLpro, RdRp and PLpro [13] [14] [15] . Some screening studies [12, [16] [17] [18] [19] have also shown that SARS-CoV-2 helicase Nsp13 is a good target for potential inhibitors. In this work, we have identified through computational approaches potential natural product inhibitors of SARS-CoV-2 helicase Nsp13 which can target the ATP-binding site. Fig. 1 Cartoon representation of the prepared crystal structure of SARS-CoV-2 helicase Nsp13 (PDB 6ZSL). The figure shows the five domains in the Nsp13 structure, namely, zinc binding domain (ZBD) colored in red, the stalk domain colored in yellow, the 1B domain colored in green, the RecA-like domains 1A and 2A colored in blue and orange, respectively. The ATP-binding site of the Nsp13 with six key residues involved in ATP hydrolysis is shown in expanded view Plant-based natural products are a major contributor to the existing approved drug space [20] . Further, medicinal plants are integral to the traditional Indian systems of medicine including Ayurveda, Siddha and Unani, and Indian herbs have been used to treat several human diseases for centuries. Recently, several herbal medicines have been proposed for treating COVID-19 [21, 22] . Previously, we have created IMPPAT [23] , the largest database on phytochemicals of the Indian medicinal plants. In the present work, we have used a small molecule library of 14,011 phytochemicals from Indian medicinal plants, mainly compiled from IMPPAT database [23] , to perform virtual screening with the goal of identification of potential phytochemical inhibitors of SARS-CoV-2 helicase Nsp13 targeting the ATP-binding site. Picrasidine M, (+)-Epiexcelsin, Isorhoeadine, Euphorbetin and Picrasidine N identified here as the top potential phytochemical inhibitors of SARS-CoV-2 helicase Nsp13, to start with, can be taken up for experimental studies.

We have earlier built the IMPPAT database [23] (https:// cb. imsc. res. in/ imppat) which is the largest resource on phytochemicals of Indian medicinal plants or herbs to date. In this work, we have used a ligand library of 14,011 phytochemicals, and this screening library was obtained by combining phytochemicals in IMPPAT [23, 24] with additional data from other literature sources [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] . Subsequently, an evaluation of these 14,011 phytochemicals based on the widelyused drug-likeness filter, Lipinski's rule of five (RO5) [41] , led to a subset of 10,510 phytochemicals that passed the R05 filter. Thereafter, we retrieved the three-dimensional (3D) structures of these 10,510 drug-like phytochemicals from PubChem [42] for virtual screening. Lastly, the 3D structures of the 10,510 phytochemicals were energy-minimized using OpenBabel [43] and converted to .pdb format from .sdf format.

In this investigation, we have used the crystal structure with 1.94 Å resolution (PDB 6ZSL) of the SARS-CoV-2 helicase Nsp13 which was downloaded from Protein Data Bank (RCSB PDB; https:// www. rcsb. org/). As of 19 February 2021, there are more than 50 crystal structures of SARS-CoV-2 Nsp13 in PDB. However, the crystal structure 6ZSL selected for this investigation has been validated by Wlodawer et al. [44] , and furthermore, is the only SARS-CoV-2 Nsp13 structure not from Pan-Dataset Density Analysis (PanDDA). In the crystal structure 6ZSL for the SARS-CoV-2 Nsp13, we have gap-filled the structural coordinates for three missing amino acid residues (339-341) by aligning to the crystal structure for the SARS-CoV Nsp13 (PDB 6JYT). Subsequently, the gap-filled structure of SARS-CoV-2 Nsp13 was energy-minimized using the minimize structure utility of UCSF Chimera [45] . Figure 1 shows the cartoon representation of this prepared crystal structure of SARS-CoV-2 Nsp13 wherein six important residues in the ATP-binding site have been highlighted [10, 12] . The SARS-CoV-2 Nsp13 adopts a triangular pyramid shape with five domains similar to the SARS-CoV Nsp13 (Fig. 1) . The triangular base of the pyramid is formed by three domains in SARS-CoV-2 Nsp13, the 2B domain and the two RecA-like domains 1A and 2A. The remaining two domains in SARS-CoV-2 Nsp13, the N-terminal zinc binding domain (ZBD) and the stalk domain, are directed toward the apex of the pyramid.

In this work, we have performed molecular dynamics (MD) simulations of SARS-CoV-2 Nsp13 protein structures using GROMACS 5.1.5 [46] with GROMOS96 54a7 force field [47] . Specifically, MD simulations were performed for: (a) the prepared crystal structure of SARS-CoV-2 Nsp13 (uncomplexed protein), and (b) the protein-ligand complexes of SARS-CoV-2 Nsp13 with predicted top five phytochemical inhibitors. The topology of the top five phytochemical inhibitors of SARS-CoV-2 Nsp13 was generated using the automated topology builder (ATB) version 3.0 (https:// atb. uq. edu. au/) [48] .

To prepare the system for MD simulation, the uncomplexed protein or the protein-ligand complex was placed at the center of a cubic box with periodic boundary conditions, and thereafter, the system was solvated by adding simple point-charge (SPC) water and neutralized by adding thirteen chloride (Cl − ) ions. Next, the system was energy-minimized using the steepest descent algorithm. Then, the solvated and neutralized system was subjected to a NVT simulation of 1 ns with 2 fs time step, wherein the number of particles, volume and temperature of the system is kept constant. Next, the pressure of the system was equilibrated to 1 bar during a NPT simulation of 1 ns with 2 fs time step, wherein the number of particles, pressure, and temperature of the system is kept constant. The position of both the protein and ligand were restrained during the above-mentioned NVT and NPT simulations.

Thereafter, a production MD run was performed, after removing the position restraint on protein and ligand, for a period of: (a) 100 ns with 2 fs time step for the SARS-CoV-2 Nsp13 uncomplexed protein, and (b) 50 ns with 2 fs time step for the SARS-CoV-2 Nsp13 protein in complex with each of the top five phytochemical inhibitors. During the production MD run, the structural coordinates of the system were written to the trajectory file after every 2 ps. Further, the system temperature and pressure were maintained at 310 K and 1 bar using the v-rescale temperature [49] and Parrinello-Rahman pressure coupling method [50] , respectively. Other parameters for the production MD simulation were fixed as per our recently published study [24] .

Finally, using GROMACS scripts and the trajectories obtained from the MD simulations, we have computed the following quantities: (a) radius of gyration (R g ) of the protein, (b) root mean square deviation (RMSD) of the C α atoms of the protein backbone, (c) root mean square fluctuation (RMSF) of the C α atoms of the protein backbone, (d) root mean square deviation (RMSD) of the heavy atoms of the ligand, and (e) distance between the center of mass of the ligand and center of mass of the key binding site residues of the protein.

The trajectory from the MD simulation of the prepared crystal structure of SARS-CoV-2 Nsp13 protein without any bound ligands was used to generate an ensemble of protein structures that capture the conformational diversity of the ATP-binding site residues of the helicase. The MD simulation trajectory from 60 to 100 ns was subjected to geometric clustering using the Daura algorithm [51] based on the conformations of the ATP-binding site residues within 8 Å from the two phosphate ions in the crystal structure of SARS-CoV-2 Nsp13. Using a clustering cutoff of 1.2 Å, we were able to identify 23 clusters from the SARS-CoV-2 Nsp13 protein simulation, from which we selected ten representative mid-point structures corresponding to the top ten populated clusters accounting for 97.6% of the sampled conformations of the ATP-binding site residues. Thus, along with the prepared crystal structure of SARS-CoV-2 Nsp13, ten structures based on clustering from the MD simulation trajectory of the protein were selected for virtual screening.

Molecular docking of the energy-minimized 3D structures of phytochemicals in the ligand library against the structures of SARS-CoV-2 Nsp13 was performed using AutoDock Vina [52] . The structures of SARS-CoV-2 Nsp13 considered for molecular docking were: (a) the prepared crystal structure of Nsp13, and (b) ten selected structures from geometric clustering of the MD simulation trajectory between 60 and 100 ns for the Nsp13 protein.

To perform docking, the structures of small molecule ligands and protein in.pdb file format were converted to .pdbqt format using prepare_ligand4.py and prepare_receptor4.py scripts, respectively, from AutoDock-Tools [53] . Thereafter, the search space center and dimension for the protein-ligand docking in AutoDock Vina for the above-mentioned eleven structures of SARS-CoV-2 Nsp13 were manually determined by considering the residues in the ATP-binding site. Subsequently, protein-ligand docking was performed using AutoDock Vina by setting the exhaustiveness parameter to 8.

Using pdb-tools [54] , we have prepared a combined structure of the best-docked conformation of a phytochemical obtained from AutoDock Vina (with the lowest binding energy) and the prepared crystal structure of SARS-CoV-2 Nsp13. Thereafter, the combined structure file of a protein-ligand complex was analyzed for non-covalent interactions between ligand and protein residues using a custom script as described in our earlier work [24] . The non-covalent interactions identified were hydrogen bonds and hydrophobic interactions. The atoms of protein and ligand are reported to have hydrophobic interactions if the distance between a carbon atom of protein (or ligand) and a carbon, halogen or sulfur atom of ligand (or protein) is ≤ 4 Å. Note that the hydrogen, chalcogen or halogen bonds identified between the atoms of protein and ligand are not reported again as hydrophobic interactions [24] .

Based on the computed Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties using SwissADME [55] and vNN-ADMET [56] , we have evaluated the pharmacokinetic properties and potential toxicity of the predicted phytochemical inhibitors for SARS-CoV-2 Nsp13 from this study.

In this study, the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) method was used to compute the binding energy of only the top five phytochemical inhibitors of SARS-CoV-2 Nsp13. From the 50 ns MD simulation trajectories for each of the protein-ligand complex, trajectories were extracted at an interval of 1 ns between 20 and 50 ns after equilibration. These were used to calculate the binding energy using g_mmpbsa [57, 58] of the phytochemical to SARS-CoV-2 Nsp13.

The helicase Nsp13 is one of the 15 non-structural proteins coded by the SARS-CoV-2 genome and shares a high sequence conservation across the family Coronaviridae [9, 12] . This protein is a critical component of the replication and transcription complex and it unwinds the RNA helices in a ATP-dependent manner similar to the SARS-CoV helicase Nsp13 [9] [10] [11] .

In this work, we have identified, using the largest small molecule library of phytochemicals from Indian herbs, potential phytochemical inhibitors of SARS-CoV-2 helicase Nsp13 that can target the ATP-binding site (Fig. 2) . From the small molecule library of 14,011 phytochemicals which was compiled based on information in our IMPPAT database [23, 24] , we identified 10,510 phytochemicals as drug-like (Methods). Figure 1 shows the prepared crystal structure of SARS-CoV-2 Nsp13 with emphasis on the six key residues in the ATP-binding site. Previously, Jia et al. [10] had identified six residues, namely, K288, S289, D374, E375, Q404 and R567, in SARS-CoV Nsp13 to be crucial for ATP hydrolysis. These six key ATP-binding site residues are also conserved in SARS-CoV-2 Nsp13. SARS-CoV-2 Nsp13 and SARS-CoV Nsp13 share 99.8% sequence identity with only one varying residue which is away from the ATP-binding site and RNA-binding site [11, 12] .

In the first stage of filtering of ligands for potential inhibitors, we have considered the binding energy of the phytochemicals with the ensemble structures, as described in the Methods section, of SARS-CoV-2 Nsp13 which represented the conformational landscape of the binding site. To decide on the binding energy cutoff for this filtering, we used the known experimental information on binding from the twelve PanDDA co-crystallized structures of SARS-CoV-2 Nsp13 with ligands bound to the ATP-binding site, namely, PDB 5RM2, 5RM7, 5RLW, 5RL9, 5RLO, 5RLY, 5RLJ, 5RLI, 5RLV, 5RLR, 5RLN and 5RLS. The binding energy of the co-crystallized ligand in the experimentally observed binding pose in the respective co-crystallized structure of SARS-CoV-2 Nsp13 was estimated using score_only option of AutoDock Vina (Supplementary Table S1 ). From the analysis of the twelve PanDDA co-crystallized structures of SARS-CoV-2 Nsp13, the best binding energy value of − 5.14 kcal/mol obtained for ligand UXG in the structure PDB 5RM2 was used as the cutoff to filter phytochemicals in the protein-ligand docking. Overall, this first stage of filtering of ligands based on known experimental evidence led to a filtered set of 5260 phytochemicals.

In the second stage of filtering of ligands for potential inhibitors, we used the interactions with the residues in the ATP-binding site based on molecular docking performed with the crystal structure conformation. Specifically, we checked that the phytochemical either binds to or forms non-covalent interactions with all six key ATP-binding site residues. At the end of the second stage of filtering, we identified 368 phytochemicals (H1-H368) as potential inhibitors of SARS-CoV-2 Nsp13 targeting the crucial ATP-binding site (Supplementary Table S2 ).

Interestingly, several of the Indian medicinal plants or herbs that produce the 368 phytochemicals H1 to H368 identified here as potential inhibitors of SARS-CoV-2 Nsp13 are mentioned for potential antiviral use in traditional medicine (Supplementary Table S3 ). Furthermore, an analysis of the best-docked pose with the prepared crystal structure of SARS-CoV-2 Nsp13 shows the extent of non-covalent interactions between the binding site residues of the protein and the complexed inhibitor (Supplementary Table S4 ). These are supplemented by the predicted physicochemical and ADMET properties of the 368 potential inhibitors of SARS-CoV-2 Nsp13 identified here (Supplementary Table S5 ).

The PubChem chemical identifier, chemical name, IUPAC name and chemical structure in SMILES format have been provided for the 368 potential phytochemical inhibitors (H1-H368) of SARS-CoV-2 Nsp13 identified here (Supplementary Table S2 ). In this section, a detailed description of the top ten potential phytochemical inhibitors (H1-H10) that have protein-ligand docking binding energies < − 8.5 kcal/ mol with the prepared crystal structure of SARS-CoV-2 Nsp13 (Tables 1 and 2 ) are discussed. Table 1 also provides the list of Indian herbs which can produce these top ten potential phytochemical inhibitors of SARS-CoV-2 Nsp13, whereas Table 2 provides the list of protein residues which form hydrogen bond or hydrophobic interaction with the top potential phytochemical inhibitors H1 to H10. The two-dimensional (2D) chemical structures and the hydrogen bond interactions in the protein-ligand complex of these top ten inhibitors are shown in Figs. 3 and 4, respectively. Of these only the top five, namely, H1 (Picrasidine M), H2 ((+)-Epiexcelsin), H3 (Isorhoeadine), H4 (Euphorbetin) and H5 (Picrasidine N), were examined in more detail through further MD simulations to estimate their binding energies using MM-PBSA (Table 3) .

Phytochemicals H1 (Picrasidine M) and H5 (Picrasidine N) are dimeric β-carboline-type alkaloid produced by the herb Picrasma quassioides. The herb Picrasma quassioides has been reported to have antiviral, antifungal and antiparasitic activities [25, 59, 60] . Additionally, the β-carboline alkaloids from Picrasma quassioides have been experimentally found to inhibit the RNA replication of the plant pathogen Tobacco mosaic virus (TMV) [61] . From Fig. 4a it is seen residue S289 forms 2 hydrogen bonds with H1 (C-H···O type and C-H···N type), and residues Q404 and R567 form 1 N-H···O type hydrogen bond each with H1. In case of H5, residues S289, Q404 and R567, which are among the six key ATP-binding site residues, form 1 C-H···O type, 1 N-H···O type and 2 N-H···O type hydrogen bonds with H5, respectively (Fig. 4e) . Phytochemical H2 ((+)-Epiexcelsin) is a lignan produced by the herb Litsea verticillata [62] . Extract of Litsea verticillata has been experimentally found to have antiviral and anti-HIV activities [63] . From Fig. 4b it is seen, the residues K288 and S289 form 2 hydrogen bonds (C-H···O type and C-H···N type) and 1 hydrogen bond (N-H···O type) with H2, respectively.

Phytochemical H3 (Isorhoeadine) is a rhoeadine alkaloid produced by the herb Papaver rhoeas [59] . H3 forms 3 C-H···O type hydrogen bonds with S289, E375 and Q404, 2 C-H···N type hydrogen bonds with R443 and R567, and 2 N-H···O type hydrogen bonds with K288 and Q404 (Fig. 4c) .

Phytochemical H4 (Euphorbetin) is a bicoumarin produced by the herb Euphorbia lathyris [64] . Figure 4d shows the extensive hydrogen bond network between H4 and the protein residues. Residues K288, S289, D374, E375 and R567 form 2 hydrogen bonds (C-H···O type and N-H···O type), 4 hydrogen bonds (2 of N-H···O type, 1 of C-H···O type and 1 of O-H···O type), 1 hydrogen bond (O-H···O type), 1 hydrogen bond (C-H···O type) and 2 hydrogen bonds (N-H···O type) with H4, respectively. Phytochemicals H6 (Ovigerine) and H8 (Hernandonine) are both produced by the herbs Hernandia guianensis and Hernandia nymphaeifolia. Figure 4f shows the hydrogen bonds H6 forms with residues of Nsp13. In case of H8, residues K288, S289, Q404, and R567 form 1 hydrogen bond (C-H···N type), 4 hydrogen bonds (2 of C-H···O type, 1 of N-H···O type and 1 of C-H···N type), 1 hydrogen bond (N-H···N type) and 1 hydrogen bond (C-H···N type) with H8, respectively (Fig. 4h) .

Phytochemical H7 (Cassamedine) is an oxoaporphine alkaloid produced by the herb Cassytha filiformis [25, 59] . From Fig. 4g it is seen, the residues S289, Q404 and R567 form 1 hydrogen bond (O-H···O type), 2 hydrogen bonds (N-H···O type and C-H···N type) and 1 hydrogen bond (C-H···N type) with H7, respectively.

Phytochemical H9 (Picriside A) is a glycoside produced by the herb Picris hieracioides [33] . H9 forms an extensive network of 17 hydrogen bonds with the residues of Nsp13 (Fig. 4i) .

Phytochemical H10 (Convolvidine) is a tropane alkaloid produced by the herb Convolvulus prostratus [26] . Residues S289, D374, E375 and Q404 form 2 hydrogen bonds (C-H···O type), 1 hydrogen bond (C-H···O type), 1 hydrogen bond (C-H···O type) and 2 hydrogen bonds (N-H···O type) with H10, respectively (Fig. 4j) .

In order to use the knowledge-based information from the PanDDA co-crystallized structures used for our analysis Table 2 The protein residues of SARS-CoV-2 Nsp13 that are involved in hydrogen bond and hydrophobic interaction with the top ten potential phytochemical inhibitors H1 to H10 in the best-docked pose For each protein-ligand complex the table lists the docking binding energy in kcal/mol, number of hydrogen bonds, and the residues forming hydrogen bond and hydrophobic interactions with the ligand atoms. Note that the hydrophobic interactions listed here are between the carbon atom of the protein residue and the carbon, halogen or sulfur atom of the ligand Protein-ligand complex Docking binding energy (kcal/mol)

Hydrogen bond interaction residues Hydrophobic interaction residues Nsp13-H1  − 10.2  9  P284, G285, S289, S310, Q404, R443,  R567   G285, G287, K288, S289, A312, A313,  E375, M378, G538  Nsp13-H2  − 9  7  G287, K288, S289, S310, R443, S535  P284, K288, A312, A313, as mentioned earlier, we have examined the ten PanDDA co-crystallized structures of SARS-CoV-2 Nsp13 with ligands bound near the ATP-binding site, namely, PDB 5RM2, 5RM7, 5RLW, 5RL9, 5RLO, 5RLY, 5RLJ, 5RLV, 5RLN and 5RLS, in relation to the top five phytochemicals identified as potential inhibitors of SARS-CoV-2 Nsp13. The ligands from PanDDA co-crystallized structures PDB 5RLI and 5RLR were not considered for comparison as they were found to be structurally similar to the ligand present in PDB 5RLJ. The comparison of the binding modes of the top five phytochemical inhibitors of Nsp13 and the ligands from the PanDDA co-crystallized structures reveals distinct modes of binding for the top inhibitors that are different from the ligands of the PanDDA structures (Supplementary Figure S1a-d) . Specifically, the binding mode of the top phytochemical inhibitor H1 in the ATP-binding site of Nsp13 doesn't overlap with that of the ligand N0E from the PanDDA structure PDB 5RM7, which has the highest 2D structural similarity with H1 (Supplementary Figure S1a , d). The 2D structural similarity was calculated using Tanimoto coefficient along with ECFP4 fingerprint as described in Vivek-Ananth et al. [65] . The above observations suggest the need for further experimental examination through cocrystallized structures and in vitro binding studies for the phytochemical inhibitors uncovered here. Figure  S3a-e) . This is also seen by the RMSD of the top inhibitors (H1-H5) (Fig. 5d) . The binding of the inhibitors to the protein is good as characterized by the distance between the center of masses of the inhibitors (H1-H5) and the six key ATP-binding site residues namely, K288, S289, D374, E375, Q404 and R567 (Supplementary Figure S4) . Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) method has been reported to be more accurate in predicting the protein-ligand binding energy in comparison with molecular docking [66] . The MM-PBSA based binding energy of the top five phytochemical inhibitors of SARS-CoV-2 Nsp13 (Table 3) Table 2; Supplementary  Table S4 ). In the MD simulation, Picrasidine M has exhibited stable interactions with SARS-CoV-2 Nsp13 (Fig. 5 ; Supplementary Figure S4 ). Therefore, the phytochemical Picrasidine M is likely to be a good potential inhibitor of SARS-CoV-2 Nsp13 which must be taken forward as a lead compound for experimental and clinical studies.

More than 3.7 million people have died due to COVID-19 which is likely to become endemic past the pandemic phase. Computational approaches such as molecular docking and molecular dynamics simulation can be used to accelerate the identification and development of anti-COVID drugs. SARS-CoV-2 Nsp13 is one of the key components of the viral replication and translation complex, and has been identified as one of the key targets for development of anti-COVID drugs [8, 11] . In this study, we have computationally screened a small molecule library of 14,011 phytochemicals against the ATP-binding site of the SARS-CoV-2 helicase Nsp13.

Natural products have directly or indirectly contributed to many FDA approved drugs [20] and the therapeutic activity space of the diverse natural product space is still largely unmapped and unexplored. This motivated this virtual Fig. 4 Cartoon representation of the hydrogen bond interactions in the best-docked pose of the top ten potential phytochemical inhibitors of SARS-CoV-2 helicase Nsp13. In this figure, hydrogen bond interactions are shown as yellow colored dotted lines between the residues of Nsp13 and the atoms of a H1, b H2, c H3, d H4, e H5, f H6, g H7, h H8, i H9 and j H10. The carbon atoms of the ligand are colored in green and the carbon atoms of the residues in Nsp13 are colored in cyan. The Nsp13 residues forming hydrogen bond interactions with the ligand are labeled with their one letter amino acid code and their residue number ◂ screening study involving the largest small molecule library of phytochemicals from Indian medicinal plants and the key drug target SARS-CoV-2 helicase Nsp13. We have identified 368 phytochemicals as potential inhibitors of SARS-CoV-2 Nsp13, targeting the ATP-binding site.

Notably, among the top five potential inhibitors, both Picrasidine M (H1) and Picrasidine N (H5) are found to be produced by the herb Picrasma quassioides, which has been used in traditional medicine as antiviral, antifungal and antiparasitic. The herb Litsea verticillata which can produce (+)-Epiexcelsin (H2), has been reported to have antiviral and anti-HIV activities [63] . It should be emphasized that, further experimental studies are required to validate the inhibitory potential of the phytochemical inhibitors of SARS-CoV-2 Nsp13 identified in this study. In conclusion, the potential phytochemical inhibitors of SARS-CoV-2 helicase Nsp13 identified here constitute a potential chemical library for the development of anti-COVID drugs.

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s11030-021-10251-1.

A novel coronavirus from patients with pneumonia in China

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

Identification of inhibitors of SARS-CoV-2 3CL-pro enzymatic activity using a small molecule in vitro repurposing screen

Drugs that inhibit TMEM16 proteins block SARS-CoV-2 spike-induced syncytia

US Food and Drug Administration (2020) FDA Issues emergency use authorization for potential COVID-19 treatment

Genome composition and divergence of the novel coronavirus (2019-nCoV) originating in China

Coronavirus biology and replication: implications for SARS-CoV-2

Therapeutic options for the 2019 novel coronavirus (2019-nCoV)

SARS-coronavirus-2 Nsp13 possesses NTPase and RNA helicase activities that can be inhibited by bismuth salts

Delicate structural coordination of the severe acute respiratory syndrome coronavirus Nsp13 upon ATP hydrolysis

Should we try SARS-CoV-2 helicase inhibitors for COVID-19 therapy

Discovery of COVID-19 inhibitors targeting the SARS-CoV-2 Nsp13 helicase

Computational drug discovery and repurposing for the treatment of COVID-19: a systematic review

A revisit to the research updates of drugs, vaccines, and bioinformatics approaches in combating COVID-19 pandemic

Structure-based identification of SARS-CoV-2 main protease inhibitors from anti-viral specific chemical libraries: an exhaustive computational screening approach

Structural elucidation of SARS-CoV-2 vital proteins: computational methods reveal potential drug candidates against main protease, Nsp12 polymerase and Nsp13 helicase

Structure-based virtual screening identifies multiple stable binding sites at the RecA domains of SARS-CoV-2 helicase enzyme

An in-silico study on selected organosulfur compounds as potential drugs for SARS-CoV-2 infection via binding multiple drug targets

Exploring the potential of novel phenolic compounds as potential therapeutic candidates against SARS-CoV-2, using quantum chemistry, molecular docking and dynamic studies

Natural products as sources of new drugs from 1981 to 2014

Traditional Chinese medicine for COVID-19 treatment

COVID-19 pandemic: a pragmatic plan for ayurveda intervention

IMPPAT: a curated database of Indian medicinal plants

In silico identification of potential natural product inhibitors of human proteases key to SARS-CoV-2 infection

New Delhi 26. The Wealth of India (2000) The wealth of India: a dictionary of Indian raw materials & industrial products. First supplement series

The treatise on Indian medicinal plants

The treatise on Indian medicinal plants

The treatise on Indian medicinal plants

Publications & Information Directorate

The treatise on Indian medicinal plants

Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings

PubChem 2019 update: improved access to chemical data

Open Babel: an open chemical toolbox

Ligand-centered assessment of SARS-CoV-2 drug target models in the protein data bank

UCSF Chimera-a visualization system for exploratory research and analysis

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

Definition and testing of the GROMOS force-field versions 54A7 and 54B7

Automated topology builder version 3.0: prediction of solvation free enthalpies in water and hexane

Molecular dynamics with coupling to an external bath

Polymorphic transitions in single crystals: a new molecular dynamics method

Peptide folding: when simulation meets experiment

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

AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility

pdb-tools: a swiss army knife for molecular structures. F1000Re-search

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

vNN web server for ADMET predictions

Electrostatics of nanosystems: application to microtubules and the ribosome

Open Source Drug Discovery Consortium, Lynn A (2014) g_mmpbsa-A GROMACS tool for high-throughput MM-PBSA calculations

Indian medicinal plants: an illustrated dictionary

Phytochemistry, traditional use and pharmacological activity of Picrasma quassioides: a critical reviews

Tobacco mosaic virus (TMV) inhibitors from Picrasma quassioides Benn

Natural anti-HIV agents-part I: (+)-demethoxyepiexcelsin and verticillatol from Litsea verticillata

Litsea species as potential antiviral plant sources

Euphorbetin: a new bicoumarin from Euphorbia lathyris L

MeFSAT: a curated natural product database specific to secondary metabolites of medicinal fungi

The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities

Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

Acknowledgements The authors acknowledge the use of the super-