24小时热门版块排行榜    

查看: 1662  |  回复: 6

floon

银虫 (初入文坛)

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

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

游子8921

铁杆木虫 (知名作家)

游子8921

我现在刚开始学习amber,你知道锌酶中的锌力场参数吗?谢谢
为了自己的目标而奋斗,一切皆有可能
2楼2012-11-26 16:21:32
已阅   回复此楼   关注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的回帖

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的回帖

floon

银虫 (初入文坛)

送鲜花一朵
引用回帖:
5楼: Originally posted by lujunyan1118 at 2012-11-27 22:46:51
amber里面有个叫MTK++的组件,在MTK++的manual里面有一章讲的是MCPB,有个tutorial可以跟着做,不过也是比较麻烦。它的原理就是从gauss的log文件里提取频率等信息,转化成键能参数。手动提到话我也不太会,可能要用 ...

谢谢,最后问下MCPB那个组的邮箱写哪一个,是Martin B. Peters (martin.b.peters@gmail.com)这个吗?
是分子模拟论坛上的amber板块吗?还是Center for Structural Biology的AMBER Archive ?  问的有点多,不要介意啊.
独自漫步
6楼2012-11-28 10:27:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lujunyan1118

金虫 (小有名气)

【答案】应助回帖

没事~
http://archive.ambermd.org/ 这个上面找或者问

邮箱应该就是那个,他们组有中国学生
7楼2012-11-28 12:41:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 floon 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 9/450 2026-08-02 16:32 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:32 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:31 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 13/650 2026-08-02 16:27 by jHcrHMhbKCFV
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:24 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:23 by jHcrHMhbKCFV
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见