ChimeraX中如何修复缺失的侧链原子？

Question

在使用ChimeraX进行蛋白质结构可视化与分析时，常遇到晶体结构或冷冻电镜密度图中部分氨基酸侧链原子缺失的问题。这会影响后续的分子对接、结构修饰或相互作用分析。用户常问：如何利用ChimeraX自动修复缺失的侧链原子？特别是当PDB结构中存在不完整的赖氨酸、精氨酸或组氨酸等带复杂侧链的残基时，能否通过内置工具如“Rotamers”或“AddH”功能补全原子坐标，并确保其符合合理的空间构象和氢键网络？此外，如何结合实验密度图或能量优化判断所补原子的合理性？这是结构生物学预处理流程中的关键步骤。

Answer 1

利用ChimeraX修复蛋白质结构中缺失的侧链原子：从基础操作到高级验证

1. 问题背景与常见场景分析

在结构生物学研究中，PDB数据库中的晶体结构或冷冻电镜（cryo-EM）密度图常存在氨基酸侧链不完整的问题，尤其是赖氨酸（Lys）、精氨酸（Arg）和组氨酸（His）等具有长且可电离侧链的残基。这些缺失直接影响分子对接、氢键网络分析以及后续的分子动力学模拟。

常见的缺失类型包括：

• 末端氮/氧原子缺失（如Lys的NZ原子）
• 侧链χ角未定义导致构象错误
• 质子化状态不确定（如His的δ vs ε质子化）
• 缺乏氢原子，影响静电相互作用建模

2. ChimeraX内置工具概览

ChimeraX提供多个模块用于结构修复与优化，主要包括：

工具名称	功能描述	适用阶段
AddH	添加氢原子，基于pKa预测质子化状态	预处理
Rotamers	替换侧链为常见rotamer构象	修复缺失
Minimize	局部能量优化，减少空间冲突	后处理
Fit in Map	将原子拟合至cryo-EM密度图	验证合理性

3. 分步操作流程：自动补全侧链原子

1. 启动ChimeraX并加载目标结构：open 7abc
2. 识别缺失侧链：select ligand or not protein 排除配体，聚焦蛋白主链
3. 使用“Rotamers”工具补全侧链：
   • 菜单路径：Tools → Structure Editing → Rotamers
   • 选择目标残基（如Lys120），点击“Replace”使用最可能rotamer
4. 执行AddH添加氢原子：addh selection 或通过界面调用
5. 设定pH值以影响质子化状态：addh pH 7.4
6. 对组氨酸特别处理：使用swapaa his命令指定HISD（δ质子化）或HISE（ε质子化）
7. 运行能量最小化：minimize steps 100 interval 10
8. 加载对应密度图进行验证：open emd_1234.map
9. 拟合侧链至密度：fitmap #1 inMap #2
10. 评估拟合质量：查看相关系数（CC）是否 > 0.8

4. 结合实验密度图判断合理性

当存在高质量cryo-EM密度图时，可通过空间一致性验证补全原子的可信度。关键步骤如下：

# 示例脚本：自动化侧链修复与密度验证
from chimerax.core.commands import run

run(session, "open 7abc.pdb")
run(session, "open emd_1234.map")
run(session, "sel /A:LYS@CA,CB,CG,CD")
run(session, "rotamers sel")
run(session, "addh sel pH 7.4")
run(session, "minimize selectedOnly true")
run(session, "fitmap #1 inMap #2")
    

5. 能量优化与构象合理性评估

仅依赖rotamer数据库可能忽略局部环境冲突。建议结合以下策略提升可靠性：

• 使用AMBER力场进行短程能量最小化
• 检查范德华冲突（Clash Score）
• 分析氢键网络变化（Tools → Structure Analysis → H-bonds）
• 对比修复前后RMSD（根均方偏差）

6. 高级技巧：批处理与脚本化修复

对于含多个缺失