基于苯基三唑啉酮衍生物作为原卟啉原氧化酶(PPO)抑制剂的综合分子建模研究

《Chemistry & Biodiversity》：A Comprehensive Molecular Modeling Study of Phenyltriazolinone Derivatives as Protoporphyrinogen Oxidase (PPO) Inhibitors

　　

抑制原卟啉原IX氧化酶(PPO)的除草剂构成了一类高效且应用广泛的农用化学品。本研究采用综合计算策略，研究了一系列含有五元杂环的苯基三唑啉酮类除草剂的PPO抑制活性。该方法结合了定量构效关系(QSAR)建模、分子对接、分子动力学(MD)模拟和泊松-玻尔兹曼表面积(PBSA)计算。最近合成的该类化合物旨在修饰已知PPO抑制剂的化学结构并增强其除草性能。通过对现有生物活性数据的QSAR建模，研究人员提出了新的衍生物，其中化合物P7的预测pIC50值为7.11。另外两个候选物P4和P5也被预测显示出比参比化合物磺草酮(sulfentrazone)更高的除草活性。分子对接和分子动力学分析结果进一步支持了这些发现，证实了其与PPO酶的结合模式有利且相互作用稳定。此外，还通过PBSA方法计算吉布斯自由能(Gibbs free energy)评估了结合亲和力。最后，研究人员为最有前景的化合物提出了可行的合成路线，以支持其作为下一代PPO抑制剂的潜在开发。

该研究聚焦于开发新型原卟啉原氧化酶(PPO)抑制剂以应对日益严峻的杂草抗性挑战。PPO是叶绿素和血红素生物合成的关键酶，抑制该酶会导致原卟啉IX在细胞质中积累，在光照下产生活性氧，最终导致植物细胞死亡。然而，现有PPO抑制剂的长期单一使用已导致杂草产生抗性，因此开发更具选择性和可持续的化合物迫在眉睫。为此，研究人员设计了一项综合计算研究，旨在通过分子建模手段指导新型苯基三唑啉酮类衍生物的理性设计，相关成果发表在《Chemistry》期刊上。

为实现这一目标，研究人员构建了包含33种苯基三唑啉酮衍生物的数据库，并采用了多尺度计算方法。关键技术手段包括：基于图像分析定量构效关系(MIA-QSAR)建模用于预测新衍生物的pIC₅₀值；分子对接研究用于阐明配体与PPO酶的结合模式；分子动力学(MD)模拟用于评估复合物在200纳秒时间尺度上的稳定性；以及分子力学-泊松-玻尔兹曼表面积(MM-PBSA)计算用于精确量化结合自由能(ΔG_binding)。此外，还利用DeepPK服务器对设计出的新化合物进行了毒理学终点预测。

1. QSAR分析与新衍生物设计

研究人员首先利用MIA-QSAR模型分析了33种已知化合物的构效关系。通过威廉姆斯图(Williams plots)确认数据集中无异常值后，建立了基于范德华半径(r_vdW)和鲍林电负性(ε)两种描述符的模型。统计分析表明，基于r_vdW的模型预测能力略优，显示出空间位阻效应在PPO结合位点的重要性。通过变量重要性投影(VIP)和偏最小二乘(PLS)回归系数图分析，研究人员识别出R¹位的-CF₃和-F基团以及R²位的-CH₃基团对提高生物活性至关重要。基于此，研究人员设计了8种新衍生物(P1-P8)，其中P4、P5和P7显示出高于6.50的pIC₅₀预测值，特别是P7的预测值高达7.11，优于参考药物磺草酮(6.70)。

2. 分子对接研究

为了验证计算预测，研究人员对高活性化合物(20、P4、P5、P7)及参考化合物进行了分子对接。对接结果显示，磺草酮与PPO活性口袋中的精氨酸(Arg98)、甘氨酸(Gly175)等残基形成碳氢键和π-π堆积作用。化合物20和提出的P4、P5、P7均能与关键氨基酸残基形成稳定的疏水相互作用和卤键。值得注意的是，P4、P5和P7获得了比化合物20和磺草酮更负的MolDock评分，从能量角度证实了含氟取代基在增强配体-受体结合效率方面的优势。

3. 分子动力学与MM-PBSA计算

分子动力学模拟结果表明，所有蛋白-配体复合物在200纳秒模拟期间均表现出良好的稳定性。均方根偏差(RMSD)分析显示，配体P4和20的波动较小，稳定在约0.075纳米，而P7虽在初期有波动但最终也趋于稳定。均方根波动(RMSF)分析表明蛋白质骨架柔性较低。随后的MM-PBSA计算揭示了结合自由能的细节：P4表现出最负的结合自由能(-27.27 kcal mol^-1)，优于参考化合物20 (-16.53 kcal mol^-1)。能量分解显示范德华力(ΔE_vdW)是所有复合物的主要稳定驱动力，这进一步验证了QSAR分析中关于空间效应的结论。

4. 毒理学预测

利用DeepPK进行的毒性预测显示，新设计的化合物P4、P5和P7在AMES致突变性、鸟类毒性和生物降解性方面均被预测为“安全(Safe)”。特别值得关注的是，与磺草酮被预测为“有毒(Toxic)”不同，所有新提案对蜜蜂均显示出“安全”的预测结果，这表明新化合物可能在保护传粉媒介方面具有生态优势。

综上所述，该研究通过整合QSAR、对接、MD和PBSA等多种计算技术，成功设计出三种具有高预测活性和良好热力学稳定性的新型PPO抑制剂(P4、P5、P7)。研究证实，在R¹位引入强吸电子基团(如-CF₃、-F)以及在R²位引入甲基是提高活性的关键策略。这项工作不仅提供了具有应用前景的候选分子，也为未来除草剂的理性设计提供了可靠的计算范式。