key: cord-263576-pn2zieek
authors: Das, Sourav; Sarmah, Sharat; Lyndem, Sona; Singha Roy, Atanu
title: An investigation into the identification of potential inhibitors of SARS-CoV-2 main protease using molecular docking study
date: 2020-05-13
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1763201
sha: 
doc_id: 263576
cord_uid: pn2zieek

A new strain of a novel infectious disease affecting millions of people, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has recently been declared as a pandemic by the World Health Organization (WHO). Currently, several clinical trials are underway to identify specific drugs for the treatment of this novel virus. The inhibition of the SARS-CoV-2 main protease is necessary for the blockage of the viral replication. Here, in this study, we have utilized a blind molecular docking approach to identify the possible inhibitors of the SARS-CoV-2 main protease, by screening a total of 33 molecules which includes natural products, anti-virals, anti-fungals, anti-nematodes and anti-protozoals. All the studied molecules could bind to the active site of the SARS-CoV-2 protease (PDB: 6Y84), out of which rutin (a natural compound) has the highest inhibitor efficiency among the 33 molecules studied, followed by ritonavir (control drug), emetine (anti-protozoal), hesperidin (a natural compound), lopinavir (control drug) and indinavir (anti-viral drug). All the molecules, studied out here could bind near the crucial catalytic residues, HIS41 and CYS145 of the main protease, and the molecules were surrounded by other active site residues like MET49, GLY143, HIS163, HIS164, GLU166, PRO168, and GLN189. As this study is based on molecular docking, hence being particular about the results obtained, requires extensive wet-lab experimentation and clinical trials under in vitro as well as in vivo conditions. Communicated by Ramaswamy H. Sarma

The recent outbreak of a novel coronavirus disease , occurring from a severe acute respiratory syndrome (SARS) like coronavirus started in Wuhan, China, is spreading rapidly in humans, which is now considered as a global pandemic . WHO, on 30 January 2020, declared this outbreak as Public Health Emergency of International Concern and later on 11 February 2020 named it as ) events as they happen," 2020). According to the data of WHO (https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports), currently (as on 24 th April 2020), more than 2.6 million COVID-19 cases have been confirmed globally, and the number of death toll rose to 1,81,938 which are increasing exponentially. To date no specific medications have been developed, and therefore considering the risk factors associated with this disease, there is an urgent need for a treatment method to treat this disease in order to limit the transmission (Enayatkhani et al., 2020) . Among all the treatment strategies, the uses of already known viral medications have great advantages as the pharmacokinetics, pharmacodynamics, and safety profiles of these drugs have already been established. This drug repurposing strategy has gained significance as they are expected to be faster with less investment (Fitchett et al., 2014; Law et al., 2013; Pant et al., 2020; Park, 2019) . Preliminary studies revealed Lopinavir/ritonavir combined therapy as a potential inhibitor of the virus. Along with these two drugs, many other antiviral drugs have also been screened (Lu, 2020; Xu et al., 2020) . Recently, anti-malarial drug chloroquine (M. and hydroxychloroquine (Gautret et al., 2020) have been reported to have certain curative effect on COVID-19. Due to diverse bioactivity and less toxicity, naturally occurring bioactive flavonoid molecules have also been screened for their potential therapeutic effects (Khaerunnisa et al., 2020; Thuy et al., 2020) . Numerous reports to identify the potential inhibitors of novel coronavirus have recently been published using computational methods (Aanouz et al., 2020; Elmezayen et al., 2020; Gupta et al., 2020; Muralidharan et al., 2020; Sarma et al., 2020) . The transcription of the genomic RNA inside the cell generates a polypeptide which is proteolytically cleaved by the main protease (M Pro which is also known as 3CL Pro or chymotrypsin-like protease) at 11 different sites to produce proteins necessary for the replication of the virus. Inhibition of the activity of the main protease is expected to block the viral replication as no human proteases are known having similar cleavage specificity (Boopathi et al., 2020; Zhang et al., 2020) . Therefore M Pro is one of the most explored drug target among the corona viruses that plays a major role in the translation and replication of new generation viruses from the viral genomic RNA (Morse et al., 2020; Wu et al., 2006) . Many crystal structures of M Pro of SARS-CoV-2 have already been submitted to protein data bank (PDB) which are in both liganded and unbounded form (Jin et al., 2020; Owen et al., 2020; Zhang et al., 2020) . The 3 D structure of SARS-CoV-2 M pro (PDB ID. 6Y84), is made up of three domains, where domain I, II and III spans across 8-101, 102-184 and 201-306 amino acid residues, repectively. Domain I and II mainly consists of b-barrels while domain III mainly consists of a-helices (S. A. . The substrate/active site of binding is located at the cleft of domain I and II, which consists of two catalytic residues, namely, HIS41 and CYS145 ( Figure 1 ).

