24小时热门版块排行榜    

查看: 1685  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 人气不行了 +9 fansofjerry 2026-08-21 9/450 2026-08-24 17:48 by zhaozhengxin
[基金申请] 没有任何消息-是不是就凉了 +6 图啦图啦 2026-08-24 6/300 2026-08-24 17:44 by 小小红帆
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见