SILCS-Kinetics: High-Throughput Ligand Dissociation Kinetics Predictions
Yu, Kumar, MacKerell Jr. et al. — J. Chem. Theory Comput. 2025, 21, 4964–4978 | DOI: 10.1021/acs.jctc.5c00265
🎯 物理ベースGFEプロファイル+MLでリガンドkoffを高スループット予測し、創薬での居留時間設計を加速する
① 背景と課題

薬物の解離速度(koff)と residence time は結合親和性(Kd)より有効性と相関する系が多く、kinetics-driven drug designへの関心が高まっている。CCR2拮抗薬など複数の臨床事例でkoffの設計が成否を分けた。

強化サンプリングMD法(メタダイナミクス・WE・MSM)は原子詳細情報が豊富だが計算コストが膨大で、多数のリガンドを早期創薬段階でスクリーニングできない
分子力学記述子ベースのML法はタンパク質ごとに専用モデルが必要で汎化性に乏しく、物理的解釈も困難

→ 物理シミュレーション(SILCS)で自由エネルギープロファイルを算出し、それをMLへの入力として使う「ハイブリッド戦略」で高速・汎用・解釈可能なkoff予測を実現

② 手法: SILCS-Kineticsワークフロー
GCMC/MD(10並列)→ GFE FragMap生成

SILCS-Pathway(A*探索)→ 解離経路クラスタ

円柱スラブ配置(DTW平均経路)

SILCS-MC slab docking → LGFE/IEプロファイル

ML予測(RF / RNN+MLP)→ -log(koff)

SILCSシミュレーション(GPU数日)完了後は1分子あたり数分でkoff予測。

② 手法: MLアーキテクチャ

決定木系(RF/GB/XGB/LGBM): スカラー5特徴(ΣLGFE障壁・ΣIE障壁・LGFESMS・MW・NRot)を入力。RFが最良。

RNN+MLP(ハイブリッド): Bi-LSTM 2層でLGFE/IEプロファイル(可変長→pad40)を処理し、スカラー3特徴をMLPで処理→連結→-log(koff)予測。過学習が少なく外部タンパク質で安定。

5特徴
LGFE障壁・IE障壁・LGFESMS・MW・NRot
③ 本研究で示したこと(要点)
  • 物理ベースGFEプロファイルがkoffの汎用ML特徴量として機能することを実証
  • 13タンパク質329リガンドの汎用モデルが未知5タンパク質94リガンドでも有効
  • 原子単位のGFE分析でどの官能基がkoffに寄与するかを定量化
  • E4+P2(LGFE+IE+LGFESMS、最低障壁経路)が最も一貫した相関を示す
  • コンセンサス予測(RF+RNN+MLP)でさらに精度向上
④ 主要結果 (a) MLモデル性能比較
モデルTest RMSETest PR
RF(最良)0.760.76
GB0.830.71
LGBM0.840.70
XGB0.900.65
RNN+MLP0.890.65

訓練/テスト=190/45化合物(8タンパク質、BindingDB)

④ 主要結果 (b) 外部タンパク質での検証
PR 0.8 0.6 0.4 0.2 RIP1 Cka TrkA FAK sEH RF RNN+MLP
④ 主要結果 (c) エネルギー障壁相関
Energy Barrier vs -log(koff) koff Energy Barrier (kcal/mol) E1(LGFE障壁のみ) E4+P2(最良)

E4(ΔLGFE+ΔIE+LGFESMS)×P2(最低障壁経路)が8タンパク質で最も一貫した相関

④ 主要結果 (d) 処理速度比較
計算時間(1リガンドあたり・対数スケール) メタダイナミクス μs〜日 WE / MSM 時間〜日 SILCS(初期シミュ) GPU数日 SILCS-Kinetics(予測) 数分 ※SILCSシミュ完了後の1分子あたり予測時間
⑤ テイクホームメッセージ
高スループット化を達成
SILCSシミュ完了後は1分子数分でkoff予測。MDベース手法のμs計算を代替し、多数リガンドの kinetics評価が可能。
🧩 物理+MLハイブリッドが汎化性の鍵
タンパク質固有モデル不要。13タンパク質訓練モデルが未知5タンパク質でも機能する汎用性を実証。
🔬 原子解像度の設計情報
GFEスコアを官能基単位で分解し「どの基がkoffに寄与するか」を定量化。リード最適化の方向性を直接指示。
📐 E4+P2が最良の組み合わせ
LGFE障壁+IE障壁+LGFESMS(表面メタステーブル状態)を最低エネルギー経路で評価することが重要。
ケムインフォマティクスへの応用
適用先ユースケース
lib/mdRMSDAnalyzer拡張: koff/居留時間予測モジュール追加
lib/mdHBondAnalyzer: 結合ポーズ安定性→koff相関分析
lib/dockingUniDockRunner後段でRFベースkinetics scorer追加
lib/dockingProLIFCalculator: 解離経路上のH-bond消失プロファイル抽出

MWとNRotはRDKit Descriptorsで即座に取得可能。エネルギー障壁はVina scoreで代理実装可能。

本研究のインパクト
  • 創薬初期段階でresidence timeを最適化するワークフローを初めて実現(高スループット)
  • koff設計の計算コストをMDベース手法比で劇的に削減
  • 物理ベースの解釈可能性が従来のブラックボックスMLより医薬化学者に受け入れられやすい
13タンパク質 329リガンド
開発・検証データセット規模(外部+94リガンド)