3. 结构质量评估与选择决策

在蛋白靶点结构准备流程中,拿到一个 PDB 条目只是起点,真正决定后续对接、口袋分析或分子动力学是否可信的,是结构质量评估与多结构取舍。同一个靶点往往有几十甚至上百个实验结构,它们的分辨率、R 因子、缺失残基、突变和构象状态各不相同。选错结构,后续所有计算都建立在不可靠的坐标上;选对结构,则能在源头控制误差。本节围绕“看哪些指标、怎么读这些指标、多个结构如何取舍”展开。

一、先分清两类质量指标

结构质量指标可以分成两类。第一类是全局质量指标,描述整个模型与实验数据的吻合程度,包括分辨率、R 因子(Rwork 与 Rfree)、以及整体 B 因子。第二类是局部质量指标,描述特定区域是否可信,包括电子密度质量、B 因子分布、缺失残基、突变、配体占据与构象状态。全局指标决定“这个结构整体能不能用”,局部指标决定“我要用的那个口袋能不能用”。

分辨率是衍射数据本身的性质,不是模型的性质。它反映电子密度图能分辨的最小细节:

  • ≤1.5 Å:通常可分辨单个原子位置,水分子和部分氢原子位置可辨,侧链构象可靠。
  • 1.5–2.0 Å:主链和大多数侧链清晰,是药物设计常用的高质量区间。
  • 2.0–2.5 Å:主链可靠,部分表面侧链构象不确定,配体位置基本可信。
  • 2.5–3.0 Å:只能看到主链走向和大致侧链,配体构象需谨慎。
  • >3.0 Å:主链可追踪,侧链和配体细节不可靠,一般不用于精细口袋分析。

R 因子衡量模型计算结构与实验衍射数据的偏差。Rwork 用参与精修的数据计算,Rfree 用未参与精修的 5%–10% 数据计算,因此 Rfree 更能反映过拟合。经验上,高分辨率结构 Rfree 常在 0.20 以下,2.5 Å 左右结构在 0.25 以下可接受。若 Rwork 很低而 Rfree 明显偏高(差值大于 0.07),说明模型可能过拟合,需要警惕。

二、局部指标:真正决定口袋可用性

电子密度是判断局部坐标可信度的直接证据。在 2Fo−Fc 图上,主链和侧链应被密度包裹;若某段侧链没有密度支撑,说明该构象是建模推测。配体尤其要看 2Fo−Fc 与 Fo−Fc 图:配体应完整落在正密度中,且 Fo−Fc 图上无显著残余正峰或负峰。若配体密度残缺,其结合模式不可信。

B 因子(温度因子)反映原子位置的柔性或不确定性。同一结构内,B 因子显著高于全局平均值的区域通常更不可靠。口袋残基若 B 因子普遍偏高,说明该口袋柔性大,对接结果对构象敏感。需要注意:B 因子受晶体堆积和精修策略影响,不同结构之间不能直接比较绝对值,应看同一结构内的相对分布。

缺失残基是另一个关键点。PDB 文件中常出现残基编号跳跃,表示这些残基在电子密度中不可见、未被建模。若缺失残基位于结合口袋内,口袋形状就是不完整的,对接或口袋体积计算会失真。此外,突变(如工程化突变、同源替换)会改变侧链甚至主链构象,必须确认目标口袋是否被突变影响。

三、多结构取舍策略

面对多个结构,建议按以下优先级决策:

  1. 先筛分辨率与 Rfree:优先保留分辨率 ≤2.5 Å 且 Rfree ≤0.28 的结构;若候选很少,可放宽到 3.0 Å,但需在报告中标注局限。
  2. 再看口袋完整性:检查目标口袋是否有缺失残基、突变或低密度区域。口袋完整、密度良好的结构优先。
  3. 比较构象状态:同一靶点可能处于活性构象、非活性构象或不同寡聚状态。应根据研究目的选择——做抑制剂设计通常选与抑制剂共结晶的构象,做变构调节则选对应变构口袋开放的结构。
  4. 看配体信息:若结构含共结晶配体,可将其作为口袋定位和对接验证的参照;无配体结构需借助口袋检测工具推断。
  5. 最后看 B 因子与寡聚状态:口袋区域 B 因子低、寡聚状态与生理相关者优先。

下面这张图把上述指标组织成一条从“候选结构”到“可用结构”的筛选决策链,并标出每个判断节点使用的具体指标与典型阈值。

多结构质量评估与取舍决策链 候选结构集合(同一靶点多个 PDB) 判断 1:分辨率 ≤2.5 Å 优先;2.5–3.0 Å 谨慎;>3.0 Å 一般排除 判断 2:R 因子 R_free ≤0.28;R_free − R_work >0.07 警惕过拟合 判断 3:目标口袋局部质量 电子密度完整?缺失残基是否落在口袋内? 是否存在影响口袋的突变?配体密度是否完整? 判断 4:构象状态与寡聚形式 活性/非活性构象?是否含目标配体? 寡聚状态是否与生理相关? 选定结构 → 进入结构准备 记录选择依据与已知局限 分辨率不足 → 排除 或仅作低精度参考 R_free 过高 → 降级 检查过拟合与建模误差 口袋缺失/突变 → 排除 或改用同源结构补全 构象不匹配 → 换结构 按研究目的重选构象 评估指标速查 ● 分辨率:数据细节上限 ● R_work / R_free:拟合与过拟合 ● 电子密度:局部坐标可信度 ● B 因子:柔性/不确定性 ● 缺失残基:口袋完整性 ● 突变:侧链/主链改变 ● 构象状态:活性/非活性 ● 配体占据:口袋定位参照 ● 寡聚状态:生理相关性 注:B 因子仅在同一结构内 比较相对分布,不跨结构比绝对值

四、一个可复用的选择记录表

为了让取舍过程可追溯,建议为每个候选结构填写一张简表。以下为示例字段,实际使用时按项目需要增删:

PDB ID 分辨率 (Å) R_free 口袋缺失残基 关键突变 构象/配体 结论
1ABC 1.8 0.21 无 无 活性构象,含抑制剂 首选
2DEF 2.4 0.26 Leu85–Ser88 缺失 无 非活性构象,无配体 备选(口袋不完整)
3GHI 3.1 0.31 无 Tyr102Ala 活性构象,含底物类似物 排除(分辨率低且突变影响口袋)

表中数据为教学示例,用于说明判断逻辑,并非真实实测结果。实际项目中,缺失残基和突变信息可从 PDB 文件 SEQRES 与 ATOM 记录比对、以及结构注释中获得;电子密度和 B 因子可用 PyMOL、Coot 或 CCP4 系列工具查看。

实用提示:当高质量结构只有一个时,不要因为“没有选择”就跳过评估。应在报告中明确写出该结构的已知局限(如分辨率、缺失残基、突变),并在后续对接或动力学中通过柔性处理、多构象采样等方式补偿,而不是默认结构完美。
常见误区:只看分辨率选结构。分辨率高不代表口袋区域一定可靠——局部 B 因子高、侧链密度差、配体建模错误都可能出现在高分辨率结构中。全局指标用于初筛,局部指标用于终选。

总结来说,结构质量评估的核心逻辑是:全局指标筛掉明显不可用的结构,局部指标确认目标口袋是否可信,构象与配体信息决定结构是否匹配研究目的。把这三层判断固化为可追溯的记录,后续的结构准备、口袋分析和对接结果才有可靠的基础。