key: cord-1034231-aqb9zh6h
authors: Kumar, Prateek; Bhardwaj, Taniya; Kumar, Ankur; Gehi, Bhuvaneshwari R.; Kapuganti, Shivani K.; Garg, Neha; Nath, Gopal; Giri, Rajanish
title: Reprofiling of approved drugs against SARS-CoV-2 main protease: an in-silico study
date: 2020-11-12
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1845976
sha: d3b1f9961df95364f4be816708cd1c8190b2f4a4
doc_id: 1034231
cord_uid: aqb9zh6h

Given the COVID-19 pandemic, currently, there are many drugs in clinical trials against this virus. Among the excellent drug targets of SARS-CoV-2 are its proteases (Nsp3 and Nsp5) that plays vital role in polyprotein processing giving rise to functional nonstructural proteins, essential for viral replication and survival. Nsp5 (also known as M(pro)) hydrolyzes replicase polyprotein (1ab) at eleven different sites. For targeting M(pro), we have employed drug repurposing approach to identify potential inhibitors of SARS-CoV-2 in a shorter time span. Screening of approved drugs through docking reveals Hyaluronic acid and Acarbose among the top hits which are showing strong interactions with catalytic site residues of M(pro). We have also performed docking of drugs Lopinavir, Ribavirin, and Azithromycin on SARS-CoV-2 M(pro). Further, binding of these compounds (Hyaluronic acid, Acarbose, and Lopinavir) is validated by extensive molecular dynamics simulation of 500 ns where these drugs show stable binding with M(pro). We believe that the high-affinity binding of these compounds will help in designing novel strategies for structure-based drug discovery against SARS-CoV-2. Communicated by Ramaswamy H. Sarma

Traditional drug development, despite of technical advancement is an expensive, laborious, and time-consuming process. In the pandemic situation like Coronavirus Disease 2019 , where no proper cure is available, drug repurposing or repositioning of an approved drug is a desirable choice. Repurposing is a progressive strategy that helps in identifying a new use for already approved drugs. With undeniable advantages of shorter time requirements and lower costs for development, use of an approved drug against a novel target can fulfill the principle need in this current pandemic (Ashburn & Thor, 2004; Breckenridge & Jacob, 2019; Pushpakom et al., 2019) . Therefore, drug repurposing against different targets of newly emerged Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) may deliver faster results.

Our previous report on SARS-CoV-2 has elaborately discussed the identified or probable functions and presence of disorder in its proteome . Among the known protein targets of SARS-CoV-2, non structural protein 3 and 5 (also called as M pro or 3CL-protease) proteases are excellent drug targets due to their roles in processing of polyprotein into the constituent nonstructural proteins needed for viral replication and survival (Wu et al., 2020) . M pro cleaves the replicase polyprotein 1ab at eleven different sites and is considered to be the best-characterized enzyme targets for drug development against coronaviruses . Recently, one of the crystal structures of SARS-CoV-2 M pro (PDB ID: 6LU7) in complex with N3 (asynthetic inhibitor) is reported by Zhenming et al. (Anand et al., 2003; Jin et al., 2020) .

Previously, we have established several studies to identify some of the existing drugs as inhibitors of other viral proteomes (Kumar et al., , 2019 Kumar, Sharma et al., 2020; Sharma et al., 2020) . One such report discussed the inhibitory action of hydroxychloroquine against Zika virus protease . Here, in this study, we made an attempt to identify druggable sites in M pro and also investigated the binding interactions between M pro and approved drug molecules. We have screened a library of drugs approved by USA, Canada and European authorities, fetched from the Drugbank database, against active-site of SARS-CoV-2 M pro . As revealed from docking and simulation results, top hit compounds and approved drugs (which are currently in clinical trials) (Cao et al., 2020; 'COVID-19 ring-based prevention trial with lopinavir/ritonavir', 2020; 'Evaluation of the efficacy of the hydroxychloroquine-azithromycin combination in the in the prevention of COVID-19 related SDRA', Khalili et al., 2020;  'Safety and efficacy of hydroxychloroquine associated with azithromycin in SARS-CoV2 virus', 2020) are showing stable and robust binding with active site residues of enzyme M pro . High-affinity binding interactions of these compounds could be a basis to design or improve the novel strategies for further clinical trials and structure-based drug discovery against SARS-CoV-2.

