g_mmpbsa工具使用

g_mmpbsa的使用教程

废话不多说,因为网上有许多教程写过g_mmpbsa的使用教程,具体见最下方参考资料[1-2]。此处我所要讲的是g_mmpbsa网站自带的一个做氨基酸突变的python脚本。因为官网没有给出具体的使用实例,所以导致很多人可能不知道或者也没有使用过,我也侥幸使用了一番,还挺有意思,因此记下使用过程。

脚本功能介绍

根据脚本中所提及的,该脚本的功能主要有两个:其一,丙氨酸扫描;其二,氨基酸突变

脚本使用方式

python环境配置

py脚本当然是使用python。此脚本使用python2写的,因此你必须安装python2,这里我建议用anaconda环境,安装和环境切换十分的方便。初次使用python mutate_traj_tpr.py -h查看帮助信息会出错,提示from modeller import *出错,不难看出这个脚本需要调用modeller来处理氨基酸突变。解决方法也十分简单,安装modeller,官网可以通过学术机构邮箱获得秘钥,安装方式我就不讲了,网上有很多。安装完成以后。再次使用python mutate_traj_tpr.py -h可以看到下图所示:

一眼就可以看出上半部分是设置氨基酸突变/丙氨酸扫描用的参数,下半部分用过GROMACS软件的人一眼就能看出是能量极小化的mdp参数部分。

辅助软件

因为该脚本是通过modeller将突变完成以后的蛋白直接进行能量最小化以备GROMACS软件模拟用,因此有几点要注意的地方:

  • 蛋白体系中不能含有其他配体或者辅酶等物质,因为GROMACS的pdb2gmx无法直接处理这样的体系。
  • 本文使用的是GROMACS 4.X版本,与脚本中出现的命令保持一致。如果你要使用最新的GROMACS软件,就需要将所有使用4.x版本的命令加上gmx前缀。另外,其中有一个4.x版本的trjcat命令可能需要删除,因为如果你直接换成新版本的gmx trjcat命令,会报错对于-f选项无法使用.gro文件.

具体实例

成功安装和配置好环境以后你就可以使用该脚本做自己的研究了。在分子动力学模拟或者生化领域,做氨基酸突变主要是为了考察该氨基酸的功能,并且对突变以后的蛋白结构和功能造成怎样的影响。本文以1AKI(鸡蛋白溶菌酶)蛋白为例,考察定点突变。
比如,我想考察一号残基位的赖氨酸(LYS)突变成半胱氨酸(CYS)以后对整个蛋白体系的结构稳定性和功能的影响。那么我们可以输入以下命令:

1
python mutate_traj_tpr.py -drWT C:\Users\liuyujie714\Desktop\test -drMT C:\Users\liuyujie714\Desktop\MT -rr "A:1"  -rnm CYS  -ff 6

这里,几个参数需要解释一下:
-drWT指的是你蛋白pdb文件所在的文件夹位置。
-drMT指的是你需要存储输出文件的文件夹位置,即突变以后的氨基酸文件。
-rr指定的是需要突变的氨基酸编号,其中 "A:1"表示A链的1号氨基酸,本例中指的是赖氨酸,如果含有其他链可以相应更改。
-rnm指定的是突变成哪种氨基酸,其中CYS指的是半胱氨酸。
-ff指定GROMACS力场。此处6表示amber99sb-ildn.ff力场。
其余参数保持默认,包括em.mdp中的参数。
好了,等待完成即可,屏幕输出如下图:

这样你会在-drMT指定的文件夹下得到如下文件:

其中,1aki_em.gro为突变后经过能量极小化以后的结构文件,你可以用它进行下一步的动力学模拟。

结构检查

过程虽然很顺利,但是还是应该检查一下突变的位置是否正确,通过可视化可以看出突变前后的变化,如下图:

看起来还是正确的突变。。
先记录到这里,以后有添加再补充。

参考资料

[1] g_mmpbsa 安装与使用笔记
[2] GROMACS教程:使用GROMACS计算MM-PBSA结合自由能