PLQM-VM2: QM-Based Mining Minima for Protein-Ligand Binding Free Energies

Computation of Protein−Ligand Binding Free Energies with a QM-Based Mining Minima Algorithm (PLQM-VM2)

Schlinsog, Sattasathuchana, Xu, Guidez, Gilson, Potter, Gordon, Webb — J. Chem. Theory Comput. 2025 (DOI: 10.1021/acs.jctc.4c01707)

🎯 MM-VM2コンフォーマーをDFTB3-D3/PCMでQM補正し、多様3系でMM比較を大幅改善。マルチプロテインスクリーニングへの展望を示す

① 背景と課題

タンパク質-リガンド結合自由エネルギー計算において、MM エンドポイント法（MM-PBSA/MM-GBSA/MM-VM2）は古典力場の限界（π-スタッキング・ハロゲン相互作用・分極・電荷移動の記述不全）から精度が一貫しない。FEP/TI は厳密だが計算コストが高く、MD の時間スケール制約もある。

MM コンフォーマー探索（Mining Minima / VM2）
↓ QM ポストプロセス（DFTB3-D3(BJ)H/PCM via GAMESS）
↓ QM 補正配置積分 → 化学ポテンシャル → 結合 ΔG
→ PLQM-VM2

② 2 種の QM モデル

カットアウトモデル: 結合ポケット周辺原子を切り出しH キャップ → 縮小系に QM 直接適用

30–45 min

シングルポイント (32 CPU コア 1 ノード)

FMO (Fragment Molecular Orbital): 全タンパク質を断片分割して QM 計算 → 最高精度

~35 h

幾何最適化込み (32 CPU コア 1 ノード)

③ Mining Minima の特徴

MM-VM2 の強力な探索能力＋ QM の精度 = PLQM-VM2

④ HIV-1 プロテアーゼ（38 リガンド）

全 4 手法で相関が大幅改善。FMO + 幾何最適化が最高精度。カットアウト単点でも実用的精度。

④ TNKS2（27 リガンド、最難関）

MM: 相関≈0

電荷混在リガンドで MM-VM2 が完全に失敗

PLQM-VM2 全手法で中程度の相関に改善。電荷-中性リガンド間の不一致を QM が解消。

さらなる改善には高次 QM 理論 + 大基底関数が必要と示唆

④ 系統誤差（DFTB3 の過大結合）

MSE -24〜-29 kcal/mol

全手法共通の系統的過大結合誤差（DFTB3 由来）

MSE オフセット補正後の非系統的誤差（RMSE）は実用的レベル。

← 重要: 相対ランク付けには大きな影響なし。絶対値の直接利用は要注意。

④ マルチプロテインスクリーニングへの展望

QM 補正により多様系間での絶対自由エネルギースケールが一貫。

系間比較が可能に

MM-VM2 では不可能だった多様タンパク質間の直接比較

オフターゲット活性評価（セレクティビティスクリーニング）への応用が期待される

⑤ テイクホームメッセージ

⚗️ QM 補正で一貫精度向上
古典 MM-VM2 を超える相関改善を 3 多様系で実証。π-スタッキング・電荷相互作用の記述が改善。

⏱️ カットアウト単点が実用的
30-45 min で FMO (35h) に匹敵する精度。ルーチン利用への現実的な入口。

🔍 系間スケール一貫性
多様ターゲット間の絶対 ΔG 比較が可能に。マルチプロテインオフターゲット評価への展望。

⚠️ DFTB3 系統誤差に注意
MSE -24〜-29 kcal/mol の過大結合。相対ランキングには実用的だが、絶対値利用は補正必須。

ケムインフォmaticへの応用

適用先	提案
lib/fep	MDAnalysis カットアウト + xTB QM ポストプロセス実装
lib/fep	MMGBSAEngine に QM 補正フックを追加
lib/docking	QM スコアをセレクティビティフィルタとして統合

xTB/DFTB+ は MIT 近傍ライセンスで利用可能。VM2 不要の軽量代替スタックとして活用できる

限界点

VM2 は VeraChem LLC 商用ソフト（要ライセンス）

DFTB3 系統的過大結合 (-24〜-29 kcal/mol) はオフセット補正が必要

プロトン化状態変化を陽に扱えていない

TNKS2 系の電荷混在には高次 QM が必要