分子对接预测

我第一次做 分子对接预测 是把一个客户给的 200 个化合物批量 dock 到靶点，按打分排了个序交差。结果实验验证只命中 3 个，客户问”为什么排第一的反而不中”，我才意识到打分函数给的是”可能性”不是”结论”。那次让我把 分子对接预测 的定位从”算命”改成”初筛 + 假设生成”。

一、为什么分子对接预测先做结合位点

对接前必须确认活性口袋，否则配体会乱跑到非功能腔。我用文献位点、保守残基或 SiteMap/FPocket 预测三者交叉定中心，半径取 10~12 Å 覆盖诱导契合空间。曾经一个激酶，我按默认网格中心放偏了 3 Å，top 构象全落在变构区，漏掉了真正的 ATP 口袋。做 分子对接预测 时，位点确认比参数调优更该先花时间。

二、核心原理：构象搜索与打分

分子对接预测 的 backbone 是采样（MC、遗传算法、RD) + 打分（经验/力场/知识）。Vina 用经验打分加局部优化，快但缺乏对溶剂和熵的处理；Glide 的 SP/XP 更细但慢。我一般先用 Vina 广筛，再对 top 50 用 Glide XP 重打分。关键认知：打分是相对序不是绝对亲和力，ΔG 误差常达 2~3 kcal/mol，不能直接当 Kd 报。

三、关键技术要点：水桥与柔性

很多口袋里有保守水桥，忽略它对接构象会错。我做 分子对接预测 时，对关键水用 WaterMap 或显式保留，让配体抢位或利用水。蛋白柔性方面，默认刚性受体是最大近似，重要残基我会做侧链柔性（Vina 不支持，用 induced-fit 或 MD 补）。曾有一个配体靠诱导契合才进袋，刚性对接打分垫底，柔性/MD 才翻上来。

四、实操流程：从准备到验证

我的步骤：受体去水/加氢/赋 AMBER 力场 → 定义位点 → 配体 3D 构象生成与质子化 → Vina 广筛 → top 候选 Glide XP 重打分 → 目检关键氢键与疏水接触 → 用短 MD（50~100 ns）验证稳定性。做分子对接预测最容易被省的是最后一步 MD 验证，但它是区分”假阳性”的关键，我默认保留。

实算一个激酶抑制剂筛选：Vina 广筛 500 个，top10 里 3 个 Glide XP 重打分上位，但 MD 一跑，其中 2 个在 30 ns 内漂出口袋，只有 1 个稳定结合。如果只交 Vina+Glide 结果，客户会拿那 3 个去湿实验，命中率 33%；补了 MD，命中率提到 100%（在验证集上）。这一轮验证多花两天算力，省下的是三轮失败的化合物合成。我现在对接项目默认把 MD 验证写进交付，不再把它当可选增值项，因为省掉它的代价往往在湿实验阶段才显现。

五、常见踩坑：电荷与质子化态

配体在生理 pH 下的质子化态错了，整个结合模式都变。我算过一个羧酸配体被当成中性，氢键网络全错。受体 His/Asp/Glu 的质子化也要按口袋环境判断，不能全默认。另一个坑是金属中心，对接软件大多不处理配位键，含 Zn/Mg 的靶点要特殊力场或手动约束，否则构象乱飞。

六、复盘：对接是”假设生成器”

回过头看，分子对接预测 最值钱的是快速缩窄候选范围、给出可检验的结合假设，而不是替实验下结论。我接项目时默认先要靶点信息和用途（筛选还是机制），再定工具链。把位点、打分、MD 验证三件事串好，对接结果才真能指导湿实验。被证明有用的，是愿意把 top1 先泼盆冷水看 MD 的清醒。下次有人只要“排个序”，我会先问——要不要我顺手用 MD 泼盆冷水。