24小时热门版块排行榜    

查看: 3237  |  回复: 6

qingtian510

银虫 (初入文坛)

[求助] gromacs中金属离子的处理

在我的蛋白里有一个锌离子与3个半胱氨酸和1个组氨酸形成配位键,由于生成金属参数挺难的,所以我想通过限制锌离子和氨基酸残基中硫和氮原子之间的距离把锌离子限制在一个固定的范围内,现在我知道锌离子与硫氮原子之间的距离在2-2.6埃之间,在gromacs或amber中,这样限制怎么实现呢,请大侠指教?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

奋斗吧,少年
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangchclove

金虫 (小有名气)

顶起!
Good
2楼2013-09-24 09:19:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
月只蓝: 金币+2, 感谢指导! 2013-10-30 08:39:45
gromacs 的Position restraints ? 你先看看手册的这一节内容;

“所以我想通过限制锌离子和氨基酸残基中硫和氮原子之间的距离把锌离子限制在一个固定的范围内” ,那么可以试试把氮原子作为“reference positions Ri”,然后你的 锌离子 就在氮原子的附近了。 建模的时候,你可以把锌离子的初始位置就放在 硫氮原子之间的距离在2-2.6埃之间 的地方,这样的初始状态更容易达到平衡。
3楼2013-09-24 10:35:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingtian510

银虫 (初入文坛)

引用回帖:
3楼: Originally posted by jiaoyixiong at 2013-09-24 10:35:02
gromacs 的Position restraints ? 你先看看手册的这一节内容;

“所以我想通过限制锌离子和氨基酸残基中硫和氮原子之间的距离把锌离子限制在一个固定的范围内” ,那么可以试试把氮原子作为“reference positio ...

那个pdb文件中含有锌离子,生成.gro文件的时候,它就是在那个位置吧,完了之后我想改top文件,加个限制,我看了那个了,还是没整明白,怎么写书写限制,还有,我曾看到说在锌离子和硫氮原子直接加harnonic pot,就是设置[bond],选择func.type = 6,我设置完之后,后边还有一个gb_x,(x代表数字),对于配位键这个数字选哪个,你知道吗?我随便选了2个试了下,能量最小化的时候需要很久,完了升温的时候就运行部下去了
奋斗吧,少年
4楼2013-09-24 19:52:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingtian510

银虫 (初入文坛)

顶顶顶啊!
奋斗吧,少年
5楼2013-09-27 09:11:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fhh2626

木虫 (正式写手)

【答案】应助回帖

★
月只蓝: 金币+1, 感谢指导! 2013-10-30 08:40:17
既然是配位键 默认的参数就可以啊 在MD里面配位键应该是视为范德华力和静电力的加合的
6楼2013-09-27 11:06:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

泊小夏

新虫 (初入文坛)

顶一个
7楼2013-10-29 09:05:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 qingtian510 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见