First crystal structure of SARS-CoV-2 main protease (PDB ID: 6LU7) (Jin et al., 2020) in complex with N3 inhibitor is made available quickly for exploring structure-based drug designing opportunities. The structure is prepared and refined using Protein Preparation Wizard of Schrodinger LLC (Madhavi Sastry et al., 2013) . The improper bond parameters and sidechain disarrangement in M pro structure are corrected in Maestro and minimized by OPLS 2005 forcefield. The M pro structure is used to generate active sites from Sitemap program (Halgren, 2009) . A grid is generated at the inhibitor binding site to analyze the interactions of approved drugs with the active site of M pro . A library of 2454 approved drugs obtained from Drugbank database (Wishart et al., 2018) is prepared in LigPrep module by generating at most 4 stereoisomers and tautomers (total 5025 structure in output) per ligand at pH 7 using Epik program.

Docking studies of m pro with approved drugs Prepared library of approved drugs is docked against prepared M pro structure using Glide XP (Extra Precision method) scoring algorithm embedded in Schrodinger. OPLS 2005 forcefield is used for calculating bonded and non-bonded parameters such as van der Waal interactions (Friesner et al., 2006; Halgren et al., 2004) . For the detailed protocol of Glide, see our previously reported articles (Kumar et al., 2019; Kumar, Saumya et al., 2020; Sharma et al., 2017; Yadav et al., 2020) . The interaction of docked compounds with M pro is analyzed in the 4 Å cutoff of active site. Further, binding energy calculation is done using Prime module (based on Molecular mechanics-Generalized Born Surface Area: MM-GBSA) Jacobson et al., 2004) .

To study the structural dynamics and binding stability of inhibitors with M pro at a molecular level, we performed MD simulation in SPC water models using GROMOS54A7 forcefield in Gromacs upto 500 ns (Berendsen et al., 1995) . Charge neutralization is done by adding 3 Na 1 ions and 0.15 M salt in the cubic simulation box. For energy minimization, 50,000 steps of steepest descent algorithm are executed with Verlet cutoff scheme to calculate the neighboring interactions. Further, equilibration process under NPT and NVT conditions is done for 1 ns. Parrinello-Rahman and V-rescale methods are availed for pressure and temperature coupling, respectively. The inhibitor topologies are generated using PRODRG server. Finally, the production runs for all four systems (M pro with Hyaluronic acid (HA), Acarbose, Amikacin, and Lopinavir) is performed in periodic boundary conditions for 500 ns each. LINCS algorithm is used for calculating bond parameters (Hess et al., 1997) . Particle Mesh Ewald for longrange electrostatics with fourier spacing of 0.16 is used for production MD simulations.

The analysis consisting root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg) for C-a atoms, are calculated using gmx, rms, rmsf, and gyrate commands, respectively. The number of hydrogen bonds are calculated in VMD using Hydrogen bonds plugin (Humphrey et al., 1996) . Further, we have used covar and anaeig commands of gromacs to interpret the cumulative motions occurring in protein by virtue of ligand binding describing the principle component analysis (PCA). We have also measured the solvent accessible surface area (SASA) using sasa command in Gromacs to determine the exposed surface area of protein to the solvents.

