Abstract
This study is carried out to find novel active drug candidates which can effectively bind to key residues of main protease (Mpro) of SARS-CoV-2. We performed molecular docking of fifty-seven (57) ligands from two classes: vanillylacetone and its derivatives and beta-hydroxy ketone derivatives against Mpro of SARS-CoV-2. We also docked three antiviral drugs as reference/benchmark drugs including remdesivir (RDV), chloroquine (CQ), and hydroxychloroquine (HCQ) against Mpro for comparison of inhibition tendencies of selected ligands. Binding energies of our reference drugs are as: CQ = -5.1 kcal mol-1 (with predicted inhibition constant (Ki pred) = 177 mu mol), HCQ = -5.7 kcal mol-1 (Ki pred = 64.07 mu mol) and RDV -6.3 kcal mol-1 (Ki pred = 13.95 mu mol). We got remarkable results for our docked ligands as 79% of total ligands indicated binding energies better than CQ, 39 % better than both HCQ and CQ, and 19 % better than all reference drugs. More interestingly interaction analysis of eight best-docked ligands showed that they interacted with desired key residues of Mpro. We further selected the four best-docked ligands L1 = -6.6 kcal mol-1 (Ki pred=13.95 mu mol), L6 = -7.0 kcal mol-1 (Ki pred = 7.08 mu mol), L34 = -6.0 kcal mol-1 (Ki pred = 38.54 mu mol), and L50 = -6.6 kcal mol-1 (Ki pred=13.95 mu mol) for further analysis by quantum chemical study, molecular dynamic (MD) simulations and ADMET analysis. We have also carried out MD-simulations of six more docked ligand L2, L14, L20, L36, L46 and L48 some of which were showing weak binding affinities and some average binding affinities to check their simulation behavior. Their RMSD, RMSF and binding free energy results were also quite satisfying. We believe the current investigation will evoke the scientific community and highlights the potential of selected compounds for potential use as antiviral compounds against Mpro of SARS-CoV-2.