為愛停留 发表于 2019-1-3 10:00:08

Gromacs中,小分子力场的生成

大家好http://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/victory.gif,初学GROMACS,在分子力场方面遇到了一个棘手的问题http://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/sad.gif。我做的是小分子药物与另一分子对接之后的动力学模拟http://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/rolleyes.gif。首先我采用autodock将两分子对接,对接结果采用Pymol保存成复合物pdb格式的输入文件;然后用Gromscs进行动力学模拟。GMX中第一步是pdb2gmx,之后会提示选择力场,我选择gromos力场(其实各个力场都试了),出现报错(http://muchongimg.xmcimg.com/data/bcs/2016/0120/w689h1901538_1453305333_363.png),经询问查找,是对应力场下残基数据库没有相应的残基,但同时我的小分子也没有相应的力场。经过大家的帮助,我得知大概有两种方法可以解决这个问题,一个是用acpype来帮助生成力场,另一个是用prodrg来生成.itp文件,#include到top文件中(说的不详细)。查找资料已经接近一周了,期间得到了各位老师同仁的帮助,但还是没有解决。本人愚钝http://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/cry.gif学的不是很明白,希望大家能在gromacs中生成小分子力场的方法或者实例(实例的链接也可以)or解决残基fetal error方法方面给出帮助。Sample Text,鄙人在此感谢大家http://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/hand.gifhttp://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/hand.gifhttp://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/hand.gifhttp://muchongimg.xmcimg.com/data/emuch_bbs_images/smilies/hand.gif

http://muchongimg.xmcimg.com/data/bcs/2016/0120/w628h1901538_1453304750_673.png
输入命令.png

心软脾气硬 发表于 2019-1-3 10:00:55


http://sobereva.com/266

她不美 发表于 2019-1-3 10:01:08

你先查找一下对接进去的复合物有没有残基缺失,或者不明原子等情况,如果有的话在pdb2gmx之前就应该修正好。然后再运行pdb2gmx加力场试试看。还有小分子应该先从复合物里抠出来,再对蛋白进行力场的赋予。

為愛停留 发表于 2019-1-3 10:01:13

我不是做的蛋白质呀,我做的是两个小分子的对接。

她不美 发表于 2019-1-3 10:01:30

两个小分子也可以对接?我们要不做蛋白与小分子的对接,要不就是蛋白和蛋白的对接。还没听说过小分子之间的对接呢。不吝赐教

為愛停留 发表于 2019-1-3 10:01:33

谢谢你,我又查了查资料,pdb2gmx工具,是gromacs通用工具,只能用于蛋白质,DNA、RNA等的生物分子,其他的要用其他的方法。

她不美 发表于 2019-1-3 10:01:48

不客气,祝你好运

為愛停留 发表于 2019-1-3 10:01:50

嘿嘿,前期不太懂。不过基本上已解决……
页: [1]
查看完整版本: Gromacs中,小分子力场的生成