Open Babel 是一个化学数据互操作工具箱，被称为“通用化学翻译器”，其 Python 封装为 Pybel。它让研究者可以在 SMILES、PDB、SDF、CIF、Gaussian 输入/输出等格式之间无缝迁移数据。
用途
- 在不同化学文件格式之间转换，例如 SDF 转 PDB。
- 从 2D 结构或 SMILES 生成 3D 坐标。
- 使用 SMARTS 模式搜索子结构。
- 进行快速的分子对接准备，例如加氢、计算电荷。
- 运行基础力场最小化，支持 UFF、GAFF、MMFF94。
- 计算分子指纹与 Tanimoto 系数。
- 批量处理大型化学数据库（数百万分子）。
主要能力
- 转换引擎：Open Babel 以“输入格式 -> 内部 OBMol -> 输出格式”的流水线方式工作，可在转换过程中添加过滤与变换，例如去盐、加氢。
- 两层接口：OpenBabel（SWIG）是低层接口，直接访问 C++ 类（OBMol、OBAtom），复杂但功能最强；Pybel 是高层封装，更符合 Python 风格，推荐用于约 90% 的任务。
- 分子属性与原子访问：可读取分子式、分子量，遍历原子及其坐标，访问 SDF 标签/元数据。
- 子结构搜索：使用 pybel.Smarts 与 SMARTS 模式查找官能团。
- 3D 结构与力场：先加氢，再用 make3D 生成初始 3D（距离几何），随后用 optimize 进行局部优化，可选 uff、gaff、mmff94、ghemical，并读取能量。
- 低层 API：可直接操作 OBMol，手动添加原子、设置原子序数与坐标、遍历键信息。
- 描述符与相似度：支持 FP2（基于路径）、FP3（SMARTS 模式）、FP4（官能团）等指纹，并可用 | 运算符计算 Tanimoto 系数。
使用方式
- 安装：推荐通过 conda/mamba 安装：conda install -c conda-forge openbabel。
- 标准导入：from openbabel import openbabel as ob；from openbabel import pybel。
- 基本格式转换：用 pybel.readstring 或 pybel.readfile 读取分子，设置 mol.title 等元数据，再用 mol.write 写出为其他格式或文件。
- 批量转换：使用 pybel.Outputfile 与 pybel.readfile 迭代处理，可选加氢。
- 对接配体准备：读取 SMILES、加氢、make3D、用 gaff 优化，再写出 pdbqt 格式。
- 去盐：通过 OBMol.Separate() 分离片段，按原子数取最大片段。
- 性能优化：先用 FP2 指纹与 Tanimoto 筛选候选，再执行昂贵的 3D 或 SMARTS 匹配；若只需格式转换，可直接使用 OBConversion 类，比创建 pybel.Molecule 对象更快。
适用场景
- 化学数据格式互转与数据清洗。
- 从 SMILES 或 2D 结构生成 3D 起始结构。
- 官能团与子结构筛选。
- 分子对接前的配体准备。
- 大规模化学数据库的批量处理与相似性筛选。
限制与注意事项
- 不要混用 RDKit 与 Open Babel 对象，二者不兼容，需通过 SMILES/SDF 传递数据。
- Open Babel 出错时较为安静，读取后应检查分子是否为 None 或原子数是否为 0。
- 不含 @ 或 / 的 SMILES 会在转换中丢失立体中心信息。
- 复杂描述符与高级 QSAR 一般更推荐 RDKit；Open Babel 更擅长转换与 3D 工作。
- 生成 3D 前必须加氢，否则结构会失真或不正确。
- 处理大文件应使用迭代器，避免将数百万分子载入列表导致内存崩溃。
- 某些系统上可能找不到力场参数数据文件，必要时可手动设置 BABEL_DATADIR。
- make3D 有时会反转立体中心，可改用 gen3d 操作以更稳健地保留立体化学。
- 优化后应检查能量是否合理，异常能量通常意味着几何结构有问题。
- 文件写入器应显式关闭，或使用上下文管理器确保数据写入。
参考文档
- 官方文档：https://openbabel.org/docs/dev/
- Python API（Pybel）：https://openbabel.org/docs/dev/UseTheLibrary/Python_Pybel.html
- 搜索模式：openbabel.OBMol、pybel.readfile、pybel.readstring、mol.make3D