Comparative Assessment of Free Energy Computational Methods for Revealing the Interactions Driving PARP1 Selective Inhibition
J. Chem. Inf. Model. | 2026年4月 | DOI: 10.1021/acs.jcim.6c00083
PARP1/PARP2 選択性を MM/PBSA・ABFE・US の 3 手法で系統比較。ABFE が ΔΔG＜1 kcal/mol の微差まで再現し、接触連結性が選択性の物理的根拠と判明。
(1) 背景と課題：PARP1 選択性問題

PARP1 阻害薬は BRCA1/2 変異がんに対して合成致死を利用した有効な治療戦略だが、配列相同性の高い PARP2 を同時に阻害すると造血毒性が生じ治療域が狭まる。saruparib・NMS-P118 等の PARP1 選択的阻害薬 はこの課題を解決する候補だが、計算による事前評価が困難。

PARP1 と PARP2 の触媒ドメインはアミノ酸配列相同性が極めて高く、選択性 ΔΔG が 1 kcal/mol 以下 の微差にとどまるケースが多く、再現が最難関。
複数の自由エネルギー手法を 同一の 8 化合物セット で系統比較したベンチマークはこれまでなく、手法選定の客観的基準が無い。

→ 3 つの代表的自由エネルギー計算手法を統一プロトコルで比較し、選択性予測の信頼性と物理的駆動因子を定量的に明らかにする。

(2) 手法の概要：3 計算手法 + 接触解析
  • MM/PBSA：Amber mmpbsa.py によるエンドポイント評価。MD トラジェクトリのスナップショットから ΔG を分解。
  • ABFE：熱力学サイクルでリガンド消滅を 16〜20 λ ウィンドウで実行。TIP3P + ff14SB、長時間 MD。
  • US (Umbrella Sampling)：解離経路に沿うハーモニック拘束ウィンドウで PMF を WHAM 推定。
  • 接触解析：MD から残基-リガンド接触頻度を計算しグラフ化、接触連結性 を選択性指標として提案。
統一評価プロトコル 8 PARPi PARP1 / PARP2 MD (ff14SB+TIP3P) 接触頻度マップ MM/PBSA (endpoint) ABFE (16–20 λ) US + WHAM (PMF) → 実験 ΔΔG (PARP1 vs PARP2) と相関比較 接触連結性: G863 等のキー残基を同定
(3) 本研究で示したこと
  • 3 手法を 同一 8 PARPi で初めて系統比較し、定量精度・コストのトレードオフを明示
  • ABFE は ΔΔG＜1 kcal/mol の微差ケースでも選択性トレンドを再現可能
  • MM/PBSA は定性傾向は捉えるが定量精度に限界、エンドポイント感度が誤差源
  • US は ABFE と同程度の精度、解離 PMF からも一致した結論を得る
  • 「接触連結性」という新指標で選択性の 物理的駆動因子 を残基レベルで説明
  • PARP1 特有残基 G863 がキー接触として浮上、選択性設計の指針に
(4a) 手法別 ΔΔG 相関
実験 ΔΔG vs 計算 ΔΔG (8 PARPi) 実験 ΔΔG (kcal/mol) 計算 ΔΔG -1 0 1 2 3 y=x ABFE (高相関) US (高相関) MM/PBSA (散布大)
ABFE と US は実験 ΔΔG と直線関係に乗り、MM/PBSA は散布が大きい。8 PARPi すべてで一貫した傾向。
(4b) 微差 ΔΔG＜1 kcal/mol の符号再現
|ΔΔG| ≤ 1 kcal/mol ケースの符号一致率 100% 75% 50% 25% 0% ~50% MM/PBSA ~88% ABFE ~75% US
最難関の |ΔΔG|≤1 kcal/mol 帯で ABFE は選択性符号を最も高頻度で再現。MM/PBSA はほぼコイントス相当で実用性に欠ける。
(4c) 接触連結性グラフ：選択性のキー残基
残基-リガンド接触ネットワーク PARP1 PARP2 LIG Y889 G863 H862 S904 D766 5 接触 (G863 含む) LIG Y473 H447 S471 3 接触 (G 等価残基なし) → G863 が PARP1 選択性のキー 太線=H結合, 細線=van der Waals 接触
(4d) コスト vs 精度の手法選定マトリクス
手法定量精度微差再現計算コスト用途
MM/PBSA低〜中不十分低 (数 CPU 時間)初期スクリーニング
ABFE最高可能 (≤1 kcal/mol)大 (数 GPU 時間)リード選択性検証
US (PMF)高概ね可能大 (数百 CPU 時間)解離経路解析と併用
8 PARPi × 2 タンパク
同一プロトコルでの系統評価規模 (saruparib・NMS-P118 等を含む臨床/研究化合物)
ABFE/US は 1 リガンドあたり数百 CPU 時間〜数 GPU 時間で大規模 VS 適用は現状困難。タンパク質の多ドメイン柔軟性が大きい系では追加サンプリングが必須。
(5) テイクホームメッセージ
ABFE が選択性予測のゴールドスタンダード
ΔΔG＜1 kcal/mol の極小エネルギー差でも符号トレンドを再現。saruparib 級のリード化合物検証に直結する精度。
MM/PBSA は予測ではなく仮説生成ツール
スループットは高いが微差選択性の定量予測には不向き。フィルタリング・ポーズスクリーニングに用途を限定すべき。
G863 が選択性の構造的駆動因子
PARP1 特有のグリシン残基周辺の接触連結性が PARP2 より高く、ここを意識した修飾が選択性最適化の指針となる。
接触連結性が新しい選択性指標
単なるエネルギー値だけでなく接触ネットワークの構造を見ることで、選択性の物理的根拠を残基レベルで説明可能。
パイプライン応用
  • lib/fep の MMGBSAEngine 拡張: ABFE プロトコル (16〜20 λ, ff14SB+TIP3P) をプリセット化
  • 接触連結性スコアラを lib/md に新規追加: HBondAnalyzer + 接触頻度グラフから選択性指標を出力
  • PARP1/PARP2 のような パラログ選択性 評価テンプレートを lib/fep 用例に登録
  • WHAM ベース PMF 計算ユーティリティを US ワークフローとして整備
インパクト
  • PARP1 選択的阻害薬（合成致死療法）の 計算的設計指針 を提示
  • パラログ選択性問題に汎用的な手法選定基準を提供
  • 接触連結性という構造解釈ツールがリード最適化に活用可能
Comparative Assessment of Free Energy Computational Methods for PARP1 Selective Inhibition · J. Chem. Inf. Model. 2026 · DOI: 10.1021/acs.jcim.6c00083