key: cord-1032074-hwe3nc85
authors: Badavath, Vishnu Nayak; Kumar, Akhil; Samanta, Pralok K.; Maji, Siddhartha; Das, Anik; Blum, Galia; Jha, Anjali; Sen, Anik
title: Determination of potential inhibitors based on isatin derivatives against SARS-CoV-2 main protease (m(pro)): a molecular docking, molecular dynamics and structure-activity relationship studies
date: 2020-11-17
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1845800
sha: 1b182e2945dab1c005bdd0b6f1bf741e7c100aca
doc_id: 1032074
cord_uid: hwe3nc85

SARS-COV-2, the novel coronavirus and root of global pandemic COVID-19 caused a severe health threat throughout the world. Lack of specific treatments raised an effort to find potential inhibitors for the viral proteins. The recently invented crystal structure of SARS-CoV-2 main protease (M(pro)) and its key role in viral replication; non-resemblance to any human protease makes it a perfect target for inhibitor research. This article reports a computer-aided drug design (CADD) approach for the screening of 118 compounds with 16 distinct heterocyclic moieties in comparison with 5 natural products and 7 repurposed drugs. Molecular docking analysis against M(pro) protein were performed finding isatin linked with a oxidiazoles (A2 and A4) derivatives to have the best docking scores of −11.22 kcal/mol and −11.15 kcal/mol respectively. Structure-activity relationship studies showed a good comparison with a known active M(pro) inhibitor and repurposed drug ebselen with an IC(50) value of −0.67 μM. Molecular Dynamics (MD) simulations for 50 ns were performed for A2 and A4 supporting the stability of the two compounds within the binding pocket, largely at the S1, S2 and S4 domains with high binding energy suggesting their suitability as potential inhibitors of M(pro) for SARS-CoV-2.

During December 2019, a novel coronavirus named "Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2), caused an outbreak of a respiratory disease in the city of Wuhan, capital of the Hubei province in China and since has been spreading globally (Wu et al., 2020b; Zhou et al., 2020) and the pulmonary disease caused by SARS-CoV-2 is known from human-to-human led to an exponential growth in the number of patients with world-wide spreading of the virus. According to the World Health Organization (WHO), the virus infects most of the patients with mild symptoms such as fever, cough and difficulty in breathing that are cured with proper care. However, in older patients or in patients with underlying health conditions, the disease can progress into fatal pneumonia and acute respiratory failure Zhou et al., 2020) . The progress of the virus was due to its novel transfer from human to human through cough, sneeze, touch from an infected person i.e. through the common droplet infection (Chen et al., 2020c; Wu et al., 2020b; Zhou et al., 2020) . As of July 2nd, 11,000,000 cases and over 525,000 deaths worldwide exerting an enormous social impact and severely effecting world economic status, as previous pandemics (Keogh-Brown & Smith, 2008) .

Currently, many researchers have been working to develop specific drugs through repurposing of drugs for diseases like HIV, Ebola, malaria etc. (Fischer et al., 2020; Wang et al., 2016; Xu et al., 2020b; Zhang et al., 2020a) . Still these significant developments have yet to achieve prevention of disease spreading. A key strategy for treatment viral infection has been inhibiting its main protease, such inhibitors are in clinical use as effective treatments for human immunodeficiency virus (HIV) and hepatitis C (Ghosh et al., 2016; Yang et al., 2006) .

The SARS-CoV-2 virion (Scheme 1) has four structural proteins (Wu et al., 2020a) spike proteins (S), Envelope (E), Membrane protein (M) and nucleocapsid proteins (N) which holds the RNA genome, similar to the other coronaviruses. The Spike protein and M pro (also called as 3CL pro ) are the two main sites for drug development for SARS-CoV-2. M pro plays a key role in mediating viral replication and transcription Xu et al., 2020b; Jin et al., 2020b; Xu et al., 2020b) . M pro operates at 11 cleavage sites on the large polyprotein 1ab (replicase 1ab, $790 kDa), inhibiting the activity of this enzyme directly blocks viral replication (Jin et al., 2020b) . As there are no human proteases with a homolog of M pro , it is an ideal target for drug design as inhibitors are less likely to be toxic for humans . Thus, several studies describing M pro inhibitors have been recently published Fan et al., 2020; Jin et al., 2020b; Raugi et al., 2016; Xu et al., 2020b) .

The crystal structure of SARS-CoV-2 M pro was recently solved and deposited to the protein data bank (PDB-6LU7) enabling rational design of specific inhibitors (Bz owka et al., 2020) . There is a strong resemblance between SARS-CoV-2 and SARS-CoV-1, thus, many researchers were motivated to develop drugs based on their activities of the SARS-CoV-1 (Dai et al., 2020; Xu et al., 2020b) . While some of the M pro inhibitors are in preclinical or early clinical stages, there have been no approved treatment yet Pillaiyar et al., 2016) . Hence, we were encouraged to perform a virtual screen of inhibitors that have a wide structural diversity and include a variety of structures to generate a large pool of potential M pro inhibitors.

