📚 课程目录第3章基本命令与工作流

第3章：基本命令与工作流

从零开始跑Open Babel，命令行就是你的“主战场”。大多数时候你都窝在终端里，靠几条子命令搞定配体处理。
我做计算化学这么多年，最深的感受就是：花10分钟把命令行练熟，后面能省下100小时的“手动改格式”时间。

这一章我们聚焦三个核心子命令：obabel、obconform 和 obminimize。另外，输入输出重定向也是天天用的技能，我一次性讲透。

3.1 三个核心子命令

先看一张总图，理清它们的分工和协作关系。

说白了，obabel 是瑞士军刀，obconform 是构象专家，obminimize 是能量修理工。多数管线就是这三者的组合。

3.1.1 obabel – 你可以理解为“格式翻译官”

obabel 是我用得最多的命令。它的基本语法很直观：

obabel input_file -O output_file [选项]

举个例子，把 PDB 转成 SDF：

obabel ligand.pdb -O ligand.sdf

如果你想一次转多个文件，用 *.pdb 通配符：

obabel *.pdb -O all.sdf -m

这里 -m 代表把每个分子输出成独立文件。我曾经有个项目需要处理 2000 个配体，-m 配合文件列表，一行命令搞定。

小技巧： 加氢操作几乎每次都用。我个人习惯加 -h 参数：
obabel ligand.pdb -O ligand_h.sdf -h
这会自动加氢并修补不完整结构。

3.1.2 obconform – 生成“扭来扭去”的构象

做虚拟筛选之前，你常常需要生成配体的多个低能构象。obconform 就是用来确定“这个分子能扭成哪些姿势”。

obconform ligand.sdf -c 100 -o conformers.sdf

-c 100 指生成最多 100 个构象。实际输出可能少于这个数，因为分子自由度不够时它会自动停。

我踩过一个坑：默认输出是 多分子 SDF，每个 model 代表一个构象。如果你想要单独的 PDB 文件，需要加上 -m：

obconform ligand.sdf -c 50 -m -o conf_

结果会生成 conf_1.pdb、conf_2.pdb…… 非常好用。

3.1.3 obminimize – 把结构“掰直”到能量最低

构象生成后，不优化直接拿去做对接，常常出诡异结果。obminimize 通过力场让结构达到局部能量最小值。最简用法：

obminimize -ff MMFF94 -c 1e-4 ligand.sdf > minimized.sdf

-ff 选力场，我常用 MMFF94 或 UFF；-c 1e-4 是收敛标准（梯度）。你知道的，默认 UFF 有时候算小分子不够准，换成 MMFF94 结果更靠谱。

注意： 如果你对多个构象做优化，必须用 -n 指定线程数，否则默认单线程跑超慢。
比如 obminimize -ff MMFF94 -n 4 multi_conf.sdf -O opt.sdf

3.2 输入输出重定向与管道

命令行里 >、<、| 三个符号，在你写自动脚本时能省大量磁盘 IO。我经常直接管道连起来，不写临时文件：

obabel ligand.pdb -O ligand.sdf -h | obconform -c 100 -o confs.sdf

这其实绕过了 ligand.sdf 这个中间文件，直接把加氢后的内存对象喂给构象生成器。不过要注意，管道输入输出对某些格式要求 stdin/stdout 是 SDF 格式，如果你需要特定格式，最好还是用临时文件加 -m 稳妥。

还记得我刚开始用 Open Babel 时，老是记不住重定向符号顺序。其实很简单：

符号	含义	示例
`>`	标准输出重定向到文件（覆盖）	`obabel a.pdb -O a.sdf 2> log.txt`
`>>`	标准输出重定向到文件（追加）	`obabel a.pdb -O a.sdf 2>> log.txt`
`2>`	错误流重定向	`obabel bad.pdb -O out.sdf 2> err.log`
`\|`	管道，前一个命令的输出作为后一个的输入	`obabel *.sdf -O temp.sdf \| obminimize -ff MMFF94 -c 1e-3 -O opt.sdf`
`<`	从文件读取输入（较少用，因为 obabel 通常用文件名参数）	`obabel -isdf -Oopdb < list.sdf`

  经验谈： 写自动化脚本时，我喜欢把错误输出单独存下来：

  obabel *.pdb -O all.sdf 2> conversion_errors.txt

  这样一来，哪些文件转换失败一目了然，不用在终端里一行行找。

3.3 组合案例：完整配体预处理工作流

这个流程是我在一个真实项目中总结的，处理 500+ 个配体，从头到脚优化完，然后给 Autodock Vina 用。命令也就三四行：

# 1. 批量转换所有 PDB 为 SDF，同时加氢
obabel *.pdb -O processed_h.sdf -m -h --gen3d

# 2. 生成最多 200 个构象，输出为单独 PDB
for f in *_h.sdf; do
  obconform "$f" -c 200 -m -o "${f%.sdf}_conf"
done

# 3. 对所有构象做能量最小化
obminimize -ff MMFF94 -c 1e-4 -n 4 *_conf*.pdb -O final_opt.sdf

# 4. 最后转成 PDBQT（Vina 专用格式）
obabel final_opt.sdf -O docking_input.pdbqt -xr

你可能会问，为什么中间要用 --gen3d？因为输入 PDB 可能只有 2D 坐标甚至无坐标，--gen3d 会用 open babel 内部算法生成一个初始 3D 结构。另外，-xr 能把残基信息保留到 PDBQT，方便对接时识别。

真的要注意： 配体预处理时千万别跳过“加氢”和“优化”两步。我有过一次教训：直接拿没加氢的配体跑 docking，结果好几个氢键供体位点因为缺失氢而完全错判，浪费了三天跑结果。

3.4 避坑与习惯

文件列表太长？ 用 xargs 或 parallel。比如 ls *.pdb | xargs -I {} obabel {} -O {}.sdf。
奇怪的中断报错？ 试试 --errorlevel 1 降低日志等级，只看关键错误。
想确认转换是否成功？ 加 --verbose 看详细进度，但生产脚本里慎用，输出文件会膨胀。
obminimize 卡住不跑？ 检查力场是否支持当前分子。MMFF94 对小有机分子好，UFF 覆盖更广但精度弱一些。

好了，命令行这部分就这些。你试试在终端敲几条上面给的例子，立刻就能感受到 Open Babel 的威力。觉得不过瘾？自己找一个小的配体库，跑一遍完整的预处理管道，看看能不能跳过任何一步还保持结果靠谱。