计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2187|回复 Reply: 19
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 适用于GROMACS的AMBER19SB力场分享

[复制链接 Copy URL]

86

帖子

3

威望

1751

eV
积分
1897

Level 5 (御坂)

本帖最后由 Dempey 于 2025-5-16 21:40 编辑

适用于GROMACS的AMBER19SB力场分享



GROMACS 2025版本一个特性就是支持了按残基匹配CMAP项(见GROMACS Release Note),这就可以支持AMBER19SB了,故本人尝试制作了一个AMBER19SB,其中DNA力场使用AMBER25推荐的OL24,RNA力场使用OL3,内置了三种水模型——OPC、OPC3、TIP4P-ice,并加入了Merz参数化的一系列离子的参数(OPC、OPC3)。现在分享给大家,有问题反馈给我,持续更新。



一些碎碎念:

1. 使用AMBER19SB力场就需配套使用OPC水,想减小水的开销就用OPC3水,尽量不要使用TIP3P。

As a result, we strongly recommend using ff19SB with OPC, and we recommend against use with TIP3P.

2. AMBER19SB力场的键参数同AMBER14SB一样,缺失了组氨酸的两个参数,已经补全在本力场中了。
3. 力场包中ions.itp中是适用于OPC水的离子参数,ions_opc3.itp是适用于OPC3水的参数,可以非常方便的切换;但这两个文件中并没有写全全部离子,可以查看ffnonbonded.itp中的离子种类自行添加。
4. 后续有时间会尝试加入GLYCAM_06j和lipid21的参数,但这两个工作量有点大
5. 如有大家在使用时发现问题请一定要反馈给我。

amber19sb.tar.gz (205.46 KB, 下载次数 Times of downloads: 138)

评分 Rate

参与人数
Participants 15
威望 +1 eV +56 收起 理由
Reason
大村驴 + 5 谢谢
mtcldys + 4 谢谢
RandomError + 5
limalvis + 3 谢谢分享
chenbq18 + 5 好物!
csm + 3 谢谢分享
爆旋陀螺 + 3 赞!
jackshoulder + 1 好物!
jamesl + 5 好物!
student0618 + 5 好物!
丁越 + 5 赞!
AxiEJohn + 5
sobereva + 1
lmch + 5
牧生 + 2 好物!

查看全部评分 View all ratings

44

帖子

0

威望

890

eV
积分
934

Level 4 (黑子)

20#
发表于 Post on 2025-7-30 08:58:43 | 只看该作者 Only view this author
感谢

7

帖子

0

威望

167

eV
积分
174

Level 3 能力者

19#
发表于 Post on 2025-7-28 09:27:41 | 只看该作者 Only view this author
Dempey 发表于 2025-7-26 22:55
2023.5不行哦,得2025.1或2025.2,或者可以用适用于Gromacs2025之前版本的Amber19SB力场及其配套相关力场 ...

好的,感谢

86

帖子

3

威望

1751

eV
积分
1897

Level 5 (御坂)

18#
 楼主 Author| 发表于 Post on 2025-7-26 22:55:28 | 只看该作者 Only view this author

2023.5不行哦,得2025.1或2025.2,或者可以用适用于Gromacs2025之前版本的Amber19SB力场及其配套相关力场中提供的力场文件

7

帖子

0

威望

167

eV
积分
174

Level 3 能力者

17#
发表于 Post on 2025-7-26 21:47:04 | 只看该作者 Only view this author
本帖最后由 杨淋锋 于 2025-7-28 09:27 编辑
Dempey 发表于 2025-7-26 19:20
你用的GROMACS版本是2025.1及更高版本吗

不是,2023.5

86

帖子

3

威望

1751

eV
积分
1897

Level 5 (御坂)

16#
 楼主 Author| 发表于 Post on 2025-7-26 19:20:52 | 只看该作者 Only view this author
杨淋锋 发表于 2025-7-26 17:00
请问gmx grompp 报错“ERROR 1 [file cmap.itp, line 60]:  Unknown atomtype C-*”?

你用的GROMACS版本是2025.1及更高版本吗