According to our literature survey, the presence of heterocyclic moieties was found in most of the potential drugs designed inhibitors for SARS-CoV-2 M pro . Here, we have screened a wide variety of heterocyclic compounds including 118 previously studied anticancer, anti-bacterial and antimicrobial compounds in comparison to the published M pro inhibitor N3 along, 5 natural products and 7 repurposed drugs.

The 7 repurposed drugs are presently used by various researchers and medical institutes for the treatment of COVID-19, like Chloroquine (Cortegiani et al., 2020; Mallikarjuna et al., 2020) , Hydroxychloroquine (Gautret et al., 2020) , Remdesivir (Scavone et al., 2020) , Favipiravir (Shiraki & Daikoku, 2020) , Niclosamide , Ebeselen (Jin et al., 2020a) and Eidd-2801 (Sheahan et al., 2020) . Five natural compounds like Ursodeoxychloric acid, Quercetin, Kaempferol, Pinocembrin and Rutin are also included in this study. Chloroquine and Hydroxychloroquine are known drugs for treatment of malaria (Ben-Zvi et al., 2012; Bhattacharjee, 2016; Cortegiani et al., 2020) . Remdesivir and Favipiravir are antiviral drugs potentially against several viruses (Scavone et al., 2020; Shiraki & Daikoku, 2020) .

The broad spectrum antiviral drug Remdesivir, developed by the Gilead Sciences and the anti malaria drug Chloroquine were recently studied in an in vitro research work by G. Xiao et al., to be effective inhibitors for the 2019-nCoV . Chloroquine was also previously predicted to be an inhibitor for the SARS-CoV-1 (De Clercq, 2006) . The other malaria drug Hydrochloroquine was also found to be effective in inhibiting SARS-CoV-2 infection in vitro and attenuate inflammatory response as given in a recent study ). An experimental antiviral drug developed for the treatment of influenza Eidd2801 was also repurposed for SARS, MERS and SARS-CoV-2, still in pre-clinical trials though very fruitful results are not expected (Sheahan et al., 2020) . Niclosamide, a drug designed to treat tapeworm infestations, is known for its broad antiviral activity, recent studies are performed to understand its activity against SARS-CoV-2 . The natural compound Ursodeoxychloric acid is a mammalian bile acid found in bear and is used in treatment of cholesterol (Einarsson, Scheme 1. Schematic representation of the structure of SARS-CoV-2. 1994). Quercetin, Kaempferol, Pinocembrin are natural flavonoids found in several fruits and vegetables and used as antioxidants (Wang et al., 2018) . Pinocembrin is also found in honey (Rasul et al., 2013) . Rutin is a bioflavonoid, found in many citrus food and have powerful antioxidant properties (Enogieru et al., 2018) .

Many of the screened drug molecules for our docking studies were taken from previously studied drugs, which were designed for anticancer activities. Some of them were also known for their antimicrobial and antibacterial activities (Afifi et al., 2017; Agarwal et al., 2016; Beedie et al., 2015; Bhatt et al., 2020; Blair et al., 2007; Fischmann et al., 2008; Gerova et al., 2016; Ivanova et al., 2007; Kok et al., 2008; Lane et al., 2001; Mendel et al., 2003; Pessoa et al., 2010; Schelman et al., 2011; Yadagiri et al., 2015; Zhang et al., 2011) . Many of these drugs are already FDA approved drugs like Sunitinib (Bhatt et al., 2020) and many are in their preclinical trials like Zibotentan (Schelman et al., 2011) .

Compounds based on oxidiazoles and isatin moieties, well known for vigorous research for anticancer, antiparasitic activities and such type of compounds were also screened on the basis of their structure similarities (Bhatt et al., 2020; Ivanova et al., 2007; Lakshmithendral et al., 2019; Mendel et al., 2003; Serafim et al., 2017; Venkatarao et al., 2019; Zhang et al., 2011) . Other heterocyclic moieties with azo based structures and antibacterial, antitumor or anticancer activities were also screened (Bhatt et al., 2018; Jha & Ramarao, 2018; Tamokou et al., 2016) . Antimicrobial compounds based on benzylpyrrolidine groups have also been screened for our cause (Sreekanth & Jha, 2020) . As for comparison our in-silico screen included the very recent virtually screened 11 compounds that showed inhibitors of SARS-CoV2 (Fischer et al., 2020) . In addition we selected compounds with a new heterocyclic ring system, piperdin-4-one semicarbazones, that exhibits moderate antibacterial and antifungal activities (Anand et al., 2019) . A recent study by Zhang et al. (2020b) showed that a-keto amide inhibitors are very important inhibitors for MERS and predicted it to be better for SARS-CoV-2 as well. We have chosen the best structure (A118) for our docking analysis.

Our docking analyses showed high interaction of 12 novel heterocyclic drug moieties with SARS-CoV-2 main protease M pro that can serve as potential novel drug candidates. Further Molecular dynamic (MD) analyses of the two best docked compounds A2 and A4, both derivatives of isatin and oxadiazole, showed that the interactions for these compounds are mostly governed by the hydrophobic and electrostatics interactions. The electrostatic charges and Molecular electrostatic potentials were also calculated for these compounds using Density Functional Theory (DFT) to understand their charge electrostatic charge, shape and the surface of the ligands. These calculations were done to predict the electronic distribution of the ligand which is one of the most important factor for analysing the binding between the protein and the ligand. Drug like properties of the compounds, ADME profile analysis, were also performed for these compounds.