The M pro crystal structure (PDB ID: 6LU7) is used in this study to screen the approved drug molecules for the purpose of identifying a potential inhibitor candidate. To proceed with the screening process, identification of a druggable site is the primary requirement for a potent outcome. As illustrated in Figure 1 , five different druggable sites are predicted in M pro structure by employing a SiteMap module. Among these five predicted sites, the most appropriate site based on the size and volume of the pocket, hydrophobicity and hydrophilicity of amino acid, site score, and druggability score (Dscore) is selected for docking and simulation analyses (Table 1) . Out of all five sites, site 1 obtains the highest site score and Dscore. Importantly, it consists of a catalytic dyad having -His41 and -Cys145, essential for the protease activity of M pro . Also, a comparison between the predicted site 1 and inhibitor binding site of N3 (co-crystallized N3 inhibitor with M pro ; 6LU7) reveals a similar type of residue distribution in both sites (Table 1) . Therefore, we have considered site 1 for further docking studies. A grid was generated at site 1 with the coordinates of À12.49 Å, 15.06 Å and 74.04 Å for X, Y and Z, respectively.

2454 approved drug molecules are selected from DrugBank database for molecular docking at predicted site 1 of M pro of SARS-CoV-2. These selected drugs are first subjected to the ligand preparation module, Ligprep, to generate at most four stereoisomers for each drug molecule (produced 5025 structures). Resulted structures are then subjected further in molecular docking protocol, which generated 9263 binding poses from 5025 input structures based on binding energies. Based on their interactions and docking scores, we have shortlisted top 100 hits which are tabulated in Table 2 . Most of these listed compounds are antibiotics and antioxidants or have anti-inflammatory as well as anti-cancer properties ( Table 2 ). In addition to top 100 compounds we have also shortlisted Ribavirin, Lopinavir, and Azithromycin (Table 3) for further analysis as these approved drugs are currently being tested on SARS-CoV-2 infected patients. The investigation of stability and binding interactions of these compounds at M pro active site 1 will help in designing novel strategies for therapeutics against COVID-19. Detailed description of docking parameters of top hits and selected compounds are given in Tables 2 and 3 .

Hyaluronic acid (HA): It is an anionic non-sulfated linear glycosaminoglycan polymer of two saccharide units, D-glucuronic acid and N-acetyl-d-glucosamine. HA is stabilized at the active site 1 of M pro via hydrogen bonds (H-bonds) with Ser46, Thr24, Thr25, Thr26, Asn142 and Glu166 residues ( Figure 2 ). HA shows docking score and binding energy equivalent to À13.54 kcal/mol and À49.46 kcal/mol, respectively (Table 2) . It is majorly found in connective tissues, umbilical cord, vitreous fluid in eyes, and synovial fluid of joints (Fraser et al., 1997) . Due to its presence in the extracellular matrix, HA is biocompatible and biodegradable, and therefore, functions as a drug delivery carrier (Widjaja et al., 2014) . Using cell culture studies, HA has been demonstrated to work as an antiviral compound against both DNA and RNA viruses such as Coxsackievirus B5, Herpes Simplex Virus Residue distribution in predicted site 1 and experimental N3 binding site is observed to be very similar. Residues of active site in crystal structure and predicted active site, in bold fonts, show similarity with each other. (Cermelli et al., 2011) .

Acarbose: Also known as acarbosa and acarbosum, it is an approved oral drug prescribed to Diabetes Mellitus type II patients. Acarbose is observed to interact with Thr25, Thr26, Ser46, Gly143, Glu166 and Gln189 residues of M pro (Figure 3) . These residues are a part of domains I and II which constitute b-barrels enclosing the active site. It binds to Mpro with a significant docking score of À12.181 kcal/mol and binding energy of À45.54 kcal/mol. It reduces the amount of glucose in blood by delaying the absorption of carbohydrates.

Acarbose is a complex oligosaccharide that works by inhibiting the a-amylase and glycosidase enzymes of pancreatic and intestinal cells, thereby preventing the hydrolysis of poly-carbohydrates to glucose (DiNicolantonio et al., 2015) .

