DynAgent / DynaMate: A Modular Multi-Agent Framework for Autonomous Protein-Ligand Molecular Dynamics Simulations
自然言語からタンパク質-リガンド MD を自律設計・実行・解析する 3 エージェント枠組み ｜ arXiv:2512.10034（EurIPS 2025 Workshop, EPFL）
🎯 計画(PrepAgent)・実行(DynAgent)・解析(Analyzer)の 3 LLM エージェントが、自然言語指示から PDB 取得〜パラメータ化〜MD〜解析〜MM-PB(GB)SA までを自己修正しながら無人実行する。
① 背景と課題

MD は創薬・タンパク質工学の基盤だが、系のパラメータ化・入力準備・ソフト設定がエラー多発で専門知識の障壁が高い。DynAgent はドメイン知識を持つ LLM エージェントにツール群を操作させ、外部 retrieval（RCSB/PubChem/PaperQA）と Reflexion 型自己修正で障壁を除去。高レベル推論と低レベル実行を意図的に分離する。

CHAPERONg/PyAutoFEP は補助的だがパイプラインが硬直的・領域特化で、拡張に手動修正と専門知識が必要
MDCrow(OpenMM)/NAMD-Agent(CHARMM-GUI) はタンパク質単独のみでリガンド系を扱えない
先行 LLM 系はシミュレーションエラーから適応的回復できず、プラットフォーム間で汎化しない

→ リガンド系対応＋外部DB/文献 retrieval＋適応的ツール選択・自己修正を統合し、AmberTools と GROMACS 双方に対応

② 手法: 3 エージェント（計画・実行・解析）
3 エージェント・ワークフロー PrepAgent（計画） NL要件抽出・PDB取得・RAGパラメータ選定 DynAgent（実行） ツール呼出→検証→内省 + 自己修正 Analyzer（解析） RMSD/RMSF/Rg/H-bond + 自動テキスト option: MM-PB(GB)SA 結合自由エネルギー連鎖
呼出 → 検証 → 内省（Reflexion 型）
ff14SB/GAFF2/TIP3P・AmberTools+GROMACS 双方対応・LiteLLM でモデル抽象化

web_search + PaperQA で系固有設定を接地（例: ヒト由来 310K / 中温性酵素 298K）。GitHub: schwallergroup/DynaMate。

③ 本研究で示したこと（要点）
  • リガンドを含む系を扱える自律 MD マルチエージェント（先行 MDCrow/NAMD-Agent はタンパク質単独）
  • 外部 DB(RCSB/PubChem)＋文献(PaperQA) retrieval でパラメータをエビデンス接地
  • 「呼出→検証→内省」ループの Reflexion 型自己修正で頻出ランタイムエラーを自動修復
  • AmberTools と GROMACS のクロスプラットフォーム対応＋MM-PB(GB)SA 連鎖
④ 主な結果 (a) 系別精度
100 3PTB 100 5UEZ 100 4GIH 100 3HTB 80 4W52 リガンド系 精度（必須タスク完了率,%）

PL 系 5/5 で本計算 MD 完走（4 系 100%、4W52 は自己修正で 80% 完走）。リガンド無 3 系は欠損残基で失敗。

④ 主な結果 (b) BRD4 親和性
-23.86 cmpd17 (IC50 2436nM) -26.06 cmpd25 (IC50 49.5nM) BRD4 BD1 計算ΔΔG (kcal/mol)

自律 10ns MD + MM-PBSA の計算ΔΔGが実験活性序列と整合し強力阻害剤 25 を判別。

④ 主な結果 (c) 成分分解
34.1 33.2 VdW 16.6 10.8 ΔEEL cmpd25 cmpd17 結合エネルギー成分 (|kcal/mol|)

成分分解で静電(ΔEEL)差が結合差の主因（25: −16.57 vs 17: −10.75）。

④ 主な結果 (d) バックエンド
Opus 4.1 / Sonnet 4.5
高効率・PL 系 4 つで精度 100%（最良バックエンド）
Haiku 3.5
ツール順序選択を誤りやすく効率・精度が劣化
8 系で検証
5 protein-ligand（3PTB/5UEZ/4GIH/3HTB/4W52）+ 3 ligand-free
⑤ テイクホームメッセージ
🤖 自律MD
自然言語から PDB 取得〜MD〜解析まで 3 エージェントが無人実行。
🔁 自己修正
Reflexion 型「呼出→検証→内省」で原子名欠落・トポロジ誤りを自動修復。
🔗 親和性連鎖
ドッキング→短時間MD→MM-PB(GB)SA で結合序列を判別（BRD4 で実証）。
⚠️ 前処理課題
欠損残基/未知原子名の PDB では下流ツールが失敗（リガンド無 3 系）。
自律MD枠組み比較
枠組みリガンド系エンジン自己修正
MDCrow×OpenMM×
NAMD-Agent×CHARMM-GUI×
CHAPERONg等一部固定的×
DynAgent（本研究）○Amber+GROMACS○
本研究のインパクト
  • lib/md: RMSDAnalyzer/HBondAnalyzer をツールスキーマ化し、計画エージェントに PDB取得→クリーニング→パラメータ化→溶媒和→平衡化→本計算→解析を無人実行させ、Reflexion 自己修正でセットアップ失敗率を削減
  • lib/fep: MMGBSAEngine を gmx_MMPBSA エージェントツールに登録し、UniDockRunner(lib/docking)→短時間 MD→MM-PB(GB)SA を自律連鎖。計算ΔΔG が活性序列を判別する初期スクリーナに
  • lib/md・lib/docking: RAG パラメータ推薦で温度・プロトン化・計算時間を RCSB/PubChem/文献に接地し自動決定、専門家以外でも妥当設定で実行
DynAgent / DynaMate: A Modular Multi-Agent Framework for Autonomous Protein-Ligand Molecular Dynamics Simulations. arXiv:2512.10034 (EurIPS 2025 Workshop). 2025. https://github.com/schwallergroup/DynaMate