第13章 距离反演与不确定度

第三部分 · 第13章

第13章　距离反演与不确定度

从复数时间迹到距离分布、构象比例和结构限制

13.1 数据对象与不可逆操作

数据层级符号允许操作需要保存
原始复数回波\(V_I(t)+iV_Q(t)\)相位校正、积分窗检查I/Q、每个scan及元数据
归一化时间迹\(V(t)/V(0)\)零点、异常点和噪声估计归一化因子及删除点
form factor\(F(t)\)背景模型联合或分步处理背景参数和残差
距离分布\(P(r)\)正则化、参数模型或全局拟合网格、约束和不确定区间
结构解释构象、比例、距离约束标签模型与多状态比较结构版本和标签库
不可逆操作

只保存实部、平滑后时间迹或背景扣除后的form factor会丢失诊断信息。正式报告应保留原始复数数据及完整参数文件。

13.2 复数相位、零点和噪声

常数相位可通过最小化虚部能量或联合拟合复数模型确定。若虚部呈时间相关结构，应优先检查接收相位、相位瞬变、脉冲重叠和频率漂移，不能用单一相位角强制消除。

步骤计算对象检查
相位校正\(V(t)e^{-i\phi}\)虚部残差应接近零均值噪声
零点确定\(t\rightarrow t-t_0\)正负时间或多路径重聚位置
异常点处理脉冲重叠区、放大器恢复区删除规则预先记录
噪声估计虚部、重复scan或后段残差异方差和相关噪声
归一化零点附近模型值或回波参考不得用受伪影影响的单点

13.3 统一前向模型

$$V(t)=A\,B(t;\boldsymbol\beta)\prod_{p=1}^{N_p}\left[1-\lambda_p+\lambda_p\int K(t_p,r;\boldsymbol\eta_p)P(r)\,dr\right]+\epsilon(t).$$(13.1)

\(A\)为幅度，\(B(t)\)为分子间背景，\(p\)表示偶极路径，\(\lambda_p\)和\(t_p\)分别为路径幅度与重聚时间，\(\boldsymbol\eta_p\)包含取向选择、谐波和有限脉冲参数。单路径、各向同性、\(S=1/2\)自旋对只是式(13.1)的简化情况。

$$K(t,r)=\int_0^1\cos\!\left[2\pi\nu_{dd,\perp}(r)(1-3u^2)t\right]du.$$(13.2)

式(13.2)对应均匀粉末取向和点偶极近似。刚性标签、金属中心、高场或明显谱位选择需要引入取向权重。

13.4 分子间背景模型

$$B(t)=\exp\left[-(kt)^{d/3}\right].$$(13.3)

三维均匀稀释体系常取 \(d=3\)。膜、二维限制、排除体积、聚集和纵向弛豫可能改变有效指数。背景与长距离分量在有限时间窗内高度相关，固定一个看似合理的维数会低估模型不确定度。

背景策略采用条件检查指标风险
固定三维指数均匀稀释体相样品浓度序列与后段残差把结构长距并入背景
自由维数伸展指数环境维度未知参数轮廓及区间维数与距离分布相关
基于浓度的物理背景有效pump浓度已知拟合速率与浓度比例标记率和排除体积误差
多数据集全局背景浓度或场位序列共享参数后的联合残差错误共享造成系统偏差

教材建议在条件允许时使 \(t_{max}>2t_{dec}\)，更充分的背景段可取 \(t_{max}>4t_{dec}\)。低调制深度、宽分布或有限\(T_m\)样品经常达不到该条件，此时应报告背景模型敏感性。

13.5 非参数反演与Tikhonov正则化

$$\hat{\mathbf P}=\arg\min_{\mathbf P\ge0}\left\{\|\mathbf W(\mathbf V-\mathbf M\mathbf P)\|_2^2+\alpha^2\|\mathbf L\mathbf P\|_2^2\right\}.$$(13.4)

\(\mathbf M\)是离散偶极核与背景模型，\(\mathbf W\)是噪声权重，\(\mathbf L\)定义平滑阶数，\(\alpha\)控制数据拟合与平滑约束的权衡。\(\alpha\)过小产生噪声驱动的窄峰，\(\alpha\)过大使真实多峰合并并加宽。

