通过先减后除法实现纵向顺磁性弛豫增强的精确测量

《Magnetic Resonance Letters》：Accurate measurement of longitudinal paramagnetic relaxation enhancement through division after subtraction

　　

顺磁性弛豫增强（Paramagnetic Relaxation Enhancement，PRE）可为生物大分子的结构与动态表征提供关于未成对电子的宝贵距离约束。虽然横向 PRE（Γ2）已被广泛应用，但其定量解释常受化学交换及其他因素干扰。相比之下，纵向 PRE（Γ1）尽管量级远小于 Γ2，却是电子?核距离的更忠实报告者。在这项工作中，研究人员引入了先减后除法（Division After Subtraction，DAS）策略，利用同步采样方案在分析上抵消共模伪影。研究人员利用位点定向共价标记与溶剂 PRE 系统验证了该方法，并证明其与标准分别拟合流程相比提高了准确性。对于易溶残基，传统非线性回归因溶剂介导效应而失效，而 DAS 尤为有效。DAS 在 Γ1 测量中提升的精度使得对生物分子构象景观的更严谨探究成为可能，从而为基于 Γ2 的 PRE 方法学提供了有力补充。

研究背景方面，顺磁性弛豫增强（Paramagnetic Relaxation Enhancement，PRE）源于未成对电子与核（通常为质子）之间的偶极?偶极相互作用，其强度比质子间核奥弗豪泽效应（Nuclear Overhauser Effect，NOE）高出600倍以上，可提供约40 ?范围内的距离限制，并适于阐明生物大分子的构象涨落与低占有态信息。常用顺磁探针包括氮氧自旋自由基、Mn²⁺或 Gd³⁺螯合复合物等，后者因金属离子位于探针中心、亲水性更强且可同时测量横向（Γ₂）与纵向（Γ₁）PRE而具优势。横向 Γ₂测量虽有单点法、两点法及优化脉冲序列等发展，但易受化学交换、同核标量耦合及伪接触位移（Pseudo?Contact Shift，PCS）影响；纵向 Γ₁虽量级较小，却对化学交换不敏感，是更忠实的电子?核距离报告者，尤其适用于Ni²⁺、Fe³⁺等快电子弛豫金属离子体系。然而，Γ₁定量需拟合纵向衰减曲线，对参数要求更高，且传统分别拟合（Separate Fitting，SEP）会因顺磁与抗磁样品平衡磁化强度差异、饱和不完全、NOE（核奥弗豪泽效应）及溶剂交换等引入系统误差，在柔性及溶剂可及残基处尤为严重，甚至得出非物理负值，因此亟需更稳健的测量方法。

研究人员开展了基于饱和恢复（Saturation Recovery）的纵向弛豫测量新方法研究，提出先减后除法（Division After Subtraction，DAS），通过顺磁与抗磁样品同步采样、时域差分与归一化相除，在数学上消去平衡磁化强度 I_eq、初始残余磁化 I₁′ 等未知常数及共模的 NOE、溶剂交换贡献，并推导了 Γ₁ 解析表达式与不确定度公式，给出 τ≈0.5 T₁、Δt 不使信号衰减低于首差强度的30%等参数选取经验准则，及 Δt=τ 的等距三点加速采样策略。研究人员利用位点定向共价标记的 Rpn10 UIM1?泛素（Ubiquitin，Ub）相互作用体系与溶剂 PRE（solvent PRE，sPRE）的 Ub＋Gd?DTPA?BMA 体系，在 Bruker 600 MHz 与 700 MHz 谱仪上采集饱和恢复数据，分别用传统 SEP 非线性回归与 DAS 处理，并结合修正所罗门方程（Solomon Equations）的数值模拟，比较二者在 NOE、溶剂介导及柔性残基下的精度与误差分布。

