ARTICLE DETAIL

资讯详情

深耕编程入门与网站建设的一线实战洞察。

AlphaFold 置信度完全指南:pLDDT 与 PAE 怎么读才靠谱

AlphaFold 置信度完全指南:pLDDT 与 PAE 怎么读才靠谱 AlphaFold 置信度完全指南pLDDT 与 PAE 怎么读才靠谱【免费下载链接】alphafoldOpen source code for AlphaFold 2.项目地址: https://gitcode.com/GitHub_Trending/al/alphafold拿到一份 AlphaFold 预测结果你很难第一眼分辨哪些结构可信、哪些只是噪声。这篇文章教你把 pLDDT 和 PAE 读明白建立一套从快速判断到深度分析的蛋白质结构预测置信度评估流程。一、3 分钟速览拿到预测结果后先看什么别急着把预测图拖进分子查看器。先按这个顺序过一遍先看 pLDDT 条形图的均值和分布→再看 PAE 矩阵的块状结构→最后检查两者是否打架。均值 70 直接存疑均值 ≥ 70 且条形图平坦才有资格往下读。PAE 矩阵用来定位哪里可能塌。两个指标不一致时才是真正需要动脑的时刻。二、pLDDT它到底在回答什么问题 它回答的问题是第 i 个残基的局部结构你该信几分拿两个具体场景对比某残基 pLDDT 95位于稳定 α 螺旋核心区Cα 原子位置预测误差通常在 1 Å 内骨架取向可以直接用于对接、突变分析。同一链另一残基 pLDDT 40落在 N 端或 loop 区主链方向大概率对但原子坐标已经不能拿来当真实位置用只能当拓扑参考。四档分级怎么落地用分数区间颜色含义你该怎么用90–100深蓝非常高原子级可信可直接做对接、位点设计70–90浅蓝中等偏上骨架可信侧链取向建议复核50–70浅绿偏低拓扑可参考别信坐标0–50橙色无序 / 不可信视为柔性或无序区勿当刚性处理为什么有的残基 pLDDT 天然偏低三个高频原因内在无序区IDR蛋白本来就没定型模型只能给一个低分。柔性 loop结构存在但自由度大模型主动压低置信度。MSA 覆盖不足该区域缺少同源序列支撑信号本身就弱。def compute_plddt(logits): # logits: [N_res, N_bins] 模型原始输出 num_bins logits.shape[-1] bin_centers np.arange(0.5/num_bins, 1.0, 1.0/num_bins) probs scipy.special.softmax(logits, axis-1) return (probs * bin_centers[None]).sum(-1) * 100 # 期望值 × 100 → 0~100 分三、PAE残基对之间的距离误差地图 它回答的问题是若把预测结构和真实结构做最优叠合残基 i 和 j 之间的距离误差大约是多少PAE 是 N×N 矩阵对角线附近 近邻残基对远离对角线 远距离残基对。看三种典型模式多亚基复合物的 3 句话读法 对角块亚基内 PAE低 → 单链骨架可信非对角块亚基间 PAE低 → 界面相对取向可信某非对角块整体高值 → 该亚基对相对位置不确定界面结论要谨慎。def compute_pae(logits, breaks): probs scipy.special.softmax(logits, axis-1) # [N_res, N_res, N_bins] bin_centers breaks (breaks[1] - breaks[0]) / 2.0 # bin 中心 bin_centers np.append(bin_centers, bin_centers[-1] (breaks[1] - breaks[0])) return (probs * bin_centers[None, None]).sum(-1) # 期望距离误差 [Å]四、组合判读pLDDT 和 PAE 打架时怎么办 ⚠️冲突场景最可能的解释下一步该做什么pLDDT 均值高但某残基对 PAE 极高该残基对相对取向建模失败可能是 loop 或界面单独看该残基 pLDDT若也低确认是柔性区而非 bugpLDDT 均值低但 PAE 对角块整体低全局拓扑对但局部原子级细节未定型骨架方向可参考原子坐标不可信考虑提高 MSA 深度复跑pLDDT 呈高-低-高块状PAE 出现对应大方块典型多结构域蛋白低 pLDDT 区是结构域连接 loop把低 pLDDT 区当柔性连接子PAE 方块边界的残基号即结构域边界复合物某非对角块整体高值对应亚基 pLDDT 正常两个亚基骨架各自 OK但相对取向不确定换更多种子复跑若仍高用其他实验数据如 FRET约束界面底层逻辑pLDDT 管单链局部精度PAE 管残基对相对可靠性。一个说原子级一个说相对级两者一致才算真的信。五、技术原理速览给好奇的你如果你不关心实现细节直接跳到下一节。pLDDT 的计算链路PAE 的数学表达一行$$\text{PAE}{ij} \sum{k1}^{K} p_{ij,k} \cdot c_k$$三个置信度头在模型里的分工头名称输出形状作用PredictedLDDTHead[N_res, N_bins]每残基 pLDDT 的 softmax 分布PredictedAlignedErrorHead[N_res, N_res, N_bins]每残基对 PAE 的 softmax 分布ExperimentallyResolvedHead[N_res, N_atoms]每个原子实验上是否可解析的概率六、实战检查清单 常见问题 ✅从预测结果到写进论文的完整操作链先扫 pLDDT 条形图均值 ≥ 70 且无大面积橙色打开 PAE 矩阵画圈标出所有异常方块。用 pLDDT 逐残基标注低置信区 70记录残基号区间。复合物分别核对每个非对角块确认界面相对定位。对存疑区域换 seed / 提高 MSA 深度复跑 1 次看 pLDDT 是否稳定。若需交叉验证用 ESMFold 或 RoseTTAFold 跑同序列看低置信区是否一致。论文里报告时同时给出 pLDDT 均值 PAE 热图缺一不可。所有写入结论的结构其 pLDDT 必须 ≥ 70 且对应 PAE 块低值。高频 FAQ问题现象最可能的原因你该做什么整体 pLDDT 50序列内在无序或 MSA 覆盖极弱用 IUPred 等工具先判无序别硬信结构PAE 出现对角线外高值大方块对应残基对相对取向建模失败定位方块边界的残基号结合 pLDDT 确认区域pLDDT 均值正常但某残基对 PAE 极高局部构象不确定可能是 loop / 界面单独看该残基 pLDDT若也低确认柔性区复合物界面非对角块高值亚基间相对取向不确定换种子复跑用实验数据如 FRET约束界面pLDDT 稳定但跨 run 变化大随机种子 / MSA 深度不稳定固定种子并记录参数多次复跑取中位数想再深挖一步用 ESMFold 或 RoseTTAFold 做交叉验证是判断模型特有偏差还是蛋白本来就不确定的最快路径。下次拿到 AlphaFold 结果按这个顺序看pLDDT 均值定基调PAE 方块定区域两者冲突时再动手——不会漏。【免费下载链接】alphafoldOpen source code for AlphaFold 2.项目地址: https://gitcode.com/GitHub_Trending/al/alphafold创作声明:本文部分内容由AI辅助生成(AIGC),仅供参考
返回列表