使用AMOEBA揭示胶原模拟肽中的变形与水合各向异性

《The Journal of Physical Chemistry B》：Using AMOEBA to Uncover Deformation and Hydration Anisotropy in Collagen-Mimetic Peptides

　　

胶原是结缔组织中超分子纤维的主要成分。为了能够预测胶原的性质，需要在分子尺度上进行严谨的表征。本工作聚焦于由(PPG)三肽重复序列折叠成三螺旋的胶原模拟肽（collagen-mimetic peptides, CMPs）。研究人员使用AMOEBA极化力场在生理温度下对这些CMPs进行分子动力学（molecular dynamics, MD）模拟。研究表明，AMOEBA能够捕捉三螺旋的长度依赖性 fraying（末端解链），与实验观察一致。研究人员还引入了新的度量指标来量化三螺旋的变形和水合动力学，从而能够系统量化多种条件下的CMPs，填补了该领域此前缺失的部分。研究人员将该方法应用于突变CMPs以及存在d-葡萄糖的CMPs中，发现平移扩散各向异性是CMP结构完整性的分子特征。

论文解读：基于AMOEBA力场的胶原模拟肽变形与水合各向异性研究

一、研究背景与问题

胶原蛋白是哺乳动物体内最丰富的结构蛋白，构成肌腱、韧带、皮肤和骨骼等结缔组织的主要成分。其分子层面的基本结构单元为三螺旋：三条左手聚脯氨酸II型多肽链缠绕形成右手超螺旋。由于天然胶原每条链超过一千个氨基酸，长度可达300 nm以上，实验和计算研究常采用基于(Pro-Pro-Gly)_n重复序列的胶原模拟肽（collagen-mimetic peptides, CMPs）作为简化模型。这些CMPs保留了三螺旋的关键几何特征，同时降低了序列异质性和层级复杂性，便于研究氢键、水合、链堆积以及局部化学扰动对三螺旋结构稳定性的影响。

然而，现有研究存在若干关键问题。实验和晶体学研究表明，(PPG)_n三螺旋的熔解温度强烈依赖链长，n ≤ 12时熔解温度低于生理温度310 K，而许多传统分子动力学模拟在310 K下却无法再现短CMPs的变形行为。例如，CHARMM36+TIP3P（C3P）模拟使短CMPs保持三螺旋结构，与实验变形研究矛盾。此外，传统非极化力场可能过度估计蛋白质-蛋白质相互作用而低估蛋白质-水相互作用，且难以准确描述羰基之间的非共价n → π*相互作用。这些局限性促使研究人员探索使用显式极化力场AMOEBA来模拟CMPs，以期更真实地捕捉其结构动态和水合行为。

二、研究内容与主要结论

研究人员利用AMOEBA极化力场，对三种长度（(PPG)₅、(PPG)₁₂、(PPG)₂₅）的CMPs进行了100 ns的分子动力学模拟，并与C3P、CHARMM36+TIP4P（C4P）、CHARMM36mGP+TIP4P（Cm4P）、AMBER+TIP4P（A4P）等非极化力场结果进行系统比较。同时，研究人员开发了新的量化指标κ_i,m(t)和ζ_i(t)，分别用于表征链内和链间变形，并分析了水分子在三螺旋周围平动和转动扩散的各向异性。此外，还构建了Gly→Ala和Pro→Lys突变体，并在不同浓度d-葡萄糖（0–2000 mM）条件下考察序列和化学扰动对三螺旋的影响。

研究得出以下主要结论：AMOEBA能够再现实验观测的长度依赖性三螺旋末端解链现象，而C3P、C4P和Cm4P则无法做到；AMOEBA对脯氨酸吡咯烷环puckering的预测与Cm4P、A4P相当，且均优于C3P/C4P；AMOEBA对n → π*相互作用的描述相对更好。水合分析发现，CMPs周围水分子的平动和转动扩散均表现出显著各向异性，沿三螺旋主轴方向的扩散慢于垂直方向，这种各向异性在自由水中不存在。Gly→Ala突变引起局部结构变形，而Pro→Lys突变通过带正电的NH₃⁺侧链增强与水及d-葡萄糖的相互作用。d-葡萄糖虽不与CMP直接长时间结合，但通过占据空间、限制链分离，间接稳定三螺旋结构。

三、关键技术方法

本研究主要采用以下关键技术方法：（1）基于THeBuScr脚本生成(PPG)_n三螺旋初始结构；（2）使用AMOEBA极化力场进行分子动力学模拟，并以C3P、C4P、Cm4P、A4P作为参照；（3）开发HeliXplore代码计算链内变形指标κ（基于rise和radius的偏差）和链间变形指标ζ（基于三链Cα三角形面积和形状的偏差）；（4）通过Bürgi–Dunitz距离和角度分析n → π*相互作用；（5）采用Impey方法计算水分子和d-葡萄糖的停留时间；（6）基于速度自相关函数计算水分子平动和转动扩散张量及其各向异性参数。样本均为计算模拟体系，无实验队列来源。

四、研究结果

1. (PPG)_n单链的构象柔性

