《JACS Au》：Solvent Reorganization in Stabilized Protein–Polymer Conjugates Visualized by Two-Dimensional Infrared and Nuclear Magnetic Resonance Spectroscopy
编辑推荐：
本文通过整合二维红外（2D IR）光谱、多维核磁共振（NMR）光谱及分子动力学（MD）模拟技术，深入探究了聚乙二醇化（PEGylation）如何通过改变蛋白质溶剂化壳层动力学而非脱水作用来增强蛋白质稳定性。研究发现，较长PEG链（5 kDa）能形成“帷幕”（shroud）状构象，减缓溶剂重排动力学，从而显著提升人半乳糖凝集素-3碳水化合物识别域（Gal3C）的热稳定性，为理性设计具有可预测性能的生物制剂提供了关键指导。
背景与意义
蛋白质类药物（生物制剂）因其高特异性和有效性，在医药领域占据重要地位，但其固有的复杂性也使其在体内环境中易降解、清除，从而限制了疗效。聚乙二醇化是提升生物制剂稳定性的关键策略之一，它通过共价连接聚乙二醇（PEG）来延长药物的血液循环时间。然而，PEG化如何从分子机制上增强蛋白质稳定性仍不完全清楚，这阻碍了具有可预测性能的蛋白质-聚合物偶联物的理性设计。其中，溶剂化作为蛋白质折叠与稳定的关键驱动力，被假设在PEG化介导的蛋白质稳定中扮演核心角色。
传统上存在两种解释PEG-蛋白质相互作用的模型：哑铃模型和帷幕模型。在哑铃模型中，PEG链大部分时间伸展在溶液中，与蛋白质表面仅有有限的直接接触，两者更像是相连但独立的单元。在帷幕模型中，PEG链则紧密贴合在蛋白质表面，形成一个动态的、与水结合的“帷幕”层。这两种模型预测了不同的蛋白质-聚合物相互作用程度，进而可能对蛋白质的溶剂化环境产生不同影响。
研究系统与方法
本研究以人半乳糖凝集素-3的碳水化合物识别域（Gal3C）为模型蛋白。Gal3C具有β-夹层折叠结构，在生物学和医学中具有代表性。研究人员制备了未修饰的Gal3C以及分别连接了两种不同分子量PEG的偶联物：Gal3C-PEG 550 Da（短链）和Gal3C-PEG 5 kDa（长链）。两者通过马来酰亚胺功能化PEG与Gal3C第243位点突变的半胱氨酸进行特异性偶联。
为了从多时间尺度和多维度揭示机制，研究整合了三种核心技术：
- 1.
二维红外光谱（2D IR）：用于探测皮秒至纳秒时间尺度的蛋白质主链动力学和溶剂动力学。
- 2.
多维核磁共振光谱（NMR）：提供了较慢时间尺度的蛋白质主链构象和动力学信息，并用于测量偶联物的旋转相关时间（τc），以推断其扩散性质和整体构象模型。
- 3.
分子动力学模拟（MD）：用于可视化PEG链的构象，并全局表征溶剂化环境和蛋白质主链柔性，为实验结果提供原子层面的解释。
研究结果与发现
1. 热稳定性与构象模型的差异
热稳定性测试表明，Gal3C-PEG 550 Da与未修饰的Gal3C类似，呈现单一的协同去折叠转变，而Gal3C-PEG 5 kDa则显示出两个独立的去折叠转变并形成中间态，表明其热稳定性更高且去折叠路径被重塑。
NMR旋转相关时间测量为区分构象模型提供了关键证据。Gal3C-PEG 550 Da的τc值与未修饰的Gal3C相似，而Gal3C-PEG 5 kDa的τc值则显著增加，这表明长链PEG显著减慢了偶联物的整体旋转扩散。这一结果支持Gal3C-PEG 5 kDa更倾向于形成“帷幕”构象，其中PEG链与蛋白质表面发生更广泛的相互作用；而Gal3C-PEG 550 Da则更接近“哑铃”模型。
进一步的NMR化学位移扰动和线宽分析也证实了这一点。与Gal3C-PEG 550 Da相比，Gal3C-PEG 5 kDa在许多残基上表现出更显著的化学位移扰动和更广泛的线宽增加，表明长链PEG对蛋白质的构象和动力学产生了更深远的影响，特别是在134-140、193-200、219和249等区域的残基附近。
2. 超快时间尺度上的蛋白质环境变化
2D IR光谱通过对酰胺I带的分析，揭示了PEG化在超快时间尺度上对蛋白质主链环境的影响。光谱显示，Gal3C-PEG 5 kDa的酰胺I正峰（1630和1680 cm-1）比未修饰的Gal3C和Gal3C-PEG 550 Da表现出更显著的沿对角线拉长，表明长链PEG诱导了更大的结构异质性和/或溶剂暴露增加。
中心线斜率（CLS）分析用于量化频率-频率相关函数（FFCF）的衰减，这反映了主链周围环境的波动速率。结果显示，Gal3C-PEG 5 kDa的CLS衰减时间常数几乎是Gal3C和Gal3C-PEG 550 Da的两倍，意味着其局部动力学（包括氢键重排）显著减慢。这表明长链PEG通过改变溶剂化网络的动力学，稳定了蛋白质周围的溶剂环境。
3. 分子动力学模拟的可视化与验证
MD模拟直观地展示了短链和长链PEG与Gal3C相互作用的差异。短链PEG（600 Da）仅覆盖连接位点附近的一小片区域，像一个“帽子”。而长链PEG（5 kDa）则更广泛地覆盖了蛋白质表面（特别是背面和侧面），形成了一个像“背包”一样的“帷幕”层。
接触概率分析表明，长链PEG的原子与蛋白质残基（特别是131-141和191-202区段）发生持续接触的概率远高于短链PEG，这与NMR观测到的化学位移扰动区域高度吻合。同时，所有PEG化蛋白的残基波动性均降低，表明PEG化对蛋白质施加了全局性的构象约束。
4. 蛋白质表面溶剂可及性与氢键动力学
一个重要的发现是，PEG化并未导致蛋白质表面脱水。通过计算溶剂可及表面积（SASA），研究发现三种蛋白（Gal3C、Gal3C-PEG 600 Da、Gal3C-PEG 5 kDa）的总SASA以及蛋白主链与水形成的氢键数量均无明显差异。这意味着PEG化稳定蛋白质的机制并非通过“脱水”或屏蔽蛋白质表面来实现。
然而，MD模拟计算的关键氢键寿命揭示了本质区别。虽然蛋白质主链与水分子形成的氢键总数相近，但这些氢键的寿命在Gal3C-PEG 5 kDa中显著延长了约50%。这直接解释了2D IR中观察到的溶剂动力学减慢现象。进一步分析表明，这种氢键寿命的延长同时源于酰胺-水（分子间）氢键和酰胺-酰胺（分子内）氢键寿命的增加，但前者的动态变化幅度更大，对皮秒时间尺度的FFCF贡献可能更强。
讨论与结论
本研究整合多尺度、多技术手段，挑战了“PEG化通过脱水蛋白质表面来稳定蛋白质”的传统观点，提出了一个新的分子机制：PEG化，特别是当使用足够长的聚合物并形成“帷幕”状构象时，是通过稳定并重组蛋白质周围的溶剂化壳层，减缓溶剂重排动力学，从而增强蛋白质的热稳定性。
具体而言，较长的PEG链（如5 kDa）能够更有效地与蛋白质表面相互作用，形成动态的“帷幕”。这种相互作用改变了溶剂化网络的结构和动态，延长了蛋白质主链与周围水分子之间氢键的寿命，减缓了局部的环境涨落。这种“稳定化的溶剂环境”提升了蛋白质的能垒，使其更不容易去折叠。
这一机制具有重要的指导意义。它强调了聚合物长度（达到形成“帷幕”的阈值）对于实现有效稳定的重要性。对于单一位点偶联，需使用足够长的PEG；而对于多位点偶联，则需达到足够的修饰度以实现充分的表面“覆盖”。同时，蛋白质自身的性质（如大小、表面特性、连接位点）也会影响PEG的构象行为和最终的稳定效果。
综上所述，本研究通过可视化溶剂重组过程，为理解蛋白质-聚合物偶联物的稳定性机制提供了深刻的分子见解，并为未来理性设计性能可控的生物制剂提供了坚实的理论和实验基础。
