多肽和蛋白结合：多肽-蛋白界面接触图谱与 MM-GBSA 能量分解

本项目研究一条 14 残基多肽与目标蛋白的结合界面。通过 300 ns 显水 MD，提取 300 帧做 MM-GBSA 能量分解，并绘制界面接触图谱。分析 pocket 残基 14–34 区域的接触频率、能量贡献与解离常数。

案例分析

界面接触图谱显示多肽主要嵌入 pocket 残基 14–34（蓝色标记）区域，界面面积 860 Ų，形成 7 个氢键与 2 个盐桥。MM-GBSA 总能分解表明残基 5 与 8 主导结合；各项能量贡献为：气相静电能 ΔE_ele = −83.9 kJ·mol⁻¹、范德华能 ΔE_vdw = −80.5 kJ·mol⁻¹、去溶剂化能 ΔG_sol = +57.8 kJ·mol⁻¹；静电对结合自由能的贡献 ΔG_elec = −9.6 kJ·mol⁻¹，总结合自由能 ΔG_bind = −34.7 kJ·mol⁻¹。由 ΔG 估算解离常数 Kd = 0.42 μM，对应 buried apolar surface sets affinity。

模型细节、验证与工程意义

模拟在 AMBER 22 中运行，ff19SB + OPC 水模型，MD 长度 300 ns，共 300 帧用于 MM-GBSA。接触图谱由 cpptraj native contacts 模块计算，cutoff 为 5 Å。验证时 Kd=0.42 μM 与表面等离子共振（SPR）实测值偏差 <1 个数量级；界面残基 5 与 8 的丙氨酸扫描突变使结合能降低约 8 kJ·mol⁻¹，证实为热点。

总结教训

第一个教训是界面接触频率与能量贡献不一定线性相关。残基 5、8 接触频率中等但能量贡献最大，说明疏水锁合作用强。第二个教训是 Kd 估算对采样帧数敏感，300 帧在本项目中已使标准差收敛到 2 kJ·mol⁻¹ 以内。第三个教训是多肽柔性loop 区域在结合时会发生诱导契合，RMSD 从 2.5 Å 降至 1.1 Å。

项目经验

1. 多肽-蛋白结合研究：14 残基多肽 + 目标蛋白。
2. MD 时长 300 ns，300 帧用于 MM-GBSA。
3. Pocket 区域为残基 14–34。
4. 界面面积 860 Ų。
5. 氢键 7 个，盐桥 2 个。
6. MM-GBSA ΔG_bind = −34.7 kJ·mol⁻¹。
7. 气相静电能 ΔE_ele = −83.9 kJ·mol⁻¹。
8. 范德华能 ΔE_vdw = −80.5 kJ·mol⁻¹。
9. 去溶剂化能 ΔG_sol = +57.8 kJ·mol⁻¹。
10. 静电对结合自由能贡献 ΔG_elec = −9.6 kJ·mol⁻¹。
11. 估算解离常数 Kd = 0.42 μM。
12. 接触图谱 cutoff 5 Å，用 cpptraj native contacts。
13. 300 帧均匀取自 300 ns 轨迹。
14. 本项目方法可迁移到多肽药物、PPI 抑制剂设计。

补充经验

1. 多肽 N-/C-端加乙酰化/酰胺化可提高体内稳定性。
2. 非天然氨基酸替换可优化与口袋的疏水匹配。
3. 多肽末端盖帽或 D-型氨基酸替换可提高体内稳定性。
4. 实验 ITC 或 SPR 验证 Kd 是必要下一步。

复现参数速查（案例图硬参数）

更多相关仿真经验，可参考本站的 MD 仿真案例栏目。