选择方法评价对象报告内容
L-curve残差范数与正则项范数曲线、拐点和相邻\(\alpha\)结果
交叉验证留出数据预测误差分割方式和误差曲线
AIC／BIC参数化模型复杂度候选模型和相对差值
残差白噪声检查时间相关结构残差图、自相关或频谱

13.6 参数化距离模型

高斯、偏斜高斯、worm-like chain或有限状态混合模型适合具有明确结构假设的样品。参数化模型给出峰位、宽度和群体比例的直接区间，也会把模型形状假设写入结果。

模型适用结构描述需检查
单高斯单一主构象和近似对称标签分布偏斜残差与肩峰
高斯混合少数离散构象状态峰数辨识和比例相关性
偏斜分布受限rotamer或非对称柔性链与背景参数的相关性
聚合物链模型无序链段或连接臂轮廓长度与持久长度范围

13.7 不确定度的三个来源

来源实例评价方法输出
统计噪声回波噪声和有限平均bootstrap、重复测量或似然区间点态或参数置信区间
参数相关背景速率、调制深度、零点和峰宽profile likelihood、后验或二维轮廓参数联合区间
模型选择背景维数、路径数、谐波系数和峰数候选模型重复拟合模型间结果范围

距离分布的点态区间不等于峰位置和峰面积的区间。构象比例应由积分或参数模型获得，并把背景和路径参数的不确定性一同传播。

13.8 多路径、多自旋与高自旋谐波

情况必要模型误用单路径的表现
5p／7p DEER多个重聚时间及路径幅度复制峰、峰宽异常和振荡相位错误
多标记或寡聚体自旋拓扑、标记概率和多距离组合调制深度过高、ghost peaks
Gd(III) RIDME基频及高次谐波权重距离按谐波比例缩放
取向选择g/A/ZFS张量与空间方向联合权重不同场位出现不一致的一维分布

高自旋RIDME的overtone coefficients应由标尺、混合时间序列或明确的物理模型确定。把系数作为任意自由参数会增加距离分布与谐波权重的不可辨识性。

13.9 取向选择与全局分析

取向选择数据应联合拟合多个场位或频率位置。每个数据集具有独立幅度、背景和选择函数，距离及结构方向参数可按物理关系共享。只把多场位距离分布取平均会丢失方向信息，也可能掩盖带宽失配。

输入共享参数数据集特有参数
多场位DEER距离分布、相对取向模型场位选择、调制深度、背景
二维SIFTER偶极张量和结构几何直接维相位、谱位响应
金属—自由基RIDME距离和标签几何场位相关\(T_1\)、核调制和背景

13.10 电子自旋距离到结构距离

$$P_{calc}(r)=\sum_{i,j}w_i^{(A)}w_j^{(B)}\,\delta\!\left(r-|\mathbf r_i^{(A)}-\mathbf r_j^{(B)}|\right).$$(13.5)

\(i,j\)表示两个位点的标签rotamer，\(w\)是相应权重。实验峰包含分子构象、标签构象、选择权重和分析分辨函数。电子自旋中心距离不能直接替换为Cα、磷酸骨架或金属配位位点距离。

报告对象最低说明
实验距离分布前向模型、背景、正则化及区间
标签预测分布结构文件、标签库、位点和温度假设
构象比例状态模型、共享参数及置信区间
骨架限制电子中心到骨架的映射方式

13.11 结果报告的最低项目

类别最低项目
原始数据复数时间迹、时间单位、零点定义、删除点和噪声估计
实验参数序列、场位、频率、脉冲、延迟、温度、SRT、shots和scans
前向模型背景、路径、核、谐波、取向及调制深度定义
反演距离网格、正则化、非负约束和模型比较
不确定度统计区间、参数相关和模型敏感性
结构解释标签模型、结构版本、状态比例及解释限制

13.12 本章结论

距离反演应围绕完整前向模型展开。背景、路径、谐波和取向选择与距离分布同时影响时间迹；仅改变正则化参数不能修正错误的物理模型。