📅 5月第2週論文レビュー — 🔬 計算化学2026-05-08〜2026-05-14

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

スキャフォールドホッピングは、中核骨格を置換することで特許を回避しながら新規薬物候補を生成するメディシナルケミストリーの重要戦略である。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種の機械学習分類モデルをコンテナ内で完全自動訓練する。

gmx ffconv はランカスター大学のAaltonenが開発したGROMACS専用の力場変換ツールで、既存のGROMACSシミュレーション系を異なる全原子力場（例: CHARMM36m ↔ AMBER/Lipid21）に変換する半自動CLIツールである。既存のMDシミュレーション系を別の力場でも並列評価したい場面や、力場間のベンチマーク比較を行う前処理として活用される。特に膜シミュレーション領域でCHARMM形式とAMBER形式の相互変換ニーズが高く、これに対応したツールとして設計されている。なお本レビューはSI（Supporting Information）PDFのみを入力として作成したため、主論文の詳細（ベンチマーク結果・比較評価）については限定的な記述となっている。

スペイン・Chemotargetsとジローナ大学のCabot-March、Jalencas、Mestresが開発したフラグメントベース創薬（FBDD）支援のための合成的3Dフラグメントリポジトリ「SAFR（Structurally Augmented Fragment Repository）」を報告している。PDB結晶構造をテンプレートとした分子場類似性アラインメント（MIMIC）とrDock拘束ドッキングの組み合わせで、ChEMBLおよびBindingDBに収録された活性化合物（pAct≥5）に対して高信頼度の結合ポーズを計算により付与した。

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

plotXVGはGROMACS分子動力学シミュレーションの標準出力形式である.xvgファイルを入力とし、出版品質のグラフを一括生成するPythonツールである。Upssala大学のRosenbaum・van der Spoelらが開発し、Alexandria Chemistry Toolkit（ACT）プロジェクトの可視化コンポーネントとして実際の出版物で使用実績を持つ。約1,000行のPythonで実装されており、CLIとPython APIの両方を提供する。

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）。