从TM-score到DockQ:蛋白结构预测模型评价指标全景解析与应用指南
1. 蛋白结构预测评价指标概述蛋白结构预测领域近年来取得了突破性进展AlphaFold系列模型的问世彻底改变了结构生物学的研究范式。但如何客观评价预测结果的准确性这需要一套科学、全面的评价指标体系。不同于传统实验结构解析如X射线晶体学、冷冻电镜主要依靠电子密度图质量如R因子评估计算预测模型的评价更加复杂多元。评价指标的选择取决于预测目标是单体蛋白还是复合物是全局结构还是局部构象不同指标各有侧重有的关注整体拓扑相似性如TM-score有的侧重局部原子精度如lDDT还有的专门评估相互作用界面如DockQ。理解这些指标的计算原理和适用场景对正确解读预测结果至关重要。2. 单体蛋白结构评价指标2.1 TM-score全局拓扑的金标准TM-score由张阳教授团队于2004年提出解决了传统RMSD均方根偏差的两个关键局限首先它对小距离误差赋予更大权重使其对全局折叠更敏感其次通过长度相关的归一化处理使评分不受蛋白大小影响。计算公式为TM-score max[1/L * Σ(1/(1(d_i/d0)^2))]其中L为蛋白长度d_i为对应残基距离d0是归一化因子。TM-score取值0-10.17随机相似0.5-1.0相同折叠实操中使用官方C工具计算wget https://zhanggroup.org/TM-score/TMscore.cpp g -static -O3 -ffast-math -lm -o TMscore TMscore.cpp ./TMscore model.pdb native.pdb2.2 lDDT局部原子精度的标尺局部距离差异测试lDDT是一种无需结构叠加的评估方法特别适合存在结构域运动的情况。它比较模型与参考结构中所有原子对的距离计算匹配比例。lDDT-Cα特指Cα原子的lDDT值范围0-100≥90极高置信度70-90可信50-70较低50极低Python实现可参考开源工具from biotite.structure import get_chains from biotite.structure.io import load_structure model load_structure(model.pdb)[0] native load_structure(native.pdb)[0] # 计算lDDT...3. AlphaFold特有评价指标3.1 pLDDT残基级置信度AlphaFold将预测的lDDT-Cα值称为pLDDT存储于PDB文件的B因子列。研究发现pLDDT与实验测定的lDDT-Cα高度相关Pearson r0.76。低pLDDT区域50往往对应固有无序区。3.2 pTM全局置信度预测pTM是TM-score的预测版本通过加权可能被实验解析的残基计算得出。在独立测试集上pTM与实际TM-score相关性达0.85。当pTM0.05时短链常见建议优先参考PAE或pLDDT。3.3 PAE相对位置误差预测对齐误差PAE矩阵表示若将结构在残基y处对齐残基x的预期位置误差Å。PAE特点非对称PAE(x,y)≠PAE(y,x)上限31.75Å低值表示结构域相对位置可靠4. 复合物结构评价体系4.1 DockQ蛋白-蛋白相互作用指标DockQ综合三个子指标评估复合物预测Fnat界面残基命中率LRMS长链对齐后短链RMSDiRMS界面10Å内原子RMSD评分标准0-0.23错误0.23-0.49可接受0.49-0.80中等≥0.80高质量4.2 ipTM界面特异性评分AlphaFold-Multimer引入ipTM界面pTM分别评估不同链间的相互作用。实际使用时采用pTM与ipTM的加权组合作为最终置信度。研究表明DockQ与ipTM存在函数关系二者具有一致性。5. 指标选择与应用指南5.1 根据预测目标选择指标单体蛋白TM-score全局 lDDT局部多结构域蛋白PAE分析结构域相对位置蛋白-蛋白复合物DockQ ipTM蛋白-配体复合物配体RMSD 接触分析5.2 结果解读注意事项置信度阈值pLDDT70的区域才适合功能分析矛盾情况当pTM与pLDDT冲突时建议检查PAE矩阵复合物评估需同时关注界面质量DockQ和单体准确性pLDDT动态区域低分区域可能是真实的无序状态而非预测错误5.3 实际案例解析以AlphaFold-Multimer预测的抗原-抗体复合物为例pTM0.65中等置信度ipTM0.58界面质量尚可抗原pLDDT均值82高抗体CDR区pLDDT45低这表明虽然整体折叠可靠但抗体互补决定区CDR可能需实验验证。通过生成多个预测模型并比较一致性可进一步提高可靠性。