第6章 · 6.8

6.8 变温、EXSY 与线形分析

Variable Temperature, EXSY, and Lineshape Analysis
本页范围与学习目标
变温谱、交换谱和线形拟合提供不同的动态证据。本节说明温度对平衡、位移和速率的共同影响,推导一个可核验的双状态 EXSY 增长模型,区分交换与 NOE,再将单向速率换算为指定方向的活化自由能。实验条件、模型检验与不确定度共同限定动力学结论。
先读图解和核心解释;需要细节时再展开。
先看图 · 阅读入口

变温、交叉峰与线形共同约束动态

变温改变多个输入,动态报告要同时说明占比、速率和模型。

变温、交叉峰与线形共同约束动态 步骤 1 改变温度 速率和平衡都可能变化 分别保存状态占比 步骤 2 记录时间转移 EXSY 随混合时间变化 检查交换与弛豫 步骤 3 联合拟合 线形与状态模型对账 再计算相应自由能 变温、交叉峰与线形共同约束动态 步骤 1 改变温度 速率和平衡都可能变化 分别保存状态占比 步骤 2 记录时间转移 EXSY 随混合时间变化 检查交换与弛豫 步骤 3 联合拟合 线形与状态模型对账 再计算相应自由能
图 6.8-1 本页概念关系图。箭头表示本页采用的检查顺序。

读图时保留这项条件:活化自由能描述越过势垒,平衡自由能描述状态相对稳定性。

6.8.1 变温实验同时改变速率与平衡

前节把交换速率与静态频率差联系起来。本节进一步讨论怎样从一组实验取得动力学证据。变温核磁共振(variable-temperature NMR, VT NMR)在控制样品条件的前提下改变温度,观察位移、占比和线形。升温常使某些构象互变加快,但位移、弛豫和氢键也会随温度变化;一个峰向另一个峰移动,不能单凭移动方向归因于交换。

展开说明与推导

开始测量前,先确认样品在所选温度范围内保持溶解、稳定且密闭条件适宜。温控气流的设定值与样品实际温度可能有差别,动力学分析应使用经校准的样品温度,并记录校准方法及平衡等待。每个温度点检查不参与目标交换的信号、溶剂和基线,尽量保持可比的激发与处理条件。温度改变后重新匀场可能必要,不能把全谱普遍变宽作为局部构象交换的证据。实验操作应遵循所在实验室的仪器流程。

对可逆的两状态 A、B,保留前节 (k_{AB})、(k_{BA}) 与 (k_{\mathrm{ex}}=k_{AB}+k_{BA}) 的约定。在平衡且没有其他出口时,(p_Ak_{AB}=p_Bk_{BA})。若两状态各有可独立积分且单位核响应可比的报告信号,慢交换谱可以估计状态占比。对于同一分子内的两个构象,这个比例描述分子处在哪个构象;它不能解释成分子中两类原子的数量发生变化。若报告区域重叠、弛豫响应不同或小峰低于检测能力,比例的误差须另行估计。

\[K=\frac{p_B}{p_A}=\frac{k_{AB}}{k_{BA}},\qquad \Delta G_{B-A}=-RT\ln K.\]
(6.8-1)

式中 (K) 是此两状态模型的无量纲平衡常数,(R) 是摩尔气体常数,(T) 用 K,\(\Delta G_{B-A}\) 用 \(\mathrm{J\,mol^{-1}}\)。这里是同一分子的状态平衡,若涉及缔合或不同分子数的反应,还须写清标准态与浓度条件。平衡差描述 B 相对 A 的稳定程度;越过过渡态的能垒则决定转换速度。两个状态可以等占比而互换很慢,也可以占比差很大却互换很快。

定义 · 三种不同信息

位移差说明两环境在当前条件下的频率区别,占比说明平衡分布,速率说明环境转换的动力学。变温时三者都可能改变;拟合应允许有物理依据的温度依赖,不能固定一切参数后把全部变化交给交换速率。

展开说明与推导

可以把一组变温观察整理为检验过程:在低温能分开 A、B 时,对每个报告核记录两列位移与比例;温度升高时跟踪原子标签,避免按峰的左右位置重新命名。若两个区域的比例同步改变,可能支持同一个全分子平衡;若只有一个区域异常,应检查重叠与局部过程。最后回到起始温度,比较完整谱与原谱是否恢复。不能恢复时,溶剂损失、沉淀或分解都可能参与,整组数据不能继续按单纯可逆构象模型解释。

