第2章距离分布与观测模型

第一部分 · 第2章

第2章　距离分布、背景与观测模型

由分子构象和自旋标签到时域信号及反演结果

2.1 章节范围与核心结论

PDS实验直接记录的是回波幅度随扫描时间的变化。距离分布通过偶极核进入该时域信号，同时受到标签构象、激发带宽、取向选择、相干衰减、分子间背景、残留偶极路径及噪声影响。最终得到的距离分布依赖观测模型和反演方法，不能直接等同于蛋白骨架或核酸骨架的几何距离分布。

本章建立从结构系综到实验信号的分层模型，明确各层变量及其适用条件。正则化算法和软件实现只介绍与观测模型有关的必要部分，完整的数据处理流程留在后续分析章节。

解释原则

距离分布的峰位、宽度和面积具有不同含义。峰位主要反映被观测自旋中心的典型间距；宽度同时包含分子构象、标签构象和实验分辨因素；峰面积只有在激发、标记、弛豫及路径权重得到合理处理时，才可用于估计构象群体比例。

2.2 从结构系综到反演结果的层级

距离实验可分为五个连续层级。每一层都引入新的变量或权重。

图2.1　距离信息的分层关系。结构与标签产生几何分布，实验选择产生有效分布，脉冲序列与背景把分布映射到时域信号，反演算法产生报告中的距离结果。

早期概述常把这些影响写成多个一维分布的卷积。该写法只能作直观示意。标签构象通常依赖分子构象和局部环境，实验权重也可能依赖距离和取向，因此一般模型应使用条件概率和观测权重。

2.3 分子构象、标签构象与自旋距离

2.3.1 条件几何分布

以 $X$ 表示分子构象，$L_A$ 和 $L_B$ 表示两个标签的构象，$E$ 表示溶剂、膜环境、温度和局部相互作用等条件。两个电子自旋代表点的位置写为 $\mathbf R_A(X,L_A)$ 与 $\mathbf R_B(X,L_B)$。理想的自旋距离分布为

$$P_{\mathrm{spin}}(r)=\int p(X)\,p(L_A,L_B\mid X,E)\,\delta\!\left[r-\left|\mathbf R_A(X,L_A)-\mathbf R_B(X,L_B)\right|\right]dX\,dL_A\,dL_B.$$ (2.1)

式(2.1)允许标签构象与分子构象相关，也允许两个标签之间存在位阻或静电相关。独立rotamer近似对应

$$p(L_A,L_B\mid X,E)\approx p(L_A\mid X,E)p(L_B\mid X,E).$$ (2.2)

式(2.2)便于快速计算，但在短距离、刚性双功能标签、标签之间直接接触及强局部相互作用条件下可能不充分。

2.3.2 不同距离定义

距离名称	几何对象	主要用途	常见偏差来源
电子自旋距离	两个磁性自旋中心或有效自旋密度中心	PDS和CW直接关联的量	自旋密度离域、取向模型
标签中心距离	nitroxide N-O中点、Cu或Gd离子等代表点	标签系综预测	代表点定义、配位结构
连接位点距离	两个被标记残基或核苷酸位点	位点设计与结构比较	侧链、连接臂和局部构象
骨架距离	Cα、磷酸骨架或指定主链原子	结构模型约束	标签到骨架的几何变换
构象群体距离	多个结构状态下上述距离的混合	构象平衡和异质性	状态数、群体权重和交换时间尺度

2.3.3 实验有效分布

激发、弛豫和相干路径会给不同构象赋予不同权重。定义综合观测权重 $W$，实验有效分布可写为

$$P_{\mathrm{eff}}(r)=\frac{1}{Z}\int p(X)p(L_A,L_B\mid X,E)W(X,L_A,L_B;\boldsymbol\eta)\,\delta[r-r_{AB}]\,dX\,dL_A\,dL_B,$$ (2.3)

其中 $\boldsymbol\eta$ 表示磁场位置、微波频率、脉冲形状、温度、序列延迟及检测窗口等实验设置，$Z$ 是归一化因子。取向选择实验中的 $W$ 对分子取向具有明显依赖；不同构象具有不同 $T_m$ 时，长延迟测量还会改变构象权重。

2.4 从距离分布到时域信号

2.4.1 单路径双自旋模型

