波动β-折叠：骨架N-氨化对DS119热稳定性影响甚微的结构动力学解析

《Journal of Peptide Science》：Fluctuating β-Sheet Secondary Structure in DS119 Explains the Small Effects of Backbone N-Amination on Thermal Stability

为了解析骨架N-氨化对平行β-折叠稳定性影响不符预期的原因，研究人员以微型蛋白DS119为模型，结合化学合成、圆二色光谱（CD）及受核奥弗豪泽效应（NOE）约束的分子动力学（MD）模拟开展研究。结果发现DS119的β-折叠区域具有高度波动性，这解释了N-氨化未能显著提升稳定性的现象。该研究强调了MD模拟在揭示蛋白质结构动态特征中的重要性。

　　

在生物化学与药物设计领域，科学家们一直在探索如何通过修饰蛋白质的骨架来赋予其新的功能或提升其稳定性。其中，骨架N-氨化（将肽键中的氮原子替换为带氨基的氮原子，形成α-氨酰肼酸残基）是一种颇具潜力的策略。这种修饰不仅能改变氢键模式、主链构象偏好，还能增强对蛋白酶解的稳定性，同时保留天然侧链的多样性，因此在理性设计新型蛋白质模拟物（proteomimetics）方面前景广阔。此前的研究表明，在反平行β-发卡肽中引入N-氨化的残基能够增加其折叠比例。那么，一个自然而然的问题是：这种稳定效应能否推广到更复杂的平行β-折叠结构，特别是在一个具有三级结构的微型蛋白质中呢？

为了回答这个问题，一个国际研究团队选择了人工设计的微型蛋白质DS119作为模型系统。DS119具有一个βαβ的三级折叠结构，包含一个两亲性螺旋和由N端与C端区域形成的两个平行β-链，是一个研究蛋白质折叠、稳定性和工程化的经典模型。研究人员最初假设，在平行β-折叠的外缘引入N-氨基取代基能够增强其热稳定性。然而，实验结果却出乎意料。

为了验证假设，研究人员合成了DS119的野生型及其一系列衍生物。这些衍生物在平行β-链外缘的色氨酸9（Trp9）或苯丙氨酸33（Phe33）位点引入了N-氨化修饰，并在相同位点合成了N-甲基化衍生物作为对比。此外，还在β-链内侧氢键结合的色氨酸34（Trp34）位点引入了N-氨化修饰，预期这会破坏结构。他们利用化学合成方法制备了这些修饰后的多肽，并采用圆二色光谱（Circular Dichroism, CD）分析了它们的二级结构和热稳定性。

关键实验技术主要包括：固相肽合成（SPPS）用于构建目标多肽序列；圆二色光谱（CD）用于表征蛋白质的二级结构并监测其热变性过程，通过测量平均残基椭圆度（MRE）随温度的变化来推算熔解温度（T_m）；分子动力学（Molecular Dynamics, MD）和随机动力学（Stochastic Dynamics, SD）模拟，利用GROMOS软件和不同力场（54A7, 54A8, 54B7），在显式溶剂和隐式溶剂条件下，对蛋白质结构进行模拟，其中部分模拟引入了从核磁共振（NMR）结构（PDB: 2KI0）中获得的核奥弗豪泽效应（NOE）距离约束，以精修结构并分析其动态特性。

3.1 在DS119中引入α-氨酰肼酸残基

研究人员选择了Trp9、Phe33和Trp34（位于色氨酸拉链基序附近）作为初始修饰位点。他们合成了野生型DS119和六个变体，包括在Trp9、Phe33位点的N-氨化和N-甲基化变体，以及在Trp34位点的N-氨化变体。CD光谱分析表明，化学合成的野生型DS119保持了重组蛋白的βαβ结构。N-氨化变体对CD谱图扰动较小，而N-甲基化变体显示出向无规卷曲结构转变的迹象。热稳定性实验显示，野生型和Trp9、Phe33的N-氨化变体在无变性剂条件下在95°C以下均未显示明显的去折叠转变，表明它们都具有高稳定性。相比之下，Trp34的N-氨化变体熔解温度（T_m）为82.5°C，而两个N-甲基化变体稳定性最低（T_m分别为74.4°C和78.7°C）。为了更精确量化稳定性差异，在1M盐酸胍（Gdn-HCl）存在下进行了耦合热/化学变性实验。结果显示，Trp9 N-氨化变体与野生型稳定性相似（T_m分别为61.0°C和62.0°C），而Phe33 N-氨化变体稳定性较低（T_m为51.0°C）。这与N-氨化能稳定平行β-折叠的初始假设不符。

