告别手动操作在Win10上用批处理脚本一键搞定AutoDock Vina批量分子对接在生物信息学研究中分子对接是药物发现和蛋白质-配体相互作用分析的核心技术之一。AutoDock Vina作为一款开源的分子对接工具凭借其高效的算法和友好的使用体验已成为众多科研人员的首选。然而当面对数百甚至上千个配体分子需要对接时手动逐个操作不仅耗时费力还容易因人为失误导致结果不一致。本文将带你用Windows批处理脚本实现全流程自动化让计算机替你完成这些重复劳动。1. 准备工作与环境配置在开始编写批处理脚本之前需要确保所有必要的组件已正确安装和配置。以下是完整的准备工作清单AutoDock Vina安装从官网下载最新版本解压后建议将vina.exe、vina_split.exe等可执行文件放入专用工作目录文件格式转换确保所有配体和受体文件已转换为.pdbqt格式可使用OpenBabel批量转换工作目录结构/DockingProject ├── vina.exe ├── vina_split.exe ├── conf.txt ├── receptor.pdbqt └── ligands/ ├── ligand1.pdbqt ├── ligand2.pdbqt └── ...提示建议使用7-Zip或类似工具创建自解压包将Vina程序和环境配置打包方便在不同计算机间迁移项目。配置文件conf.txt是对接参数的核心典型内容如下receptor receptor.pdbqt center_x 15.2 center_y 22.7 center_z 18.3 size_x 25 size_y 25 size_z 25 num_modes 10参数说明参数名说明center_x/y/z对接盒子中心坐标(Å)size_x/y/z盒子各边长度(Å)num_modes输出构象数量2. 基础批处理脚本编写对于从PubChem下载的配体文件通常以Conformer3D_开头我们可以编写如下批处理脚本echo off setlocal enabledelayedexpansion :: 创建结果总览文件 echo 对接结果汇总 summary.csv echo 配体名称,最佳亲和力(kcal/mol),RMSD lower,RMSD upper summary.csv :: 批量对接循环 for %%f in (Conformer3D_*.pdbqt) do ( echo 正在处理配体: %%~nf if not exist %%~nf mkdir %%~nf :: 执行对接命令 vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt :: 提取关键结果数据 for /f tokens1,2 delims: %%a in (findstr Affinity %%~nf\log.txt) do ( set affinity%%b ) for /f tokens1,2,3 delims, %%a in (findstr RMSD %%~nf\log.txt) do ( set rmsd_l%%b set rmsd_u%%c ) :: 写入汇总文件 echo %%~nf,!affinity!,!rmsd_l!,!rmsd_u! summary.csv ) echo 批量对接已完成 start excel summary.csv脚本优化点添加了进度显示实时反馈当前处理的配体自动提取关键结果指标亲和力、RMSD值生成结构化的CSV汇总文件可直接用Excel分析使用enabledelayedexpansion解决变量延迟问题3. 高级功能扩展3.1 多受体批量对接对于需要测试多个受体的情况可扩展脚本如下echo off setlocal enabledelayedexpansion :: 受体列表 set receptorsreceptor1.pdbqt receptor2.pdbqt receptor3.pdbqt :: 为每个受体创建结果目录 for %%r in (%receptors%) do ( if not exist results_%%~nr mkdir results_%%~nr ) :: 双重循环处理 for %%r in (%receptors%) do ( echo 正在处理受体: %%~nr copy %%r conf_temp.txt nul (echo receptor %%~nxr) conf_temp.txt findstr /v receptor conf.txt conf_temp.txt for %%f in (ligands/*.pdbqt) do ( echo 处理配体: %%~nf if not exist results_%%~nr\%%~nf mkdir results_%%~nr\%%~nf vina --config conf_temp.txt --ligand %%f --out results_%%~nr\%%~nf\out.pdbqt --log results_%%~nr\%%~nf\log.txt ) ) del conf_temp.txt echo 多受体批量对接完成3.2 结果可视化与筛选添加以下脚本自动生成可视化报告echo off :: 生成HTML报告 echo ^html^^body^^h1^对接结果报告^/h1^ report.html echo ^table border1^^tr^^th^配体^/th^^th^亲和力^/th^^th^RMSD^/th^^/tr^ report.html for /d %%d in (*) do ( if exist %%d\log.txt ( for /f tokens2 delims: %%a in (findstr Affinity %%d\log.txt) do set affinity%%a for /f tokens2,3 delims, %%b in (findstr RMSD %%d\log.txt) do set rmsd%%b,%%c echo ^tr^^td^%%d^/td^^td^!affinity!^/td^^td^!rmsd!^/td^^/tr^ report.html ) ) echo ^/table^^/body^^/html^ report.html start report.html4. 错误处理与性能优化4.1 常见错误处理在批处理脚本中添加错误检查机制:: 在主要对接命令后添加错误检查 vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt if errorlevel 1 ( echo [错误] 配体 %%~nf 对接失败 error.log echo 错误详情 error.log type %%~nf\log.txt error.log echo ---------------------------- error.log ) else ( echo 配体 %%~nf 对接成功 )4.2 性能优化技巧并行处理利用start命令实现简单并行set /a max_threads4, count0 for %%f in (*.pdbqt) do ( set /a count1 if !count! gtr !max_threads! ( wait set /a count1 ) start Docking_%%~nf /min cmd /c vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt )资源监控添加系统资源检查:: 检查内存使用 for /f tokens2 delims: %%a in (tasklist /fi imagename eq vina.exe /fo list ^| findstr 内存使用) do ( set mem%%a set mem!mem:MB! if !mem! gtr 2000 ( echo 警告内存使用过高(!mem!MB)暂停新任务 timeout /t 30 ) )断点续接记录已完成任务if not exist completed.txt ( type nul completed.txt ) for %%f in (*.pdbqt) do ( findstr /x %%~nf completed.txt nul if errorlevel 1 ( vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt echo %%~nf completed.txt ) )在实际项目中我发现将配体分批处理每批50-100个能更好地平衡性能和稳定性。对于特别大的数据集可以考虑使用Python等脚本语言结合批处理实现更复杂的任务调度。
告别手动操作！在Win10上用批处理脚本一键搞定AutoDock Vina批量分子对接
告别手动操作在Win10上用批处理脚本一键搞定AutoDock Vina批量分子对接在生物信息学研究中分子对接是药物发现和蛋白质-配体相互作用分析的核心技术之一。AutoDock Vina作为一款开源的分子对接工具凭借其高效的算法和友好的使用体验已成为众多科研人员的首选。然而当面对数百甚至上千个配体分子需要对接时手动逐个操作不仅耗时费力还容易因人为失误导致结果不一致。本文将带你用Windows批处理脚本实现全流程自动化让计算机替你完成这些重复劳动。1. 准备工作与环境配置在开始编写批处理脚本之前需要确保所有必要的组件已正确安装和配置。以下是完整的准备工作清单AutoDock Vina安装从官网下载最新版本解压后建议将vina.exe、vina_split.exe等可执行文件放入专用工作目录文件格式转换确保所有配体和受体文件已转换为.pdbqt格式可使用OpenBabel批量转换工作目录结构/DockingProject ├── vina.exe ├── vina_split.exe ├── conf.txt ├── receptor.pdbqt └── ligands/ ├── ligand1.pdbqt ├── ligand2.pdbqt └── ...提示建议使用7-Zip或类似工具创建自解压包将Vina程序和环境配置打包方便在不同计算机间迁移项目。配置文件conf.txt是对接参数的核心典型内容如下receptor receptor.pdbqt center_x 15.2 center_y 22.7 center_z 18.3 size_x 25 size_y 25 size_z 25 num_modes 10参数说明参数名说明center_x/y/z对接盒子中心坐标(Å)size_x/y/z盒子各边长度(Å)num_modes输出构象数量2. 基础批处理脚本编写对于从PubChem下载的配体文件通常以Conformer3D_开头我们可以编写如下批处理脚本echo off setlocal enabledelayedexpansion :: 创建结果总览文件 echo 对接结果汇总 summary.csv echo 配体名称,最佳亲和力(kcal/mol),RMSD lower,RMSD upper summary.csv :: 批量对接循环 for %%f in (Conformer3D_*.pdbqt) do ( echo 正在处理配体: %%~nf if not exist %%~nf mkdir %%~nf :: 执行对接命令 vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt :: 提取关键结果数据 for /f tokens1,2 delims: %%a in (findstr Affinity %%~nf\log.txt) do ( set affinity%%b ) for /f tokens1,2,3 delims, %%a in (findstr RMSD %%~nf\log.txt) do ( set rmsd_l%%b set rmsd_u%%c ) :: 写入汇总文件 echo %%~nf,!affinity!,!rmsd_l!,!rmsd_u! summary.csv ) echo 批量对接已完成 start excel summary.csv脚本优化点添加了进度显示实时反馈当前处理的配体自动提取关键结果指标亲和力、RMSD值生成结构化的CSV汇总文件可直接用Excel分析使用enabledelayedexpansion解决变量延迟问题3. 高级功能扩展3.1 多受体批量对接对于需要测试多个受体的情况可扩展脚本如下echo off setlocal enabledelayedexpansion :: 受体列表 set receptorsreceptor1.pdbqt receptor2.pdbqt receptor3.pdbqt :: 为每个受体创建结果目录 for %%r in (%receptors%) do ( if not exist results_%%~nr mkdir results_%%~nr ) :: 双重循环处理 for %%r in (%receptors%) do ( echo 正在处理受体: %%~nr copy %%r conf_temp.txt nul (echo receptor %%~nxr) conf_temp.txt findstr /v receptor conf.txt conf_temp.txt for %%f in (ligands/*.pdbqt) do ( echo 处理配体: %%~nf if not exist results_%%~nr\%%~nf mkdir results_%%~nr\%%~nf vina --config conf_temp.txt --ligand %%f --out results_%%~nr\%%~nf\out.pdbqt --log results_%%~nr\%%~nf\log.txt ) ) del conf_temp.txt echo 多受体批量对接完成3.2 结果可视化与筛选添加以下脚本自动生成可视化报告echo off :: 生成HTML报告 echo ^html^^body^^h1^对接结果报告^/h1^ report.html echo ^table border1^^tr^^th^配体^/th^^th^亲和力^/th^^th^RMSD^/th^^/tr^ report.html for /d %%d in (*) do ( if exist %%d\log.txt ( for /f tokens2 delims: %%a in (findstr Affinity %%d\log.txt) do set affinity%%a for /f tokens2,3 delims, %%b in (findstr RMSD %%d\log.txt) do set rmsd%%b,%%c echo ^tr^^td^%%d^/td^^td^!affinity!^/td^^td^!rmsd!^/td^^/tr^ report.html ) ) echo ^/table^^/body^^/html^ report.html start report.html4. 错误处理与性能优化4.1 常见错误处理在批处理脚本中添加错误检查机制:: 在主要对接命令后添加错误检查 vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt if errorlevel 1 ( echo [错误] 配体 %%~nf 对接失败 error.log echo 错误详情 error.log type %%~nf\log.txt error.log echo ---------------------------- error.log ) else ( echo 配体 %%~nf 对接成功 )4.2 性能优化技巧并行处理利用start命令实现简单并行set /a max_threads4, count0 for %%f in (*.pdbqt) do ( set /a count1 if !count! gtr !max_threads! ( wait set /a count1 ) start Docking_%%~nf /min cmd /c vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt )资源监控添加系统资源检查:: 检查内存使用 for /f tokens2 delims: %%a in (tasklist /fi imagename eq vina.exe /fo list ^| findstr 内存使用) do ( set mem%%a set mem!mem:MB! if !mem! gtr 2000 ( echo 警告内存使用过高(!mem!MB)暂停新任务 timeout /t 30 ) )断点续接记录已完成任务if not exist completed.txt ( type nul completed.txt ) for %%f in (*.pdbqt) do ( findstr /x %%~nf completed.txt nul if errorlevel 1 ( vina --config conf.txt --ligand %%f --out %%~nf\out.pdbqt --log %%~nf\log.txt echo %%~nf completed.txt ) )在实际项目中我发现将配体分批处理每批50-100个能更好地平衡性能和稳定性。对于特别大的数据集可以考虑使用Python等脚本语言结合批处理实现更复杂的任务调度。
