为什么在PyMOL中使用cealign进行结构比对时提示target选择过短

报错原因 #

CEalign-Error: Your target selection is too short. 报错的核心触发逻辑是CEalign算法无法在你指定的目标选择集中找到≥8个可用于对齐的Cα原子，除了原子无坐标的问题外，更常见的触发原因是未限定有效肽链范围、多链结构选择逻辑异常、结构未完全加载就执行对齐命令。

解决方案 #

按以下步骤操作即可解决问题：

• 第一步：加载结构时禁用异步加载，确保结构完整加载后再执行后续操作

fetch 2mi1, async=0
fetch 3zkt, async=0

• 第二步：执行对齐前限定选择范围为蛋白质的有效肽链Cα原子，同时指定需要对齐的链（如果不清楚结构的链信息，可先运行chains命令查看所有已加载结构的链标识）

# 示例为对齐2mi1的A链和3zkt的A链，可根据实际需求替换链标识
cealign 2mi1 and name CA and chain A, 3zkt and name CA and chain A

异常排查 #

如果操作后仍有报错，按以下规则排查：

• 检查待对齐的两条肽链长度，单条残基数少于8的结构无法使用CEalign对齐，可替换为align或super命令执行普通结构对齐
• 清除结构的冗余构象，避免多构象导致的坐标识别错误

remove (not alt ''+A)

• 若需要做全结构所有蛋白质链的对齐，先合并选择集再执行对齐

select mobile, 2mi1 and polymer.protein and name CA
select target, 3zkt and polymer.protein and name CA
cealign mobile, target

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