对于单一调制路径、孤立自旋对和均匀粉末取向，分子内偶极演化函数为

$$D(t)=\int_0^\infty P_{\mathrm{eff}}(r)K(r,t)dr.$$ (2.4)

若泵浦或等效翻转过程只影响比例 $\lambda$ 的伙伴自旋，归一化分子内信号可写为

$$V_{\mathrm{intra}}(t)=1-\lambda+\lambda D(t).$$ (2.5)

常规四脉冲DEER的基础观测模型为

$$V(t)=V_0\,V_{\mathrm{intra}}(t-t_0)\,V_{\mathrm{inter}}(t-t_0)+\varepsilon(t).$$ (2.6)

其中 $V_0$ 是尺度因子，$t_0$ 是实验时间零点，$V_{\mathrm{inter}}$ 是分子间背景，$\varepsilon(t)$ 表示噪声和未建模误差。部分序列还需加入相干包络、未调制路径和额外偶极路径。

2.4.2 离散前向模型

在距离网格 $r_j$ 和时间点 $t_i$ 上，式(2.4)可写为

$$\mathbf D=\mathbf K\mathbf P,\qquad K_{ij}=K(r_j,t_i)\Delta r.$$ (2.7)

矩阵 $\mathbf K$ 的列在有限时间窗内高度相关，使距离反演成为病态问题。时域噪声的微小变化可能对应多种明显不同的 $P(r)$。正则化或参数化模型用于限制解空间，不能补回实验中未记录的信息。

模型成分	物理含义	主要实验来源	误设后的典型影响
$K(r,t)$	特定序列的偶极响应	哈密顿量、相干路径、粉末平均	峰位、峰数或分布宽度失真
$P_{\mathrm{eff}}(r)$	被实验加权的自旋距离分布	结构、标签和选择函数	结构群体解释偏差
$\lambda$	调制路径的有效幅度	翻转效率、标记率、路径概率	群体比例或多自旋判断偏差
$V_{\mathrm{inter}}$	分子间偶极背景	自旋浓度、维度和空间分布	低频及长距离成分偏差
$t_0$	偶极演化的实验零点	脉冲中心、有限脉冲和时序定义	相位扭曲及短距离伪成分
$\varepsilon$	随机噪声和模型残差	接收噪声、漂移、脉冲不稳定	假峰、不确定区间过窄

2.5 分子间偶极背景

2.5.1 稀释自旋簇近似

目标分子内的自旋构成有限尺寸自旋簇。其他分子上的自旋虽然距离较远，数量却很多，会产生平滑的时域衰减。Goldfarb与Stoll第19章在簇内和簇间贡献可分离时写出

$$V(t)=F(t)B(t).$$ (2.8)

对于均匀三维分布和常规DEER条件，背景常表现为单指数衰减。更一般的经验形式为

$$B(t)=\exp\left[-(k|t|)^{d/3}\right],$$ (2.9)

其中 $k$ 与被泵浦自旋浓度、激发概率及偶极常数有关，$d$ 是有效空间维度。均匀三维分布对应 $d=3$。膜、表面、排除体积和非均匀浓度可能产生不同指数。

2.5.2 背景模型的适用范围

背景类型	适用条件	模型特征	主要检查方式
三维均匀DEER背景	体相冻结溶液、自旋随机分布	近似单指数	浓度系列、单标样品、长时间段
低维或非均匀背景	膜、表面、聚集体、排除体积	伸展指数或显式空间模型	改变浓度、样品环境和拟合区间
SIFTER背景	单频率相干转移与过滤并存	调制部分和未调制部分可具有不同背景形状	单自由基、双自由基及SIDRE对照
RIDME背景	伙伴自旋在混合期发生纵向跃迁	依赖混合时间、弛豫及跃迁概率	混合时间系列和单标对照
多路径背景	五／七脉冲DEER及其他多路径序列	各路径具有自己的重聚时间和背景因子	路径模型、时间原点及残差结构

2023年的SIFTER理论研究指出，调制信号所乘的分子间背景与未调制SIFTER贡献具有不同形状。适当处理顺序包括先拟合并减去未调制项，再除以另一背景函数。由此可见，DEER背景公式不能直接作为全部PDS方法的通用背景。

背景与长距离

