24小时热门版块排行榜    

查看: 1372  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 一志愿211,0703化学310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-21 22:21 by peike
[考研] 考研化学学硕调剂,一志愿985 +5 张vvvv 2026-03-15 7/350 2026-03-21 19:23 by ColorlessPI
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 317求调剂 +9 申子申申 2026-03-19 15/750 2026-03-21 17:31 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
信息提示
请填处理意见