24小时热门版块排行榜    

查看: 1749  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:12 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:00 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 QUjhNVAcOSff 2026-09-13 7/350 2026-09-14 21:21 by vZfe6xYu34yj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:49 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
请填处理意见