This study suggests two compounds that bind M pro extremely potently in-silico and may serve as the basis for potent inhibitors of COVID-19. Our calculations also shed light on the shape, electrostatic charge and surface of the binders showing adequate ADME properties.

For our virtual screen of inhibitors of M pro we selected structures of a library of newly designed compounds, A1-A118, in addition to natural products used as antiviral, anticancer, antibacterial drugs and COVID-19 drugs in clinical trials ( Figure S1 ), as discussed in the introduction. All the structures were optimized with Gaussian 16 software at B3LYP/6-31 þ G(d) level of theory (Gibbs free energy and the zero point energies are also given in Figure S1 ). Drugs used in clinical trials with specific crystal structures were taken directly without optimizing. All the structures were subjected to the docking analysis with the main protease M pro of the SARS-CoV-2, PDB id: 6LU7. The free energy of binding (kcal/ mol) and inhibitory constant (Ki) for these compounds are given in Table 1 . The active site of M pro is consists of five sub-sites (S1, S2, S1', S2', S4). S1 site consists of residues (Phe140, Leu141, Gln166, His163, Met165 and His172). The bulky hydrophobic S2 subsite (Met49 and Asp187, Glu189) (Wang et al., 2016) , while S1' (Thr25, Leu27, Cys38, Pro-39, Val42, and Cys145) is essential for catalysis (Xue et al., 2008) , S2' (Thr-26, Asn-28, Tyr-118, Asn-119, and Gly-143) and flexible S4 site (Leu167, Gln192) in nature (Mittal et al., 2020; Wang et al., 2016) .

Our obtained results were also compared with the recently published inhibitors for possible interaction with 6LU7-the main protease (M pro ) of the SARS-CoV-2. The docking analyses were also performed for the native drug N3 (À8.88 kcal/mol) and ebselen (À6.45 kcal/mol) as given in Table 1 . Furthermore, we compared these to the FDA approved drugs Remdesivir (À8.42 kcal/mol), Chloroquine (À7.01 kcal/mol); Hydroxychloroquine (À7.39 kcal/mol); Niclossamide (À7.60 kcal/mol) and EIDD-2801 (À6.93 kcal/ mol). The natural flavonoids Quercetin (À7.94 kcal/mol), Kaempferol (À7.96 kcal/mol), Pinocembrin (À7.76 kcal/mol) and Rutin (À8.12 kcal/mol) also showed very similar scores like the FDA approved drugs (Table 1) . Ursodeoxychloric acid (À9.42 kcal/mol) showed a higher docking score than all the other natural products and the FDA approved drugs featured (Table 1) . Our calculations (À8.42 kcal/mol) with AutoDock4.2 showed comparable interactions with the main protease of the SARS-CoV-2 with the Remdesivir drug by Ji et al., using Schrodinger docking suites (À7.215 kcal/mol) (Hall Jr & Ji, 2020 ).

The two best inhibitors found from our docking studies were A2 (À11.22 kcal/mol) and A4 (À11.15 kcal/mol) with calculated Ki of 5.94 nM and 6.66 nM, respectively, see structures in Table S1 . The oxodiazole moiety of A2, oriented toward the S1 site and establishes a hydrogen bond with residue Glu166, and weakens the salt bridge between positively charged N-terminal Ser1 (one monomer) of the negatively charged Glu166 (another monomer), thus destablilizes the M pro pocket (Huynh et al., 2020) this may be the reason for increasing the potency of compound A2. The Ebselen a Mpro inhibitor FDA approved drugs Remdesivir, Chloroquine and Niclosamide were also formed a similar type of hydrogen bonding with Glu166. The the Phenol segment of A2 has formed two H-bond with Ser144 by fully occupying S1' subsite. The chlorobenzene segment orients toward the S4 pocket, while the isatin segment is oriented towards the bulky hydrophobic S2 site ( Figure 1 ). Compound A4 showed a different picture with its benzo[d]oxazol-2(3H)-1 segment, the benzyl moiety orienting towards the S2' subsite while the oxazole moiety has accommodated well in S1'subsite by forming a strong hydrogen bond between the ring oxygen with the catalytic Cys145 residue; The chlorobenzene segment near oxazole moiety orients toward the S1 site and form a p-p interaction with the residue His172. The connecting 2 amine segment of the A4 is oriented towards the bulky hydrophobic S1 site and forms a hydrogen bond with the residue His164. The chlorobenzene segment near oxadiazole moiety orients towards the S4 pocket ( Figure 1 ). In Figure 1 , the other higher docked compounds like A2, A4, A8 and A20 are also represented. The best docked compounds A2 and A4 have been selected for further analysis with MD simulations using GROMACS.

In this article, we performed screening of 16 distinct heterocyclic skeletons (Table 2 ) with a total of 118 compounds, 5 natural products and 8 repurposed drugs (including Ebselen, which has a similar structure to isatin (Jin et al., 2020b) ) ( Figure 2 ). The docking analysis of several isatin compounds A2, A4, A8, A20, showed the best interaction with the M pro (above À10.84 kcal/mol) as presented in Table 1 and Figure 1 . In addition, highly potent compounds with Ki under 10 nM, A2, A4 and A8 were comprised of three separate ring systems, branched from a central carbon enabling multiple binding interactions with M pro .

