TACK: PROTAC分解活性の統計的ML評価と新規ベンチマークデータセット
Ribes, Dunlop & Mercado（Chalmers）— KDD AI for Sciences Track '26 | arXiv:2605.19579
🎯 3,514 PROTACsの標準化データセット構築と厳密な統計的ML評価 — XGBoostが特化GNNを凌駕し、校正済みUQで実験優先度付けを支援
① 背景と課題

PROTAC（Proteolysis-Targeting Chimera）は標的タンパク質をユビキチン-プロテアソーム系で触媒的に分解する新しい創薬モダリティ。従来の阻害剤では困難な「undruggable」標的に対応できるが、ウォーヘッド・リンカー・E3リガンドの組み合わせ空間が膨大で合理的設計が困難。

公開データセット: 3ソース（PROTAC-DB, PROTACpedia, TPDdb）が異なる単位・閾値・アノテーションで管理されており横断比較が不可能
既存モデルはランダムスプリットを採用 → スキャフォールド類似化合物が訓練・テストに混在し過楽観的な性能評価
統計的有意性の欠如: 多重比較補正なしで手法比較 → p-hackingリスク

→ 標準化データセット + 厳密な統計評価プロトコル + 再現可能なMLベンチマークの三位一体が必要

② TACKデータセット統計
3,514
ユニークPROTAC化合物
6,561
分解エンドポイント
164
標的タンパク質 (POI)
155
細胞株

主要POI: AR(26.7%), SMARCA2(12.0%), BTK(8.9%) | 主要E3: CRBN+VHL = 98.7% | Active率: 55.1%（DC50≤100nM & Dmax≥80%）

ホールドアウトセット(10%)はTanimoto距離最大化で構造的多様性確保
② キュレーションパイプライン
PROTAC-DB 3.0 (6,110) + PROTACpedia (1,189) + TPDdb (21,429)
↓ SMILES正規化（RDKit）
↓ POI/E3 → UniProt ID統一
↓ 細胞株 → Cellosaurus検証
↓ DC50 nM統一・比較演算子別保存
↓ 有効エンドポイントフィルタ
→ TACK: 3,514 PROTACs / 6,561 endpoints
④ 主要結果 (a) 予測性能比較（ROC-AUC）

活性二値分類タスク（Tukey HSD検定, p<0.05で全ペア有意差あり）

0.70 0.75 0.80 0.85 0.90 PROTAC-STAN MLP XGBoost 0.741 0.799 0.851 ★ ★ p<0.001（全ペア）Tukey HSD Δ(XGB–STAN) = +0.11 | Δ(MLP–STAN) = +0.06
④ 主要結果 (b) pDC50 vs Dmax 予測困難性の非対称性

ホールドアウトセットでのR²比較（低いほど予測困難）

0.0 0.5 1.0 0.66 0.36 R²=0.66 pDC₅₀ （効力） R²=0.36 Dmax （最大分解率） DC50 ← 結合親和性支配 | Dmax ← 細胞文脈依存
④ 主要結果 (c) アンサンブルUQ校正

std（不確実性）vs |絶対誤差|のSpearman相関

不確実性（std）→ |予測誤差|→ pDC50 ρ=0.207 (p<0.01) Dmax ρ=0.543 (p<0.001)
④ 主要結果 (d) 特徴セット選定結果
タスク最良モデルCellProtein分子Best Val.
pDC50MLP†OneHotOneHotDescRMSE 0.770
pDC50XGBTextESM-S-PCADescRMSE 0.781
DmaxXGB†TextOneHotDesc+FPRMSE 18.34
活性XGB‡OneHotESM-S-PCADesc+FPAUC 0.851
活性MLP‡OneHotOneHotDescAUC 0.799
⚠ SimpleなOneHot特徴がESM-Sと統計的同等 — 4K件ではアーキテクチャより特徴工学が支配的
③ 提案手法の要点
  • TACK dataset: 3,514 PROTAC / 6,561エンドポイントを3DBから統合・標準化
  • Scaffold-based 5×5 CV: Murckoスキャフォールドでグループ化、化学空間外挿を厳密評価
  • 階層的統計検定: Friedman検定 + Wilcoxon + BH補正でp-hacking防止
  • 特徴アブレーション: 10特徴セット × 2アーキテクチャを統計的同等性で選定
  • Caruana greedy ensemble: 500候補から最適サブセットを選択 + bagging(10×50%)
  • 校正済みUQ: std/entropyが予測誤差と有意正相関 → active learning支援
⑤ 限界点・残る議題
OneHot特徴は未知POI/E3を表現不可 — 新規ターゲット汎化に課題
3D情報（三次複合体構造）なし — 構造予測ツール成熟後に再評価が必要
CRBN/VHL偏重(98.7%) — 新興E3リガーゼへの汎化未検証
UQカバレッジ過小(1σ: 46.5% vs 理論68.3%) — aleatoric不確実性混入
Dmax予測R²=0.36 — 細胞コンテキスト依存性がML特徴で捉えられない
⑥ ケムインフォマティクスパイプラインへの応用
  • lib/docking: XGBoost活性スコアラー（ROC-AUC 0.85）をUniDockRunnerの後処理フィルタとして統合 → CRBN/VHL系PROTAC仮想スクリーニング
  • lib/molgen: pDC50回帰 + UQ(std)をMolgenYaml scorerとして登録 → 生成分子の分解ポテンシャル評価 + active learning
  • 汎用評価: scaffold-based GroupKFold + BH補正統計フレームワークを任意ML評価ユーティリティとして整備
  • lib/fep × MD hybrid: DmaxにMD由来E3発現量・局在特徴を追加 → 細胞文脈を補完

コード/データ全公開: github.com/ribesstefano/TACK / huggingface.co/datasets/ailab-bio/TACK / zenodo.org/uploads/15691822

arXiv:2605.19579 | KDD AI for Sciences Track '26 | Ribes, Dunlop & Mercado — Chalmers University of Technology & Univ. of Gothenburg | 2026年5月