key: cord-0858402-f03m8aod
authors: EL Khatabi, Khalil; Aanouz, Ilham; Alaqarbeh, Marwa; Ajana, Mohammed Aziz; Lakhlifi, Tahar; Bouachrine, Mohammed
title: Molecular docking, molecular dynamics simulation, and ADMET analysis of levamisole derivatives against the SARS-CoV-2 main protease (M(Pro))
date: 2021-10-02
journal: Bioimpacts
DOI: 10.34172/bi.2021.22143
sha: a31e149942674d5aac33e6089c216f9f86c106f4
doc_id: 858402
cord_uid: f03m8aod

[Image: see text] Introduction: The new species of coronaviruses (CoVs), SARS-CoV-2, was reported as responsible for an outbreak of respiratory disease. Scientists and researchers are endeavoring to develop new approaches for the effective treatment against of the COVID-19 disease. There are no finally targeted antiviral agents able to inhibit the SARS-CoV-2 at present. Therefore, it is of interest to investigate the potential uses of levamisole derivatives, which are reported to be antiviral agents targeting the influenza virus. Methods: In the present study, 12 selected levamisole derivatives containing imidazo[2,1-b]thiazole were subjected to molecular docking in order to explore the binding mechanisms between these derivatives and the SARS-CoV-2 M(pro) (PDB: 7BQY). The levamisole derivatives were evaluated for in silico ADMET properties for wet-lab applicability. Further, the stability of the best-docked complex was checked using molecular dynamics (MD) simulation at 20 ns. Results: Levamisole derivatives and especially molecule N°6 showed more promising docking results, presenting favorable binding interactions as well as better docking energy compared to chloroquine and mefloquine. The results of ADMET prediction and MD simulation support the potential of the molecule N°6 to be further developed as a novel inhibitor able to stop the newly emerged SARS-CoV-2. Conclusion: This research provided an effective first line in the rapid discovery of drug leads against the novel CoV (SARS-CoV-2).

In late December 2019, the entire world has witnessed an unexpected outbreak of novel coronavirus disease-19 . The virus started in Wuhan city, China has rapidly turned out to be a global pandemic and got spread worldwide very fast. 1 The virus is recognized to be a considerable threat to the global health systems. As of the 13 th of July 2020, the novel coronavirus has infected about 4.5 million patients all over the world with a mortality rate reaching approximatively 6.7%. 2 SARS-CoV-2 is one of the three single strand-RNA beta-coronaviruses in addition to the Middle East respiratory syndrome (MERS) virus and severe acute respiratory syndrome (SARS) virus identified as highly developing viruses. 3 Coronaviruses (CoVs) are an etiological factor of severe infections belonging to a large family of viruses causing disorders in the respiratory tract infections in both humans and other mammals. 4 In addition to the respiratory tract, the virus can also cause disorders in the digestive tract. The virus is transmitted primarily from an infected human to another through saliva, viral droplets after sneezing and coughing, or touching an infected surface. Severe pneumonia, fever, sore throat, sneezing, cough, loss of taste and smell are the symptoms of coronavirus infection. [5] [6] [7] [8] As the most study. The chemical structure of the levamisole is shown in Fig. 1 . The selected levamisoles containing imidazo[2, 1-b] thiazole moiety were prepared by Gürsoy et al 15 and evaluated against diverse DNA and RNA viruses, using different cell-based assays cells infected with influenza B virus or influenza A/H3N2 and A/H1N1. The general chemical structures reported as antiviral agents are depicted in Fig. 2 , while chemical groups representing n, R 1, and R 2 for 12 Levamisole derivatives are available in Table 1 .

The 3D structures of ligands were built using SYBYL-X 2.0 (Tripos Inc., St. Louis, USA). They were energetically minimized using Tripos force field, 19 Gasteiger-Huckel charges, and Powell method. 20 Receptor The crystal structure of SARS-CoV-2 M Pro (PDB: 7BQY, resolution: 1.7Å), 20 which reacts in such a way to affect the host cells like the parabronchial epithelial and pulmonary cell, was obtained from the protein databank.

