第一次做分子对接，很多人都把时间花在装软件、转格式、调参数上。

等Vina终于跑出一个“-8.5 kcal/mol”，总算松了一口气。

但做多了会发现：分子对接最容易出问题的地方，往往不是软件，而是软件运行之前。

第一步：先把蛋白选对

蛋白结构不是随便从PDB下载一个就能用。

同一个蛋白可能有很多结构：有的带配体，有的不带；有的是单体，有的是二聚体；有的口袋附近还有缺失残基。

做小分子对接时，如果有合适的共晶复合物，通常优先参考带配体的结构，因为它直接提供了真实结合位点的信息。

分辨率很重要，但不是唯一标准。

一个“更像真实结合状态”的结构，往往比单纯追求最高分辨率更重要。

第二步：别把蛋白上的“最大洞”当成结合位点

有共晶配体，先看共晶配体的位置。

没有共晶结构，再结合文献报道的关键残基、突变结果和口袋预测软件综合判断。

口袋预测软件能告诉你“哪里像个洞”，但不能自动证明这个洞就是功能位点。

Vina官方也提醒，搜索空间并不是越大越好。搜索区域过大，会增加采样难度，需要相应提高搜索强度。

第三步：受体和配体都要处理

蛋白一般需要检查缺失残基、氢原子、质子化状态、金属和辅因子。

配体则要注意3D构象、质子化状态、互变异构体和可旋转键。

尤其是金属离子、辅因子和关键水分子，不能机械地“一键全删”。

有些水只是结晶环境中的“路人”，有些却真正参与了蛋白—配体识别。

第四步：有共晶配体，先做一次redocking

这一步很值得做。

把共晶配体从蛋白中取出，再按照自己的参数重新对接回去，看能否基本恢复实验观察到的结合姿势。

这相当于先拿一道“有标准答案的题”测试自己的对接流程。

很多研究会以配体重原子RMSD约2 Å以内作为成功复现结合姿势的常用参考，但这个阈值不是所有体系都必须机械套用。

第五步：先看构象，再看分数

结果出来以后，不要第一眼只找最负的数字。

先打开pose看看：

配体是不是真的进入预期口袋？
有没有明显空间冲突？
疏水基团是不是待在合理的疏水区域？
关键残基能不能形成合理相互作用？
配体有没有被扭成很不自然的构象？

Vina官方也明确指出，不同靶点上的预测准确性差异很大，最好利用已知活性分子或共晶配体先评估自己的体系。

所以，一套比较稳妥的思路其实很简单：

选结构 → 定位点 → 处理受体和配体 → redocking → 正式对接 → 分析pose → 实验验证。

至于分子动力学模拟，要不要做取决于研究目的。

它可以进一步观察复合物在动态环境中的行为，但并不是每次分子对接都必须机械地再跑一遍MD。

总结：

分子对接其实很“危险”。

不是因为它难，而是因为——它太容易给你一个结果。

只要输入文件没问题，软件通常都会给出分数、pose和漂亮图片。

真正决定这张图能不能写进论文的，从来不是“跑出来了没有”，而是：

这个结合模式有没有依据？

参考文献

1. Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31:455–461.

2. AutoDock Vina Documentation. The Scripps Research Institute.

3. Hevener KE, et al. Validation of molecular docking programs for virtual screening against dihydropteroate synthase. J Chem Inf Model. 2009;49:444–460.