【摘要】 做 GROMACS 簇分析前，如何处理周期性边界、跨盒跳跃和轨迹取帧？本文给出 trjconv 命令与检查顺序。

科研人必备干货｜全套科研知识库

常用科研资料汇总，点击领取即可获取。

 

GROMACS 簇分析前，轨迹预处理主要解决三个问题：分子是否被周期性边界切开、是否存在跨盒跳跃、分析帧数是否过多。目标是让 RMSD 反映构象变化，而不是轨迹坐标问题。

 

需要准备哪些文件

基本文件是轨迹 md.xtc 和对应的 md.tpr。如果要选择特定原子组，还需要 index.ndx。这些文件应来自同一套模拟任务，以保证原子编号、拓扑和轨迹对应。

 

三步处理周期性边界

恢复完整分子

gmx trjconv -s md.tpr -f md.xtc -o md_whole.xtc -pbc whole

去除跨盒跳跃

gmx trjconv -s md.tpr -f md_whole.xtc -o md_nojump.xtc -pbc nojump

居中并选择盒子表示

gmx trjconv -s md.tpr -f md_nojump.xtc -o md_center.xtc -center -pbc mol -ur compact

交互选择时，应根据体系和分析目标选择合适的分组。处理后用 VMD 播放 md_center.xtc，检查分子是否完整、是否突然跳跃、中心位置是否符合预期。

 

取帧范围与间隔

gmx cluster 需要计算帧之间的 RMSD 关系。帧数过多会增加计算量，因此可限定分析窗口：

-b 20 -e 100 -dt 0.1

这里表示从 20 ns 到 100 ns，每 0.1 ns 读取一帧。实际时间单位应与轨迹设置保持一致。不同体系需要比较时，应固定起止时间、取帧间隔和分析原子组。

拟合组和 RMSD 计算组

拟合组用于去除整体平移和旋转，RMSD 计算组用于衡量研究对象的结构差异。两者可以相同，也可以按研究问题分别设置。默认分组不合适时，可创建索引：

gmx make_ndx -f md.tpr -o index.ndx

自定义索引时，要检查每帧原子编号和顺序一致。多分子体系还要明确比较的是单个分子、骨架、配体还是整个复合体系。

如果正在咨询科学指南针·唯理计算的理论计算相关服务，可以在沟通时先提供研究对象、轨迹文件状态、目标原子组和预期输出，使后续方案讨论更聚焦。

 

常见问题

不处理 PBC 就做 RMSD 会怎样？

如果分子被边界切开或跨盒跳跃，RMSD 可能把人为坐标变化当成构象差异，导致簇边界和代表性结构失真。具体处理方式仍需结合体系和分析对象确认。

取帧越多越准确吗？

帧数增加能提供更密集的采样，但会增加两两 RMSD 计算成本，也可能让高度相近的帧占比过高。应根据时间窗口、采样频率和研究问题取帧。

为什么不同体系要固定取帧间隔？

时间范围和取帧间隔会改变样本数量与构象出现频率。比较簇大小、簇数量或状态转换时，分析条件应尽量一致。

继续阅读：GROMACS 簇分析教程 和 gmx cluster 参数与输出文件。