酶分子动力学模拟

我做过一个酶-底物复合物的分子动力学项目，客户想理解底物结合后酶活性位点周边残基的柔性变化，为后续突变设计提供依据。我这个项目的关键图就是左侧的酶-底物复合物示意图和右侧的残基柔性 RMSF 曲线。

左侧图展示了酶活性中心的三残基催化 triad：His、Asp、Ser。底物分子位于三个残基围成的口袋里，用橙色棍状表示。我做 MD 之前先通过分子对接把底物放到活性位点，然后做 100 ns 显式溶剂 MD。对接时用 AutoDock Vina，格盒中心放在催化 Ser 附近，搜索空间 25×25×25 Å。

右侧图是每个残基的 RMSF（根均方涨落）。整体看，大部分残基 RMSF 在 0.8-1.3 Å 之间波动，但在残基索引约 210 的位置出现一个尖锐高峰，RMSF 接近 1.6 Å。这个位置是一个 loop 区，远离活性中心。高峰说明这段 loop 在底物结合后变得特别柔，可能是诱导契合机制的一部分：loop 运动帮助底物进入并稳定过渡态。客户后续正是针对这个 loop 做了丙氨酸扫描突变，验证了我的预测。

模拟流程：蛋白用 AMBER ff19SB，底物和关键水分子用 GAFF2， RESP 电荷，TP3P 水盒子，12 Å 缓冲。先最小化、升温、NPT 平衡 2 ns，然后 NVT 生产 100 ns。RMSF 分析去掉前 20 ns，取平衡后的 80 ns。我还做了动态互相关矩阵 DCCM 和主成分分析 PCA，发现前两个主成分主要对应这个柔性 loop 的开合运动。

踩坑经验：第一个是底物质子化状态。酶催化通常涉及质子转移，如果底物的质子化态设错，整个 MD 都会跑偏。我用 PROPKA 和 H++ 做了预测，并结合文献选了最可能的状态。第二个是结合模式：对接可能给出多个构象，我选了催化 Ser 羟基氧到底物羰基碳距离最近的 pose，这样更接近反应起始态。第三个是金属离子：这个酶不含金属，如果有 Zn²⁺/Mg²⁺，力场参数要单独处理。

工程经验：酶 MD 不要只给一张 RMSF 图，要把柔性区域和已知功能位点、实验突变数据对照。如果 RMSF 高峰正好在活性口袋入口或变构位点，就非常有价值；如果高峰在远离功能区的尾巴，可能只是本身无序。我这个项目我还测了底物与催化 triad 的距离分布，确认 His-Ser 氢键和 Ser-底物亲核攻击距离在模拟中都保持稳定。

更多酶分子动力学和分子对接联合分析的内容，我汇总在 [酶分子动力学模拟](https://www.keyanxueshu.com/category/md/) 栏目里。[酶分子动力学模拟](https://www.keyanxueshu.com/)

除了 RMSF，酶-底物项目我还会补做一个 MM-GBSA 结合自由能分解。虽然 MD 本身不直接给结合能，但把轨迹喂给 gmx_MMPBSA 或 Amber 的 MMPBSA.py，可以得到每个残基对结合的贡献。通常催化 triad 的 His 和 Ser 会贡献较大的极性相互作用，而口袋入口的疏水残基贡献疏水作用。把 RMSF 高峰残基和 MM-GBSA 能量贡献对照，能判断柔性 loop 是“功能性柔性”还是“无关柔性”。如果高 RMSF 残基同时有较大的负结合能贡献，那它就是诱导契合的关键；如果能量贡献接近零，可能只是背景噪声。这种联合分析比单看 RMSF 更有说服力。