📚 2026年5月 月次論文レビュー — 🔬 計算化学

分子フレキシビリティは創薬や物性予測で頻繁に語られる概念でありながら、明確に定量化することが難しい。著者らは本研究で、柔軟性を「分子の全トーション角の運動範囲、すなわち取りうる全トーション状態のアンサンブル」と定義し、これを離散的に表現する新しい記述子 Torsion Angular Bin String (TABS) を提案した。TABS は 1 つのコンフォマーを、その各回転可能結合の二面角を離散ビンに量子化した整数ベクトル（文字列）として表す。さらに、ある分子について構成可能な distinct な TABS の総数を nTABS と定義し、これをコンフォメーション空間サイズの上限推定を与える 2D フレキシビリティ記述子として位置づけた。

医薬品の溶解度は経口バイオアベイラビリティを左右する重要物性であり、製剤開発から大規模製造まで創薬の全段階で精度の高い予測が求められる。COSMO-RS（Conductor-like Screening Model for Real Solvents）は量子化学計算から得られる分子表面の分極電荷密度（σ-profile）を統計熱力学と統合し、溶解度や溶媒和自由エネルギーを a priori に予測できる代表的手法だが、その精度は QM 計算レベルや基底関数のみならず入力として与える分子のコンフォマーに強く依存する。とりわけ分子内水素結合の有無は σ-profile を大きく変え、予測値を左右することが知られている。

本研究は、オープンソースかつ MIT ライセンスで開発された Open Free Energy（OpenFE）エコシステムの第一弾となる相対結合自由エネルギー（RBFE）計算プロトコルを、15 社の製薬企業が分担して大規模に独立検証した前競争的（precompetitive）協働ベンチマークである。アルケミカル自由エネルギー法は化合物の相対親和性をランク付けする創薬の標準ツールとなっているが、これまで産業利用に耐えるオープンソース実装は大規模検証・長期保守・高スループット最適化のいずれかを欠いていた。OpenFE はこの空隙を埋めるべく、学界（特に Chodera 研の Perses 法）と産業界の結節点で、…

構造ベース創薬におけるドッキングは、実験的に決定された高分解能タンパク質構造の有無に律速されてきた。本研究は、配列から直接 3 次元構造を予測する AlphaFold(AF)モデルがこの律速をどこまで埋められるかを、広く使われる Glide ドッキングを用いて体系的に評価したものである。Boston 大学の Maidanik、Joseph-McCarthy らのグループは、フラグメント・大型リガンドの双方について、AF モデルへのドッキング精度を実験的な結合(holo)構造・非結合(apo/unbound)構造と比較した。重要なのは、結合リガンドの座標を一切使わずに FTMap のホットスポットマッピングでドッキングボックスを定義する「非バイアス(ligand-agnostic)」手法を採った点で、…

OTMol は、二つの分子構造の重ね合わせ（superimposition）と RMSD 計算において本質的に重要な「原子間の対応付け（atom-to-atom correspondence）」を、最適輸送（Optimal Transport, OT）理論によって原理的に解く手法である。RMSD は柔軟リガンドの conformer 比較から分子クラスタ配置の比較まで広く使われる構造類似性指標だが、原子順序がソフトウェアや実験パイプラインごとに食い違う問題、水ネットワークや希ガス集合体のように等価な原子割当が多数存在する対称・不可識別配置の問題、そしてアライメント時の鏡像操作によってキラリティが反転してしまう問題を抱えている。これらはいずれも「正しい原子対応をどう決めるか」に帰着する。

ドッキングベースの仮想スクリーニング(VS)は創薬の初期ヒット探索で主役になったが、偽陽性率が高く、計算で優先順位づけしたヒットを実際に購入・実験検証する最終段は1化合物あたりのコストが他工程より桁違いに高い。著者らはこれを物流・通信の「ラストマイル問題」になぞらえ、この最終濾過にどの結合自由エネルギー(BFE)予測法が最もコスト効率よく真の活性化合物を濃縮できるかを問うた。本研究は、その問いに対して物理ベース法とAI/ML法を統一プロトコルで横並び比較した、著者らの知る限り初の大規模ベンチマークである。評価対象はアルケミカル絶対結合自由エネルギー(ABFE)、MM/PBSA、MM/GBSA の物理3法に加え、…

qMol は、創薬の初期ヒット探索における類似性検索を、フラグメントベースの reduced graph という表現を軸に再構成した Web サーバである。分子をそのまま原子レベルの指紋やグラフとして扱うのではなく、化学的に意味のあるフラグメント（主にリングシステムと適度な長さの側鎖・リンカー）をノード、それらをつなぐ共有結合をエッジとする縮約グラフへ変換する点が中核にある。各ノードにはフラグメント識別子（FragID）が割り当てられ、ライブラリ中の全化合物をこの reduced graph のトポロジー単位でグルーピングしておく。検索時には、query 分子と同一トポロジーを共有する化合物を厳密に取得することも、…

本研究は、構造ベースバーチャルスクリーニング (SBVS) の最大の弱点である「受容体構造の変動に対する脆弱性」を、分子表面パッチ表現によって克服する手法 PL-PatchSurfer3 (PL-PS3) を提案している。AutoDock Vina や DOCK6 に代表される従来のドッキングプログラムは、原子間距離に基づくスコアリング関数で結合自由エネルギーを推定するため、受容体のわずかな配座変化でも予測が大きく崩れる。実際、論文の引用する先行研究では、受容体を holo 形から apo 形に変えただけでリガンド RMSD が 1.34Å から 3.65Å に増大し、EF1% が 10.0 から 3.

本研究は、深層生成モデルが創薬パイプライン全体ではなく「ヒット同定」という一工程を直接代替できるかを問うた、著者らの主張では初めての系統的研究である。従来のヒット同定は HTS や仮想スクリーニングに依存し、数か月から数年と多大なコストを要する律速工程である。著者らは問いを「生成モデルは妥当でドラッグライクな分子を作れるか」という一般的なものから、「生成モデルはヒット様の特徴を持つ化合物を直接生成してヒット同定を加速できるか」というより具体的で要求の厳しいものへと再設定している。そのために、構造・物理化学・生物活性関連の指標を多段フィルタリングパイプラインに統合した評価フレームワークを提案し、…

芳香族複素環は医薬品の中核をなす骨格であり、近年 FDA 承認低分子薬の 82%（2013–2023 年）が窒素複素環を含むほど重要度が増している。複素環の入れ替え（バイオアイソステア置換・スキャフォールドホップ）は、活性を保ちながら物性や ADMET を調整する常套手段だが、合成済み・既知の置換空間に偏りがちで、承認薬が含むユニーク芳香族骨格は約 120 種に過ぎないとされる。Oxford の Brennan らは、未合成のものを含む芳香族複素環空間を網羅的に探索できるオープンソース Python ツール Heterocycle Isostere Explorer (HCIE) を開発した。