key: cord-0720142-9cblzv98
authors: Chandra, Anshuman; Gurjar, Vaishali; Qamar, Imteyaz; Singh, Nagendra
title: Identification of potential inhibitors of SARS-COV-2 endoribonuclease (EndoU) from FDA approved drugs: a drug repurposing approach to find therapeutics for COVID-19
date: 2020-06-09
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1775127
sha: f8a48fbd787524736acc93b488201aa182190d20
doc_id: 720142
cord_uid: 9cblzv98

SARS-CoV-2 is causative agent of COVID-19, which is responsible for severe social and economic disruption globally. Lack of vaccine or antiviral drug with clinical efficacy suggested that drug repurposing approach may provide a quick therapeutic solution to COVID-19. Nonstructural protein-15 (NSP15) encodes for an uridylate-specific endoribonuclease (EndoU) enzyme, essential for virus life cycle and an attractive target for drug development. We have performed in silico based virtual screening of FDA approved compounds targeting EndoU in search of COVID-19 drugs from commercially available approved molecules. Two drugs Glisoxepide and Idarubicin used for treatment for diabetes and leukemia, respectively, were selected as stronger binder of EndoU. Both the drugs bound to the active site of the viral endonuclease by forming attractive intermolecular interactions with catalytically essential amino acid residues, His235, His250, and Lys290. Molecular dynamics simulation studies showed stable conformation dynamics upon drugs binding to endoU. The binding free energies for Glisoxepide and Idarubicin were calculated to be –141 ± 11 and –136 ± 16 kJ/mol, respectively. The IC(50) were predicted to be 9.2 µM and 30 µM for Glisoxepide and Idarubicin, respectively. Comparative structural analysis showed the stronger binding of EndoU to Glisoxepide and Idarubicin than to uridine monophosphate (UMP). Surface area calculations showed buried are of 361.8Å(2) by Glisoxepide which is almost double of the area occupied by UMP suggesting stronger binding of the drug than the ribonucleotide. However, further studies on these drugs for evaluation of their clinical efficacy and dose formulations may be required, which may provide a quick therapeutic option to treat COVID-19. Communicated by Ramaswamy H. Sarma

Novel coronavirus (nCoV) originated from Wuhan, China, is a spherical or pleomorphic shaped, single-stranded RNA virus with characteristic crown-shaped glycoproteins on its surface (Graham Carlos et al., 2020) . nCOV causes coronavirus disease-2019 (COVID-19) characterized by pneumonia and high fever (Corkery et al., 2015; Majid Rezaei, 2020) . The nCOV is also named as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) upon the similarity of its symptoms to those induced by the severe acute respiratory syndrome (SARS). The first case of COVID-19 was reported in late December 2019 and was declared pandemic on March 11, 2020. As of June 3, 2020, more than 6.28 million cases of SARS-CoV-2 infection have been confirmed with 379,941 global fatalities (https://www.who.int/emergencies/diseases/novel-coronaviru s-2019/situation-reports/). The genomic sequences of SARS-CoV-2 viruses isolated from a number of patients share sequence identity higher than 99.9%, suggesting a very recent host shift into humans (Kim, Chung, Jo, et al., 2020) .

SARS-CoV-2 belongs to the beta CoVs category and has a diameter of approximately 60-140 nm. Like other CoVs, it is sensitive to UV rays and higher temperature (Cascella et al., 2020) . Diagnosis of COVID-19 is based on the real-time polymerase chain reaction (RT-PCR) technique of isolated RNA from the patients. The SARS-CoV-2 genome consists of a long replicase gene encoding nonstructural proteins (NSPs), along with structural and accessory genes (Vijgen et al., 2005) . The replicase gene encodes for two ORFs named as rep1a and rep1b by ribosomal frameshifting, which are translated into pp1a and pp1ab polyproteins (Cui et al., 2019) . These two large polyproteins are processed by two viral proteases known as 3C-like protease (3CLpro) and papain-like protease (PLP) that are encoded by NSP5 and NSP3, respectively. Sixteen viral NSPs are produced by cleavage (Baez-Santos et al., 2015) , which gets assembled into a membranebound large replicase complex having several functions. NSPs have been exhibiting many enzymatic activities along with RNA replication and modification of subgenomic RNAs. One of these proteins, NSP15 a nidoviral RNA uridylate-specific endoribonuclease (NendoU) containing carboxy-terminal catalytic domain linked to the EndoU family (Deng & Baker, 2018) . Previous study has suggested that Nsp15 has crucial role in evasion of host defense mechanism and loss of EndoU function results in protective immune response (Deng & Baker, 2018; Liu et al., 2019) . Therefore, Nsp15 is considered critical in coronavirus biology. Detection technique(s) also remains a critical factor for treatment and may still need improvements Elasnaoui & Chawki, 2020) . There is no existing vaccine or proven drug for this disease, but various treatment options such as utilizing medicines effective in other viral ailments, are being attempted against SARS-COV2 (Babadaei, Hasan, Vahdani, et al., 2020; Elmezayen et al., 2020; Khan, Jha, Amera, et al., 2020; Mahanta et al., 2020; Muralidharan et al., 2020) . The use of de-risked compounds, with potentially lower overall development costs and shorter development timelines is being widely accepted against drug targets (Adeoye et al., 2020; Lobo-Galo et al., 2020) . Novel and natural compounds as inhibitors of SARS-CoV-2 are also being explored using regress computational approaches (Abdelli et al., 2020; Bhardwaj et al., 2020; Elfiky, 2020a; Enmozhi et al., 2020; Islam et al., 2020; Kumar et al., 2020; Umesh et al., 2020; Wahedi et al., 2020) .

