核磁共振(NMR)波谱法能够以残基特异性的方式表征蛋白质的结构动力学。本文提供了一套实用的实验方案,用于采集对皮秒到纳秒时间尺度敏感的NMR 15N R1 和 R2 弛豫以及{1H}-15N异核Overhauser效应(hetNOE)实验数据。
核磁共振(NMR)波谱法能够以残基特异性的方式表征蛋白质的结构动力学。本文提供了一套实用的实验方案,用于采集对皮秒到纳秒时间尺度敏感的NMR 15N R1 和 R2 弛豫以及{1H}-15N异核Overhauser效应(hetNOE)实验数据。
核磁共振(NMR)波谱技术可用于在生理温度下研究溶液中的蛋白质。通常,蛋白质主链的酰胺基团或侧链中的甲基基团被用作报告蛋白质结构动力学的探针。对于经15N标记且完全质子化的球状蛋白质样品,其主链结构动力学研究通常适用于分子量不超过50 kDa的蛋白质。当结合侧链氘代与横向弛豫优化波谱(TROSY)技术时,该上限可扩展至球状蛋白质200 kDa,若研究重点为侧链,则可进一步扩展至1 MDa。当研究对象为内在无序蛋白质(IDPs)或含有内在无序区域(IDRs)的蛋白质时,这些分子量限制不再适用,研究范围可远超上述数值。其原因是IDPs或IDRs具有高度的内部灵活性,通常在动力学上相互解耦。多种NMR方法可在从皮秒到小时的广泛时间尺度上,提供原子分辨率的蛋白质结构动力学信息。标准的15N弛豫测量可概述蛋白质的内部灵活性,并表征蛋白质主链在快速的皮秒至纳秒时间尺度上的动力学行为。本文提供了一套实用的操作指南,用于设置和记录NMR的15N R1、R2以及异核Overhauser效应(hetNOE)实验。我们展示了典型数据,并解释了如何在进行更复杂分析之前对其进行简单的定性解读。
蛋白质的功能不仅由其三维结构决定,还取决于其结构动力学,包括蛋白质内部的柔韧性以及在不同构象之间发生的结构转变。核磁共振(NMR)波谱技术可用于研究溶液中蛋白质的结构动力学1,2,3。近年来质子检测的固态NMR技术的发展,也使得在溶解度较低的状态下(例如脂双层膜)对蛋白质动力学进行表征成为可能4,5,6。在溶液NMR中,可以研究蛋白质主链和侧链的结构动力学。对于球状蛋白质,只要蛋白质经过15N同位素标记,便可对其主链结构动力学进行研究,分子量上限可达50 kDa。若结合侧链氘代和横向弛豫优化谱(TROSY)技术,该上限可扩展至200 kDa7,8。当研究重点为侧链动力学时,可研究的蛋白质及复合物的分子量范围可进一步扩展至1 MDa2,9。
所述的分子量限制不适用于内在无序蛋白(IDPs),这些蛋白通常表现出较高的内在动态性。真核生物蛋白质组中超过30%由IDPs或内在无序区域(IDRs)组成10,11,12,13。它们在许多细胞过程中发挥核心作用,例如信号转导和转录1,并经常参与细胞内相分离14,15,16,17。IDPs在生理条件下缺乏明确的三维(3D)天然结构,其能量景观呈弱漏斗状或崎岖不平17,18。由于IDPs或IDRs主链上疏水性较低且存在强烈的静电排斥,导致其缺乏向刚性结构折叠的驱动力19。IDPs在与其他结合伴侣形成复合物时,常采用折叠构象10,20,21。此外,翻译后修饰(PTMs)可扩展IDPs或IDRs的折叠可能性22,23。IDPs的错误折叠已被确认为多种疾病的原因,包括神经退行性疾病15,24,25,26。
固有无序蛋白(IDPs)和固有无序区域(IDRs)具有高度的内部灵活性21,27,28。通过分子动力学模拟以及来自实验数据的约束条件,获得了展示原子位置和二面角变化的构象集合29,30,31,32。由于动态性以及在冷冻状态下导致的无序性,弥散的电子密度使得利用结构生物学中的前沿方法(如冷冻电镜或X射线晶体学)对其进行结构表征变得困难。此外,用于低温实验的结晶条件或样品制备技术可能影响IDPs所经历的构象空间。然而,溶液核磁共振(NMR)适用于高度动态的蛋白质,因此非常适合研究IDPs16,20,22,28,29,30,31,32,33,34,35,36,37,38。
如上所述,溶液核磁共振(NMR)提供了多种技术,可用于研究蛋白质内部动力学在宽时间尺度上的变化(图1),这些技术主要基于自旋弛豫31,33,38,39,40,41,42。
蛋白质主链酰胺基团中15N原子核的自旋弛豫由蛋白质内部动力学和集体运动(在相关情况下包括旋转扩散)引起的1H-15N键角取向变化所诱导27,43,44,45,46,47,48,49,50,51。在短于旋转相关时间τR(分子旋转一个弧度所需的时间,也称为整体翻滚相关时间)的时间尺度上,化学位移各向异性(CSA)和偶极耦合(D)具有活性,且不会被蛋白质的旋转扩散所平均。蛋白质主链的内部动力学,包括键角的变化、键的重新取向以及旋转翻滚,会引起CSA和偶极耦合张量的随机波动,导致局部磁场的变化,最终引发核磁共振(NMR)自旋弛豫47,48,52,53。这些波动可通过一个总体相关函数来描述。总体相关函数的傅里叶变换称为谱密度函数。在半经典的Redfield弛豫理论中,NMR弛豫速率常数可由这些谱密度函数的线性组合来描述54。
20世纪90年代初发展的主链15N核磁共振(NMR)弛豫实验包括15N R1、R1ρ以及{1H}-15N核Overhauser效应实验,这些实验对快于蛋白质旋转相关时间τR的皮秒(ps)至纳秒(ns)时间尺度的运动具有敏感性45,55,56,57。为了表征慢于旋转相关时间τR的主链动力学,通常采用所谓的弛豫弥散实验,如R1ρ和Carr-Purcell-Meiboom-Gill(CPMG)实验,这些实验对微秒(µs)至毫秒(ms)时间尺度的动力学过程敏感44,46,58,59,60,61。慢于微秒的动力学过程可通过15N化学交换饱和转移(CEST)核磁共振62、交换谱学(EXSY,毫秒至秒级)或实时(RT)核磁共振(秒至小时级)进行捕捉63,64。顺磁探针的顺磁弛豫增强(PRE)效应以及残余偶极耦合(RDCs)可用于评估从皮秒到毫秒整个时间范围内的动力学行为65,66,67,68。

