Mahmud Shafi, Islam Md Jahirul, Parves Md Rimon, Khan Md Arif, Tabussum Lamiya, Ahmed Sinthyia, Ali Md Ackas, Fakayode Sayo O, Halim Mohammad A
Division of Computer Aided Drug-Design, The Red-Green Research Center, BICCB, Tejgaon, Dhaka, Bangladesh.
Department of Genetic Engineering and Biotechnology, University of Rajshahi, Rajshahi, Bangladesh.
J Biomol Struct Dyn. 2022;40(19):9403-9415. doi: 10.1080/07391102.2021.1930159. Epub 2021 Jun 1.
The multidrug transporter P-glycoprotein is an ATP binding cassette (ABC) exporter responsible for resistance to tumor cells during chemotherapy. This study was designed with computational approaches aimed at identifying the best potent inhibitors of P-glycoprotein. Although many compounds have been suggested to inhibit P-glycoprotein, however, their information on bioavailability, selectivity, ADMET properties, and molecular interactions has not been revealed. Molecular docking, ADMET analysis, molecular dynamics, Principal component analysis (PCA), and binding free energy calculations were performed. Two compounds D1 and D2 showed the best docking score against P-glycoprotein and both compounds have 4-thiazolidinone derivatives containing indolin-3 one moiety are novel anti-tumor compounds. ADMET calculation analysis predicted D1 and D2 to have acceptable pharmacokinetic properties. The MD simulation discloses that D1-P-glycoprotein and D2-P-glycoprotein complexes are in stable conformation as apo-form. Hydrophobic amino acid such as phenylalanine plays significant on the interactions of inhibitors. Principal component analysis shows that both complexes are relatively similar variables as apo-form except planarity and Columbo energy profile. In addition, Quantitative Structural Activity Relationship (QSAR) of the ligand candidates were subjected to the principal component analysis (PCA) for pattern recognition. Partial-least-square (PLS) regression analysis was further utilized to model drug candidates' QSAR for subsequent prediction of the binding energy of validated drug candidates. PCA revealed groupings of the drug candidates based on the similarity or differences in drug candidates QSAR. Moreover, the developed PLS regression accurately predicted the values of the binding energy of drug candidates, with low residual error of prediction.Communicated by Ramaswamy H. Sarma.
多药转运蛋白P-糖蛋白是一种ATP结合盒(ABC)转运体,在化疗期间负责肿瘤细胞的耐药性。本研究采用计算方法设计,旨在鉴定P-糖蛋白的最佳强效抑制剂。尽管已提出许多化合物可抑制P-糖蛋白,然而,它们在生物利用度、选择性、ADMET性质和分子相互作用方面的信息尚未揭示。进行了分子对接、ADMET分析、分子动力学、主成分分析(PCA)和结合自由能计算。两种化合物D1和D2对P-糖蛋白显示出最佳对接分数,且两种化合物均为含有吲哚啉-3-酮部分的4-噻唑烷酮衍生物,是新型抗肿瘤化合物。ADMET计算分析预测D1和D2具有可接受的药代动力学性质。分子动力学模拟表明,D1-P-糖蛋白和D2-P-糖蛋白复合物呈脱辅基形式的稳定构象。疏水性氨基酸如苯丙氨酸在抑制剂的相互作用中起重要作用。主成分分析表明,除平面性和哥伦布能量分布外,两种复合物作为脱辅基形式的变量相对相似。此外,对配体候选物的定量构效关系(QSAR)进行主成分分析(PCA)以进行模式识别。进一步利用偏最小二乘(PLS)回归分析对候选药物的QSAR进行建模,以便随后预测经过验证的候选药物的结合能。PCA根据候选药物QSAR的相似性或差异揭示了候选药物的分组。此外,所开发的PLS回归准确预测了候选药物结合能的值,预测残差较低。由Ramaswamy H. Sarma通讯。