key: cord-304660-w7rs2dvt
authors: Bharadwaj, Shiv; Azhar, Esam Ibraheem; Kamal, Mohammad Amjad; Bajrai, Leena Hussein; Dubey, Amit; Jha, Kanupriya; Yadava, Umesh; Kang, Sang Gu; Dwivedi, Vivek Dhar
title: SARS-CoV-2 M(pro) inhibitors: identification of anti-SARS-CoV-2 M(pro) compounds from FDA approved drugs
date: 2020-11-05
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1842807
sha: 
doc_id: 304660
cord_uid: w7rs2dvt

Recent outbreak of COVID-19 pandemic caused by severe acute respiratory syndrome-Coronavirus-2 (SARS-CoV-2) has raised serious global concern for public health. The viral main 3-chymotrypsin-like cysteine protease (M(pro)), known to control coronavirus replication and essential for viral life cycle, has been established as an essential drug discovery target for SARS-CoV-2. Herein, we employed computationally screening of Druglib database containing FDA approved drugs against active pocket of SARS-CoV-2 M(pro )using MTiopen screen web server, yields a total of 1051 FDA approved drugs with docking energy >−7 kcal/mol. The top 10 screened potential compounds against SARS-CoV-2 M(pro) were then studied by re-docking, binding affinity, intermolecular interaction, and complex stability via 100 ns all atoms molecular dynamics (MD) simulation followed by post-simulation analysis, including end point binding free energy, essential dynamics, and residual correlation analysis against native crystal structure ligand N3 inhibitor. Based on comparative molecular simulation and interaction profiling of the screened drugs with SARS-CoV-2 M(pro )revealed R428 (−10.5 kcal/mol), Teniposide (−9.8 kcal/mol), VS-5584 (−9.4 kcal/mol), and Setileuton (−8.5 kcal/mol) with stronger stability and affinity than other drugs and N3 inhibitor; and hence, these drugs are advocated for further validation using in vitro enzyme inhibition and in vivo studies against SARS-CoV-2 infection. Communicated by Ramaswamy H. Sarma

Upper and lower respiratory infections caused by viruses are very common in temperate climates (Hing et al., 2006) and cause rhinitis, pharyngitis, sinusitis, bronchiolitis and pneumonia (Canducci et al., 2008; Jev snik et al., 2012) . Coronaviruses (CoVs) are a large group of RNA viruses with single-stranded RNA genomes that cause 30% of upper and lower respiratory tract infections in humans. However, in the last two decades, four human coronaviruses (HCoVs) have been identified, including severe acute respiratory syndrome (SARS)-CoV in 2003 (Fouchier et al., 2003) , HCoV-NL63 in 2004 (van der Hoek et al., 2004 , HCoV-HKU1 in 2005 (Vabret et al., 2006) , and MERS-CoV in 2012 (de Groot et al., 2013) . Recent COVID-19 pandemic caused by the fifth HCoVs considered as novel coronavirus (2019-nCoV), also known as SARS-CoV-2; suggested to emerged from Wuhan, Hubei Province of China (Giwa & Desai, 2020) , which is already announced as pandemic by the WHO Zhu et al., 2020) . Generally, person infected with SARS-CoV-2 exhibits clinical symptoms, such as fever, cough, shortness of breath and dyspnea; in severe case, infection can cause pneumonia, kidney failure, SARS and even death . However, there is currently no specific therapeutic drug or vaccine has been approved for the treatment of HCoVs; hence, under present conditions, HCoVs outbreak is a potential danger to the human health and world economy (Kirtipal et al., 2020) .

SARS-CoV-2 is a medium-sized, enveloped, positive-strand RNA virus ($30 kb) of genus Betacoronavirus and characterized by a crown-shaped (corona) appearance in electron micrographs of negatively stained preparations (Chang et al., 2016; Kim et al., 2020) . The viral genome deciphers several structural and non-structural proteins that assist virion in the infection to reproduce in a conserved linear arrangement (Kirtipal et al., 2020) . Since, HCoVs, including SARS-CoV-2, encoded the polyprotein which further releases the functional polypeptides through proteolytic digestion by 33.8 kDa main protease (M pro ), also known as 3-chymotrypsin-like cysteine protease, 3CL pro (Kirtipal et al., 2020; Yang et al., 2003) .

