24小时热门版块排行榜    

查看: 1943  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

3115321

铁虫 (正式写手)

[求助] gromacs pull的问题!10金币 已有1人参与

设定了一个离子通道蛋白A,和一个离子B,想让B 通过pull的umbrella模式通过通道A
我的办法1,设定蛋白A是参考组,pull B,但是离子B停留在蛋白A中心不再向下穿透。
办法2,在通道A另一端设置一个参考离子C,把B向离子C 拉动,但是离子没有向蛋白方向运动而是从另一侧反向接近C。
请问怎么回事,我mdp设置砍了好几遍也看不出解决问题的方法
有没有人知道啊?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
月只蓝: 金币+2, 感谢指导! 2013-12-12 08:30:44
3115321: 金币+10, ★★★★★最佳答案 2013-12-22 22:31:35
引用回帖:
5楼: Originally posted by 3115321 at 2013-12-11 21:39:46
用direction也跨不过中间,只能分两步走,先到中心位置,然后再另外设置一个参考点跨越后半部分。是不是应该这样操作?...

我觉得最正确的操作是:
仔细阅读手册上有关pull  的有关命令,再结合你要模拟的内容,搞明白如何设置参数。

用gromacs 的 pull code 我也做过一些模拟, 根据“设定了一个离子通道蛋白A,和一个离子B,想让B 通过pull的umbrella模式通过通道A” 这个目标,要是我的话,“pull_geometry   = direction” 肯定行的。只是你在3楼贴的”pull-vec1 “ 命令也是要做调整的,我给你提示了,可能你还没有理解手册上的内容。还有一个”pull-rate1 = 0.001“,我不知道你这个值是根据什么设置的?

现在的提示已经非常明确了,多多思考吧。
祝好运!
6楼2013-12-11 22:09:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★
感谢参与,应助指数 +1
月只蓝: 金币+1, 鼓励交流! 2013-12-11 18:21:31
你不把你的 MDP文件的pull部分贴出来,恐怕只有神仙知道,
还有你仔细确认过你设置的 A B吗?
2楼2013-12-11 17:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3115321

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2013-12-11 17:00:03
你不把你的 MDP文件的pull部分贴出来,恐怕只有神仙知道,
还有你仔细确认过你设置的 A B吗?

pull = umbrella
pull-geometry = distance
pull-dim =N N Y
pull-ngroups = 1
pull-start = yes
pull-pbcatom1 = 0
pull-group0 = protein
pull-group1 = CL
pull-vec1 = 0 0 0
pull-init1 = 0
pull-rate1 = 0.001
pull-k1 = 1000
freezegrps = na
freezedim = Y Y Y
3楼2013-12-11 17:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★
月只蓝: 金币+2, 感谢指导! 2013-12-12 08:30:39
1、你真的理解了你贴的那些命令的意思吗???
     建议你把手册这段命令再仔细阅读理解一下,读完也许你就知道为你这些命令运行完,会是你在一楼说的那样。

2、“设定了一个离子通道蛋白A,和一个离子B,想让B 通过pull的umbrella模式通过通道A”
      如果是我,我会选择“pull_geometry   = direction” ,提示到此,其他的设置请仔细阅读手册吧。

3、祝好运!
4楼2013-12-11 19:53:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:14 by TgvkbZZl9tqB
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:09 by TgvkbZZl9tqB
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-11 16:41 by TgvkbZZl9tqB
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 3/150 2026-11-11 15:43 by TgvkbZZl9tqB
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-11 15:24 by TgvkbZZl9tqB
[基金申请] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-11 15:13 by TgvkbZZl9tqB
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 YHiw1gou1Eda 2026-11-10 4/200 2026-11-11 14:13 by cdqGD6ftHmbl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 14:02 by cdqGD6ftHmbl
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 11:41 by gRND4NVAZ8VG
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 11:38 by gRND4NVAZ8VG
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 10:46 by 1Nb0IVlmV6v0
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 10:10 by 1Nb0IVlmV6v0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 07:06 by Z7qfCsdIiD1T
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
信息提示
请填处理意见