6.8.2 EXSY 跟踪混合时间内的环境转换

交换谱(exchange spectroscopy, EXSY)把一个核在混合时间之前和之后的环境频率关联起来。若核先在 A 中编码频率,在混合时间内转到 B,随后以 B 的频率被检测,就可能出现 A—B 交叉峰。适用的通常是两个状态的共振仍可分辨、而混合时间内又有足够转换的范围。没有交叉峰可能因为交换过慢、混合时间不适合、小状态过弱或信号重叠,不能直接证明两状态不交换。Jeener 等(1979)的原始研究建立了这种二维交换分析思路。

展开说明与推导

常见质子 EXSY 与 NOESY 使用相同类型的脉冲实验,名称强调其分析任务。交叉峰既可能来自同一核的环境交换,也可能来自不同核间的偶极交叉弛豫。判断前先完成原子指认:属于两个构象的同一个氢,应与整套构象信号的对应关系相容;空间 NOE 则连接不同核的位置。必要时结合选择性一维实验、ROESY、温度变化和多组报告核,检验这些解释。峰的相位可作辅助依据,但相位处理、分子运动范围和传递伪影会改变观察,不能只凭一个颜色宣布交换。

同一核由环境 A 转到 B,可在 A 编码与 B 检测处形成交叉峰演化:记录 A 频率核的初始环境混合时间 tmA → B 的转换检测:记录 B 频率同一核的新环境二维位置(示意,非实测)A—A、B—B:仍在原环境的对角峰A—B、B—A:需检验来源的交叉峰ABAB
图 6.8-2 EXSY 的环境关联示意。两轴频率标签须结合实际处理方向阅读;图中交叉点只表示可能的交换关联,并不提供距离或速率的实测数值。
展开说明与推导

一次 EXSY 实验给出一个混合时间的传递结果。若目标是求速率,应采集多个混合时间,包括接近零的参考,比较交叉峰增长和对角峰衰减,检查是否符合一个共同的动力学模型。长混合时间可提高某些交换峰的可见性,同时加大弛豫损失和间接传递。定量时使用经过基线和重叠检查的峰体积,重复测量评估误差。仅有一张高对比度等高线图,尚不足以反演两个方向的速率。

6.8.3 从传递概率理解交叉峰增长

先用可逐步验算的简化模型。假定两状态等占比、单向速率均为 \(k\),纵向弛豫速率相同为 \(R_1\),没有 NOE、其他交换路径或额外磁化源。跟踪初始标记在 A 的磁化,环境概率满足一级方程。总概率为一,B 的概率从零开始,以 \(2k\) 的时间常数趋向二分之一:

展开说明与推导
\[P_{A\to B}(t_m)=\frac{1-e^{-2kt_m}}2,\qquad P_{A\to A}(t_m)=\frac{1+e^{-2kt_m}}2.\]
(6.8-2)

本节把 \(I_{AB}\) 定义为先在 A 编码、后在 B 检测的峰体积,\(I_{AA}\) 为对应对角峰。二者在所设模型中带有相同的初始占比、尺度因子与 \(e^{-R_1t_m}\) 衰减,取比值时这些共同因子消去,得到:

\[r=\frac{I_{AB}}{I_{AA}}=\tanh(kt_m),\qquad k=\frac{1}{2t_m}\ln\frac{1+r}{1-r}.\]
(6.8-3)

这里 \(t_m\) 用 s,\(0\le r<1\),返回的 \(k\) 是单向速率,仍须乘二才是总交换速率。此关系来自本节假设下的概率模型,并非任意 EXSY 谱的通用公式。不同 \(R_1\)、不等占比、NOE、非选择性激发或第三状态会破坏共同因子的消去;真实数据应使用带弛豫和交换项的磁化传递矩阵,联合拟合对角峰与交叉峰。

教学设定:在 \(t_m=0.100\ \mathrm{s}\) 时测得已经正确归属的体积比 \(r=0.197375\)。先代入比值式,\(\ln[(1+r)/(1-r)]\) 约为 0.400,得到 \(k=2.00\ \mathrm{s^{-1}}\),\(k_{\mathrm{ex}}=4.00\ \mathrm{s^{-1}}\)。再把同一速率用于 0.050 s,预测比值约 0.099668;用于 0.200 s,预测约 0.379949。第二、第三个时间点提供模型检验,如果偏差远超积分误差,不能只保留第一个点的结果。

方法 · 用多时间点检验假设

