基于QSAR引导设计、DFT、分子对接及ADMET预测的新型噻唑烷酮-吲哚衍生物作为强效α-葡萄糖苷酶抑制剂的研究

《Letters in Drug Design & Discovery》：QSAR-guided design, DFT, molecular docking, and ADMET prediction of novel thiazolidinone-based indole derivatives as potent α-glucosidase inhibitors

　　

糖尿病是一种复杂的代谢性疾病，主要特征为因胰岛素分泌不足或作用障碍导致的高血糖水平。该疾病分为两种主要类型：1型糖尿病，即机体无法产生足够的胰岛素；以及2型糖尿病，即胰岛素利用效率低下。控制糖尿病的关键在于维持最佳血糖水平，这对于预防与该疾病相关的并发症至关重要，因为该疾病具有高死亡率和沉重的医疗负担。 对于2型糖尿病而言，一个有前景的治疗策略是抑制α-葡萄糖苷酶（α-glucosidase），这是一种关键的肠道酶，负责水解α-D-糖苷键并将葡萄糖释放入血液。抑制该酶可延缓碳水化合物消化并降低餐后血糖水平，对于以碳水化合物为主食的人群是一种有效策略。目前获批的α-葡萄糖苷酶抑制剂（AGIs）包括阿卡波糖（acarbose）、米格列醇（miglitol）和伏格列波糖（voglibose），以缓解葡萄糖吸收而闻名。尽管具有潜力，但AGIs的使用受到显著缺陷的制约，如口服生物利用度差、胃肠道副作用和疗效有限。这些挑战推动了开发具有更优药代动力学性质的高效小分子α-葡萄糖苷酶抑制剂的研究。 噻唑烷酮（Thiazolidinone）是一种含硫和氮的杂环化合物，以其结构多样性和广泛的生物活性谱而著称。噻唑烷酮的衍生物和杂合物表现出强效的抗菌、抗真菌、抗病毒、抗炎、抗癌和抗结核作用。特别相关的是，带有4-噻唑烷酮部分的噻唑烷酮衍生物已显示出作为葡萄糖苷酶抑制剂的潜力，并具有有记录的抗糖尿病活性。 吲哚（Indole）是一种由苯环和吡咯环稠合而成的双环杂芳烃骨架，因其多样的生物和治疗活性而被认为是药物化学中最重要的优势结构之一。天然和合成的吲哚衍生物表现出显著的抗癌、抗菌、抗抑郁和抗焦虑活性。基于吲哚的分子如褪黑素（melatonin）具有抗氧化和昼夜节律调节功能，并调控神经退行性疾病，包括阿尔茨海默病和帕金森病。 这两种生物活性骨架的融合由Khan等人（2022）探索，他们合成了一系列20个基于噻唑烷酮的吲哚衍生物并评估了其α-葡萄糖苷酶抑制活性，其中吲哚环5位带有氟取代基的衍生物（化合物5）表现出最强的α-葡萄糖苷酶抑制活性，IC50为2.40 μM，优于阿卡波糖（IC50 = 11.70 μM）。尽管该实验数据集和初步构效关系具有价值，但该研究缺乏详细的定量构效关系（QSAR）分析，这可以进一步阐明活性的结构决定因素并指导新类似物的合理设计。QSAR方法通过能够在合成前对候选化合物进行计算机模拟优先排序来简化药物开发，从而降低与实验筛选相关的成本和时间。这些预测模型通过分子描述符将化学结构与生物活性相关联，有助于设计具有预测抗糖尿病特性的新化合物。当与分子对接和ADMET预测相结合时，QSAR成为合理设计具有改进效力和有利药代动力学的新化合物的强大工具。 本研究通过以下方式弥补这一空白：（1）通过GA-MLR方法在20个基于噻唑烷酮的吲哚衍生物数据集上开发和验证QSAR模型，以评估分子描述符对α-葡萄糖苷酶抑制活性的影响；（2）基于最有效的母体化合物（化合物5）设计具有预测α-葡萄糖苷酶抑制潜力的新类似物；（3）利用分子对接分析有前景的类似物与α-葡萄糖苷酶活性位点的结合相互作用；（4）使用B3LYP/6–31G**水平的DFT计算表征电子反应性特征；以及（5）评估所设计类似物的类药性和ADMET性质。这一综合计算策略为识别和优先排序基于噻唑烷酮的吲哚衍生物作为潜在α-葡萄糖苷酶抑制剂的进一步研究提供了合理框架。

