等离子体衍生活性物种经植物水通道蛋白通道的选择性转运：分子动力学研究

《Scientific Reports》：Selective transport of plasma-derived reactive species through the plant aquaporin channels: a molecular dynamics study

　　

研究人员探究了冷大气等离子体（CAP）生成的活性氧氮物种（RONS）经植物水通道蛋白（AQP）的选择性通透规律，以识别可进入细胞内的等离子体衍生物种。采用全原子分子动力学（MD）与增强采样方法，研究人员量化了嵌入脂质双分子层中的PIP2;1水通道内八种RONS（HNO3、HO2、顺式-HNO2、反式-HNO2、N2O4、NO、NO2与O3）的自由能谱。通透行为受两个互补特征调控：芳香族/精氨酸（ar/R）收缩区入口处的自由能垒，以及孔道内部的自由能极小值深度。强极性、具氢键能力的HNO3呈现最高入口势垒（≈16 kJ·mol?1），故通透性最低；而弱极性物种如NO以极小阻力跨通道。值得注意的是，若干物种（最突出为N2O4）兼具低入口势垒与深内部极小值，表明其易被吸入孔道却在内部被强烈滞留而非自由转运。结果表明，PIP2;1对RONS的选择性由分子尺寸、极性与氢键能力共同决定，并为可持续等离子体农业技术中CAP生成物种的理性筛选提供了机制基础。

研究背景与立项依据

冷大气等离子体（CAP）技术在可持续农业中被认为是一种有前景的环保手段，可用于促进种子萌发、刺激植物生长与提高作物产量。CAP作用时会释放短寿命与长寿命混合的活性氧氮物种（RONS），包括H₂O₂、HO₂、·OH、NO、NO₂、O₃、N₂O₄、HNO₃及顺式/反式亚硝酸（HONO）等，外部施加于种子或组织后可能影响植物生理。实验已观察到等离子处理可加速萌发、促进早期发育、缩短成熟并减少灌溉需求，但其底层分子机制尚不清楚。一种假说认为等离子体RONS与细胞转运系统相互作用，而水通道蛋白（AQP，一类介导水及小溶质跨膜被动转运的整合膜蛋白）可能是RONS进入植物细胞的中介。植物中质膜内在蛋白（PIP）亚家族（尤其是PIP2;1）在水分稳态、营养吸收与胁迫适应中关键作用，其单体由六次跨膜螺旋及Asn-Pro-Ala（NPA）基序与芳香族/精氨酸（ar/R）选择性滤筛构成。已有研究提示部分AQP可转运H₂O₂等中性小溶质，但针对植物AQP传导CAP源RONS尚缺系统计算研究；既往Cordeiro与Yusupov等虽用MD考察过哺乳动物AQP1及相关脂质双层，但未专门面向植物促生长RONS转运。该论文发表于《Scientific Reports》，旨在从原子层面定量PIP2;1对典型CAP源RONS的被动通透自由能谱，为等离子体农业装置优化提供机制依据。

主要关键技术方法

研究人员以PDB ID 1Z98的四聚体PIP2;1为模型，嵌入单组分1,2-二油酰-sn-甘油-3-磷酸胆碱（DOPC）脂质双层，两侧显式水层并加Cl^?中和，采用GROMOS 53A6联合原子力场与SPC/E水模型，RONS参数取自Cordeiro等复现水-烷烃分配系数的GROMOS兼容集。先用GROMACS 2022做50,000步最速下降能量最小化，NPT系综下303.15 K（Nose-Hoover）与1.01325 bar（Parrinello-Rahman半各向同性）平衡，继以1 μs生产MD（2 fs步长）确认骨架RMSD约0.32–0.34 nm平台与膜厚稳定。通透自由能谱用伞状采样（US）：沿通道z轴加2000 kJ·mol^?1·nm^?2谐波偏置，xy平面0.5 nm平底限制；利用四孔并行，每孔8窗×11套错位0.1 nm组合成81窗（间距0.1 nm，跨z=?4至+4 nm），每窗5 ns（弃1 ns），总计2816个US模拟（≈14 μs），经加权直方图分析法（WHAM）周期边界重构自由能谱，bootstrap（50样本）估误差。

研究结果

全局结构稳定性：1 μs MD显示PIP2;1骨架RMSD于~80 ns后达0.32–0.04 nm平稳区，Z向盒长稳于~10 nm，跨膜螺旋残基RMSF<0.20 nm，仅N端、胞外环与C尾柔性高，导电通路刚性保持。

亲水RONS自由能谱：HO₂、顺式-HNO₂、反式-HNO₂与HNO₃均可低阻进入胞外前庭，但在ar/R区分化显著。HNO₃的ΔG_max≈16 kJ·mol^?1最高；HO₂与顺式-HNO₂几乎无上升势垒（ΔG_max≈0.1与1.0 kJ·mol^?1），反在孔内出现≈?11与?10 kJ·mol^?1深阱；反式-HNO₂居中（ΔG_max≈4.6 kJ·mol^?1伴浅阱）。说明强极性HNO₃因水相保留倾向被ar/R排斥，弱极性自由基则在极性孔衬残基处热力学稳定。

疏水/弱极性RONS自由能谱：NO谱线近平坦（ΔG_max≈3.5 kJ·mol^?1），近似自由扩散；O₃ΔG_max≈6.7 kJ·mol^?1伴?9.6 kJ·mol^?1深阱；NO₂ΔG_max≈8.3 kJ·mol^{?1；N₂O₄入口垒仅≈4 kJ·mol?1却落入≈?17 kJ·mol?1的全样本最深阱，表明其易入孔但被内部强滞留，归因于其体积大、四极矩几何与弱水合。}

势垒-阱深综合比较：ΔG_max排序为HNO₃ > NO₂ > O₃ > 反式-HNO₂ > N₂O₄ ≈ NO > 顺式-HNO₂ > HO₂。低ΔG_max不等于高通量，N₂O₄、HO₂、O₃、顺式-HNO₂皆“易进难出”。与无AQP的硝基氧化脂质膜相比，极性物种在PIP2;1中势垒更高，证实ar/R施加额外立体/静电选择。

讨论与结论翻译

讨论指出，疏水弱极性种（NO、O₃）可经PIP2;1高效入胞调控生长或抗逆；强极性种（HNO₃）受动力/立体阻碍，除非依赖孔瞬扩或脂质介导扩散否则胞内到达受限；CAP尾气组分可据此向高通透物种倾斜以放大农艺收益。结论部分表明：本研究通过全原子MD与US量化八种RONS在PIP2;1中的被动通透，发现选择性由ar/R入口垒与孔内极小值深度共同决定——强极性氢键种（最典型HNO₃）入口垒最高、最不易被动通透；弱极性种（如NO）几乎无罚穿越；N₂O₄、HO₂、O₃等易入却被深阱滞留而非自由转运。该“势垒排斥”与“深阱滞留”的二分法及二者对尺寸、极性、氢键能力的依赖，构成PIP2;1对RONS选择性的中心能量特征。上述结论限于中性被动非反应转运框架，未涵盖质子转移、氧化或自由基复合，也未含生理盐浓度与混合脂质/多亚型效应，但在该简化模型内为判别CAP物种经AQP入胞提供了预测框架，支撑可持续农业中等离子体技术的理性优化。