AIQM-PBSA — MLIP-Enhanced MMPBSA Binding Free Energy

AIQM-PBSA: Integrating ML Interatomic Potentials with MMPBSA

Wei, Chen, Yang, Xu, Dral & Chen — J. Phys. Chem. B 130, 4982-4993 (2026) | DOI: 10.1021/acs.jpcb.6c00306

🎯 古典力場の代わりにΔ-learning型MLIP「AIQM3」をONIOMでMMPBSAへ統合し、QMレベル精度を高速のまま実現する結合自由エネルギー計算

① 背景と課題

MMPBSAはFEPより圧倒的に高速な end-point 結合自由エネルギー（BFE）計算法だが、その信頼性は古典力場のポテンシャルエネルギー曲面（PES）と統計熱力学的近似に本質的に縛られている。古典力場は分極・電荷移動・複雑な分子内相互作用を捉えられず、内部誘電率εintで経験的に補償するしかない。

古典MMPBSAはεintに極端に敏感：データセットAでR=−0.40(εint=1.0)→0.55(εint=5.0)と激変

純MLIP（ANI-2x, MACE）は学習分布外データで精度が落ちやすい

→ SQMベースラインにNN補正を加えるΔ-learning型 AIQM3 をONIOM(MLIP:MM)でガス相エネルギー項に導入

② AIQM-PBSA 手法

結晶構造 → AMBER24で系調製
↓
10 ns MD（1000スナップショット）
↓
ONIOM二層: E=E_AIQM(QM) + E_MM(全系) − E_MM(QM)
QM領域=リガンド+6Å、境界Hキャップ
↓
PBSA溶媒和（PB極性 + SPT非極性 SASA+SAV）
↓
全フレーム平均 → ΔG_bind

AIQM3

GFN2-xTB* + 8 ANI NN補正 + D3(BJ)。CCSD(T)/CBSへ二段fine-tuning。H,C,N,O,F,S,Cl対応

③ 計算化学パイプライン連携 (lib/fep)

lib/fep: MMGBSAEngineのガス相エネルギー計算をMLIP/QMバックエンドに差し替え可能に抽象化
lib/fep: ONIOM(MLIP:MM)二層スキーム + QM領域自動切り出し（リガンド+6Å, Hキャップ）
lib/fep: SPT非極性溶媒和（SASA+SAV併用）とPB極性項
lib/fep: IE/QHA/NMA/FE エントロピー推定（IE推奨）とεintスキャン

実装ギャップ: MLIPアダプタ(MLatom/TorchANI/MACE)・QM領域切り出し・SPT項・IEエントロピーがlib/fepに未実装

④ 主な結果 (a) JACS set 相関 (εint=2.1)

AIQM-PBSAはMMPBSAより+0.31、ANI-2xより+0.16高いPearson R。Spearman/Kendallも0.58/0.42で最良。

④ 主な結果 (b) εint への感度

古典MMPBSAはεintでR=−0.40→0.55と激変。MLIPは分極・電荷移動を本質的に捉えεintにほぼ鈍感。

④ 主な結果 (c) 精度・効率の総括

手法	Pearson R (JACS)	Type
MMPBSA	0.27	MM
PBSA_E	0.47	MM
MACE-OFF23(S)	0.53	MLIP
MACE-OFF24(M)	0.56	MLIP:MM
AIQM-PBSA	0.59	MLIP:MM

+2〜10%

64コア+GPUの実用パイプラインでの総計算時間増（古典MM置換比）

データセットA/Bで Pearson R = 0.84 / 0.80。荷電リガンド60系も係数再フィットで PBSA_E に匹敵。

④ 主な結果 (d) 物理的起源とエントロピー

外れ値「4GIH_ejm31」: C-Cl···O ハロゲン結合のσ-holeを古典力場は等方球で誤り、AIQM3は第一原理的に正しく記述
エントロピーはIE推奨（1.6秒/系で全4ターゲット安定改善）。NMAはCDK2で0.73→0.60と劣化・1.8時間/系で高コスト
AIQM3はデータ記憶でなく物理を学ぶため未知ターゲットへの高汎化性

限界: 長距離相互作用とimplicit solvation近似が残課題。荷電系は係数再フィット必須。エントロピー推定は依然系依存。