1. PAE 概述:定义、意义与应用场景

PAE,全称是 Predicted Aligned Error。这东西说白了,就是给蛋白质模型打的“质量分”。

但我说的可不是那种笼统的好或不好。它告诉你的是:模型里任意两个残基之间的位置可靠性。你想想看,一个蛋白几百个残基,有的区域可能算得准,有的区域可能完全是瞎猜。PAE 干的就是这个——把这种不确定性量化出来。

1.1 PAE 的定义

PAE 是一个二维矩阵。横轴和纵轴都是残基编号,从 1 到 N。矩阵中每个值表示:如果你把第 i 个残基的位置固定住,那么第 j 个残基的位置误差期望值,单位是 Ångström(埃)。

怎么理解?举个例子:

PAE[i][j] = 5.2 Å

意思是,当 i 这个残基不动时,j 残基预测位置和真实位置之间的误差,平均在 5.2 埃左右。值越小,说明这两个残基的相对位置越可信。

我个人习惯:拿到 PAE 矩阵,先看对角线附近的值。自比对区域通常 PAE 很低(<2 Å),说明局部结构稳定。如果某个区域的对角线出现高值,那这段骨架大概率是 disordered。

1.2 为什么 PAE 比整体分数更「有用」?

整体指标比如 pLDDT 或 TM-score,它们给你的是一个「总分」。但蛋白质结构的实际应用,往往只关心特定区域:

  • 活性位点是不是靠谱?
  • 两个 domain 之间的相对位置对不对?
  • 某个 loop 会不会影响结合界面?

整体高分不代表局部可靠。我在项目中遇到过多次:pLDDT 整体 >85,但 PAE 矩阵显示 domain 间的旋转方向完全不确定。这个结构要是拿去对接,结果肯定是错的。

场景 整体高分就可以? 更依赖的指标
核心折叠验证 基本够用 pLDDT
结构域划分 不够 PAE 矩阵 + 聚类
蛋白-蛋白对接 完全不够 PAE (间域相对位置)
口袋识别 偶尔够 pLDDT + PAE 局部过滤
我曾经犯过的错:接手一个 AlphaFold2 预测结果,pLDDT 全局 >90,我直接拿去做分子替换。结果结构解析死活跑不动。回去一查 PAE,发现两个 domain 之间的 PAE 值高达 12 Å。原来 domain 间相对排布完全不可信。从那以后,我给自己定了个规矩:不看 PAE 矩阵,绝对不乱用结构。

1.3 PAE 矩阵怎么「读」?

PAE 矩阵通常用热图表示。颜色从蓝色(低误差)到红色(高误差)。读图的核心思路是看模式:

  1. 对角线颜色深蓝:每个残基和它自己的相对位置是确定的——这是 trivial 的。
  2. 局部方块颜色浅蓝或绿:说明这个区域内残基之间的相对位置可靠。通常对应稳定的二级结构元素。
  3. 远距离区域出现红色或橙色:两个 domain 之间相对位置不确定。这是多 domain 蛋白的典型特征。
  4. 整行/整列高值:某个残基和其他所有残基的 PAE 都高。这个残基可能是柔性 linker 或者末端尾巴。

嗯,这里要注意:PAE 矩阵反映的是「相对误差」,不是「绝对误差」。两个 domain 如果 PAE 高,不代表它们各自折叠错了,只说明它们的空间排布不确定性高。

1.4 应用场景:你什么时候需要 PAE?

  • 多 domain 蛋白结构分析:判断 domain 间相对位置的可信度。这是 PAE 最经典的应用。
  • 结构域划分:用 PAE 矩阵做层次聚类,自动识别结构域边界。
  • 对接前的筛选:蛋白-蛋白复合物预测时,用 PAE 过滤掉界面不可信的模型。
  • 实验结构验证:把 PAE 高的区域标记为「不敢直接用」,实验设计时重点避开。
  • Loop 建模可靠性:某个 loop 的 PAE 值如果是 8 Å 以上,那它大概率是盲目猜的。
我的一个习惯操作:拿到一个新蛋白的结构预测结果后,先做两个操作:一是看 pLDDT 直方图,二是看 PAE 矩阵的 domain 块是否清晰。如果 PAE 矩阵里没有明显的低 PAE 方块,那这个蛋白大概率是 intrinsically disordered 或者多构象的。这种情况不要强行用一个状态去表示。

1.5 PAE 和 TM-score 之间的关系

简单说:TM-score 衡量「整体形状相似度」,PAE 衡量「内部相对位置的确定性」。它们是互补的。TM-score 高不一定 PAE 好,PAE 好也不一定能得到高 TM-score。两者结合使用才能做出合理判断。

PAE 与 TM-score 知识体系结构 蛋白质结构可靠性评估 PAE(局部间域相对可靠性) TM-score(全局形状相似度) 定义:残基间位置误差 读图:热图的模式 定义:结构比对分值 阈值划分等级 多domain结构分析 结构域自动划分 对接模型筛选 实验验证参考 校验结论:通过 PAE 矩阵 + TM-score 联合判定 核心思维:整体与局部不可偏废,PAE 补足 TM-score 的「区域盲区」

你经常需要判断两个 domain 之间的方向是不是对的。PAE 矩阵如果显示出清晰的低 PAE 区域,那说明这个排布值得信赖。如果整个矩阵一片橙红,那最好别拿着这个结构去做分子替换或对接——大概率是幻觉。

避坑指南:我曾经在分析一个 800 残基的受体蛋白时,发现 TM-score 高得离谱——0.95。当时想这个预测真完美。结果细看 PAE 矩阵,N 端 200 残基和 C 端 200 残基之间的 PAE 全红,12 Å 以上。也就是说,两个 domain 的相对方向是瞎蒙的。用这个构象去做对接,前 500 个 pose 全部集中在错误界面上。所以,TM-score 再高也得看过 PAE 再下结论。