molecular-docking

分子ドッキング — パイプラインスキル（薄い足場）

概要

受容体／リガンドの調製ツールの選択、グリッドボックスのサイズ設定、探索の網羅性の調整、検証プロトコルに関する判断ロジックを追加します。誤ったポーズが気付かれないまま生成される原因となる落とし穴に重点を置いています。

使用方法

• ドッキングキャンペーンで調製ツール（obabel と ADFR／meeko）の選択時に有効化します
• グリッドボックスの寸法や配置を定義する際に有効化します
• 再ドッキングによるドッキング設定の検証時に有効化します

基本概念

判断ロジック

Receptor preparation tool:
├── Quick exploration / single ligand → obabel (-d, reduce, -xr)
└── Production virtual screen → prepare_receptor (ADFR) or meeko
    └── Handles altlocs, metals, non-standard residues correctly

Ligand preparation:
├── Single ligand from SMILES/SDF → obabel --gen3d -h
├── Large library (>100 ligands) → meeko (mk_prepare_ligand.py)
└── Stereocenters undefined? → Enumerate explicitly BEFORE docking

Grid box placement:
├── Co-crystal ligand available → Center on ligand centroid
└── No co-crystal → Use pocket predictor (fpocket, P2Rank, SiteMap)

Grid box sizing:
└── Ligand longest diameter + 10–15 Å padding per axis
    └── Typical: 20–30 Å per side for drug-like molecules

Exhaustiveness:
├── Quick exploration → 8 (Vina default)
├── Production / ranking → 32
└── Large/flexible ligands (>8 rotatable bonds) → 64+

検証ゲート（スクリーニング結果を信頼する前に必ず合格する必要があります）： 共結晶リガンドを再ドッキングし、結晶構造のポーズに対する RMSD が 2 Å 未満であることを確認します。失敗した場合は、スクリーニングの前に調製またはグリッドを修正します。

重要なパラメーター

パラメーター	推奨値	変更する場合
exhaustiveness	32	回転可能な結合が 8 本を超える場合は 64 以上
num_modes	10	アンサンブル解析では増やします
energy_range	3 kcal/mol	多様なポーズを捉えるために広げます
ボックスのパディング	+10–15 Å	アロステリック部位では大きくします
プロトン化 pH	7.0（デフォルト値から下げます）	酸性ポケットでは PROPKA／H++ を使用します（例：アスパラギン酸プロテアーゼ、pH 4.5）

スコアの解釈：

• Vina のスコアを比較できるのは、同じ受容体とグリッドの設定内に限られます
• 創薬向け化合物に典型的なヒットのスコア：-7 ～ -12 kcal/mol
• スコアは、異なる標的やグリッド設定間では比較できません

よくある間違い

• 誤り： 結晶構造の水を受容体に残す 正しい方法： PDBQT に変換する前に HOH／WAT 残基を削除します（明示的な水分子を含めたドッキングを行う場合を除く） 理由： 水が結合部位を塞ぎ、物理的に不自然なポーズを生みます

• 誤り： すべての結合部位でデフォルトの pH 7 のプロトン化状態を使う 正しい方法： 酸性／塩基性ポケットでは PROPKA または H++ を使用します 理由： 誤ったプロトン化状態は水素結合ネットワークを変化させ、誤ったポーズを生みます

• 誤り： グリッドボックスを結合部位ぎりぎりの大きさにする 正しい方法： リガンドの直径に、各軸で 10–15 Å のパディングを加えます 理由： パディング不足は部位を切り取り、ポーズをボックスの中心付近に偏らせます

• 誤り： 大型または柔軟なリガンドにデフォルトの exhaustiveness = 8 を使う 正しい方法： 本番環境では 32 以上、回転可能な結合が 8 本を超える場合は 64 以上を使います 理由： 探索の網羅性が低いと大域的最小値を見逃し、ランキングの信頼性が低くなります

• 誤り： スクリーニングの前に再ドッキングによる検証を省略する 正しい方法： 共結晶リガンドを再ドッキングし、RMSD が 2 Å 未満であることを確認します 理由： RMSD が 2 Å 以上の場合、調製またはグリッドが誤っています。スクリーニング結果は信頼できません

• 誤り： 本番環境のバーチャルスクリーニングに obabel -xr の PDBQT を使う 正しい方法： 本番環境では prepare_receptor（ADFR）または meeko を使用します 理由： obabel は altloc、金属、または非標準残基を見落とす場合があります

• 誤り： 定義されていない立体中心を含む SMILES に対して obabel --gen3d を実行する 正しい方法： 異性体を明示した SMILES を指定するか、立体異性体を明示的に列挙します 理由： 未定義の立体中心は任意に決定されるため、誤った鏡像異性体がドッキングされる可能性があります

• 誤り： 異なる受容体またはグリッド設定間で Vina の親和性を比較する 正しい方法： 同じ受容体とグリッドの設定内でのみ比較します。どちらかを変更した場合は再ランキングします 理由： スコアは条件に依存し、異なる条件間では比較できません

応答形式

• まずユーザーが必要とするコマンドまたはコードを示し、その後に説明します
• 次の構成にします：入力の確認 → 動作するコード → 主要なパラメーターの説明 → 注意点
• 各タスクにつき、完全に動作する例を 1 つ示します。あらゆる代替案は示しません
• コードコメントは最小限かつ機能的なものにします（存在理由ではなく、内容を説明します）
• 目安：コード 50～100 行と、簡潔な補足説明