cpptraj 中 pdb 格式选项：完整指南

在使用 cpptraj 将轨迹数据输出为 pdb 文件时，您可以通过各种选项来控制输出 pdb 文件的结构。这些选项包括链 ID、残基名称、原子类型、B 因子和占用率数据。

基本选项

[dumpq | parse | vdw]：选择输出的电荷和半径数据类型。
- dumpq：输出 PQR 格式，将电荷（单位为 e-）写入占用率列，将 GB 半径写入 B 因子列。
- parse：输出 PQR 格式，将电荷写入占用率列，将 PARSE 半径写入 B 因子列。
- vdw：输出 PQR 格式，将电荷写入占用率列，将 vdW 半径写入 B 因子列。
[pdbres]：使用 PDB V3 残基名称。
[pdbatom]：使用 PDB 原子名称。
[pdbv3]：强制使用 PDB V3 格式。
[teradvance]：将终止符（TER）记录添加到每个残基的末尾。
[terbyres | pdbter | noter]：控制终止符记录的行为。
- terbyres：每当遇到新残基时，都添加终止符记录。
- pdbter：使用 PDB 文件中的终止符记录。
- noter：不添加终止符记录。
[model | multi]：控制输出 pdb 文件中模型的数量。
- model：仅输出当前模型。
- multi：输出所有模型。
[chainid ]：指定输出 pdb 文件中的链 ID。
[sg ]：指定要输出的原子组。
[include_ep]：包含额外属性数据。
[conect]：输出 CONNECT 记录。
[keepext]：保留输入 pdb 文件中的扩展信息。
[usecol21]：使用第 21 列来存储原子类型。

B 因子和占用率选项

[bfacdefault <#>]：指定默认的 B 因子值。
[occdefault <#>]：指定默认的占用率值。
[bfacdata ]：指定用于 B 因子的数据源。
[occdata ]：指定用于占用率的数据源。
[bfacbyres]：根据残基类型设置 B 因子。
[occbyres]：根据残基类型设置占用率。
[bfacscale]：缩放 B 因子。
[occscale]：缩放占用率。
[bfacmax ]：设置 B 因子的最大值。
[occmax ]：设置占用率的最大值。
[adpdata ]：指定用于 ADPT 记录的数据源。

例子：将每帧输出为单独的 pdb 文件

要将轨迹中每一帧都输出为一个 pdb 文件，您可以使用以下命令：

# 从输入轨迹文件中读取帧
trajin input.nc

# 设置参考结构，如果需要的话
reference reference.pdb

# 设置拓扑文件
parm topology.prmtop

# 使用 dump 命令将每一帧输出为 pdb 文件
for i in {1..100}
do
trajout frame$i.pdb pdb nobox
go
done

在这个例子中，我们首先使用 trajin 命令将输入轨迹文件 (input.nc) 加载到 cpptraj 中。然后，我们使用 reference 命令设置参考结构（如果需要的话）。接下来，使用 parm 命令加载拓扑文件 (topology.prmtop)。最后，使用 for 循环将每一帧使用 trajout 命令输出为 pdb 文件，文件名为 frame$i.pdb，其中 $i 是循环变量，范围从 1 到 100。

请确保将上述命令中的输入文件名和输出文件名替换为您实际使用的文件名。

提示：

使用 pdbset 命令可以修改输出 pdb 文件的格式。
使用 write 命令可以将轨迹数据输出为其他格式，例如 dcd 或 netcdf。

希望这篇指南可以帮助您在使用 cpptraj 时更好地理解和使用 pdb 格式选项。