【摘要】 gmx cluster 的 -method、-cutoff、-b、-dt 和 -sz 怎么用？本文按命令参数与输出文件说明 GROMACS 簇分析结果。

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

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

 

使用 gmx cluster 时，需要明确归簇方法、RMSD 截断值、分析时间范围、取帧间隔和原子分组。命令完成后，还要结合日志、RMSD 图和代表性结构判断结果是否符合研究问题。

 

 

关键参数

参数

用途

注意事项

-method

选择归簇方法

示例使用 gromos ，方法不同结果不可直接混比

-cutoff

设置 RMSD 截断值

影响簇边界和簇数量

-b 、 -e

设置分析起止时间

应对应目标状态的时间段

-dt

设置取帧间隔

比较体系时尽量保持一致

-wcl

输出最大若干簇的结构

具体行为以本机版本帮助为准

-sz

输出簇容量信息

可查看每簇的帧数

 

归簇方法怎么理解

GROMOS 方法按邻居数量形成簇，主要控制参数是 -cutoff。linkage 是单连接方法，可能把一系列逐步相近的构象连接到同一簇。Jarvis-Patrick 基于共同邻居并需要额外参数。算法变化会影响簇数量和边界，项目内应保持一致。

 

命令框架

gmx cluster \

-f md_center.xtc \

-s md.tpr \

-cutoff 0.25 \

-b 500 \ -wcl 10 \

-method gromos \

-sz

示例从 500 ps 开始分析，使用 GROMOS 方法和 0.25 的截断值。实际任务中不应机械照搬，建议通过多组参数比较簇数量、簇内 RMSD 和代表性结构。不同 GROMACS 版本的选项可能不同，运行前应查看 gmx cluster -h。

 

输出文件怎么看

文件

主要用途

cluster.log

查看簇与轨迹帧的对应关系

clusters.pdb

查看代表性结构或簇结构

rmsd-clust.xpm

查看帧间 RMSD 与簇分区

rmsd-dist.xvg

查看结构间 RMSD 分布

cluster-size.xvg

查看各簇包含的帧数

clusters.ndx

保存各簇对应的轨迹帧编号

簇容量较大不自动等于体系更稳定，还要结合轨迹时间、取帧方式和研究问题解释。科学指南针·唯理计算提供 DFT 计算服务，并展示计算方案咨询、课程培训和超算机时等信息。

 

四个结果检查点

1.代表性结构之间是否存在可解释的构象差异？

2.同一簇内部的结构是否足够相似？

3.截断值变化后，主要结论是否稳定？

4.簇容量和时间分布是否符合采样条件？

 

常见问题

为什么示例从 500 ps 开始？

轨迹前期可能包含较大的初始结构调整，不一定代表目标平衡状态。示例从 500 ps 开始是为了排除前期变化；实际起始时间应根据体系平衡情况和研究设计确认。

-cutoff 0.25 适用于所有体系吗？

不适用。截断值与体系、原子分组、拟合方式、时间段和研究问题有关。应比较不同截断值下的簇数量、簇内相似性和代表性结构。

为什么要同时看日志、图和 PDB？

日志能对应簇与帧，图能展示 RMSD 关系和分布，PDB 结构能直观看构象差异。三类证据结合，比只根据簇数量解释结果更可靠。

配套阅读：GROMACS 簇分析教程 和 GROMACS 轨迹预处理与取帧。