使用 CHARMMff 文件计算蛋白质-配体结合自由能（单轨迹法）

信息

该示例可以在存储库文件夹中的examples/Protein_ligand_CHARMMff目录中找到。如果您之前没有使用 gmx_MMPBSA_test，请使用downgit从 gmx_MMPBSA GitHub 存储库下载特定文件夹。

要求¶

在这种情况下，gmx_MMPBSA需要：

需要输入文件	类型	描述
输入参数文件	`in`	输入文件包含有关将要执行的计算类型的所有规范
MD 结构 + 质量 (db) 文件	`tpr` `pdb`	包含系统坐标的结构文件
受体和配体组	`integers`	索引文件中的受体和配体组编号
轨迹文件	`xtc` `pdb` `trr`	最终的 GROMACS MD 轨迹，已拟合且没有 pbc。
拓扑文件	`top`	考虑到属于拓扑文件的 *.itp 文件也应该存在于该文件夹中
参考结构文件	`pdb`	具有所需链 ID 和残基编号分配的复杂参考结构文件（不含氢）

-> 必须定义 ---> 可选，但推荐 --> 可选

在此处查看gmx_MMPBSA 命令行中所有标志的详细列表

命令行¶

也就是说，一旦进入包含所有文件的文件夹，命令行将如下所示：

连续剧与MPIgmx_MMPBSA_test

gmx_MMPBSA -O -i mmpbsa.in -cs com.tpr -ci index.ndx -cg 1 13 -ct com_traj.xtc -cp topol.top -o FINAL_RESULTS_MMPBSA.dat -eo FINAL_RESULTS_MMPBSA.csv

其中mmpbsa.in输入文件是包含以下行的文本文件：

PB 计算的示例输入文件
Sample input file for PB calculationThis input file is meant to show only that gmx_MMPBSA works. Althought,we tried to used the input files as recommended in the Amber manual,some parameters have been changed to perform more expensive calculationsin a reasonable amount of time. Feel free to change the parameters according to what is better for your system.&generalsys_name="Prot-Lig-CHARMM",startframe=1,endframe=4,# In gmx_MMPBSA v1.5.0 we have added a new PB radii set named charmm_radii. # This radii set should be used only with systems prepared with CHARMM force fields. # Uncomment the line below to use charmm_radii set#PBRadii=7,/&pb# radiopt=0 is recommended which means using radii from the prmtop file for both the PB calculation and for the NP# calculationistrng=0.15, fillratio=4.0, radiopt=0/

记住

radiopt = 0建议使用prmtop文件中的半径

请在此处查看gmx_MMPBSA输入文件中所有选项的详细列表以及几个示例

注意事项¶

在这种情况下，遵循单一轨迹 (ST) 近似，这意味着受体和配体结构和轨迹将从复合物的结构和轨迹中获得。为此，需要一个 MD Structure+mass(db) 文件 ( com.tpr)、一个索引文件 ( index.ndx)、一个轨迹文件 ( com_traj.xtc) 以及索引文件 ( ) 中的受体和配体组编号1 13。mmpbsa.in输入文件将包含 MM/PB(GB)SA 计算所需的所有参数。在这种情况下，还需要一个拓扑文件（强制性的）来生成琥珀色格式的拓扑文件，其中包含 CHARMM 力场的所有项。

将保存包含所有统计信息的纯文本输出文件（默认值FINAL_RESULTS_MMPBSA.dat：）和包含每个计算中每个帧的所有能量项的 CSV 格式输出文件。'-eo' 标志中的文件名将强制以 [.csv] 结尾（FINAL_RESULTS_MMPBSA.csv在本例中）。仅当在命令行上指定时才写入此文件。

笔记

计算完成后，可以分析结果gmx_MMPBSA_ana（如果-nogui未在命令行中使用标志）。请查看gmx_MMPBSA_ana部分以获取更多信息

源文件链接：https://valdes-tresanco-ms.github.io/gmx_MMPBSA/dev/examples/Protein_ligand_CHARMMff/

上一篇gromacs文件中关于top的解读下一篇pdb格式解读

最新评论

请先登录才能进行回复登录

使用 CHARMMff 文件计算蛋白质-配体结合自由能（单轨迹法）

要求¶

命令行¶

注意事项¶

电话直呼