24小时热门版块排行榜    

查看: 1258  |  回复: 5

xulinan

木虫 (小有名气)

[求助] 什么原因可能会引起MD时蛋白构象没有变化?

前几天用gromacs做了个两个蛋白相互作用的MD。初始的pdb是由单体蛋白MD的预实验中提取的,然后把两个蛋白的坐标文件合并到了一个gro文件。然后用这个合并的结构文件跑了20ns的MD,过程中没有报错。可提数据的时候问题就出现了,我发现在整个过程中那两个蛋白基本没有动过(mdp中也没做冻结或限定)。所有我觉得是不是初始结构文件做的有问题?具体步骤是这样的(之前发帖提到过)
1、用pdb2gmx分别生成两个已知pdb的top、gro和itp文件(为方便描述,两个蛋白下面分别用a、b表示)
2、修改gro文件。将其中b的gro的坐标复制到a的gro文件中(粘贴到a坐标之后就可以),同时修改原子总数。
3、修改top文件。在a的top文件中加入下列内容
;Include ligand topology
# include “b.itp”
最后在[molecules]结尾加:protein  1(分子的名字和数量)
4、之后的就是正常的MD了(建盒子、加溶剂等等)


其实我的目的是将一个蛋白加到另一个蛋白体系中,然后MD观察相互作用,请问上述方法是否欠妥?若有错误请高手指出,提供些其他的方法也很欢迎。不打算做对接,因为文献中是直接随机加入另一个蛋白的,加完后蛋白距离大于cutoff就行,而对接的结构是直接接触的,反应不了过程的变化吧。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snowinheart

银虫 (正式写手)

你好,我是做蛋白质模拟的新手。请问你是用的MS做模拟吗?
2楼2011-06-20 22:26:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulinan

木虫 (小有名气)

引用回帖:
Originally posted by snowinheart at 2011-06-20 22:26:09:
你好,我是做蛋白质模拟的新手。请问你是用的MS做模拟吗?

你好,我用的是gromacs
3楼2011-06-21 08:39:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqingwen

木虫 (正式写手)

独行者

【答案】应助回帖

★ ★
xulinan(金币+2): 不是这个问题,谢谢回复 2011-06-27 09:26:15
ghcacj(金币+2): 谢谢 2011-07-04 12:56:39
你是不是在加完溶剂水以后,没有扣除蛋白质那部分的水啊?
gromacs我听说加水的时候是把蛋白质放到水里,然后在扣除蛋白质所占位置的水的。不好意思,我是做NAMD的,以前好像也出现过这个问题:蛋白质感觉没动....
没有雨伞的孩子只能奔跑!
4楼2011-06-25 09:14:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bridgeqbf

新虫 (初入文坛)


御剑江湖(金币+1): 谢谢 2011-06-29 10:37:20
what kind of properties did you measure to say that the proteins don't move at all?
5楼2011-06-28 17:13:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hhwinhuanghe

木虫 (小有名气)

【答案】应助回帖

★ ★
ghcacj(金币+2): 谢谢 2011-07-04 12:57:01
xulinan(金币+3): 谢谢回复 2011-07-05 10:40:05
引用回帖:
Originally posted by xulinan at 2011-06-20 14:47:18:
前几天用gromacs做了个两个蛋白相互作用的MD。初始的pdb是由单体蛋白MD的预实验中提取的,然后把两个蛋白的坐标文件合并到了一个gro文件。然后用这个合并的结构文件跑了20ns的MD,过程中没有报错。可提数据的时候 ...

建议现将两个蛋白做成一个pdb文件,然后再pdb2gmx等等。。。
6楼2011-07-04 04:38:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xulinan 的主题更新
信息提示
请填处理意见