背景衰减和长距离偶极振荡都集中在低频区域。短时间窗内二者相关性较强。背景形式、拟合区间和时间零点的变化需要作为距离分布稳定性检查的一部分。

2.6 激发、取向和弛豫选择

实验测得的 $P_{\mathrm{eff}}(r)$ 包含序列选择权重。权重可拆成若干相互关联的组成：

$$W=W_{\mathrm{label}}W_{\mathrm{exc}}W_{\mathrm{orient}}W_{\mathrm{relax}}W_{\mathrm{path}}W_{\mathrm{detect}}.$$ (2.10)

权重	来源	可能与结构相关的情形	对照或校正
$W_{\mathrm{label}}$	标记效率和标签完整性	不同构象的位点可及性不同	CW自旋计数、质谱或化学标记对照
$W_{\mathrm{exc}}$	脉冲带宽、频偏和B₁	不同构象对应不同谱位或取向	反转轮廓、nutation和频率扫描
$W_{\mathrm{orient}}$	g、A各向异性及磁场位置	刚性标签与骨架方向相关	多场位测量、全局取向模拟
$W_{\mathrm{relax}}$	$T_1$、$T_m$及瞬时扩散	局部运动和溶剂暴露随构象变化	延迟系列、温度系列及弛豫测量
$W_{\mathrm{path}}$	相干路径和翻转概率	多自旋拓扑与路径概率相关	路径模型和相位循环对照
$W_{\mathrm{detect}}$	回波积分窗口和接收相位	不同谱区相位或线形不同	保存复数瞬态并统一积分定义

柔性nitroxide标签在常规低场选择条件下常接近各向同性权重。刚性标签、Cu(II)、Gd(III)、光生三重态及高场测量更容易出现显著取向选择。取向选择既可提供结构方向信息，也会改变一维距离分布的权重。

2.7 时间零点、有限时间窗与噪声

2.7.1 时间零点

实际核函数使用 $t-t_0$。矩形脉冲常以脉冲中心定义有效翻转时刻；整形脉冲的有效作用时刻还与扫频方向、频偏和动态相位有关。若 $t_0$ 设置错误，时域偶极振荡获得额外相位，频域Pake谱形失去预期对称关系。

2.7.2 有限时间窗

实验信号等于无限时域信号与窗口函数的乘积：

$$V_{\mathrm{obs}}(t)=W_t(t)V_{\mathrm{ideal}}(t).$$ (2.11)

时域截断在频域对应与窗口函数傅里叶变换的卷积。较短的 $t_{\max}$ 首先限制低偶极频率，因此长距离峰容易变宽并与背景相关。过大的时间增量 $\Delta t$ 会产生快速偶极振荡的欠采样，主要影响短距离。

2.7.3 噪声和残差

最简单模型假设 $\varepsilon\sim\mathcal N(0,\sigma^2\mathbf I)$。实际数据还可能含有时间相关噪声、基线漂移、相位误差、周期性核调制和脉冲幅度漂移。随机噪声决定统计区间，结构化残差通常指向模型缺项或实验伪影。

数据表现	优先检查	可能影响
残差呈随机分布且方差稳定	噪声模型及权重	适合统计重采样
残差含规则振荡	核调制、额外偶极路径、时间零点	可能形成假距离峰
末端持续单向偏离	背景形式和基线漂移	长距离及分布宽度偏差
零时间附近尖峰或折点	脉冲重叠、死时间和相位瞬变	短距离伪成分
重复扫描间整体尺度变化	温度、频率、谐振器和样品稳定性	调制深度及群体比例偏差

2.8 多自旋、多路径和高自旋谐波

2.8.1 独立翻转近似下的多自旋信号

一个observer自旋与多个伙伴自旋耦合时，在独立翻转和单一observer路径近似下，分子内信号可写为多个自旋对因子的乘积：

$$F(t\mid\mathcal C)=\prod_{j=1}^{N}\left[1-\lambda_j+\lambda_jK(r_{0j},t)\right],$$ (2.12)

其中 $\mathcal C$ 表示一个具体空间构型。样品信号还需对构型系综平均。乘积中的交叉项会产生组合频率；若使用线性双自旋模型处理，可能出现没有对应实际自旋对的峰。

2.8.2 多路径实验