Here in this report, a total of 33 compounds belonging to different classes (i) control drugs, (ii) natural products, (iii) anti-virals, (iv) anti-fungals, (v) anti-nematodes and antiprotozoal were identified to investigate their inhibitory effects towards M Pro by molecular docking analysis. Ritonavir, hydroxychloroquine, penciclovir and lopinavir are known drugs that have been screened for the treatment of COVID-19. The reason for them being used as control drugs in the study is that, as mentioned earlier, lopinavir/ritonavir combination, and anti-malarial drug hydroxychloroquine has recently been found to be a potential inhibitor of the virus, while penciclovir could also inhibit the virus to some extent (M. . Considering the diverse pharmaceutical properties of the natural products including anti-viral properties (Jo et al., 2020; Thayil & Thyagarajan, 2016; Zakaryan et al., 2017) , a total of 17 compounds: curcumin, demethoxycurcumin, (-)-Epigallocatechin gallate (EGCG), (-)-epigallocatechin, (EGC), hesperidin, myricitrin, puerarin, scutellarin, rutin, quercitrin, capsaicin, ursolic acid, glabridin, apiin, rhoifolin, glycyrrhizin, and vitexin were screened for their inhibitory potential against SARS-CoV-2 M pro . Azidothymidine, indinavir, tipranavir, and saquinavir are the well-known anti-viral drugs (De Clercq & Li, 2016) , while diethylcarbamazine, primaquine, mepacrine, artemisinin and niclosamide are the anti-nematodes (Panic et al., 2014) and emetine, an anti-protozoal drug (Akinboye & Bakare, 2011) , which have also been screened for their inhibitory potentials. Further, two anti-fungals, namely, fluconazole and itraconazole (Vardakas et al., 2005) were also docked to the main protease of COVID-19. This study is focused on the repurposing/identification of compounds with a sole aim to expedite the identification of specific drugs for the treatment of COVID-19.

Molecular docking methods are the best tools to predict the drug interactions with macromolecules. The blind docking method involves a search throughout the entire surface of the macromolecule for binding sites. Therefore, blind molecular docking analyses of some known drugs along with some bioactive natural compounds, have been screened with the main protease of COVID-19, SARS-CoV-2 M pro . 

and niclosamide (4322) and emetine (9802) were obtained from ChemSpider (www.chemspider.com) as a .mol file and further geometry optimization was carried out in ArgusLab using PM3 methods (Thompson, 2004) . The optimized structures of the ligands are shown in Figure S1 . The code within the parenthesis signifies the ChemSpider ID of the respective molecules.

The drug likeliness of a molecule is indicated by the Lipinski's rule of five parameters (molecular weight <500 Da, no more than 5 hydrogen bond donors, no. of hydrogen bond acceptors should be less than 10 and logP should not be greater than 5). The Lipinski's rule of five parameters were obtained from the SWISSADME server (www.swissadme. ch/index.php) (Daina et al., 2017) . The chemical structures, chemical formula and the Lipinski's rule parameters of the ligands are listed in Table S1 (Supplementry Information).

The crystal structure of M Pro (PDB ID: 6Y84) was obtained from Protein Data Bank (PDB) (Owen et al., 2020) . The obtained structure was a dimer of two homologous amino acid chains (A and B). The chain A was used for molecular docking. Water molecules were removed from the PBD using PyMOL (Schrodinger LLC, 2017 ).

The final PDB file of the protease and geometry optimized ligand structures were selected for molecular docking analyses in the SwissDock server (http://www.swissdock.ch/docking). SwissDock is based on EADock DSS docking software which estimates the CHARMM (Chemistry at HARvard Macromolecular Mechanics) energy on a grid and the most favourable binding modes are clustered in the server (Grosdidier et al., 2011) . The docked pose with the minimum fullfitness score is considered for futher analysis. UCFS Chimera (Pettersen et al., 2004) , a molecular visualization tool, was used to visualize the results obtained from the server. Discovery Studio Visualizer was further employed in order to prepare the docked poses and 2 D interaction plots (Biovia, 2019) . The changes in the accessible surface area of the SARS-CoV-2 M pro protease upon ligand binding were determined using an online server (http://cib.cf.ocha.ac.jp/bitool/ASA/ ). The electrostatic surface potential was determined using the APBS plugin of PyMOL (Schrodinger LLC, 2017) .

As mentioned earlier, in addition to the control drugs for the treatment of COVID-19, here we have performed the molecular docking studies of various naturally occurring compounds, anti-virals, anti-fungals, anti-nematodes and anti-protozoal drug with SARS-CoV-2 M pro . According to the SwissDock server, the binding modes of the ligands are ranked according to their fullfitness score, therefore the conformers having the minimum fullfitness score were analyzed for further evaluations. Table 1 list the fullfitness score and the free binding energy associated with each corresponding docked pose, this Table 1 will be further utilized in the following sections which rationalize the results obtained from the docking studies one by one.

The docked pose of the minimum energy (fullfitness score) conformers of the four control drugs, namely, penciclovir, ritonavir, hydroxychloroquine and lopinavir along with their corresponding 2 D interaction plots are depicted in Figure 2 .

The docked poses clearly demonstrate that the drugs molecules bind within the active site of the SARS-CoV-2 M pro macromolecular structure. Figure 2a shows that penciclovir binds firmly through three conventional hydrogen bonds with residues HIS41, CYS44 and GLU166 besides other interactions such as carbon-hydrogen bonds (ASN142 and MET165) and van der Waals interactions. Ritonavir binds to the active site of SARS-CoV-2 with different types of noncovalent interactions, as depicted in Figure 2b . Ritonavir stabilizes in the active site through hydrogen bonding with THR25 and SER46, p-sigma with THR26, alkyl hydrophobic with MET165, p-alkyl hydrophobic with LEU27, CYS145 and PRO168, and van der Waals interactions with other residues as depicted in the 2 D plot. Hydroxychloroquine, a promising candidate for the treatment of the current pandemic due to SARS-CoV-2 (Gautret et al., 2020) , has been found to bind within the active site of the protease through p-sulphur interaction with MET165, p-sigma with GLN189, alkyl hydrophobic with LEU167 and PRO168, and van der Waals interactions with other residues as shown in Figure 2c . Similarly, lopinavir interacts with the residues in the active site through hydrogen bonding with THR26, p-sulphur interaction with CYS145 and alkyl hydrophobic with PRO168, besides van der Waals interaction with other residues as depicted in the 2 D plot (Figure 2d ). The fullfitness score along with the estimated free energy (DG) of binding of the above mentioned drugs are listed in Table 1 . If we compare in terms of the DG value, we can see that ritonavir has the most negative DG value with -9.52 kcal/mol followed by lopinavir, then hydroxychloroquine and finally penciclovir. This can be directly correlated with the number of non-covalent interactions that these drugs undergo with the surrounding residues within the active site of SARS-CoV-2. Moreover, the stability of a particular drug within the active site is also associated with the number of p-interactions that it undergoes with the surrounding residues (Arthur & Uzairu, 2019) . The electrostatic surface potential of the binding site along with the simultaneous presence of the four drugs are shown in Figure S2 .

The docked pose of the minimum energy (fullfitness score) conformers of the 17 natural products, namely, curcumin, demethoxycurcumin, EGCG, EGC, hesperidin, myricitrin, puerarin, scutellarin, rutin, quercitrin, capsaicin, ursolic acid, glabridin, apiin, rhoifolin, glycyrrhizin, vitexin along with their corresponding 2 D interaction plots are depicted in Figure 3 . Curcumin, a potent bioactive molecule binds in the active site of SARS-CoV-2 M pro (Figure 3a ) through hydrogen bonding with GLY143 and GLN192, p-sulphur, p-sigma interactions with CYS145 and PRO168, respectively, along with other non-covalent interactions such as van der Waals interactions with other residues as shown in the 2 D plot (Figure Table 1 . Fullfitness score and estimated change in free energy of the minimum docked pose of the respective possible inhibitors with SARS-CoV-2 M pro . 3a). EGC, a tea polyphenol, binds to the active site through hydrogen bonding with THR26, HIS41 and ASN142, and van der Waals forces with other residues (Figure 3b ). Demethoxycurcumin binds in the active site (Figure 3c ) through hydrogen bonding with CYS44 and PRO168, p-sigma with PRO168 and p-alkyl with MET49. Hesperidin interacts through hydrogen bonding with THR24, THR25, THR45, HIS4, SER46, CYS145, amide-p stacked interaction with THR45, p-alkyl interactions with MET49 and CYS145 (Figure 3d ). EGCG interacts with PHE140, GLU166, GLN192 through hydrogen bonding, CYS145 through p-sulphur, GLN189 through p-sigma and other non-covalent forces as depicted in the 2 D plot of Figure 3e . Myricitrin stabilizes in the active site mainly through hydrogen bonding with the THR24, THR25, THR26, ASN119, ASN42 residues as shown in Figure  3f , while puerarin gets stabilized by hydrogen bonding with HIS41, CYS44, GLY143 and GLU166 residues (Figure 3g) . Quercitrin, on the other hand, gets stabilized by hydrogen bonding with THR25, GLY143 and GLU166, and amide-p stacking with THr45 along with other interactions as depicted in Figure 3h . Scutellarin (Figure 3i ) forms hydrogen bonds with THR26, GLY143 and CYS145 along with a p-sulphur interaction with CYS145. Capsaicin forms hydrogen bonds with THR190 along with alkyl hydrophobic (with CYS145) and p-alkyl interactions with HIS163 and PRO168 (Figure 3j ). Rutin undergoes several non-covalent interactions with the residues within the active site (Figure 3k ), it gets stabilized through hydrogen bonding with HIS41, LEU141, ASN142, GLU166, THR190 and GLN192, further it undergoes p-sulphur interaction with CYS145 and p-alkyl with PRO168. Ursolic acid undergoes only van der Waals interactions with the surrounding residues as shown in Figure 3l . Vitexin forms hydrogen bonds with THR26, THR45 and GLY143, and p-sigma interaction with ASN142 ( Figure 3m ). Glabridin interacts through hydrogen bonding with GLY143, alkyl hydrophobic with LEU27 and CYS145, and p-alkyl with HIS41 and CYS145 ( Figure 3n ). Glycyrrhizin forms a hydrogen bond with THR190, alkyl hydrophobic with MET49, CYS145 and MET165, and p-alkyl with HIS41 ( Figure 3o ). Rhoifolin stabilizes in the active site ( Figure 3p ) mainly through hydrogen bonding with HIS41, CYS44, ASN119 and GLU166, alkyl, and p-alkyl hydrophobic interactions with CYS145 and HIS41, respectively. Apiin undergoes several non-covalent interactions at the active site, it forms hydrogen bonds with HIS41, SER46, GLN189, p-sigma interaction with PRO168, and p-alkyl interaction with MET165 ( Figure 3q ). In all the above docked pose, van der Waals play a major role in the binding processes as could be observed from the respective 2 D plots from Figure 3 . On comparison of the DG value, which gives the estimated free energy of binding, we can infer that the binding affinity of the natural products towards SARS-CoV-2 protease follows the following order rutin (À9.55 kcal/mol).> hesperidin (À9.02 kcal/mol) > apiin (À8.74 kcal/mol) > puerarin (À8.43 kcal/mol) > Rhoifolin (À8.37 kcal/mol) > scutellarin (À8.32 kcal/mol) > quercitrin (À8.26 kcal/mol) > curcumin (8.15 kcal/mol) > capsaicin (À8.12 kcal/mol) > EGCG (7.96 kcal/ mol) > demethoxycurcumin (À7.87 kcal/mol) $ myricitrin (À7.85 kcal/mol) $ ursolic acid (À7.84 kcal/mol) $ vitexin (À7.82 kcal/ mol) $ glycyrrhizin (À7.81 kcal/mol) > glabridin (À7.59 kcal/ mol) > EGCG (À7.30 kcal/mol). The electrostatic surface potential of the binding site along with the simultaneous presence of the 17 natural products is shown in Figure S3 .

The docked pose of the minimum energy (fullfitness score) conformers of the two anti-fungals, namely, fluconazole and itraconazole along with their corresponding 2 D interaction plots are depicted in Figure 4 . It can be seen that the two ligands bind within the active site of SARS-CoV-2. Fluconazole interacts through hydrogen bonding with GLY163 and GLU166, p-alkyl hydrophobic interactions with CYS145 and MET165, halogen (fluorine) bonding with CYS144 and van der Waals forces with other residues as depicted in the 2 D plot (Figure 4a ). Another anti-fungal drug itraconazole is found to form hydrogen bonds with SER46 and GLY143, p-sulphur interaction with CYS145, alkyl hydrophobic and p-alkyl interactions with PRO168, along with van der Waals interactions as shown in the 2 D plot (Figure 4b ). it can be seen from Table 1 , that the DG value for itraconazole is more negative (À9.05 kcal/mol) as compared to that of fluconazole binding to SARS-CoV-2, which can also be seen from the number of interactions that each molecule undergoes with the active site residues. The electrostatic surface potential of the binding site along with the simultaneous presence of the two drugs is shown in Figure S4 .

The docked pose of the minimum energy (fullfitness score) conformers of the four anti-virals, namely, saquinavir, azidothymidine, indinavir, and tipranavir with SARS-CoV-2 protease along with their corresponding 2 D interaction plots are depicted in the Figure 5 . It can be seen that the anti-virals bind within the active site of the SARS-CoV-2 M pro . It could be observed from Figure 5a , that saquinavir interacts with the residues in the active site through hydrogen bonding with SER46, ASN142 and GLU166, p-p T-shaped interaction with HIS41, p-alkyl interaction with CYS145, and van der Waals forces of attraction with other residues as depicted in the 2 D plot. Azidothymidine forms hydrogen bonds with GLY143, CYS145 and GLU166, p-p stacked interaction with HIS41, p-alkyl hydrophobic interaction with CYS145 and van der Waals forces with other residues as shown in the 2 D plot (Figure 5b ). Indinavir forms hydrogen bonds with ASN142 and CY145, p-sigma and p-sulphur interactions with MET49, p-p T-shaped interaction with HIS41, p-alkyl interaction with PRO168, and van der Waals interactions with other residues as depicted in 2 D plot (Figure 5c ). Tipranavir forms hydrogen bond and halogen bond with THR24 residue, p-sulphur with MET165, p-alkyl with PRO168, alkyl hydrophobic interaction with LEU27, and van der Waals forces with other residues as vivid from the 2 D plot ( Figure  5d ). On comparison of the estimated DG values from Table 1 for the four anti-virals, we observed that indinavir has the most negative DG value (À8.84 kcal/mol), followed by tipranavir (À8.42 kcal/mol), saquinavir (À8.35 kcal/mol) and azidothymidine (À7.40 kcal/mol). The electrostatic surface potential of the binding site along with the simultaneous presence of the four drugs is shown in Figure S5 . 

The docked pose of the minimum energy (fullfitness score) conformers of five anti-nematodes, namely, diethylcarbamazine, primaquine, mepacrine, artemisinin and niclosamide, and one anti-protozoal drug, emetine with SARS-CoV-2 protease along with their corresponding 2 D interaction plots are depicted in Figure 6 . It can be seen that the drugs bind to the active site of the SARS-CoV-2. Diethylcarbamazine ( Figure 6a ) interacts with LEU27, MET49 and CYS145 through alkyl hydrophobic interactions and van der Waals interaction with other residues. Primaquine (Figure 6b ) forms hydrogen bond with CYS44 residue along with p-sulphur interaction with CYS145 and van der Waals interaction with other residues. Mepacrine (Figure 6c ) undergoes alkyl hydrophobic interactions with CYS145 and MET165 residues, p-alkyl hydrophobic interactions with HIS163, MET165, MET49, and van der Waals interaction with other residues. Similarly, artemisinin ( Figure 6d ) undergoes alkyl hydrophobic interaction with MET49 and CYS145 residues along with p-alkyl interactions with HIS163. Niclosamide (Figure 6e) , forms a hydrogen bond with GLY143, p-sigma and p-sulphur interactions with GLN189 and MET165, respectively, alkyl hydrophobic interactions with LEU167 and PRO168, p-alkyl interaction with CYS145. The anti-protozoal drug, emetine undergoes attractive electrostatic (salt bridge) interaction with GLU166, alkyl hydrophobic interaction with MET49, and Met165, p-alkyl interaction with PRO168 ( Figure 6f ). On comparison of the estimated DG values from Table 1 for the five anti-nematodes and one anti-protozoal, we observe that emetine has the most negative DG value (À9.07 kcal/mol), followed by mepacrine (À7.80 kcal/mol), primaquine(À7.45 kcal/mol), artemisinin (À7.15 kcal/mol) and finally niclosamide (À6.77 kcal/mol). The electrostatic surface potential of the binding site along with the simultaneous presence of the six drugs is shown in Figure S6 . Herein, the blind molecular docking studies carried out for a total of 33 compounds under different categories reveal that they might act as potential anti-SARS-CoV-2 drugs, as these molecules bind to the active site of the main protease of SARS-CoV-2 M pro . Inhibition of this enzyme/protease is essential for blocking the viral replication . This inhibitor binding site within the protease of SARS-CoV-2 consists of residues HIS41, MET49, GLY143, CYS145, HIS163, HIS164, GLU166, PRO168, and GLN189 as evident from the recent study on the a-ketoamide inhibitors for SARS-CoV-2 main protease . We observe from the blind docking study of all the 33 molecules with the SARS-CoV-2 protease that the molecules are generally surrounded by the above mentioned residues, which clearly suggests that this molecule can prevent the viral replication of SARS-CoV-2. The distance of the ligands along with the change in accessible area of the two important catalytic residues (HIS41 and CYS145) within the active site of the protease is shown in Table 2 . Although the blind docking studies indicates that all the molecules can act as potential agents for COVID treatments, but from the estimated free energy of binding (DG) values (Table 1) , we can infer that rutin with the highest negative minimum DG value (À9.55 kcal/mol) among all the studied compounds can be the best possible SARS-CoV-2 inhibitor, followed by ritonavir (À9.52 kcal/mol), emetine (À9.07 kcal/mol), hesperidin (À9.02 kcal/mol), and indinavir (À8.84 kcal/mol).

The emergence of novel viruses requires efficient therapeutic strategies to develop or identify drugs that can treat the disease or control the further spread of the virus. In recent years, molecules that influence the life cycle of viruses have been identified, and a broad range of anti-viral agents have been proposed, which includes inhibiting the entry and stopping of viral replication or targeting intracellular signal transduction pathways (Elfiky & Azzam, 2020; Hasan et al., 2020) . In this context, repurposing of already known drugs and molecules is an essential concept (R. J. , especially in the pandemic that has resulted due to the novel coronavirus, as it is cost-effective in terms of research and development of new novel drug molecules. Sir James Black, the recipient of the 1988 Nobel Prize in Physiology and Medicine, had famously stated that "The most fruitful basis for the discovery of a new drug is to start with an old drug". In addition to the FDA approved drugs, here we observed that rutin and hesperidin, two dietary polyphenols have significant potential to function as inhibitors of SARS-CoV-2 M pro . Dietary polyphenols have low systemic toxicity and are highly beneficial for human health (Cory et al., 2018) . Rutin, also known as vitamin P, has several pharmacological properties including anti-viral activity (Ganeshpurkar & Saluja, 2017) . Similarly, hesperidin has also been observed to act as an anti-viral compound (Parvez et al., 2019) . Therefore, this finding also indicates a promising potential for the use of dietary polyphenols for the treatment of COVID-19.

This blind molecular docking study proposes a potential approach for the utilization of natural products, anti-virals, anti-fungals, anti-nematodal, and anti-protozoal drugs as prospective inhibitors of SARS-CoV-2 M pro , a pandemic and a global threat that is currently affecting millions of people, leading to death in extreme cases. We observe that the studied bind within the active site of the 6Y84, the main protease of SARS-CoV-2, which is essential for preventing the protein mutarotation of the virus by inhibiting the viral replication, and further spread of the infection. The estimated free energy of binding varied in the range of -6.57 to -9.55 kcal/mol, suggesting the molecules can spontaneously interact within the binding site of SARS-CoV-2 M pro . Among all the molecules, the inhibition activity of rutin was found to be the highest as its DG value (À9.55 kcal/mol) is the minimum, followed by ritonavir (À9.52 kcal/mol), emetine (À9.07 kcal/mol), hesperidin (À9.02 kcal/mol), and indinavir (À8.84 kcal/mol). As this study has been carried out using computational methods, therefore being particular about these results would require further wet-lab experiments to be carried out under in vivo as well as in vitro conditions.

No potential conflict of interest was reported by the author(s). Table 2 . Distance and the change in accessible surface area of the active site residues (HIS41 and CYS145) of SARS-CoV-2 M pro on interaction with the potential inhibitors. 

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

Biological activities of emetine

Molecular docking studies on the interaction of NCI anticancer analogues with human Phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit

Discovery studio visualizer

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

COVID-19) events as they happen

The role of polyphenols in human health and food systems: A mini-review

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

Approved antiviral drugs over the past 50 years

Novel Guanosine Derivatives against MERS CoV polymerase: An in silico perspective

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

Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: An in silico study

Funding infectious disease research: A systematic analysis of UK research investments by funders 1997-2010

The pharmacological potential of rutin

Hydroxychloroquine and azithromycin as a treatment of COVID-19: Results of an open-label non-randomized clinical trial

SwissDock, a proteinsmall molecule docking web service based on EADock DSS

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

A review on the cleavage priming of the spike protein on coronavirus by angiotensin-converting enzyme-2 and furin

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

Inhibition of SARS-CoV 3CL protease by flavonoids

Potential Inhibitor of COVID-19 Main Protease (M pro) from Several Medicinal Plant Compounds by Molecular Docking Study. Preprints (Www.Preprints.Org) j NOT PEER-REVIEWED j Posted

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

Drug repurposing: A better approach for infectious disease drug discovery?

Drug treatment options for the 2019-new coronavirus (2019-nCoV)

Outbreak of pneumonia of unknown etiology in Wuhan, China: The mystery and the miracle

Learning from the Past: Possible urgent prevention and treatment options for severe acute respiratory infections caused by 2019-nCoV

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

RCSB Protein Data Bank (PDB) ID

Repurposing drugs for the treatment and control of helminth infections

Peptide-like and small-molecule inhibitors against Covid-19

A review of computational drug repurposing

Plant-derived antiviral drugs as novel hepatitis B virus inhibitors: Cell culture and molecular docking study

UCSF Chimera-A visualization system for exploratory research and analysis

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

The PyMOL Molecular Graphics System, Version 2.0, Schr€ odinger, LLC

Pa-9: A flavonoid extracted from plectranthus amboinicus inhibits HIV-1 protease

Molecular docking using ArgusLab, an efficient shape-based search algorithm and the AScore scoring function

Investigation into SARS-CoV-2 resistance of compounds in garlic essential oil

Fluconazole versus itraconazole for antifungal prophylaxis in neutropenic patients with haematological malignancies: A meta-analysis of randomisedcontrolled trials

Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro

Comparative effectiveness of combined favipiravir and oseltamivir therapy versus oseltamivir monotherapy in critically ill patients with influenza virus infection

Antiviral drug discovery against SARS-CoV

Nelfinavir was predicted to be a potential inhibitor of 2019-nCov main protease by an integrative approach combining homology modelling, molecular docking and binding free energy calculation

Flavonoids: Promising natural compounds against viral infections

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

The authors are grateful to National Institute of Technology Meghalaya for providing research platform. SD and SS thank TEQIP III, National Institute of Technology Meghalaya for their fellowship. The study was not supported by any funding authority. The authors thank the reviewers for their valuable suggestions to increase the scientific quality of the manuscript.The authors declared that no competing conflict of interest exists. All authors have read and approved this version of the article.