AutoDock Vina分子对接从0到1完整实战流程、参数解读与避坑指南【免费下载链接】AutoDock-VinaAutoDock Vina项目地址: https://gitcode.com/gh_mirrors/au/AutoDock-VinaAutoDock Vina 是计算化学与药物研发领域被引用最广的开源分子对接引擎之一它的任务可以浓缩成一句话判断一个候选小分子配体会以什么姿态、多大结合亲和力进入蛋白质受体上的口袋。无论你是刚接触结构生物学的学生还是想搭建虚拟筛选管线的研发人员本教程都从一张真实运行的“成绩单”倒推整套流程带你一步步跑通分子对接并把新手最容易卡住的环节提前拆开讲明白。一次真实对接的终端输出长什么样先认识 mode、affinity 与 RMSD很多人第一次跑通 Vina盯着终端却不知道自己在看什么。一次对接顺利结束后屏幕上会出现这样一张表mode | affinity | dist from best mode | (kcal/mol) | rmsd l.b.| rmsd u.b. ------------------------------------- 1 -13.23 0 0 2 -11.29 0.9857 1.681 3 -11.28 3.044 12.41 4 -11.15 3.813 12.24三个核心名词先记在脑子里名词人话解释mode算法找到的一个“候选结合姿势”pose按得分从好到差排序affinity预测的结合亲和力单位 kcal/mol负值代表结合过程放能越负通常意味着结合越紧RMSD该姿势与参考姿势的空间差异l.b./u.b. 是两种计算口径单位 Åmode 1 和它自己比当然是 0mode 1 通常是所有候选姿势中最优的那个也是后续做可视化、做分析时最常关注的起点。但请记住这张表只是终点难的是走到终点之前的那几步。下面我们倒推一下拿到这张成绩单需要经历什么。跑通分子对接的前置准备工具链、PDBQT 输入文件与搜索盒子一次对接背后其实是一条流水线从配体/受体的原始结构出发经过预处理变成 Vina 认识的格式再划定搜索区域最后才轮到引擎上场。整条管线可以用项目仓库里这张流程图一眼看全用一条 pip 命令装齐 AutoDock Vina 和 Meeko 预处理工具Vina 本身负责“算”但输入文件的准备工作推荐交给配套的 Meeko 工具包。建议在 conda 或 venv 虚拟环境里执行pip install -U numpy scipy rdkit vina meeko这条命令装入了计算引擎vina和三个预处理命令行工具mk_prepare_ligand.py、mk_prepare_receptor.py、mk_export.py。装完可以敲一下vina --version确认引擎可用。想直接看官方示例的话可以克隆整个项目仓库git clone https://gitcode.com/gh_mirrors/au/AutoDock-Vina仓库里的example/basic_docking/目录自带完整的输入与预期输出非常适合边看边练。为什么不能把 PDB 文件直接丢给 Vina这是新手第一道坎Vina 的输入格式是 PDBQT而不是你从数据库下载的 PDB。PDB 文件里没有键连关系用它描述小分子会丢信息很多情况下直接导致对接结果不可信。PDBQT 相当于分子界的“统一制服”它在 PDB 坐标基础上补上了原子类型、部分电荷、可旋转键等信息Vina 只认这套格式。一句话总结预处理省掉的功夫都会变成结果里的坑。配体与受体的 PDBQT 制作命令以仓库自带示例1iep抗癌药伊马替尼与 c-Abl 激酶为例受体准备mk_prepare_receptor.py -i 1iep_receptorH.pdb -o 1iep_receptor -p -v \ --box_size 20 20 20 --box_center 15.190 53.903 16.917这条命令把受体 PDB 转成1iep_receptor.pdbqt-p同时额外输出一个盒子参数文件-v加盒子参数。配体准备mk_prepare_ligand.py -i 1iep_ligand.sdf -o 1iep_ligand.pdbqt注意配体输入推荐用 SDF 而非 PDB——前者自带键连信息。如果你的小分子结构里还没有氢原子建议先用加氢工具处理后再走这一步成败有时就差在一个氢原子上。对接盒子中心与尺寸怎么定、最容易错在哪“盒子”就是 Vina 允许配体探索的三维搜索区域可以理解成用手电筒在黑暗房间里扫一块区域灯光照不到的地方一概不算。你需要给出中心坐标--box_center和三个方向的边长--box_size单位是 Å。这里有两个经典翻车点提前打预防针别把格点数当埃用。AutoDock4 时代盒子的尺寸常用格点0.375 Å/格表示Vina 直接按 Å 算混用会得到一个巨大无比的搜索区。盒子不是越大越好。搜索区超过 27000 ų约 30×30×30会触发警告算法要在更大的空间里碰运气结果反而更不稳定。原则是“刚好罩住结合口袋不贪大”。第一次对接的三种启动姿势直传参数、配置文件与 AutoDock4 力场输入文件备齐后启动对接有三种常见姿势按你的使用习惯任选。姿势一所有参数堆在命令行适合快速试跑一次性把条件写全vina --receptor 1iep_receptor.pdbqt --ligand 1iep_ligand.pdbqt \ --center_x 15.190 --center_y 53.903 --center_z 16.917 \ --size_x 20 --size_y 20 --size_z 20 \ --exhaustiveness 32 --out 1iep_ligand_vina_out.pdbqt--exhaustiveness控制搜索的“卖力程度”默认是 8像伊马替尼这种比较难缠的体系官方教程建议调到 32结果会更稳定代价是更长的计算时间。姿势二把条件写进配置文件参数多的时候用配置文件更清爽也好归档复现。新建一个config.txtreceptor 1iep_receptor.pdbqt ligand 1iep_ligand.pdbqt center_x 15.190 center_y 53.903 center_z 16.917 size_x 20 size_y 20 size_z 20 exhaustiveness 32 cpu 8 out 1iep_ligand_vina_out.pdbqt然后一句命令搞定vina --config config.txt姿势三切换 AutoDock4 力场走“预生成亲和力图谱”路线Vina 默认自带 Vina 力场会内部算好所有图谱但如果你想用 AutoDock4 力场例如做水合对接就得先让 AutoGrid4 预生成亲和力图谱autogrid4 -p 1iep_receptor.gpf -l 1iep_receptor.glg然后带着图谱跑vina --ligand 1iep_ligand.pdbqt --maps 1iep_receptor --scoring ad4 \ --exhaustiveness 32 --out 1iep_ligand_ad4_out.pdbqt提醒Vina 力场和 AutoDock4 力场的分数没有可比性别拿两边的能量直接对比那是拿苹果比橘子。对接成绩单的读法能量、RMSD 与“结果对不对”的判断标准跑完只是第一步会读结果才是关键能力。能量怎么读结合亲和力是负值越负越好但这个“好”是有上下文的。在 1iep 这个示例体系里接近晶体构象的 pose 用 Vina 力场打分大约在 -13 kcal/mol 附近AutoDock4 力场约 -14 kcal/mol。如果你的体系和教程一致跑出的分数和这个量级差距很大就该回头检查输入了。怎么判断 pose 对没对最硬核的办法是把对接结果和实验测定的晶体构象叠在一起算 RMSD通常 2 Å 就算复现得不错。如果你有已知活性的配体或天然配体结构先做一次“回对接”把已知配体放回原口袋验证流程比盲目批量跑要靠谱得多。结果转格式PDBQT 不是通用格式很多可视化软件比如某些版本的 PyMOL显示起来很怪。用 Meeko 转成 SDF能保留正确的键级与形式电荷mk_export.py 1iep_ligand_vina_out.pdbqt -s 1iep_ligand_vina_out.sdf进阶玩家还可以直接调打分函数的权重例如下面这条命令会把所有氢键的强度翻倍vina ... --weight_hydrogen -1.2 ...进阶玩法一柔性对接让指定残基的侧链也动起来标准对接把受体当成一尊雕像配体在雕像表面找位置。但真实蛋白在结合时往往会有侧链“让位”这就是柔性对接想解决的问题——允许口袋里的某几个残基侧链跟着配体一起动。以仓库示例1fpu为例把第 315 号苏氨酸链 A设为柔性残基mk_prepare_receptor.py -i 1fpu_receptorH.pdb -o 1fpu_receptor -p -v \ -f A:315 -a --box_size 20 20 20 --box_center 15.190 53.903 16.917-f A:315就是“指定柔性残基”的关键参数。这条命令会把受体拆成两个文件1fpu_receptor_rigid.pdbqt不动的刚性部分1fpu_receptor_flex.pdbqt会动的柔性侧链对接时把两个文件一起喂给 Vinavina --receptor 1fpu_receptor_rigid.pdbqt --flex 1fpu_receptor_flex.pdbqt \ --ligand 1iep_ligand.pdbqt --config 1fpu_receptor.box.txt \ --out 1fpu_ligand_flex_vina_out.pdbqt注意柔性残基数量不要贪多侧链越多搜索空间越复杂跑得越慢、越容易发散。详细教程在 docs/source/docking_flexible.rst。进阶玩法二水合对接把水分子也摆进结合口袋蛋白口袋里的水分子不是摆设有些水被牢牢锁在保守位点配体结合时未必能把它们全部挤走这些“钉子户”水会直接影响配体怎么摆。水合对接的思路是给配体装饰一圈可移动的“假想水”能放得下就加分和受体撞了就自动丢弃。实操分三步给配体加氢、加水生成带水标记的 PDBQTscrub.py 1uw6_ligand.sdf -o 1uw6_ligandH.sdf mk_prepare_ligand.py -i 1uw6_ligandH.sdf -o 1uw6_ligand.pdbqt -w-w是关键告诉 Vina 配体上带有一圈显式水W 原子。按 AutoDock4 力场预生成亲和力图谱水合对接是基于 AutoDock4 力场标定验证的所以绕不开autogrid4那一步。对接结束后用仓库example/autodock_scripts/里的辅助脚本收尾mapwater.py分析水分子在结果中的保留情况dry.py去掉多余的水并重新打分。对于片段级小分子比如碎片药物发现场景显式考虑水通常能明显改善构象预测。完整的套路见 docs/source/docking_hydrated.rst。进阶玩法三用 Python 脚本把对接升级成批量虚拟筛选流水线单打独斗跑一个配体只是热身真正的生产力在批量。Vina 提供 Python 绑定让你在脚本里完成“加载受体 → 算图谱 → 对接 → 导出”的全过程。仓库里的 example/python_scripting/first_example.py 是一个极好的入门模板from vina import Vina v Vina(sf_namevina, verbosity1) v.set_receptor(1iep_receptor.pdbqt) v.set_ligand_from_file(1iep_ligand.pdbqt) v.compute_vina_maps(center[15.190, 53.903, 16.917], box_size[20, 20, 20]) v.dock(exhaustiveness32, n_poses9) v.write_poses(1iep_ligand_py_out.pdbqt, n_poses5, overwriteTrue)其中sf_name可以换成ad4或vinardo来切换力场。想换成optimize()还能做局部能量优化用来消除手动摆放姿势里的空间冲突很顺手。批量筛选的窍门如果先调用compute_vina_maps()再设置配体Vina 会把 Vina 力场定义的全部 22 种原子类型图谱一次性算好。这样不管后面的配体库里出现什么原子类型图谱都齐了不必逐个配体重算——这正是虚拟筛选能跑起来的关键import glob from vina import Vina v Vina(sf_namevina) v.set_receptor(receptor.pdbqt) v.compute_vina_maps(center[15.0, 53.0, 16.0], box_size[20, 20, 20]) for path in glob.glob(library/*.pdbqt): v.set_ligand_from_file(path) v.dock(exhaustiveness16, n_poses1) v.write_poses(results/ path.split(/)[-1], overwriteTrue)跑完把每个配体的最优得分汇总就是一张可供排序的初筛榜单。脚本细节的逐步讲解见 docs/source/docking_python.rst。新手高频翻车点与排查思路速查表把社区里被问得最多的问题整理成一张速查表遇到症状直接对号入座症状原因药方报错can not open conf.txt文件明明存在文件浏览器把扩展名藏了实际文件名是conf.txt.txt打开“显示文件扩展名”确认文件名按老视频教程敲命令报 usage errors老版本选项是--all新版改成了--out使用vina --help核对当前选项在集群上报boost thread resource错误集群配置禁止创建线程Vina 没法并行联系系统管理员或换一台能开线程的机器同样参数跑两次结果不一样搜索算法带随机性这不是 bug想要复现就固定随机种子--seed出现搜索空间超过 27000 ų 的警告你很可能把盒子尺寸按格点0.375 Å填了确认尺寸单位是 Å把盒子改小改了配体部分电荷结果毫无变化Vina 用自己内部的静电处理忽略你喂的电荷想调静电就去改打分权重改电荷没用结果构象合理但氢原子的位置很奇怪Vina 是 united-atom 打分只关心重原子正常现象真正要管好的是输入结构的质子化状态指定了--num_modes 20输出却不足 20 个输出数量还受energy_range截断限制适当调大energy_range结果在 PyMOL 里显示得歪七扭八PDBQT 不是标准结构格式用mk_export.py转成 SDF 再看更完整的问题库在 docs/source/faq.rst遇到罕见报错先去翻一翻大概率能找到答案。下一步行动指南用好项目自带的“教材”目录学到这你已经不是对接小白了。接下来的成长路线项目仓库本身就替你铺好了对照练习example/basic_docking/solution/里放着一整套预期输出PDBQT、图谱、日志你跑出的结果可以直接和它对答案这种自带参考答案的练习在开源项目里相当难得。按需进阶基础掌握后按“柔性对接 → 水合对接 → 大环分子 → 锌金属蛋白 → 多配体”的顺序逐个攻破每个示例都有配套文档和现成数据。看源码想知道--exhaustiveness背后到底在干什么可以去 src/ 源码目录逛一圈入口在 src/main/main.cpp搜索算法在monte_carlo.cpp、quasi_newton.cpp这几个文件里。最后给你一个可以立刻执行的小任务克隆仓库进入example/basic_docking/data/照着本文的“姿势二”配好config.txt把exhaustiveness从 8 试到 32对比两次成绩单的差异。这一步做完你对“搜索强度”这个概念的理解会超过绝大多数只看教程不动手的人。分子对接是一门实践出真知的功夫理论看得再多也不如亲手跑通一次对接来得踏实。打开终端开始你的第一次 Vina 之旅吧。【免费下载链接】AutoDock-VinaAutoDock Vina项目地址: https://gitcode.com/gh_mirrors/au/AutoDock-Vina创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考