[论文解读] Free-Energy Calculations in Protein Folding by Generalized-Ensemble Algorithms

Y. Sugita, Y. Okamoto|arXiv (Cornell University)|Feb 16, 2001

Protein Structure and Dynamics被引用 4

一句话总结

本文提出了两种先进的广义系综算法——回火交换多宏观算法（REMUCA）和多维回火交换伞形采样（REUS），以提升蛋白质折叠过程中的自由能计算。通过结合回火交换的高效性与多宏观方法的广泛采样能力，该方法能够对复杂能量景观进行充分探索，实现宽温度和反应坐标空间内的平坦能量分布，并获得精确的自由能分布。

ABSTRACT

We review uses of the generalized-ensemble algorithms for free-energy calculations in protein folding. Two of the well-known methods are multicanonical algorithm and replica-exchange method; the latter is also referred to as parallel tempering. We present a new generalized-ensemble algorithm that combines the merits of the two methods; it is referred to as the replica-exchange multicanonical algorithm. We also give a multidimensional extension of the replica-exchange method. Its realization as an umbrella sampling method, which we refer to as the replica-exchange umbrella sampling, is a powerful algorithm that can give free energy in wide reaction coordinate space.

研究动机与目标

• 解决传统蒙特卡洛和分子动力学模拟在蛋白质折叠中采样不足的问题。
• 通过在势能景观上实现随机行走，克服能量势垒陷阱。
• 通过将回火交换与多宏观方法结合，提升自由能计算的精度与效率。
• 实现在包括溶剂化体系在内的宽反应坐标空间内的系统性自由能剖面分析。
• 开发鲁棒的、自动化的多宏观采样权重因子确定方法，基于回火交换模拟。

提出的方法

• 提出回火交换多宏观算法（REMUCA），利用短时回火交换模拟通过WHAM方法确定最优多宏观权重因子。
• 应用加权直方图分析方法（WHAM）从回火交换数据中提取多宏观权重，避免迭代调优。
• 引入多维回火交换方法（MREM），实现温度与势能参数的交换。
• 通过将MREM扩展至包含偏置势参数交换，实现一种新型伞形采样变体——回火交换伞形采样（REUS）。
• 使用确定的多宏观权重进行长时间生产运行，以实现广泛且平坦的能量分布。
• 采用单直方图与多直方图重加权技术，仅通过一次模拟即可计算不同温度下的正则系综平均值。

实验结果

研究问题

• RQ1 回火交换能否在无需迭代调优的情况下高效确定多宏观权重因子？
• RQ2 回火交换与多宏观采样相结合是否能提升蛋白质折叠模拟中的采样效率？
• RQ3 多维回火交换方法在多大程度上将回火交换的应用范围从温度交换扩展至其他参数？
• RQ4 回火交换伞形采样方法在气相与溶剂化肽类中对复杂崎岖能量景观的采样效果如何？
• RQ5 溶剂化对蛋白质能量景观崎岖性的影响如何？自由能分布揭示了何种结果？

主要发现

• 回火交换多宏观算法（REMUCA）成功地仅通过一次短时回火交换模拟即确定了多宏观权重，消除了对迭代权重优化的需求。
• 在水溶液中的Met-Enkephalin模拟中，实现了跨越2,500 kcal/mol能量范围的平坦概率分布，证实了对整个能量景观的广泛采样。
• 通过REMUCA与传统多宏观方法计算的势能正则系综平均值完全一致，验证了新方法的准确性。
• 与气相相比，水溶液中Met-Enkephalin的能量景观显著更平滑，300 K时二面角分布更宽且更弥散。
• 回火交换伞形采样（REUS）通过结合参数交换与增强采样，实现了在宽反应坐标空间内的高效自由能计算。
• 仅通过一次使用优化权重的模拟运行，该方法成功获得了从200 K到1000 K温度范围内的可靠自由能剖面。