UniDock-Pro: Unified GPU-Accelerated Platform for High-Throughput Virtual Screening
Journal of Chemical Information and Modeling, 66, 2735-2752 | 2026 | DOI: 10.1021/acs.jcim.5c02587
SBVS・LBVS・Hybrid を 1 つの GPU エンジンに統合し、EF1% を 2.42 倍に押し上げ、1 日あたり数百万化合物のウルトララージ VS を実用化する。
① 背景と課題

創薬初期のバーチャルスクリーニング(VS)では、レセプター構造を使う SBVS と既知活性化合物の形状/ファーマコフォアを使う LBVS が相補的な情報源となる。しかし両者は別々のソフトウェア(AutoDock Vina、ROCS、Phase など)で運用されることが多く、エンジン間でスコアリング基準も座標系も異なるため、ハイブリッド戦略の構築には人手のパイプライン結合が必要だった。

既存の AutoDock-SS(先代 UniDock の LBVS モジュール)は再帰的ポテンシャルモデルを使うため、勾配が滑らかでなく GPU 上での共役勾配最適化と相性が悪かった。
数十億規模の make-on-demand ライブラリ(Enamine REAL 等)に対する 1 ノードでのスループットは、依然として SBVS 単独のドッキング速度がボトルネック。

→ SBVS・LBVS・Hybrid を同一 GPU カーネル上で動かし、参照リガンド情報をオンザフライでスコアに混ぜ込めば、早期濃縮と速度を同時に引き上げられるはず。

② 手法の概要 — 3 モード統合エンジン
  • SBVS:レセプター由来のポテンシャルグリッドを用いた従来型ドッキング
  • LBVS:参照リガンドのファーマコフォアグリッドを生成し、L-J 様の滑らかなポテンシャルへ置換
  • Hybrid:両グリッドを重み付き和でオンザフライ評価し、共役勾配最適化を 1 パスで実行
  • FFCA(Force Field Complementarity Analysis):受容体場とリガンド場のコサイン類似度を格子点ごとに算出し相乗領域を可視化
Hybrid モードのデータフロー 受容体格子 参照リガンド ファーマコフォア 重み付き和 w_R V_R + w_L V_L GPU 共役 勾配最適化 FFCA で領域ごとの補完性を後段解析
③ 本研究で示したこと
  • LBVS モード単独で EF1% を AutoDock-SS の 2.42 倍(DUDE-Z)
  • Hybrid は DUDE-Z/VSDS-vd TrueDecoy の双方で競争力ある早期濃縮
  • 1 GPU・1 日で数百万化合物の処理を実演(ウルトララージ VS 実用域)
  • FFCA により "力場が補完しあう領域" とドッキング精度向上の対応を定量化
  • 滑らかな L-J 様ポテンシャルが共役勾配と整合し収束ステップ数を削減
④-a EF1% 比較(DUDE-Z, LBVS)
早期濃縮率 EF1% — AutoDock-SS vs UniDock-Pro LBVS 0 x1 x2 x3 x1.00 AutoDock-SS 先代 LBVS x2.42 UniDock-Pro LBVS(本研究)
④-b モード別早期濃縮(DUDE-Z + VSDS-vd)
SBVS / LBVS / Hybrid のベンチマーク位置づけ low mid high → 偽陽性に対する厳しさ(DUDE-Z 厳格 / VSDS-vd TrueDecoy) SBVS LBVS Hybrid レセプタ単独 x2.42 改善 両ベンチで競争力
④-c スループット — 1 GPU・1 日
スクリーニング規模ファネル(1 GPU / 1 day) 数百万化合物 投入 make-on-demand ライブラリ想定 Hybrid フィルタ通過 受容体+ファーマコフォア両方を満足 EF1% Top x2.42 濃縮(vs AutoDock-SS) ※ GPU 型番は論文で詳細未記載
④-d FFCA — 力場補完性の解析軸
解析軸計算内容解釈
受容体場V_R(x)レセプタ由来ポテンシャル
リガンド場V_L(x)参照ファーマコフォア
FFCAcos(V_R, V_L) /格子点補完性の高い空間領域を特定
相関FFCA × ドッキング精度高補完領域で Hybrid 優位
FFCA は予測自体ではなく "なぜ Hybrid が効くか" の説明ツール。リード最適化時のフラグメント置換場所のヒントとして使える。
⑤ テイクホームメッセージ
3 モード単一エンジン化

SBVS・LBVS・Hybrid を 1 本の GPU カーネルにまとめたことが最大の貢献。座標系・スコアの整合が取れ、ハイブリッド戦略をユーザー側で自前実装する必要がなくなる。

滑らかな LBVS ポテンシャル

L-J 様への置換で共役勾配最適化と整合。AutoDock-SS の再帰モデル時代と比べ EF1% が 2.42 倍に達した実装上の鍵。

参照リガンドが必須という制約

Hybrid は既知活性化合物(参照)を要求する。完全 de novo の標的では SBVS 単独に頼るしかなく、運用フェーズの選択が必要。

ウルトララージ VS の実用ライン

1 GPU・1 日で数百万処理という指標は、Enamine REAL クラスのライブラリを単一ノードで回す現実解として機能する(GPU 型番は要追試)。

応用補足 — lib/docking への取り込み
  • UniDockRunner を Pro 仕様に拡張:mode={sbvs, lbvs, hybrid} と reference_ligand を受け付けるインタフェースに改修
  • ファーマコフォアグリッド生成をサブモジュール化し、ProLIFCalculator の interaction fingerprint を参照リガンド側のシードに使う
  • FFCA を解析 API として公開し、相互作用の "肝" 領域をフラグメント設計に橋渡し
  • 滑らかポテンシャルの実装はオプション化し、旧 AutoDock-SS 互換モードも残せる
インパクト
  • lib/docking の UniDockRunner を直接的に上位互換化できる(実装優先度 high)
  • Hybrid モードの導入で社内のリード最適化ループに参照リガンド主導のリスコアリングを組み込める
  • FFCA を MMGBSAEngine と組み合わせれば、結合親和性向上の物理的根拠を可視化する報告書が作りやすくなる