通过Ramachandran图分析发现，所有力场给出的平均二面角ψ、φ都在10°以内，但AMOEBA模拟中ψ、φ的标准差是其他力场的两倍以上，表明AMOEBA下肽链骨架具有更大的柔性。ω角在所有力场下均为约176°，与实验值178°一致。脯氨酸吡咯烷环puckering分析显示，Cm4P、A4P和AMOEBA能正确预测Xaa位exo:endo为50:50、Yaa位为75:25的群体比例，而C3P和C4P则高估exo群体。对于n → π*相互作用，AMOEBA使δd_BD和δθ_BD同时满足参考条件的帧比例为16%（(PPG)₁₂），高于Cm4P的9.7%和A4P的13%，表明AMOEBA能更灵活地容纳CMPs的变形。

2. 三螺旋变形

通过主成分分析得到三条链的主轴方向偏差δθ和轴向位移s(t)。AMOEBA模拟中，(PPG)₅的?δθ?_t = 0.7 ± 0.2，(PPG)₁₂为0.2 ± 0.1，(PPG)₂₅接近0，呈现明显的长度依赖性；而C3P、C4P、Cm4P对所有长度均给出接近零的变形，不符合实验预期。新指标κ_i,m(t)和ζ_i(t)进一步显示：(PPG)₅所有三肽单元均发生链内和链间变形；(PPG)₁₂前四个N端三肽发生链内变形，仅N端第一个三肽发生链间变形；(PPG)₂₅只有N端第一、二个三肽发生链内变形。变形主要集中于N端，与实验观察的三螺旋解链方向一致。

3. 氢键与水合分析

AMOEBA模拟中，Gly(NH)-Xaa(CO)链间氢键在(PPG)₅、₁₂、₂₅中平均每三肽分别为0.2、0.75和0.8个，断裂主要发生在N端。Xaa(Cα)-CO氢键则会在模拟中动态形成。三肽与水形成的氢键数在末端三肽增多，变形区域暴露更多水分子。水合壳层水分子的平动扩散系数约为1.1–1.2 × 10^-5 cm²/s，比体相水（3 × 10^-5 cm²/s）慢2.5–2.7倍。平动弛豫时间常数τ₁^trans≈0.3 ps，τ₂^trans≈3.15 ps；转动弛豫时间τ^rot≈0.65–0.80 ps，均在蛋白质水合水的文献范围内。关键发现是，水分子沿三螺旋主轴（z轴）的平动扩散D_zz^trans显著低于垂直方向，各向异性参数c_aniso^trans在(PPG)₁₂和(PPG)₂₅中分别为0.44和0.38，而在完全变形的(PPG)₅中接近零。转动扩散同样表现出围绕水偶极轴（r₃轴）较慢的各向异性。

4. 胁迫因素对(PPG)_n CMPs的影响

在(PPG)₁₂中引入Gly→Ala和Pro→Lys突变，并加入0–2000 mM d-葡萄糖。结果显示，d-葡萄糖浓度不影响非突变(PPG)₁₂的氢键数目。Ala突变在0 mM时导致CMP弯曲成L形，但d-葡萄糖存在时弯曲消失，链间变形指标ζ下降；Lys突变显著增加突变位点附近的CMP-水氢键数（从3–5增加到6–10），但水合壳层水分子数因侧链空间排斥而减少。所有CMP在d-葡萄糖存在下，高ζ群体（链间分离态）的比例从0 mM时的5.13%降至2000 mM时的0.07%，表明d-葡萄糖通过空间限制间接稳定三螺旋。Lys的NH₃⁺使d-葡萄糖的最大停留时间延长至5.1–5.7 ns，可能促进非酶糖化反应。平动扩散各向异性c_aniso^trans在健康CMPs中较高，在弯曲的Ala突变体中降低，而转动各向异性基本不受影响。

五、讨论与结论

讨论部分指出，AMOEBA力场能够同时再现CMPs的构象柔性、长度依赖性变形和正确的氢键动态，这归因于其显式极化处理以及更平衡的蛋白质-水相互作用。研究新提出的变形指标κ和ζ可以区分拉伸、压缩和链分离等不同变形模式，具有普适性，可用于其他多链螺旋体系。水合各向异性的发现将三螺旋构象与水动力学联系起来，提示各向异性水扩散可作为胶原样有序结构的敏感分子标志物。在糖化相关研究中，d-葡萄糖通过空间占用稳定三螺旋，而Lys侧链提供了持久结合位点，可能为糖化早期反应创造有利微环境。

结论部分总结：研究人员使用AMOEBA力场系统研究了(PPG)_n CMPs，验证了AMOEBA在预测脯氨酸puckering、长度依赖性变形和n → π*相互作用方面的优势。新指标揭示了变形的位置和类型，发现水分子在CMPs周围存在平动和转动各向异性。突变研究显示Gly→Ala破坏局部结构，Lys引入额外相互作用，而d-葡萄糖整体上稳定三螺旋。该工作为理解胶原稳定性及糖化相关结构变化提供了分子层面的控制模型，未来可扩展至更接近天然胶原的序列和复杂细胞外环境。