Computationally Driven Design of Small-Molecule Reversible Covalent Inhibitors for Targeted Therapeutics
熱力学×速度論×耐性を統合する可逆共有結合阻害剤のマルチスケール設計（批判的レビュー） ｜ MSc Essay (2026, Laurentian Univ.)
🎯 可逆共有結合阻害剤(RCI)を、ケムインフォ→QSAR/DL→ドッキング/MD→FEP/QM/MM→ADMET と縦に積むファネル型パイプラインで合理設計し、ΔG_bind・K_i・ΔG‡・K_cov で統一解釈する。
① 背景と課題

共有結合阻害剤は持続的標的占有を可能にするが、不可逆型は累積毒性リスクを抱える。RCI は過渡的な共有アダクトで持続占有と毒性回避を両立する。だが結合親和性・共有結合速度論・耐性進化が複雑に絡み合い合理設計が難しい。本稿は単一手法依存(siloed)を批判し、計算化学を縦に積むマルチスケールファネルとして統合する。

従来設計はドッキング/古典 MD で「姿勢が安定か」止まりで、共有結合形成という化学反応の実現可能性を問わない
姿勢の安定性は共有結合形成の必要条件だが十分条件ではない（空間的近接だけでは遷移状態を有利にしない）
求核標的がシステイン一辺倒、単一共有結合は monovalent vulnerability（耐性）を抱える

→ 選択性は共有結合状態の安定化で決まる／ΔG‡ は active-site 内 QM/MM で計算／二重ワーヘッド bidentate で耐性ハードニング

② 手法: マルチスケール・ファネル（FEP + QM/MM）
マルチスケール設計ファネル ケムインフォ/ML トリアージ 臨床標的・耐性マッピング QSAR / 深層学習 potency 予測 ドッキング / MD 姿勢・占有評価 FEP / QM/MM 精密生物物理 ΔG_bind・ΔG‡ in silico ADMET / off-target
二状態 FEP: 非共有(要拘束) vs 共有アダクト(拘束不要)
ΔG‡ は QM/MM で in situ 評価・Eyring K_inact=(kBT/h)exp(−ΔG‡/RT)

選択性は共有結合状態の安定化で決まる／二重ワーヘッド bidentate で耐性克服。批判的レビュー（独自データなし）。

③ 本研究で示したこと（要点）
  • 熱力学(ΔG_bind/K_i)・速度論(ΔG‡/K_cov)・耐性を統合するファネル型マルチスケール設計枠組み
  • 選択性は初期非共有認識でなく共有結合状態の安定化で決まるという認識転換（二状態 FEP）
  • ΔG‡ を QM/MM で active-site 内計算する設計パラメータとして提示（治療レンジ 15–22 kcal/mol）
  • 単一共有結合の脆弱性を二重ワーヘッド bidentate で克服する耐性ハードニング戦略
④ 主な結果 (a) SrtA リシン標的
17.0 W7(親) 4.6 P3(APBA) 1.3 P5(RMR1) SrtA リシン標的 IC50 (µM,↓良)

RMR1 ワーヘッドで ~13倍向上、k_on 6.3(遊離比~300倍)・滞留~34h・Lys173 単一アダクト(LC-MS)。

④ 主な結果 (b) カルパイン選択性
0.87 共有アダクト FEP カルパイン選択性 vs 実験 (Pearson r)

非共有状態のみでは弱い逆相関、共有アダクト状態を含めると r≈0.87。系統誤差が相殺。

④ 主な結果 (c) 耐性ハードニング
50× ZNL-0056(bidentate) 300× オシメルチニブ EGFR L858R/C797S potency シフト (×,↓良)

二重ワーヘッドが C797S 耐性を緩和（50× vs 300×の活性喪失）。

④ 主な結果 (d) 速度論・機構
ΔG‡ 治療レンジ 15–22 kcal/mol
~47+ kcal/mol で共有捕捉せず（蛋白質環境依存）。1 kcal/mol 低下→K_inact ~4倍
KRASG12C: 反応性が potency 支配
非共有が一定なまま共有反応性が2桁超変動（GSH 反応性と相関）
ERK2: MD 安定でも QM/MM 高障壁
機構誘導の再設計で低障壁経路を作り共有捕捉を回復
⑤ テイクホームメッセージ
🧬 共有状態で選択性
選択性は非共有認識でなく共有アダクト安定化で決まる（FEP r≈0.87）。
⚗️ ΔG‡ が設計変数
QM/MM で active-site 内に計算、治療レンジ 15–22 kcal/mol。
🔁 可逆共有
過渡的アダクトで持続占有と毒性回避を両立。
🛡️ 耐性ハードニング
二重ワーヘッド bidentate で C797S 耐性を緩和（50× vs 300×）。
4 ケーススタディ
系手法知見
カルパイン二状態 FEP共有状態で選択性 r≈0.87
KRASG12CQM/MM 速度論共有反応性が potency 支配
SrtA (Lys)ジアザボリンP5 IC50 1.3µM・滞留 34h
EGFR (bidentate)covalent docking耐性 50× vs 300×
本研究のインパクト
  • lib/fep: DockFEP/MMGBSAEngine に CovalentFEP ラッパを追加し、二状態（非共有レグ=flat-bottom 距離拘束 / 共有レグ=結合形成・カスタム結合パラメータ）を λ-windows で並列実行、熱力学サイクルで ΔΔG_cov を組み立てアイソフォーム対をランキング（系統誤差相殺で相対選択性を信頼）
  • lib/md: QM/MM 反応座標マッピングで共有結合形成の活性化障壁 ΔG‡ を active-site 内計算し Eyring 式で K_inact に変換。治療レンジ(15–22 kcal/mol)で候補をトリアージし、MD で姿勢安定でも QM/MM 高障壁の偽陽性を排除
  • lib/docking・lib/molgen: covalent docking + bidentate スペーサ長/ワーヘッド配置最適化を UniDockRunner/MolgenYaml に組込み、求核標的(Cys/Lys)と前組織化幾何で耐性ハードニング設計を支援
Shah, F. Computationally Driven Design of Small-Molecule Reversible Covalent Inhibitors for Targeted Therapeutics. MSc Critical Review Essay, Laurentian University. 2026.