如何用Biopython修改PDB文件中配体氢原子的名称？

你遇到的问题根源在于：Biopython中Atom对象的id属性本质是二元元组（格式为(原子名称, 替代位置标识符)），而非单纯的字符串。直接赋值带空格的字符串会导致内部结构不兼容，保存PDB时程序仍会读取初始化时的原子名元组，而非你修改后的字符串。

以下是两种可行的解决方案：

方案一：正确修改Atom.id的元组结构 #

直接将Atom.id设置为包含新原子名和空替代位置的元组，这样会同步更新full_id，保存时即可生效：

from Bio.PDB import PDBParser, PDBIO

io = PDBIO()
target_pdb_path = 'mypdb.pdb'
pdb = PDBParser(QUIET=True).get_structure('target', target_pdb_path)[0]

hydrogens = []
for atom in pdb.get_atoms():
    if atom.parent.id[0].startswith('H_') and 'H' in atom.name:
        hydrogens.append(atom)

# 重命名配体氢原子
for h_num, h in enumerate(hydrogens, 1):
    # 正确设置id为(新原子名, altloc)元组
    h.id = (f'H{h_num}', ' ')

io.set_structure(pdb)
io.save('test.pdb')

修改后查看h.full_id，会发现原子名部分已更新为新的元组（比如('H29', ' ')），输出的PDB文件也会使用新的氢原子命名。

方案二：自定义Select类控制输出 #

如果方案一在特定Biopython版本中仍有问题，可以通过继承Select类，在输出阶段动态修改原子名：

from Bio.PDB import PDBParser, PDBIO, Select

class HydrogenRenameSelect(Select):
    def __init__(self, hydrogen_map):
        self.hydrogen_map = hydrogen_map
    
    def accept_atom(self, atom):
        if atom in self.hydrogen_map:
            # 动态修改原子id元组
            atom.id = (self.hydrogen_map[atom], ' ')
        return True

# 读取PDB结构并筛选配体氢原子
pdb = PDBParser(QUIET=True).get_structure('target', 'mypdb.pdb')[0]
hydrogens = [
    atom for atom in pdb.get_atoms()
    if atom.parent.id[0].startswith('H_') and 'H' in atom.name
]

# 建立氢原子到新名称的映射
hydrogen_map = {h: f'H{i+1}' for i, h in enumerate(hydrogens)}

# 使用自定义规则保存PDB
io = PDBIO()
io.set_structure(pdb)
io.save('test.pdb', HydrogenRenameSelect(hydrogen_map))

这种方式通过拦截原子输出的逻辑，确保只有目标氢原子被重命名，其余原子不受影响。

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