24小时热门版块排行榜    

查看: 1751  |  回复: 3
当前主题已经存档。

panda_wendao

金虫 (小有名气)

[交流] 在gromacs中使用Amber力场

Gromacs是有名的分子动力学模拟软件,重写商业的Gromos,免费,而且号称是世界上最快的动力学代码。
Amber是生物分子模拟里的最好的力场之一,甚至可以不说之一。
但是由于种种原因,Gromacs中没有直接包含Amber力场。
参考 http://folding.stanford.edu/ffamber/ 的介绍,下面给出在gromacs中使用Amber力场的方法:

Install FF:
1.下载力场 :
http://folding.stanford.edu/ffamber/ffamber_v3.3.1.tar.gz
(490 kB)
http://folding.stanford.edu/ffamber/ffamber_v3.3.1-doc.tar.gz
(6.54 MB)带文献例子
2.解压缩
3.把aminoacides.dat、vdwradii.dat拷到top文件夹下覆盖,实际上一股脑考过去就完了。
4.子文件夹下的东西都拷到外面,也就是top里面。
5.修改FF.dat,修改第一行的数字,加上新添的力场数。
6.为每个新的力场增加新行,如:
     ffamber94 AMBER94 Cornell protein/nucleic forcefield
     ffamber99 AMBER99 Wang protein/nucleic acid forcefield
     ffamber99p AMBER99p protein/nucleic forcefield
     ffamber03 AMBER03 Duan protein/nucleic forcefield

Edit PDB:
1.DNA残基名改成DA、DT、DC、DG,分子为A、B。
2.所有的(*)字符换成(')。
3.端基的残基名改为DX5/DX3,或NXXX/CXXX。
4.DT中:O2->O
5.DT中:C5M->C7
6.DT中:H5M->H7
7.DC中:O2->O
8.比较麻烦:两个等位的H,形如nHX,->HXn
9.蛋白质的话,LYS->LYP

残基中的原子修改方法是将标准的DNA文件和gromacs里Amber的rtp文件逐一比较得到的,希望会对这方面的研究者有所帮助。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

panda_wendao

金虫 (小有名气)

然后就可以用pdb2gmx转pdb文件为gro和top,并继续标准步骤了。
2楼2007-07-23 12:54:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈硕真

不错!

作结构生物学和蛋白质折叠很需要!
3楼2007-07-23 19:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

itismineok

捐助贵宾 (正式写手)

收藏了!有空一定要试试!
4楼2007-08-07 15:37:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 panda_wendao 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见