本人用Pmemd做分子动力学模拟,可能是起初in文件并未设置固定蛋白中心参数的原因,在随后跑的过程中蛋白分子的部分区域跑出水盒。由于周期边界条件设置跑出水盒的部分又从另一侧出现,因此在对C@进行RMSD分析的图谱出簇升簇降现象,翻看Amber Tool说明发现RMSD分析有将蛋白固定中心的center语句,所以重新对轨迹进行RMSD分析时在导入轨迹后特意加上了:
center @CA mass origin
image origin center
两句但是RMSD分析数值和未加以上两句并未发生变化,图谱依旧出簇升簇降。请问是我的center语句相关参数设置不太恰当仍可通过分析时加上固定中心的语句弥补轨迹的缺陷,还是这条轨迹就直接不能用了需要重新跑?如果重新跑设置什么语句可以在MD过程中不跑出边界?
谢谢!