ParametrizANI: Fast, Accurate and Free Parametrization for Small Molecules
Arantes, Sinha & Palermo — ChemRxiv 2025 · DOI: 10.26434/chemrxiv-2025-m26n9 · UC Riverside · github.com/palermolab/ParametrizANI
🎯 ANI NNP で QM を代替し GAFF/OpenFF 二面角パラメータ化を6分・無料・DFT精度で実現
① 背景と課題 — 二面角パラメータ化のボトルネック

MD シミュレーションの精度は力場の品質に直結し、特に小分子の二面角パラメータ化が鍵となる。従来は DFT/CCSD(T) スキャン計算が必須で商用ソフト（Gaussian等）とHPC資源を要する。これが中小規模の研究室や創薬プロジェクトでの障壁となっていた。

QM 計算：数時間〜数日/分子・商用ライセンス必須・HPC 必須
既存 GAFF パラメータ：汎用であり新規スキャフォールドでの精度が不安

→ ANI NNP（TorchANI/AIMNet2/MACE-OFF/GFN2-xTB）を参照に置き換えることで「6分・無料・Google Colab」を実現

② 手法: 7ステップワークフロー
7ステップ・パラメータ化フロー ① SMILES → RDKit コンフォマー ② NNP で参照 E 計算 (4 種) ③ OpenMM で ELE+VdW のみ最小化 ④ ΔE を Fourier 級数フィット (n=1,2,3,4,6 / 線形最小二乗) ⑤ 新パラメータで OpenMM 再最小化 ⑥ 参照 NNP と比較・検証 ⑦ GAFF/OpenFF topology 出力
② 手法: NNP モデル選択ガイド
モデル元素精度速度
ANI-2xH,C,N,O,F,Cl,SCCSD(T)/CBS 相当◎
ANI-1ccxH,C,N,OCCSD(T)/CBS 相当◎
AIMNet2-wB97M有機全般DFT 相当○
MACE-OFF生体有機DFT 相当○
GFN2-xTB広範半経験的◎◎

F/Cl/S 含有リガンドは ANI-2x または AIMNet2 を選択

③ 主な利点
~6分
1二面角/CPU Colab
無料
Google Colab 利用
DFT精度
ANI-2x参照
4種
NNP 参照モデル
✓ HPC / 商用ソフト不要
✓ GAFF・OpenFF 両対応（4 ノートブック）
✓ OSS 公開（github.com/palermolab/ParametrizANI）
ANI-2x 対応元素 (7 種) と参照モデル数 H C N O F Cl S F/Cl/S 含有時は ANI-2x または AIMNet2 を選択 参照 NNP モデル 4 種 提供ノートブック 4 本 ワークフロー工程 7 ステップ
④ 主な結果 (a) 速度比較（定性）
計算時間比較（1 二面角・対数スケール） 1分 10分 1時間 10時間 100時間 ~6分 ParametrizANI (ANI-2x) 30〜60分 B3LYP/DFT (HPC) 10〜100時間 CCSD(T) (HPC)
④ 主な結果 (b) ANI vs 従来手法 精度

Rezaee et al. (2024) による小有機分子の配座エネルギー比較：

手法OPLS比B3LYP比
ANI-2x✅ 優位✅ 優位
ANI-1ccx✅ 優位✅ 優位
B3LYP/DFT○ 良好—
OPLS FF—△ 劣る

ANI モデルは vdW・H 結合等の分子内非共有相互作用を適切に捉える

二面角フィットに用いる Fourier 級数の次数 (n=1,2,3,4,6) n=1 n=2 n=3 n=4 n=6 線形最小二乗で各次数の係数を ΔE にフィット → 二面角パラメータ決定 二面角 0° 360°
⑤ lib/md への統合提案
  • SMILES → GAFF/OpenFF topology 自動生成：lib/md の前処理ステップとして ParametrizANI ワークフローを組み込み、新規リガンドの力場パラメータ自動生成を実現
  • ANI-2x エネルギー評価モジュール：lib/md の構造最適化・コンフォマー評価に TorchANI を統合し DFT 品質の高速評価を実現
  • HBondAnalyzer/RMSDAnalyzer との連携：ANI パラメータ生成→GROMACS/OpenMM MD→既存の HBond/RMSD 解析という統合パイプラインを構築
依存: torchani · openmm · rdkit · ambertools · openff-toolkit
⑥ 限界
ChemRxiv プレプリント（未査読）— 多化合物クラスでの系統的ベンチマークは今後の課題
金属含有リガンド非対応（ANI-2x は金属元素をカバーしない）
Google Colab の実行時間制限（12h/無料）が多二面角分子で問題になる可能性
  • MCPB.py との連携（メタロプロテイン）は将来課題
  • 論文内の定量ベンチマーク（RMSE 等の数値）が不十分
⑦ 総括
  • QM 不要・無料・6分/二面角 — 二面角パラメータ化の民主化
  • GAFF/OpenFF 両対応 × 4種 NNP 参照（ANI-2x/AIMNet2/MACE-OFF/GFN2-xTB）
  • DFT/CCSD(T) 精度を Google Colab CPU で達成
  • lib/md の SMILES → topology 前処理統合の最有力候補

OSS: github.com/palermolab/ParametrizANI

ParametrizANI · ChemRxiv 2025 · DOI: 10.26434/chemrxiv-2025-m26n9 · category: computational_chemistry · 生成日: 2026-05-11