GAMA：一种稳健且自动化的基于片段的量子化学方法用于生物分子系统

《The Journal of Physical Chemistry Letters》：GAMA: A Robust and Automated Fragment-Based Quantum Chemistry Method for Biomolecular Systems

　　

基于片段的量子化学方法处理共价键合生物分子仍是一项主要挑战，因为跨共价键的片段化会破坏必要的电子关联和长程极化。重要的是，先前开发的几种基于片段的方法中，少数能够准确且高效地处理非共价和共价键合分子系统，这凸显了该领域的一个显著空白。网格自适应多体分析（GAMA）框架的一个关键新颖之处在于克服了这一局限。在研究人员先前为非共价系统建立GAMA的工作基础上，研究人员将该框架扩展到共价键合生物分子，并开发了GAMA2，这是一种全自动协议，整合了简单的基于网格的片段化方案、在二体阶截断的具有重叠片段的多体展开以及多层低层校正。在各种肽（从柔性生物活性基序到结构化18聚体螺旋）中，当使用HF作为低层理论时，GAMA2对柔性中小型肽系统的全系统MP2/6-311G(d,p)能量的无符号绝对误差约为0.01–4 kcal/mol；当使用M06-2X/6-311G(d,p)作为低层方法时，对复杂螺旋型肽结构的误差约为2–5 kcal/mol，显示出使用精确DFT基方法相对于HF的显著改进。除了这些高精度结果外，GAMA2还展示了与参考全MP2计算相比，使用HF作为全系统低层方法时显著的计算加速，从而将GAMA2确立为一种可扩展、高效且可系统改进的途径，用于生物分子系统的关联量子化学计算。

**论文解读**

**研究背景**
基于片段的量子化学（FBQC）方法是扩展关联量子力学到远超出传统MP2或耦合簇计算能力范围系统的重要手段。这些方法通过将大分子分解为可管理的子系统，独立计算其能量，并利用多体展开重建总能量。尽管取得了显著进展，但可靠处理共价键合系统仍是一个长期挑战，因为跨化学共价键的片段化会破坏轨道离域、局域极化和关联效应，这些效应必须被准确恢复。目前，少数先前开发的FBQC方法能够高效且准确地同时处理非共价和共价键合分子系统，这凸显了该领域的显著空白。为解决这一问题，研究人员发展了网格自适应多体分析（GAMA）框架，并进一步扩展至共价键合生物分子系统。

**研究内容与意义**
本研究中，研究人员将GAMA框架推广至共价键合肽系统，开发了全自动协议GAMA2。该方法结合了基于网格的片段化方案、在二体阶截断的具有重叠片段的多体展开以及多层低层校正，能够在显著降低计算成本的同时提供MP2级别的能量精度。研究采用了一系列化学多样的肽系统进行基准测试，包括短生物活性肽（Segglu、Segtrp、Tuftsin）、中等长度甘氨酸低聚物（Gly12）、结构化肽1YJP以及三种丙氨酸18聚体构象异构体（β-折叠、α-螺旋和3₁₀-螺旋）。结果表明，GAMA2对柔性小中型肽系统的MP2/6-311G(d,p)能量误差约为0.01–4 kcal/mol（HF低层），对复杂螺旋型肽结构误差约为2–5 kcal/mol（DFT低层）。此外，GAMA2展现了显著的计算加速：α-螺旋-Ala₁₈的全MP2计算需约1178 CPU小时，而GAMA2仅需约77 CPU小时（约15倍减少）。该工作发表在《The Journal of Physical Chemistry Letters》上，意义在于建立了一种可扩展、高效且可系统改进的通道，用于生物分子系统的关联量子化学计算。