论文解读文章

一、研究背景与问题提出

糖尿病是一种以血糖水平升高为主要特征的复杂代谢性疾病，其根本原因在于胰岛素分泌不足或胰岛素作用障碍。该疾病分为1型糖尿病（机体无法产生足够胰岛素）和2型糖尿病（胰岛素利用效率低下）两种主要类型。控制血糖水平是预防糖尿病相关并发症的关键，该疾病具有高死亡率和沉重的医疗负担。

在2型糖尿病的治疗策略中，抑制α-葡萄糖苷酶（α-glucosidase）是一种有前景的方法。α-葡萄糖苷酶是肠道中负责水解α-D-糖苷键并将葡萄糖释放入血液的关键酶。抑制该酶可延缓碳水化合物消化并降低餐后血糖水平。目前已获批的α-葡萄糖苷酶抑制剂（AGIs）包括阿卡波糖（acarbose）、米格列醇（miglitol）和伏格列波糖（voglibose），然而这些药物存在口服生物利用度差、胃肠道副作用显著和疗效有限等缺陷，限制了其临床应用。因此，开发具有更优药代动力学性质的新型小分子α-葡萄糖苷酶抑制剂成为研究热点。

噻唑烷酮（thiazolidinone）是一种含硫和氮的杂环化合物，具有结构多样性和广泛的生物活性，其衍生物表现出抗菌、抗真菌、抗病毒、抗炎、抗癌和抗结核等多种药理作用。带有4-噻唑烷酮部分的衍生物已显示出作为葡萄糖苷酶抑制剂的潜力。吲哚（indole）则是由苯环和吡咯环稠合而成的双环杂芳烃骨架，是药物化学中最重要的优势结构之一，天然和合成的吲哚衍生物具有显著的抗癌、抗菌、抗抑郁和抗焦虑活性。

Khan等人（2022）合成了20个基于噻唑烷酮的吲哚衍生物并评估了其α-葡萄糖苷酶抑制活性，其中吲哚环5位带有氟取代基的化合物5表现出最强的抑制活性（IC₅₀ = 2.40 μM），优于阿卡波糖（IC₅₀ = 11.70 μM）。然而，该研究缺乏详细的定量构效关系（QSAR）分析，未能充分阐明活性的结构决定因素。QSAR方法能够在合成前对候选化合物进行计算机模拟优先排序，从而降低成本和缩短时间。当与分子对接和ADMET预测结合时，QSAR成为合理设计新化合物的强有力工具。本研究旨在弥补这一空白，通过综合计算策略为基于噻唑烷酮的吲哚衍生物作为潜在α-葡萄糖苷酶抑制剂的开发提供合理框架。

二、研究内容与主要结论

研究人员基于20个噻唑烷酮-吲哚衍生物数据集，通过遗传算法-多元线性回归（GA-MLR）方法建立了QSAR模型，并基于该模型设计新类似物，结合分子对接、DFT计算和ADMET预测对候选化合物进行综合评估。研究得出以下主要结论：

（1）建立的GA-MLR QSAR模型包含三个分子描述符：氢键受体数目（nHBAcc）、价链序7（VCH-7）和沿X/Y轴的转动惯量（MOMI-XY），模型方程表达为：pIC₅₀ = 9.4039 - 0.7484nHBAcc - 11.659VCH-7 + 1.3869*MOMI-XY。模型验证显示R²_train = 0.911，Q²_cv = 0.832，R²_test = 0.697，各项验证指标均达到或超过标准阈值，Y-随机化分析确认模型不存在偶然相关性（cR²_p = 0.821），Williams图显示所有化合物均位于适用域内。

（2）基于化合物5作为结构模板，研究人员设计了三个新类似物D1、D2和D3。其中D3表现出最优的预测活性（pIC₅₀ = 5.566，IC₅₀ = 2.72 μM），几乎与母体化合物5相当（IC₅₀ = 2.40 μM），且明显优于阿卡波糖（IC₅₀ = 11.70 μM）。

