[論文レビュー] Identification and Analysis of Transition and Metastable Markov States

Linda Martini, Adam Kells+3 more|arXiv (Cornell University)|May 13, 2016

Protein Structure and DynamicsBiochemistry, Genetics and Molecular Biology参考文献 40被引用数 18

ひとこと要約

本論文は、分子動力学（MD）トラジェクトリから遷移状態および準安定状態を自動的に同定するための新規なアルゴリズムを提案する。マコフ状態モデル（MSM）アプローチを用い、集合座標のセットから内在する最も遅い緩和レートを捉えることで、小規模ペプチド（Ala5）およびより大きなタンパク質（EGFR）において自由エネルギー障壁と動的経路を正確に特定する。本手法は、複雑なバイオ分子系における希な構造的遷移を研究するための堅牢で一般化可能なツールを提供する。

ABSTRACT

We present a new method that enables the identification and analysis of both transition and metastable conformational states from atomistic or coarse-grained molecular dynamics (MD) trajectories. Our algorithm is presented and studied by using both analytical and actual examples from MD simulations of the helix-forming peptide Ala5, and of a larger system, the epidermal growth factor receptor (EGFR) protein. In all cases, our method identifies automatically the corresponding transition states and metastable conformations in an optimal way, with the input of a set of relevant coordinates, by capturing accurately the intrinsic slowest relaxation rate. Our approach provides a general and easy to implement analysis method that provides unique insight into the molecular mechanism and the rare but crucial rate limiting conformational pathways occurring in complex dynamical systems such as molecular trajectories.

研究の動機と目的

• 分子動力学トラジェクトリにおける遷移状態および準安定状態を、反応座標の事前知識なしに一般化可能で自動化された手法で同定すること。
• 複雑なバイオ分子系において、希だが動的に重要な構造的遷移経路を同定する課題に取り組むこと。
• 高次元のトラジェクトリデータから自由エネルギー障壁および最も遅い緩和モードを正確に同定することを改善すること。
• 原子的スケールおよび粗粒度シミュレーションに適用可能な、計算的に効率的で実装可能なフレームワークを提供すること。

提案手法

• 本手法は、トラジェクトリデータからの構造的ダイナミクスを分析するため、マコフ状態モデル（MSM）フレームワークを採用する。
• 動的挙動および最も遅い緩和モードに基づいて、構成をクラスタリングすることで準安定状態を同定する。
• アルゴリズムは、状態空間を定義し、遷移確率を計算するために、集合座標のセットを入力として使用する。
• 遷移状態は、主要な反応座標上に位置し、フラックスの蓄積が最小限である構成を特定することで検出する。
• 遷移行列の固有値特性を活用して、遅い緩和レートを抽出し、これが制限段階に対応することを確認する。
• 本手法は、解析的モデルおよび実際のAla5およびEGFRのMDシミュレーションを用いて検証され、堅牢性と正確性が示された。

実験結果

リサーチクエスチョン

• RQ1 反応座標の事前知識なしに、分子動力学トラジェクトリにおける遷移状態および準安定状態を信頼性高く同定する方法は何か？
• RQ2 複雑なバイオ分子系において、最も遅い緩和モードが制限段階となる構造的遷移経路をどのように定義するか？
• RQ3 本手法は、既存の手法と比較して、自由エネルギー障壁および動的中間体を同定する上でどのように優れているか？
• RQ4 Ala5 や EGFR などの高い構造的複雑性を示す系において、本手法は遷移状態を正確に同定できるか？
• RQ5 本手法の性能は、集合座標の選択にどの程度依存するか？

主な発見

• 本手法は、既知のヘリックス-コイル遷移挙動と一致するように、Ala5ペプチドにおける遷移状態および準安定状態を正確に同定した。
• EGFRタンパク質においては、活性化経路に明確な遷移状態が検出され、実験的および先行シミュレーション結果と整合的であった。
• 同定された最も遅い緩和レートは、内在する動的障壁とよく一致しており、本手法が制限段階のプロセスを正確に捉えていることを確認した。
• 異なる集合座標の選択に対しても本手法は高い頑健性を示し、一般化可能性を裏付けた。
• アルゴリズムは手動によるバイアスなしに構造的状態を自動的に同定し、状態割り当てにおける主観性を低減した。
• 本手法は明確な動的経路と well-defined な中間状態を明らかにし、EGFR活性化機構におけるメカニスティックな知見を提供した。