先建立两套状态信号的对应表,再检查零混合时间参考、交叉峰增长、对角峰衰减和独立报告核。拟合残差随混合时间系统变化时,应考虑弛豫、间接传递、第三状态或积分问题。不同机制需要不同控制实验,不能通过增加拟合参数自动证明机制。

习题 6.8-1 · 单向速率的反演

采用本节简化模型,混合时间为 0.200 s,A—B 与 A—A 的体积比为 0.379949。求单向和总交换速率。若另一报告氢给出明显不同的比值,应如何处理?

参考解答

\(k=\ln[(1+0.379949)/(1-0.379949)]/(2\times0.200)\approx2.00\ \mathrm{s^{-1}}\),因此总交换速率约为 \(4.00\ \mathrm{s^{-1}}\)。同一全分子两状态转换应能由共同速率描述;不一致时先查峰体积、重叠、归属和 NOE 等额外传递,再检查两状态及相同弛豫假设。不能平均两个不相容的数值后声称模型通过检验。

展开说明与推导

短混合时间下,交叉传递概率近似为 \(kt_m\),这使增长初始斜率容易与单向速率关联;时间很长时,比值趋向一,速率差异在比值上变得很小。即使实验峰仍很强,这一饱和区域也未必能精确确定速率。选择时间点时要兼顾初始增长、足够的交叉峰信噪比与弛豫损失,不能只重复最显眼的长时间谱。误差估计也应重新考虑这个非线性变换:比值接近一时,小的积分误差可被对数反演明显放大。

6.8.4 用整段线形约束交换模型

线形分析(lineshape analysis)比较实验谱与包含交换、位移、弛豫及必要耦合的预测谱。前节的等占比双单峰模型适合作为起点;真实谱若有多重峰,就必须保留耦合及自旋系统的影响。观测峰顶的间隔在中间交换时已经偏离静态频率差,不能把该间隔重新输入合并公式,再将输出当作独立测量。

展开说明与推导

可以从交换很慢的低温谱估计静态位移和耦合,但低温值未必等于合并温度下的静态值。应利用有依据的温度趋势或联合数据约束它们。处理窗函数也会改变线宽,模拟比较需使用相同处理作用,或明确把它计入模型。每个温度点可允许合理的振幅、相位和基线参数,避免仪器变化被迫解释为动力学;同时不能给每一点任意曲线,使交换参数失去约束。

参数相关是这类分析的实际困难:未知静态频率差、占比和内禀线宽,可能与速率相互补偿。应检查不同初值是否收敛到同一解释,比较多个报告信号,报告参数不确定度和残差。若数据只包含一个很窄的快交换峰,往往无法同时确定所有参数;EXSY、独立弛豫测量或更宽温度范围可以增加约束。软件返回一个数值只说明优化程序结束,尚未检验参数是否被数据确定。

6.8.5 活化自由能与平衡自由能分开计算

在过渡态理论及合适的一级过程描述下,Eyring 方程(Eyring equation)把速率与活化 Gibbs 能(Gibbs energy of activation)联系起来。IUPAC 的活化自由能与过渡态理论词条给出了相应定义。本节采用传递系数 \(\kappa=1\) 的常规近似,写成:

展开说明与推导
\[k_{AB}=\frac{k_BT}{h}\exp\!\left(-\frac{\Delta G^{\ddagger}_{A\to B}}{RT}\right),\qquad \Delta G^{\ddagger}_{A\to B}=RT\ln\!\left(\frac{k_BT}{hk_{AB}}\right).\]
(6.8-4)

\(k_B\) 是 Boltzmann 常数,\(h\) 是 Planck 常数,\(R\) 是气体常数;右侧对数内为无量纲比值。注意这里用的是对应方向的速率,并非 \(k_{\mathrm{ex}}\)。若把总速率代入对称两状态的单向式,所得能垒会降低 \(RT\ln2\)。采用其他传递系数、状态简并或不同标准态时,应说明能垒定义随模型怎样变化。

平衡比较两个最低点,动力学比较各最低点到过渡态的高度状态 A状态 B过渡态从 A 起算的能垒从 B 起算的能垒反应或构象坐标(示意)G
图 6.8-3 两状态自由能曲线示意,未按实测能量绘制。两最低点高度差决定平衡偏向,两个方向的活化能垒从各自状态起算;状态占比相等不会使能垒为零。
展开说明与推导