五脉冲和七脉冲DEER可同时产生多个偶极路径。每条路径具有重聚时间 $T_p$、幅度 $\lambda_p$、频率缩放或谐波因子 $n_p$，并可能带有独立背景。数据模型需要直接描述这些路径。把时域数据强行变换为理想四脉冲DEER曲线会丢失路径信息。

2.8.3 高自旋RIDME

对于 $S>1/2$ 的伙伴自旋，混合期可能发生多个 $\Delta m_S$ 跃迁，产生偶极基频及高次谐波。前向模型可写成若干谐波核的加权和：

$$K_{\mathrm{RIDME}}(r,t)=\sum_{n=1}^{n_{\max}}c_nK(r,nt),\qquad \sum_n c_n=1.$$ (2.13)

系数 $c_n$ 与自旋量子数、温度、混合时间、弛豫网络和实验条件有关。使用四脉冲DEER核处理此类信号会把谐波误换算为较短距离。

2.9 距离反演、模型选择与不确定度

2.9.1 正则化目标函数

非参数距离分布常通过带非负约束的正则化求解：

$$\hat{\mathbf P}=\underset{\mathbf P\ge0}{\operatorname{argmin}}\left\{\|\mathbf V_{\mathrm{exp}}-\mathbf M(\boldsymbol\theta)\mathbf P\|_2^2+\alpha^2\|\mathbf L\mathbf P\|_2^2\right\}.$$ (2.14)

$\mathbf M(\boldsymbol\theta)$包含偶极核、背景、调制深度、时间零点和路径参数，$\alpha$控制数据拟合与平滑约束的平衡，$\mathbf L$定义平滑阶数。DeerLab论文展示了背景参数与非参数 $P(r)$ 同时拟合的单步方法，并指出短时间窗内先拟合背景、再反演距离的两步处理可能不稳定。

2.9.2 参数化与非参数分布

方法	表示方式	适合情形	主要限制
高斯或多高斯模型	有限个峰位、宽度和面积	状态数有结构依据、信号质量较高	峰数和峰形假设影响结果
非参数正则化	距离网格上的非负向量	分布形状缺少先验	平滑强度影响峰宽和细节
贝叶斯模型	参数或函数的后验分布	需要显式先验和模型比较	结果依赖先验及计算设置
结构系综模型	构象及标签的加权组合	多距离对和已有结构模型	可辨识性取决于距离对数量与标签模型

2.9.3 不确定度的分层报告

拟合区间通常只描述给定模型内的噪声传播。背景形式、路径模型、正则化强度和标签模型的差异属于模型相关变化，需要另行比较。

层级	来源	推荐检查	报告形式
重复性	重复扫描和样品重复制备	独立重复测量	峰位、均值或分布差异
噪声传播	有限S/N	bootstrap、协方差或噪声重采样	置信区间或置信带
背景敏感性	背景函数与拟合区间	多种合理背景及时间截断	结果包络和残差比较
反演敏感性	$\alpha$、$\mathbf L$、峰形和峰数	正则化扫描、参数模型比较	稳定峰及不稳定细节
实验模型	取向、多路径、谐波、交换和有限脉冲	替代核与独立实验	条件化结论
结构模型	标签rotamer和分子构象	多标签位点、多个结构模型	结构群体范围

2.10 距离分布的结构解释与报告规则

2.10.1 峰位、宽度和面积

分布特征	可直接陈述的内容	需要附加条件的解释
峰位或众数	有效自旋距离的主要位置	映射到骨架距离需要标签模型
均值或中位数	指定距离区间和归一化条件下的统计量	长尾和截断会改变数值
峰宽	观测分布的宽度	分子柔性、标签柔性和实验分辨不能仅凭一条分布分离
多峰	存在多个可拟合距离成分	离散构象状态需要独立结构信息或多距离对验证
峰面积	当前模型下的相对分布权重	群体比例需要校正标记、激发、弛豫和路径权重

2.10.2 最小报告内容

距离结果应同时保存原始复数时域数据、处理后的信号、拟合曲线、残差、背景模型、时间零点、调制深度、距离网格、反演设置和不确定区间。应用结论同时注明自旋距离与结构距离的区别。

结论强度

单一距离对通常适合判断结构模型是否相容或检测明显构象变化。多个距离对、标签模型和独立结构方法可以约束构象群体。平滑分布中的细小峰不宜单独作为新增构象状态的依据。