The 7BQY protein is classified as a viral protein and involved two chains, A and B, which subordinate to form a homodimer by a crystallographic 2-fold symmetry axis. Chain A was employed for macromolecule preparation of a sequence of 306. The original ligand for 7BQY is (E,4S)-4-azanyl-5-[(3S)-2-oxidanylidenepyrrolidin-3-yl]pent-2enoic acid. Molecular docking studies Molecular docking is powerful in silico approach, effective way to strengthen the immune defense against the SARS-CoV-2, prophylactic vaccines should have essential roles. 9 Live attenuated or killed whole virus vaccines, recombinant protein vaccine, Ii-Key peptide COVID-19 as subunit protein vaccine, multi-epitope vaccine, DNA-and mRNA-based vaccines, convalescent sera or plasma therapy, and S-trimer are examples of vaccine candidates to be considered for more rapid development of prophylactics and reinforce immunity system against SARS-CoV-2 virus with minimum side-effects. [10] [11] [12] [13] At this time, no specific drugs or treatments are available for this virus reported. Repurposing medication seems to be the only way to develop potent drugs to control the pandemic and stop viral infection. Levamisole is a drug having important immunomodulatory properties and was reported in the literature to be one of the most potent drug families in the fight against influenza virus antiviral agents is considered to be an attractive target for an antiviral drug. 14 Gürsoy et al evaluated the levamisole derivatives against DNA and RNA viruses, using diverse cell-based assays infected with influenza B virus or influenza, A/ H3N2 and A/H1N1. 15 Therefore, this work evaluates the levamisole derivatives as a potential inhibitor of the SARS-CoV-2 M Pro through computational studies such as molecular docking and ADMET prediction. The molecular docking results of the studies compounds were compared to the docked molecules chloroquine and mefloquine which have been originally developed for viral inhibition. Chloroquine is at present experiencing for the treatment of SARS-CoV-2. 16 It is an amino quinolone derivative indicated to treat malaria as well as rheumatoid arthritis, HIV, and prophylaxis of Zika virus, 17 while mefloquine is an antimalarial drug acting as a blood schizonticides that has the potential to be applied to inhibit the erythrocytic forms of plasmodium species and some sorts of influenza viruses. 18 The docking of the levamisole, chloroquine, and mefloquine was performed using the same protease (PDB: 7BQY) to provide molecular insights into the stability of the complexes.

Consequently, one levamisole derivative with higher affinity was prioritized and proceeded to a detailed analysis through molecular docking, then compared to the docked molecules Chloroquine and Mefloquine. MD simulations were applied through RMSD, RMSF, and Rg properties to identify the stabilization of the protein-ligand complexes. The prioritized molecule showed more promising results as a potential inhibitor which provided a suitable basis for in vivo and in vitro anti-viral studies for further drug discovery research against SARS-CoV-2 disease.

Ligands A series of 12 selected Levamisole derivatives previously assessed for their antiviral activities were collected from literature, 15 and were subjected to molecular docking Computational approaches to identify potent antiviral agents against SARS-CoV-2 protease BioImpacts, 2022, 12(2), 107-113 109 performed using Surflex-Dock module in SYBYL-X v2.0, 21 to investigate the binding interaction of compounds and forecast the optimized conformation of a ligand and protein target. The preparation of the protein structure was done by removing water molecules and other atoms and adding polar hydrogen atoms. The 12 molecules were docked in the active site of the protein, and the ligandreceptor interactions were studied. 22, 23 The Total-score of Surflex-Dock was used to represent binding affinities, which was selected as the docking result. The total Surflex-Dock score was expressed in -log10 (Kd) units.

The pharmacokinetics ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties of the 12 selected levamisole derivatives were predicted using pKCSM and admetSAR predictor. 24, 25 MD simulations To further confirm the reliability of molecular docking and reveal the binding mode and conformational changes during the interaction between the ligand (compound 6) and receptor protein, molecular dynamics simulations of the complex was carried out using GROMACS 5.1.4 software. 26, 27 Protonation and minimization steps were applied for the system using GROMOS96 43A1 force field to avoid space crashes. The molecular topology files parameters for the ligand were generated using PRODRG server. 28 The docked complex was solvated and immersed in water cubic boxes with 12 Å as a margin distance. Four ions of Na + were added to neutralize the systems. The energy minimization was carried on the system through 5000 steps of the steepest descent minimization to eliminate the weak van der Waals contacts. The system was further equilibrated to carry out 20 ns MD simulations while the pressure and temperature were set at 1atm and 300 K, respectively. All covalent bonds including heavy atom-H were constrained by the LINCS algorithm. MD simulation was run for 20 000 ps for the docked complex, writing coordinates every two ps interval. Finally, the resulting trajectories of simulated systems were saved for detailed analysis. The root-mean-square deviation (RMSD), rootmean-square fluctuation (RMSF), and radius of gyration (Rg) were analyzed throughout the trajectory using the gmxrms, RMSF, and gyrate, respectively, the built-in function of the GROMACS software. The graphs plotting was made by the QtGrace.22 program.