Amikacin: According to our results, docking score of drug Amikacin is À10.966 kcal/mol, and binding energy of À52.25 kcal/ mol. Drug-interaction view exposes the interaction of Amikacin with His41 residue of M pro which is one of the two residues of catalytic dyad conserved in catalytics site of coronaviruses ( Figure  4 ). It also interacts with Phe140, Asn142, Gly143 and Glu166 residues present in domain II. It is a well-known broad-spectrum antibacterial drug prescribed for the treatment of life-threatening gram-negative bacterial infections (Ramirez & Tolmasky, 2017; Tamma et al., 2012) . Recently in 2018, Amikacin liposome inhalation suspension (ALIS, ArikayceV R ) is approved for treating patients suffering from mycobacterium avium complex lung diseases therapy (Shirley, 2019) .

Ribostamycin: Docking score of Ribostamycin on SARS-CoV-2 protease is estimated to be À11.845 kcal/mol with a glide binding energy of À36.45 kcal/mol ( Table 2 ). The M pro -Ribostamycin complex is stabilized by five H-bonds formed between Phe140, Leu141, His164 and Glu166 residues of protease and the compound (Figure 5(a) ). It is an aminoglycosidic anti-bacterial compound but is demonstrated to have an antiviral activity against HIV-1 as well (Ennifar et al., 2006) . Our previous in-silico study has also reported nsp2 cysteine protease of chikungunya virus as one of the targets of Ribostamycin (Kumar et al., 2019) .

Octreotide: The long-acting octapeptide has a cyclic structure mimicking to that of natural hormone somatostatin. Due to similar pharmacologic properties to endogenous hormone, it is a potent inhibitor of growth hormone, luteinizing hormone, insulin and glucagon (Battershill & Clissold, 1989; McKeage et al., 2003) . Our docking analysis has given it a very significant binding score of À11.162 kcal/mol. The drug form H-bonds with Thr26, Ser46, Gly143, Leu141, Glu166 and Gln189 residues and a salt bridge interaction with Glu166 of protease with a binding energy of À66.41 kcal/mol ( Figure  5 (b)). It is recommended for the treatment of acromegaly and thyrotrophinomas (Battershill & Clissold, 1989; McKeage et al., 2003) . As an analogue of somatostatin, Octreotide has been reported to bind with somatostatin receptor and exert its influence on downstream signaling (Hofland & Lamberts, 1996) . In a direct way, this drug has not been reported against any of the viral infection, as per our knowledge. However, several octapeptide substrates have been functionally characterized against Nsp3 helicase protein of all serotypes of Dengue virus (Li et al., 2005) .

Paromomycin: It is prescribed as a first line drug for the treatment of intestinal amebiasis and visceral leishmaniasis (Davidson et al., 2009; Shalev-Benami et al., 2017) . In our docking study, it is observed to interact via H-bonds with Phe140, Leu141, Asn142, Glu166, Gln189 and Thr190 residues of M pro and shows a powerful docking score of À11.159 kcal/ mol ( Figure 5(c) ). Furthermore, the stability of complex is increased by a salt bridge interaction between NH3 1 and Glu166 residue (Details of binding parameters are given in Table 2 ).

Lopinavir: According to our docking results (docking score of À6.119 kcal/mol), Lopinavir-protease complex is stabilized by two H-bonds formed by Thr25, and Glu166 amino acids of protease with the drug (Table 3 ). In conjunction with Hbonds, a Pi-Pi stacking interaction with His41 further increases the stability of complex upto À46.98 kcal/mol of binding energy (Figure 6(a) ). It is an FDA-approved peptidomimetic compound containing a hydroxyethylene scaffold that mimics the peptide bond. Co-formulated capsules with Ritonavir, marketed as Kaletra V R are widely used for the treatment of HIV (Chandwani & Shuter, 2008) . As an anti-retroviral drug, Lopinavir in combination with Ritonavir is demonstrated to suppress viral replication.