主要关键技术方法方面，研究人员采用两类样本：共价标记样本为与人 Ub 瞬态相互作用的人 Rpn10 UIM1 域（E232C 突变并共轭 EDTA?Mn²⁺），溶剂 PRE 样本为在 2 mM Gd?DTPA?BMA 存在下的 ¹⁵N 标记 Ub；利用饱和恢复脉冲序列在 298 K 下采集多延迟（如 0.1–1.7 s 八点、0.01–1.6 s 六点）的 ¹H?¹⁵N HSQC 数据；传统分析对抗磁与顺磁强度分别做单指数拟合得 R_1,dia 与 R_1,para，相减得 Γ₁ 并以蒙特卡洛迭代估不确定度；DAS 分析取同步时点的强度差 ΔI，计算归一化差之比解析得出 Γ₁ 并按误差传播公式算不确定度；数值模拟在 Python 中用矩阵形式所罗门方程与矩阵指数模拟目标自旋与环境自旋的交叉弛豫（σ）对表观衰减的影响，对比两种处理在不同 σ、溶剂饱和/保留状态下的误差热图；所有谱图用 NMRPipe 处理，余弦窗函数与零填充后提取峰强。

研究结果部分，在 2.1 Derivatization of the DAS measurement approach 中，研究人员指出饱和恢复优于反转恢复（Inversion Recovery），但因脉冲不完善、场不均匀及饱和不完全，实际衰减含未知初值与时序偏移。传统分别拟合需从两条曲线各提取 R₁ 再相减，受限于有限重复延迟 D₁ 导致的平衡态差异与残余磁化而引入系统误。DAS 通过对同步采样点定义间隔 τ 的差分信号 ΔI_t=I_t+τ?I_t，再在 Δt 后取 ΔI_t+Δt 并归一化，顺磁与抗磁分别得含 e^{?(R₁+Γ₁)τ 与 e?R₁τ 的表达式，相除后 R₁ 相消，留下仅含 Γ₁ 的解析关系 Γ₁=?(1/Δt)ln[(ΔI_para,t+Δt/ΔI_para,t)/(ΔI_dia,t+Δt/ΔI_dia,t)]，从而在数学上免疫初始饱和缺陷与样品间平衡磁化差异；不确定度由原始强度方差与偏导直接算出；τ 过短则 ΔI 小、误差放大，过长则信号趋近 I_eq 也降信噪比，实践取 τ≈0.5 T₁，Δt 不使次窗信号低于首窗 30%；采用 Δt=τ 的等距三点采样（如 0.1、0.6、1.1 s）可复用重叠数据点，将采集时间减半且保持伪影消除。}

在 2.2 NOE contributions are minimized with the DAS scheme 中，研究人员用引入哑变量 E=1 的修正所罗门方程将非齐次边界转为齐次矩阵形式，模拟目标自旋 H_z^G 与环境自旋 H_z^SS 的纵向弛豫与交叉弛豫 σ。结果显示，随负 NOE（σ 负）增大，分别拟合的 Γ₁ 误差可超 60%，而 DAS 的低误差区域显著扩大，因内部 NOE 被当作共模因子在归一化相除中消去；溶剂介导的 NOE 与交换亦类似，若在饱和恢复前用选择性水翻转回脉冲（water flip?back）保留溶剂磁化，可进一步减小溶剂对 H_z^G 恢复的影响与 Γ₁ 误差，DAS 结合保留溶剂磁化在强弱 NOE 下均维持高精度，而分别拟合在溶剂饱和时误差剧增。

在 2.3 Experimental demonstration of improved measurement accuracy with the DAS scheme 中，共价标记体系下，DAS 所得 Γ₁ 在多数结构残基上系统性高于分别拟合，差异微小（刚性核平均基线偏移 0.08±0.06 s^?1）源于 I₀<0 的不完全饱和，分别拟合敏感于此基线误而 DAS 相减消去；差异最大的残基 T9（位于 Ub β1?β2 间溶剂可及环）与 G76（高度柔性 C 端）在分别拟合中决定系数（Coefficient of Determination，R²，COD）明显下降，偏离单指数行为：T9 的抗磁曲线初斜率被单指数过估、I_eq 低估，导致 R_1,dia 虚高、相减后 Γ₁ 严重低估甚至负值；G76 的顺磁初段加速恢复使 R_1,para 过估，相减后 Γ₁ 虚高；机制在于溶剂水交叉弛豫 σ 与局部 R₁ 或 R₁+Γ₁ 相近时，多指数成分与单指数回归相构破坏或相长干涉，而 DAS 同步相减将非指数交叉弛豫作共模解析消去。sPRE 体系（Ub＋Gd?DTPA?BMA）中同样在 β1?β2 环与 C 端残基见最大分歧，分别拟合的 R² 跌落指明溶剂 NOE 扰动致多指数偏离，DAS 通过同步相减消去多指数畸变，给出更准确一致的 Γ₁，反映未被局域水交换动力学污染的 Outer?sphere 顺磁增强，适于定量表面作图。

