AutoDock Vina 分子对接快速上手3 步跑通你的第一次对接【免费下载链接】AutoDock-VinaAutoDock Vina项目地址: https://gitcode.com/gh_mirrors/au/AutoDock-Vina如果你第一次打开 AutoDock Vina 的终端大概率会经历这样一段事故现场盯着--receptor、--exhaustiveness这些参数发呆照抄一段教程命令却报出cannot find box好不容易跑起来又看不懂满屏的数字。别急——AutoDock Vina 是当今全球使用最广、速度最快的开源分子对接引擎它的设计初衷就是开箱即用你只需要准备好两个分子文件、圈出结合位点剩下的构象搜索与打分交给它。这篇文章要做的就是陪你把从装不上到跑通、再到跑得又快又好的这条路完整走一遍。第一幕 装不上、跑不通先搞清楚你卡在哪一环新手最容易犯的错是把分子对接当成一个黑箱扔进两个分子吐出结合能。于是第一次失败时他们毫不犹豫地把锅扣在 Vina 头上。但你看一眼上方的流程图就会明白全流程其实拆成了三个独立阶段结构预处理 → 对接输入准备 → 对接计算。实测下来你遇到的报错里九成发生在第二阶段之前——配体还是 SDF 格式没转成 PDBQT、受体带着水分子和辅因子没清理、盒子坐标根本没圈住结合位点。换句话说Vina 本身很少出错出错的是喂给它的文件还不够格。这个项目由 Scripps 研究所的 Forli 实验室持续维护以 Apache 2.0 开源协议免费分发。它内置 Vina、Vinardo、AutoDock4 三套打分函数支持大环分子、水合对接、多配体同时对接和批量虚拟筛选还附赠 Python 3 绑定。看到这里你应该换个心态它不是玄学黑箱而是一台需要你正确装料的精密仪器。第二幕 破除误区网格图的事真不用你操心很多人以为做分子对接之前必须先跑一遍 AutoGrid手工算好每种原子的亲和力网格图那一堆.map文件再让 Vina 去查表。这是 AutoDock4 时代留下的肌肉记忆而 AutoDock Vina 的设计哲学恰恰相反。项目官方文档里写得很直白用户不需要理解实现细节、不需要调整晦涩的搜索参数、不需要聚类结果、更不需要懂四元数代数。试着把 Vina 想象成一对磨合多年的双人舞搭档配体是那个不断变换步伐的领舞者受体则是舞台上纹丝不动的布景。Vina 一轮轮地尝试新的舞步组合构象、朝向、扭转打分函数像裁判一样即时给出这段合不合拍的分数然后它自动朝更协调的方向继续试探。你完全不需要亲手把舞台地板一块块擦亮——选用 Vina 力场时它会在内部默默准备好自己需要的一切。那网格图是不是彻底没用了也不是。当你切换 AutoDock4 力场或者做水合对接这种特殊协议时才需要autogrid4出面生成地图。理解这条分界线能帮你省下不少白折腾的时间多数入门场景跳过 AutoGrid 完全没问题。第三幕 10 分钟跑通第一次 AutoDock Vina 分子对接闲聊到此为止跟着做一遍。示例数据都在example/basic_docking/data/里官方逐步教程见docs/source/docking_basic.rst。我们用抗癌药伊马替尼对受体 c-AblPDB 条目 1iep做一次标准的刚性对接。第 1 步装齐工具链推荐用 conda 建一个干净环境把 Vina 本体和文件准备工具 Meeko 一次装齐conda create -n vina python3 -y conda activate vina conda install -c conda-forge numpy scipy rdkit vina meeko -y这一步解决装得上的问题Meeko 负责把结构文件翻译成 Vina 认得的 PDBQTVina 负责实际计算缺了任何一环后面都会卡在半路。第 2 步把受体和配体翻译成对接语言PDB 格式不带键连接信息而 PDBQT 才是 Vina 的母语。用 Meeko 的两条命令完成翻译cd example/basic_docking/data 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 mk_prepare_ligand.py -i 1iep_ligand.sdf -o 1iep_ligand.pdbqt这一步解决格式对的问题第一条命令生成受体 PDBQT同时借-v参数写出盒子坐标的配置文件和一个可视化用的 box PDB第二条把配体的 SDF 转成 PDBQT。20 Å 的盒子足够罩住伊马替尼的结合口袋——如果盒子画偏了后面一切努力都会白费。第 3 步跑起来读结果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 8 --out 1iep_out.pdbqt --log 1iep.log这一步解决跑得动、看得懂终端会逐行打印每个构象的亲和力分数单位 kcal/mol负值越大结合越强所有构象一并写入1iep_out.pdbqt日志存进1iep.log。看到分数刷出来的那一刻你就正式入门了。顺带记一个官方明说的技巧如果结果时好时坏、偶有找不到正确结合姿势把--exhaustiveness从默认的 8 提到 32代价只是多等一会儿。这个参数控制的就是 Vina 愿意下多大功夫去搜索构象空间。第四幕 进阶与避坑从能跑到跑得又快又好批量虚拟筛选几百个配体一次跑完适用场景你手里攒了一批候选分子想做虚拟筛选一个一个跑既慢又容易手滑漏文件。关键命令把--ligand换成--batch再用--dir指定输出目录输出文件名会自动按配体名命名vina --receptor 1iep_receptor.pdbqt --batch ligands/*.pdbqt \ --config 1iep_receptor.box.txt --dir poses常见坑与对策别把批量对接和多配体对接混为一谈——前者是串行筛选一批分子后者是把多个配体同时塞进一个结合位点算协同作用两者语义完全不同详细说明见docs/source/docking_in_batch.rst和docs/source/docking_multiple_ligands.rst。另外在 Windows 上通配符展开要格外小心建议用$(ls .\ligands\*.pdbqt)这类显式写法。水合对接把结合腔里的钉子户水算进去适用场景结合口袋里存在强烈保守的水分子它们实际上参与并稳固了配体结合尤其在片段类药物发现里忽略水分子可能导致姿势预测跑偏。关键命令先用example/autodock_scripts/mapwater.py给配体装上显式水W 原子完成对接后再用同目录的dry.py去掉未参与相互作用的水并重新打分。完整协议见docs/source/docking_hydrated.rst现成示例在example/hydrated_docking/。常见坑与对策水合对接依赖 AutoDock4 力场和autogrid4生成的改进版网格图别拿默认的 Vina 力场硬套否则显式水的存在意义会完全丢失。用 Python 脚本把对接塞进自动化管线适用场景参数扫描、重复实验、批量打分或者想把对接嵌入自己已有的分析流程。关键代码项目自带的最小示例在example/python_scripting/first_example.py核心逻辑只有几行from vina import Vina v Vina(sf_namevina) 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_poses20) v.write_poses(1iep_ligand_vina_out.pdbqt, n_poses5, overwriteTrue)常见坑与对策Python 绑定目前只覆盖 Linux 和 macOSWindows 用户请走 WSL。还有一个很容易踩的坑可执行文件vina和 Python 包vina是两条互相独立的安装路径——命令行能跑通不代表pip里能 import 成功需要时得单独pip install vina。结尾现在的你已经走完了从对着报错发呆到批量跑完几百个分子的完整旅程。下次当你盯着屏幕上不断刷新的亲和力分数时大概会想起最初连 box 都找不到的自己——那感觉就像第一次学会骑车后回头看见还不敢上车的朋友。去动手吧打开终端clone 一份源码git clone https://gitcode.com/gh_mirrors/au/AutoDock-Vina从example/basic_docking开始你的第一次分子对接就差这一条命令的距离。【免费下载链接】AutoDock-VinaAutoDock Vina项目地址: https://gitcode.com/gh_mirrors/au/AutoDock-Vina创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考