Azithromycin: It shows a docking score of À5.636 kcal/mol and is observed to form a single H-bond with Leu141 in addition with a salt bridge interaction between NH 1 and Glu166 residue (Figure 6(b) ). The M pro -Azithromycin complex exhibits a binding energy equivalent to À46.25 kcal/mol. Approved by FDA in 1991, it is a broad-spectrum antibiotic drug administered orally during respiratory, enteric, and genitourinary infections (Peters et al., 1992) . As a macrolide having a 15-membered ring, it shows bacteriostatic activity against Gram-negative as well as Gram-positive bacteria. Although considered safe for use during pregnancy, it occasionally causes diarrhea in breast-feeding infants (McMullan & Mostaghim, 2015; Peters et al., 1992) . As mentioned above, presently, Azithromycin is in clinical trials along with Hydroxychloroquine for the treatment of COVID-19 infection ('Evaluation of the efficacy of the hydroxychloroquine-azithromycin combination in the in the prevention of COVID-19 related SDRA', 2020; 'Safety and efficacy of hydroxychloroquine associated with azithromycin in SARS-CoV2 virus', 2020). However, FDA has approved it for emergency use against SARS-CoV-2.

Ribavirin: In our study, docking analysis of Ribavirin revealed a docking score equivalent to À6.813 kcal/mol (Table 3 ). The drug establishes H-bonds with His164, Glu166, Gln189 and Thr190 residues of SARS-CoV-2 main protease (binding energy À35.63 kcal/mol) (Figure 6(c) ). It is a broad-spectrum FDA-approved antiviral drug used against a number of DNA and RNA viruses including respiratory syncytial virus, parainfluenza virus and hepatitis C virus (Krilov, 2001; Reddy et al., 2009) . It is a synthetic guanosine analogue which interfers with viral RNA synthesis. Administered orally, Ribavirin is prescribed in combination with pegylated interferon-a2 for treating hepatitis C virus infection (Loustaud-Ratti et al., 2016; Reddy et al., 2009 ). On the basis of its in-vitro efficacy against SARS-CoV and MERS-CoV, it is being clinically tested on COVID-19 patients (Khalili et al., 2020) . The catalytic site of M pro and nearby residues in contact such as Phe140, Leu141, His164, Glu166, Glu189 are observed to be important for binding of these drugs. Therefore out of 100 shortlisted drugs, we have chosen HA, Acarbose, and Amikacin on the basis of their high docking scores and molecular interactions with catalytic site and nearby residues of M pro . Additionally, binding stability of the drug Lopinavir is also investigated through molecular simulation.

The stability of protein-ligand complexes with respect to protein conformation is determined by MD simulations. In general, a small deviation in protein conformation during the course of simulation indicates a stable protein structure. We analyzed the atomic distance and fluctuation in protein upon binding of drug molecule at the active site 1 of M pro . Acarbose, HA, Amikacin and Lopinavir complexes with M pro are monitored for 500 ns using our in-house high-performance cluster facility.

Assessment of root mean square deviation (RMSD), root mean square fluctuation (RMSF) and radius of gyration (Rg)

