蛋白结合位点预测

项目过程：螺旋蛋白的可药性口袋识别

客户拿到一条新解析的螺旋蛋白结构，不知道小分子该往哪里结合。他希望能先预测出潜在的配体结合口袋，并评估每个口袋的“可药性”（druggability），为后续虚拟筛选定向。

我采用了几何+能量相结合的口袋探测方法。先用 SiteFinder 或 fpocket 扫描蛋白表面凹陷，计算每个候选口袋的几何体积、疏水性、极性原子比例、以及口袋开口程度；再用 druggability 评分函数给出 0–1 之间的可药性指数。口袋倾向得分 cutoff 设为 0.6。

蛋白是一条约 250 残基的单链螺旋蛋白，没有配体共晶结构，所以预测完全基于 apo 结构。

案例分析：P1 口袋得分最高，体积 482 ų，可药性 0.83

第一张案例图是沿蛋白序列的口袋倾向得分曲线。可以看到在残基约 110 的位置有一个非常高的尖峰，标注为 pocket 1（P1），得分接近 1.0；在残基约 200 的位置有第二个峰 P2，得分约 0.7。两条虚线分别标出了口袋 cutoff 0.6 和几个关键残基 D、H、E 的位置。

第二张案例图是 P1 口袋的三维示意图。口袋位于螺旋束之间的缝隙中，估计体积 482 ų，可药性评分 0.83。这个数值通常被认为是高可药性，适合小分子药物设计。图中用颜色深浅表示口袋深度，红色区域是口袋最深、最容易形成强相互作用的位置。

总结教训：apo 结构的口袋可能和 holo 结构不同

第一，apo 结构预测出的口袋有时候会在配体结合后发生诱导契合而变大或变深。我一般会建议客户如果后续有共晶结构，把 holo 结构拿来再验证一次。

第二，口袋体积不是越大越好。482 ų 对一个小分子药物来说比较合适；如果体积超过 1000 ų，往往是蛋白表面浅沟，不适合高亲和力小分子。

第三，druggability 评分要结合序列保守性一起看。P1 如果同时位于高保守区，那它很可能是功能位点，结合小分子可能带来生物学效应。

项目经验：口袋预测是筛选漏斗的第一步

这个项目的输出不是一个最终配体，而是一个优先口袋列表。P1 应该作为虚拟筛选的首要靶点，P2 可以作为备选或别构位点。客户拿到结果后，可以优先对 P1 做分子对接和片段筛选，比盲目全蛋白扫描高效得多。

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