FE-ToolKit: A Versatile Software Suite for High-Dimensional Free Energy Surfaces and Alchemical Free Energy Networks
Giese, Snyder, Piskulich, Barletta, Zhang, McCarthy, Ekesan, York | J. Chem. Inf. Model. 2025, 65, 5273–5279 | DOI: 10.1021/acs.jcim.5c00554
🎯 ndfes（多次元FES解析）・EdgeMBAR（FEPネットワーク解析）・fetkutils（λスケジュール最適化）の3コンポーネントで自由エネルギー計算の解析を自動化・統合する。
① 背景と課題

FEPシミュレーションが大規模化するにつれ（数千ウィンドウのFES、数百エッジのFEPネットワーク）、手動解析が限界に達した。サイクルクロージャーエラー・未平衡サンプリング・不良な位相空間オーバーラップの自動検出・補正が不可欠となっている。

従来MBARnetはメモリ大・計算コスト高、新データ追加時に全体再最適化が必要
最小長サイクル以外の大サイクルのクロージャーが保証されない問題が残存

→ EdgeMBAR法でエッジ事前計算による効率的ネットワーク解、実験値拘束にも対応。

② FE-ToolKit の3コンポーネント
コンポーネント機能
edgembarEdgeMBAR法でFEPネットワーク全体解析・サイクルクロージャー補正
ndfesMBAR/vFEP/wTP/gwTP でN次元FES構築 + SASM最小自由エネルギーパス
fetkutilsλスケジュール最適化 + 動力学同位体効果（KIE）計算

MIT License, AMBERTools同梱、独自ファイル形式でMD非依存

③ EdgeMBAR法の要点
  • 各エッジを独立MBAR解析 → unconstrained ΔΔG*とforce constant k取得
  • quadratic近似でエッジ目的関数を事前計算（F̃(ab)(x)）
  • グラフ全体: F(c) = Σ F̃(ab)(cb-ca) を最小化
  • 線形方程式 M·c = KᵀXᵀg で N_lig-1 の相対自由エネルギーを一括解
  • 実験値をLinear constraintsとして任意に組み込み可能
④ 主な結果 (a) EdgeMBAR vs MBARnet 比較
EdgeMBAR vs MBARnet 比較 特性 MBARnet EdgeMBAR メモリ効率 低（全状態同時） 高（エッジ分離） 新データ追加 全体再計算 差分のみ 大サイクル閉包 保証なし 保証あり 実験値拘束 非対応 Linear constraint グラフレポート 限定的 HTML自動生成 🎯 EdgeMBAR = 効率的・確実なFEPネットワーク解析の決定版 差分更新・全閉包保証・実験値統合で産業応用レベルに対応
④ 主な結果 (b) ndfes: FES解析フロー
ndfes 多次元FES解析フロー ① メタファイル: ウィンドウ設定 + サンプルデータ ② MBAR / vFEP / wTP / gwTP でN次元FES構築 ③ SASM: 全サンプル集積でパス最適化（反復） ④ XML出力: ビン毎の自由エネルギー・標準誤差 ⑤ 未平衡自動検出 (CheckEquil) + トラブルシュート指標
④ 主な結果 (c) fetkutils λスケジュール最適化
λスケジュール最適化の概念 λ (0 → 1) dU/dλ 均一間隔（非最適） 最適化後（均一オーバーラップ） 急峻領域に多くのλ点を配置してオーバーラップを均一化
④ 主な結果 (d) 自動品質管理機能
自動 未平衡サンプリング検出・除外（ndfes-CheckEquil）
HTML エッジレポート・グラフレポートを自動生成
Lagrange マルチプライアー・dU/dλ・位相空間オーバーラップ指標
QM/MM + ML ポテンシャル（QDπ, GNN）との連携対応
⑤ テイクホームメッセージ
EdgeMBAR = FEPネットワーク標準解析へ
差分更新・全サイクル閉包保証・実験値拘束の三拍子。産業規模FEPキャンペーンに対応。
多次元FES解析を自動化
数千ウィンドウのアンブレラサンプリングをndfes+SASMで一括処理。QM/MM対応も完備。
λスケジュール最適化で精度向上
バーンイン後に本番λ値を自動決定。HRE効率改善とMBAR信頼性向上を同時達成。
MIT License / AmberTools統合
GitLabで全ソース公開。既存AMBERワークフロー（ProFESSA）とシームレスに連携。
lib/fep への適用ポイント

DockFEP のFEPネットワーク解析を EdgeMBAR に置き換え。fetkutils でλ最適化を前処理に組み込む。

  • edgembar-WriteGraphHtml.py でHTML結果レポート自動化
  • 実験IC50値をΔΔG拘束として組み込み精度向上
  • fetkutils λ最適化を DockFEP 実行前処理に統合
本研究のインパクト
  • FEPネットワーク解析の標準ツールとして実運用実績（AMBER ProFESSA統合）
  • EdgeMBAR で大規模ネットワークのサイクルクロージャー問題を解決
  • 実験値（IC50等）拘束付き自由エネルギー予測が可能に
Giese T.J. et al. FE-ToolKit: A Versatile Software Suite for Analysis of High-Dimensional Free Energy Surfaces and Alchemical Free Energy Networks. J. Chem. Inf. Model. 2025, 65, 5273–5279. https://doi.org/10.1021/acs.jcim.5c00554 | Code: https://gitlab.com/RutgersLBSR/fe-toolkit (MIT)