在农业害虫防治的漫长博弈中，传统神经毒性杀虫剂正逐渐陷入“杀敌一千，自损八百”的困局。这不仅体现在对蜜蜂、家蚕等非靶标生物的广谱毒性上，更体现在害虫抗性的爆发式增长。目前全球小菜蛾种群对氯虫苯甲酰胺的抗性已飙升至数千倍，2019年侵入我国的草地贪夜蛾更是让多地玉米减产超30%。在这样的背景下，昆虫生长调节剂（IGR，Insect growth regulator）凭借“靶向昆虫独有生长发育过程、对人畜低毒、抗性风险低”的独特优势，成为绿色农药研发的核心赛道。然而，IGR的研发绝非坦途，如何在提升活性的同时降低成本、并克服单靶标药剂易产生抗性的问题，是当前科研痛点。
引言
IGR的作用逻辑非常聪明，其不直接毒杀成虫，而是精准打击昆虫蜕皮、变态发育过程中的关键蛋白，让害虫在蜕皮时卡壳、畸形，最终无法正常发育为成虫或蛹。目前研究最成熟的两类IGR靶标，都与昆虫蜕皮过程直接相关。
第一类是蜕皮激素受体（EcR，ecdysone receptor）。蜕皮激素是昆虫蜕皮的总开关，而EcR则是接收这个信号的信号接收器。非甾体型二苯甲酰肼类（DBH，dibenzoylhydrazine）IGR正是模拟蜕皮激素的结构，抢占EcR的结合位点，导致昆虫内分泌紊乱，提前蜕皮并最终死亡。目前已商品化的DBH类IGR结构如下图所示，其中虫酰肼（Tebufenozide，112410-23-8）是应用最广的代表。然而，这类药剂结构骨架单一，长期大面积使用后，小菜蛾、甜菜夜蛾等害虫已对其产生高水平抗性，应用空间被不断压缩。
图1 已商品化的非甾体型二苯甲酰肼类（DBH）昆虫生长调节剂结构
第二类是几丁质酶（Chitinase）。几丁质是昆虫表皮和围食膜的核心成分，相当于昆虫的盔甲。昆虫每次蜕皮都需要几丁质酶来降解旧表皮、合成新表皮。如果能抑制几丁质酶的活性，害虫就会因无法完成蜕皮而死亡。目前虽有茶碱（Theophylline）、咖啡因（Caffeine）、己酰胍（Hexanoylguanidine）等几丁质酶抑制剂的报道，但至今没有任何一款几丁质抑制剂类IGR获批商品化，要么活性不足，要么选择性差，对作物和环境存在潜在副作用。
为突破单靶标IGR的局限性，中国农业大学的科研团队此前通过骨架跃迁（Scaffold hopping）策略，发现了首个能同时作用于EcR和亚洲玉米螟几丁质酶（OfChtI，Ostrinia furnacalischitinase）的双靶标先导化合物D-27。该化合物是七元环的吡唑酰胺结构，虽然对OfChtI的抑制活性优异，但对EcR的结合活性偏弱。更棘手的是，其合成原料环庚酮价格昂贵，关键中间体收率仅为50%，原子经济性差，难以规模化应用（笔者认为这并不是问题，合成大咖多了去了）。此外，团队此前并不清楚连接酰胺键和苯环的桥链长度，究竟会如何影响双靶标结合活性？能不能通过简化结构、优化桥链，得到活性更好、成本更低的候选化合物？
带着这些问题，团队在D-27的基础上开展了本次研究。设计策略如Scheme 1所示：将七元环替换为更廉价、更易合成的六元环，同时在酰胺键和取代苯环之间插入1~2个亚甲基柔性间隔基。最终合成了24个六环吡唑酰胺衍生物（其中23个为新化合物），系统评价了它们的杀虫活性、双靶标结合能力，号称找到了活性显著优于D-27的潜力候选化合物6j，为双靶标IGR的开发提供了扎实的实验依据。Scheme 1 六环吡唑酰胺类衍生物的设计策略示意图
1. 合成路线
本次合成的核心思路是“降本+增效”，完整路线如Scheme 2所示。团队摒弃了昂贵的环庚酮，改用价格仅为其1/5的环己酮作为起始原料，通过克莱森缩合、环化、皂化、酰胺缩合四步核心反应，构建目标化合物。Scheme 2 目标六环吡唑酰胺类衍生物的合成路线
原料的替换带来了显著的效益提升。合成关键中间体3（2-氧代-2-(2-氧代环己基)乙酸乙酯）的收率达75%，远高于此前七元环中间体的68%。随后的环化反应中，关键中间体4（2-苯基-4,5,6,7-四氢-2H-吲唑-5-羧酸乙酯）收率为58%，也高于七元环中间体的50%。最终目标化合物的酰胺缩合收率普遍比对应七元环化合物高5%~10%。全流程不仅原料成本低廉，且原子经济性显著提升，非常适合后续的产业化放大。所有目标化合物的结构均通过¹H NMR、¹³C NMR和高分辨质谱（HRMS）完成了确证。
2. 生物活性
研究团队选取农业生产中危害最严重的三种鳞翅目害虫，小菜蛾（P. xylostella）、草地贪夜蛾（S. frugiperda）、亚洲玉米螟（O. furnacalis），按照行业标准开展了系统杀虫活性测试，并以已商品化的IGR标杆产品虫酰肼作为阳性对照。实验中用“n”表示酰胺键和取代苯环之间的亚甲基间隔基数，n=0表示无间隔基，n=1表示1个亚甲基，n=2表示2个亚甲基，测试数据汇总如表1。
表1 | 目标化合物对三种鳞翅目农业害虫的杀虫活性（校正死亡率，%）
2.1 对小菜蛾的活性
在500 mg/L浓度下，多数目标化合物都表现出中等至良好的杀虫活性。其中最亮眼的是6j（n=1，苯环对位溴取代），其校正死亡率达100%，与虫酰肼的效果相当，且远高于先导化合物D-27的65.8%。当浓度降低到200 mg/L时，6j的死亡率仍达40.3%，是所有测试化合物中最高的。
构效关系分析揭示了有趣的现象。将七元环改为六元环并未显著提升活性（如图2所示，六元的6f和七元的D-27活性几乎一致），但桥链长度对活性影响极大。当苯环取代基为溴、乙基、甲氧基时，带1个亚甲基间隔基（n=1）的化合物活性远高于无间隔基（n=0）或带2个亚甲基（n=2）的同类物。例如，6j（n=1，4-Br）的活性远高于6b（n=0，4-Br）和6u（n=2，4-Br）。此外，苯环对位取代的化合物活性普遍高于间位取代，如对位三氟甲氧基取代的6l活性远高于间位取代的6k。
图2 六元环与七元环先导化合物的杀虫活性对比
桥链亚甲基数目对小菜蛾杀虫活性的具体影响规律如图3(a)所示。值得注意的是，处理后的小菜蛾幼虫出现了典型的蜕皮异常症状：体色发黑、后肠脱出、表皮粗糙、无法正常蜕皮，与虫酰肼处理后的症状类似。但细微之处见真章，6j处理的幼虫头部大小和颜色与虫酰肼处理的存在差异（图3b），这提示6j的作用机制可能与传统DBH类IGR有所不同。
图3 (a) 桥链亚甲基数目对小菜蛾杀虫活性的影响；(b) 6j、虫酰肼处理后小菜蛾幼虫的症状对比
2.2 对草地贪夜蛾的活性
草地贪夜蛾是2019年侵入我国的入侵害虫，对玉米等作物危害极大。本次测试中，所有目标化合物都表现出一定的杀虫活性，其中6h、6i、6j、6t在500 mg/L浓度下的死亡率为68.7%~83.3%，远高于D-27的29.3%。6j的活性依然位居前列：500 mg/L浓度下死亡率达83.3%，200 mg/L浓度下死亡率达40.3%。
构效关系同样显示n=1的化合物活性最优（图4a），且苯环对位连吸电子取代基（如氟、溴）的化合物活性，显著高于连给电子取代基（如甲基、乙基、甲氧基）的同类物。例如，对氟取代的6h活性远高于对甲基取代的6n，对溴取代的6j活性远高于对乙基取代的6o。处理后的草地贪夜蛾幼虫生长发育被显著抑制，同样出现蜕皮异常、畸形的症状，但其体色与虫酰肼处理的个体有明显差异（图4b），进一步印证了6j具有独特的作用模式。
图4 (a) 桥链亚甲基数目对草地贪夜蛾杀虫活性的影响；(b) 6j、虫酰肼处理后草地贪夜蛾幼虫的症状对比
2.3 对亚洲玉米螟的活性
在200 mg/L浓度下，多数目标化合物对亚洲玉米螟的活性偏低，但6j依然表现出33.5%的死亡率，远高于D-27的11.8%。桥链亚甲基数目对亚洲玉米螟杀虫活性的影响规律见图5a。处理后的幼虫同样出现后肠脱出、表皮褶皱等蜕皮异常症状，只是体毛颜色比虫酰肼处理的个体更深（图5b）。
综合三类害虫的测试结果，6j是本次研究中综合活性最优的化合物，不仅活性高，更重要的是，其引发的生理反应与传统药剂既有相似之处又有细微差别，这暗示了其双靶标作用机制的独特性。
图5 (a) 桥链亚甲基数目对亚洲玉米螟杀虫活性的影响；(b) 6j、虫酰肼处理后亚洲玉米螟幼虫的症状对比
3. 靶标验证
6j处理后害虫出现的蜕皮异常症状，与D-27、虫酰肼的症状高度相似，但这仅是表象。为从分子水平验证其作用机制，团队开展了蛋白水平的结合/抑制实验，结果如下图所示：
图6 | (a) 6j、D-27、虫酰肼对Plutella xylostella EcR的结合活性；图6(b) 6j、D-27、(GlcN)₅对OfChtI的抑制活性
EcR结合实验：以氚标记的蜕皮激素PonA为竞争配体，测试化合物对Plutella xylostellaEcR/USP异源二聚体的结合能力。结果显示，40 mg/L浓度下，6j的结合率达85.5%，远高于D-27的32.5%，仅略低于虫酰肼的93.7%。这表明6j对EcR的亲和力极强，是其发挥杀虫活性的重要基础。
OfChtI抑制实验：以4-甲基伞形酮-β-D-N,N'-二乙酰壳二糖苷为底物，测试化合物对OfChtI的抑制活性。结果显示，40 mg/L浓度下，6j的抑制率达66.0%，与D-27的67.6%相当，仅略低于阳性对照五-N-乙酰氨基葡萄糖（(GlcN)₅）的75.6%。这说明6j对OfChtI的抑制能力同样出色。
两项实验共同证实：6j确实可以同时靶向EcR和OfChtI，是典型的双靶标IGR候选物。这种一把钥匙开两把锁的特性，既能通过干扰蜕皮信号通路，又能通过阻碍表皮降解，双重机制协同作用，不仅能提升杀虫效率，更能大幅降低害虫产生抗性的风险。
4. 分子对接
为从原子层面解析“为什么六元环比七元环活性差不多？为什么1个亚甲基的间隔基活性最好？”等问题，团队开展了分子对接研究，将微观结构差异与宏观活性数据联系起来。
4.1 与EcR的结合模式
将6j、D-27和虫酰肼分别对接到棉铃虫EcR/USP异源二聚体的晶体结构（PDB ID：3IXP）中，发现三者的结合模式高度相似，都通过酰胺键的羰基氧与EcR的关键氨基酸Tyr408、Asn504形成氢键。然而，氢键的长度是决定结合强度的关键。对接结果如下图所示，6j的羰基氧与Tyr408的氢键距离为3.0 Å，短于D-27的3.2 Å，说明结合更紧密，这完美解释了6j对EcR的结合活性高于D-27的原因。而虫酰肼的氢键距离仅为2.8 Å，结合更紧密，因此活性略高于6j。
图7 6j(a)、D-27(b)、虫酰肼(c)与EcR的结合模式；6j酰胺羰基与Tyr408的氢键距离(d，3.0 Å)
对于相同苯环取代基的化合物，桥链长度会直接影响氢键的形成。n=0的6b和n=2的6u都只能与Asn504形成氢键，无法与Tyr408形成氢键，因此结合活性远低于n=1的6j。这从分子动力学角度解释了为何n=1是活性最优的桥链长度。
4.2 与OfChtI的结合模式
将6j、D-27和底物(GlcN)₅分别对接到亚洲玉米螟I组几丁质酶的晶体结构（PDB ID：3WL1）中，结合模式叠加结果如图8(a) 所示。6j和D-27都能与关键氨基酸Trp107形成氢键，也能与Trp34、Phe309、Tyr272形成π-π相互作用。但6j有一个独特的优势，即可以与另一个关键氨基酸Trp372形成额外的π-π相互作用（图8d），而D-27不具备这个相互作用。这一额外的分子间作用力，是6j对OfChtI抑制活性与D-27相当甚至更优的重要原因。
图8 (a) 6j、D-27、底物(GlcN)₅与OfChtI的结合模式叠加；图8(b) D-27与OfChtI的相互作用；图8(c) 6j与OfChtI的相互作用；图8(d) 6j与Trp372的额外π-π相互作用
同样，桥链长度会影响与疏水残基的相互作用。n=1的6j与疏水残基Phe309、Tyr272的相互作用，远强于n=0的6b和n=2的6u，这也从分子层面解释了为什么1个亚甲基间隔基的化合物活性最优。
分子对接的结果完美呼应了生物活性的构效关系，为后续的化合物优化提供了明确的方向。后续可以通过调整苯环取代基的电子效应、优化间隔基的柔性和长度，进一步提升6j的双靶标结合活性。
结论与展望
这项研究通过“七元环改六元环+插入柔性亚甲基间隔基”的策略，成功开发出了24个六环吡唑酰胺类衍生物。其中化合物6j的综合杀虫活性优于先导D-27，且同时具备EcR结合活性和OfChtI抑制活性，是极具潜力的双靶标IGR先导化合物。更重要的是，研究明确了“适当增加分子柔性可以提升双靶标结合活性”的规律，为后续绿色IGR的创制提供了可复制的设计思路。