[论文解读] Molecular dynamics, Langevin, and hybrid Monte Carlo simulations in multicanonical ensemble

Ulrich H. E. Hansmann, Yuko Okamoto+1 more|ArXiv.org|Oct 18, 1997

Protein Structure and Dynamics参考文献 1被引用 10

一句话总结

本文将多基准系系方法从传统的蒙特卡罗模拟扩展至分子动力学、朗之万动力学及混合蒙特卡罗算法。结果表明，仅通过一次模拟运行，即可高效采样甲硫氨酰脑啡肽的完整能量景观，从而计算任意温度下的正则系综性质与热力学量，包括全局能量最低构象。

ABSTRACT

We demonstrate that the multicanonical approach is not restricted to Monte Carlo simulations, but can also be applied to simulation techniques such as molecular dynamics, Langevin, and hybrid Monte Carlo algorithms. The effectiveness of the methods are tested with an energy function for the protein folding problem. Simulations in the multicanonical ensemble by the three methods are performed for a penta peptide, Met-enkephalin. For each algorithm, it is shown that from only one simulation run one can not only find the global-minimum-energy conformation but also obtain probability distributions in canonical ensemble at any temperature, which allows the calculation of any thermodynamic quantity as a function of temperature.

研究动机与目标

• 将多基准系系方法从标准蒙特卡罗扩展至动态模拟技术。
• 实现对肽类等生物分子系统中复杂能量景观的高效采样。
• 实现仅通过一次模拟运行即可计算任意温度下的正则系综性质与热力学量。
• 以甲硫氨酰脑啡肽为测试体系，验证该方法在蛋白质折叠问题上的有效性。

提出的方法

• 通过修改动力学以优先采样未充分采样的能量状态，将多基准系系框架适配至分子动力学。
• 在多基准系系框架内实现朗之万动力学，通过引入阻尼力与随机力，同时保持偏置势能。
• 通过将确定性动力学与随机更新相结合，将混合蒙特卡罗与多基准系系方法集成。
• 采用依赖于系统能量的偏置势能，以增强对罕见状态（尤其是低能量构象）的采样。
• 使用特定的蛋白质折叠势能函数，将该方法应用于五肽（甲硫氨酰脑啡肽）。
• 通过重加权方法，从多基准系系数据中提取正则系综分布与热力学性质。

实验结果

研究问题

• RQ1 多基准系系方法能否成功应用于分子动力学模拟？
• RQ2 朗之万动力学能否与多基准系系方法有效结合以增强采样？
• RQ3 在多基准系系系综中，混合蒙特卡罗模拟能否获得准确的热力学性质？
• RQ4 一次多基准系系系综模拟能否提供所有温度下的完整正则系综信息？
• RQ5 该方法在复杂生物分子系统中能否可靠地定位全局能量最低构象？

主要发现

• 多基准系系方法成功扩展至分子动力学、朗之万动力学及混合蒙特卡罗模拟。
• 一次多基准系系系综模拟即可实现任意温度下正则系综性质的计算。
• 使用三种动态方法均能可靠地定位甲硫氨酰脑啡肽的全局能量最低构象。
• 可从一次模拟中计算出内能与熵等热力学量随温度的变化关系。
• 该方法显著增强了对低能量状态的采样，包括蛋白质折叠中至关重要的罕见构象。
• 该方法可实现精确重加权，从而获得正则系综分布，而无需额外模拟。