According to a literature survey, isatin based compounds are well known for their anticancer (Bhatt et al., 2020; Popp, 1969) and antibacterial activities (Pandeya et al., 2000; Varma & Nobles, 1967) . They have also been known for their antiviral activities for a various pathogen virus (Bauer & Sadler, 1960; Burger, 1979) 89, 90 including HIV (Pauwels et al., 1988; Teitz et al., 1994) and SARS-CoV-1 (Selvam et al., 2008) . Oxadiazole, based derivatives have also been found very important for potent biological functions. Especially antiviral, Table 1 . The calculated Ki, Free binding energy and bond length of hydrogen bonds and p-p interaction of selected FDA approved drugs, natural products containing heterocyclic moieties compared with the top 12 compounds with the best binding affinity to SARS-Co-V-2 M pro . anti-HIV, anti-HCV, anti-HBV, anti-HSV activities, etc (Dong et al., 2016; Li et al., 2011) . Triazole derivatives also occupy a pivotal position in modern medicinal chemistry and several derivatives are well known for their applications in medicine especially anticancer and antiviral drugs (Kharb et al., 2011; Xu et al., 2019) . The specific reason behind the isatin; oxadiazole and triazole scaffold acting as potential anticancer, antiviral drugs are their strong tendency to accept multiple hydrogen bonds from their target protein. The Oxidiazole based compounds A1, A2, A4, A5, A7-A9, A11 & A12 and triazole based compound A20 showed good binding energies and Inhibition constants. Previously designed Cyclindependent kinases inhibitors known for their anticancer activity, A38 and A40 also showed high binding interaction with the M pro of the SARS CoV-2, though their inhibition constants were higher than the isatin oxidiazoles/triazoles based compounds (Table 1) . Compound A68, also showed high scores in our calculation (À10.82 kcal/mol), also comprised of a three separate ring system. A68 inhibitor was already reported by Fischer et al. 10 . Other compounds with very efficient docked scores that were subjected to further analysis, from our docking calculations are A1 (À10.67 kcal/ mol); A5 (À10.25 kcal/mol); A7 (À10.75 kcal/mol); A8 (À10.97 kcal/mol); A9 (À10.44 kcal/mol); A11 (À10.43 kcal/ mol); A12 (À10.87 kcal/mol); A20 (À10.84 kcal/mol); A38 (À10.00 kcal/mol) and A40 (À10.41 kcal/mol). All other compounds with lower docked scores are given in the supporting information (Table S2) . A2 (À11.22 kcal/mol) and A4 (À11.15 kcal/mol) compounds showed the highest docked scores and were subjected to further analysis. 

Molecular recognition between a protein and small molecules (ligand) occurs at their surfaces. The extent of binding forces in the protein-ligand complex largely depends on the electronic distribution of the ligand which can be predicted by the electrostatic charges of the ligand molecule. The electrostatic charges of the electronegative atoms of the A2 and A4 ligands are given in Figure 3 . Such molecular details of the ligands helps predicting the binding probabilities of the ligands in the S1, S2, S1', S2' and S4 pockets as discussed above. Our charge analysis shown that more electronegative atoms are present in A2 predicting it to interact better with the protein than A4. The ChelPG charges of the electronegative atoms for the other potential ligands are given in the supporting information ( Figure S1 ). As observed from the charge analysis, A1-an analogue of A2, showed less interaction to M pro due to the presence of chlorine atom in place of the hydroxide. A5, an analogue of A4, showed low interactions despite of the presence of hydroxide ion. This may be due to the steric hindrance between the isatin and the chlorobenzene occurring, owing to the flexibility of bonds. A7, A8, A9, A11 and A12 having the similar structures showed very comparable interactions. The ChelPG charge analysis could not be performed for A7, A8 and A9 ligands due to restrictions in the Gaussian format with presence of bromine atom. The ligand A20, having a distinct structure with a hydrophobic end, showed higher interaction to the protein. The ChelPG charges of the other important ligands (A38, A40 and A68) with good interactions to the M pro are also given in the Figure S1 in the supporting information. MESP analysis shows A2 can form at least three hydrogen bonds, while A4 can form two (red patches). The electron density of the A2 was found to be higher than the A4 (green to yellow). The incorporation of an OH group in A2, rather than a Cl as in A4, is the reason behind the increasing electron density of A2, which should increases its binding interactions. Such enhancement of the electron density is favourable for the p-staking interaction. The Cl atom for both compounds also occupies a large electron density surface, showed in green, with negative values which is also suitable for multiple p-staking interactions. The small red, yellow and blue patches on the large green surface of compounds are balancing the hydrophilic and hydrophobic parts which are essential for good binding to the protein. The MESP analysis of the other ligands as given in Table 1 are also calculated and are given in the Figure S2 , in the supporting information. The red, yellow and blue patches on the large green surface of compounds in these ligands are observed and the interactions patterns are clear from the MESP.

Further calculations using MD simulation tools were performed for A2 and A4 to check the interaction analysis and binding ability with the M pro protein of the SARS-CoV-2 to evaluate the results of the docking and the DFT charge analyses.