In order to explore anitviral compounds several reports of using drug repurposing methods are also recently been made available (Aanouz et al., 2020; Al-Khafaji et al., 2020; Chetri et al., 2019; Elfiky, 2020b; Gupta et al., 2020; Gyebi et al., 2020; Joshi et al., 2020; Mittal et al., 2020; Sarma et al., 2020; Sinha et al., 2020) . In the present study, in silico methods have been used to identify potential drugs which may act as specific inhibitor for EndoU enzyme of SARS-CoV-2. Food and drug administration (FDA) approved drug database has been used for structure based virtual screening (SBVS) and molecular docking. Top compounds were then re-docked and the protein complexes were analyzed. To remove false-positive hits further validation was performed using molecular dynamics (MD) simulation studies. Furthermore, the molecular mechanics-Poisson-Boltzmann surface area (MM/PBSA) based binding free energies between the EndoU-ligand complexes were calculated. Glisoxepide, sulfonylurea agent used for the treatment of type 2 diabetes mellitus and Idarubicin, anticancer drug were found most potent inhibitors of SARS-COV2 EndoU.

Crystal structure of the EndoU of SARS-Cov-2 was downloaded from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank. The crystal structure (PDB ID: 6W01) has the highest resolution among all EndoU crystal structures and hence was used in the study. 6W01 is a complex of EndoU with citrate ion and was solved at 1.9 Å resolution . The bound citrate ion was used as a control for docking and simulation studies. For ligand retrieval, the FDA approved drug database was downloaded from the ZINC 12 database (Irwin & Shoichet, 2005; Irwin et al., 2012) in mol2 file format.

Protein and ligand were prepared for virtual screening. For protein preparation, only single chain 3D structure of EndoU was retained from the PDB and all other heteroatoms were removed. The coordinates were then energy minimized to avoid close contact within atoms using Chimera (Pettersen et al., 2004) . Amber ff99SB force field was employed for the minimization process with 100 steepest descent steps. Then, polar hydrogen along with Kollman charges were added to the protein structure. Finally, AD4 atom type was assigned and structure was saved in pdbqt (Protein Data Bank, Partial Charge (Q), & Atom Type (T)) format using Autodock Tools (ADT)) (Steffen et al., 2010) . For ligand preparation, prepareli-gand4.py script provided by Autodock developers was used. First, the non-polar hydrogens were added and merged, then the atoms types were set, and the Gasteiger charges were added. Finally, each ligand from the FDA drug database (n ¼ 2895) was converted to required pdbqt format.