（3）分子对接研究（PDB ID: 3W37）显示，D3的结合亲和力为-6.2 kcal/mol，与化合物5（-6.4 kcal/mol）相当。D3通过卤键（含-CF₃基团与ILE233和ASP232）、Pi-硫相互作用（与TRP329）、Pi-阴离子相互作用（与ASP469）、三个Pi-Pi T型相互作用（与PHE601和PHE236）以及Pi-烷基相互作用（与ALA234）维持结合，采用了一种不依赖常规氢键的替代结合策略。

（4）DFT计算（B3LYP/6–31G^**水平）表明，D3具有较小的能隙（ΔEg = 3.07 eV）和较高的亲电性指数（ω = 4.23 eV），显示出优于化合物5和阿卡波糖的化学反应性和动力学稳定性。

（5）ADMET预测显示，D3无Lipinski规则违例，胃肠吸收率为92.95%，表现出良好的类药性和药代动力学特征。

三、主要技术方法

本研究采用的主要技术方法包括：（1）GA-MLR QSAR建模：使用遗传算法结合多元线性回归进行变量选择和模型构建，数据集通过Kennard-Stone算法按75%:25%划分训练集和测试集；（2）分子对接：使用MOE 2015.10软件，采用Triangle Matcher放置算法、London dG评分函数和GBVI/WSA dG精化，对α-葡萄糖苷酶（PDB ID: 3W37）进行对接模拟；（3）DFT计算：使用Spartan 2014软件在B3LYP/6–31G^**水平进行几何优化和电子结构计算；（4）ADMET预测：使用pkCSM在线服务器评估类药性和毒性。样本数据来源于Khan等人（2022）报道的20个化合物及其实验活性数据。

四、研究结果

QSAR模型开发与验证： 通过GA-MLR方法建立了包含三个描述符的QSAR模型，模型具有良好的拟合优度（R²_train = 0.911）和预测能力（Q²_cv = 0.832，R²_test = 0.697）。描述符分析表明，氢键受体数目增加会降低抑制活性，分子外围分支（VCH-7）增加同样降低活性，而分子伸长或各向异性质量分布（MOMI-XY）则增强活性。

新类似物设计： 基于QSAR模型和SAR分析，以化合物5为模板设计了D1、D2、D3三个新类似物。D3（含-CF₃取代基）预测IC₅₀ = 2.72 μM，活性最接近母体化合物，被选为进一步研究对象。

分子对接验证与结合模式分析： 通过将共结晶配体阿卡波糖重新对接至α-葡萄糖苷酶活性位点进行方法验证，RMSD = 1.91 ?（低于2.0 ?阈值），确认对接方案可靠。D3通过卤键、Pi-硫、Pi-阴离子、Pi-Pi T型和Pi-烷基等多种非氢键相互作用维持结合，采用与化合物5不同的替代结合策略。

DFT电子结构表征： D3的HOMO-LUMO能隙（ΔEg = 3.07 eV）小于化合物5和阿卡波糖，亲电性指数（ω = 4.23 eV）更高，表明其具有更强的化学反应活性和电子转移能力。

ADMET性质评估： D3无Lipinski规则违例，分子量、Log P、氢键供体和受体数目均符合口服药物标准，胃肠吸收率高达92.95%，具有良好的类药性和口服生物利用度潜力。

五、讨论与结论

讨论部分指出，D3虽然未能形成常规氢键，但通过卤键和疏水相互作用网络实现了与化合物5相当的结合亲和力，这种非氢键介导的抑制策略为未来优化提供了新方向。D3的预测活性（IC₅₀ = 2.72 μM）优于阿卡波糖（IC₅₀ = 11.70 μM），且具有良好的ADMET特征，弥补了阿卡波糖口服生物利用度差和胃肠道副作用等缺陷。

研究结论表明，所建立的QSAR模型可靠且具有预测能力，成功指导了新型类似物D3的设计。D3在预测活性、结合亲和力、化学反应性和药代动力学性质方面均表现出色，与母体化合物5活性相当且优于标准药物阿卡波糖，可作为有前景的α-葡萄糖苷酶抑制剂候选化合物，值得进一步的实验合成和体外活性验证。该研究为基于噻唑烷酮-吲哚骨架的抗糖尿病药物开发提供了系统的计算设计框架。