《Biophysical Journal》：Molecular dynamics reveals how calcium drives hetero- versus homo-trimerisation of type I collagen
编辑推荐：
本研究针对I型胶原异常同源三聚化与多种年龄相关疾病的关联机制不明的问题，通过分子动力学模拟揭示了钙离子如何通过稳定链间界面、调控二硫键形成及赋予热力学优势来驱动异源三聚化。结果表明，钙离子缺失会破坏氢键网络，增大半胱氨酸间距，阻碍异源三聚体形成，从而促使同源三聚体形成。该发现为理解相关疾病的分子病理提供了新视角。
在人体这个精密的生命体中，I型胶原扮演着“脚手架”的角色，是构成骨骼、皮肤、肌腱等组织的主要结构蛋白，对维持机体的结构完整性和力学性能至关重要。这种蛋白质在合成时，需要由两条proα1α(I)链和一条proα2α(I)链精确地组装成一个异源三聚体分子。这个过程并非自发完成，钙离子的参与是其关键一环。然而，在衰老或某些病理状态下，机体有时会产生一种由三条proα1α(I)链组成的异常同源三聚体。这种异常的胶原分子与多种严重的年龄相关疾病密切相关，包括癌症、组织纤维化、肌肉骨骼系统疾病以及心血管疾病。但是，科学家们长期以来一直困惑：在什么情况下，正常的异源三聚化过程会被打破，转而形成这种可能有害的同源三聚体？钙离子在这一分子“抉择”中究竟扮演了怎样的角色？为了回答这些基础且重要的生物学问题，研究人员展开了一项深入的研究。
为了揭示钙离子在I型胶原三聚化过程中的分子机制，本研究主要采用了以下几种关键的计算机模拟与理论计算方法：
- 1.
分子动力学模拟：对I型胶原C-端肽(C-propeptide)的蛋白质结构进行了长时间的动态模拟，以此观察在钙离子存在和缺失条件下，蛋白质分子内和分子间相互作用的动态变化。
- 2.
拉伸模拟：通过施加外力，模拟钙离子从单体胶原链上解离的过程，用以评估钙离子与不同胶原链（proα1α(I)与proα2α(I)）的结合强度差异。
- 3.
增强采样方法：采用先进的采样技术来克服能垒，更有效地探索不同三聚体组合（异源三聚体与同源三聚体）形成的可能性与相对稳定性。
- 4.
热力学积分法：这是一种基于“炼金术”思想的精确自由能计算方法。通过计算机模拟，逐步地将一个体系中的某些原子或残基“变成”另一个体系中的对应物，同时精确计算过程中的自由能变化。本研究利用此方法计算了proα2α(I)链特定残基在钙离子存在下，为异源三聚体形成所带来的相对结合自由能优势，从而从热力学角度量化了钙离子的驱动作用。
钙离子缺失导致链间界面去稳定化
通过分子动力学模拟，研究人员发现钙离子的缺失会显著影响胶原C-端肽链内和链间的氢键网络。这种氢键网络的破坏导致三条肽链的稳定性下降，尤其体现在链与链相互接触的界面区域。模拟结果显示，在缺乏钙离子的环境中，维持三聚体结构的关键相互作用被削弱。
钙离子缺失阻碍链间二硫键形成
二硫键是稳定蛋白质三级结构的重要共价键。本研究特别关注了由半胱氨酸(Cysteine)残基形成的链间二硫键。模拟分析表明，钙离子的丢失会导致形成这些二硫键的半胱氨酸残基之间的空间距离增大。距离的增加直接阻碍了二硫键的有效形成，从而从结构上阻止了稳定的异源三聚体的组装。
钙离子与α1(I)链的结合更强
为了探究钙离子对不同类型胶原链的偏好性，研究人员进行了拉伸模拟和钙离子解离建模。结果清晰地显示，钙离子与proα1α(I)链的结合力要强于与proα2α(I)链的结合力。这意味着在竞争环境下，proα1α(I)链对钙离子有更高的亲和力。
α2(I)链在钙离子存在下对三聚体形成具有更高亲和力
尽管钙离子更“喜欢”与proα1α(I)链结合，但利用增强采样方法进行的模拟却得出了一个有趣的结论：在结构钙离子存在的情况下，proα2α(I)链对形成三聚体的亲和力实际上高于第三条proα1α(I)链。这表明，虽然钙离子本身与proα1α(I)结合更牢，但它创造了一个更有利于proα2α(I)链参与并形成异源三聚体的微环境。
α2(I)链界面残基赋予异源三聚化热力学优势
研究最关键的部分是通过热力学积分法对三聚体组装的结合自由能进行了精确的量化计算。结果表明，在proα2α(I)链的链间界面处存在一些特异性残基，这些残基在钙离子存在的条件下，能够为异源三聚体的形成提供可测量的热力学优势。也就是说，从能量上看，形成由两条proα1α(I)链和一条proα2α(I)链组成的异源三聚体，比形成三条都是proα1α(I)链的同源三聚体更加有利（能量更低，更稳定）。而这种优势依赖于钙离子的存在。
综上所述，本研究通过多层次的分子模拟，清晰地描绘了钙离子驱动I型胶原三聚化选择的分子图景：在正常的钙离子浓度下，钙离子通过稳定链间界面、促进二硫键形成，并特别通过增强proα2α(I)链界面残基的贡献，为异源三聚体的组装提供了显著的热力学优势，从而确保了正常胶原的形成。然而，当钙离子可利用性降低时，这种优势消失。由于proα1α(I)链对钙离子具有更强的结合能力，它们可以“抢占”有限的钙离子资源，进而促使三条proα1α(I)链聚集形成同源三聚体。这项研究首次在分子动力学层面揭示了钙离子作为“分子开关”调控I型胶原异源与同源三聚化的详细机制，不仅深化了对胶原蛋白生物合成基础生物学过程的理解，更重要的是为解释与异常同源三聚体相关的多种年龄性疾病的发病机制提供了新的理论依据，并为未来开发针对这些疾病的干预策略指明了潜在靶点。该研究成果发表于《Biophysical Journal》。
