蛋白分子模拟对接

项目过程：蛋白-蛋白对接与能量漏斗分析

客户需要预测两个已知结构的蛋白如何形成复合物，没有实验复合物结构。他的目标是得到一个近天然（native-like）的对接构象，并用量化指标说明这个构象的可信度。

我用的是蛋白-蛋白对接流程：先用 ZDOCK 快速生成 9000 个 decoy，再用 ZRANK 对这些 decoy 重新打分，最后按界面 RMSD 聚类，分析每个 cluster 的能量和几何质量。

打分指标包括 ZRANK 能量、界面 RMSD（iRMSD）、天然接触比例 F_nat，以及 CAPRI 质量标准 Q。

案例分析：cluster 1 近天然，iRMSD 1.9 Å，F_nat=0.82

第一张案例图是 900 个 decoy 的 ZRANK 能量与界面 RMSD 散点图。横坐标 iRMSD 越低表示越接近天然结构，纵坐标 ZRANK 能量越低表示打分越好。左下角绿色阴影区是 CAPRI acceptable 区域。可以看到大量 decoy 聚集在 iRMSD 10–15 Å、能量 −50 到 −70 kcal·mol⁻¹ 的区域；而真正的近天然构象用黑色五角星标出，位于 cluster 1，iRMSD 很低。

第二张案例图是聚类排名。cluster 1 包含约 220 个 decoy，F_nat=0.82，iRMSD=1.9 Å，Q 评级达到 medium/high。随着 cluster rank 增加，F_nat 迅速下降，说明排名第一的 cluster 明显优于其他。

总结教训：对接打分必须配合聚类才有意义

第一，只看最低能量构象很容易选到错误的最小值。很多打分最低的 decoy 其实是 iRMSD 很大的非天然构象。聚类后取最大 cluster 的中心，通常更接近天然结构。

第二，F_nat 和 iRMSD 是判断对接质量的金标准。F_nat>0.5 且 iRMSD<5 Å 一般认为是可接受的，这里 cluster 1 的 0.82 和 1.9 Å 已经属于高质量预测。

第三，蛋白柔性是蛋白-蛋白对接的最大挑战。ZDOCK 把蛋白当刚性体，如果结合伴随明显构象变化，需要结合 MD 或 HADDOCK 的柔性对接来修正。

项目经验：能量漏斗是说服审稿人的利器

把 900 个 decoy 画成能量-几何图，天然的构象会形成一个明显的“漏斗”——能量低且 iRMSD 小。这个图比单纯文字更有说服力。客户如果要把结果放到论文里，这张图是核心证据。

这篇文章把我在这个项目里用到的计算流程、参数取舍和图上读数做了完整复盘。如果你也在做类似方向，可以参考我的思路：蛋白分子模拟对接。