[論文レビュー] Weighted Ensemble Path Sampling for Multiple Reaction Channels

Bin W. Zhang, David Jasnow+1 more|ArXiv.org|Feb 16, 2009

Advanced Chemical Sensor Technologies被引用数 5

ひとこと要約

この論文は、遷移経路サンプリング（TPS）のような従来の手法とは異なり、特定の経路に閉じ込められることなく、分子系における複数の反応経路を効果的にサンプリングできる重み付きアンサンブル（WE）経路サンプリング法を示している。本手法はアラニンジペプチドを500Kおよび300Kで高品質な遷移経路アンサンブルを生成し、特に常温ではブルートフォースシミュレーションが非現実的であるにもかかわらず、経路の相対的重要性を定量的に分析可能である。

ABSTRACT

Finding and sampling multiple reaction channels for molecular transitions remains an important challenge in physical chemistry. Here we show that the weighted ensemble (WE) path sampling method can readily sample multiple channels. In a first test, both the WE and transition path sampling methods are applied to two-dimensional model potentials. The comparison explains why the weighted ensemble approach will not be trapped in one channel. The WE approach is then used to sample the full transition path ensemble in implicitly solvated alanine dipeptide at two different temperatures. The ensembles are of sufficient quality to permit quantification of the fractional importance of each channel, even at T=300K when brute-force simulation is prohibitively expensive.

研究の動機と目的

• 分子のコンformational遷移における複数の反応経路をサンプリングする課題に対処すること。
• 重み付きアンサンブル（WE）手法が、TPSや類似の経路サンプリング手法に共通する閉じ込め問題を回避することを示すこと。
• ブルートフォースシミュレーションが極めて高コストとなる常温（300K）における、統計的に厳密で高品質な遷移経路アンサンブルを生成すること。
• 遷移経路アンサンブル内における異なる反応経路の相対的寄与度を定量化すること。

提案手法

• WE手法は、ブルートフォース経路の無偏見な複製と結合を用いて、特定の経路に偏らない経路空間の探索を実現する。
• 重み付き経路の集団を維持し、重みを動的に調整するとともに、経路の分裂・統合を制御することで統計的正確性を確保する。
• WEとTPSの両手法が複数の経路をどのようにサンプリングするかを比較するため、二次元モデルポテンシャルを用いる。
• アラニンジペプチドに対しては、イオン化溶媒モデルを用いて500Kおよび300KにおけるC7eqとC7ax状態間の遷移をWE法でサンプリングする。
• ΦおよびΨ二面角の回転方向に基づいて、異なる経路を同定・分類できるように、拡張された二面角平面で経路アンサンブルを分析する。
• WEの結果をブルートフォースシミュレーションと比較することで、反応速度および経路分布の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 WE手法は、TPSとは異なり、特定の経路に閉じ込められることなく、複数の反応経路を効果的にサンプリングできるか？
• RQ2 アラニンジペプチドのような現実的なバイオ分子系において、WE手法は高温および常温での遷移経路のサンプリングにどの程度有効か？
• RQ3 300Kにおける遷移経路アンサンブルにおいて、各々の異なる反応経路が全体の経路アンサンブルに占める相対的寄与度は何か？
• RQ4 ブルートフォースシミュレーションが非現実的である状況下でも、WEは正確な反応速度と経路アンサンブルを同時に提供できるか？

主な発見

• 500Kでは、WE手法が1.4×10⁻¹/nsの反応速度を生成し、ブルートフォース結果の1.5×10⁻¹/nsと非常に良好に一致した。
• 300Kでは、WE手法が1.6×10⁻³/nsの反応速度を算出し、500Kと比較してエネルギー障壁が100倍に増加していることが示された。
• 300Kにおける経路アンサンブルは、ΦおよびΨの時計回り・反時計回り回転に対応する4つの主要な経路タイプを明らかにした。特に「右下」のC7ax状態が最も支配的であった。
• 「右下」経路は、Ψ=0およびΦ=0付近に高い障壁によって2つのサブ経路に分かれていることが判明した。
• 経路分布の分析から、低温になるほど支配的経路の相対的寄与度がさらに高くなることが示され、その動的優位性が裏付けられた。
• 本研究は、統計的に厳密なアラニンジペプチドの遷移経路アンサンブルを報告した初の例であり、経路の相対的寄与度の定量的評価を可能にした。