芳香环翻转揭示晶体和复合物中蛋白质动力学的重塑

《Nature Chemistry》：Aromatic ring flips reveal reshaping of protein dynamics in crystals and complexes

　　

蛋白质构象能量景观不仅受分子内相互作用影响，还受其环境影响。在蛋白质晶体和蛋白质-蛋白质复合物中，分子间接触改变了这一能量景观，但这种改变的确切性质难以解读。理解晶格如何影响蛋白质动力学对于基于晶体学的运动研究至关重要，但其对集体运动的影响仍不清楚。疏水核心中的芳香环翻转是此类动力学的敏感探针。本研究结合先进同位素标记与定量核磁共振（NMR）方法，比较了GB1蛋白在晶体、与其结合伙伴免疫球蛋白G（IgG）的复合物以及溶液中的芳香环翻转动力学。研究人员发现，晶体中核心环的翻转频率比溶液中低近三个数量级。基于本研究报道的GB1变体晶体结构的增强采样分子动力学模拟重现了这些升高的能垒，并揭示了晶体如何约束运动。

### 论文解读：芳香环翻转揭示晶体和复合物中蛋白质动力学的重塑

#### 研究背景、问题与意义

蛋白质的构象能量景观不仅由分子内相互作用决定，还受其环境（如晶体晶格或结合伙伴）的显著影响。在蛋白质晶体学中，长期以来存在一个关键问题：晶体堆积对蛋白质动力学的约束程度相对于溶液状态究竟如何？尽管早期研究表明晶格可改变蛋白质的结构与动力学，但其对集体运动（尤其是慢时间尺度运动）的影响仍不明确。芳香环的翻转——苯丙氨酸（Phe）和酪氨酸（Tyr）侧链绕Cβ–Cγ轴（χ²角）的180°旋转——被认为是集体动力学的敏感探针，因为埋藏环的翻转需要蛋白质发生大的协同“呼吸运动”。然而，晶体学无法直接检测这些翻转，而核磁共振（NMR）波谱学可有效探测。本研究聚焦于免疫球蛋白结合蛋白G（GB1）的芳香环翻转，系统比较其在溶液、晶体及与天然结合伙伴免疫球蛋白G（IgG）复合物中的动力学，旨在揭示分子间接触如何重塑蛋白质内部运动的时间尺度。该工作发表在《Nature Chemistry》。

#### 关键技术方法

研究人员采用了整合实验与计算的方法。关键实验技术包括：（1）先进同位素标记：利用α-酮酸前体对Phe和Tyr的(CH)^?位点进行位点特异性¹³C¹H标记，并结合均匀²H和¹⁵N标记；（2）定量NMR波谱学：包括魔角旋转（MAS）NMR（用于晶体与复合物样品）和溶液NMR，运用交换谱（EXSY）、旋转回波双共振（REDOR）、纵向和旋转坐标系¹³C自旋弛豫测量（R₁和R_1ρ）以及检测器（detectors）分析方法；（3）单晶X射线衍射（XRD）：分别在100 K和295 K解析了GB1变体GB1^QDD（T2Q, N8D, N37D三重突变）的晶体结构（空间群C121，不对称单元含两个分子）。计算方面，采用增强采样的分子动力学（MD）模拟（基于多行走子温控元动力学扩展自适应偏置力（WTM-eABF）算法）获取自由能垒。所有模拟均基于AMBER ff19SB力场和OPC水模型，温度维持在298 K。

#### 研究结果

**GB1^QDD在晶体和IgG:GB1^QDD复合物中的堆积**
通过单晶XRD解析了GB1^QDD的晶体结构（1.08 ?分辨率，100 K；1.69 ?，295 K），发现晶体中每个Cα原子10 ?范围内的非氢原子数在两条链中相似，但与IgG复合物模型中的接触分布不同，表明局部堆积差异可能调节蛋白动力学。