3.2 DS119的MD和SD模拟

为了理解上述意外结果，研究人员对DS119结构进行了NOE约束的MD和SD模拟精修。模拟从已发表的DS119 NMR结构（PDB: 2KI0）开始。无约束的显式溶剂MD模拟显示出相当程度的去折叠，且与实验NOE数据存在大量严重冲突，表明起始结构中的芳香侧链堆积可能不够稳定。因此，他们进行了包含NOE约束的模拟，这些模拟很好地满足了实验NOE数据。在所有NOE约束的显式溶剂模拟中，残基14-26形成了至少持续90%模拟时间的α-螺旋。然而，残基6-10和30-34的β-折叠二级结构占比很低（0%-25%）。虽然这些残基在大部分轨迹中采取伸展构象，但链间氢键的占比很低。模拟时间序列显示β-折叠在整个模拟过程中是瞬时存在的。

对氢键的详细分析表明，在原始的20个NMR结构 ensemble 中，大多数结构只存在一两个链间氢键。而在NOE约束的MD模拟中，观察到了平行β-折叠预期的六个氢键，其中8NH-30O、30NH-6O和34NH-10O的占比最高（24%–77%），其余三个占比很低（1%–8%）。隐式溶剂模拟显示了类似的氢键模式，但还存在不属于典型平行β-折叠的额外氢键（如12NH-35O）。因此，DS119的NOE精修结构显示其平行β-折叠比原始NMR结构观察到的更为规则，但具有高度波动的特性。此外，Trp9和Trp34的侧链始终保持平行排列，这与NMR观测到的平行相互作用而非典型的边对面（edge-to-face）相互作用一致。

讨论与结论

该研究通过综合实验与模拟手段，揭示了DS119微型蛋白中一个未被充分认识的结构动态特征。最初的设计旨在通过色氨酸拉链基序稳定一个平行β-折叠，但本研究显示，尽管该基序成功地将蛋白质的N端和C端区域拉近，却并未形成一个持久、稳定的平行β-折叠。相反，β-折叠区域的主链氢键处于高度波动状态，其占比从低于10%到最高77%不等。这种固有的结构柔性（或称之为“波动性”）意味着，在此区域引入旨在稳定特定骨架构象（如伸展的β-折叠构象）的N-氨化修饰，其效果必然有限。因为修饰所倾向的构象只在部分时间存在，其稳定化效应被整个结构的动态平均所稀释。这完美地解释了为何在Trp9和Phe33位点进行N-氨化未能如预期那样显著提升蛋白质的热稳定性。其中，aPhe33变体稳定性反而下降，可能是由于骨架构象的改变影响了苯丙氨酸侧链的朝向及其与其他疏水残基的相互作用。

这项研究具有双重重要意义。首先，在应用层面，它提醒蛋白质设计者需要全面考虑目标结构的动态特性。旨在引入特定相互作用的理性设计，必须评估该结构区域是刚性的还是柔性的，因为动态区域会削弱局部修饰的预期效果。其次，在方法论上，该研究有力展示了分子动力学模拟，特别是结合了实验约束（如NMR的NOE数据）的模拟，在精修蛋白质结构和揭示其重要动态特征方面的强大能力。这些动态特征在标准的、基于“单一确定构象”假设的蛋白质结构解析流程中很容易被忽略。该论文发表在《Journal of Peptide Science》上，其工作强调了构象无序性和熵贡献对设计蛋白质稳定性的重要性，为未来开发功能更优的蛋白质模拟物提供了关键的动力学视角。