Biopython的PDB模块疑似存在Bug：某残基坐标获取错误

坐标差异对比 #

原始PDB文件中4IP残基前4个原子坐标：

HETATM 5636 C1 4IP A 405 80.967 85.113 26.680 1.00 22.42 C
HETATM 5637 O1 4IP A 405 82.327 85.039 27.129 1.00 23.40 O
HETATM 5638 C2 4IP A 405 80.917 85.791 25.309 1.00 22.60 C
HETATM 5639 O2 4IP A 405 81.463 87.121 25.385 1.00 20.92 O

Biopython处理后保存的坐标：

HETATM 1 C1 4IP A 405 30.570 61.217 -13.415 1.00 22.42 C
HETATM 2 O1 4IP A 405 29.672 60.422 -14.201 1.00 23.40 O
HETATM 3 C2 4IP A 405 30.182 61.120 -11.938 1.00 22.60 C
HETATM 4 O2 4IP A 405 28.836 61.592 -11.740 1.00 20.92 O

排查步骤 #

• 检查Biopython版本：不同版本的PDB模块解析逻辑可能存在差异，执行以下代码查看当前版本，尝试升级至最新稳定版：

  import Bio
  print(Bio.__version__)
  

• 确认PDB文件完整性：确保使用的是直接获取的原始PDB文件，未经过其他工具预处理。
• 简化测试流程：排除其他代码干扰，仅执行加载-保存操作，验证问题是否复现：

  from Bio.PDB import PDBParser, PDBIO
  
  parser = PDBParser(QUIET=True)
  structure = parser.get_structure("2r09", "2r09.pdb")
  io = PDBIO()
  io.set_structure(structure)
  io.save("test_2r09.pdb")
  

• 排除对称操作影响：部分结构包含空间群或对称操作记录，Biopython可能自动应用变换导致坐标偏移。尝试禁用相关解析逻辑：

  parser = PDBParser(QUIET=True, PERMISSIVE=True)
  # 或遍历结构时跳过对称衍生的链/残基
  

• 检查自定义代码逻辑：确认是否在解析过程中无意应用了坐标变换（如平移、旋转），或对杂残基进行了特殊过滤处理。

内容的提问来源于stack exchange，提问作者TheEivill