| 查看: 1539 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
活性中心含双铜的蛋白酶的分子动力学模拟（Amber）
|
|||
| 求助大侠们：蛋白为双铜中心，每个铜配位3个His，现在想用amber要跑一个空蛋白的动力学，如何开发铜中心的力场参数，就是frcmod文件的准备？ |
回复此楼» 收录本帖的淘帖专辑推荐
 amber |
» 猜你喜欢
366求调剂
已经有9人回复
-
材料工程085601，270求调剂
已经有37人回复
-
279学硕食品专业求调剂院校
已经有18人回复
-
290调剂生物0860
已经有31人回复
-
一志愿085802 323分求调剂
已经有13人回复
-
277求调剂
已经有23人回复
-
322求调剂，08工科
已经有4人回复
-
材料工程281还有调剂机会吗
已经有30人回复
-
化学070300 求调剂
已经有16人回复
-
化工学硕294分，求导师收留
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求，铜离子及EDTA对蛋白酶活性的影响与原因
已经有1人回复
lujunyan1118
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1306.9
- 帖子: 73
- 在线: 20.4小时
- 虫号: 871655
- 注册: 2009-10-14
- 性别: GG
- 专业: 生物大分子结构与功能
【答案】应助回帖
★ ★ ★
zh1987hs: 金币+3, 谢谢 2012-11-28 10:37:29
zh1987hs: 金币+3, 谢谢 2012-11-28 10:37:29
|
amber里面有个叫MTK++的组件，在MTK++的manual里面有一章讲的是MCPB，有个tutorial可以跟着做，不过也是比较麻烦。它的原理就是从gauss的log文件里提取频率等信息，转化成键能参数。手动提到话我也不太会，可能要用到专门的脚本。 不同的体系关系我觉得不太大，因为你算得时候本来就只把活性中心附近的残基抠出来的，没有考虑其他残基的作用，优化以后可能就可pdb结构不一样了 不过还是越接近越好，起码也要是6个his和两个cu配位的结构（也许不太好找）。。。你可以发邮件给发展MCPB方法的那个组，或在amber论坛上求助，看有没有现成的参数 |
floon5楼2012-11-27 22:46:51
游子8921
铁杆木虫 (知名作家)
游子8921
- 应助: 12 (小学生)
- 金币: 8915.2
- 散金: 132
- 红花: 5
- 帖子: 6402
- 在线: 750.9小时
- 虫号: 685694
- 注册: 2009-01-01
- 性别: GG
- 专业: 理论和计算化学
2楼2012-11-26 16:21:32
lujunyan1118
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1306.9
- 帖子: 73
- 在线: 20.4小时
- 虫号: 871655
- 注册: 2009-10-14
- 性别: GG
- 专业: 生物大分子结构与功能
【答案】应助回帖
★ ★ ★
chaizhm: 金币+3, 谢谢~ 2012-11-26 22:33:26
chaizhm: 金币+3, 谢谢~ 2012-11-26 22:33:26
|
如果是锌的参数，可以参考这篇文献，有很多现成的http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2941202/ 铜的话比较少，电荷可以自己用高斯算，推荐用R.E.D 键长键角参数最好是去文献里面找，实在没办法了才选自自己发展。也可以用上面那个MCPB的方法 |
» 本帖已获得的红花（最新10朵）
3楼2012-11-26 22:13:18
4楼2012-11-27 21:52:38
10