AutoDock软件家族：AD4与Vina的版本演变
新手刚开始接触AutoDock时，最容易在哪件事上犯迷糊？选版本。打开官网一看，AutoDock4和AutoDock Vina两套工具摆在那，很多人直接懵了。我刚开始也是这样，心想“不都是AutoDock吗，有啥区别？”后来踩了几次坑才明白——这两兄弟虽然同根同源，但设计理念完全不同。
这章咱们就聊聊这两个版本的核心差异，以及它们各自该怎么选。我个人习惯把AD4比作“手术刀”，把Vina比作“冲锋枪”。前者精准但需要耐心，后者快速但不那么精细。
1. 两个版本的核心定位
| 对比维度 | AutoDock4 (AD4) | AutoDock Vina |
|---|---|---|
| 引擎类型 | 经典遗传算法 + 局部搜索 | MC随机扰动 + 梯度下降 |
| 构象搜索 | 拉马克遗传算法（LGA） | Broyden-Fletcher-Goldfarb-Shanno（BFGS） |
| 评分函数 | 半经验自由能力场 | 基于X-Score的经验评分 |
| 输出格式 | DLG文件，需要ADT解析 | 直接输出PDBQT和TXT |
| 适用场景 | 学术研究、蛋白-配体细节分析 | 高通量虚拟筛选、快速评估 |
我的个人建议：如果你在写文章时需要详细分析结合模式，或者研究的蛋白有特殊金属离子、水分子参与，选AD4。如果你只是做高通量筛选，或者项目时间紧需要快速出结果，直接上Vina。
2. 算法差异：一个求深度，一个求速度
AD4用的拉马克遗传算法，说白了就是把“遗传”和“拉马克进化”结合起来。每一代配体构象会经历交叉、变异，然后通过局部优化调整。我上学那会儿用AD4跑一个配体要半小时，但结果出来是真的精细，RMSD能到0.5Å以内。
Vina则不同。它用的是MC随机扰动配合BFGS梯度下降。什么意思呢？每次随机动一下配体，看能量降不降，降就收下，不降就按一定概率接受。这个方法的优势是快，非常快。我记得有回过一个虚拟筛选项目，用Vina跑了1000个配体，才花了不到三小时。
为什么要快？因为Vina把搜索空间做了简化处理。它对二面角旋转的处理不像AD4那么全面，对ligand柔性区域的描述也更粗粒度。所以你在Vina上跑出来的结果，趋势是对的，但细节不一定能复现。
关键区别：AD4擅长“精搜”，适合做1-2个配体的深入分析；Vina擅长“粗筛”，适合从几十万个小分子中挑出有潜力的候选物。
3. 评分函数的差异
AD4的评分函数是基于经验自由能力场。它把结合自由能拆成几个项：范德华力、静电力、氢键、去溶剂化能和扭转熵。你想想看，每个项都有独立的系数，需要大量拟合。这些系数来自对已有晶体结构的回归分析，所以AD4的评分体系本身就有一定的局限性——它更偏向于它训练集中出现过的蛋白类型。
Vina的评分函数则简单很多，它主要沿用了X-Score的思路，把范德华和静电作用打包成一套经验项，配合一个简单的疏水项。这个评分函数鲁棒性更好，但精度上不如AD4细致。我在项目中遇到过特别有意思的事：同一个配体-蛋白复合物，AD4给的打分是-9.2 kcal/mol，Vina给的是-8.7，实验测得是-9.0。嗯，AD4更准，但Vina也能用。
4. 输入输出与操作习惯
AD4的整个流程离不开AutoDock Tools（ADT）。你需要用ADT准备pdbqt文件、设置盒子、生成GPF和DPF参数文件。整个过程比较繁琐，但每一步你都能看到具体的参数。
# AD4的典型命令行
autodock4 -p ind.dpf -l ind.dlgVina就简单多了，你只需要准备好受体和配体的pdbqt文件，写一个配置文件，一行命令就能跑。
# Vina的典型命令
vina --receptor receptor.pdbqt --ligand ligand.pdbqt --config conf.txt --out out.pdbqt配置文件长这样：
# conf.txt
center_x = 15.0
center_y = 20.0
center_z = 10.0
size_x = 20
size_y = 20
size_z = 20
exhaustiveness = 8
num_modes = 9我曾经踩过的坑：在AD4中，如果你不小心把配体的起始构象放得离盒子太远，遗传算法可能会找不到正确的结合位点。即使跑了100代，结果也可能完全跑偏。后来我学乖了，每次跑之前都会先用Vina做一个快速粗搜，确认一下大致的结合位置，再用AD4精搜。
5. 版本图谱：一图看懂全局
下面我画了一张图，帮你理清AutoDock家族的进化路线和各自的任务场景。
6. 版本选择的实战建议
到底什么时候用AD4，什么时候用Vina？我给你列个快速判断指南：
- 你做的是虚拟筛选，配体数量超过100个：别犹豫，直接上Vina。AD4跑100个配体可能需要一整天，Vina两小时搞定。
- 你的蛋白有特殊金属离子（Zn、Fe、Cu）：用AD4。Vina对金属配位的处理非常粗糙，容易出假阳性。
- 你需要发表文章，评审可能会看RMSD和结合模式的细节：用AD4做完后，最好再用已知晶体结构做一个re-docking验证。
- 你只是想快速看看某个化合物有没有潜力，后续还有大量实验验证：Vina足够了，别在计算上浪费时间。
避坑指南：我曾经在一次项目里用Vina跑了一批配体，结果出来的打分排序跟实验活性完全相反。后来查了半天，发现是受体准备时加氢方向设错了。所以不管用哪个版本，受体和配体的预处理都是最关键的步骤。如果你在AD4上跑不通，可以先在Vina上试试，看是不是同一个问题。
7. 版本未来的趋势
最近几年AutoDock团队推出了更新的版本，比如AutoDock-GPU和AutoDock Vina 1.2.x。这些版本在底层做了不少优化：
- AutoDock-GPU：把AD4的评分函数和遗传算法搬到了GPU上跑，速度提升了50倍以上。我试过用一张RTX 4090跑AD4，100个配体十几分钟就跑完了。
- Vina 1.2.x：加入了更完善的柔性残基支持和新的评分项，精度比老版Vina提升了不少。但代价是速度稍微慢了一点。
嗯，这些新版本本章不展开，后面的章节会单独讲。你先把AD4和Vina这两个基础版本搞明白，后面的进阶版本学起来就顺手了。
好了，这一章你就记住了：AD4是手术刀，Vina是冲锋枪。选对武器，对接工作就成功了一半。