Exploring the Capabilities of Machine-Learned Potentials for Biomolecular Simulations (NNP/MM)
Galvelis & De Fabritiis (Acellera) | arXiv:2201.08085 | 2022
🎯 ANI-2x をリガンド内ポテンシャル、AMBER/GAFF をその他相互作用に使う NNP/MM ハイブリッドを OpenMM/ACEMD で実装し、複合体長時間 MD を実用化。
① 背景と課題

古典 MM 力場のリガンド内ポテンシャルは単純フーリエ展開で量子化学的エネルギー面と乖離し、配座サンプリングにバイアスが入る。一方 QM/MM はサブナノ秒サンプリングしかできない。

GAFF/CGenFF の二面角ポテンシャルは QM 面と乖離、配座バイアスの原因
DFT/MM は 1 ステップ計算が膨大で長時間 MD が困難

→ ANI-2x をリガンドの内分子項に限定的に使う NNP/MM で QM 精度と MM スループットを両立。

② 手法の概要: NNP/MM 分割

E_total = E_NNP(リガンド内) + E_MM(タンパク内 + リガンド-タンパク + リガンド-溶媒)

  • ANI-2x: H, C, N, O, S, F, Cl の 7 元素対応
  • DFT/CCSD(T) 訓練データから内分子項を学習
  • その他相互作用は AMBER/GAFF 古典力場
② 手法の概要: 実装と高速化
  • OpenMM 上の Custom TorchForce
  • NNPOps GPU 最適化 CUDA fused カーネル
  • TorchANI 比 5.9× 推論高速化
  • 古典 MD の 5〜10× コストで長時間 MD 可能
  • ACEMD MD エンジンとも統合
③ 本研究で示したこと
  • 4 タンパク質-リガンド複合体で各 ~250 ns NNP/MM 軌道
  • 骨格 RMSD 1.8〜2.8 Å の安定軌道
  • リガンド配座サンプリングで古典 MD と差異を観測
  • 実用的なコストで QM 精度のリガンド動力学を実現
④ 主な結果 (a) NNPOps 推論高速化
推論時間（相対） 5.9× TorchANI 1.0× NNPOps CUDA fused kernel で大幅高速化
④ 主な結果 (b) 検証4系
RMSD ranges (Å) per system 3 1.5 0 1AJV 1HPO 2P95 3BE9 Lig 1.7
④ 主な結果 (c) 既存手法比較
手法リガンド精度速度
GAFF低1× (高速)
QM/MM★高~1000× 遅い
NNP/MMQM相当5-10×
④ 主な結果 (d) 系の特徴
~250 ns 各複合体での NNP/MM 軌道時間
4 系 1AJV (Xa), 1HPO (HIV-1Pr), 2P95 (CDK2), 3BE9 (thrombin)
5-10× 古典 MD に対するコスト倍率
⑤ テイクホームメッセージ
NNP/MM はバランスの良い選択
QM 精度（リガンド内）+ MM スループット。
NNPOps が実用化の鍵
5.9× 推論高速化で長時間 MD が現実に。
配座サンプリングに差
NNP/MM は古典 MD より広い空間を探索。
ANI-2x の元素制限
金属・ハロゲン (Br/I) は対象外、適用前に要確認。
lib/md + lib/fep 応用

RMSDAnalyzer/HBondAnalyzer のバックエンドに OpenMM TorchForce + NNPOps を採用。

  • NNP/MM ヘルパ
  • FEP のハイブリッドハミルトニアン
  • 軌道再現性チェック
本研究のインパクト
  • タンパク質-リガンド複合体長時間 MD への NNP/MM 系統的適用の先駆け
  • NNPOps 共有による NNP-MD コミュニティの加速
  • 創薬 MD ワークフローへの ML 力場の標準的組込み準備
Galvelis, R.; De Fabritiis, G. Exploring the Capabilities of Machine-Learned Potentials for Biomolecular Simulations. arXiv 2022, arXiv:2201.08085. https://arxiv.org/abs/2201.08085