[論文レビュー] Optimized parallel tempering simulations of proteins

Simon Trebst, Matthias Troyer+1 more|JuSER (Forschungszentrum Jülich)|Feb 6, 2006

Advanced Data Storage Technologies被引用数 16

ひとこと要約

本稿では、エネルギー障壁に相当する地点に温度レプリカを動的に集中させることで、タンパク質の並列温度法シミュレーションの収束を最適化する適応的フィードバックアルゴリズムを提案する。36アミノ酸からなるビリリンヘッドピece HP-36に適用した結果、従来の高性能分子動力学研究に比べて計算コストの1%で、実験的PDB構造と3.7 Å RMSD以内の最低エネルギー構造を同定した。

ABSTRACT

We apply a recently developed adaptive algorithm that systematically improves the efficiency of parallel tempering or replica exchange methods in the numerical simulation of small proteins. Feedback iterations allow us to identify an optimal set of temperatures/replicas which are found to concentrate at the bottlenecks of the simulations. A measure of convergence for the equilibration of the parallel tempering algorithm is discussed. We test our algorithm by simulating the 36-residue villin headpiece sub-domain HP-36 wherewe find a lowest-energy configuration with a root-mean-square-deviation of less than 4 Angstroem to the experimentally determined structure.

研究の動機と目的

• 粗いエネルギー・ランドスケープと低温での遅い収束性に起因するタンパク質折りたたみシミュレーションにおける非効率的サンプリングの課題に対処する。
• 標準的な並列温度法では、最適でない温度設定がコンformation空間の効率的探索を妨げるという制限を克服する。
• 収束性と収束の改善を図るため、温度レプリカ分布を体系的に最適化するフィードバック駆動型アルゴリズムを開発する。
• サンプリング効率の向上と収束基準の妥当性を検証することで、真のネイティブ構造または低エネルギー構造を信頼性を持って同定することを保証する。
• 従来のHP-36のシミュレーションが、温度設定が不適切であったために長時間実行されていても熱的平衡に達していなかった可能性を示す。

提案手法

• 温度空間におけるレプリカの拡散をモニタリングすることで、エネルギー・ランドスケープのボトルネックを特定する適応的フィードバックループを実装する。
• 観察された温度空間における拡散を基に、反復的にレプリカ温度の分布を調整し、緩和が遅い領域に集中させる。
• 最低温度への統計的に独立した訪問回数に基づく収束基準を用い、ターゲット温度での収束を保証する。
• メトロポリス受容確率を用いるレプリカ交換（並列温度法）フレームワークを採用：$ p = \min(1, \exp(-\Delta\beta \Delta E)) $、ここで $ \Delta\beta = 1/T_{i+1} - 1/T_i $ で、$ \Delta E $ はエネルギー差である。
• フィードバックを繰り返し適用することで温度セットを最適化し、受容率の向上と緩和時間の短縮を図る。
• シミュレーション全体を通して最低エネルギー構造を追跡し、実験的PDB構造（1vii）とRMSDを用いて比較する。

実験結果

リサーチクエスチョン

• RQ1 適応的フィードバックアルゴリズムは、タンパク質折りたたみにおける並列温度法シミュレーションの効率性と収束性を向上させることができるか？
• RQ2 複雑なエネルギー・ランドスケープにおいて、緩和時間を最小化し、サンプリング効率を最大化する最適なレプリカ温度分布は何か？
• RQ3 従来のHP-36の並列温度法シミュレーションは、熱的平衡に達していたのだろうか、それとも最適でない温度セットによって制限されていたのか？
• RQ4 計算コストを著しく削減した状況でも、ネイティブ構造とRMSD < 4 Å以内の低エネルギー構造を信頼性を持って同定できるか？
• RQ5 現在の力場とイオン化溶媒モデルは、サンプリングの向上にもかかわらず、依然としてシミュレートされたタンパク質構造の正確性を制限している程度はどの程度か？

主な発見

• 最適化された並列温度法アルゴリズムでは、約10回の往復時間後に最低温度での完全収束が達成され、収束には最低100,000 MCスイープが必要である。
• 2回のフィードバック反復後、半径の gyration が11 Å未満の構造の99%がサンプリングされた。これは、コン pact で低エネルギー状態の有効なサンプリングを示している。
• 最適化されたシミュレーションで得られた最低エネルギー構造は、実験的PDB構造（1vii）と3.7 Å RMSDを示し、従来の結果（5.8–6.0 Å）よりも顕著に低い。
• 本手法ではネイティブに類似した構造がグローバルエネルギー最小点として同定されたが、従来のシミュレーションでは長時間実行されていても、誤折りたたみ状態と区別できなかった。
• 計算リソースは、類似したRMSD値を報告した大規模分子動力学シミュレーションの約1,000 CPU年分の1%（約10 CPU年）で、最良の構造的一致が達成された。
• RMSDが低いにもかかわらず、ネイティブに類似した構造とPDB構造とのエネルギー差はわずか約10 kcal/molであり、現在のエネルギー関数に依然として限界があることが示唆された。