The average RMSD of 6LU7-HA, 6LU7-Acarbose, 6LU7-Amikacin and 6LU7-Lopinavir complex is observed to be 0.42 nm, 0.31 nm, 0.37 nm and $0.4 nm, respectively. 6LU7-HA complex attains a stable conformation till 15 ns, after which deviation starts increasing and the complex become unstable near 400 ns. For 6LU7-Acarbose, RMSD keeps on fluctuating between 0.2 nm and 0.3 nm during the entire simulation period. Nonetheless, these RMSD values are still much lower than that of other two compounds 6LU7-HA and 6LU7-Amikacin (Figure 7) . The RMSF plot identifies flexible regions in protein-ligand complexes. Residues with low RMSF values indicate its contribution to the structured regions, while higher residual fluctuation indicates the presence of unstructured regions such as loops and turns in protein. The binding of an inhibitor or druggable compound through one or more stable bonds, however, reduces the fluctuation in residues. As illustrated in Figure 8 , Lopinavir bound protease complex is observed to have higher fluctuating (0.4-0.6 nm) peaks showing less stable binding. Contrary to this, 6LU7-Acarbose, 6LU7-HA and 6LU7-Amikacin complexes are showing less fluctuations (between 0.1 to 0.3 nm). Similarly, Rg for 6LU7-HA and 6LU7-Acarbose are observed to be 2.18 nm and 2.17 nm, respectively. For 6LU7-Amikacin, and 6LU7-Lopinavir complexes, respective mean Rg values are calculated to be 2.2 nm and 2.16 nm. Graphs in Figure 9 illustrate the Rg for complex of M pro (6LU7) with HA, Amikacin, Acarbose, and Lopinavir. The time-dependent trace of Rg is stable for HA, and Acarbose throughout the simulation. In case of Amikacin, the complex is showing higher fluctuation between 2.15 nm to 2.25 nm beyond 275 ns. Although, 6LU7-Lopinavir, at initial stage of simulation acquires a high Rg upto 2.25 nm, however, beyond 60 ns itdecreases and becomes stable till 425 ns.

Binding of ligand induces local as well as global fluctutaions in protein structure that are hard to distinguish. In order to identify the simultaneously occurring conformational changes, PCA is performed for C-a atoms of protease. It is a multivariate statistical technique which builds a co-varinace or correlation matrix for methodically calculating eigenvectors and eigenvalues. PCA generates eigenvectors which represents the reduced collective atomic motions in protein structure (David & Jacobs, 2014) .

The effect of ligand stabilization on protein structure is shown in Figure 10 where PC1 and PC2 majorly contributes to global conformations of 6LU7-ligand complexes. As observed in the analysis, dense clusters are formed in all drug bound forms of protease but Amikacin and HA have higher range from À10 to 10 nm and À12 to 6 nm, respectively. However, Lopinavir has shown stable clusters with a range of À5 to 5 nm and similarly, the second top hit drug, Acarbose is observed to form clusters of stable positions (À7 to 2 nm). This indicates that binding of Acarbose has reduced the global conformational changes occurring in M pro while HA has slightly higher values, may be due to little fluctuation in the last 100 ns period of simulation.

Hydrophobic force deriving protein folding results in formation of hydrophobic cores with hydrophilic residues assembling near the surface of a protein structure. Binding of ligand can significantly affect the SASA of a protein, in which overall differences related to ligand-induced stability in complex are frequently estimated. Therefore, we also investigated the SASA parameter in this study. Figure 11 depicts the SASA values of C-a atoms of M pro complexes upon binding with screened drugs. In correlation with RMSD, RMSF and Rg analysis, 6LU7-Lopinavir has a gradually increasing surface area (188 nm 2 ) till 200 ns illustrating higher protein surface exposure. The complex is then observed to be fluctuating a little which again decreases after 400 ns to 178 nm 2 . However, Acarbose and HA has a stable trend of SASA values than Amikacin bound M pro with an average value of 184 nm 2 but Amikacin experiences higher fluctuations after 250 ns. This also confirms that the top screened drugs Acarbose and HA have gained firm binding with M pro .

Numerous attempts of drug repurposing have been made against SARS-CoV-2 M pro enzyme. In one of the reports, Sang et al. explored inhibitory effect of approved anti-HIV drugs, which have shown stable interactions with residues Met49, Leu141, Cys145, Met165, Pro168 and Glu189 of M pro (Sang et al., 2020) . In another study, synergistic effect of Lopinavir, Oseltamivir and Ritonavir are investigated. These drugs are stabilized by interaction with Lue287, Tyr101, Asp33, His41, Asn142 and Glu166 residues of M pro (6LU7) (Muralidharan et al., 2020) . Khan et al. showed stable binding of two compounds Paritaprevir and Raltegravir with SARS-CoV-2 M pro (shows multiple interactions with Thr24, Ser46, Leu50, Asn142, Cys145 and Pro168) (Khan et al., 2020) . Another approach of using combination of drugs is frequently observed to be highly effective. Lopinavir-Ritonavir complex shows a stable binding affinity of À10.6 kcal/mol with SARS-CoV-2 M pro . The combination drug is found to interact with one of the catalytic dyad residue His41, and further His164, Glu166, Arg188 and Gln192 of M pro active site .