SBVS is a technique used in the drug discovery and development methods to explore compound library in search of novel bioactive molecules against a target structure (Shukla & Singh, 2020) . SBVS was performed by using Autodock Vina (Trott & Olson, 2009 ) and then redocked with Idock (Li et al., 2012) and Smina (Koes et al., 2013) . The previously prepared protein coordinates were used to prepare the grid. The position of co-crystallized compound citrate was used as a center of grid. The active site of EndoU is comprised of His235, His250, and Lys290 residues. The citrate bound with the His235, Gln245, His250, Lys290, Val292, Thr341 and Tyr343 in the co-crystallized structure. The parameters for the grid were Center_X ¼ À63.98, Center_Y ¼ 72.101, Center_ Z ¼ 28.841 and Size_X ¼ 40 Å, Size_Y ¼ 40 Å and Size_Z ¼ 40 Å. The number of binding mode and energy range were set to maximum with exhaustiveness at 50. The prepared ligands were screened against the EndoU enzyme. To avoid any false-positive identification, we further performed SBVS with Smina tool. Same parameter and ligands were used in the process. Top five drug molecules (Dihydroergotamine, Glisoxepide, Idarubicin, Ergotamine and Tasosartan) with good binding affinity to NSP15 were selected for further analysis.

In order to calculate binding affinity of the drugs Dihydroergotamine, Glisoxepide, Idarubicin, Ergotamine and Tasosartan, the half-maximal inhibitory concentration (IC 50 ) against EndoU was predicted using Autodocktools package version (Steffen et al., 2010) . The grid box spacing was 0.375 Å with Center_X ¼ À63.98, Center_Y ¼ 72.10, Center_Z ¼ 28.84 and Size_X ¼ 40 npts, Size_Y ¼ 40 npts and Size_Z ¼ 40 npts. Number of individuals in population and maximum number of energy evaluations in each run was kept 150 and 25,000,000, respectively to determine best docking conformer and other docking parameters were set to default. Overall, 100 conformation of each compound was generated with IC 50 values. Finally, Glisoxepide, Idarubicin and Tasosartan complexes with EndoU were used for MD simulation.

To understand the conformational dynamics upon drugs binding to EndoU, MD simulation studies were performed. All atom simulation method was used to gain insights by solving Newton's equation of motion. MD simulations of the EndoU-drugs complexes were performed with the GROMACS 4.6.5 package using gromos53a6 force field (Pronk et al., 2013; Van Der Spoel et al., 2005) . The topology of EndoU was generated using pdb2gmx modules of GROMACS. The model was placed in a dodecahedron box using the editconf module keeping protein EndoU at the center and with 10 Å distance from the edges. The system was then solvated with simple point charge water model (SPC216) to attain real dynamics. Solvated systems were then energy minimized using the Steepest Descent algorithms with 5000 steps cutoff. Further, individual systems were slowly warmed up using a V-rescale thermostat with a coupling constant to reach 310 K to perform equilibration in NVT (number of atoms, volume of the system, and temperature of the system) ensemble. Later, the solvent density was maintained using a Parrinello-Rahman barostat with a coupling constant at 1 bar and 310 K to perform equilibration in NPT (number of atoms in the system, the pressure of the system and temperature of the system) ensemble. The systems were equilibrated for about 100 ps under both ensemble processes (NVT and NPT) with position restrained dynamics. The LINCS algorithm was used for constraining all the bonds. Finally the systems were submitted to molecular dynamics simulation for 100 ns to observe stability of each EndoU-ligand complex Shukla et al., 2018; Sk et al., 2020) .

