适用于GROMACS的AMBER19SB力场分享

Dempey · 发表于 Post on 2025-5-16 20:59:48

本帖最后由 Dempey 于 2026-1-21 10:44 编辑

适用于GROMACS的AMBER19SB力场分享

注意：GROMACS 2026.0版本已经集成了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: 347)

喵星大佬 · 发表于 Post on 2025-10-7 06:36:42

ASURA 发表于 2025-10-5 21:09
您好请问您解决这个问题了吗可以单纯用适用于Gromacs2025之前版本的Amber19SB力场及其配套相关力场[/ba ...

重新生成一遍拓扑文件，里面的原子类型名做了调整

ASURA · 发表于 Post on 2025-10-5 21:09:37

本帖最后由 ASURA 于 2025-10-5 21:20 编辑

杨淋锋发表于 2025-7-26 17:00
请问gmx grompp 报错“ERROR 1 [file cmap.itp, line 60]: Unknown atomtype C-*”？

您好请问您解决这个问题了吗可以单纯用适用于Gromacs2025之前版本的Amber19SB力场及其配套相关力场这个帖子中的cmap.itp替换现在力场的文件吗还是要从头生成蛋白拓扑

youknowdcf · 发表于 Post on 2025-7-30 08:58:43

感谢

杨淋锋 · 发表于 Post on 2025-7-28 09:27:41

Dempey 发表于 2025-7-26 22:55
2023.5不行哦，得2025.1或2025.2,或者可以用适用于Gromacs2025之前版本的Amber19SB力场及其配套相关力场 ...

好的，感谢

Dempey · 发表于 Post on 2025-7-26 22:55:28

杨淋锋发表于 2025-7-26 21:47
是的，2023.5

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

杨淋锋 · 发表于 Post on 2025-7-26 21:47:04

本帖最后由杨淋锋于 2025-7-28 09:27 编辑

Dempey 发表于 2025-7-26 19:20
你用的GROMACS版本是2025.1及更高版本吗

不是，2023.5

Dempey · 发表于 Post on 2025-7-26 19:20:52

杨淋锋发表于 2025-7-26 17:00
请问gmx grompp 报错“ERROR 1 [file cmap.itp, line 60]: Unknown atomtype C-*”？

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

杨淋锋 · 发表于 Post on 2025-7-26 17:00:48

请问gmx grompp 报错“ERROR 1 [file cmap.itp, line 60]: Unknown atomtype C-*”？

Dempey · 发表于 Post on 2025-6-29 17:57:51

SPHK1 发表于 2025-6-29 17:39
你好，我还想再请教一下。我在19sb的ffnonbonded.itp和ions.itp中发现有镁离子的信息，为啥蛋白离子复合物 ...

需要你在任何一个rtp文件里加入

[ MG ]
[ atoms ]
MG MG 2.00000 1

复制代码

SPHK1 · 发表于 Post on 2025-6-29 17:39:43

你好，我还想再请教一下。我在19sb的ffnonbonded.itp和ions.itp中发现有镁离子的信息，为啥蛋白离子复合物就识别不出呢？
图片分别是ffnonbonded.itp和ions.itp文件中镁离子信息的截图

SPHK1 · 发表于 Post on 2025-6-29 17:25:31

Dempey 发表于 2025-6-29 16:25
是的，我没把离子和水写进rtp里，你可以自己补充或在pdb2gmx得时候去掉Mg离子，生成拓扑后手动加入。

...

原来如此，怪我没看仔细

CAMP我没接触过，不敢瞎说；modifide是软件名的后缀，这个modified版本是sob老师编译的windows版本

Dempey · 发表于 Post on 2025-6-29 16:25:00

SPHK1 发表于 2025-6-29 15:27
我在对含有Mg2+的蛋白质生成拓扑时，发现这个力场下Gmx不能识别出镁离子

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

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

SPHK1 · 发表于 Post on 2025-6-29 15:27:36

我在对含有Mg2+的蛋白质生成拓扑时，发现这个力场下Gmx不能识别出镁离子

fhh2626 · 发表于 Post on 2025-6-10 11:32:33

看具体的问题，比如GAFF、lipid21都是在TIP3P的环境下拟合的，和这些联用的情况下OPC就不一定比TIP3P好，使用OPC3就更不合适了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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

评分 Rate

浏览过的版块