《Biology》：Ligand-Specific Remodeling of the Anopheles gambiae RDL GABA Receptor by the A296S Mutation Differentiates Recognition of Fipronil and Fluralaner
编辑推荐：
γ-氨基丁酸(GABA)门控氯离子通道抗狄氏剂(RDL)受体是杀虫剂的主要分子靶点，然而与抗性相关的突变如何差异性影响化学结构不同的配体的结构基础尚未完全阐明。本研究利用同源建模、分子对接、200 ns分子动力学(MD)模拟（包括载脂蛋白受体及三个独立重复的每
γ-氨基丁酸(GABA)门控氯离子通道抗狄氏剂(RDL)受体是杀虫剂的主要分子靶点，然而与抗性相关的突变如何差异性影响化学结构不同的配体的结构基础尚未完全阐明。本研究利用同源建模、分子对接、200 ns分子动力学(MD)模拟（包括载脂蛋白受体及三个独立重复的每种配体结合复合物）以及包含相互作用熵(IE)项的MM/PBSA自由能计算，探究了经典的A296S置换对冈比亚按蚊(Anopheles gambiae)RDL受体与苯基吡唑类氟虫腈(fipronil)和异噁唑啉类氟雷拉纳(fluralaner)相互作用的影响。载脂蛋白状态下的分子动力学模拟显示，A296S保持了受体的整体结构稳定性，同时局部改变了孔道内衬M2跨膜螺旋的环境。尽管分子对接在两种受体变体中均对氟雷拉纳的打分优于氟虫腈，但重复的分子动力学模拟揭示了配体和基因型依赖的行为。在野生型受体中，氟雷拉纳保持稳定结合，其MM/PBSA自由能估算值为?18.3 ± 4.2 kcal/mol，而氟虫腈的配体RMSD逐渐增加，其估算值接近零(+1.2 ± 9.5 kcal/mol)。在A296S模型中，氟虫腈的估算值相似(+1.6 ± 4.3 kcal/mol)，但氟雷拉纳的估算值变为不利(+10.0 ± 1.8 kcal/mol)，且Ser296本身对两种配体的贡献均为轻微不利。综上所述，这些发现表明A296S并非作为配体识别的通用决定因素，而是与建模互作景观中配体特异性的变化相关联，其中对氟雷拉纳的预测效应更大。这些预测源自同源模型，未经电生理学或结合试验验证，也不衡量功能性抗性。尽管如此，它们可能有助于指导交叉抗性模式的解读以及针对GABA受体类杀虫剂的基于结构的研发工作。
一、研究背景与研究目的
疟疾仍然是全球公共卫生的重大负担，尤其是在撒哈拉以南非洲地区，冈比亚按蚊是恶性疟原虫的主要传播媒介。目前对蚊虫的防控高度依赖杀虫剂，而杀虫剂靶点产生抗药性构成了严重威胁。γ-氨基丁酸（GABA）门控氯离子通道中的抗狄氏剂（RDL）受体是昆虫中枢神经系统中主要的突触后GABA受体，也是多类杀虫剂的重要分子靶点。苯基吡唑类代表药物氟虫腈作为非竞争性拮抗剂（NCA）作用于该受体，而异噁唑啉类代表药物氟雷拉纳则被归类为GABA门控氯离子通道的别构调节剂。位于M2跨膜螺旋2′位置的丙氨酸到丝氨酸的保守置换（A296S）是一个经典抗性突变，最初在环戊二烯类抗性果蝇中发现，随后在多种蚊媒中被报道。该突变对环戊二烯类的抗性作用强且明确，但对氟虫腈敏感性的影响却表现出高度可变性，而对异噁唑啉类的影响证据相对有限。因此，一个关键的未解机制问题是：一个单一的保守置换如何对作用于部分重叠但非完全相同结合位点的不同配体产生显著不同的效应？表型抗性测定无法解释这种差异的分子基础。本研究旨在通过同源建模、分子对接、分子动力学（MD）模拟和MM/PBSA自由能计算等计算手段，比较A296S突变如何差异性地调控冈比亚按蚊RDL受体与氟虫腈和氟雷拉纳的相互作用，从而从原子层面揭示配体特异性效应的机制。
二、主要技术方法
研究人员首先利用Phyre2.2服务器以人β3δ GABAA受体冷冻电镜结构（PDB 7QN8，链D）为模板构建冈比亚按蚊RDL受体的单体同源模型，再通过GalaxyHomomer组装成五聚体，并经GalaxyRefine优化。随后在PyMOL中引入A296S突变。将野生型和突变型受体嵌入POPC磷脂双层膜中，使用AMBER FF19SB力场和GAFF2力场，通过GROMACS 2023进行200 ns的载脂蛋白状态分子动力学模拟。利用TTClust对轨迹进行构象聚类，选取代表性结构用于AutoDock Vina v1.2.7的分子对接。将打分最高的对接构象再次进行三组独立的200 ns配体结合态MD模拟。最后采用包含相互作用熵（IE）项的MM/PBSA方法计算结合自由能并进行残基能量分解，并使用HOLE程序分析孔道半径。
三、研究结果
3.1 结构建模揭示A296S置换后RDL整体架构得以保持
通过同源建模成功获得了野生型和A296S突变型受体的单体和五聚体结构，两者均展现出典型的半胱氨酸环配体门控离子通道超家族折叠。A296S的引入并未破坏整体结构，但丝氨酸侧链的额外羟基朝向孔道腔，导致预测的孔道直径相比野生型略有缩小。
3.2 载脂蛋白模拟揭示野生型和突变型受体全局动力学大体保守但存在局部差异
200 ns的载脂蛋白MD模拟显示，两种受体的骨架均方根偏差（RMSD）最终收敛至约3 ?，回转半径几乎相同，表明突变不影响整体结构稳定性。但A296S突变以链特异性方式重新分布了局部柔性，且突变受体平均分子内氢键数（400）略低于野生型（420）。
3.3 对接分析揭示配体特异性差异
分子对接显示，在所有受体构象中，氟雷拉纳的预测结合亲和力（Vina打分）均优于氟虫腈。在突变型受体中，氟雷拉纳与残基形成广泛的疏水接触、氢键和卤键网络，而氟虫腈的互作网络较为稀疏。在野生型受体中，氟雷拉纳同样比氟虫腈拥有更丰富的相互作用。
3.4 A296S突变型复合物的分子动力学和MM/PBSA分析
在三组重复的MD模拟中，两种配体结合的突变型受体骨架动态相似。配体RMSD值在最后150 ns内约为3.5–4.5 ?，表明两者均停留在结合区域。MM/PBSA计算表明，突变型-氟虫腈复合物的总结合自由能为+1.590 ± 4.316 kcal/mol，不具优势；而突变型-氟雷拉纳复合物则为+10.023 ± 1.848 kcal/mol，更为不利。残基分解显示，Ser296自身对两种配体均贡献轻微不利的能量。
3.5 分子动力学和MM/PBSA分析揭示野生型受体中氟虫腈和氟雷拉纳稳定性的差异
在野生型受体中，氟雷拉纳在整个模拟过程中保持稳定（配体RMSD约2 ?），而氟虫腈的平均配体RMSD逐渐增加至约6–7 ?，显示出更高的位置移动性。MM/PBSA分析确认氟雷拉纳的结合自由能非常有利（?18.317 ± 4.208 kcal/mol），而氟虫腈的估算值接近零（+1.177 ± 9.480 kcal/mol）。残基分解显示，Leu297对野生型-氟雷拉纳复合物的贡献最大（?2.663 ± 0.208 kcal/mol）。
3.5.3 跨基因型比较
综合来看，氟虫腈的自由能估算值在野生型（+1.177 ± 9.480 kcal/mol）和A296S突变型（+1.590 ± 4.316 kcal/mol）中相似。而对于氟雷拉纳，自由能估算值从野生型的?18.317 ± 4.208 kcal/mol变为突变型的+10.023 ± 1.848 kcal/mol，变化显著。这表明在该模型中，A296S与氟雷拉纳预测能量学的较大改变相关，而非氟虫腈。
四、讨论与结论
讨论部分指出，A296S并不作为一个统一的位阻或静电屏障起作用，而是局部重构了拮抗剂结合腔的物理化学环境，且这种改变被两种配体不对称地利用。载脂蛋白模拟证实突变不损害整体结构，这与A296S等位基因在田间种群中的持续存在一致。对于氟虫腈，其自由能估算值在两种基因型中无差异，与该突变对氟虫腈有限且可变的交叉抗性报道相符。对于氟雷拉纳，其在突变型中变得不利的估算值与异噁唑啉类占据部分重叠但非相同的跨膜腔这一药理特性并不直观吻合，可能反映了突变对局部口袋几何形状的间接效应、非等效的对接起始姿态或方法本身的局限性。研究强调，平衡结合估算值不能等同于功能敏感性或抗性预测，后者取决于状态依赖性通道阻断、门控动力学等多个因素。
结论翻译：
在所建立的RDL受体模型中，A296S置换与拮抗剂结合口袋的配体特异性变化相关联。该突变并未作为一个通用决定因子起作用，它未能检测性地改变氟虫腈微弱的MM/PBSA自由能估算值，且Ser296并未充当有利的能量锚点；而氟雷拉纳的估算值则从野生型受体中的有利（?18.3 ± 4.2 kcal/mol）转变为A296S受体中的不利（+10.0 ± 1.8 kcal/mol）。在野生型受体中，氟虫腈显示出比保持稳定的氟雷拉纳更大的位置移动性。这些发现表明，一个单一的保守孔道内衬置换可以对占据部分重叠结合环境的化学结构不同的杀虫剂产生不同的分子后果。这种静电和疏水接触的差异性重塑与报道的A296S对氟虫腈敏感性的可变表型效应相符；对于异噁唑啉类，表型证据有限，本研究中氟雷拉纳的结果对抗性而言尚无定论。这些计算分析描述了一个未经实验验证的同源模型的平衡互作景观，不应被解读为功能性抗性的直接度量。通过解析残基水平的能量贡献，本研究提供了A296S如何在RDL口袋差异性塑造结构不同杀虫剂识别的原子层面见解。研究结果并未确定异噁唑啉类是否对携带A296S的冈比亚按蚊种群保持活性，这需要电生理学或表型测试，但它们可能有助于指导未来针对RDL化合物的基于结构的研发工作，以评估其对经典2′孔道内衬残基的依赖性。
