《Biology》：DUOX-Generated Intestinal ROS Restrict Spore Adhesion and Invasion of Nosema bombycis in Silkworm Bombyx mori
编辑推荐：
低频正态模分析（NMA）广泛用于预测蛋白质的集体运动，但直接笛卡尔缩放随着振幅增大会逐步扭曲α-碳（Cα）骨架。研究人员提出骨干角兼容尽可能刚性正态模分析（BA-ARAP-NMA），一种非线性Cα结构生成方法，在各向异性网络模型（ANM）计算中对相邻及次邻Cα
低频正态模分析（NMA）广泛用于预测蛋白质的集体运动，但直接笛卡尔缩放随着振幅增大会逐步扭曲α-碳（Cα）骨架。研究人员提出骨干角兼容尽可能刚性正态模分析（BA-ARAP-NMA），一种非线性Cα结构生成方法，在各向异性网络模型（ANM）计算中对相邻及次邻Cα距离的一阶变化施加约束，并在结构生成阶段应用共旋转尽可能刚性（ARAP）校正。在35对实验解析的双态蛋白结构上的测试表明，几何约束保留了大部分从输入态到第二态的位移信息；在35对中的34对中，重算约束模态使保留的过渡信息在主导低频模态中更集中。在匹配位移下，BA-ARAP-NMA在每一对中都大幅降低了局部Cα距离与角度扭曲。在蛋白对层面，BA-ARAP-NMA相较线性Cα-ANM通常给出更低的配对态根均方偏差（RMSD）和更高的过渡覆盖率，同时维持改进的虚拟角度精度。因此，BA-ARAP-NMA将Cα正态模结构生成扩展至直接线性缩放之外，在保留与过渡相关的集体信息的同时，提供带帧级几何度量的振幅有序结构，用于后续重建与精炼。
研究背景与问题提出
蛋白质功能依赖其特定三维构象，而构象变化中的铰链位点和柔性区域编码了功能态转换的动态信息。正态模分析（NMA）结合弹性网络模型（ENM）是解码集体运动的基础方法，其中Cα粗粒化各向异性网络模型（Cα-ANM）通过低频模态捕捉域级重排方向。然而标准NMA基于平衡构型附近的无穷小位移线性近似，直接笛卡尔缩放低频模态至数埃振幅时，会渐进破坏Cα虚拟键长与虚拟键角，导致局部骨架几何退化，阻碍全原子重建与后续模拟。已有非线性方案如刚体块旋转（NOLB）或内坐标表述（iMODS）虽缓解该问题，但或预设刚体块，或脱离笛卡尔输出。研究人员据此开展本研究，旨在不预设刚体块、保留Cα笛卡尔表示的前提下，在模态计算与结构生成两阶段同时引入局部几何控制，提出骨干角兼容尽可能刚性正态模分析（BA-ARAP-NMA）。论文发表于《Biology》。
主要关键技术方法
研究人员选取35对实验解析的双态蛋白结构对（覆盖域闭合、伴侣循环、交替接入转运、折叠切换、机械敏感通道门控等运动类，每对常见Cα残基80–1200个、序列一致度≥95%、态间Cα-RMSD≥3 ?），以15 ?为基准ANM截断（另测12/18 ?鲁棒性）。方法上：先移除六维刚体自由度，在模态计算前对连续段内相邻（i,i+1）与次邻（i,i+2）Cα距离施加一阶零变分约束，将ANM黑塞矩阵投影至约束零空间重对角化得BA约束模态；结构生成时沿单位Cα-RMS模态方向线性预测后，对局部边向量做正交Procrustes拟合得共旋转尽可能刚性（ARAP）修正，经组件质心规整的图拉普拉斯求解输出最终Cα坐标。以线性Cα-ANM为基线，设四路消融（线性缩放、线性+ARAP、BA约束+线性缩放、完整BA-ARAP-NMA），在0–20 ?实际位移网格上匹配实际Cα位移评估，并以NOLB作外部非线性对照（12对子集、等候选数、无配对结构引导）。统计单元为蛋白对，采用中位数效应与BCa自助CI及精确符号检验。
研究结果
3.1 局部Cα约束在35对基准中保留过渡相关信息
通过累积重叠分析，K=10时约束基在35对中的31对累积重叠高于线性Cα-ANM（对中位差+0.00952），约束空间保留中位94.3%实验位移，且34对中前导模态浓度（η10）提升（中位+0.0625）。表明约束重算模态并未剔除集体运动，而是将其重新集中于低频端，该性质在12/15/18 ?截断下稳定。
3.2 BA约束与ARAP对局部Cα几何贡献不同
四路消融在匹配实际位移1–10 ?下显示：完整BA-ARAP-NMA在每对、每位移均降低相邻Cα距离MAE与虚拟角MAE。BA约束模态独立贡献虚拟角控制（10 ?时中位降3.19°），ARAP独立贡献有限振幅邻边距校正（10 ?时完整法距MAE降0.665 ?，RMSD降0.527 ?）。二者互补而非互换。
3.3 局部几何控制改善蛋白对级生成结构集
匹配位移下完整法对配对态RMSD降低、过渡覆盖率提升，32–35对有利，BCa区间均居优侧；几何控制可达位移（无>1 ?键距异常、无>10°角异常）中位提升1.340 ?（BA约束单独贡献1.063 ?）。等量1600结构回溯最小RMSD中，完整法中位再降0.346 ?，30/35对有利，PPDK与MscS等5对未更低但不改局部几何优势。
3.4 HCV解旋酶示构级效应
以HCV解旋酶同模域为例，同一模态支上完整法最低RMSD候选自1.50 ?降至1.28 ?，虚拟角MAE自1.53°降至0.15°，消除线性缩放产生的6个>10°角异常点，展示振幅有序帧中几何退化被抑制。
3.5 无目标NOLB对照
12对多样子集、等20结构候选、匹配实际位移下，完整法最小RMSD候选局部距/角误差全对低于NOLB，中位最小RMSD有利，3 ?后方向一致性显著；PPDK在≥3 ?仍NOLB有利，显示配对态质量存蛋白级异质性，但局部几何优势均匀。
3.6 计算成本
单输入态端到端中位时间：线性Cα-ANM 14.6 s，BA约束14.2 s，线性+ARAP 48.0 s，完整47.0 s；峰值内存中位94.7 MB。约束模态求解因子空间维数中位235（约非约束内维1/3）而未显著增本。
讨论部分总结
讨论指出局部几何获益在酶、核酸结合蛋白、伴侣、转运体、折叠切换及MscS中一致出现，强于配对态结果。核心结论不是“几何变好”而是“在保留过渡集体信息下几何变好”：BA约束在对角化前投影重算模态，使方向与新序重排，保留信息向低频集中；ARAP在有限振幅补偿残差旋转漂移。胶原酶N端连接肽α→β转换、PPDK中心域111°摆动等跨多模态事件，单模态分支难全复现，说明局限来自单模态NMA表达力而非约束剔除。方法输出振幅有序可检Cα帧，便于结合独立实验证据优选后重建。与NOLB对照显示在同等信息/候选数下局部几何更优，PPDK例外反映运动分布异质。未来可允许多约束模态联合生成或局部松弛拓展。
结论部分翻译
BA-ARAP-NMA旨在解决直接放大Cα模态导致的局部骨架几何渐进损失。通过在一级局部距离约束下重算低频模态并施以共旋转ARAP，该方法在降低局部扭曲的同时保留了模态空间承载的大部分过渡相关构象信息。在35对实验双态结构中，BA约束与ARAP映射对局部几何贡献分明；匹配实际位移下完整法在每对降低相邻Cα距离与虚拟角误差，提升几何控制可达位移并改善蛋白对级配对态候选集。约束并非简单删除常规Cα-ANM基中的运动，而是改变保留构象信息在低频空间的分配。局部几何效应对双输入结构与12/15/18 ? ANM截断均鲁棒。在等结构数无目标NOLB受控协议下，回溯选出的BA-ARAP-NMA最小RMSD候选在12对比较中所有对局部Cα误差更低，1–10 ?中位集级最小RMSD更低（PPDK仍为配对态RMSD的蛋白级反例）。该框架在集体低频运动主导构象变化时最有效，宜视为可检振幅分辨Cα结构生成器，经独立证据优选后供重建精炼。