Hence, because of functional importance of M pro in the viral life cycle, viral M pro has been suggested as potential target in the antiviral development against SARS-CoV-2 (Cui et al., 2019; Du Toit, 2020) . Although, lack of therapeutic agents against SARS-CoV-2 posing a great difficulty to clinicians in disease management, but recent studies have highlighted the significance of existing drugs and their repurposing for disease management. For instance, some of noted recommendations are the synergistic use of antimalarial drugs, such as Chloroquine-Hydroxychloroquine (Colson et al., 2020) and Remdesivir-Favipiravir (Dong et al., 2020) . Additionally, Huanzhu Lu has suggested neuraminidase inhibitors, Remdesivir, Peptide (EK1), abidol, RNA synthesis inhibitors (TDF and 3TC), anti-inflammatory drugs (hormones and other molecules), and Chinese traditional medicine (ShuFengJieDu Capsule and Lianhuaqingwen Capsule) as possible treatment against SARS-CoV-2 (Lu, 2020) .

Provided with rapid surge in the information about the viral structure, mechanism and pathology, targeted drug discovery processes have been accelerated against SARS-CoV-2. However, considering the conventional timeframe for new drug development and alarming SARS-CoV-2 proliferation, drug repurposing has been established as better solution against SARS-CoV-2 as recommended by WHO (Kupferschmidt & Cohen, 2020) . As limited research is conducted to screen the known FDA approved drugs against SARS-CoV-2 M pro , therefore, we investigate the potential FDA approved drugs using structure-based virtual screening, all atoms molecular simulation, and post-simulation analysis against SARS-CoV-2 M pro , as shown in Figure 1 .

Structure-based virtual screening was performed against SARS-CoV-2 M pro (PDB ID: 6LU7) (Jin et al., 2020) to identify the potential inhibitors from Druglib database-containing FDA approved drugs using MTiopen screen web server (Labbe et al., 2015) . The viral protease contained three structural domains; domain I (Phe 8 -Tyr 101 ), domain II (Lys 102 -Pro 184 ), and domain III (Thr 201 -Val 303 ), interconnected by loop of residues Phe 185 to Ile 200 . The active pocket with catalytic dyad (Cys 145 and His 41 ) was defined in the cleft between domains I and II of SARS-CoV-2 M pro (Jin et al., 2020) . For virtual screening, SARS-CoV-2 M pro was prepared under default parameters using Dock prep tool in Chimera-1.14 (Pettersen et al., 2004) . Besides, residues (His 41 , Phe 140 , Gly 143 , Cys 145 , His 163 , His 164 , Glu 166 , Gln 189 , and Thr 190 ) interacting with native ligand N3 inhibitor in the crystal structure of viral protease were considered as active residues for virtual screening against Druglib database via MTiOpenScreen virtual screening web server (Labbe et al., 2015) .

Following, top 10 docked poses of protein-drug complexes were retrieved from virtual screening and re-docked again in the pocket covering same active residues using Chimera-1.14-AutoDock Vina plugin setup, as reported earlier Bharadwaj, Lee, Dwivedi, Yadava, et al., 2020) . Briefly, protein and screened drugs were minimized under default parameters in structure minimization tool in USCF Chimera-1.14 (Pettersen et al., 2004) and saved in PDB format. Following, these 3D structures were prepared with Dock prep tool in Chimera-1.14, where polar hydrogen atoms and charges were added. Later, molecular docking was conducted using AutoDock Vina (Trott & Olson, 2010) , as plugin in USCF Chimera-1.14 under default parameters on the active residues covered by docking grid box of size 60 Â 40 Â 40 Å along the three (X, Y and Z) axes while center at À8.918, 17.918, and 62.905 Å region. Finally, top poses with the highest docking score and least RMSD (0 value by default) were selected for intermolecular interaction profiling using ligand-receptor interaction module in academic version of Maestro version 12.3 (Schr€ odinger Release 2020-1: Maestro, Schr€ odinger, LLC, New York, NY, 2020). Herein, various intermolecular interactions between ligands and active residues of protein, i.e. hydrogen bonding, hydrophobic, p-cation, p-p interaction, interaction, negative, positive, glycine, polar, and salt bridges formation, were calculated at cutoff radius of 4 Å under default conditions. Likewise, same docking procedure was also employed for the selected reference compound, i.e. Michael acceptor inhibitor or N3 inhibitor, which was reported as native ligand in the crystal structure of SARS-CoV-2 M pro (Jin et al., 2020) , to justify the docking methodology and for relative analysis against screened FDA approved drugs. Both 3D and 2D interaction images were rendered using academic version of Maestro version 12.3 package (Schr€ odinger Release 2020-1:Maestro, Schr€ odinger, LLC, New York, NY, 2020).