**关键技术方法**
研究人员采用了以下主要技术方法：
1. **网格基片段化方案**：将系统置于三维空间盒子中，划分为立方网格单元；通过切割C_α-C_羰基单键（避免切割极性肽键和N-C_α键）将肽骨架分段，并用标准氢链接原子饱和断裂价键，将每个“原子组”映射到网格单元，生成主要片段和自然重叠片段。
2. **多体展开（MOBE）**：在二体阶截断，利用重叠片段恢复电子离域和极化效应。
3. **多层低层校正**：采用类似ONIOM的两层校正，片段间相互作用在MP2/6-311G(d,p)水平计算（高层），全系统低层计算（HF/6-311G(d,p)或DFT方法）提供互补的长程贡献。

**研究结果**
**1. GAMA2精度评估**：使用2 ?盒尺寸和5 ?截止半径，以HF/6-311G(d,p)为低层方法，GAMA2对柔性小中型肽的误差为0.01–4 kcal/mol，对Gly12和β-折叠-Ala₁₈误差分别为3.87和4.0 kcal/mol。对于紧凑且氢键丰富的螺旋结构，α-螺旋-Ala₁₈误差为13.53 kcal/mol，3₁₀-螺旋-Ala₁₈误差为7.66 kcal/mol，表明HF低层方法难以捕捉长程协同极化。
**2. 低层方法的影响**：将HF替换为B3LYP或M06-2X可显著降低误差，尤其对α-螺旋-Ala₁₈（从13.53降至约2–5 kcal/mol）和3₁₀-螺旋-Ala₁₈，表明精确低层校正是实现稳健GAMA2能量的关键。
**3. 网格盒尺寸效应**：固定截止半径5 ?，增大网格盒从1 ?至3 ?一致提升精度，因为更大盒子能更有效纳入长程片段相互作用并减少多体截断误差。
**4. 截止半径效应**：固定2 ?盒尺寸，增大截止半径从5 ?至9 ?系统降低误差，在5 ?后改善趋于平缓，支持5 ?作为最优二体截断距离。
**5. 计算时间对比**：对α-螺旋-Ala₁₈，全MP2需约192小时（壁钟时间），GAMA2仅需约1.17小时（约164倍加速）；对Gly12也观察到类似趋势。CPU时间方面，GAMA MP2较全MP2显著减少（α-螺旋-Ala₁₈：1178 vs 77 CPU小时），而GAMA HF的CPU时间可能高于全HF，但由于高效并行化，壁钟时间仍更低。

**讨论与结论**
总结讨论部分指出，GAMA2的精度受网格盒尺寸、截止半径和低层方法选择的可系统调优，实现了对共价键合生物分子系统的可控且高效处理。研究结论翻译如下：总之，研究人员将GAMA框架推广至共价键合肽，实现了稳健、自动且准确的片段基量子化学。GAMA2协议采用简单网格基片段化和二体截断的重叠片段多体展开并配以多层校正，为多样生物分子结构提供了关联MP2级别能量，兼具显著效率。其性能关键在于能量随网格尺寸和截止半径的系统收敛以及低层方法的关键作用，其中DFT基方法（B3LYP和M06-2X）即使对较大螺旋肽系统也提供良好精度。特别地，GAMA2设计用于相对能量、构象偏好和定性至半定量能量趋势为主要关注的应用，而非亚kcal mol^?1热化学精度。对于柔性中等大小生物分子系统（如肽和蛋白质片段），1–4 kcal mol^?1的能量差异通常足以进行可靠的构象排序和结构分析。对于更紧凑且高度协同的系统，约2–5 kcal mol^?1的误差仍能进行有意义的相对稳定性和能量趋势评估，尤其当全MP2或更高级别计算在计算上不可行时。在此背景下，GAMA2提供了精度与计算效率的有利平衡，使得全系统常规方法无法触及的关联量子力学处理成为可能。这种自动化、精度与计算效率的结合使GAMA2成为中等大小生物分子系统关联波函数计算的实用且可系统改进的框架。更先进的静电嵌入变体（以点电荷嵌入替代全系统低层校正）将大幅降低计算成本，并为GAMA扩展到更大生物分子系统提供路径，这将在未来工作中探索。