如何将RDKit的Conformer对象保存为坐标一致的SDF文件

问题成因 #

1. 调用AllChem.EmbedMultipleConfs(mol, numConfs=8)生成构象后，mol对象本身已经内置了8个构象，构象ID依次为0~7，你提取的conformer_i就是ID为0的构象，不需要再额外添加回mol对象。
2. 你执行mol.AddConformer(conformer_i)时，RDKit会自动为新加入的构象分配不重复的新ID（避免和已有的ID=0冲突），不会覆盖原有构象，这一步属于多余操作，反而增加了构象管理的混乱。
3. SDWriter.write()方法如果不指定第二个参数confId，默认会写入非你预期的构象，最终导致输出SDF的坐标和你打印的conformer_i坐标不匹配。

解决方法 #

你可以选择以下任意一种方式解决问题：

方式1：直接写入指定ID的构象 #

不需要额外添加构象，直接在写入时指定你要保存的构象ID即可，代码如下：

w = Chem.SDWriter('conformer.sdf')
# 直接指定写入ID为0的构象
w.write(mol, confId=0)
w.close()

方式2：单独提取构象保存 #

如果你需要单独保存该构象，避免和原分子的其他构象混淆，可以创建分子副本后写入：

# 创建分子的深拷贝
mol_single = Chem.Mol(mol)
# 清除副本中所有原有构象
mol_single.RemoveAllConformers()
# 添加目标构象
mol_single.AddConformer(conformer_i)
# 写入SDF
w = Chem.SDWriter('conformer.sdf')
w.write(mol_single)
w.close()

如果需要批量导出所有构象到同一个SDF文件，直接遍历所有构象ID写入即可：

w = Chem.SDWriter('all_conformers.sdf')
for conf_id in range(mol.GetNumConformers()):
    w.write(mol, confId=conf_id)
w.close()

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