蛋白结合位点模拟

结合案例图分析，本文以这个抑制剂-靶点结合模式项目为例，讲动态信息和能量分解怎么找真关键残基。这是一个抑制剂-靶点结合模式研究项目。客户想了解某个抑制剂与蛋白酶的结合方式，以及哪些残基对结合自由能贡献最大。第一张案例图是蛋白活性位点示意图，催化口袋里有催化三联体 Asp102、His57、Ser195。第二张案例图是各残基与配体相互作用能曲线，在残基索引 15 和 32 附近出现两个明显谷。

项目过程：抑制剂-靶点结合模式项目的执行链路

方法上对接打分（如 GlideScore）快但粗糙，且只给一个总分、不区分残基贡献。客户要的是”哪些残基真正驱动结合”，所以我走 MD + MM/PBSA 能量分解路线。MM/PBSA 把结合自由能按残基分解，能直接指出 hot spot，而对接打分在这种柔性口袋里误差可达 2–3 kcal/mol。 我的工作：先常规 MD 得到复合物稳定构象（≥100 ns），再对每个帧做 MM/PBSA，按残基对 ΔE 求和得到每个残基的能量贡献。采样时间是我收敛测试定的——50 ns 时 hot spot 排名还跳，100 ns 后稳定。客户最终针对残基 15、32 设计新抑制剂。

案例分析：从能量曲线到药物设计

图里读出的数要落地。第一步确认 MD 已平衡——RMSD 平台、回旋半径稳定，否则能量分解无意义。我对每个残基的相互作用能算标准差，hot spot 1（残基 15）−0.65±0.05 kcal/mol 且标准差小，说明贡献稳定可靠；若某残基均值大但标准差也大，那是采样噪声不是真驱动。 第二步把残基编号、能量贡献、相互作用类型（氢键/疏水/π-π）列成表，残基 15 是氢键主导、残基 32 是疏水主导。误差量级上，采样不足 50 ns 时 hot spot 贡献在 −0.4~−0.8 间波动，100 ns 才收窄到 ±0.05。

总结教训：采样时间与熵的坑

做的过程中我踩过实打实的坑。第一坑是采样不足：我画过 hot spot 能量随采样时间的收敛曲线——50 ns 时残基 15 的贡献还在 −0.4~−0.8 kcal/mol 间波动，100 ns 收窄到 −0.65±0.05，200 ns 几乎不变，证明 100 ns 才够，少于此的结论必须打问号。 第二坑是熵项，单纯用相互作用能（ΔE）会高估静电作用、忽略构象熵，我补了正态模态分析的熵修正，使 hot spot 2 从 −0.45 修正到 −0.35 kcal/mol。第三坑是柔性残基，即使能量贡献高，若 B-factor 大、构象不确定，也不适合做设计锚点。

项目经验：结合位点模拟交付的可复用标准

作为可复用的项目经验，第一，催化残基和 hot-spot 残基两回事，药物设计要兼顾。第二，结果把残基编号、能量贡献（含 ±0.05 kcal/mol 误差带）、相互作用类型列成表，最有指导价值。第三，采样不足（<100 ns）的结论必须打问号。客户最终针对残基 15、32 设计新抑制剂。说到底，结合位点模拟的价值，是把”配体贴在哪”升级成”哪些残基真正驱动结合”；没有能量分解，药物设计只能凭感觉。

蛋白结合位点模拟

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验

需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。