The resultant trajectories were sampled using GROMACS for analysis (Van Der Spoel et al., 2005) . Protein dynamics and ligand interaction were analyzed for the whole trajectories. Respective frames of the trajectories were visualized and analyzed using Chimera (Pettersen et al., 2004) . Initially, aligning and rotational fitting of the trajectory was performed to exclude simulation artifacts if any. The root mean square deviation (RMSD) of the C a backbone atoms of the protein was calculated for each system. The root mean square fluctuation (RMSF) values of protein C a atoms were also measured to identify the fluctuating residues during simulation. To evaluate the compactness of the protein structure upon ligand binding, radius of gyration (R g ) was measured. The correlated motion of protein was also measured using principal component analysis (PCA). In PCA, covariance matrix is generated and studied by diagonalizing all the motions of protein in form of eigenvector and eigenvalues. The principal component PC1 versus PC2 was plotted to examine the cluster in phase space. Finally, the binding free energy was calculated for each protein ligand complex. In addition, Poisson-Boltzmann surface area continuum solvation (MM/PBSA) (Fu et al., 2018; Hu et al., 2017; Kumari et al., 2014; Moesgaard et al., 2020; Xue et al., 2018) was also performed to calculate the free binding energy of the ligand. Distance of 3.2 Å was kept as cutoff for bond length for salt bridge and hydrogen bond interactions (Mohammad et al., 2020; . All graphical presentation was prepared using Origin 6.0.

The structure-based virtual screening was carried out for the identification of inhibitors against EndoU using Autodock vina, Idock and Smina tools (Koes et al., 2013; Li et al., 2012; Trott & Olson, 2009) . Each molecule was scored on the basis of binding affinity. To validate screening and redocking protocol, bound citrate molecule was used as control. The binding affinity of resulted compounds ranged from -5.1 to -9.8 kcal/mol and -5.2 to -10.23 kcal/mol and -5.2 to -10.5 kcal/mol from Autodock vina, Idock and Smina, respectively (Table 1) . Idarubicin was found to be the lowest energy molecule with binding affinity -10.5 kcal/mol and Tasosartan was found having lowest binding affinity (Table 1 ). All the compounds showed a favorable binding affinity towards EndoU indicating their competency of inhibiting enzyme (Table 1) . To avoid toxicity, selected compounds were filtered by removing the drugs targeting human kinase receptors in apoptosis.

Finally, Dihydroergotamine, Glisoxepide, Idarubicin, Ergotamine and Tasosartan that showed good binding affinity for EndoU were selected (Table 2) .

The IC 50 values of selected molecules were observed in the range from 0.0092 mM to 0.52 mM. The IC 50 of Dihydroergotamine, Glisoxepide, Idarubicin, Ergotamine and Tasosartan were predicted to be 0.52 mM, 0.0092 mM, 0.039 mM, 0.19 mM and 0.09 mM, respectively (Table 1) . From the runs of prediction by Autodock, Glisoxepide was observed to be the best possible inhibitory molecule against EndoU of SARS-Cov-2 with the IC 50 value of 9.22 mM.

In order to provide dynamics information, all the four models including citrate as control (EndoU-Idarubicin, EndoU-Glisoxepide, and EndoU-Tasosartan) were subjected to MD simulations for 100 ns each. Further, to determine stability of all the systems, RMSD values of backbone atoms were analyzed. During 100 ns simulations, all the systems achieved stability after 10 ns and average RMSD values were found to be 0.372 nm, 0.393 nm, 0.385 nm and 0.4 nm, respectively (Figure 1(a) ). The RMSD profiles indicated that the binding of the all selected compounds significantly stabilized the EndoU structure. The RMSD value for individual drugs was also calculated. Average RMSD value of Citrate, Glisoxepide, Idarubicin and Tasosartan was 0.126 nm, 0.203 nm, 0.106 nm and 0.174 nm, respectively. Average RMSD value of Idarubicin was found to be lowest whereas rest of the ligand showed comparable RMSD values (Figure 1(b) ). RMSD of EndoU and its all complexes got stabilized at 50 ns of simulation. Further analysis of complex was carried over for after 50 ns (Table 3) .

To provide more detailed information of motions in EndoU upon binding of selected compounds, RMSF as function of the residue number for the last 50 ns trajectories was plotted (Figure 1(c) ). The RMSF values of EndoU-Citrate, EndoU-Idarubicin, EndoU-Glisoxepide, and EndoU-Tasosartan were observed to be 0.224, 0.228, 0.19 and 0.212 nm, respectively. As shown in RMSF plot several residual fluctuations were observed in the RMSF profile of each complex. Upon the binding of Glisoxepide to EndoU, the RMSF of binding site was significantly reduced suggesting gain in thermodynamic stability of the enzyme. Intramolecular hydrogen bonds in a protein are considered to play vital roles in molecular recognition, stability and overall conformation. The number of hydrogen bonds was analyzed to get insight into the protein-ligand interaction and stability. The possible number of hydrogen bonds in the complexes of EndoU with Citrate, Idarubicin, Glisoxepide, and Tasosartan were observed to 4, 7, 8 and 1, respectively (Figure 1(d) ). Drugs also gained more van der Waals interactions with EndoU in comparison with the control citrate compound, which has provided base for stronger binding of drugs to the enzyme. Furthermore, the R g is defined as overall conformational shape of a protein and utilized to know about protein compactness. We observed the compactness of EndoU-complexes by plotting their R g (Figure 2(a) ). The R g values of EndoU-Citrate, EndoU-Idarubicin, EndoU-Glisoxepide, and EndoU-Tasosartan were found to be 2.35, 2.41, 2.3, and 2.98 nm, respectively. The average R g values of each compound were comparable while the binding of Glisoxepide showed the decrease in R g values. The solvent accessible area was also calculated for each complex. The average solvent accessible area was 196, 199, 198 and 198 nm 2 for Idarubicin, Glisoxepide, Tasosartan and Citrate complexes, respectively (Figure 2(b) ).

PCA investigate the total expansion of protein during the MD simulation. The PCA was performed for EndoU-Citrate, EndoU-Idarubicin, EndoU-Glisoxepide, and EndoU-Tasosartan via determining their collective motions and dynamics were illustrated utilizing gmx covar and gmx anaeig modules with reference to the backbone. PCA reveals the average mobility of a protein. The first ten eigenvectors for EndoU complexes with Citrate, Idarubicin, Glisoxepide, and Tasosartan accounted for 81%, 80%, 70% and 77% of the motions observed for the last 50 ns, respectively (Figure 3(a) ). The results show that EndoU-Tasosartan showed fewer motions as compared to other predicted hits and the reference compound, whereas reference compound citrate showed higher motions among all ligands. EndoU-Idarubicin, EndoU-Glisoxepide, and EndoU-Tasosartan showed fewer motions during ligand binding, indicating their stabilizing effects on the enzyme. Since the first few eigenvectors play an important role in overall motions, we considered only the first two eigenvectors from each EndoU complex for a clear depiction of the results (Figure 3(b) ). The PC1 versus PC2 principal component graph showed that EndoU-Tasosartan showed a much defined and stable cluster as compared to all other selected ligands, while EndoU-Citrate did not show a stable cluster in phase space. So, the twodimensional PCA result analysis suggested that the predicted drugs are better than the reference citrate ion for recognizing endonuclease enzyme.

Conformational behavior and degree of folding of a protein depends on its secondary structure. We have calculated the secondary structure elements of EndoU during the course of simulation to identify the stability of EndoU upon ligand binding. a-helix, b-sheet, bridge and turn of EndoU were analyzed considering individual residues per step, and were plotted as a function of time with the total secondary structure with respect to time. The total structural elements of all EndoU complexes remain stable, whereas, increment was observed in the secondary structure contents of selected ligands (Figure 4 ). This increase in secondary structure is due to the conversion of coils into b-sheet as a result of structural stabilization of enzyme due to drugs binding. Overall, no reduction in secondary elements was observed after drugs binding to EndoU, suggesting the strong stability of the complexes.

The MM/PBSA tool is used to evaluate the protein and ligand binding energy during MD simulation. The free binding energy was determined using polar and apolar solvation energy. The free binding energy was investigated as electrostatic energy, polar solvation energy, van der Waals energy, SASA energy and average binding energy (Table 4) . MM/PBSA based free binding energy of all EndoU-ligand complexes was calculated for the last 50 ns trajectories. The reference compound citrate showed free binding energy of -54 kJ/mol whereas Glisoxepide, Idarubicin and Tasosartan showed binding energy of -141, -134 and -103 kJ/mol suggesting Glisoxepide having highest affinity for EndoU. In order to elucidate the energy contribution prediction of potentially important residues of EndoU participating in ligand binding during the course of simulation, the energy decomposition plot was calculated. The plot showed that Gly240, Gly241, Gly247, Gly248, His250, Asn278, Ser294, Trp333, Thr341, Tyr343, Pro344, Lys345 and Leu346 were consistently involved in stabilization of the EndoU complexes during simulation ( Figure 5 ).

NSP15 is shown to have endoribonuclease activity by hydrolyzing RNA with specificity towards uridylate (Bhardwaj et al., , 2008 . The catalytic site of NSP15 is located at the C-terminal NendoU domain in threedimensional crystal structure . The active site is formed by a shallow groove between the two b-sheets, with six conserved key residues His235, His250, Ser294, Lys290, Thr341, and Tyr343 (Figure 6(a) ) . Based on sequence similarities His235, His250, Lys290 have been known as catalytic triad of EndoU (Ricagno et al., 2006) .

As Idarubicin and Glisoxepide are shown to have maximum affinity for EndoU based on binding energies and IC 50 predictions, the intermolecular interactions between EndoU and both the drug molecules were analyzed. The drug complexes were stabilized by noncovalent interactions including hydrogen bonds and van der Waals interactions between drug molecules and EndoU. Both the drug molecules occupied active site of the enzyme. In EndoU-citrate complex (PDB: 6W01), citrate formed four hydrogen bond with the side chains of Gln245, Gly248, His250 and Lys290 amino acid residues of the active site. Idarubicin and Glisoxepide occupied same place as citrate ion in the active site of EndoU. Glisoxepide formed maximum eight number of hydrogen bonded interactions with Gly248 (NH), His235 (Ne2), His250 (Ne2), Asn278 (NH2), Lys290 (Nf), Ser294 (Oc), Thr341 (Ok) and Leu346 (NH), whereas Val292, Trp333, Tyr343 and Glu340 of EndoU were involved in vander Waals interaction with the Glisoxepide (Figure 6(b) ). Idarubicin bound EndoU less strongly than Glisoxepide with seven number of hydrogen bonded interactions. It formed hydrogen bonds with Gln245 (NH2), Gly248 (NH), Val292 (O) and Ser294 (NH, Oc), Lys290 (Nf), Tyr343 (OH) and also made van der Waals contacts with the side chains of Val292, Trp333, Thr341, Tyr343, and Lys345 residues of the enzyme (Figure 6(c) ). 3.8. Comparison of the binding of drugs and ribonucleotide to EndoU

Recently, crystal structure of SARS-COV2 EndoU complex with uridine mono phosphate (UMP) has been solved (6WLC.pdb) . Idarubicin and Glisoxepide both share same place with UMP in the active site of EndoU.

UMP formed only four hydrogen bonded interaction with His250 (Ne2), Ser294 (Oc, NH) and Lys290 (Nf) (Figure 6(d) ). UMP also formed van der Waal contacts with Tyr343 and Val292 residues of EndoU. Comparatively, analysis showed that both Glisoxepide and Idarubicin formed more number of attractive interactions with EndoU and hence bound more strongly than UMP. Buried surface area calculations of EndoU showed 361.8 Å 2 , 310.5 Å 2 and 186.6 Å 2 area occupied by Glixoxepide, Idarubicin and UMP, respectively further indicating more contact area for the drugs than the UMP.

We aim to contribute to tackle the therapeutic crisis arisen due to SARS-COV-2, by using commercially available FDA approved drugs for the development of novel class of antivirals. As EndoU is found to be involved in host envision mechanism and thus can be utilized as an attractive drug target for the development of anti-COVID-19 therapeutics.

We have utilized in-silico based approach by using structure based virtual screening of the FDA approved database (d) Binding of UMP to NSP15 is shown. Interacting residues of NSP15 are shown in light yellow and ligand molecules are colored brown. Some of the interacting residues have been omitted for the sake of clarity. The figure was drawn using Chimera (Pettersen et al., 2004) .

against EndoU to identify its potent inhibitors. Two drugs Idarubicin and Glisoxepide, were selected out of 2895 compounds based on their binding affinity towards EndoU. Idarubicin inserts itself into DNA and prevents DNA unwinding by interfering with the enzyme topoisomerase II in order to perform its function as antileukemic drug (Fukushima et al., 1993) . Glisoxepide is an anti-diabetic drug belongs to second-generation sulfonylureas (Seino, 2012) . Both the drugs recognized active site of EndoU suggesting them to be competitive inhibitors of the enzyme. MD simulation studies confirmed that selected drugs efficiently bind to EndoU and formed stable complexes with least structural perturbations. Glisoxepide has shown to be most potent inhibitors of EndoU with predicted IC 50 of 9.2 mM whereas Idarubicin binds enzyme with IC 50 of 30 mM. A comparative analysis also showed stronger binding of drugs than UMP and citrate ion to EndoU. We assume that Glisoxepide and Idarubicin both may have a good potential of controlling SARS-COV2 infection and are needed to be further validated by in vitro and in vivo studies. Furthermore, these compounds may also act synergistically to develop novel class of inhibitor against SARS-CoV2.

Moroccan Medicinal plants as inhibitors of COVID-19: Computational investigations

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

Repurposing of chloroquine and some clinically approved antiviral drugs as effective therapeutics to prevent cellular entry and replication of coronavirus

Using integrated computational approaches to identify safe and rapid treatment for SARS-CoV-2

The expression level of angiotensin-converting enzyme 2 determine the severity of COVID-19: Lung and heart tissue as targets

Development of remdesivir repositioning as a nucleotide analog against COVID-19 RNA dependent RNA polymerase

The SARS-coronavirus papain-like protease: Structure, function and inhibition by designed antiviral compounds

The severe acute respiratory syndrome coronavirus Nsp15 protein is an endoribonuclease that prefers manganese as a cofactor

Structural and functional analyses of the severe acute respiratory syndrome coronavirus endoribonuclease Nsp15

RNA recognition and cleavage by the SARS coronavirus endoribonuclease

Identification of bioactive molecules from tea plant as SARS-CoV-2 main protease inhibitors

Features, evaluation and treatment coronavirus (COVID-19)

Identification and characterization of glyceraldehyde 3-phosphate dehydrogenase from Fasciola gigantica

Connecting the speckles: Splicing kinases and their role in tumorigenesis and treatment response

Origin and evolution of pathogenic coronaviruses

An "Old" protein with a new story: Coronavirus endoribonuclease is important for evading host antiviral defenses

Using X-ray images and deep learning for automated detection of coronavirus disease

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

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

Drug repurposing for coronavirus (COVID-19): In silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes

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

Exploring the binding mechanism of metabotropic glutamate receptor 5 negative allosteric modulators in clinical trials by molecular dynamics simulations

Action mechanism of idarubicin (4-demethoxydaunorubicin) as compared with daunorubicin in leukemic cells

Novel Wuhan (2019-NCoV) coronavirus

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

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

Exploration of the binding affinities between ecdysone agonists and EcR/USP by docking and MM-PB/GBSA approaches

ZINC -A free database of commercially available compounds for virtual screening

ZINC: A free tool to discover chemistry for biology

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

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

Marine natural compounds as potents inhibitors against the main protease of SARS-CoV-2. A molecular dynamic study

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

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

Center for Structural Genomics of Infectious Diseases (CSGID). (2020). RCSB PDB -6WLC: Crystal structure of NSP15 endoribonuclease from SARS CoV-2 in the complex with uridine-5'-monophosphate

Identification of coronavirus isolated from a patient in Korea with covid-19

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

Lessons learned in empirical scoring with smina from the CSAR 2011 benchmarking exercise

Understanding the binding affinity of noscapines with protease of SARS-CoV-2 for COVID-19 using MD simulations at different temperatures

g_mmpbsa-a GROMACS tool for highthroughput MM-PBSA calculations

Idock: A multithreaded virtual screening tool for flexible ligand docking

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

Potential anti-viral activity of approved repurposed drug against main protease of SARS-CoV-2: An in silico based approach

Theory about treatments and morbidity prevention of corona virus disease (Covid-19)

Identification of potential molecules against COVID-19 main protease through structure-guided virtual screening approach

Modeling the sterol-binding domain of aster-A provides insight into its multiligand specificity

Virtual screening approach to identify high-affinity inhibitors of serum and glucocorticoid-regulated kinase 1 among bioactive natural products: Combined molecular docking and simulation studies

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

UCSF chimera -A visualization system for exploratory research and analysis

GROMACS 4.5: A high-throughput and highly parallel open source molecular simulation toolkit

Crystal structure and mechanistic determinants of SARS coronavirus nonstructural protein 15 define an endoribonuclease family

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

Cell signalling in insulin secretion: The molecular targets of ATP, cAMP and sulfonylurea

Activity loss by H46A mutation in Mycobacterium tuberculosis isocitrate lyase is due to decrease in structural plasticity and collective motions of the active site

Structure-based discovery of phenyl-diketo acids derivatives as Mycobacterium tuberculosis malate synthase inhibitors

Virtual screening, pharmacokinetics, molecular dynamics and binding free energy analysis for small natural molecules against cyclin-dependent kinase 5 for Alzheimer's disease

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

Elucidating biophysical basis of binding of inhibitors to SARS-CoV-2 main protease by using molecular dynamics simulations and free energy calculations

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

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

Identification of new anti-nCoV drug chemical compounds from Indian spices exploiting SARS-CoV-2 main protease as target

GROMACS: Fast, flexible, and free

Complete genomic sequence of human coronavirus OC43: Molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event

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

What contributes to serotonin-norepinephrine reuptake inhibitors' dual-targeting mechanism? The key role of transmembrane domain 6 in human serotonin and norepinephrine transporters revealed by molecular dynamics simulation

We are thankful to Bioinformatics Resources and Applications Facility (BRAF), CDAC, Pune, for providing access to the computational facility.

No potential conflict of interest was reported by the authors. 

http://orcid.org/0000-0003-0419-0684