图1:蛋白质主链动力学的时间尺度及不同核磁共振动力学实验的敏感时间窗口。核磁共振(NMR)提供了多种方法,可用于表征跨越广泛时间尺度的蛋白质主链动力学。图中在相应的时间尺度上标示了蛋白质主链所经历的不同类型的运动。蛋白质的旋转相关时间τR是指蛋白质完成一次整体旋转(一个弧度)所需的时间。比蛋白质旋转相关时间τR更快的运动可归因于蛋白质内部的柔性。各种核磁共振实验及其对相应时间尺度的敏感性在箭头下方标出。请点击此处查看该图的放大版本。
以下方案描述了Lakomek等人69和Stief等人70建立的核磁共振(NMR)弛豫实验,采用灵敏度增强的异核单量子相干(HSQC)检测方案。在进入实验操作之前,将简要概述NMR自旋弛豫及NMR弛豫实验的基本原理。由于篇幅限制并为保证本方案的可理解性,该概述只能保持简略(因而并不完整)。
纵向弛豫或自旋-晶格弛豫,由T1 时间或 R1 = 1/T1 速率常数,描述磁化矢量恢复至玻尔兹曼平衡的过程。在平衡状态下,磁化矢量沿外加磁场方向排列,该方向定义为实验室坐标系的 z 轴。高频下的谱密度1小(H)15N) 拉莫尔频率(核磁共振共振频率,例如600 MHz用于¹H) 1H 对于14.1特斯拉磁体)以及这些拉莫尔频率的线性组合对 15N R1 舒张,其特征为 15N R1 以 rad·s 为单位测得的速率常数-1. 运动的时间尺度与这些拉莫尔频率成反比;因此,皮秒到纳秒时间尺度上的运动对弛豫速率常数R有贡献1对于表现出整体翻滚运动且可定义旋转相关时间的分子,其 R1 (T1曲线在 ωτ 处显示一个极大值(极小值)R = 1,旋转相关时间τR 以及所考虑自旋的拉莫尔频率 ω。如果多个拉莫尔频率均有贡献,则频率最低者占主导地位,例如 ωN 在……的情况下 15N R1快速运动区(ωτR 远小于1)适用于在低磁场和低黏度条件下快速翻滚的小分子。慢运动 regime(ωτR 远大于1)适用于翻滚较慢的大分子、高磁场和高黏度条件。
球状折叠蛋白在溶液中表现出整体翻滚运动,可分配一个旋转相关时间。然而,对于内在无序蛋白而言,整体翻滚的概念不再适用,且通常无法用单一的旋转相关时间来描述。此时,残基特异性的内禀相关时间变得更为重要。
所描述的脉冲序列用于测量15N R1弛豫速率(图2),其基于一种灵敏度增强的HSQC检测实验,并采用Echo/Anti-Echo检测方式实现四象限检测69,70,71。该序列使用强度和时长可变的短梯度,用于相干选择并增强水信号抑制70。在此期间,15N纵向磁化将发生弛豫。较长的弛豫时间将导致该准3D谱中相应2D平面的信号强度降低(延迟数据点记录在第三维中)。为实现更长的弛豫时间,下文所述的循环单元将被重复执行的次数逐步增加。由于在弛豫延迟期间,15N化学位移各向异性(CSA)与1H和15N偶极耦合(D)之间的交叉弛豫效应同时存在,因此需要在中心位置施加一个选择性作用于酰胺质子的I-BURP-2 180°脉冲72,以重新聚焦交叉弛豫的贡献(若不重新聚焦,将导致15N R1速率常数出现偏差和错误)。

图2:测定NMR弛豫速率常数的NMR脉冲序列示意图。(A) 15N R1ρ,(B) 15N R1,以及(C) 采用灵敏度增强型HSQC检测方案的hetNOE实验69,70。90°(x)脉冲用窄矩形表示,180°(x)脉冲用宽矩形表示,除非另有说明。应用以下相位循环:φ6 = y, y, -y, -y;φ7 = y, -y;φrec = y, -y, -y, y。通过反转梯度G5的极性和φ7的相位循环实现四象限检测(回波/反回波检测)。(A) 15N R1ρ 实验:黑色矩形代表自旋锁定,其持续时间变化以获得不同的弛豫延迟。自旋锁定前后三角形表示绝热形状脉冲,用于将磁化矢量沿有效磁场轴Beff对齐。G10是可选梯度,用于在演化阶段防止水磁化产生的辐射阻尼。(B) 15N R1 实验:括号部分表示序列的循环单元,重复n次以匹配所需的弛豫延迟。(C) hetNOE脉冲序列与R1和R1ρ脉冲序列的后半部分相似,即t1演化时间和HSQC检测单元。然而,15N磁化是直接激发的,无需INEPT过程。通过施加至少5秒的一系列180(1H)脉冲实现质子磁化的饱和(以实现1H与15N之间的交叉弛豫)。在参考实验中,则施加相同长度(此处为5秒)的无脉冲空闲延迟。G5是可选梯度,用于防止辐射阻尼;梯度G4极性的反转结合相位φ7 = y, -y, -y, y实现四象限检测。以积算符表示的磁化转移步骤以红色标出。请点击此处查看此图的放大版本。
弛豫速率常数 R2 描述了横向极化(位于与外加磁场垂直的 xy 平面)由于自旋之间相位相干性的丧失而导致的弛豫过程,从而引起可检测磁化信号的衰减53,54。在高频和低频处的谱密度函数对 R2 有贡献,这一点与 R1 类似。然而,R2 的最大贡献来源于零频处的谱密度。因此,R2 对旋转翻滚运动非常敏感,该运动由旋转相关时间 τR, 描述,在室温下对于一个小的球状蛋白而言,其数量级为几纳秒(ns)。因此,主链在数百皮秒(ps)到低纳秒(ns)范围内的较慢运动贡献最大。15N 原子核化学位移张量各向同性部分的调制所引起的主链交换动力学,原则上会向 R2 弛豫速率常数添加一个交换贡献项 R2ex43,44,49,60,61。在所述实验中,通过施加自旋锁定来抑制 R2ex 的贡献,该自旋锁定可重聚焦比自旋锁定圆频率倒数更慢的动力学过程。自旋锁定是一种长时间的连续波射频脉冲,其作用是使磁化矢量保持沿有效磁场轴 Beff 方向排列(该方向为自旋锁定场 ω1 与 15N 载波频率的化学位移偏移量的矢量和(见下文))。沿 B1,eff 轴排列的磁化弛豫称为 R1ρ 弛豫,其包含 R1 分量和 R2 分量。利用公式 (1),可由 R1ρ 和 R144,73 计算出 R2:
(1)。
有效磁场B轴之间的夹角eff 以及外加磁场 B0 是
. ω1 自旋锁定的射频振幅,Ω 为自旋锁定与化学位移之间的偏移量 15相应残基的N化学位移以及 15N载波频率 44,73.
R1ρ脉冲方案(图2A,70)与15N R1方案非常相似,区别仅在于弛豫延迟部分。为了测量15N R1ρ弛豫速率,必须在磁化矢量通过一个与自旋锁定具有相同射频(RF)幅度的绝热脉冲被调整至有效场轴Beff方向后,启动自旋锁定。通过改变自旋锁定的持续时间,可以获得不同的弛豫延迟。
稳态 {1H}-15N核核奥弗豪泽效应(1H-15N NOE),以下称为hetNOE,是交叉弛豫速率与 15N纵向弛豫速率。这会导致稳态极化的降低 15由于质子极化饱和时与质子发生交叉弛豫而导致的N45,53,54,74,75交叉弛豫依赖于拉莫尔频率之和与差的谱密度函数 1H 和 15N 的拉莫尔频率。因此,hetNOE 对快速的皮秒级动力学(< 100 ps)以及皮秒-纳秒级动力学(由于其R1 依赖性)。该序列69 (图2C基于灵敏度增强的HSQC检测,采用梯度回波/反梯度回波(Echo/Anti-Echo)梯度实现四相检测。为实现质子磁化饱和及由此产生的异核NOE,需先将平衡态的质子磁化反转,随后通过快速施加约5倍于该弛豫时间的180°脉冲进行饱和。 15 N T1在对照实验中,恢复延迟时间等于饱和延迟时间,但不包含饱和过程 1H 180°脉冲序列。在参考实验及另一实验中额外加入 D1 = 2 s 的延迟。 1氢饱和。两项实验连续进行,仅在是否施加 1H 180°脉冲(饱和)或不施加(参考)。实验中记录的谱峰强度比值 1H饱和度除以参照实验(无180°质子脉冲序列)的信号强度,得到{1H}-15N NOE(hetNOE)值。
以下方案描述了Lakomek等人69和Stief等人70建立的核磁共振(NMR)弛豫实验。我们重点介绍采用灵敏度增强型HSQC检测方案的NMR脉冲序列。15N R1和R1ρ实验的实施方法详见Stief等人70的描述,而异核NOE(hetNOE)实验则由Lakomek等人69描述。
1. 核磁共振样品制备
注意:对蛋白质进行同位素标记是为了进行高维核磁共振(NMR)和高级NMR实验。当在富含培养基(例如Luria-Bertani [LB] 或 2x酵母提取物胰蛋白胨培养基 [2YT])中实现大肠杆菌(E. coli)表达并完成蛋白质纯化,且每升培养液获得数毫克蛋白质后,制备同位素标记的NMR样品通常相对简单。
2. 核磁共振弛豫实验上机前的准备工作
注意:所述的核磁共振弛豫实验专用于布鲁克(Bruker)波谱仪。这些实验已在配备低温和室温1H、15N和13C三共振探头的Avance III与Avance Neo控制台(使用Bruker软件Topsin 3.6或更高版本操作)上进行过测试。
3. 实施核磁共振弛豫实验
注意:NMR 弛豫脉冲序列(图 2)可在以下网址获取: https://www.ipb.hhu.de/en/teams/team-lakomek/pulsesequences ,或在扩展的生物核磁共振数据库(BMRB)存储库(bmrbig102)中获取。
4. 记录的核磁共振实验的处理与分析
注意:谱图使用 Bruker 系统记录,数据处理在 Unix 或 Linux 操作系统上进行。谱图处理与数据分析使用 NMRPipe80 和 python3 完成。NMRPipe 软件可从 https://www.ibbr.umd.edu/nmrpipe/index.html 下载。 基于 NMRPipe 的处理脚本可从以下网站获取:https://www.ipb.hhu.de/en/teams/team-lakomek/pulsesequences 或扩展的生物核磁共振数据库(BMRB)存储库(bmrbig102)。 建议使用 NMRPipe。若无法使用或不希望使用 NMRPipe,也可选用其他替代软件,如 CCPN81 或 SPARKY(SPARKY 3 及其后续版本 NMRFAM-SPARKY82 或 POKY83)。
以下展示的是在囊泡SNARE蛋白突触小泡相关膜蛋白2(Synaptobrevin-2,1-96位氨基酸,常称为VAMP2,即囊泡相关蛋白2)上记录的一些典型的核磁共振(NMR)弛豫数据。NMR数据采集所用样品为171 µM的15N标记Synaptobrevin-2(1-96)(下文简称为Syb-2),溶于含150 mM NaCl、0.1 mM TCEP和1 mM EDTA的50 mM MES缓冲液(pH 6.0)中。所有实验数据均在278.15 K条件下采集,样品体积为250 µL,装入3 mm NMR样品管。实验在配备Bruker 5 mm QCI 1H、15N、13C、31P四共振低温探头的Bruker 600 MHz AVANCE III HD波谱仪上进行。
图3 展示了根据核磁共振(NMR)弛豫实验获得的NMR弛豫曲线拟合得到的残基特异性15N R1和R2弛豫速率常数。纵轴表示沿蛋白质序列(横轴)各个残基对应的弛豫速率常数。此外,还显示了从1H-15N异核NMR测量中饱和与非饱和(参考)实验的强度比值提取得到的1H-15N异核NOE值(hetNOE)。
异核NOE数据和15N R2弛豫速率常数比15N R1数据更容易解释,因此我们将从它们开始。
异核NOE:
异核{1H}-15N NOE对快速皮秒时间尺度(<100 ps,因其依赖于交叉弛豫)的内部动力学敏感,同时也对皮秒至纳秒(ps - ns)的动力学过程敏感(因其依赖于R1),对于折叠蛋白,其值通常介于0与(略小于)1之间。在较高磁场下,该值会系统性升高(但始终低于1)。
它们描述了蛋白质内部的柔性。当比较相同磁场强度下的hetNOE值时,规律如下:{1H}-15N NOE值越小,相应残基的动态性越高。对于内在无序蛋白(IDPs),通常观察到较小的残基特异性{1H}-15N NOE值,例如在600 MHz下约为0.5或更小。然而,对于IDPs,甚至可观察到负的{1H}-15N NOE值,最低可达500 MHz下的-2。对于单体Syb-2,在600 MHz下测得的{1H}-15N NOE值从N端的-0.5到C端连接结构域的0.5不等。其数值从N端到C端逐渐增加,表明柔性逐渐降低/刚性逐渐增强(图3C)。
15N R2弛豫速率常数
15N R2弛豫速率常数对调节N-H键矢量取向的皮秒至纳秒尺度动力学过程敏感,这些过程包括蛋白质整体翻滚运动、肽平面运动以及主链二面角的扭转。与{1H}-15N核奥弗豪泽效应(NOE)值和R1不同,15N R2弛豫速率常数主要由与蛋白质整体(旋转)翻滚相关的J(0)谱密度所主导。因此,对于球状折叠蛋白,15N R2弛豫速率常数的大小会根据蛋白质的整体翻滚速率和分子尺寸进行缩放。而对于真正的内在无序蛋白(IDPs),无法定义全局的旋转相关时间。尽管如此,仍可拟合出如聚合物体系所预期的、依赖于链长的慢速相关时间84,但残基特异性的内部分子运动对15N R2弛豫速率常数的影响更为显著。在同一蛋白质中,15N R2弛豫速率常数可用于区分具有不同残基特异性内部运动性的结构区域与非结构区域。观测到的15N R2弛豫速率常数越小,相应残基的流动性越高。较高的15N R2弛豫速率常数是二级结构元件(如α-螺旋或β-折叠)的特征,而较低的R2弛豫速率常数则典型地出现在环状区域和内在无序区域。对于单体形式的Syb-2 (1-96),我们观察到整体较低的15N R2弛豫速率常数,这与其内在无序的特性一致。从N端开始至C端连接结构域,内部刚性逐渐增强(流动性降低),随后在柔性C末端区域再次减弱(图3B)。
15N R1弛豫速率常数
R1弛豫速率对皮秒(ps)及更快的纳秒(ns)尺度的分子内运动具有敏感性。然而,与之前的方法相比,其对内部动力学的直接解释更为复杂。当将T1(即R1的倒数)作为分子旋转相关时间τR的函数作图时,会出现一个T1的极小值(对应R1的极大值),该极小值出现在τR, min = 1/ω处。对于快速翻滚的分子,若其τR < τR, min,例如小分子或小型球状蛋白,较低的R1速率常数表明蛋白主链动力学增强(与R2的情况相同)。而对于较大的蛋白质,若其τR > τR, min,情况则相反,较高的R1弛豫速率常数代表更强的动力学行为(与R2相反)。

图 3:核磁共振弛豫数据。 数据按照上述描述的实验方案与工作流程进行采集与分析。(A)Syb2(1-96)蛋白主链上各残基特异性的15N R1弛豫速率常数。(B)15N R2弛豫速率常数,通过公式(1)由 R1ρ 和 R1 数据计算得到。(C)Syb2(1-96)主链上的异核NOE转移速率。请点击此处查看该图的放大版本。
补充图1:导入NMR参数。 显示宏 rpar 的界面。(A) 所有参数集的概览。红色框标出搜索栏。红色数字对应实验方案中的相应步骤。双击以选择一个参数集。(B) 可导入当前实验的参数概览。执行 getprosol 或保留当前实验中特定参数的选项以红色框标出。请点击此处下载该文件。
补充图 2:为实现最佳水信号抑制而定位 1H 载波频率。 显示的是宏 gs 的界面。在实验循环运行期间,可以更改实验的不同参数。此处可使用标记在采集过程中调节所选参数。单击标记上方可使参数值增加“灵敏度”下方指定的数值。灵敏度开关可用于调整标记的步进大小。在右侧可见实时自由感应衰减(FID)信号。校准 O1 的目标是最小化用红色框标出的 FID 区域,因为 FID 通常以水信号为主。FID 区域的减小意味着水信号降低,水抑制效果得到改善。红色数字对应于方案中的相应步骤。请点击此处下载该文件。
补充图3:对 1H-15N HSQC 脉冲序列及其采集参数 (A脉冲序列参数 1H-15N HSQC 实验 HSQCFPF3GPPHWG。点击“E”以在 B 中打开脉冲程序编辑器窗口。B脉冲程序编辑器 1H-15N HSQC实验 HSQCFPF3GPPHWG。带注释的行(第91行)包含一个 15N 90°脉冲显示原始谱线 1H-15N HSQC实验。在线92中, 15N 90°脉冲被更改为一个 15N 180°脉冲。这一改变使得 1H-15N HSQC实验:90°脉冲的校准脉冲程序 15N 硬脉冲校准。如需编辑脉冲序列,请先将序列另存为新名称。使用以下工作流程:点击 文件,另存为,输入新名称,单击 好的,然后单击 设置 PULPROG 添加到数据集。红色数字对应于方案中的相应步骤。 请点击此处下载此文件。
补充图4:设置梯度。 显示了用于更改梯度形状、长度和强度的 gpnam 宏界面。要打开所使用的脉冲序列参数界面,必须输入所有梯度形状 GPNAM。红色数字对应于方案中的相应步骤。请点击此处下载此文件。
补充图5:设置弛豫延迟。 界面 vplist 宏观(与 vclist 宏相同)。(A)选择 vplist 的界面。点击三个点后,界面 B 将打开。B)通过双击选择 vplist,或通过点击创建新列表 文件,以及 新. (C) 界面通过创建新列表或编辑现有列表打开。可通过点击保存列表 文件 和 保存红色数字对应于方案中的相应步骤。 请点击此处下载此文件。
补充图6:15N R1ρ实验的脉冲序列参数。 (A-D)协议中必须修改的所有关键参数均用红色框标出,并附有相应步骤。 请点击此处下载该文件。
补充图7:shapetool宏的界面。 实验中使用的I-BURP2被作为形状脉冲打开。(A)所选形状脉冲的参数。点击用红框标记的Start NMR Simulation按钮后,将弹出模拟参数窗口,其中包含模拟参数:脉冲时长[µs](= 2000 µs)和旋转角度[°](= 180°)。点击Start NMR-SIM按钮后,打开(B)模拟窗口,显示I-BURP2脉冲的反相带宽。红色数字对应于实验方案中的相应步骤。请点击此处下载该文件。
补充图8:calcpowlev宏的界面。(A)用于输入待计算功率脉冲的脉冲长度。(B)参考脉冲(通常为90°硬脉冲)的输入窗口。(C)结果窗口显示新脉冲长度与参考脉冲长度之间脉冲强度的dB差值。注意:仅当两个脉冲长度对应相同形状和旋转角度的脉冲时,该差值才有效。红色数字对应于实验方案中的相应步骤。请点击此处下载该文件。
补充图9:核磁共振谱图处理:NMRPipe转换工具界面。点击读取参数可自动读取实验参数。右侧显示转换脚本。所有参数均可单独设置。使用清除脚本更改参数后,建议更新脚本。点击保存脚本可将脚本保存为fid.com文件,如右侧所示。红色数字对应实验方案中的相应步骤。请点击此处下载该文件。
补充图10:用于谱图处理的nmrproc.com文件。 nmrproc.com 文件中用于处理 NMRPipe 格式原始数据的代码如图所示。在进行处理之前,必须设置 tauValues 项和 idxExpmt 项。tauValues 对应 ms 中的 vplist 项,单位为 ms 15N R1ρ 实验(vclist 条目以毫秒为单位重新计算) 15N R1 实验)。关于它们的弛豫延迟长度,idxExpmt 条目对应于 vplist 条目(对于 15N R1 实验为 vclist 条目)的顺序。在脚本中,用于更改谱图相位的位置用红色方框标出。红色数字对应于方案中的相应步骤。 请点击此处下载此文件。
补充图11:核磁共振谱图显示与峰识别。 展示了打开谱图时nmrDraw的概览界面。使用+和-按钮调整Factor值,该值会伴随等高线水平的变化而变化。更改Z维度可切换带有不同弛豫延迟的HSQC实验。使用Draw按钮可使用修改后的参数重新绘制谱图。P0滑块和P1滑块用于调整谱图的相位。这些相位调整不会保存在处理后的数据中;必须手动将相位参数输入nmrproc.com文件,并重新执行该文件。使用Peak Detection宏可在设定的等高线水平下识别谱图中的峰。红色数字对应实验方案中的相应步骤。请点击此处下载该文件。
补充图12:使用NMRDraw进行半自动化的核磁共振峰归属。 显示了ass.com文件的代码,必须执行该代码以启动核磁共振信号的归属。将- assName更改为所生成的目标峰列表。该峰列表的格式必须与文件夹peaklist中附带的峰列表ass.syb2.hn_1200.tab一致(可在https://www.ipb.hhu.de/en/teams/team-lakomek/pulsesequences下载的zip文件夹中找到)。红色数字对应于实验方案中的相应步骤。请点击此处下载该文件。
补充图13:核磁共振峰的可视化及化学位移指认检查。 运行 ass.com 可打开化学位移指认宏的界面。使用鼠标左键指认峰,使用右键取消指认。点击 Next 进入下一个氨基酸,点击 Prev 返回上一个氨基酸。通过 + 和 - 按钮可调整等高线水平,再使用 Draw 重新绘制谱图。按照 Quit、Save 和 Exit 的流程操作,可保存指认结果并退出程序。图中红色数字对应实验方案中的相应步骤。请点击此处下载该文件。
补充图14:从 R1ρ 和 R1 数据中提取 15N R2 速率常数。步骤1。 显示了 t1rho.tab 文件的内容。在根据 R1ρ 速率计算 R2 速率之前,必须添加用红色框标出的代码。红色数字对应于方案中的相应步骤。请点击此处下载该文件。
补充图15:从R1ρ和R1数据中提取15N R2速率常数。步骤2。 显示了crtR1p.tcl文件的内容。在执行之前,必须按照实验方案中的说明修改红色框中的条目。红色数字对应于方案中的相应步骤。请点击此处下载该文件。
本方案描述了Lakomek等人69和Stief等人70建立的NMR 15N弛豫实验方法。我们重点采用了基于灵敏度增强型HSQC检测方案的NMR脉冲序列。15N R1和R1ρ实验的实施方法详见Stief等人70的描述,而异核NOE(hetNOE)实验则由Lakomek等人69进行了详细说明。
在设置核磁共振弛豫实验时,15N R1ρ 实验最为关键,因为选择过长的最大弛豫时间或过高的自旋锁定射频幅度可能会损坏探头。对于弛豫时间,理想情况是从100%强度到25%强度(相对于参考实验)均匀采样完整的衰减曲线。对于弛豫较慢的蛋白质(如内在无序蛋白IDP),由于其具有较长的15N T2时间,通常难以对衰减曲线进行全程采样,因为需要限制探头内的功率沉积。
对于15N射频(RF)幅度,理论上RF幅度νspin-lock越高越好,因为可以重聚焦更多的化学交换过程。一般经验法则是,交换速率慢于1 / (2π νspin-lock)的过程可以被重聚焦。然而,实际中受到射频功率的限制。例如,在我们的实验中,为确保安全,当施加2 kHz的RF幅度时,自旋锁定时间不会超过65 ms。同时,在600 MHz核磁共振谱仪上,自旋锁定RF幅度应高于1 kHz,例如600 MHz时使用1.4 kHz,该条件已被证明足以满足小球状蛋白gb369的实验需求。为了覆盖更宽的谱宽并避免在更高磁场强度下的失谐效应,相应的RF幅度应随磁场强度成比例增加,例如在900 MHz时应提高至1.5倍(νspin-lock
2 kHz)。但该幅度必须保持在探头允许的功率限制范围内。(请参考本地谱仪功率限制清单,并联系您的核磁共振设施管理人员以获取建议。)
在选择参考实验(最短延迟)与最长延迟时间(由上述因素决定)之间的弛豫延迟间隔时,我们倾向于对弛豫延迟的数据点采用等间距设置,并选择8个数据点(在分配的测量时间和信噪比允许的情况下)。第一个增量(以及后续增量)主要由信号最强(且通常弛豫最慢)的残基主导。因此,根据蛋白质的动态范围,在选定的最大弛豫时间后,蛋白质序列中多个残基的核磁共振信号可能已经衰减。通过等间距设置延迟点,可确保快速弛豫残基的弛豫曲线也能被足够数量的数据点充分采样。
此外,选择合适的恢复延迟时间对于获得准确的弛豫数据并避免探头长期损坏至关重要。我们选择至少 D1 = 2 s 的恢复延迟,以防止对探头造成任何损伤。为了获得准确的结果,蛋白质的1H 和 15N 纵向极化也必须有足够的时间恢复至平衡状态,或至少达到一个稳定、稳态的状态。在 600 MHz 条件下,小至中等尺寸蛋白质(分子量不超过 20 kDa)的1H T1 通常介于 1 s 到 4 s 之间;对于更大的蛋白质和更高的磁场强度,1H T1 还将进一步增加。对于小球状蛋白 gb3,我们发现 D1 = 3.5 s 已足以满足15N R1 和 R1ρ 实验的要求。对于本方案中介绍的小型固有无序蛋白 Syb-2,我们近期的研究表明 D1 = 2 s 已经足够70。
由于激发起始于15N,且脉冲序列开始时未应用1H-15N INEPT,因此15N hetNOE的灵敏度最低。此外,为了获得准确的数据,1H饱和期和15N恢复延迟应设置为15N T1的5倍。对于本文所展示的固有无序蛋白Syb-2,5秒的饱和/恢复延迟加上额外2秒的D1恢复延迟(总计7秒)已足够;但对于更大的蛋白质和/或氘代蛋白质,该时间可能显著更长,例如10秒或更长。然而,如果定性数据已足够,则可接受轻微的系统性偏差。在此情况下,当延迟时间为15N T1的3倍时即可获得合理结果,因为此时已实现95%的饱和。这种情况下,所得的hetNOE值将表现出系统性偏高的偏差。
在本实验方案中,我们使用了完全质子化的15N标记的Syb-2,这是最具成本效益的选择。建议进行侧链氘代,以减少弛豫路径的数量,从而获得最孤立的1H-15N骨架酰胺自旋对。然而,通过比较15N2H标记的Syb-2与15N标记的完全质子化Syb-2,我们发现氘代样品与完全质子化样品之间的差异极小70。出于实际目的,核磁共振(NMR)弛豫数据通常在双标记或三标记样品(15N13C或2H15N13C)上采集。所讨论的NMR脉冲序列能够重聚焦13C与1H或15N核之间的J耦合演化。显然,引入13C标记后,弛豫路径也会增加。然而,据我们经验,从实际应用角度看,由此引入的误差(如果存在)小于5%,在仅需定性数据的情况下可能是可接受的。
在检测不同大小的蛋白质时,常常需要考虑选择灵敏度增强的HSQC检测实验70还是TROSY检测实验69。对于小分子量蛋白质,灵敏度增强的HSQC检测核磁共振实验可提供两倍的信噪比提升(前提是能够实现良好的水信号抑制)。而对于大分子量蛋白质,则应优先选择TROSY检测实验。对于折叠蛋白,两者的临界点大约在20 kDa左右。对于内在无序蛋白(IDPs),由于其R2弛豫速率常数较低,且大致随肽链长度线性增加,该临界点可能向更高的分子量偏移。在实际应用中,若难以确定选择,建议同时实施灵敏度增强的HSQC检测和TROSY检测实验,然后根据首个增量数据中获得的信噪比优劣来选择更合适的实验方法。
关于选择灵敏度增强的HSQC检测实验还是TROSY检测实验的讨论,仅在信噪比为主要考量因素时适用。当需要最窄的谱线宽度和最高分辨率时,TROSY始终更具优势。TROSY还能提供最佳的水信号抑制效果,因其所需的质子脉冲更少,且无需进行质子去耦。有关更详细的讨论,请参见先前发表的文章69。
关于水峰抑制问题,还有一点需要说明。虽然在Bruker Avance NEO谱仪上将接收增益(rg)设为16可能已足够,但对于(较旧的)Bruker Avance谱仪,通常建议接收增益在128–256范围内。可实现的rg值(由Bruker命令rga确定)通常根据残余水峰信号强度与噪声水平的比值获得。水峰抑制效果不佳总会给谱图引入额外噪声。因此,即使在Bruker Avance NEO谱仪上能够达到rg = 16,残余的水峰信号仍可能显著增加噪声水平,从而导致信噪比下降。因此,我们建议将rg值设定在128–256范围内。文中提供的脉冲序列70已在600 MHz和1200 MHz条件下测试。两种磁场强度均使用脉冲序列中所描述的相同梯度设置。我们至少推荐这些参数作为进一步优化的良好起点。
如引言中所述,所讨论的核磁共振(NMR)弛豫实验仅涵盖皮秒到纳秒时间尺度的动力学过程。对于微秒至毫秒时间尺度的动力学研究,可参考弛豫弥散实验、CPMG实验,以及通过改变自旋锁定幅度(或载波频率)扩展的R1ρ实验44,46,58,59,60,61,85,86。慢于微秒的动力学过程可通过15N CEST(化学交换饱和转移)NMR62、交换谱学(EXSY,毫秒至秒级)或实时(Real-time, RT)NMR(秒至小时级)63,64进行检测。顺磁探针引起的顺磁弛豫增强(PRE)效应以及残余偶极耦合(RDCs)可用于评估从皮秒到毫秒整个时间范围内的动力学行为65,66,67,68。
所有这些实验仅涵盖主链动力学。为了研究侧链动力学,Kay 及其同事引入了另一系列核磁共振(NMR)实验,可提供分子量高达 1 MDa 的蛋白质侧链动力学信息2,87。
为了更定量地了解蛋白质主链的动力学特性,有必要对全部三种核磁共振(NMR)参数——15N R1、R2 和 hetNOE——进行联合分析,有时甚至需要引入更多参数,例如交叉弛豫速率常数。目前已发展出多种复杂的分析方法。其中所谓的“无模型”(model-free)方法已在生物大分子NMR中得到广泛应用43,88,89,90,91。此外,谱密度映射也可提供有价值的信息55,92,93。尽管这些方法对于具有明确旋转相关时间的球状蛋白质非常有效,但其在内在无序蛋白(IDP)中的适用性受到限制。定义旋转相关时间需要一个折叠的结构域39。当所涉及的不同时间尺度能够良好分离时,可通过将相关函数 C(t) 建模为多个指数衰减函数之和的方式部分克服这一限制23,31,84。在缺乏明确旋转相关时间的情况下,可采用涉及相关时间分布的替代方法,例如 IMPACT 分析94 或 DETECTOR 分析95,96。此外,将NMR弛豫数据与分子动力学(MD)模拟相结合的分析方法对IDP也具有前景97;然而目前在较大IDP的全原子模拟中,我们仍面临采样和力场方面的局限性31,98。
作者声明,他们不存在任何已知的可能被认为会影响本文报道工作的竞争性经济利益或个人关系。
感谢Melinda Jaspert和Kevin Bochinsky提供的有益讨论。N.L.感谢德国科学基金会通过海森堡计划(DFG资助编号433700474)提供的经费支持。本工作还得到了亥姆霍兹协会倡议与网络基金资助项目"新冠病毒致病机制的病毒学与免疫学决定因素——为未来大流行病做好准备的经验教训(KA1-Co-02 "COVIPA")的支持。我们感谢获得由于利希研究中心和杜塞尔多夫海因里希·海涅大学(HHU)共同运营的于利希-杜塞尔多夫生物分子核磁共振中心的大力支持。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 布鲁克 600 MHz AVANCE III HD 波谱仪 | Bruker | https://www.bruker.com/en/products-and-solutions/mr/nmr/avance-nmr-spectrometer.html | 进行的核磁共振实验 |