key: cord-0743628-2832gi7f authors: Tachoua, Wafa; Mohamed, Kabrine; Mushtaq, Mamona; Ulhaq, Zaheer title: An in-silico evaluation of COVID-19 main protease with clinically approved drugs date: 2020-09-21 journal: J Mol Graph Model DOI: 10.1016/j.jmgm.2020.107758 sha: 798b4341bd871aaaad19ad0beca4cb5f3c639bef doc_id: 743628 cord_uid: 2832gi7f A novel strain of coronavirus, namely, SARS-CoV-2 identified in Wuhan city of China in December 2019, continues to spread at a rapid rate worldwide. There are no specific therapies available and investigations regarding the treatment of this disease are still lacking. In order to identify a novel potent inhibitor, we performed blind docking studies on the main virus protease M(pro) with eight approved drugs belonging to four pharmacological classes such as: anti-malarial, anti-bacterial, anti-infective and anti-histamine. Among the eight studied compounds, Lymecycline and Mizolastine appear as potential inhibitors of this protease. When docked against M(pro) crystal structure, these two compounds revealed a minimum binding energy of -8.87 and -8.71 kcal/mol with 168 and 256 binding modes detected in the binding substrate pocket, respectively. Further, to study the interaction mechanism and conformational dynamics of protein-ligand complexes, Molecular dynamic simulation and MM/PBSA binding free calculations were performed. Our results showed that both Lymecycline and Mizolastine bind in the active site. And exhibited good binding affinities towards target protein. Moreover, the ADMET analysis also indicated drug-likeness properties. Thus it is suggested that the identified compounds can inhibit Chymotrypsin-like protease (3CL(pro)) of SARS-CoV-2. The outbreak of the novel coronavirus disease, COVID-19, caused by coronavirus 2019-nCoV officially designated as severe acute respiratory syndrome related coronavirus 2 (SARS-CoV-2), represents a pandemic threat to global public health [1, 2] . On January 30, the World Health Organization (WHO) announced a Public Health Emergency of International Concern (PHEIC) for the 2019-nCoV outbreak [3] . Coronaviruses are relatively large viruses containing a single stranded positive-sense RNA genome encapsulated within a membrane envelope. The viral membrane is studded with glycoprotein spikes that give coronaviruses their crown-like appearance. There are four classes of coronaviruses, namely alpha, beta, gamma and delta. The beta-coronavirus class include severe acute respiratory syndrome virus (SARS-CoV), the middle east respiratory syndrome virus (MERS-CoV), and the newly identified severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [4] . Although SARS-CoV-2 is classified into beta-coronaviruses group, it is diverse from MERS-CoV and SARS-CoV. It has been reported that SARS-CoV-2 genes share less than 80% identity (nucleotides) with SARS-CoV and it is more transmissible than other SARS-CoV viruses [5] [6] [7] . The SARS-CoV-2 genome consists of approximately 30000 nucleotides. It encodes several structural proteins such as the glycosylated spike protein (S), envelope protein (E), membrane protein (M), and nucleocapsid protein (N). Additionally, the viral genome also encodes numerous nonstructural proteins, including RNA-dependant RNA polymerase (RdRp), coronavirus main protease (M pro ), and papain-like protease (PLpro). Upon entrance in the host cell, the viral genome is released and subsequently translated into viral poly-proteins using host cell translation machinery [1, 6] . The poly-proteins are then cleaved into effector proteins by viral proteases PLpro and M pro [8, 9] . The M pro , also known as 3-chymotrypsin-like protease (3CL), plays a critical role in the virus replication process. It cleaves pp1a and pp1b polyproteins subsequently releasing functional proteins including RNA polymerase, endoribonuclease and exoribonuclease. Therefore, M pro is a potential target for anti-coronaviruses screening. Indeed, inhibiting M pro activity could stop the spread of infection [10, 11] . The released crystal structure of M pro (6lu7) was obtained by a co-crystallization with a peptidelike inhibitor called N3 (PRD_002214). The enzyme has a molecular weight of 33.79 kDa as determined by Mass spectroscopy. M pro forms a dimer, where each monomer comprises three domains: Domain I (residues 8-101) and II (102-184) consists of an antiparallel beta barrel, and the alpha helical domain III (residues 201-301) is required for enzymatic activity. It shares a similar structure with cysteine protease, however the active site lacks the third catalytic residue. Thus, catalytic harbors Histidine 41 (H41) and Cysteine 145 (C145) [12, 13] located in between domain I and II, while amino acids, T24, L27, H41, F140, C145, H163, M165, P168 and H172 form a hydrophobic surrounding in the pocket [14] . The concept of drug repurposing is widely used nowadays to identify potential drugs for COVID-19 disease. The practice of reusing drugs for more than one purpose can significantly reduce the cost, time and risks of the drug development process [15] . To date many studies have targeted SARS-CoV-2 M pro , out of which, the study of Khan and collaborators, based on computational drug design methods, yielded five promising hits corresponding to three antiviral drugs (Remdesivir, Saquinavir and Darunavir) and two natural compounds (Flavone and coumarine) [16] . Similarly, Aanouz and collaborators proposed three compounds (Corcin, Digitoxigenin, and β-Eudesmol) as potential inhibitors of M pro from Moroccan medicinal plants [17] . Interestingly, Talampicillin and Lurasidone are reported to possess the highest affinity towards M pro enzyme and might be repurposed against COVID-19 [11] . Although, several molecular docking studies have been established to find a potential inhibitor of M pro activity based on antiviral compounds commonly used to treat immunodeficiency virus (HIV), phytochemical, antimalarial agents, spices or marine products [18] [19] [20] [21] [22] . Unfortunately, no specific therapies for COVID-19 disease are available and investigations regarding the treatment of COVID-19 disease are still lacking. Therefore, there is a dire need to identify approved drugs which may inhibit SARS-CoV-2 virus proteins and to find at the same time their optimal association and concentration to treat COVID-19 patients. In the following study, we investigated the M pro inhibitory potential of eight clinically approved drugs belonging to four pharmacological classes: anti-viral, anti-bacterial, anti-infective and antihistamine, using a molecular docking and a molecular dynamics approach. This research is a descriptive-analytical study. In this study, the interaction of several approved compounds was investigated. A total of eight compounds were tested against COVID-19 main protease (M pro ). N3 compound was used as a docking target for comparison. In order to obtain the structure information of selected compounds, a Drugbank database (https://www.drugbank.ca/) was used [23] (Table 1 ). The Three dimensional structure of M pro ( Fig.1 ) was retrieved from the Protein Data Bank (PDB) (https://www.rcsb.org/) [24] . It corresponds to a complex between the enzyme and its inhibitor N3 [13] . The 6lu7 structure preparation consists of several steps such as deleting all water molecules, N3 inhibitor and adding hydrogens. The new file was saved for docking analysis. In order to investigate the molecular interaction between several approved compounds and COVID-19 main protease (M pro ), blind docking was performed using a SwissDock server (http://www.swissdock.ch/) under the accurate mode with no flexibility of the side chain of any amino acid of the target protein. In addition, a binding pocket was not defined so as not to bias the docking towards the active site. The protein and ligand were specified by uploading PDB and Mol2 files, respectively. SwissDock generates all possible binding modes for each ligand and the most favorable binding modes at a given pocket were clustered. All ligand clusters were saved in an output file called "prediction file". The prediction file provided; Cluster rank, Element, Fulfitness and estimated binding free energy (ΔG). A cluster corresponds to a predicted binding pocket on the target protein and the cluster rank represents the different conformations of a ligand in a certain J o u r n a l P r e -p r o o f cluster [25] . Only the lowest energy model of cluster zero was considered to be the most favorable interaction. After docking Chimera and Pymol softwares were used to visualize the receptor ligand interactions for the lowest energy model of the clusters obtained from the previous step. Each ligand cluster was inspected for amino acids interacting with the ligand, hydrogen bonds (Hbonds), the specific atoms involved. All the interacting amino acids with the target were noted for each cluster [26, 27] . 2D diagrams were generated by PoseView integrated into Proteins Plus server (https://proteinsplus.zbh.uni-hamburg.de/#dogsite) [28] , which automatically creates twodimensional diagrams of protein ligand complexes according to chemical drawing conventions. The generation of structure diagrams and their layout modifications are based on the library 2Ddraw. Interactions between the molecules are estimated by a builtin interaction model that is based on atom types and simple geometric criteria [29, 30] . Although the binding site is well characterized for N3 within many CoV M pro crystals. We have applied the DoGSiteScorer online tool (https://proteinsplus.zbh.uni-hamburg.de/#dogsite) to predict and describe binding pockets within native M pro and the complexes M pro /inhibitors obtained after docking analysis. DoGSiteScorer is a grid-based method which uses a Difference of Gaussian filter to detect potential binding pockets solely based on the 3D structure of the protein and splits them into subpockets. Global properties, describing the size, shape and chemical features of the predicted (sub) pockets are calculated. Per default, a simple druggability score is provided for each (sub) pocket, based on a linear combination of the three descriptors describing volume, surface, hydrophobicity and enclosure. Furthermore, a subset of meaningful descriptors is incorporated in a support vector machine (libsvm) to predict the (sub) pocket druggability score (values are between zero and one). The binding pockets are ranked according to their size, surface area and druggability score [31] . To analyze drug likeliness and the pharmacokinetics parameters, all the molecules were subjected to ADMET (Absorption, distribution, metabolism, excretion and toxicity) predictions with SwissADME (http://www.swissadme.ch/) and PROTOX web servers (http://tox.charite.de/protox_II/) [32, 33] . PROTOX predict the median oral lethal dose, LD 50 (Lethal dose, 50 %) based on the analysis of the two-dimensional similarity to compounds with known LD 50 values and the identification of fragments over-represented in toxic compounds [33] . The molecular dynamics simulation (MD) was performed on the best two compounds obtained from molecular docking study (Lymecycline and Mizolastine), which helped to get more insight into the protein and docked complexes in biological conditions. To gain atomic level insight into the binding interaction of Lymecycline and Mizolastine with SARS-Cov2 M pro , all atom-molecular dynamics simulations were carried out using PMEMD.cuda module implemented in AMBER 18 with explicit water model under periodic boundary conditions [34] . The parameters for both systems (M pro -Lymecycline and M pro -Mizolastine) were prepared using the antechamber and tleap modules. Amber ff14SB [35] force field was assigned to protein whilst Amber Generalised Force Field (GAFF) [36] was selected for the two compounds. Each system was solvated using the TIP3P water model in a periodic box, with a minimum spacing of 10 Å from the solute/protein atoms. After that neutralization of the system with counter ions (Na + ) replacing solvent molecules at the position of electrostatic favorable potential. The Particle Mesh Ewald Method (PMEM) and the SHAKE algorithm were used to calculate the long-range electrostatic interactions and to constrain the hydrogen bonds respectively [37, 38] . To remove bad contacts and possible steric clashes both complexes were subjected to energy minimization. Briefly systems were relaxed by adjusting hydrogen position, whereas for the Mizolastine complex the first 2500 cycles were carried out using steepest descent algorithm and the next 2500 steps involved the conjugate gradient algorithm. However, for Lymecycline complex the number of cycles were increased to achieve minimized system. The minimization steps were repeated multiple times with gradual decrease in force applied to restrain heavy atoms. It was followed by 5000 (50000 for lymecycline complex) additional minimization steps in absence of any restraint. Afterwards each system was equilibrated under NVT (gradual heating from 0 K to 300 K ) and NPT (pressure = 1 atm) ensembles. The final production run of 120 ns was performed at constant temperature (300 K) and pressure (1 atm) with integration time steps of 2 fs. The obtained results were analyzed via CPPTRAJ modules. Chimera and VMD were used for visualization and the graphs were plotted using Xmgrace tool [27, 39] . To estimate the binding free energies of the selected inhibitors against SARS-CoV-2 M pro , the widely known MM-PBSA was performed. This tool computes components of binding free energy utilizing the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA), implemented in AMBER 18. A total of 2500 frames were extracted from the most stable trajectories and of them 500 frames were used to compute various energy components involved in ligand binding. In the MMPBSA protocol, the following equation is used to calculate binding free energy (ΔG bind ) of an inhibitor. Where ΔE MM represents the alterations in molecular mechanics energy at gas phase, ΔG solv and -TΔS indicate desolvation free energy and the conformational entropic contributions upon binding with ligand, respectively. Further, ΔE MM is sum of ΔE internal (bond, dihedral and angle energies), ΔE elec (electrostatic), and ΔE vdW (Van der Waals) energies while ΔG pol and ΔG np represents polar and nonpolar solvation energy [40] . In order to investigate the possible mechanism by which selected drugs act, an silico theoretical molecular docking approach was used. During our study, we simulated the binding mode of N3 against 6lu7 crystal structure using SwissDock to ensure the effectiveness of docking results and to compare results produced by several drugs to those of N3. Indeed, this compound is a well characterized inhibitor of COVID-19 main protease. Docking results revealed that N3, Indinavir and Chloroquine had the best energies of binding -10.83, -9.81 and -9.71 kcal/mol, respectively ( Table 2 , column 5), which is consistent with three studies. The first one reported the entire complex N3/M pro crystal structure saved in the PDB database under 6lu7 accession number [13] . The second reported that Indinavir exhibited a good docking score (-7.05) when docked against 5r7z M pro structure using flexible docking with Glide and the last one revealed that Chloroquine and its derivatives can bind to M pro [18. 21] . According to a full fitness score, Lymecycline and Mizolastine had more favorable binding mode, indicated by more negative fullfitness scores -1332.56 and -1300.12 kcal/mol, respectively ( Table 2 , column 4). Azelastine can undergo an optimal binding with M pro [41] . When the clusters were analyzed it was found that N3, Indinavir, Chloroquine and Mizolastine showed that all clusters were able to fit into M pro binding pocket. However, Lymecycline revealed a binding energy of -8.87 kcal/mol and it was able to occupy 23 clusters constituting a total of 168 possible conformations within the substrate binding cavity, out of a total of 256 elements ( Table 2 , column 3). It is striking that all or more than half of the total predicted elements are docked in the substrate binding pocket. To investigate the possible reasons for differences in the binding energies, we examined the docked complexes with Pymol software and PoseView integrated in Protein Plus server. Table 2 showed the number and length of H-bonds formed between the target protease and the different compounds. Chloroquine, Nitazoxanide and Cetirizine established only one H-bond with N142, E166 and N142 residues, respectively. Otherwise, Quinine, Doxycycline and Indinavir were found to form two H-bonds with E166, E166 and (L141, Gly143) residues, respectively. Interestingly, Mizolastine and Lymecycline were found to form three H-bonds with M pro . Indeed, Lymecycline established three specific H-bounds (two H-bonds with E166 residue and one H-bond with Q189 residue) (Fig. 3A, and C) . However, Mizolastine formed three H-bonds (two H-bonds with T24 J o u r n a l P r e -p r o o f residue and one with G143) and hydrophobic interactions with N142 and C145 residues (Fig. 3B and D). N3 compound forms multiple hydrogen bonds with the main chain of residues in the substratebinding pocket [31] . However, only two hydrogen bonds were detected after SwissDock analysis because the complex of M pro with its inhibitor N3 was obtained in theoretical (in silico) not in experimental conditions. Although Lymecycline and Doxycycline belong to the same family of tetracyclines, Lymecycline bind more effectively to M pro with a minimum energy of -8.87 kcal/mol compared to Doxycycline with -7.52 kcal/mol binding energy, suggesting, the role of NH (CH2) 4 CH COOH NH2 chemical substituting group in increasing the binding affinity of Lymecycline towards M pro enzyme. To elucidate and describe binding pockets within M pro , Dogsitescorer server was used to analyze unliganded structure (6lu7) and N3/6lu7, Lymecycline/6lu7, Mizolastine/6lu7 complexes. Here, we have reported only the first three pockets in Table 3 . Results revealed that N3, Lymecycline and Mizolastione occupied the same pocket (P0) within 6lu7 structure (Fig. 4) with a high druggability score of 0.8 and a volume of 1191.74, 1061.18 and 1266.18 Å 3 , respectively. Based on Docking and Dogsitescorer studies, it is evident that Lymecycline and Mizolastine showed favorable binding with the new M pro , and the data comparable with those of N3. Drug likeness and pharmacokinetics proprieties analysis showed that Mizolastine and Chloroquine had four interesting features; a high lipophilicity, bioavailability score of 0.55, and higher gastrointestinal (GI) absorption. Otherwise, Lymecycline is soluble compared to Chloroquine, Indinavir and N3. The absorption is low and presents a bioavailability score of ≈0.1 with N3 compound. Interestingly, Lymecycline and Chloroquine are not substrates of P-glycoprotein compared to other drugs (Table 4) . Indeed, P-glycoprotein is a multidrug transporter which is apically expressed in the gastrointestinal tract, liver, kidney and brain endothelium. Consequently, P-glycoprotein plays an important role in the oral bioavailability, CNS (Central Nervous System) distribution biliary and renal elimination of drugs which are substrates of this transporter. It has been reported that Indinavir and Saquinavir displayed a reduced antiviral activity. Both drugs partially recover the ability to inhibit the replication of HIV-1 in the presence of Verapamil (Pglycoprotein Inhibitor). Thus, P-gp expression may affect the activity of antiviral drugs [42, 43] . Toxicity data analysis showed that Lymecycline (LD 50 =3000 mg/kg) and Mizolastine (LD 50 =450 mg/kg) are safe to used compared to Chloroquine (LD 50 =311 mg/kg). In addition, several studies reported that Tetracyclines had clinical utility against intracellular parasites such as malaria, Chlamydia and a-proteobacteria [44] . J o u r n a l P r e -p r o o f Molecular dynamic is a state-of-the-art simulation method for studying the physical motion and trajectory of the atoms in the presence of other molecules along with the various interactions within a system. It helps to follow and understand the structural features and conformational dynamics in the system. Thus, to validate the stability of the system and to probe ligand induced perturbations, MD simulation was performed with two best compounds as a function of time. The MD trajectories were examined based on various parameters including Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (Rg), Intermolecular hydrogen bond interaction and occupancy. Moreover, binding free energy calculations were also performed. RMSD monitors the deviations in average distance between the atoms of target protein during simulation with respect to initial docking structure/ reference frame. In short it is the deviation of the 3D structure over time. It provides insight into the system's stability, equilibrium and convergence whereas, the smaller fluctuations and constant backbone atoms (C, Cα, N, and O) RMSD, is indicative of the stable system. As described in Fig. 5A after an initial period of fluctuation both systems attained equilibrium during the last 50 ns of the simulation run. In general, the M pro and Lymecycline system displayed slightly higher fluctuation, whilst in comparison the lowest deviations were observed for M pro -Mizolastine complex. For Lymecycline complex during the initial frames continuous increase in RMSD value was observed in the range of < 2 -4Å however, in the last 50 ns trajectories the system obtained stability with the deviation of < 3 Å whereas no sharp fluctuations were observed during this time frame. In comparison, for Mizolastine complex, after gradual increase in fluctuation during the initial 45 ns time period, the system attained equilibrium state in the remainder MD trajectories except the frames in between 60-65ns where sharp fluctuation peaks yet in acceptable range (< 3.8 Å) were observed. The average RMSD for M pro -Lymecycline and M pro -Mizolastine complex was maintained at 3.10± 0.43 and 3.66± 1.77 Å, which implies that both systems attained a more stable structure compared to initial structure. Additionally, there was not much deviation between average and observed RMSD of protein at the end of the120 ns simulation and for both systems the RMSD during the whole run was < 4 Å which is in acceptable range. RMSF analysis is highly suitable to calculate the time dependent fluctuations in each residue, in turn is an essential parameter to determine the protein's flexibility. The RMSF of backbone atoms was calculated for M pro with 306 residues with the two potential inhibitors. As illustrated in Fig. 5B , both complexes had almost similar fluctuations trend. For most of the protein residues the observed RMSF value was lower than 2.5 Å. Whereas, the residues interacting with the ligands in the active site were found stable and displayed little fluctuations over time indicating the stable nature of compounds with target protein. In contrast the main fluctuations corresponded to the region that were distant from the ligand active site and others were found around the flexible loop region. Furthermore, the RMSF fluctuations in the Mizolastine complex was observed to be lower than the Lymecycline complex, suggesting that it had relatively lower structural mobility than the Lymecycline. The average RMSF values were 1.46± 0.93 and 1.32± 0.57 for Lymecycline and Mizolastine respectively. The Rg is a parameter to assess the compactness and overall dimension of the protein which in turn signifies folding and unfolding of the protein. The lower the gyration values the more folded the protein is and vice versa. Therefore, Rg was calculated to determine whether the tested drugs maintained the folding of the system. The gyration graph for backbone carbon atoms of the protein relative to time is presented in the Fig. 5C . The Rg values for both the systems range from 22-22.5 Å with an average score of 22.31± 0.142 Å for Lymecycline and 22.43± 0.28 Å for Mizolastine. As evident from the plot, the Lymecycline complex attained more compacted form as the MD simulation progressed, indicating a well converged system. Free energy calculations are a collection of methods used to envisage ligand binding affinities by considering several atomic level interactions that could be responsible for the affinity of ligand toward targeted protein. Thus, to understand the biophysical basis of interaction between Lymecycline and Mizolastine with main protease in further detail MM/PBSA calculations were carried out. For each complex non-polar solvation energy (G nonpolar ), polar solvation (G PB ), electrostatic (E elec ), Van der waals (E vdw ) and binding free energy (∆G binding ) were calculated. A summary of individual components involved in the binding free energy of tested inhibitors with main protease is shown graphically in Fig. 6 with the data listed in Table 5 . It is evident from the data that the binding of Lymecycline to the target protein is favored mostly by intermolecular electrostatic interactions which is also evident from the numbers of hydrogen bond contacts formed between inhibitor and targeted enzyme. Van der Waals interactions and nonpolar solvation energy also contributed in binding affinity. In contrast, for Mizolastine Van der Waals interactions were mainly responsible for complexation along with electrostatic and nonpolar interactions. However, the complexation in each system is disfavored by polar solvation energy. The calculated ∆G binding energy for Lymecycline and Mizolastine with target protein was found to be -22.19± 5.23 and -27.87±3.91 kcal/mol, respectively. Generally, Protein-ligand interaction is stabilized by different types of weak interactions with hydrogen bonding contacts being the most important interaction. Hydrogen bond analysis was performed for both the complexes, with the occupancies reported in Table 6 . In general H-bond patterns were observed with various active site residues during simulation. The inter-molecular interactions were plotted utilizing the structural ensembles from stable trajectories. In the M pro -Lymecycline complex the active residues engaged in H-bond by Lymecycline were L141, G143, S144, catalytic dyad residue C145, E166. The highest occupancy was observed for E166 (36.1%), followed by S144 (34.7%). The molecular interaction diagram (Fig. 7A ) indicated that carbonyl oxygen on cyclohexene ring from Lymecycline contributed to the formation of hydrogen bonds with S144, and G143 amide hydrogen and also with the carbonyl oxygen of S144. Moreover, E166 displayed two H-bond contacts with the side chain of Lymecycline. Other H-bonds were observed between carbonyl oxygen of methylated cyclohexane ring from Lymecycline and amide hydrogen of C145 and with the L141 residue. Mizolastine was found to make H-bond interaction with G143 (18.8%), N142 (1.17%) and C145 (0.78%) and T190 (0.33%) from the active site. Furthermore, T24, T26, T45 and S46 were also involved. To investigate the stability of interactions and to explore binding mode, interactions analysis was also carried out utilizing the same information used to calculate MMPBSA. As depicted in Fig. 7B three hydrogen bonds were observed. The Imidazole ring from benzimidazole scaffold contributes to the hydrogen bond with C145 and G143. Another H-bond was formed between the pyrimidinone ring and T26 (0.83%). It is important to note here that from the H-bond interactions analyzed in docking simulation only the interaction with G143 was maintained till the end of the simulation. For T24 the H-bond occupancy found was 11.2 and 3.45% however, the distance increased from 2.1 before MD simulation to 3.6 Å in the most stable structure during simulation. This blind molecular docking study proposes potential available approved drugs: Lymecycline and Mizolastine as prospective inhibitors of SARS-CoV-2 M pro . The estimated free energy of binding of these two drugs is in the range of -8.87 and -8.71 kcal/mol, proffering the spontaneous and energetically favored production of protein ligand complex. Molecular dynamic study concludes that both of the tested compounds, can serve as potential M pro inhibitors based on their stable nature during a long term molecular dynamic simulation. Though the maximum number of hydrogen bonds were formed by Lymecycline during the simulation the Mizolastine was stabilized in the pocket mainly via non-bonded interactions as indicated from energy component analysis. In summary, the data from RMSD, RMSF, Rg and MMPBSA indicated that there was very insignificant/minor difference between the two in imparting stability to the system. Therefore, these are potential candidates for further in vitro testing and Anti SARS-CoV-2 therapeutics development. The species Severe acute respiratory syndromerelated coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2 Will novel virus go pandemic or be contained? Going global-Travel and the 2019 novel coronavirus Global patterns in coronavirus diversity Early dynamics of transmission and control of COVID-19: a mathematical modelling study A new coronavirus associated with human respiratory disease in China Full-genome evolutionary analysis of the novel corona virus (2019-nCoV) rejects the hypothesis of emergence as a result of a recent recombination event Coronavirus main proteinase (3CLpro) structure: basis for design of anti-SARS drugs Research and development on therapeutic agents and vaccines for COVID-19 and related human coronavirus diseases Structures of the Middle East respiratory syndrome coronavirus 3C-like protease reveal insights into substrate specificity Drug repurposing for coronavirus (COVID-19): in silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes Novel 2019 Coronavirus Structure, Mechanism of Action, Antiviral drug promises and rule out against its treatment Structure of M pro from SARS-CoV-2 and discovery of its inhibitors The crystal structures of severe acute respiratory syndrome virus main protease and its complex with an inhibitor Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 Protease against COVID-19 Ul-Haq, Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach Moroccan Medicinal plants as inhibitors of COVID-19: Computational investigations In silico studies on therapeutic agents for COVID-19: Drug repurposing approach A molecular modeling approach to identify effective antiviral phytochemicals against the main protease of SARS-CoV-2 Marine natural compounds as potents inhibitors against the main protease of SARS-CoV-2, A molecular dynamic study Screening of Chloroquine, Hydroxychloroquine and its derivatives for their binding affinity to multiple SARS-CoV-2 protein drug targets Identification of new anti-nCoV drug chemical compounds from Indian spices exploiting SARS-CoV-2 main protease as target DrugBank: a comprehensive resource for in silico drug discovery and exploration The Protein Data Bank: A computer-based archival file for macromolecular structures SwissDock, a protein-small molecule docking web service based on EADock DSS Pymol: An open-source molecular graphics tool UCSF Chimera a visualization system for exploratory research and analysis Proteins Plus: a web portal for structure analysis of macromolecules Automated drawing of structural molecular formulas under constraints Molecular complexes at a glance: automated generation of two-dimensional complex diagrams Combining global and local measures for structure-based druggability predictions SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules ProTox: a web server for the in silico prediction of rodent oral toxicity The Amber biomolecular simulation programs Simmerling, ff14SB: improving the accuracy of protein side chain and backbone parameters from ff99SB Development and testing of a general amber force field Particle mesh Ewald: An N⋅ log (N) method for Ewald sums in large systems Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes VMD: visual molecular dynamics MMPBSA. py: an efficient program for endstate free energy calculations Molecular docking and dynamics simulation of FDA approved drugs with the main protease from 2019 novel coronavirus P-glycoprotein expression affects the intracellular concentration and antiviral activity of the protease inhibitor saquinavir in a T cell line May the drug transporter P glycoprotein affect the antiviral activity of human immunodeficiency virus type 1 proteinase inhibitors? The Antibiotic and Nonantibiotic Tetracyclines The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.