使用BioPython的BioPDB修改蛋白质链ID后导出PDB失效求助

你的核心问题是没正确拆分残基到新链，而是反复修改原链的ID，导致导出时无法保留拆分后的结构。以下是问题分析和正确实现：

原代码的问题 #

• 第一种代码：遍历残基时不断修改同一链的ID，且每次只保存单个链，最终导出的是最后一次修改的链，丢失了另一部分残基
• 第二种代码：遍历残基时反复覆盖原链的ID，最终链ID会被最后一个残基的判断结果覆盖，要么全A要么全B，根本没实现拆分

正确实现代码 #

用Biopython的PDB模块，核心是创建新链B，把符合条件的残基从原链A移过去，再保存完整结构：

from Bio.PDB import PDBParser, PDBIO

# 解析原PDB结构
parser = PDBParser()
structure = parser.get_structure("heterodimer", "heterodimer/heterodimer.pdb")
model = structure[0]

# 获取原链A，创建新链B并添加到模型
chain_a = model["A"]
chain_b = model.add_chain("B")

chain_a_max_resid = 275

# 遍历原链的残基，转成列表避免遍历过程中修改原链导致的迭代异常
for res in list(chain_a.get_residues()):
    res_id = res.get_id()[1]
    if res_id > chain_a_max_resid:
        # 从链A移除残基，添加到链B
        chain_a.detach_child(res.get_id())
        chain_b.add(res)

# 保存修改后的完整结构
io = PDBIO()
io.set_structure(structure)
io.save("heterodimer_edit.pdb")

关键说明 #

• 必须创建新链B：原链所有残基属于同一个链对象，直接修改ID只会让整个链变A或变B，无法实现拆分
• 遍历残基要转成list()：遍历过程中修改原链的残基集合（移除残基）会导致迭代器异常
• 最后保存整个结构：这样才能同时保留链A和链B的所有残基

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