《Molecules》：QSAR-Based Ecotoxicological Assessment of Novel Imidazole Derivatives for Sustainable Crop Protection
编辑推荐：
咪唑已被证实是非常有效的农药，尤其对抗植物病原真菌和昆虫。由于其对环境的不利影响，仅有少数咪唑类化合物获欧洲委员会（EC）批准使用。亟需开发对多种害虫高效、广谱，同时对环境和有益生物及人类安全的新咪唑衍生物。为减少昂贵且耗时的实验，研究人员采用基于定量构效关系
咪唑已被证实是非常有效的农药，尤其对抗植物病原真菌和昆虫。由于其对环境的不利影响，仅有少数咪唑类化合物获欧洲委员会（EC）批准使用。亟需开发对多种害虫高效、广谱，同时对环境和有益生物及人类安全的新咪唑衍生物。为减少昂贵且耗时的实验，研究人员采用基于定量构效关系（QSAR）模型的计算机（in silico）方法预测新化合物或未测试化学品的毒性。本研究使用Vega与ChemFREE网络平台评价拟合成咪唑衍生物的农药相似性、环境风险属性与生态毒理效应。研究人员对金刚烷（adamantane）、烷基（alkyl）及三唑（triazole）类的酰胺（amide）、酯（ester）、氨基甲酸酯（carbamate）和酮（ketone）衍生物按上述属性筛选，其中四个烷基酰胺被凸显为潜在高效且环境安全的抗真菌、除草与杀虫剂。分子对接（molecular docking）研究指示了上述化合物抗真菌、除草、杀虫及抗菌活性的可能作用机制，并揭示与特定受体结合相关的结构特征。
研究背景方面，合成化合物在植物保护中限制了病虫害发生与发展，但害虫抗药性、农药对环境和健康危害促使需要新的植物保护产品活性成分。这类化合物须具有高特异性和环境及毒理可接受性。农药生态效应受理化性质、施用量、暴露时间和环境条件影响，其中土壤理化性质尤为关键。持久性指农药在土壤、水、空气或沉积物中保持至发生化学或生物转化前的趋势，微生物酶可介导降解。化合物结构同时影响持久性与生物降解，因此这些终点适合用基于结构的计算方法预测。现代作物保护要求化合物兼具高生物活性和对非靶标生物低毒，欧盟Regulation (EC) No. 1107/2009引入更严格准入标准，并推动非动物实验。农药实验室发现与商业化成本高、周期长，计算机辅助分子设计（CAMD）因高效被广泛用于农用化学品生态毒理建模与设计，欧洲化学品管理局（ECHA）的REACH指南允许用QSAR模型预测危险属性以避免动物试验。咪唑为五元杂环芳香环、1和3位有两个氮原子，具广谱药理活性，农用上咪唑类杀菌剂较成熟，但欧盟再评估后仅少数获批，如imazalil仍可用于特定用途，prochloraz因健康与环境安全未续批。咪唑类仍被广泛研究抗真菌、抗菌、抗病毒、杀虫和线虫活性，咪唑鎓离子液体对Fusarium和Alternaria有强抗真菌效应。现有咪唑环境毒理研究少，imazalil慢性暴露可对蚯蚓Eisenia andrei致亚致死繁殖、细胞毒性和氧化应激效应，高浓度咪唑离子液体抑制植物萌发与根茎生长。土壤微生物维持有机质矿化与养分循环，破坏其多样性会影响土壤健康和作物生产力，欧盟Green Deal与Farm to Fork Strategy也主张减少常规化学农药。因此研究人员以32个咪唑衍生物（烷基、金刚烷、三唑的酰胺、酯、氨基甲酸酯、酮）为对象，用QSAR与分子对接在合成前评估农药相似性、环境风险和生态毒理，并探明作用机制，论文发表于《Molecules》。
关键技术方法方面，研究人员将32个设计咪唑衍生物结构转为SMILES作为输入，使用ChemFREE平台基于PaDEL-Descriptor与图卷积/图注意力卷积神经网络（GACNN）预测农药相似性（以RDL表示概率）、生物降解性（28天降解60%为易生物降解）及水、土、沉积物、空气中降解半衰期（HL）；用VEGA平台调用鱼类急性毒性、胖头鱼短期毒性、大型溞（Daphnia magna）短期毒性、藻类毒性分类与急性毒性、蜜蜂急性毒性、蚯蚓毒性、活性污泥毒性、生物富集因子（BCF）等QSAR模型，并以ALOGPS 2.1算亲脂性。分子对接前用Spartan ’08以MM+力场与PM3半经验量子法优化配体，对蛋白激酶K（PDB ID: 2PWB）、甾醇14α-脱甲基酶（CYP51，PDB ID: 5EAH）、几丁质酶（PDB ID: 4TXE）、谷氨酸门控氯通道（CGluCl，PDB ID: 3RIF）、果蝇多巴胺转运体（dDAT，PDB ID: 7WGT）、乙酰辅酶A羧化酶羧基转移酶域（PDB ID: 3K8X）、肽脱甲酰酶（PDB ID: 1G2A）用iGEMDOCK做对接，以重新对接根均方偏差（RMSD）<2 ?验证方法可靠性，BIOVIA Discovery Studio可视化相互作用。
Results部分，研究人员展示咪唑衍生物农药相似性与环境风险属性，以及生态毒理预测值，并标明预测可靠性G/M/L。通过RDL筛选发现金刚烷衍生物6、7、16和烷基18、27、28、29、30抗真菌概率最高；除草以金刚烷8和烷基12–16、18、27–32最佳；杀虫以金刚烷9和烷基13–15、27、31最佳。但金刚烷类难生物降解被初筛排除。水生动物高毒预测集中于金刚烷乙酰胺、碘盐羧胺和醋酸酯类，藻类高毒与氯苯基（chlorophenyl）片段21–24有关。蜜蜂毒性低至中等，三唑基尤其1,2,3-三唑19增强蚯蚓毒性的概率。污泥毒性和BCF与亲脂性相关，多数VEGA预测可靠性低需谨慎。经环境风险与生态毒理过滤，抗真菌剩烷基酰胺29、30；环保除草剂为29、30、31、32；31为领先杀虫剂。雷达图显示候选满足农药相似性（可旋转键≤9、氢键供体≤2、受体≤6、分子量≤435、AlogP≤6），仅30有3个氢键供体、穿透蜡质角质层能力弱。分子对接表显示29、30与蛋白激酶K和几丁质酶结合能相近，30与蛋白激酶K为?87.69 kcal mol?1、几丁质酶?86.69 kcal mol?1；31与CGluCl为?72.09 kcal mol?1、GABA转运体同工型?74.77 kcal mol?1；30与羧化酶?56.00 kcal mol?1、肽脱甲酰酶?83.63 kcal mol?1，RMSD均<2。共结晶配体如38F、ivermectin（IVM）、NO711（9BC）、tepraloxydim（B89）、actinonin（BB2）结合能更低属预期。图4显示30进入蛋白激酶K、羧化酶、肽脱甲酰酶位点，31进入CGluCl/GABA转运体位点，相互作用含氢键、范德华、π–σ、π–π堆积、π–烷基。
Discussion部分，研究人员指出金刚烷咪唑有高抗真菌亲和性但难降解、亲脂导致水生高毒；三唑类抑菌但具斑马鱼发育毒；羧胺类亲水致水生毒；烷基咪唑适配为安全农药骨架。对接机制上，30抑制蛋白激酶K靠2-羟基与Asn161/Thr223/Ser224成氢键，咪唑氮与Leu133/Ala158/Ser170成氢键，芳环与Gly134 π–π、Ala172 π–π堆积。31入NO711位点，酰胺氮与Tyr124氢键，咪唑环与Tyr43/Ser327范德华，丙基与Ala48 π–烷基。30抑乙酰辅酶A羧化酶羧基转移酶（CT）域靠吡啶型氮与羧胺氮双氢键Ala1627，羟基氢键Val1733/Ser1625，咪唑环π–烷基Ile1629。30抑肽脱甲酰酶（PDF）靠羟基与Gln50/Glu133、羧胺氧与Gly45/Ile44、吡咯型氮与Glu133氢键，咪唑π–σ与His132。研究人员也说明局限：ChemFREE适用性域检查受限，VEGA多数生态毒理预测可靠性低，分子对接仅为假设结合模式需实验确认。
Conclusions部分翻译：本研究是基于化学信息学方法理性开发植物保护产品的初始阶段，可在分子水平初步测试毒性和环境效应。一组设计的咪唑衍生物（金刚烷、烷基和三唑的酰胺、酯、氨基甲酸酯和酮）按估计的农药相似性、环境风险属性和生态毒理效应筛选，四个烷基酰胺脱颖而出成为可持续作物保护的有前景药剂。分子对接研究提示了抗真菌、抗菌、除草和杀虫活性的可能作用模式，并暗示与特定受体结合的重要结构特征。研究的下一步将包括“绿色”合成与针对病原及有益生物的生物测定。
该研究意义在于以in silico流程在合成前排除高环境风险结构、聚焦烷基酰胺候选，为减少动物实验和农用化学品开发成本提供可复用框架，并显示咪唑骨架经结构修饰可兼顾靶标活性和非靶标安全性，对《Molecules》所代表的交叉农化与计算毒理方向有方法学示范价值。