The stability of each simulated model was determined by the calculations of backbone atom RMSD for ligand À M pro complex as shown in Figure 5 . The RMSD of A2 and A4 with-M pro rises till 20 ns and stabilizes after 10 ns at 2 Å. This clearly indicates that ligands are stabilized inside the pocket and don't change their orientation in the active site of M pro . Similarly protein backbone atoms for A2 showed stability around 2.5 Å. However, A4 showed higher RMSD and maintained at 3 Å, which directly resembels to the charge analyses. Furthermore, A2 and A4 hydrogen bond with M pro are plotted along the 50-ns MD simulation, shown in Figure 6 . These simulations indicated that on an average A2 and A4 show approximately two hydrogen bonds and are in accordance with our docking results and further suggested the stability of the ligand in active site.

Other than the Hydrogen bond, hydrophobic interactions were established among non-polar amino acids and provide the stability of proteins in solution by shielding the nonpolar amino acids in hydrophobic cores, away from the aqueous environment. Solvation plays an important role in ligand-protein association and has a strong impact on comparisons of binding energies for dissimilar molecules.

Change in Solvent Accessible Surface Area (SASA) directly indicates the change of hydrophobic core caused by changes in tertiary structure. The calculation of the SASA of A2 and A4 do not show any drastic change that indicate the ligand were stable during the simulation run, as shown in Figure 7 , and showed very similar hydrophobic area. This indicates that A2 and A4 are less exposed to water molecules and buried in M pro active site.

The calculated binding energy using the mmpbsa approach indicates that A2 has higher energy (À48.33 kcal/mol) in comparison to A4 (-33.10 kcal/mol), Table 3 , similar to the finding in the docking analysis, Table 1 . Individual components of energy for each system indicate that A2 and A4 both show higher van der Waals energies and very comparable SASA energies and polar salvation energies. However a significant difference was observed in the electrostatic energies. As explained in the ChelpG charge analysis (Figure 3) , the electrostatic interactions for the A2 showed much higher interactions than the A4 ligand, which effects the binding energy difference of the two systems. Thus, the binding energies are governed by the van der Waals energy and electrostatic energy.

To identify important interacting residues with A2 and A4 the contribution of each residues to the binding energy was calculated and plotted ( Figure 8 ) indicating three important regions: amino acids 26 to 56, amino acids 136 to 148, and amino acids 162 to 194. These regions showed significant energy contribution to the binding energy. Leu27 Asn51 and Pro52 showed negative energy in both systems while His41 showed positive energy in A2. In the second region (amino acid 136 to 148) A4 showed higher energy values for residues Leu141, Asn142, Gly143, Ser144, Cys145. In the third region (amino acid residue 162 to 194) A2 had higher binding energy values for interactions with residues Met165, Leu167, Pro168, Phe185, Asp187, Thr190, Ala191 and Gln192. While A2 interacted with amino acid from region 162 to 192 and residue 44, A4 showed good interaction with region 140 to 146. The M pro complex interacts with the ligands A2 and A4 with the different amino acid residues between energy À2.4 kcal/mol to À0.24 kcal/mol. A2 formed more interaction with high binding energy (up to 2.4 kcal/mol) in the active site, Figure 9 . Cys44, Met165, Gln189 and Ala191 are major residues that highly contributing to binding with M pro , this indicates the ligand fits well into the active site of M pro . However, residues participating in A4 binding are somewhat different, very few residues bind with more than À0.95 kcal/ mol, the major binding residue were found to be Gln189, Met165, Cys145 and Pro168. Several amino acids were common to binding both A2 and A4: Leu27, Cys44, Met49, His164, Met165, Leu167, Asp187, Gln189, Thr190, Ala191 and Gln192 these also contribute to the binding energy. Overall this analysis supports the better binding interactions of the A2 on the active site of M pro than A4, Figure 10 .

The RMSF captures, for each residue, the fluctuation from its average position, gives insight into the flexibility of regions of the M pro . The fluctuations per residue of COVID-19 M pro indicated that after the binding, three regions, amino acid residues from 40 to 54, 124 to 144 and 180 to 200, and in addition 272 to 306 showed higher fluctuations, Figure 11 . Moreover, these residues are also participating in ligand binding. These might suggest that these flexible regions are important in substrate or ligand binding.

