TY - JOUR
T1 - Molecular docking, dynamics, and quantum chemical study of vanillylacetone and beta-hydroxy ketone derivatives against Mproof SARS-CoV-2
AU - Amin, Saniyah
AU - Muhammad, Shabbir
AU - Iqbal, Javed
AU - Ullah, Sami
AU - Al-Sehemi, Abdullah G.
AU - Algarni, H.
AU - Alarfaji, Saleh S.
AU - Ayub, Khurshid
N1 - KAUST Repository Item: Exported on 2022-12-23
Acknowledgements: The authors extend their appreciation to the Institute of Research and Consulting Studies at King Khalid University for supporting this research through grant number 2-N-20/22 and the support of Research Center for Advanced Materials Science is highly acknowledged. For computer time, this research used the resources of the Supercomputing Laboratory at King Abdullah University of Science & Technology (KAUST) in Thuwal, Saudi Arabia.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.
PY - 2022
Y1 - 2022
N2 - 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.
AB - 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.
UR - http://hdl.handle.net/10754/686606
U2 - 10.17159/0379-4350/2022/v76a12
DO - 10.17159/0379-4350/2022/v76a12
M3 - Article
SN - 1996-840X
VL - 76
SP - 79
EP - 90
JO - SOUTH AFRICAN JOURNAL OF CHEMISTRY-SUID-AFRIKAANSE TYDSKRIF VIR CHEMIE
JF - SOUTH AFRICAN JOURNAL OF CHEMISTRY-SUID-AFRIKAANSE TYDSKRIF VIR CHEMIE
ER -