《Scientific African》：Investigating the Herbicidal Potential of Substituted Urea Derivatives as HPPD Inhibitors: A Comparative Computational Study with Standard Mesotrione
编辑推荐：
摘要
抗除草剂杂草的日益普遍存在要求开发具有更高效率、选择性和环境友好性的新型除草剂。在本研究中，研究人员采用分子对接、MM-GBSA结合自由能计算、理化性质与类除草剂评估、ADMET与毒性预测、分子动力学（MD）模拟及密度泛函理论（DFT）分析组成的联合工作
摘要
抗除草剂杂草的日益普遍存在要求开发具有更高效率、选择性和环境友好性的新型除草剂。在本研究中，研究人员采用分子对接、MM-GBSA结合自由能计算、理化性质与类除草剂评估、ADMET与毒性预测、分子动力学（MD）模拟及密度泛函理论（DFT）分析组成的联合工作流程，对一系列取代脲衍生物作为4-羟苯基丙酮酸双加氧酶（HPPD）潜在抑制剂进行了计算评估。分子对接结果确认化合物2为最有前景的候选分子，其展现出优于参考除草剂硝磺草酮（Mesotrione，MolDock评分：-118.689 kcal/mol）的结合亲和力（MolDock评分：-166.256 kcal/mol），这归因于有利的氢键、疏水作用、π-π相互作用和范德华相互作用。MM-GBSA计算进一步证实了蛋白质-配体复合物的热力学稳定性，其中化合物2在所有评估化合物中表现出最优的结合自由能。ADMET、毒性和类除草剂分析表明其具有理想的农用化学品特性和良好的预测安全性。此外，100 ns的MD模拟揭示，化合物2-HPPD复合物相对于硝磺草酮表现出更高的结构稳定性、持久的分子间相互作用和改善的构象完整性。DFT分析表明其具有增强的电子稳定性和适当的化学反应活性，支持在HPPD活性位点内的高效分子识别。综合而言，对接、MM-GBSA、MD、DFT和ADMET分析之间的高度一致性将化合物2确定为该系列中最有前景的先导化合物。这些发现凸显了取代脲衍生物作为合理设计下一代HPPD抑制剂的重要分子骨架，并为未来的实验验证和除草剂开发提供了坚实的计算基础。
取代脲衍生物作为HPPD抑制剂的计算研究：从分子识别到先导化合物鉴定
一、研究背景与意义
化学除草在现代农业生产体系中占据着不可替代的核心地位，这主要归因于其稳定的田间表现、操作效率以及对大规模农业系统的适配性。尽管包括利用根际细菌的生物防治在内的多种环境友好型杂草防除策略已被探索，但这些替代方法在实际应用中往往难以达到常规田间使用所需的稳定性和可操作性。在各类化学除草剂中，基于脲的化合物凭借其灵活可调的化学结构和潜在的选择性生物活性，作为有效的杂草治理工具引起了广泛关注。然而，杂草种群中除草剂抗性的持续演化构成了重大挑战，凸显了持续开发新型高效除草剂的迫切需求。
4-羟苯基丙酮酸双加氧酶（HPPD）是植物体必需的酶类之一，催化酪氨酸分解代谢途径中对羟苯基丙酮酸转化为尿黑酸的反应。抑制该酶会阻断质体醌和生育酚的生物合成，进而干扰类胡萝卜素的生成，最终导致叶片白化（褪绿）和植株死亡。商品化除草剂硝磺草酮（Mesotrione）正是通过这一机制发挥除草功效。然而，随着除草剂抗性问题的日益严峻，开发结构新颖、作用机制明确且环境相容性更优的HPPD抑制剂成为农药化学研究的重要方向。
计算方法的快速发展为除草剂研发提供了变革性工具。分子对接、定量构效关系（QSAR）分析和分子动力学（MD）模拟等方法使研究人员能够在实验合成之前对候选化合物进行系统的虚拟筛选和优化评估。这些计算机模拟策略显著减少了对耗时实验筛选的依赖，加速了具有理想理化性质和生物活性的候选分子的识别进程。基于此，本研究旨在通过整合多种计算方法的联合工作流程，系统评估一系列取代脲衍生物的HPPD抑制潜力，并将其与标准除草剂硝磺草酮进行比较分析。
二、研究方法概述
研究采用多重计算策略整合的工作流程。数据集来源于文献中已报道的41种取代脲衍生物，其抑制活性以半数抑制浓度（IC??）表示并转换为pIC??值进行标准化。靶蛋白结构从蛋白质数据银行（PDB ID: 1TFZ）获取。主要技术方法包括：（1）利用Molegro Virtual Docker软件进行分子对接研究，分析配体与HPPD活性位点的结合模式及相互作用；（2）采用Schr?dinger套件的Maestro 14.3进行MM-GBSA结合自由能计算，量化蛋白质-配体复合物的热力学稳定性；（3）通过GROMACS 2024进行100 ns分子动力学模拟，评估复合物的动态行为和结构稳定性；（4）运用Spartan '14软件在B3LYP/6-31G水平上进行密度泛函理论（DFT）计算，获取前沿分子轨道能量和全局量子化学描述符；（5）利用SwissADMET和ADMETlab 3.0平台预测化合物的吸收、分布、代谢、排泄和毒性（ADMET）特征。
三、主要研究结果
分子对接结果
分子对接分析揭示了41种化合物在HPPD酶活性位点内的结合亲和力存在显著差异。参考除草剂硝磺草酮的MolDock评分为-118.689 kcal/mol。多种化合物表现出更强的预测结合相互作用，最突出的是化合物2（-166.256 kcal/mol）、36（-159.167 kcal/mol）、32（-154.060 kcal/mol）和23（-151.845 kcal/mol）。这些显著更低的能量值表明这些配体在酶结合口袋内具有增强的结合亲和力和改善的构象稳定性。氢键相互作用分析显示，化合物2和21表现出更强的氢键贡献，而化合物32和36虽仅表现出中等的氢键能量，却获得了优异的整体对接评分，表明疏水力和范德华吸引力等其他非共价相互作用对稳定配体-酶复合物同样具有重要贡献。从结构角度，磺酰脲键、杂环系统和卤素及三氟甲基等吸电子基团的存在可能改善了分子在HPPD活性位点的匹配度。
交互作用分析表明，参考配体硝磺草酮在HPPD活性位点内采取紧凑稳定的构象，其β-三酮部分定位于催化中心，与GLN286形成常规氢键，并与催化Fe2?离子配位——这对β-三酮类除草剂的抑制活性至关重要。π-π堆积（PHE360）和涉及PHE398、LEU244、VAL248的疏水接触进一步增强了结合。相比之下，配体2通过与GLN358、SER246、ASN261和PHE371的常规氢键相互作用，以及与ILE273、ALA247、ILE260、LEU244和GLY399的范德华相互作用稳定结合。π-π T型相互作用和与PHE360、PHE403及MET314的烷基接触将配体稳定在疏水区域。尽管与PHE398存在轻微的不利受体-受体相互作用，但这一局部且较弱的作用并未损害复合物的整体稳定性和结合亲和力。
MM-GBSA结果
MM-GBSA结合自由能计算为对接结果提供了热力学验证。所有五个得分最高的配体（2、9、23、32和36）均表现出负的结合自由能，范围从-41.75到-60.01 kcal/mol，证实与HPPD受体的复合物形成在热力学上是有利的。配体2展现出最优的结合自由能（-60.01 kcal/mol），紧随其后的是配体23（-59.89 kcal/mol），两者仅相差0.12 kcal/mol，表明这两种化合物具有几乎相同的结合强度。值得注意的是，与文献报道的硝磺草酮MM-GBSA结合自由能（-28.34 kcal/mol）相比，所有五种选定的脲衍生物均表现出更优的热力学参数。有趣的是，对接评分与MM-GBSA结果之间存在一定差异：配体9和36获得最优对接评分但并非最低结合自由能，而配体2和23虽对接评分适中却展现出最强的MM-GBSA结合能量，凸显了两种方法互补评估的必要性。
ADMET与安全性评估
ADMET预测结果显示，所选衍生物的理化性质大体落在除草剂和药物相似性的公认范围内。分子量（401.06-475.07 g/mol）低于500 g/mol的上限，氢键供体数（2-3）在推荐限值内，Log P值（0.799-2.835）反映适度的亲脂性，有利于平衡膜渗透性和环境迁移性。化合物33和39完全满足Lipinski五规则，而化合物2存在两项违规（主要为氢键受体数较高和拓扑极性表面积较大），但其整体计算性能最优。毒性预测方面，所有评估的脲衍生物均被预测为非致癌性，且未显示皮肤致敏、眼刺激或呼吸毒性。相比之下，硝磺草酮在多个毒性终点呈阳性预测，包括致癌性、皮肤致敏、呼吸毒性和AMES致突变性。尽管设计化合物和参考除草剂均显示人肝毒性和血液毒性迹象，但脲衍生物缺乏AMES相关致突变性尤为值得关注。
分子动力学模拟分析
RMSD分析：两种复合物在100 ns模拟中均达到平衡并保持稳定，骨架RMSD值主要介于1.9-2.5 ?之间。1TFZ-化合物2复合物表现出特别一致的RMSD模式，波动幅度更小、轨迹更平滑，表明该配体与HPPD结合口袋具有高度相容性。
RMSF分析：两种复合物的波动模式总体相当，最大波动发生在N端区域。然而，硝磺草酮结合系统在残基180-250和C端区域附近表现出略高的波动（最高约4.0-4.4 ?），而化合物2复合物相应区域波动较低（约2.2-3.1 ?），表明化合物2通过限制酶的过度构象变化贡献了更稳定的蛋白质环境。
SASA分析：配体结合复合物的SASA值高于脱辅基蛋白。1TFZ-化合物2复合物记录到最高SASA值（约268-278 nm2），但轨迹在整个模拟过程中保持高度一致，无突变变化，表明配体结合诱导的结构重排是稳定性的体现而非结构完整性丧失。
半径回旋半径分析：化合物2复合物的Rg值最大（约23.1-23.6 ?），但轨迹高度一致，无重大结构破坏证据。硝磺草酮复合物的Rg值较低（约21.4-22.1 ?），在55-65 ns和75-85 ns出现中度增加但迅速恢复平衡。
氢键分析：硝磺草酮复合物在模拟早期维持约2-3个氢键，约35-42 ns达最大值，但45 ns后氢键占据显著减少。相比之下，化合物2复合物虽然初始阶段氢键形成有限，但随着模拟推进逐渐建立更一致的相互作用，30-50 ns维持约1-2个氢键，约80 ns后持续表现出一至多个以上氢键，反映出逐渐稳定的结合模式。
DFT分析结果
DFT计算显示，硝磺草酮的HOMO能量（-0.26859 a.u.）略低于配体2（-0.25270 a.u.），表明其供电子倾向稍低。配体2的HOMO-LUMO能隙（0.16071 a.u.）大于硝磺草酮（0.13051 a.u.），暗示配体2具有更高的动力学稳定性、更低的化学反应活性及更强的抗电子扰动能力。配体2的化学硬度（0.080355 a.u.）高于硝磺草酮（0.065255 a.u.），与硝磺草酮表现出更高的化学软度（15.3245 vs. 12.44478）一致。两者的亲电指数均非常低，表明亲电性不太可能是它们结合行为的主要控制因素。
四、讨论与结论
讨论部分指出，综合多种计算方法的评估结果具有高度一致性。分子对接识别出结合模式优于参考化合物的多个候选分子，MM-GBSA计算确认了这些复合物的热力学稳定性，MD模拟揭示了化合物2在长期动态条件下的稳定结合特征，DFT分析提供了电子层面的稳定性依据，ADMET预测则支持其良好的安全性。尽管化合物2在Lipinski五规则上存在两项违规，但考虑到其整体优异的计算性能，这一偏差不构成实质性障碍。研究还分析了取代脲骨架、杂环系统和吸电子取代基的组合效应在优化结合稳定性和抑制潜力中的协同作用。
结论部分明确指出，本研究通过整合分子对接、MM-GBSA结合自由能计算、分子动力学模拟、DFT分析和ADMET/毒性分析的计算策略，系统评估了设计化合物的HPPD抑制潜力。在所有研究的衍生物中，化合物2始终展现出最佳综合性能，表现为强结合亲和力、高度有利的结合自由能、稳定的蛋白质-配体相互作用、良好的理化与药代动力学性质，以及与参考除草剂硝磺草酮相比更优的预测安全性。分子动力学模拟通过稳定的构象行为和持续的分子间相互作用确认了化合物2-HPPD复合物的结构稳定性，MM-GBSA分析验证了其热力学有利的结合。DFT计算揭示了增强的电子稳定性和适当的化学反应活性，表明化合物2的结构架构促进了HPPD活性位点内的高效分子识别。替代脲骨架、杂环框架和吸电子取代基的组合效应在优化结合稳定性和抑制潜力方面发挥了关键作用。该研究为取代脲骨架作为合理除草剂设计的有效模板提供了坚实的计算基础，并为有效、选择性和环境可持续的HPPD靶向除草剂的开发以及应对现代农业中日益严峻的除草剂抗性挑战提供了重要理论依据。这些发现仍需通过合成、酶学测定和温室或田间评估进行实验验证。
