MolWorld: Molecule World Models for Actionable Molecular Optimization
分子=ノード・MMP変換=エッジの「分子世界」を世界モデルで逐次拡張する actionable 最適化 ｜ arXiv:2605.08954 [cs.LG] (2026, Emory & Merck)
🎯 生成分子が既知系列から有効な局所変換(MMP)の連鎖で到達可能であることを保証しつつ性質を最適化する（actionable なリード最適化）。
① 背景と課題

リード最適化は単に高スコア分子を見つけるだけでは不十分で、その分子が既知化合物から有効な局所変換の連鎖で到達できる（actionable な類縁体系列の妥当な改変である）ことが要る。MolWorld は分子=ノード・有効局所変換(MMP)=エッジの molecule-transfer graph を導入し、最適化をこのグラフの逐次拡張（状態=グラフ、行動=分子生成、遷移=世界モデル）として定式化する。

de novo 生成（潜在空間/RL/拡散/構造ベース）は到達可能性を全く考慮しない
ペアワイズ最適化（graph-to-graph/分子翻訳/LLM）は単一分子からの一対一編集に閉じている
複数分子に条件づけた最適化（multi-molecule-conditioned）を扱う既存研究が存在しない

→ 単一出発分子への大域近接でなく「局所変換パスで既存系列から到達可能」という緩い制約で、構造化された分子世界の逐次拡張として最適化

② 手法: MolWorld（世界モデル RL）
MolWorld 閉ループ アンカー選択 性質+探索+多様性のビームサーチ(幅Ka) context-conditioned generator GNN文脈 C_Z → Transformer 自己回帰生成 オラクル評価 + 妥当性/重複除去 世界モデル Wφ グラフ挿入(p_φ(x∼u)>τ) → 状態更新
St+1 = Wφ(St, Xt), Xt ~ πθ(·|St)
SearchState St=(Gt, Ht, Pt)・MMP エッジは mmpdb で構築

ChEMBL 事前学習。生成分子は有効リンクを持つ場合のみ到達可能グラフに挿入。

③ 本研究で示したこと（要点）
  • リード最適化を「molecule-transfer graph の逐次拡張」として定式化（到達可能性を明示保証）
  • 複数分子に条件づけた最適化（multi-molecule-conditioned）を扱う初の枠組み
  • context-conditioned generator が複数関連分子の共有変換規則に整合した分子を生成
  • 高スコアを保ちつつ孤立ノード率をほぼゼロ・平均次数を大幅向上
④ 主な結果 (a) 性質最適化 AUC
0.9473 MolWorld 0.8491 MOLLEO ZINC250 JNK3 AUC@1

DRD2 0.9993・GSK3β 0.9977 も高水準で、ほぼ全タスクで最良/高競争力。

④ 主な結果 (b) 連結性(平均次数)
54.3 10.9 JNK3 58.9 32.4 DRD2 MolWorld ベスト比較 平均 transfer degree (↑良)

到達可能グラフの連結性が他手法を大きく上回る（JNK3 54.26 vs MOLLEO 10.90）。

④ 主な結果 (c) 孤立ノード率
0.0118 GSK3β 0.0039 DRD2 孤立ノード率 (↓良)

生成分子がほぼ全て既存系列から到達可能（孤立ノードほぼゼロ）。

④ 主な結果 (d) ドッキング・解釈
ドッキング全標的で最良
A2AR/DRD3/EGFR・両カットオフ(−8/−6)で成功率最大（PyScreener Vina）
actionable な多段パス
QED O1 0.827→G9 0.936 を化学的に解釈可能な局所編集の連鎖で到達
予算限定で優位顕著
オラクル予算が限られる PMV21(1000回)で特に強い
⑤ テイクホームメッセージ
🕸️ 分子世界
分子=ノード・MMP=エッジのグラフを世界モデルで逐次拡張。
🔗 到達可能性
生成分子は局所変換パスで既存系列から到達可能（孤立ノードほぼゼロ）。
🧩 多分子条件
複数の関連分子が共有する変換規則に整合した生成。
🧪 解釈可能
QED 0.827→0.936 へ化学的に解釈できる多段デザインパス。
ZINC250（MOLLEO 比）
指標MOLLEOMolWorld
JNK3 AUC@10.84910.9473
JNK3 平均次数10.9054.26
孤立ノード率高い0.0118
本研究のインパクト
  • lib/molgen: MolgenYaml/JobManager 生成を「molecule-transfer graph を状態とする閉ループ」でラップし、孤立的な高スコア外れ値でなく入力系列から到達可能な actionable 類縁体系列を出力（グラフ状態＋アンカー選択＋context-conditioned generator を追加、MMP は mmpdb で構築）
  • lib/docking: オラクル f を UniDockRunner(Vina) に差し替え構造ベース actionable 最適化。高親和性サブセットの孤立ノード率・平均次数を併報し「出発系列から辿れるヒット」を優先選別
  • lib/fep: MMP で連結されたエッジに限定して MMGBSAEngine/DockFEP を回し、相対自由エネルギー計算を化学的に意味ある小変換に集中させ計算効率と解釈性を向上
MolWorld: Molecule World Models for Actionable Molecular Optimization. arXiv:2605.08954 [cs.LG]. 2026 (Emory & Merck). https://arxiv.org/abs/2605.08954