Our simulation analysis of 6LU7-HA, 6LU7-Acarbose, 6LU7-Amikacin and 6LU7-Lopinavir complex reveals stable binding of Acarbose, and HA at M pro active site. Interaction of Lopinavir stabilized by multiple non-covalent bonds with catalytic as well as surrounding residues (such as Asn142 and Glu166) of M pro is in line with the previously reported literature Muralidharan et al., 2020) . Antiviral effects of HA against both DNA and RNA viruses are already well described (Cermelli et al., 2011) . Although Amikacin and Acarbose are not much known to exert antiviral effects, however, they have shown significant docking and simulation results. Therefore, these drugs may have inhibitory potential and are need to be tested experimentally for their efficacy against SARS-CoV-2. A study by Victorovich et al., has observed that among more than 250 nucleotide mutations (known yet), there are more number of synonymous Cytosine to Uracil transitions in ORF 1a and 1b which have caused mutational U-pressure in RNA plus strand of SARS-CoV-2 and that can not be repaired by proofread machinery of coronaviruses (Victorovich et al., 2020) .

The current scenario of the COVID-19 pandemic urgently needs drugs to treat patients. However, at the same time, we need to understand the mechanisms or the interactions of those drugs on the molecular targets on SARS-CoV-2. In this study, we have screened the approved drugs against the M pro of SARS-CoV-2 to investigate their binding efficacies and molecular interactions. We found that HA, Acarbose, Ribostamycin, Paromomycin, and Amikacin have good binding efficacy with the target. The interaction with the catalytic dyad (-His41, -Cys145) at the active site, seen in Amikacin, and Lopinavir may be responsible for inhibiting M pro . Additionally, HA and Acarbose have also given the stable binding in 500 ns MD simulation time as observed from atomic distance, fluctuations and PCA. Further, the inhibition of M pro may lead to the inhibition of SARS-CoV-2 replication and propagation. As elucidated, these drugs may act as potential inhibitors against COVID-19 by targeting M pro of SARS-CoV-2. However, in-vitro validation of these identified drugs is essential to test the inhibition potency.

RG, NG and GN: Conception and design. RG: review and writing of the manuscript, and study supervision. PK, TB, AK, BRG, SK: acquisition andinterpretation of data, writing of the manuscript.

E7 oncoprotein of human papillomavirus: Structural dynamics and inhibitor screening study

Coronavirus main proteinase (3CLpro) structure: Basis for design of anti-SARS drugs

Drug repositioning: Identifying and developing new uses for existing drugs

Octreotide. A review of its pharmacodynamic and pharmacokinetic properties, and therapeutic potential in conditions associated with excessive peptide secretion

GROMACS: A message-passing parallel molecular dynamics implementation

Overcoming the legal and regulatory barriers to drug repurposing

A trial of lopinavir-ritonavir in adults hospitalized with severe Covid-19

In vitro evaluation of antiviral and virucidal activity of a high molecular weight hyaluronic acid

Lopinavir/ritonavir in the treatment of HIV-1 infection: A review

COVID-19 ring-based prevention trial with lopinavir/ritonavir -Full text view -ClinicalTrials.gov

Principal component analysis: A method for determining the essential dynamics of proteins BT

Acarbose: Safe and effective for lowering postprandial hyperglycaemia and improving cardiovascular outcomes

Targeting the dimerization initiation site of HIV-1 RNA with aminoglycosides: From crystal to cell