7

帖子

0

威望

167

eV
积分
174

Level 3 能力者

15#
发表于 Post on 2025-7-26 17:00:48 | 只看该作者 Only view this author
请问gmx grompp 报错“ERROR 1 [file cmap.itp, line 60]:  Unknown atomtype C-*”?

86

帖子

3

威望

1751

eV
积分
1897

Level 5 (御坂)

14#
 楼主 Author| 发表于 Post on 2025-6-29 17:57:51 | 只看该作者 Only view this author
SPHK1 发表于 2025-6-29 17:39
你好,我还想再请教一下。我在19sb的ffnonbonded.itp和ions.itp中发现有镁离子的信息,为啥蛋白离子复合物 ...

需要你在任何一个rtp文件里加入
  1. [ MG ]
  2. [ atoms ]
  3.    MG     MG           2.00000     1
复制代码

6

帖子

0

威望

85

eV
积分
91

Level 2 能力者

13#
发表于 Post on 2025-6-29 17:39:43 | 只看该作者 Only view this author
你好,我还想再请教一下。我在19sb的ffnonbonded.itp和ions.itp中发现有镁离子的信息,为啥蛋白离子复合物就识别不出呢?
图片分别是ffnonbonded.itp和ions.itp文件中镁离子信息的截图

202506291736117932..png (31.11 KB, 下载次数 Times of downloads: 112)

202506291736117932..png

202506291737167025..png (33.22 KB, 下载次数 Times of downloads: 108)

202506291737167025..png

6

帖子

0

威望

85

eV
积分
91

Level 2 能力者

12#
发表于 Post on 2025-6-29 17:25:31 | 只看该作者 Only view this author
Dempey 发表于 2025-6-29 16:25
是的,我没把离子和水写进rtp里,你可以自己补充或在pdb2gmx得时候去掉Mg离子,生成拓扑后手动加入。

...

原来如此,怪我没看仔细
CAMP我没接触过,不敢瞎说;modifide是软件名的后缀,这个modified版本是sob老师编译的windows版本

86

帖子

3

威望

1751

eV
积分
1897

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2025-6-29 16:25:00 | 只看该作者 Only view this author
SPHK1 发表于 2025-6-29 15:27
我在对含有Mg2+的蛋白质生成拓扑时,发现这个力场下Gmx不能识别出镁离子

是的,我没把离子和水写进rtp里,你可以自己补充或在pdb2gmx得时候去掉Mg离子,生成拓扑后手动加入。

另2020.6-MODIFIED支持残基分辨的CMAP吗,我不知道modified了什么

6

帖子

0

威望

85

eV
积分
91

Level 2 能力者

10#
发表于 Post on 2025-6-29 15:27:36 | 只看该作者 Only view this author
我在对含有Mg2+的蛋白质生成拓扑时,发现这个力场下Gmx不能识别出镁离子

202506291526147542..png (30.8 KB, 下载次数 Times of downloads: 112)

202506291526147542..png

1149

帖子

6

威望

6631

eV
积分
7900

Level 6 (一方通行)

9#
发表于 Post on 2025-6-10 11:32:33 | 只看该作者 Only view this author
看具体的问题,比如GAFF、lipid21都是在TIP3P的环境下拟合的,和这些联用的情况下OPC就不一定比TIP3P好,使用OPC3就更不合适了

18

帖子

0

威望

43

eV
积分
61

Level 2 能力者

8#
发表于 Post on 2025-6-8 19:45:51 | 只看该作者 Only view this author
19SB相对14SB的优势到现在为止都没看到显著的报道

129

帖子

0

威望

1495

eV
积分
1624

Level 5 (御坂)

7#
发表于 Post on 2025-6-6 20:44:12 | 只看该作者 Only view this author
Dempey 发表于 2025-6-6 09:02
是的,现在格式转换很方便,我比较熟悉tleap+parmed,但是我发现可能有点小问题,也有可能是我操作不对: ...

看来是我想简单了,以为parmed或者acpype都能够正确的转化Amber ff19SB的格式呢。。。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 21:51 , Processed in 0.261015 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list