Computational repurposing approach for targeting the critical spike mutations in B.1.617.2 (delta), AY.1 (delta plus) and C.37 (lambda) SARS-CoV-2 variants using exhaustive structure-based virtual screening, molecular dynamic simulations and MM-PBSA methods

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the contagious coronavirus disease 2019 (COVID-19) which was first identified in Wuhan, China, in December 2019. Around the world, many researchers focused their research on identifying inhibitors against the druggable SARS-CoV-2 targets. The reported genomic mutations have a direct effect on the receptor-binding domain (RBD), which interacts with host angiotensin-converting enzyme 2 (ACE-2) for viral cell entry. These mutations, some of which are variants of concern (VOC), lead to increased morbidity and mortality rates. The newest variants including B.1.617.2 (Delta), AY.1 (Delta plus), and C.37 (Lambda) were considered in this study. Thus, an exhaustive structure-based virtual screening of a ligand library (in which FDA approved drugs are also present) using the drug-likeness screening, molecular docking, ADMET profiling was performed followed by molecular dynamics (MD) simulation, and Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculation to identify compounds or drugs can be repurposed for inhibiting the wild type, Delta, Delta plus and Lambda variants of RBD of the spike protein. Based on the virtual screening steps, two FDA approved drugs, Atovaquone (atv) and Praziquantel (prz), were selected and repurposed as the best candidates of SARS-CoV-2 RBD inhibitors. Molecular docking results display that both atv and prz contribute in different interaction with binding site residues (Gln493, Asn501 and Gly502 in the hydrogen bond formation, Phe490 and Tyr505 in the π- π stacking and Tyr449, Ser494, and Phe497 in the vdW interactions) in the wild type, Delta, Delta plus and Lambda variants of RBD of the spike protein. MD simulations revealed that among the eight studied complexes, the wild type-atv and Delta-prz complexes have the most structural stability over the simulation time. Furthermore, MM-PBSA calculation showed that in the atv containing complexes, highest binding affinity is related to the wild type-atv complex and in the prz containing complexes, it is related to the Delta-prz complex. The validation of docking results was done by comparing with experimental data (heparin in complex with wild type and Delta variants). Also, comparison of the obtained results with the result of simulation of the k22 with the studied proteins showed that atv and prz are suitable inhibitors for these proteins, especially wild type t and Delta variant, respectively. Thus, we found that atv and prz are the best candidate for inhibition of wild type and Delta variant of the spike protein. Also, atv can be an appropriate inhibitor for the Lambda variant. Obtained in silico results may help the development of new anti-COVID-19 drugs.