Rafique Amina, Muhammad Shabbir, Iqbal Javed, Al-Sehemi Abdullah G, Alshahrani Mohammad Y, Ayub Khurshid, Gilani Mazhar Amjad
Department of Chemistry, University of Agriculture, Faisalabad 38000, Pakistan.
Department of Chemistry, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia.
J Mol Liq. 2023 Jul 15;382:121904. doi: 10.1016/j.molliq.2023.121904. Epub 2023 Apr 26.
In the current study, a hybrid computational approach consisting of different computational methods to explore the molecular electronic structures, bioactivity and therapeutic potential of piperidine compounds against SARS-CoV-2. The quantum chemical methods are used to study electronic structures of designed derivatives, molecular docking methods are used to see the most potential docking interactions for main protease (M) of SARS-CoV-2 while molecular dynamic and MMPBSA analyses are performed in bulk water solvation process to mimic real protein like aqueous environment and effectiveness of docked complexes. We designed and optimized piperidine derivatives from experimentally known precursor using quantum chemical methods. The UV-Visible, IR, molecular orbitals, molecular electrostatic plots, and global chemical reactivity descriptors are carried out which illustrate that the designed compounds are kinetically stable and reactive. The results of MD simulations and binding free energy revealed that all the complex systems possess adequate dynamic stability, and flexibility based on their RMSD, RMSF, radius of gyration, and hydrogen bond analysis. The computed net binding free energy ) as calculated by MMPBSA method for the complexes showed the values of -4.29 kcal.mol for P1, -5.52 kcal.mol for P2, -6.12 kcal.mol for P3, -6.35 kcal.mol for P4, -5.19 kcal.mol for P5, 3.09 kcal.mol for P6, -6.78 kcal.mol for P7, and -6.29 kcal.mol for P8.The ADMET analysis further confirmed that none of among the designed ligands violates the Lipinski rule of five (RO5). The current comprehensive investigation predicts that all our designed compounds are recommended as prospective therapeutic drugs against M of SARS-CoV-2 and it provokes the scientific community to further perform their in-vitro analysis.
在当前研究中,采用了一种混合计算方法,该方法由不同的计算方法组成,用于探索哌啶化合物对严重急性呼吸综合征冠状病毒2(SARS-CoV-2)的分子电子结构、生物活性和治疗潜力。量子化学方法用于研究设计衍生物的电子结构,分子对接方法用于观察SARS-CoV-2主要蛋白酶(M)最具潜力的对接相互作用,而分子动力学和MMPBSA分析则在大量水溶剂化过程中进行,以模拟真实的蛋白质样水环境和对接复合物的有效性。我们使用量子化学方法从实验已知的前体设计并优化了哌啶衍生物。进行了紫外可见光谱、红外光谱、分子轨道、分子静电图和全局化学反应性描述符分析,结果表明所设计的化合物在动力学上是稳定且具有反应活性的。分子动力学模拟和结合自由能的结果表明,基于均方根偏差(RMSD)、均方根波动(RMSF)、回转半径和氢键分析,所有复合系统都具有足够的动态稳定性和灵活性。通过MMPBSA方法计算的复合物的净结合自由能显示,P1为-4.29 kcal·mol,P2为-5.52 kcal·mol,P3为-6.12 kcal·mol,P4为-6.35 kcal·mol,P5为-5.19 kcal·mol,P6为3.09 kcal·mol,P7为-6.78 kcal·mol,P8为-6.29 kcal·mol。ADMET分析进一步证实,所设计的配体均未违反Lipinski五规则(RO5)。当前的综合研究预测,我们所有设计的化合物都被推荐作为针对SARS-CoV-2 M蛋白的潜在治疗药物,这促使科学界进一步进行体外分析。