例如,教学设定 \(T=298.15\ \mathrm{K}\)、\(k_{AB}=2.00\ \mathrm{s^{-1}}\),代入常数得到 \(\Delta G^{\ddagger}_{A\to B}\approx71.3\ \mathrm{kJ\,mol^{-1}}\)。这是一点温度下、指定模型的活化自由能,不能据此分别确定活化焓与活化熵。只有在适当温区内取得多个可靠速率,并有理由近似把活化焓和熵视为常数时,才可用 \(\ln(k/T)\) 对 \(1/T\) 的关系分别估计它们;拟合斜率的不确定度通常比单点自由能更敏感。

习题 6.8-2 · 比例与速度

在 298.15 K,A:B 为 4:1,测得 \(k_{AB}=2\ \mathrm{s^{-1}}\)。求 \(k_{BA}\)、总交换速率和 B 相对 A 的自由能差。采用上述过渡态近似,哪个方向的能垒更高?

参考解答

\(p_A=0.8,p_B=0.2\),通量平衡给 \(0.8\times2=0.2k_{BA}\),所以 \(k_{BA}=8\ \mathrm{s^{-1}}\),总速率为 \(10\ \mathrm{s^{-1}}\)。\(K=1/4\),\(\Delta G_{B-A}=RT\ln4\approx3.44\ \mathrm{kJ\,mol^{-1}}\),B 更高。A→B 较慢且从更低的状态起算,其能垒更高;两方向能垒相差约 3.44 kJ/mol,不能将此差当作任一方向的整个能垒。

6.8.6 从动态证据到可检查的报告

一份交换分析记录应保留状态定义、原子对应、温度、介质、场强、混合时间、处理参数和拟合假设。变温结果说明线形如何变化;EXSY 说明哪些环境在实验窗口内关联;线形拟合在模型条件下估计速率;能垒换算再依赖速率方向与理论近似。这些证据可以互相检查,但不应合并成没有条件的“发生某种运动”。

警示 · 交叉峰和合并均有替代解释

NOE、间接传递或重叠可制造看似交换的关联;温度引起的位移靠近也可降低分辨度。可靠结论需检验整套信号对应及可逆性,并比较其他实验。无法区分时保留机制候选,将下一步实验指向候选产生不同预测的位置。

习题 6.8-3 · 选择补充证据

两个信号间出现 NOESY 交叉峰,升温后二者合为一峰。能否直接用等占比合并公式报告转动能垒?列出尚需检查的条件。

参考解答

不能。先确认两信号是否代表同一核的两个环境,检查整套状态对应、相位、NOE 或间接传递;测温度序列和降温恢复,检查占比、静态频率差、内禀宽度与耦合。满足等占比、窄线双单峰近似后,合并速率也只是近似估计;完整线形和多混合时间 EXSY 可进一步检验模型。最后使用样品实际温度与单向速率换算指定方向的活化自由能,报告假设及不确定度。

关键术语

变温 NMR (VT NMR)
控制并记录温度,比较谱学参数与线形的实验方法。
状态占比 (state population)
样品处于某一指定状态的分数。
平衡常数 (equilibrium constant)
按指定状态与标准态定义的平衡比值。
交换谱 (EXSY)
关联同一核在混合时间前后所处环境的谱学实验。
混合时间 (mixing time)
环境交换或磁化传递发生并被分析的一段实验时间。
线形分析 (lineshape analysis)
用带明确交换及弛豫假设的模型比较完整谱形。
活化 Gibbs 能 (Gibbs energy of activation)
指定方向中初始状态到过渡态的自由能差。
Eyring 方程 (Eyring equation)
过渡态理论框架下将速率与活化自由能关联的关系。

参考文献与延伸阅读

  1. McConnell HM. 1958. Reaction Rates by Nuclear Magnetic Resonance. The Journal of Chemical Physics, 28: 430–431. 原始交换方程。
  2. Jeener J, Meier BH, Bachmann P, Ernst RR. 1979. Investigation of exchange processes by two-dimensional NMR spectroscopy. The Journal of Chemical Physics, 71: 4546. 二维交换原始研究。
  3. Claridge TDW. 2016. High-Resolution NMR Techniques in Organic Chemistry, 3rd ed. Elsevier. 动态效应与 EXSY 教材。
  4. Hebrew University of Jerusalem, NMR Laboratory. NOESY and EXSY. 实验目的与传递解释。
  5. IUPAC. Gibbs energy of activation, Gold Book, G02631. 活化自由能定义。
  6. IUPAC. Transition state theory, Gold Book, T06470. 过渡态理论与速率关系。