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

モレキュラードッキングは、創薬の構造ベース薬物設計（SBDD）において中核的な計算手法である。本論文は2026年4月にPreprints.orgで公開されたスコーピングレビューであり、古典的な物理ベース手法から最新の生成AI手法までを体系的に整理したものである。著者らはPRISMA-Scopingガイドラインに従い、2015年以降の文献を対象としてPubMed・Scopus・arXiv・ChemRxivなどを横断的に検索し、以下5つの研究課題への回答を構造化している。すなわち、(1)手法の種類と進化、(2)古典手法とAI手法の精度・物理的妥当性の比較、(3)HTPVSインフラの現状、(4)受容体柔軟性対応の手法、(5)ハイブリッドパイプラインの将来展望である。

共有結合阻害剤は、ウォーヘッドと呼ばれる反応性官能基がタンパク質の求核残基（主にシステインやセリン）と共有結合を形成することで、非共有結合型リガンドに比べて著しく高い結合力と選択性を発揮する。FDAが承認した薬物の約7%が共有結合型であることからも、その薬学的重要性は明確だが、従来の分子ドッキングプログラムは共有結合形成を適切にシミュレートできず、専用ツール（CovDock, GOLD, AutoDock等）もre-docking成功率40〜60%程度に留まっていた。本論文のCarsiDock-Covは、著者らが先に開発した非共有結合ドッキングツールCarsiDockのTransformerベース距離行列予測モデルを共有結合リガンドに転用し、…

スキャフォールドホッピングは、中核骨格を置換することで特許を回避しながら新規薬物候補を生成するメディシナルケミストリーの重要戦略である。ChemBounceはこれをオープンソースで実現するフレームワークで、入力SMILES分子をHierSアルゴリズム（ScaffoldGraph）で断片化して環系スキャフォールドを特定し、ChEMBLから構築した322万件の合成実証済みスキャフォールドライブラリとECFP4 Tanimoto類似度でマッチングして置換、さらにElectroShape（3D形状+電荷）フィルタリングでファーマコフォア保持を確認することで、合成アクセシブルな新規構造を体系的に生成する。

本論文は、2012年の創設論文「Solvent structure improves docking prediction in lectin-carbohydrate complexes」以来のバイアスドッキング研究10年間の成果をまとめたレビューである。MDシミュレーションから決定される溶媒サイト（WS・HS・PhS）をAutoDock4/Vinaのエネルギーグリッドマップに逆ガウスエネルギーウェルとして組み込む「知識ベースバイアスドッキング（AutoDock-Bias）」の原理・実装・適用範囲を包括的に解説する。

分子ドッキングのスコアリング関数は絶対結合親和性の予測精度に限界があり、汎用スコアリング関数（CHEMPLP・Vina等）は活性化合物の選択（スクリーニング能力）において不十分な性能を示すことが多い。SMARTDockはこの問題をターゲット特異的スコアリング関数（TSSF）の自動構築で解決するDockerベースのVSワークフローである。GOLDのCHEMPLPスコアリングから抽出したタンパク質原子別スコア寄与フィンガープリント（PADIF）と、ChEMBLの公開バイオアクティビティデータおよびDark Chemical Matterからの合成デコイを統合し、RF/XGBoost/MLPの3種の機械学習分類モデルをコンテナ内で完全自動訓練する。

スペイン・Chemotargetsとジローナ大学のCabot-March、Jalencas、Mestresが開発したフラグメントベース創薬（FBDD）支援のための合成的3Dフラグメントリポジトリ「SAFR（Structurally Augmented Fragment Repository）」を報告している。PDB結晶構造をテンプレートとした分子場類似性アラインメント（MIMIC）とrDock拘束ドッキングの組み合わせで、ChEMBLおよびBindingDBに収録された活性化合物（pAct≥5）に対して高信頼度の結合ポーズを計算により付与した。続いてBRICS断片化を適用し、フラグメントライク性フィルターを通じて157,080ユニークフラグメントと1,142タンパク質の3D相互作用情報を格納したリポジトリを構築…

Rutgers大学のYorkグループがAMBER/AmberToolsと統合したFEP（自由エネルギー計算）解析統合ツールスイート「FE-ToolKit」を報告している。主要コンポーネントは3つ：(1) **ndfes**：アンブレラサンプリングからMBAR/vFEP法で多次元自由エネルギー面を構築し、surface-accelerated string method（SASM）で最小自由エネルギーパスを探索、(2) **edgembar**：アルケミカルFEPシミュレーションのエッジを集約したネットワーク（熱力学グラフ）に対し、EdgeMBAR法でサイクルクロージャー条件を課した全体最適化を実施、…

plotXVGはGROMACS分子動力学シミュレーションの標準出力形式である.xvgファイルを入力とし、出版品質のグラフを一括生成するPythonツールである。Upssala大学のRosenbaum・van der Spoelらが開発し、Alexandria Chemistry Toolkit（ACT）プロジェクトの可視化コンポーネントとして実際の出版物で使用実績を持つ。約1,000行のPythonで実装されており、CLIとPython APIの両方を提供する。折れ線グラフ（MD時系列データ）、散布図（相関プロット）、ヒートマップ、等高線図をサポートし、Gibbs自由エネルギーランドスケープ計算（G = -k_B T ln(P/P₀)）やSciPyを用いたカーネル密度推定（KDE）も内蔵する。

RinPyはタンパク質の残基相互作用ネットワーク（RIN: Residue Interaction Network）モデルを実装したオープンソースPythonパッケージである。重原子コンタクトトポロジーに基づいてタンパク質構造をグラフとして表現し、ベットウィーネス・クロースネス・次数の3種類の中心性指標からアロステリック・機能的残基を特定する。さらにグラフラプラシアンのスペクトル解析により動的ドメインとヒンジ領域を特定し、リガンド結合・変異等の摂動前後での構造変化を定量的に比較できる。Istanbul Technical UniversityのKurkcuoglu研究室が開発し、PyPI（`pip install rinpy`）とGitHubで公開されている。

CReM-pharm は、3D ファルマコフォアモデルを制約条件として用い、フラグメントを反復的に成長させることで新規薬様分子を設計するde novoジェネレーターである。コアとなる技術は Pavel Polishchuk らが開発した CReM（Chemically Reasonable Mutations）フレームワークで、ChEMBL由来の実在フラグメントとその化学コンテキストを組み合わせることにより、化学的に妥当でかつ合成到達性の高い構造を生成する。ファルマコフォアモデルへの適合と合成容易性制御を同時に達成している点が最大の特徴であり、オープンソースとして公開されている（https://github.com/ci-lab-cz/crem-pharm）。