A molecular docking study was performed to explore the binding affinity, the binding types, and the active amino acid residues of studied molecules in the target enzyme. The twenty-two selected molecules, chloroquine and Mefloquine were docked into the binding site of SARS-CoV-2 main protease and evaluated for their affinity. The binding affinity values of the 12 selected compounds are reported in Table 2 .

The best (or lowest) energy level of molecules originating from the contribution of various interactions with the SARS-CoV-2 M Pro (PDB: 7BQY) is observed for molecule N°6 (Table 2) . This molecule could have more inhibitory potential of the studied enzyme than the complexes formed with the reference molecules Chloroquine and Mefloquine. The results of the re-docked Chloroquine and Mefloquine molecules and their position in the SARS-CoV-2 M Pro are shown in Figs. 3 and 4 , respectively. Consequently, the chosen molecule was prioritized and proceeded to a detailed analysis as shown in Fig. 5 .

It is necessary to carry out the applicability of the 12 selected compounds especially molecule N°6 as antiviral agents using virtual properties before the experiment. The ADMET properties were predicted, the results are presented in Table 3 .

To assure the structural stability of the protein-ligand complex in the binding site of SARS-CoV-2 M Pro and analyze the conformational changes of the system, the selected molecule N°6 was prioritized and subjected to 20 000 ps MD simulations. As result, several important dynamics properties of the simulated complex were generated.

The root-mean-square deviation (RMSD), the rootmean-square fluctuation (RMSF), and the radius of gyration (Rg) were evaluated to gain an insight into the overall stability of the complex during the simulation as shown in Figs. 6A, 6B, and 6C, respectively. 

Based on the further analysis of the docking results, hydrogen bonds and some non-classical hydrogen bonds could be also observed for the Chloroquine binding process. The hydrogen bond was generated with residue Cys145 at a distance of 2.67 Å, while the aliphatic chain was in carbon-hydrogen bond contact with His163, Glu166, Phe140, and Asn142. Mefloquine is involved in hydrogen bond interactions with Cys145 and Asn142. The piperidine ring of Mefloquine forms non-classical hydrogen bonds with His163 and Leu141, while halogen interaction is created only with Gln189 because of the fluorine atom (F). The unfavorable donor-donor interaction is observed with Gly143 amino acid. Fig. 5 shows the interaction pose of the active molecule (with low energy), which is mentioned as molecule N°6 in the present study. Docking study of molecule N°6 in SARS-CoV-2 M Pro shows more number and type of favorable interactions in comparison with Chloroquine and Mefloquine of the studied enzyme (H-Bond, Van der Waals, p-Donor Hydrogen Bon, alkyl, and p-alkyl interaction). The presence of hydrogen bonding and hydrophobic interactions could impact the structure and the function of the protease and play a key role in stabilizing the complex conformation.

The molecule N°6 was able to be docked deeply within the binding site of the SARS-CoV-2 M Pro , displaying as results a favorable binding interaction as well as better docking energy compared to Chloroquine and Mefloquine. The results of the molecular docking study support the potential of the molecule N°6 to be further developed as a novel inhibitor of Coronavirus (SARS-CoV-2).

The tested molecules showed a high percentage of human intestinal absorption ranged from 83.826% to 92.273%, which were labeled to be fully absorbed by the human enteric. Molecules with logVDss< -0.15 are considered to be poorly distributed in tissues rather than in plasma. The logVDss results in Table 3 show most of the tested molecules were distributed in the plasma. For metabolism, all molecules were predicted as substrates for the CYP450 3A4 subtype, which indicated may be well metabolized by CYP 3A4. Moreover, all tested molecules could inhibit the 2C9subtype, whereas they could not inhibit the 2D6 subtype. Nevertheless, some of the cytochrome P450 subtypes might be inhibited by one or more of the tested molecules. Based on the high values of total clearance, all tested compounds could be filtered by renal and hepatic tissues in a combinational way. The predicted toxicity indicated all tested molecules were non-mutagenic and non-carcinogen, indicating the safety of the molecules which is important to develop an active drug.

The tested molecules demonstrated promising results of computational pharmacokinetics and toxicity evaluations, which could serve as promising inhibitors targeting SARS-CoV-2 for further development. The molecule N°6which exhibited the good potential to be a CoV inhibitor compared to Chloroquine and Mefloquine showed a reasonably good ADMET profile. Fig. 6A displays the RMSD for heavy atoms of the protein during 20 000 ps, while Fig. 6B shows the RMSF during 20 000 ps. In the RMSD graph, the studied ligand showed fluctuation during the initial simulations till about 10 000 ps. A negligible deviation was detected between 10 000 ps and 15 000 ps, while stable dynamics were observed thereafter throughout the simulation. It can be clearly seen that the ligand reached equilibrium at a lower RMSD value of 0.27 nm, indicating the good stability of the system.