In order to screen the pharmacokinetic behaviour, the in-silico ADME parameters of the designed ligands with free energy of binding higher than À10 kcal/mol were calculated using the SwissADME software (http://www.swissadme.ch). Lipinski's rule of five is a prerequisite to ensure drug-like properties when using rational drug design. The designed derivatives follow the four properties of Lipinski's rule of five (mol. wt. 500 Da; log P o/w 5; HBD 5; HBA 10). All the ligands have showed ADME parameters in an acceptable range and possess significant drug-like characteristics based on Lipinski's rule of 5, Table 4 .

In conclusion, we have studied the binding of a series of heterocyclic moieties to the active site of M pro of the SARS-CoV2. Molecular docking analysis show that the isatin based derivatives possess excessive interactions to M pro , with binding energies of over À10 kcal/mol. The best docked ligands A2 and A4 consisted of an Isatin moiety along with an oxidiazole ring, and both moieties are highly used in anticancer, antiviral, antibacterial drug development. Further, the binding interactions were more vividly discussed with a comprehensive MD simulation analysis using GROMACS, and binding energies were calculated using the MMPBSA method on the active site of the protein. These analyses show hydrogen bonds and p-p interactions with the residues His41, Cys145, Val42, Cys44, Thr45, Ser46, Glu47, Asp48, Met49 and Glu166, Met165, Leu167 and Pro168. Our analysis also shows that hydrophobic and electrostatic interactions are most important to predict the binding energies of the ligands to the protein and that A2 interacts more strongly than the A4. This finding was reinforced by the electrostatic charge analysis and MESP analysis through DFT calculations. All the ligands have showed ADME parameters in an acceptable range and possess significant drug-like characteristics based on Lipinski's rule of 5. Our findings indicate that both isatin and oxidiazoles based drugs are indeed very good candidates as potential inhibitor for M pro of SARS-CoV-2 and A2 and A4 compounds are the best candidates with extremely high potency in silico and further in vitro analysis will be performed as future studies.

All the ligands were initially optimized using Density functional theory with Gaussian 16 suite of programs to obtain the most stable structure of each system (Ogliaro et al., 2016) except the repurposed drugs with crystal structures available. All the compounds were optimized at B3LYP/6-31 þ G(d) (Becke, 1993; Lee et al., 1988) level of theory in gas phase. Positive vibrational frequencies confirmed the structures to be minimum. PDB structures were saved using Gaussview and used for the Molecular docking studies. Electrostatic charges were calculated for the specific compounds using the ChelpG charge analysis (Martin & Zipse, 2005; Rozas, 1997; Tsuzuki et al., 1996) implemented in Gaussian software. The Molecular Electrostatic Potential (MESP) was generated for both the ligands at the same level of theory.

The crystal structure of SARS-CoV-2 has been retrieved from the rcsb.org (PDB ID: 6LU7) (Jin et al., 2020b ) and used to generate initial 3D coordinates. Co-crystallized water molecules were deleted along with addition of polar Hydrogen and compute gasteiger charges. The structures of selected inhibitors were superimposed against pre-docked ligand in the PDB, and the latter was then removed to generate initial conformation of inhibitor at the active site of SARS-CoV-2. Grid box was then determined by the native ligand position on the binding site (Thr26, Tyr54, Phe140, Asn142, Gly143, Ser144, Cys145, His163, His164, Glu166, and His172). Finally both the Auto-grid and Auto-dock were run with the default parameters (Badavath et al., 2016a) . All molecular docking analysis were performed by AutoDock4.2 using earlier reported protocol ( Gangireddy et al., 2019; Munusamy et al., 2019; Nath et al., 2018; Nayak et al., 2015; Vishnu Nayak et al., 2013) software version 4.2 (Forli et al., 2012) . After completing the molecular docking process, the generated DLG (.dlg) file was analyzed using AutoDock v4.2 tools for all the possible interaction and binding free energies following the methodology (Huey et al., 2007) . Top scoring molecules from largest cluster were considered for the interactions. used for MD simulation studies of two system 50 ns each. The PRODRG2 Server (Schuttelkopf & Van Aalten, 2004 ) was used to generate the topology of A2 and A4. Each system was placed in the centre of cubic box having distance of 1.0 nm between protein and edge of the simulation box and solvated with explicit water molecules, with the SPC solvation model (Mark & Nilsson, 2001) . Before proceeding towards energy minimization, all the systems were neutralized by adding Na þ /Cl À ions accordingly. Steepest descent method was used for energy minimization of each system. MD simulations with NVT (isochoric-isothermal) and NPT (isobaric-isothermal) ensembles (N ¼ constant particle number, V ¼ Volume, P ¼ Pressure, T ¼ Temperature) were performed for 1 ns to equilibrate the protein-ligand system for constant volume, pressure (1 atm) and temperature (300 K). To calculate electrostatic interaction, Particle Mesh Ewald (PME) algorithm (Darden et al., 1993; Essmann et al., 1995) was used with grid spacing of 1.6 Å and a cutoff of 10 Å and LINCS method (Hess et al., 1997) was used to restrain the bond length. Finally, 50 ns of production MD was performed for each system and trajectories were saved at every 2 ps.

G_mmpbsa (Kumari et al., 2014) was developed using two widely used open source software i.e. GROMACS and APBS and it has similar user interface like other GROMACS tools. The tool calculates components of binding energy using MM-PBSA method except the entropic term and energetic contribution of each residue to the binding using energy decomposition scheme. Binding free energy estimation of all the docked complexes was done by MM/PBSA method through g_mmpbsa module. For each simulated system, to calculate average binding energy by using bootstrap analysis, a total of 200 snapshots were taken from the last 10 ns of the trajectory on the intervals of 50 ps. The module calculates electrostatic interactions, van der Waals interactions, polar solvation energy and non-polar solvation energy.

ADME parameters of the designed ligands were calculated using the SwissADME (http://www.swissadme.ch) (Daina et al., 2017) accessed on-5 th May 2020. All the compounds were drawn in Chem Draw software and smiles were taken as input, before being used in SwissADME. Lipinski's rule of five (Lipinski et al., 1997) .

Structure-activity relationships and molecular docking studies of chromene and chromene based azo chromophores: A novel series of potent antimicrobial and anticancer agents

Design and synthesis of new 2,5-disubstituted-1,3,4-oxadiazole analogues as anticancer agents

Synthesis, spectroscopic investigation, computational, stereochemical and biological studies of 1, 3-dimethyl-2, 6-diphenylpiperidin-4-one N (4 0 )-cyclohexylsemicarbazone: Crystal structure and Hirshfeld surface analysis

Monoamine oxidase inhibitory activity of novel pyrazoline analogues: Curcumin based design and synthesis

Monoamine oxidase inhibitory activity of ferulic acid amides: Curcumin-based design and synthesis

Monoamine oxidase inhibitory activity of novel pyrazoline analogues: Curcumin based design and synthesis-II

Monoamine oxidase inhibitory activity of novel pyrazoline analogues: Curcumin based design and synthesis

Design, synthesis and MAO inhibitory activity of 2-(arylmethylidene)-2,3-dihydro-1-benzofuran-3-one derivatives

The structure-activity relationships of the antiviral chemotherapeutic activity of isatin beta-thiosemicarbazone

Density-functional thermochemistry. III. The role of exact exchange

Anticancer properties of a novel class of tetrafluorinated thalidomide analogues

Hydroxychloroquine: From malaria to autoimmunity

Design, synthesis and anticancer evaluation of oxa/thiadiazolylhydrazones of barbituric and thiobarbituric acid: A collective in vitro and in silico approach

Design and ultrasound assisted synthesis of novel 1, 3, 4-oxadiazole drugs for anti-cancer activity

Chemistry of antibiotics and related drugs

Structure-guided development of affinity probes for tyrosine kinases using chemical genetics

Structural and evolutionary analysis indicate that the SARS-CoV-2 Mpro is a challenging target for small-molecule inhibitor design

Synthesis and evaluation of isatin derivatives as effective SARS coronavirus 3CL protease inhibitors

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Prediction of the 2019-nCoV 3C-like protease (3CLpro) structure: Virtual screening reveals velpatasvir, ledipasvir, and other drug repurposing candidates

Diagnosis and treatment recommendations for pediatric respiratory infection caused by the 2019 novel coronavirus

A systematic review on the efficacy and safety of chloroquine for the treatment of COVID-19

Structure-based design of antiviral drug candidates targeting the SARS-CoV-2 main protease

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

Particle mesh Ewald: An NÁ log (N) method for Ewald sums in large systems

Potential antivirals and antiviral strategies against SARS coronavirus infections

Synthesis and antiviral activity of novel thioether derivatives containing 1, 3, 4-oxadiazole/thiadiazole and emodin moieties. Phosphorus, Sulfur, and Silicon and the Related Elements

Effect of ursodeoxycholic acid on hepatic cholesterol metabolism

Rutin as a potent antioxidant: Implications for neurodegenerative disorders

A smooth particle mesh Ewald method

Repurposing of clinically approved drugs for treatment of coronavirus disease 2019 in a 2019-novel coronavirusrelated coronavirus model

Potential inhibitors for novel coronavirus protease identified by virtual screening of 606 million compounds

Structure-guided discovery of cyclin-dependent kinase inhibitors

AutoDock (Version 4.2)

Design and synthesis of piperazine-linked imidazo [1, 2-a] pyridine derivatives as potent anticancer agents

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

Combretastatin A-4 analogues with benzoxazolone scaffold: Synthesis, structure and biological activity

Recent progress in the development of HIV-1 protease inhibitors for the treatment of HIV/ AIDS

Severe acute respiratory syndrome-related coronavirus: The species and its viruses-a statement of the Coronavirus Study Group

A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease

LINCS: A linear constraint solver for molecular simulations

GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation

A semiempirical free energy force field with charge-based desolvation

In silico exploration of the molecular mechanism of clinically oriented drugs for possibly inhibiting SARS-CoV-2's main protease

Cytotoxic Mannich bases of 6-(3-aryl-2-propenoyl)-2(3H)-benzoxazolones

Expeditious microwave-assisted synthesis and bio-evaluation of novel bis (trifluoromethyl) phenyl-triazole-pyridine hybrid analogues by the click chemistry approach

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

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

The economic impact of SARS: How does the reality match the predictions? Health Policy

Recent advances and future perspectives of triazole analogs as promising antiviral agents

Synthesis and anti-cancer activity of benzothia

g_mmpbsa-A GROMACS tool for high-throughput MM-PBSA calculations

Design, synthesis and biological evaluation of 2-(phenoxymethyl)-5-phenyl-1,3,4-oxadiazole derivatives as anti-breast cancer agents

A Novel cdk2-selective inhibitor, SU9516, induces apoptosis in colon carcinoma cells

Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density

Therapeutic drugs targeting 2019-nCoV main protease by high-throughput screening

1,3,4-oxadiazole: A privileged structure in antiviral agents

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

Hydroxychloroquine, a less toxic derivative of chloroquine, is effective in inhibiting SARS-CoV-2 infection in vitro

Potential inhibitors against 2019-nCoV coronavirus M protease from clinically approved medicines

Screening of chloroquine, hydroxychloroquine and its derivatives for their binding affinity to multiple SARS-CoV-2 protein drug targets

Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K

Charge distribution in the water molecule-a comparison of methods

In vivo antitumor activity of SU11248 antitumor activity of SU11248, a novel tyrosine kinase inhibitor targeting vascular endothelial growth factor and plateletderived growth factor receptors

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

Novel halogenated pyrido [2, 3-a] carbazole with enhanced aromaticity as potent anticancer and antioxidant: Rationale design and microwave assisted synthesis

Curcumin-based pyrazoline analogues as selective inhibitors of human monoamine oxidase A

Monoamine oxidase inhibitory activity of 2-aryl-4H-chromen-4-ones

Synthesis, antibacterial, antifungal and anti-HIV activities of norfloxacin Mannich bases

Rapid and automated tetrazolium-based colorimetric assay for the detection of anti-HIV compounds

Discovery of phthalimides as immunomodulatory and antitumor drug prototypes

An overview of Severe Acute Respiratory Syndrome-Coronavirus (SARS-CoV) 3CL protease inhibitors: Peptidomimetics and small molecule chemotherapy

Synthesis of potential antineoplastic agents. XX. Compounds related to the 3-o-nitrophenylhydrazone of isatin

Pinocembrin: A novel natural compound with versatile pharmacological and biological activities

Four amino acid changes in HIV-2 protease confer class-wide sensitivity to protease inhibitors

Atomic charges derived from different methods: A comparative study applied to SO2 heterocycles

Current pharmacological treatments for COVID-19: What's next?

A phase I study of zibotentan (ZD4054) in patients with metastatic, castrate-resistant prostate cancer

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

PRODRG: A tool for high-throughput crystallography of protein-ligand complexes

In vitro antiviral activity of some novel isatin derivatives against HCV and SARS-CoV viruses

Molecular modeling and structure-activity relationships studies of bioisoster hybrids of N-acylhydrazone and furoxan groups on cruzain

An orally bioavailable broad-spectrum antiviral inhibits SARS-CoV-2 in human airway epithelial cell cultures and multiple coronaviruses in mice

Favipiravir, an anti-influenza drug against life-threatening RNA virus infections

Microwave assisted synthesis and antimicrobial activity of novel pyrrolidine derivatives

Antibacterial and cytotoxic activities and SAR of some azo compounds containing thiophene backbone

Inhibition of human immunodeficiency virus by N-methylisatin-beta 4':4'-diethylthiosemicarbazone and N-allylisatin-beta-4':4'-diallythiosemicarbazone

Comparison of atomic charge distributions obtained from different procedures: Basis set and electron correlation effects

GROMACS: Fast, flexible, and free

Synthesis and antiviral and antibacterial activity of certain N-dialkylaminomethylisatin beta-thiosemicarbazones

Synthesis and biological evaluation of chalcone fused quinoline derivatives as anticancer agents

Monoamine oxidase inhibitory activity of 3,5-biaryl-4,5-dihydro-1H-pyrazole-1-carboxylate derivatives

Structure of main protease from human coronavirus NL63: Insights for wide spectrum anti-coronavirus drug design

Antitumor, antioxidant and anti-inflammatory activities of kaempferol and its corresponding glycosides and the enzymatic preparation of kaempferol

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

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

A new coronavirus associated with human respiratory disease in China

1, 2, 3-Triazole-containing hybrids as potential anticancer agents: Current developments, action mechanisms and structure-activity relationships

Broad spectrum antiviral agent niclosamide and its therapeutic potential

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

Structures of two coronavirus main proteases: Implications for substrate binding and antiviral drug design

Synthesis and evaluation of benzosuberone embedded with 1,3,4-oxadiazole, 1,3,4-thiadiazole and 1,2,4-triazole moieties as new potential anti proliferative agents

Synthesis, crystal structure, structure-activity relationships, and antiviral activity of a potent SARS coronavirus 3CL protease inhibitor

Deep learning based drug screening for novel coronavirus 2019-nCov. Interdisciplinary Sciences

a-Ketoamides as broad-spectrum inhibitors of coronavirus and enterovirus replication: Structure-based design, synthesis, and activity assessment

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

Synthesis, biological evaluation, and molecular docking studies of 1,3,4-oxadiazole derivatives possessing 1,4-benzodioxan moiety as potential anticancer agents

Isatin compounds as noncovalent SARS coronavirus 3C-like protease inhibitors

A pneumonia outbreak associated with a new coronavirus of probable bat origin

AK acknowledges Bioinformatics laboratory, Biotechnology Division of CSIR-CIMAP, Lucknow. AD, AJ and AS acknowledge to GITAM University, Visakhapatnam, India for the support. AS also acknowledges Kalyanashis Jana for fruitful discussions.

The manuscript was written through contributions of all the authors. All the authors have given approval to the final version of the manuscript.

There are no conflicts to declare. Pralok K. Samanta http://orcid.org/0000-0002-1926-7720 Galia Blum http://orcid.org/0000-0002-9374-2489 Anik Sen http://orcid.org/0000-0002-1336-5995