Quantum chemical modelling of sequence encoded energetics and electrostatic landscapes governing early peptide association in α- synuclein and amyloid β

技术摘要：序列编码的早期肽类结合中的能量与静电景观的量子化学建模

问题陈述
蛋白质错误折叠和聚集是帕金森病（PD）和阿尔茨海默病（AD）等神经退行性疾病的核心机制，这些疾病由 $\alpha$-突触核蛋白（$\alpha$-synuclein）和淀粉样 $\beta$（A$\beta$）等肽类的早期分子结合所驱动。虽然实验研究已经表征了成熟纤维的动力学和结构，但对于在瞬态、异质性的早期结合阶段中，单个分子构型对能量的具体贡献仍难以解析。目前存在一个关键的研究空白，即如何理解氨基酸的总组成之外，特定的残基线性顺序（序列）是如何塑造决定这些初始结合事件的分子属性的。本研究旨在通过将易聚集肽段与组成匹配的打乱对照组进行比较，从而分离出序列排列的贡献，进而确定残基顺序如何影响能量、结构组织和静电景观。

方法论
本研究采用了一套一致的计算工作流，整合了半经验量子化学和密度泛函理论（DFT），用于模拟单体和同源二聚体状态。

• 体系： 将易聚集肽段（$\alpha$-突触核蛋白：VTGVTAVAQK；A$\beta$：LVFFAEDVGS）与组成匹配的打乱对照组（AGATQTVVVK 和 EFLVGSDVAF）进行了对比。
• 几何优化： 使用 GFN2-xTB 方法（第二代几何、频率、非共价扩展紧束缚方法）在 xTB 中对单体和同源二聚体结构进行了优化。二聚体结构在最终优化前先使用 GFN 力场进行了弛豫。
• 能量计算： 结合能（$\Delta E_{assoc}$）通过优化后的二聚体的收敛 GFN2-xTB 电子能量与两个独立优化的单体能量之差来计算。
• 电子结构分析： 在优化的单体上使用 GPAW 进行 DFT 计算，采用 PBE 交换相关泛函和实空间网格（间距 0.20 Å）上的投影缀加平面波（PAW）形式。
• 描述符：
  • 静电学： 将分子静电势（ESP）映射到电子密度表面（0.001 e Bohr$^{-3}$），并通过均值中心化进行分析，以可视化相对景观。
  • 电荷分布： 使用 Hirshfeld 粒子分析来量化原子级电荷分布，包括统计指标（均值、中位数、标准差）和残基级总和。
  • 结构/序列： 描述符包括序列电荷装饰（SCD）、局部亲水性、疏水序列长度、回转半径（$R_g$）、端对端距离和不对称度。
  • 界面： 通过重原子距离（$\le$ 4.5 Å 和 5.0 Å）以及特定的 N/O/S 接触（$\le$ 3.5 Å）来量化分子间接触。

主要结果
研究表明，尽管氨基酸组成保持不变，残基重排会对能量和结构领域产生截然不同的、肽特异性的响应。

1. 单体能量： 残基顺序显著改变了孤立单体的电子能量。对于 A$\beta$，打乱对照组的能量比易聚集序列低约 45.33 kJ/mol。相反，对于 $\alpha$-突触核蛋白，易聚集序列比其打乱对照组低约 151.12 kJ/mol。
2. 结合能量与界面堆积：
   • A$\beta$： 易聚集序列（115.617 kJ/mol）和打乱对照组（232.552 kJ/mol）均表现出正的结合能，表明优化的二聚体在电子能量上高于独立优化的单体。然而，易聚集序列表现出较低的正结合能（低 116.935 kJ/mol），但其形成的重原子接触较少（40 对 81），且接触密度低于对照组。
   • $\alpha$-突触核蛋白： 打乱对照组表现出更具利于结合的（更负的）结合能（-494.371 kJ/mol），而易聚集序列的结合能为 -172.819 kJ/mol，差异为 321.552 kJ/mol。然而，易聚集序列形成了更广泛的界面，具有 166 个重原子接触（对比 105 个）以及更高的近距离 N/O/S 接触密度（15 对 5）。
   • 结论： 结合能和界面堆积（接触数/密度）并不呈线性变化；它们是结合过程的互补而非等效度量。
3. 静电景观 (ESP)： 残基重排重新组织了 ESP 的空间分布。虽然 A$\beta$ 和 $\alpha$-突触核蛋白对具有相似的全局电荷状态，但易聚集序列与打乱序列之间，正电势和负电势区域的空间组织发生了显著偏移。这表明序列顺序决定了电荷特征的三维呈现方式。
4. Hirshfeld 电荷分布： 原子电荷分析显示，残基重排主要是重新分配了原子电荷群体，而非统一改变全局电荷量级。在 $\alpha$-突触核蛋白中，正/负原子的比例保持几乎不变，但电荷值的概率分布发生了移动。在 A$\beta$ 中，观察到了正/负原子平衡的轻微偏移以及电荷的重新分布。
5. 序列与结构描述符：
   • A$\beta$： 重排改变了序列电荷装饰（SCD）、疏水序列长度（5 个残基对比 3 个）以及分子紧凑度，易聚集肽的 $R_g$ 约为 10.1 Å，而打乱对照组约为 11.0 Å。
   • $\alpha$-突触核蛋白： 重排主要改变了局部亲水性（对照组的最大 3 残基平均亲水性从 ~2.6 增加到 ~4.2），而 SCD、$R_g$ 和不对称度保持相当。

意义与主张
本文确立了残基顺序是早期肽类结合的基本决定因素，且独立于氨基酸组成。研究证明，相同的组成可以通过序列依赖的分子组织产生截然不同的能量、结构和静电特性。

核心主张包括：

• 序列作为决定因素： 残基的排列（而非仅仅是组成）支配着分子间相互作用的平衡，从而影响结合的热力学过程和产生的界面几何结构。
• 肽特异性响应： 残基重排的分子后果并非普适性的；$\alpha$-突触核蛋白和 A$\beta$ 表现出截然不同的响应，即序列顺序的变化会导致能量有利性和界面堆积的不同组合。
• 互补描述符： 研究强调结合能和界面接触指标捕捉到了分子结合的不同方面。一个更密集的界面并不一定意味着更具利的结合能，正如在两种不同的肽系统中所观察到的发散趋势。
• 静电重组： 由于电荷模式和三维结构的重组，相同的氨基酸组成可以产生独特的静电环境，这为序列特异性的结合行为提供了分子基础。

这项工作提供了一个识别序列编码分子决定因素的定量框架，表明早期聚集的“代码”在于残基的具体线性排列，这种排列决定了最终形成的能量和静电景观。