2.11 应用案例：双组氨酸Cu(II)-NTA标签的距离建模

研究对象与科学问题：双组氨酸Cu(II)-NTA标签可用于蛋白质结构和骨架构象研究。实验距离直接对应Cu(II)中心，结构解释需要描述标签在α螺旋和β折叠位点上的构象分布。

方法选择依据：Cu(II)中心的位置比长柔性连接臂nitroxide标签更接近蛋白骨架，双功能配位还限制标签构象。该标签适合检验rotamer library能否把实验Cu-Cu距离分布转换为结构模型约束。

实验与建模概述：Hasanbasri等构建了基于分子动力学和半经验构象采样的两组rotamer libraries，并通过chiLife双功能标签模型预测GB1、GSTA1和YopO等体系的Cu-Cu距离分布，再与实验分布比较。

主要结果：论文报告测试数据中预测与实验分布众数的绝对偏差为0.0–1.2 Å，平均值为0.6 Å。YopO案例还比较了直线形和弯曲α螺旋结构对距离分布的影响。

结果所说明的问题：该研究表明，标签构象模型可把电子自旋距离与具体骨架模型建立定量联系。双功能Cu(II)标签较受限的构象空间有利于区分相近结构模型。

适用范围与限制：所报误差来自论文使用的测试体系和特定rotamer libraries。配位几何、局部二级结构、溶剂及蛋白构象超出建库条件时，需要重新检查标签分布。实验距离与预测结果的一致性也受PDS背景、取向选择和反演设置影响。

2.12 本章结论

结论	适用条件	报告要求
结构与标签共同决定自旋距离分布	标签分布可由条件概率描述	区分电子自旋、标签中心和骨架距离
观测信号可分为分子内项和背景项	稀释自旋簇和可因子化条件	注明背景类型及拟合区间
标准Pake核适合均匀粉末取向	取向选择及带宽效应较弱	金属中心和刚性标签需检查选择函数
有限时间窗限制长距离信息	可用演化时间由相干衰减和S/N决定	进行时间截断和背景敏感性检查
多路径、高自旋和多自旋需要专用核	额外路径或谐波不可忽略	保存路径幅度、重聚时间和谐波系数
拟合区间只描述指定模型内的不确定度	模型结构固定	另行比较背景、反演和标签模型

返回第一部分目录

模型成分	物理含义	主要实验来源	误设后的典型影响
\(K(r,t)\)	特定序列的偶极响应	哈密顿量、相干路径、粉末平均	峰位、峰数或分布宽度失真
\(P_{\mathrm{eff}}(r)\)	被实验加权的自旋距离分布	结构、标签和选择函数	结构群体解释偏差
\(\lambda\)	调制路径的有效幅度	翻转效率、标记率、路径概率	群体比例或多自旋判断偏差
\(V_{\mathrm{inter}}\)	分子间偶极背景	自旋浓度、维度和空间分布	低频及长距离成分偏差
\(t_0\)	偶极演化的实验零点	脉冲中心、有限脉冲和时序定义	相位扭曲及短距离伪成分
\(\varepsilon\)	随机噪声和模型残差	接收噪声、漂移、脉冲不稳定	假峰、不确定区间过窄

权重	来源	可能与结构相关的情形	对照或校正
\(W_{\mathrm{label}}\)	标记效率和标签完整性	不同构象的位点可及性不同	CW自旋计数、质谱或化学标记对照
\(W_{\mathrm{exc}}\)	脉冲带宽、频偏和B₁	不同构象对应不同谱位或取向	反转轮廓、nutation和频率扫描
\(W_{\mathrm{orient}}\)	g、A各向异性及磁场位置	刚性标签与骨架方向相关	多场位测量、全局取向模拟
\(W_{\mathrm{relax}}\)	\(T_1\)、\(T_m\)及瞬时扩散	局部运动和溶剂暴露随构象变化	延迟系列、温度系列及弛豫测量
\(W_{\mathrm{path}}\)	相干路径和翻转概率	多自旋拓扑与路径概率相关	路径模型和相位循环对照
\(W_{\mathrm{detect}}\)	回波积分窗口和接收相位	不同谱区相位或线形不同	保存复数瞬态并统一积分定义

第2章 距离分布与观测模型