As depicted in Fig. 6B , the RMSF was also estimated to assess the flexibility of the protein residues. The results exhibited rigid flexibility of residue fluctuation profile of the protein with an average RMSF of 0.2 nm. This indicated that the binding of ligand to protein showed well stable interaction and did not result in a large movement in protein conformation. Rg plot was measured (Fig. 6C) . The Rg imparts the compactness of a structure throughout the simulation. The Rg-time fluctuations were observed almost constants within the acceptable range, mostly maintained between 2.1 and 2.18 nm, indicating steady conformation changes of the protein. The obtained results revealed that the complex could keep stable throughout the simulation and interact well with the SARS-CoV-2 main protease (M Pro ).

To summarize the above results, levamisole derivatives and especially molecule N°6 showed more promising results as a potential inhibitor able to stop the newly emerged SARS-CoV-2 disease. Overall, this work might provide valuable insights that have the potential to be applied in designing novel SARS-CoV-2 inhibitors for further drug discovery research.

Molecular docking was carried out for 12 levamisole derivatives with the SARS-CoV-2 M Pro (PDB: 7BQY). The studied compounds were evaluated for in silico ADMET prediction and the results showed a reasonably good ADMET profile. Molecular docking results indicate that the hydrogen bonds, van der Waals and hydrophobic What is the current knowledge? √ We expect that our approach; the structure-based molecular docking and MD simulations, may provide new promising avenues for optimizing and developing more potent COVID-19 inhibitors. √ Levamisole derivatives with strong binding were prioritized and proceeded to a detailed analysis through molecular docking and MD simulations.

What is new here? √ Levamisole derivatives, especially molecule N°6 showed more promising docking results compared to Chloroquine and Mefloquine. √ Good ADMET prediction results for levamisole derivatives were acquired.

interactions were proved to be of great importance in ligand-receptor binding. The docked position of the molecule N°6 showed strong interactions, which was well stable during the 20 ns MD simulation, indicating that molecule N°6 represents a good potential to be coronavirus inhibitor with high binding affinity, which could be a promising candidate for further drug discovery research against SARS-CoV-2 disease.

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

SARS-CoV-2 Vaccines: Status Report

COVID-19: towards controlling of a pandemic

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster

Clinical Characteristics of Coronavirus Disease 2019 in China

Data sharing and outbreaks: best practice exemplified

Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China

Current status of multiple drug molecules, and vaccines: an update in SARS-CoV-2 therapeutics

COVID-19 -Recent advancements in identifying novel vaccine candidates and current status of upcoming SARS-CoV-2 vaccines

A candidate multi-epitope vaccine against SARS-CoV-2

A domain-based vaccine construct against SARS-CoV-2, the causative agent of COVID-19 pandemic: development of selfamplifying mRNA and peptide vaccines

Next-generation vaccines and the impacts of state-of-the-art in-silico technologies

Levamisole and inosiplex: antiviral agents with immunopotentiating action

Design and synthesis of novel Imidazo[2,1-b]thiazole derivatives as potent antiviral and antimycobacterial agents

Breakthrough: Chloroquine phosphate has shown apparent efficacy in treatment of COVID-19 associated pneumonia in clinical studies

Repurposing of the anti-malaria drug chloroquine for Zika Virus treatment and prophylaxis

The position of mefloquine as a 21st century malaria chemoprophylaxis

Validation of the general purpose tripos 5.2 force field

A brief review and table of semiempirical parameters used in the Hueckel molecular orbital method

Surflex-Dock 2.1: robust performance from ligand energetic modeling, ring flexibility, and knowledge-based search

Discovery Studio Modeling Environment

Pymol: An Open-Source Molecular Graphics Tool

ADMET evaluation in drug discovery. 11. PharmacoKinetics Knowledge Base (PKKB): a comprehensive database of pharmacokinetic and toxic properties for drugs

pkCSM: Predicting Small-Molecule Pharmacokinetic and Toxicity Properties Using Graph-Based Signatures

Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS

GROMACS: fast, flexible, and free

PRODRG: a tool for highthroughput crystallography of protein-ligand complexes

We are grateful to the "Association Marocaine des Chimistes Théoriciens" (AMCT) for its pertinent help concerning the programs.

There is none to be clarified.

None to be declared.

The authors declare that they have no conflict of interest.Authors' contribution KE: draft preparation, data handling, data analysis, writing, and reviewing. IA: conceptualization, data analysis, reviewing. MA: data handling and reviewing. MAA: data analysis, study justification, and supervision. TL: supervision and project administration. MB: writing and reviewing study justification, and supervision.