[論文レビュー] 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

ひとこと要約

この論文は、伝統的なモンテカルロシミュレーションにとどまらず、分子動力学、ランジュバン動力学、ハイブリッドモンテカルロアルゴリズムへとマルチカノニカル集合の手法を拡張する。本手法により、1回のシミュレーション走行でメチオネイン・エンkeファリンの全エネルギー・ランドスケープを効率的にサンプリングでき、任意の温度におけるカノニカル集合の性質や熱力学的量を計算可能となる。

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.

研究の動機と目的

• 標準的なモンテカルロ手法にとどまらず、マルチカノニカル集合の手法を動的シミュレーション技法へ拡張すること。
• ペプチドのようなバイオ分子系における複雑なエネルギー・ランドスケープの効率的サンプリングを可能にすること。
• 1回のシミュレーション走行から、任意の温度におけるカノニカル集合の性質や熱力学的量を計算可能にすること。
• メチオネイン・エンケファリンをテスト系として用いて、タンパク質折りたたみ問題へのアプローチの妥当性を検証すること。

提案手法

• マルチカノニカル集合フレームワークを分子動力学に適応させ、未十分にサンプリングされたエネルギー状態を優遇するように動力学を修正する。
• 摩擦力と確率的力を取り入れつつバイアスポテンシャルを維持する形で、マルチカノニカルフレームワーク内にランジュバン動力学を実装する。
• 決定的動力学と確率的更新を組み合わせることで、ハイブリッドモンテカルロをマルチカノニカルアプローチに統合する。
• 系のエネルギーに依存するバイアスポテンシャルを用いて、とりわけ低エネルギー状態のサンプリングを強化する。
• 特定のタンパク質折りたたみエネルギー関数を用いて、ペプチド（メチオネイン・エンケファリン）に本手法を適用する。
• リウェートリングによりマルチカノニカルデータからカノニカル集合の分布や熱力学的性質を抽出する。

実験結果

リサーチクエスチョン

• RQ1 マルチカノニカル集合は分子動力学シミュレーションに効果的に適用可能か？
• RQ2 ランジュバン動力学はマルチカノニカルアプローチと効果的に組み合わせられるか？
• RQ3 マルチカノニカル集合におけるハイブリッドモンテカルロシミュレーションは、正確な熱力学的性質をもたらすか？
• RQ4 マルチカノニカル集合における1回のシミュレーション走行で、すべての温度における完全なカノニカル集合の情報が得られるか？
• RQ5 本手法は、複雑なバイオ分子系において、グローバルな最小エネルギー構造を信頼性高く特定できるか？

主な発見

• マルチカノニカルアプローチは、分子動力学、ランジュバン動力学、ハイブリッドモンテカルロシミュレーションへと効果的に拡張可能である。
• マルチカノニカル集合における1回のシミュレーション走行で、任意の温度におけるカノニカル集合の性質を計算可能である。
• すべての3つの動的手法を用いて、メチオネイン・エンケファリンのグローバルな最小エネルギー構造が信頼性高く特定された。
• 内部エネルギー、エントロピーなどの熱力学的量を、1回のシミュレーションから温度関数として計算可能である。
• 低エネルギー状態、とりわけタンパク質折りたたみに重要な希少な構造のサンプリングが強化された。
• 追加のシミュレーションを要せず、正確なリウェートリングによりカノニカル分布が得られる。