AMBER 簇分析：从轨迹提取 RMSD 拟合的 PDB 文件

第一步：准备目录

簇分析过程会创建大量文件，即轨迹中每一帧的 pdb 文件。为了避免这些文件 '污染' 我们的工作目录，我们需要创建一些子目录来进行这个过程：

mkdir clustering
cd clustering
mkdir PDBfit

第二步：使用 ptraj 从轨迹中提取 RMSD 拟合的 PDB 文件

首先，我们需要使用 ptraj 命令从轨迹文件中提取 RMSD 拟合的 PDB 文件。以下是一个示例 ptraj 命令文件（extract_pdb.ptraj）：

#orient all frames best fit to backbone of helix in NMR structurerms first mass :2-20@CA,C,N

run

这个 ptraj 命令文件将执行以下操作：

1. 读取轨迹文件 `../equil_1-10.binpos`。
2. 将所有 PDB 帧输出到 `PDBfit` 子目录下的 `trp.pdb` 文件中。
3. 对每个帧进行 RMSD 拟合，将轨迹中的所有帧与 NMR 结构的第 2 到 20 号残基的 backbone（CA, C, N）进行最佳拟合。

运行这个 ptraj 命令文件后，`PDBfit` 目录中将会生成 50,000 个 pdb 文件，每个文件代表一个已经被 RMSD 拟合到轨迹中 NMR 结构的第 2 到 20 号残基的单个快照。

**提示：**您可以根据您的具体需要修改 ptraj 命令文件中的参数，例如：

* `trajin` 参数指定了轨迹文件的路径。
* `trajout` 参数指定了输出 PDB 文件的路径和文件名。
* `rms` 参数指定了用于 RMSD 拟合的原子组和参考结构。

通过使用 ptraj 命令提取 RMSD 拟合的 PDB 文件，您可以为后续的簇分析做好准备，并更加有效地分析和理解您的模拟结果。