Hyaluronan: Its nature, distribution, functions and turnover

Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes

When darkness becomes a ray of light in the dark times: Understanding the COVID-19 via the comparative analysis of the dark proteomes of SARS-CoV-2, human SARS and bat SARS-like coronaviruses

Identifying and characterizing binding sites and assessing druggability

Glide: A new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening

LINCS: A linear constraint solver for molecular simulations

Somatostatin receptors and disease: Role of receptor subtypes. Bailli ere's

VMD: Visual molecular dynamics

A hierarchical approach to all-atom protein loop prediction

Structure of Mpro from COVID-19 virus and discovery of its inhibitors

Novel coronavirus treatment with ribavirin: Groundwork for an evaluation concerning COVID-19

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

Respiratory syncytial virus: Update on infection, treatment, and prevention

Hydroxychloroquine inhibits Zika virus NS2B-NS3 protease

Targeting the NTPase site of Zika virus NS3 helicase for inhibitor discovery

Mechanistic insights into Zika virus NS3 helicase inhibition by epigallocatechin-3-gallate

Targeting the nsp2 cysteine protease of chikungunya virus using FDA approved library and selected cysteine protease inhibitors

Identification of peptidomimetic compounds as potential inhibitors against MurA enzyme of Mycobacterium tuberculosis

In silico prediction of potential inhibitors for the Main protease of SARS-CoV-2 using molecular docking and dynamics simulation based drug-repurposing

Functional profiling of recombinant NS3 proteases from all four serotypes of dengue virus using tetrapeptide and octapeptide substrate libraries

Ribavirin: Past, present and future

Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments

Octreotide long-acting release (LAR): A review of its use in the management of acromegaly

Prescribing azithromycin

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

Drug repurposing: Progress, challenges and recommendations

Amikacin: Uses, resistance, and prospects for inhibition

Safety and efficacy of hydroxychloroquine associated with azithromycin in SARS-CoV2 virus (Coalition Covid-19 Brasil II) -Full text view -ClinicalTrials.gov

Anti-HIV drug repurposing against SARS-CoV-2

Atomic resolution snapshot of Leishmania ribosome inhibition by the aminoglycoside paromomycin

Polysaccharides like pentagalloylglucose, parishin A and stevioside inhibits the viral entry by binding the Zika virus envelope protein

Epigallocatechin gallate, an active green tea compound inhibits the Zika virus entry into host cells via binding the envelope protein

Amikacin liposome inhalation suspension: A review in Mycobacterium avium complex lung disease

Combination therapy for treatment of infections with gram-negative bacteria

Evaluation of the efficacy of the hydroxychloroquine-azithromycin combination in the in the prevention of COVID-19 related SDRA -Full Text View -ClinicalTrials.gov

Translation-associated mutational Upressure in the first ORF of SARS-CoV-2 and other coronaviruses

Hyaluronic acid-based nanocomposite hydrogels for ocular drug delivery applications

DrugBank 5.0: A major update to the DrugBank database

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

Investigating into the molecular interactions of flavonoids targeting NS2B-NS3 protease from Zika virus through in-silico approaches

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved a-ketoamide inhibitors

All the authors would like to thank IIT Mandi for providing HPC facilities. RG would like to thank the Department of Biotechnology, Govt of India (BT/11/IYBA/2018/06). TB is grateful to the Department of Science and Technology for INSPIRE fellowship. PK, SK, AK, and BRG are supported by MHRD fellowship, Govt. of India.

All authors declare that there is no financial or any other type of competing interest.

This work was supported by Department of Biotechnology, Ministry of Science and Technology.

http://orcid.org/0000-0001-7392-5045Taniya Bhardwaj http://orcid.org/0000-0002-7226-6348 Neha Garg http://orcid.org/0000-0003-2227-8292 Rajanish Giri http://orcid.org/0000-0002-2046-836X