总结讨论部分，研究人员指出横向 Γ₂ 常被化学交换卷积，而纵向 Γ₁ 较少受相干退相干伪影影响；且 Γ₁ 量级更小，快交换极限条件 k_ex>>ΔΓ 更易满足，对两态构象交换 A?B 而言 ΔΓ₁=|Γ_1,A?Γ_1,B| 远小于 ΔΓ₂=|Γ_2,A?Γ_2,B|（Γ₂ 可达数千 s^?1），故 Γ₁ 在表征复杂/异质构象景观具显著优势。DAS 概念上类比已建立的 Γ₂ 强度比法，通过归一化强度差消去影响抗磁与顺磁样品的共模系统误，对需长时间采集、谱仪漂移或脉冲不稳的 Γ₁ 测量尤为关键；方法可直接嵌入标准 NMR 脉冲序列，现有等距时间点的 ¹H 弛豫数据集也可用 DAS 重分析以修正此前因拟合不确定而误判的残基。DAS 提升的精度与实验效率将促进 PRE 实验应用于更具挑战的生物体系，包括多时间尺度构象交换过程映射、内禀无序蛋白（Intrinsically Disordered Proteins，IDPs）与内禀无序区（Intrinsically Disordered Regions，IDRs）表征，其中细微结构偏好常被高柔性与溶剂交换掩盖。该研究发表于《Magnetic Resonance Letters》。

翻译研究结论部分如下：与常受化学交换效应贡献卷积的横向顺磁性弛豫增强（Γ₂）不同，纵向分量（Γ₁）较少易受相干退相干伪影影响。此外，由于 Γ₁ 值量级上固有地小于 Γ₂，快交换极限条件即 k_ex>>ΔΓ 更易达到。具体而言，对于两态构象交换模型（A?B），定义 ΔΓ₁=|Γ_1,A?Γ_1,B| 与 ΔΓ₂=|Γ_2,A?Γ_2,B|，代表交换态间本征 PRE 速率之差。横向 PRE 值（Γ₂）通常比纵向 PRE 值（Γ₁）大一数量级，靠近顺磁中心的核 Γ₂ 易达数千 s^–1。因此，对 Γ₂ 而言，要实现真实快交换区并使观测 PRE 速率等于本征态群体加权平均值，需要极快动力学接近微秒（μs）或更快速率。就此而言，Γ₁ 测量在表征复杂或异质构象景观方面具显著优势。本研究提出的 DAS 方案在概念上类比已建立的 Γ₂ 测量强度比法。通过归一化信号强度差，DAS 方案有效消去影响抗磁与顺磁样品的共模系统误差。此误差抑制能力对 Γ₁ 测量至关重要，因其常需较长实验采集时间，谱仪漂移或脉冲不稳定性可能损害数据完整性。本方法学可轻易整合入标准 NMR 脉冲序列而无需显著修改。值得注意的是，现有具等距时间点的 ¹H 弛豫数据集可用 DAS 方案重新分析，以恢复此前因拟合不确定而被判错的残基的更准确测量。DAS 方法提升的精度与实验效率将促进 PRE 实验更广泛地应用于具挑战性的生物体系，这包括但不限于跨多时间尺度映射构象交换过程，以及表征内禀无序蛋白（IDPs）与区域（IDRs），其中细微结构偏好被高柔性与溶剂交换所掩盖。