自由エネルギー計算によるキノーム全体の選択性最適化 — Wee1 阻害剤ケーススタディ
Knight, Gerasyuto et al. (Schrödinger) — Nature Communications 16:7962 (2025) | DOI: 10.1038/s41467-025-62722-w
🎯 タンパク質残基変異FEP(PRM-FEP+)でゲートキーパー5変異だけを計算し、数百キナーゼへの結合傾向を一括予測。Wee1選択的阻害剤を効率的に発見
① 背景と課題

キナーゼはATP結合ポケットが高度に保存されており、目的キナーゼだけを選択的に阻害するのは難しい。リガンド相対結合自由エネルギー計算(L-RB-FEP+)は単一標的の効力を実験精度(平均1.0 kcal/mol = 6〜8倍以内)で予測できるが、数十〜数百のオフターゲットを個別にモデル化する選択性評価は計算量的に破綻する。

QSAR・Free-Wilson・MLモデルはスパースな訓練データに依存し、未踏の化学空間(コア変更等)へ外挿できない
オフターゲットごとに高分解能構造とFEPを用意するのは数百キナーゼでは非現実的

→ 稀少なWee1ゲートキーパー Asn376 を「選択性ハンドル」とし、残基側を摂動するPRM-FEP+でキノーム選択性を一括予測

② 手法概要(PRM-FEP+)
Wee1 WT (Asn376 ゲートキーパー)

残基376をアルケミカル変異
Asn → Thr / Phe / Met / Leu / Val

bound(複合体)−unbound(apo) の ΔG

ΔpKi (ゲートキーパーファミリー別結合傾向)
5回のFEP
1リガンドあたり5ゲートキーパー変異で全キノーム選択性を推定

OPLS4 / Desmond GPU の FEP/REST。12 λ窓・10 ns(トリアージ2 ns、荷電/大変異は36窓・25 ns)。変異残基376をREST hot regionに包含。ΔpKiがより負ほどWee1選択的、5変異すべてpKi<8で "clean" 判定。

③ 計算化学パイプライン連携(lib/fep)
  • lib/fep: 新規 PRMFEPEngine。残基変異アルケミカル摂動でΔpKiを算出(MMGBSAEngine拡張)
  • lib/fep: L-RB-FEP+ の on-target→off-target カスケードフィルタ(DockFEP拡張)
  • lib/docking: 物理化学フィルタ+ドッキング+MLで設計を事前絞り込み(UniDockRunner)
  • lib/molgen: AutoDesigner型の大規模列挙→FEP優先付け(JobManager/MolgenYaml)
実装ギャップ: lib/fep にタンパク質残基変異FEPエンジン・選択性ハンドル自動提案・ゲートキーパー5変異の選択性スコアラーが未実装
④ 主な結果 (a) レトロスペクティブ検証精度
91% 予測精度 1 偽陽性 2 偽陰性 0 100% 7スキャフォールド×5変異 vs scanMAX

PRM-FEP+ 予測がscanMAX 403キナーゼのヒットパターンを91%精度(FP 1・FN 2)で再現。promiscuous/ファミリー指向/選択的を正しく区別。

④ 主な結果 (b) 6,6-core 最適化の選択性ヒット率
82% Cpd4(出発) Thr GKヒット 0% Cpd11 22 nM 3 Cpd12 1 nM, 3キナーゼ 0 82% promiscuous → 選択的への変遷

出発点 compound 4 はThrゲートキーパーの82%をヒット。PRM-FEP+誘導でcompound 11(22 nM・サブパネル全20キナーゼ100倍超選択)、compound 12(1 nM・scanMAXで3キナーゼのみ)へ。

④ 主な結果 (c) 前向きキャンペーン成果
段階規模 / 成果
初期ヒット同定(L-RB-FEP+)9000設計を評価, 7ヶ月で80化合物合成
PLK1選択性(早期)最大 1000倍(vs PLK1)
選択性最適化(PRM-FEP+)6000+化合物を計算, 3ヶ月
合成・目標達成42合成中 22化合物が目標達成
新規スキャフォールド5,6 / 6,6 / トリシクリックの3系統

22化合物が単〜二桁nMのWee1効力と AZD1775同等以上の選択性を達成。compound 14は1 nM・scanMAXで14キナーゼのみヒット。

⑤ テイクホームメッセージ
  • 残基変異FEP(PRM-FEP+)をキノーム選択性予測に初適用
  • 選択性ハンドル1残基の5変異だけで広いオフターゲットを一括評価
  • 水素結合の有無等の定性SARを超え、脱溶媒和・エントロピー・補償効果まで捉える
  • キナーゼ以外(プロテアーゼ・GPCR・核内受容体)への拡張も提案
限界: ハンドルがゲートキーパー以外/単一点変異で表せない場合は適用限定。不確かさ約1 kcal/mol。選択性ハンドル自動同定は未統合。最終判定は実測アッセイ併用が必要。