蛋白酶体分子模拟

项目过程：20S 蛋白酶体 500 ns 全原子动力学

客户研究蛋白酶体抑制剂，需要了解 20S 蛋白酶体在溶液中的动力学行为，特别是 β 亚基催化位点 Thr1 的可及性和周围柔性。他希望模拟能解释为什么某些抑制剂对 β5 亚基比对 β1/β2 更有效。

我用 GROMACS 对 20S 蛋白酶体（α₇β₇β₇α₇）做了 500 ns 的全原子 MD 模拟。力场 AMBER99SB-ILDN，水模型 TIP3P，温度 300 K，NPT 系综。20S 蛋白酶体加溶剂后体系规模在几十万原子量级，我用 2 fs 步长、PME 处理静电、LINCS 约束氢键。

分析指标包括：核心 Cα RMSD、回转半径 Rg、β 亚基各残基的 RMSF，以及催化 Thr1 的局部环境。

案例分析：核心稳定，loop 和 N 端前肽柔性大

第一张案例图是 500 ns 内的 Cα RMSD 和 Rg 随时间变化。蓝色实线是核心 Cα RMSD，从 0 开始上升，在约 100 ns 后稳定在 3.5 Å 左右；橙色虚线是 Rg，稳定在 61 Å 附近。这说明蛋白酶体整体桶状结构保持得很好，没有明显的解聚或塌缩。

第二张案例图是 β 亚基残基的 RMSF。大部分核心残基 RMSF 低于 0.5 Å，非常刚性；但 N 端前肽和 loop 区出现两个明显的高峰，分别达到约 1.9 Å 和 3.0 Å。催化 Thr1 用红色标出，位置靠近 N 端前肽，但本身 RMSF 不高，说明催化残基相对稳定，而周围的门控区域比较灵活。

总结教训：大体系模拟需要稳定的初始结构和长轨迹

第一，20S 蛋白酶体这种大体系，初始结构的氢原子添加很关键。我先用 pdb2gmx 检查了一遍，修正了几个缺失侧链和二硫键，否则前 50 ns RMSD 会异常漂移。

第二，500 ns 对观察门控 loop 的慢运动仍然可能不够。如果客户特别关心底物进入通道的开闭，可能需要微秒量级的模拟或增强采样方法。

第三，RMSF 分析要结合 B 因子对比。我对比了晶体结构 B 因子和模拟 RMSF，整体趋势一致，说明模拟可靠。

项目经验：柔性区域是抑制剂设计的突破口

催化 Thr1 本身刚性，但它周围的 N 端前肽和 loop 柔性大。抑制剂设计可以利用这个特点：小分子不一定要直接占据催化位点，也可以通过稳定或锁定柔性门控区来阻断底物进入。这个思路可以作为客户后续抑制剂优化的方向。

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