[論文レビュー] Rejection-free kinetic Monte Carlo simulation of multivalent biomolecular interactions

Jin Yang, William S. Hlavacek|arXiv (Cornell University)|Dec 25, 2008

Gene Regulatory Network Analysis参考文献 35被引用数 3

ひとこと要約

本稿では、サイトおよび凝集体の状態に基づいて反応確率を直接サンプリングすることで、非効率な試行錯誤による拒否ステップを回避する、多価バイオ分子相互作用をシミュレートする拒否なしのキネティックモンテカルロ法を提示する。この手法は、全分子サイトではなく凝集体の数に比例してスケーリングするため、特に広範な凝集を示す系において、従来の拒否ベースの手法よりも速度とスケーラビリティに優れた優れた効率性を達成する。

ABSTRACT

The system-level dynamics of multivalent biomolecular interactions can be simulated using a rule-based kinetic Monte Carlo method in which a rejection sampling strategy is used to generate reaction events. This method becomes inefficient when simulating aggregation processes with large biomolecular complexes. Here, we present a rejection-free method for determining the kinetics of multivalent biomolecular interactions, and we apply the method to simulate simple models for ligand-receptor interactions. Simulation results show that performance of the rejection-free method is equal to or better than that of the rejection method over wide parameter ranges, and the rejection-free method is more efficient for simulating systems in which aggregation is extensive. The rejection-free method reported here should be useful for simulating a variety of systems in which multisite molecular interactions yield large molecular aggregates.

研究の動機と目的

• 多価リガンド-レセプター相互作用のキネティックモンテカルロシミュレーションにおける拒否サンプリングの非効率性、特に大規模凝集体が支配的なソル-ゲル領域において、これを解決すること。
• 試行サイトを拒否するのを避けるために、各候補サイトの正確な反応確率を計算する統計的に正確なシミュレーション手法の開発。
• 広範な凝集を示す系において、すべての個々の分子サイトではなく凝集体のみを追跡することで、計算コストを低減すること。
• 連続的解法との比較を通じて、さまざまな運動論的パラメータにおいて本手法の妥当性を検証すること。
• 多サイト・多価相互作用および大規模分子クラスターを含む複雑な生化学的系の効率的かつスケーラブルなシミュレーションを可能にすること。

提案手法

• ルールリストから反応ルールをサンプリングし、各サイトの状態と凝集体内での接続性に基づいて、そのサイトが反応に参加する正確な確率を計算する。
• サイトをランダムにサンプリングして不適合なものを拒否するのではなく、アルゴリズムは反応確率分布に従って直接サイトを選択することで、拒否イベントを完全に排除する。
• 動的な凝集体リストを維持し、各凝集体内のサイト数と接続性に基づいて反応速度を計算することで、全サイトからではなく、有効な凝集体のみに探索空間を制限する。
• 全凝集体の累積反応速度に基づいて反応イベントを選択することで、統計的正確性と詳細つり合いを保証する。
• 全ネットワークの列挙を避けるために、次の反応イベントを総速度に基づいて効率的に特定するため、優先度キューまたは類似のデータ構造を用いる。
• 循環的複合体における凝集体内結合の禁止などの非局所的適用条件を正しく計算するために、接続性情報を追跡する。

実験結果

リサーチクエスチョン

• RQ1 大規模凝集体が系を支配する場合、多価バイオ分子相互作用のキネティックモンテカルロシミュレーションは、どのようにしてより効率的に行えるか？
• RQ2 拒否なしのアプローチは、凝集プロセスにおいて高い拒否率が原因で性能が制限される従来の拒否サンプリングのボトルネックを解消できるか？
• RQ3 拒否なし手法のスケーリング特性は、分子サイトの数と凝集体の数のどちらに依存するか？
• RQ4 オフレートと結合Valenceが変化する条件下で、拒否なし手法と拒否ベース手法の性能はどのように比較されるか？
• RQ5 拒否なし手法は、低次元の多価結合モデルに対して連続的解を正確に再現できるか？

主な発見

• 拒否なし手法は、広範な運動論的パラメータの範囲で、拒否ベース手法と同等またはそれ以上の性能を示し、特に広範な凝集を示す系において顕著な優位性を示す。
• 本手法の計算コストは、全分子サイト数ではなく凝集体の数に比例するため、ソル-ゲル領域において顕著な効率的向上が達成される。
• 平衡状態では、平均凝集体数⟨NA⟩は受容体数NR=5000の四分の一未満であった。これは、凝集体リストがコンパクトで管理可能であることを示している。
• 拒否率θが1に近づく、つまり凝集傾向が強い系では、拒否なし手法が拒否手法を上回る。
• 反応速度と非局所的接続性条件（例：循環的複合体における凝集体内結合の禁止）を正しく考慮することで、統計的正確性を維持する。
• 大規模凝集体の繰り返しグラフ走査を回避することで、標準的手法の主要なボトル neck をさらに軽減できる。