蛋白质主链中曲率支撑场的拉普拉斯结构揭示的不稳定性驱动的扭转符号翻转

《The Journal of Physical Chemistry B》：Instability-Driven Torsion Sign Flips Revealed by the Laplacian Structure of Curvature Support Fields in Protein Backbones

　　

扭转符号翻转是蛋白质主链动力学中的基本离散事件，但其几何触发因素仍知之甚少。虽然传统上将这些转变与奇异临界点相关联，但研究人员表明，这些转变主要由局部几何不稳定性驱动，其量化指标为曲率派生支撑场的拉普拉斯量ΔS。通过对来自3000个PDB结构的超过130万个残基位置的大规模分析，研究人员证明了一条稳健的物理规律：翻转事件在拉普拉斯为零（局部平衡）的区域显著减少，并随着不稳定性幅度|ΔS|的增加而呈现逐渐升高的富集。这支持了一种与不稳定性相关的统计解释，即在所分析的静态结构数据集中，几何失衡与扭转反转倾向增加相关。研究人员进一步识别出一种“拓扑锁定”效应，其中环状约束系统地抑制翻转，特别是在高不稳定性区域。多变量建模证实，ΔS提供了标准局部描述符（曲率和扭转）无法获取的预测信息。研究结果建立了一个将离散微分几何与随机主链动力学联系起来的物理可解释框架，揭示了局部几何失衡与全局结构约束如何共同关联于实验确定的静态蛋白质结构中的扭转符号翻转倾向。

蛋白质主链在维持整体结构完整性的同时表现出显著构象柔性，理解局部几何转变如何发生是分子生物物理学的核心问题。扭转符号翻转是主链动力学中最基本的离散事件，但长期以来仅被归因于临界点或奇异性，缺乏对其结构背景和全局约束作用的定量解释。本文在《The Journal of Physical Chemistry B》发表，提出一种以曲率支撑场（curvature support field）的离散拉普拉斯量ΔS为核心的几何框架，旨在检验扭转符号翻转是否由局部几何不稳定性驱动。

研究人员从RCSB蛋白质数据银行（PDB）获取2,926个实验解析结构，涵盖1,335,936个残基位置，识别出570,478个翻转事件。基于Cα坐标计算单位切线、离散曲率κ、扭转τ，定义翻转；构造曲率支撑场S为相邻曲率的局部加权平均（邻域为i±1，均匀权重），并计算ΔS=S_i+1+S_i-1?2S_i。不稳定性幅度定义为|ΔS|。环约束定义为序列间隔大于10个残基且Cα距离小于6.0 ?的残基对。

**扭转符号翻转在拉普拉斯零区域耗竭**
研究人员发现，在|ΔS|<0.01的局部平衡区域，翻转条件概率为0.0298，显著低于基线0.0448，富集因子0.666，Fisher精确检验p<10^?300，优势比0.655。这表明翻转并非优先发生在临界点，而是在拉普拉斯零区域被抑制。

**不稳定性幅度驱动转变**
通过设置不同阈值δ（0.05、0.10、0.20、0.30），条件翻转概率从0.453升至0.587，富集因子从1.059升至1.374，优势比从1.144升至1.954。分箱分析也显示P(flip)随|ΔS|增大而连续上升，最高箱超过0.6，证明几何不稳定性幅度与翻转倾向存在系统性关联。

**多变量分析与局部描述符比较**
以(S,ΔS)为特征建立的逻辑回归模型AUC为0.558；在多变量模型P(flip)=σ(aκ+bτ+cΔS+d)中，ΔS的系数显著，且加入ΔS可稳定提升预测性能。说明ΔS提供了曲率κ和扭转τ所不包含的二阶空间失衡信息。

**拓扑约束抑制不稳定性相关转变**
环约束区域的翻转概率为0.4149，显著低于非环区域的0.4387（p<10^?15），绝对差异约2.4%。在高不稳定性条件下抑制增强：δ=0.10时，环区域为0.448，非环为0.498，抑制比约0.90。这表明环状拓扑约束通过限制可及构象空间，削弱了不稳定性驱动的翻转倾向。

**实值几何不稳定性的充分性**
将特征扩展为复值z=S+iΔS后，模型AUC未提升（≈0.51），性能低于实值模型（AUC=0.558），校准曲线亦显示实值模型更优。说明翻转行为可由实值几何失衡充分描述，无需引入相位型自由度。

**观察总结**
基于超过1.3×10⁶个残基位置的大规模分析，研究确认翻转概率在拉普拉斯零区域耗竭，并随|ΔS|增加而富集；环约束在统计上抑制翻转，复值扩展无增益。

讨论部分强调，该研究基于静态结构，所揭示的是统计几何偏好而非动态过程本身。结果支持不稳定性积累而非几何分岔的机制：ΔS≈0表示曲率支撑均匀分布、净驱动力为零，而大|ΔS|反映曲率诱导的有效应力，可降低翻转的等效势垒。环约束可视为对驱动力的重正化，即|ΔS|→α_loop|ΔS|（0<α_loop<1）。复值检验的负结果进一步排除了相位自由度，支持标量驱动场的图像。

总之，本研究建立了联系离散微分几何与随机主链动力学的可解释框架，证明扭转符号翻转为不稳定性偏置的随机转变，并受全局拓扑约束调节；该框架可用于识别高柔性区域，对配体结合、别构调控和构象转换等功能过程具有潜在意义。