24小时热门版块排行榜    

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

floon

银虫 (初入文坛)

[求助] 活性中心含双铜的蛋白酶的分子动力学模拟(Amber)

求助大侠们:蛋白为双铜中心,每个铜配位3个His,现在想用amber要跑一个空蛋白的动力学,如何开发铜中心的力场参数,就是frcmod文件的准备?
回复此楼
独自漫步
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lujunyan1118

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
chaizhm: 金币+3, 谢谢~ 2012-11-26 22:33:26
引用回帖:
2楼: Originally posted by 游子8921 at 2012-11-26 16:21:32
我现在刚开始学习amber,你知道锌酶中的锌力场参数吗?谢谢

如果是锌的参数,可以参考这篇文献,有很多现成的http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2941202/

铜的话比较少,电荷可以自己用高斯算,推荐用R.E.D
键长键角参数最好是去文献里面找,实在没办法了才选自自己发展。也可以用上面那个MCPB的方法

» 本帖已获得的红花(最新10朵)

3楼2012-11-26 22:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

游子8921

铁杆木虫 (知名作家)

游子8921

我现在刚开始学习amber,你知道锌酶中的锌力场参数吗?谢谢
为了自己的目标而奋斗,一切皆有可能
2楼2012-11-26 16:21:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

floon

银虫 (初入文坛)

送鲜花一朵
引用回帖:
3楼: Originally posted by lujunyan1118 at 2012-11-26 22:13:18
如果是锌的参数,可以参考这篇文献,有很多现成的http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2941202/

铜的话比较少,电荷可以自己用高斯算,推荐用R.E.D
键长键角参数最好是去文献里面找,实在没办法了才选 ...

你好,我的蛋白是PDB库里下载的,然后我是将活性中心的两个铜和配位的6个his拿出来,在高斯下固定两个铜,进行优化和频率以及ESP电荷的计算,然后写frcmod文件的时候,不知道如何从高斯输出文件中把拉伸力常数和扭曲力常数提取出来?
文献里的都是不同体系的键长键角,可以直接拿来用吗?
MCPB你指得是amber里的Metal Center Parameter Builder ?这个我不太会用,看什么能够比较快上手?
独自漫步
4楼2012-11-27 21:52:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lujunyan1118

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
zh1987hs: 金币+3, 谢谢 2012-11-28 10:37:29
amber里面有个叫MTK++的组件,在MTK++的manual里面有一章讲的是MCPB,有个tutorial可以跟着做,不过也是比较麻烦。它的原理就是从gauss的log文件里提取频率等信息,转化成键能参数。手动提到话我也不太会,可能要用到专门的脚本。

不同的体系关系我觉得不太大,因为你算得时候本来就只把活性中心附近的残基抠出来的,没有考虑其他残基的作用,优化以后可能就可pdb结构不一样了
不过还是越接近越好,起码也要是6个his和两个cu配位的结构(也许不太好找)。。。你可以发邮件给发展MCPB方法的那个组,或在amber论坛上求助,看有没有现成的参数

» 本帖已获得的红花(最新10朵)

5楼2012-11-27 22:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 14:22 by 0aDbZAghdVq9
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 10:19 by xSrQNSSjcgoX
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 09:35 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 07:56 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:33 by xSrQNSSjcgoX
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:43 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 uKOSCocCb6Fm 2026-11-12 8/400 2026-11-15 02:26 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:51 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 FChKX8meu1ds 2026-11-09 8/400 2026-11-14 22:04 by Equinoxhua
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-14 21:58 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[基金申请] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:14 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 09:21 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
信息提示
请填处理意见