Furthermore, re-docked complexes were also analyzed for binding free energy using Prime molecular mechanics (MM)/ generalized born surface area (GBSA) module in Schrodinger suite (Schr€ odinger Release 2019-2:Maestro, Schr€ odinger, LLC, New York, NY, 2019) and can be expressed using the following Equation (1):

where DG MMGBSA Bind depicted the binding free energy, DG Complex marked for the binding energy of docked receptor-ligand complex while DG Ligand and DG Receptor showed the energy of receptor and ligand, respectively.

The selected docked poses for each protein-ligand complexes were further subjected to 100 nanoseconds (ns) molecular dynamics (MD) simulation under Linux environment on HP Z280 workstation using academic version of Maestro-Desmond version 5.6 module of Schr€ odinger-Maestro version 11.8 suite (Schr€ odinger Release 2018-4, Schr€ odinger, LLC, New York, NY, 2018) (Bowers et al., 2006) . Briefly, protein was prepared using protein preparation wizard module of Schr€ odinger suite under default parameters followed by generation of orthorhombic grid box (10 Â 10 Â 10Å buffer) as simulation box using system building tool. The complete prepared system was immersed in Monte Carlo equilibrated periodic transferable intermolecular potential with 4 points (TIP4P) water bath and neutralized with suitable number of sodium counter ions; salt and ions placement were omitted at 20 Å resolution around ligand. Following, 0.002 ps time steps for anisotropic diagonal position scaling were selected to maintain constant pressure for MD simulation process. Moreover, temperature was gradually raised from 100 to 300 K coupled with 20 ps NPT reassembly at 1 atm pressure. Following, density of the complete simulation system was preserved at 1 g/cm 3 and MD calculation was conducted under default force field, i.e. optimized potentials for liquid simulations (OPLS)-2005, parameters in academic version of Desmond version 5.6-Maestro suite version 11.8. Finally, MD simulation was performed for each selected docked complex under similar conditions for 100 ns interval. Subsequently, MD simulation trajectory for each system was investigated using simulation interaction diagram tool of the Desmond version 5.6 module in Schr€ odinger-Maestro version 11.8 package (Bharadwaj, Rao, et al., 2020; Dwivedi et al., 2020 

Essential dynamics aided in the elucidation of correlated fluctuations in protein structure that are essentially required for the protein function. Hence, generated MD trajectory for each complex was considered for essential dynamics via principle components analysis (PCA) using Bio3d package (Grant et al., 2006) . Moreover, residual displacement in the docked protein structure during the simulation interval was also scrutinized for residue correlation coefficient by dynamic cross-correlation analysis in Bio3d package (Grant et al., 2006) . Both essential dynamics and DCCM analysis were conducted on all the Carbon-alpha (Ca) atoms of each complex extracted as 5000 frames from 100 ns MD trajectory and later superimposed to the initial docked pose to reduce the root mean square variances among the equivalent residues in the protein structure. These calculations were computed for each simulation trajectory under R environment (Team, 2013) with Bio3d package for simulation trajectory analysis (Grant et al., 2006) .

3.1. Virtual screening, re-docking and MM/ GBSA analysis Virtual screening is the use of high-performance computing to screen the large-small molecule databases for the identification of potential ligands against a specified drug target (Bharadwaj, Rao, et al., 2020; Dwivedi et al., 2016) . MTiOpenScreen webserver predicts the potential drug molecules based on the minimum binding energy and correspondingly generates three significant binding poses against each screened molecule (Labbe et al., 2015) . In this study, structure-based virtual screening of SARS-CoV-2 M pro against Druglib database yields a total of 1051 FDA approved drugs with significant docking score between À10.3 and À7.5 kcal/ mol (Table S1 ). Following, top 10 docked drugs were selected for re-docking analysis in AutoDock Vina (Figure 2 ) revealed improved binding score, as mentioned in Table 1 , and occupied the same site on in the protease active pocket ( Figure S1 ). Interestingly, all the selected drug molecules showed binding affinity greater than À8 kcal/mol against N3 inhibitor (7.6 kcal/mol), which is considerably significant for SARS-CoV-2 M pro docked ligands, as reported earlier ( Bharadwaj, Lee, Dwivedi, Yadava, et al., 2020; Elmezayen et al., 2020) .

The screened best poses for all the 10 drugs molecules docked with SARS-CoV-2 M pro were analyzed for intermolecular interaction which assisted in the formation of stable docked complex. All the selected molecules showed significant placement in the binding pocket of SARS-CoV-2 M pro ( Figure 3 ) and shared some common residues for intermolecular interaction (Table 1, Figure S2 ). Since, hydrogen bonding in the docked complex marked for stability of the ligand with receptor, a maximum and a minimum of four and one hydrogen bonds, respectively were observed in the respective docked complexes; in addition to other intermolecular interactions (Table 1, Figure S2 ). Based on the number of hydrogen bond formation, selected drugs can be placed in the descending order; R428, UK-432097, Etoposide, Teniposide, Orvepitant, Vs5584, MK-3207, Tadalafil, CEP-11981, and Setileuton, as potent inhibitor of SARS-CoV-2 M pro . Moreover, all the drugs showed polar and hydrophobic interactions with the catalytic dyad residues His 41 and Cys 145 , respectively, except R428 and Etoposide also formed p-p stacking with His 41 while Teniposide exhibited hydrogen bond formation with Cys 145 residue. Moreover, these redocked poses were also studied in comparison to the best pose of SARS-CoV-2 M pro -N3 inhibitor for intermolecular interactions profiling which showed occupancy of similar residues (Table 1, Figure S3 ). It is important to mention that SARS-CoV-2 M pro -N3 inhibitor complex showed intermolecular molecular interaction via formation of four hydrogen bonds with residues Cys 145 , Glu 166 , and Gln 189 (2) along with other intermolecular interactions ( However, additional residue contacts were observed between the screened drugs and viral protease M pro (Table 1 , Figures  S2 and S3) . These drug molecules, except R428, MK-3207, Orvepitant, and UK-432097, also interacted with Asp 187 residue which was suggested to enhance the catalytic efficiency of viral M pro (Zhao et al., 2008) . The molecular contacts of screened drugs with the essential residues in the protease structure suggested the inhibition of SARS-CoV-2 M pro proteolytic function, which is essentially required for the viral replication and pathogenesis. Hence, comparative molecular docking analysis of screened FDA approved drugs against N3 inhibitor suggested the potential of selected drugs to inhibit SARS-CoV-2 M pro by formation of hydrogen and non-covalent interaction with its catalytic dyad and substrate binding residues.

To further established the role of intermolecular interactions between screened drugs and SARS-CoV-2 M pro , binding free energy along with contributing energy components were assessed using Prime MM/GBSA method for each docked complex (Table S2 , Figures 4 and S4) . These results suggested the maximum contribution of DG Bind Coulomb and DG Bind vdW in the stability of the respective complexes while DG Bind Covalent and DG Bind Packing contribute in the destabilization of the docked respective complexes. Remarkably, SARS-CoV-2 M pro -Etoposide, SARS-CoV-2 M pro -Cept-11981, SARS-CoV-2 M pro -Teniposide, SARS-CoV-2 M pro -Orvepitant, and SARS-CoV-2 M pro -Setileuton docked complexes exhibited relatively higher binding free energy (>À60 kcal/mol) against other docked drug complexes and N3 Inhibitor complex (À84.75 kcal/mol). These observations further support the potential of screened drugs for the inhibition of SARS-CoV-2 M pro .

MD simulation is a widely accepted computational method in drug discovery in order to understand the physical interactions, such as structure-function relationships, intramolecular/ intermolecular interactions, and other structural properties, at the atomic-level for the biological macromolecules (Bharadwaj, Lee, Dwivedi, Yadava, et al., 2020; . Thereof, the selected SARS-CoV-2 M pro -screened drug complexes were monitored for the docked complex stability under 100 ns MD simulations conducted in Desmond suite. The produced MD trajectory for each docked complex were scrutinized in terms of last snapshot of simulated docked complex, root mean square deviation (RMSD), root mean square-fluctuations (RMSF), and protein-ligand contacts mapping with respect to 100 ns interval, as reported earlier ( Bharadwaj et al., 2019) .

The last snapshot from each simulated complex trajectory was analyzed for the occupancy in the active pocket and formation of molecular contacts with viral protease by comparison to N3 inhibitor as reference ligand (Figures 5 and S5) . Interestingly, all the drugs were relatively noted within the catalytic pocket of SARS-CoV-2 M pro and observed for formation of substantial interactions, including hydrogen bonds and other intermolecular interactions against N3 inhibitor Teniposide, VS-5584, and Setileuton with viral protease based on formation of hydrogen bonds with essential residues against N3 inhibitor, as noticed in the respective initial docked poses (Table 1) . Hence, these observations indicate the considerable stability of these protein-drug complexes against other docked drugs and reference ligand N3 inhibitor with SARS-CoV-2 M pro . Moreover, an average change in the displacement of Ca atoms (RMSD) of viral protease and protein fit drug molecule as ligand in reference to their corresponding initial frames were analyzed for each SARS-CoV-2 M pro -drug complex by comparison to SARS-CoV-2 M pro -N3 inhibitor complex (Figures 6 and S7) . Remarkably, Ca atoms in each docked complex showed acceptable deviation (<2.5 Å) during 100 ns simulation interval (Figures 6 and S7) . These variations in Ca atoms were further supported by RMSF of respective receptor in docked complexes, where residues showed deviations <2.5 Å except in N-and C-terminal of the viral protease (<4.5 Å) (Figures S8 and S9 ). As N-and C-terminal of the protein structure does not contribute substantially to the activity of the viral protease, these residual fluctuations in each complex are considered as acceptable during the molecular dynamic's simulation. Likewise, RMSD values for protein fit drugs were extracted and analyzed from 100 ns simulation interval ( Figure 6) . Although, occasional fluctuations <6.0 Å were recorded for some of protein fit drugs, i.e. R428, Etoposide, MK-3207, Teniposide, VS-5584, Tadalafil, and Setileuton, but these docked ligands in the active pocket of (Figure 6 ). However, N3 inhibitor docked with viral protease initially observed with considerable RMSD values followed by higher deviations (<18.0 Å) after 80 ns till end of MD simulation (Figures S7), suggested that N3 inhibitor required more than 100 ns interval to achieve the equilibrium with viral protease. Besides, RMSF analysis for all the screened drug atoms showed acceptable variations (<3.5 Å) except in the regions with heavy atoms (<10.5 Å) by comparison to N3 inhibitor ( Figures S10 and S11) .

Furthermore, protein-ligand intermolecular interactions for the screened drugs against N3 inhibitor were also extracted during the course of 100 ns MD simulation, including hydrogen bonds, hydrophobic contacts, ionic interactions, and water bridge formations, which reflect the binding strength and permanency of docked ligands with viral protease (Figures 7 and S12) . Interestingly, various residues were recorded for the substantial interactions with the screened drugs by comparison to N3 inhibitor during 100 ns simulation (Figures 7 and S12 ). For instance, screened drugs showed one or more than one type of interaction with active residues, viz. catalytic dyad (His 41 and Cys 145 ), substrate binding residues (Phe 140 , Leu 141 , Gly 143 , Ser 144 , His 163 , Met 165 , and Glu 166 ), and other essential substrate-binding residues (Thr 24 , Thr 25 , Met 49 , Phe 140 , Asn 142 , His 163 , Met 165 , Asp 187 , and Gln 189 ) (Figure 7) , as defined in the crystal structure of SARS-CoV-2 M pro with N3 inhibitor (Jin et al., 2020) . Remarkably, these residues were also logged in the respective docked complexes of the protease with screened drugs (Table 1) , further support the stability of selected drugs with active pocket of SARS-CoV-2 M pro . Also, protein-ligand interactions profile extracted at 30% of simulation interval for each MD trajectory demonstrates the highest number of non-covalent interactions for all the docked drugs, except Orvepitant, with catalytic dyad residues and other essential residues in the active pocket of SARS-CoV-2 M pro against N3 inhibitor (Figures S13 and S14).

Conclusively, only drugs R428, Teniposide, VS-5584, and Setileuton against N3 inhibitor docked with SARS-CoV-2 M pro were noted for significant molecular contacts via hydrogen bonding with viral protease in the extracted last pose from MD trajectories, acceptable variations values, and substantial interactions with essential residues of viral protease. Hence, based on MD simulation analysis, R428, Teniposide, VS-5584, and Setileuton drugs are concluded as potential inhibitors of SARS-CoV-2 M pro than other drug molecules and N3 inhibitor.

To compute the influence of MD simulation on the binding free energy of ligands with the active pocket of SARS-CoV-2 M pro , snapshots collected from respective MD simulations were analyzed for binding affinity using Prime MM/GBSA method (Figure 8 ). Total DG Bind and individual energy components values computed on the respective snapshots from simulation trajectories are given Table S4 . All the respective snapshots for SARS-CoV-2 M pro -FDA approved drugs exhibited no considerable change in the net binding energy, except reduction in binding affinity was noted for Etoposide, Orvepitant, Setileuton, and N3 inhibitor docked with viral protease (Figures 8 and S15) . Additionally, analysis of individual components contributing in the total binding free energy of the respective complexes, i.e. DG Bind Coulomb , DG Bind Covalent , DG Bind Hbond , DG Bind Lipo (hydrophobic interactions), DG Bind Solv GB (Generalized Born electrostatic solvation energy), DG Bind packing (p-p interactions), and DG Bind vdW (van der Waals energy) revealed significant contribution of DG Bind Coulomb and DG Bind vdW energy in the respective complex stability (Table S4 , Figures 8 and S15 ), supports the favorable enthalpy formation in the respective complexes. Similar results for the energy components DG Bind Coulomb and DG Bind vdW contributing in the docked complexes were reported earlier ( Bharadwaj et al., 2019) . Moreover, no significant contribution of ligand strain energy was observed in all the complexes, indicates the considerable affinity of screened drugs by comparison to N3 inhibitor with the active pocket of SARS-CoV-2 M pro .

Essential dynamics were performed on the MD trajectories to collect the major Eigen values, also known as principal component analysis ( Figures S16 and S17) . These values suggest a major level of conformational mobility caused by docked ligand in active pocket of viral protease. However, successive elbow point and no notable change in variations of the Eigen fraction were detected following 4 th Eigen value (Figures S16 and S17). These results intimated that substantial flexibility exist in SARS-CoV-2 M pro docked with selected compounds during initial phase of MD simulation that diminished with respect to time. Also, slow decrement in proportional contribution of the Eigen modes advised the additional localized fluctuations in SARS-CoV-2 M pro docked with each ligand to obtain the favorable stability. Hence, these variations in each complex were advised to play a key role in the stability of respective docked complexes. Furthermore, first two Eigen vectors for SARS-CoV-2 M pro docked with each compound that were obtained from respective MD trajectory as cluster groups, indicated compact and cluster motions for SARS-CoV-2 M pro during 100 ns simulation, except in SARS-CoV-2 M pro -Etoposide complex and SARS-CoV-2 M pro -MK-3207 complex (Figures 9 and S17) . Also, the generated plots suggested the variations in cluster distribution in each conformation during MD simulation via color gradient change from blue to red stands for periodic jumps among the various conformational poses of docked viral protease. In conclusion, a restricted correlated fluctuations and motions of the viral protease in all the studied systems represent the rigidity and stability of respective docked complexes, except in SARS-CoV-2 M pro -Etoposide complex and SARS-CoV-2 M pro -MK-3207 complex, during MD simulation. Furthermore, to calculate the conformational dynamics changes instituted in SARS-CoV-2 M pro structure in each complex due to inhibitory activity of docked ligands, DCCM analysis was computed based on the positions of Ca atoms. Figure 7 shows high correlated motions which varies from the light blue to cyan color (þ1) while anti-correlated motions are in the range from light purple to red brick color (À1). Analysis of residue cross correlation suggested a considerable correlated motions and dynamic changes in all the systems by comparison to N7 inhibitor (Figures 10 and S17) . The calculated results established that screened drugs and N3 inhibitors significantly induced the restricted conformation changes in the viral protease structure during 100 ns MD simulation (Figure 10) .

Based on the combinatorial computational analysis, including structure-based virtual screening, molecular docking, binding free energy calculations, MD simulation and post-MD simulation analysis for the screened FDA drugs with viral protease, suggested the drugs, R428, Teniposide, VS-5584, and Setileuton with comparatively higher stability and affinity with SARS-CoV-2 M pro against N3 inhibitor via strong intermolecular interactions formation as well disturbing the conformation of viral protease active pocket.

This study was aimed to repurpose the FDA approved drugs as potential inhibitor of SARS-CoV-2 M pro using structure-based virtual screening approach and validation by MD simulation. Initially, 1051 compounds were obtained from structure-based virtual screening of Druglib database and only top 10 drug molecules were selected as potential ligands of SARS-CoV-2 M pro for further intermolecular and complex stability. The selected drugs showed significant interactions with the catalytic (His 41 and Cys 145 ) and other substrate-binding residues in docked complexes as well as during 100 ns MD simulation. Thereof, based on collective observations, R428, Teniposide, VS-5584, and Setileuton drugs were suggested to hold the potential to significantly occupy the active pocket of SARS-CoV-2 M pro and can be used as its inhibitors. Hence, computational screening and validation of top 10 drugs recorded in this study can be further evaluated against in vitro SARS-CoV-2 M pro inhibition and viral infection for the drug formulation against SARS-CoV-2 infection.

Computational insights into tetracyclines as inhibitors against SARS-CoV-2 M-pro via combinatorial molecular simulation calculations

Density functional theory and molecular dynamics simulation support Ganoderma lucidum triterpenoids as broad range antagonist of matrix metalloproteinases

Discovery of Ganoderma lucidum triterpenoids as potential inhibitors against Dengue virus NS2B-NS3 protease

Structure-based screening and validation of bioactive compounds as Zika virus Methyltransferase (MTase) inhibitors through first-principle density functional theory, classical molecular simulation and QM/MM affinity estimation

Scalable algorithms for molecular dynamics simulations on commodity clusters

Two-year prospective study of single infections and co-infections by respiratory syncytial virus and viruses identified recently in infants with acute respiratory disease

Recent insights into the development of therapeutics against coronavirus diseases by targeting N protein

Emerging coronaviruses: Genome structure, replication, and pathogenesis

Chloroquine and hydroxychloroquine as available weapons to fight COVID-19

Origin and evolution of pathogenic coronaviruses

Middle East respiratory syndrome coronavirus (MERS-CoV): Announcement of the coronavirus study group

Discovering drugs to treat coronavirus disease 2019 (COVID-19)

Outbreak of a novel coronavirus

Anti-dengue infectivity evaluation of bioflavonoid from Azadirachta indica by dengue virus serine protease inhibition

Identification of new potent inhibitors of dengue virus NS3 protease from traditional Chinese medicine database

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

Aetiology: Koch's postulates fulfilled for SARS virus

Novel coronavirus COVID-19: An overview for emergency clinicians

Bio3d: An R package for the comparative analysis of protein structures

National ambulatory medical care survey: 2004 summary

Coronavirus infections in hospitalized pediatric patients with acute respiratory tract disease

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

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

From SARS to SARS-CoV-2, insights on structure, pathogenicity and immunity aspects of pandemic human coronaviruses

Race to find COVID-19 treatments accelerates

MTiOpenScreen: A web server for structure-based virtual screening

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

Study of the differential activity of thrombin inhibitors using docking, QSAR, molecular dynamics, and MM-GBSA

UCSF Chimera-a visualization system for exploratory research and analysis

R: A language and environment for statistical computing

Software news and update auto dock vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

Detection of the new human coronavirus HKU1: A report of 6 cases

Identification of a new human coronavirus

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

The crystal structures of severe acute respiratory syndrome virus main protease and its complex with an inhibitor

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

Structure of the main protease from a global infectious human coronavirus, HCoV-HKU1

A novel coronavirus from patients with Pneumonia in China

Authors are highly thankful to Dr. Amaresh Kumar Sahoo, Indian Institute of Information Technology, Prayagraj, India for providing his kind support in binding free energy calculation using Prime MMGBSA module of Schrodinger Suite 2019.2.

The authors declare no competing interests.

This study was supported by The Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, Saudi Arabia, under grant no. FP-3-42.

http://orcid.org/0000-0003-2832-8617 Umesh Yadava http://orcid.org/0000-0002-9127-532X Sang Gu Kang http://orcid.org/0000-0002-5216-5494 Vivek Dhar Dwivedi http://orcid.org/0000-0003-0767-1328