**GB1^QDD晶体的MAS NMR**
对晶体样品进行MAS NMR谱图归属，观察到底物Q32–D40区域信号缺失，归因于动力学导致的线宽增强；而GB1^T2Q晶体中该区域可观测。15N自旋弛豫检测器分析显示，晶体中残基28–31和41的μs–ms运动显著增强，进一步证实Q32–D40区域的微秒运动加剧。与GB1^T2Q晶体相比，GB1^QDD晶体中快速运动（ρ₁）类似，但整体ρ₁和ρ₂升高，表明突变影响了动力学。

**GB1^QDD三种状态下芳香环的NMR特征**
高分辨MAS NMR谱图显示，溶液中每个芳香残基因快速翻转（k_flip > Δω）仅出现一个平均信号；IgG:GB1复合物中同样为单峰；而晶体中Y3和Y45各出现两个等间距分裂峰，表明翻转变慢（k_flip < Δω）；Y33保持快速翻转单峰；Phe区域则出现宽峰，指示微秒运动。

**晶体堆积减慢埋藏酪氨酸环的翻转**
通过¹H–¹³C EXSY实验定量得到Y3^crystal和Y45^crystal的翻转速率分别为(55.8 ± 2.2) s^-1和(18.37 ± 0.11) s^-1（304 K），比溶液中约慢三个数量级。温度依赖实验和联合Eyring拟合表明，晶体中翻转减慢源于激活熵降低（ΔS^?,solution > ΔS^?,crystal），而激活焓相近。

**通过MAS NMR方法探测快速翻转环**
对于快速翻转的Y33和Phe残基，利用REDOR测得的偶极序参数（S）和¹³C自旋弛豫结合检测器分析方法推断：Y33^crystal的翻转时间尺度在37–600 ns之间；Phe信号（F30和/或F52）对应ns–μs运动。

**IgG:GB1^QDD中环的翻转比晶体中更快**
IgG:GB1复合物中所有酪氨酸均呈现单峰，且偶极序参数低、¹³C R₁和R_1ρ弛豫速率高，检测器分析表明翻转时间尺度在数十至数百纳秒，比晶体中快约五个数量级。

**晶体中环翻转的熵有利性降低**
综合三种状态的翻转相关时间（τ_flip）和温度依赖数据，得到Y3和Y45在晶体中的活化吉布斯自由能（ΔG^?）高于溶液，差值分别为(15 ± 5) kJ mol^-1和(17.9 ± 1.5) kJ mol^-1，与MD模拟结果(13 ± 5)和(17 ± 4) kJ mol^-1高度吻合。

**MD模拟揭示环动力学涉及复杂相互作用网络**
增强采样MD模拟在溶液和晶体中均重现了Y3和Y45的高能垒，且与实验定量一致。模拟发现，Y45的翻转与χ¹角的波动以及相邻晶胞分子间接触（如D40、E42、W43、E56）的瞬时变化相关，但未观察到与疏水核心整体膨胀/收缩（全局“呼吸运动”）的明显耦合。这些分子间接触在溶液中不存在，因此被认为是晶体中翻转能垒升高的主要来源。

#### 总结与讨论

研究结论表明，分子间接触可深刻影响内部蛋白质动力学的时间尺度。晶体堆积使埋藏酪氨酸（Y3、Y45）的翻转减慢约三个数量级，这归因于颗粒间接触限制了翻转所需构象的熵可利用性；而IgG结合则意外地加速了翻转，可能是由于GB1–IgG相互作用稳定了利于翻转的构象或复合物整体的低幅摇摆运动。此外，苯丙氨酸（F30、F52）在晶体中翻转明显快于酪氨酸，提示不同芳香残基的翻转机制存在差异，不支持单一的全局“呼吸运动”模型。这些发现强调蛋白质自由能景观的复杂性和对环境的高度依赖性，对基于晶体学的动力学研究以及细胞拥挤环境中的动力学具有重要启示。