gromacs学习

入门初学者

Posted by xiaodechao on March 29, 2021

1 环境配置

(1)构建分子模型,导出pdb,mol(生成top的工具,输入文件为pdb或mol文件)文件,使用Gaussian或者MS以及其他的一些软件。gaussian 409w MS oufc

(2)搭建盒子软件,packmol,我用的是archlinux,使用yay命令安装。win系统搜官网,下载安装。

(3)运行软件,gromacs.

(4)分析软件vmd,win10官网,linux使用命令安装即可。

(5)分析软件grace和qtgrace,打开xvg图,输出结果的分析,推荐qtgrace,提取码kkkw。

2 各类文件了解

(1) 在不修改力场的情况下,需要使用,pdb(高斯等软件构建),gro(pdb文件通过命令转化),top文件(蛋白质类型,直接命令生成,参考教程。分子通过工具生成top生成工具 ),itp文件(属于top的子文件,建议将top中的分子的信息,提取到itp文件中,include到top文件)

(2)需要添加新的力场参数,使用的文件。

  • gromacs安装目录中的top文件夹中的残疾类型(.dat文件)文件,自己定义残基名,以及类型 。
  • 打开使用的力场的目录,查阅相关文章,查找参数。
  • 打开原子类型文件,添加自定义的原子,命名前部分为元素符号,添加自定后缀,必须与前面的不一致(注意: 相关文章使用了存在原子类型,不需要定义)。
  • 打开非键相互作用文件,添加自定义原子,电荷,sigma,epsla值,是范德华力的信息.
  • 打开键相互作用文件,添加键长,键的力大小。找到键角,添加角度,和力大小。二面角(一般添加都是特殊粒子,太小,没有)修改以及其它信息一般不需要修改。
  • 定义rtp文件,描述详细的残基信息,不要在原有的rtp文件修改,新建rtp文件即可,gromacs会自动检索rtp文件。
  • 信息添加完成后,按照蛋白质教程,gmx pdb2gmx命令即可生成,top文件

3 gromacs学习网站

思家公社

中文手册

top生成网站

视频课 6t3z

4 力场选择

根据文献来

5 模拟流程

固定流程

  1. 生成gro文件,和top(以及相关的itp子文件)文件

  2. 能量最小化

  3. 温度耦合

  4. 压力耦合(由压力耦合一次,有两次的,我采用两次,不同的方式耦合)

  5. md模拟

  6. 生成文件,进行分析