计算化学公社

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

[GROMACS] 采用CGenFF产生str文件,施加力场出现错误提示“str和mol2文件残基名不匹配”?

[复制链接 Copy URL]

15

帖子

0

威望

797

eV
积分
812

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
具体提示如下:
D:\Python 27>python cgenff_charmm2gmx.py RESI JZ4.mol2 JZ4.str charmm36-mar2019.ff
NOTE1: Code tested with python 2.7.3. Your version: 2.7.17 (v2.7.17:c2f86d86e6, Oct 19 2019, 21:01:17) [MSC v.1500 64 bit (AMD64)]


NOTE2: Please be sure to use the same version of CGenFF in your simulations that was used during parameter generation:
--Version of CGenFF detected in  JZ4.str : 4.0
--Version of CGenFF detected in  charmm36-mar2019.ff/forcefield.doc : 4.1


WARNING: CGenFF versions are not equivalent!




NOTE3: In order to avoid duplicated parameters, do NOT select the 'Include parameters that are already in CGenFF' option when uploading a molecule into CGenFF.
Error in atomgroup.py: read_mol2_coor_only: no. of atoms in mol2 (22) and top (0) are unequal
Usually this means the specified residue name does not match between str and mol2 files

JZ4.mol2

2.37 KB, 下载次数 Times of downloads: 38

JZ4.str

2.45 KB, 下载次数 Times of downloads: 46

15

帖子

0

威望

797

eV
积分
812

Level 4 (黑子)

2#
 楼主 Author| 发表于 Post on 2019-10-26 15:11:02 | 只看该作者 Only view this author
已解决,谢谢

2

帖子

0

威望

7

eV
积分
9

Level 1 能力者

3#
发表于 Post on 2020-5-17 23:16:06 | 只看该作者 Only view this author
请问楼主怎么解决的?

2

帖子

0

威望

7

eV
积分
9

Level 1 能力者

4#
发表于 Post on 2020-5-17 23:21:03 | 只看该作者 Only view this author
通过改jz4.str中的残基名跑成功了这一步,但我想问的是这里提示的CGenFF的版本号不一致是否对后续分析有影响?

15

帖子

0

威望

797

eV
积分
812

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2020-5-18 08:38:48 | 只看该作者 Only view this author
lcd83389 发表于 2020-5-17 23:21
通过改jz4.str中的残基名跑成功了这一步,但我想问的是这里提示的CGenFF的版本号不一致是否对后续分析有影 ...

你好,理论我也刚开始接触,最好统一吧,我记得好像 是2017的那个力场,版本匹配不匹配都成功了,不出现错误提示,是有note的话,可以正常运行的,后续分析也基本和官网一致; 后来小分子top文件改了方法,用的sob老师推荐中的ambertools。 官方教程我续跑了70ns,自己的就不放这了

rmsd_jz4.png (22.79 KB, 下载次数 Times of downloads: 51)

rmsd_jz4.png

12

帖子

0

威望

1365

eV
积分
1377

Level 4 (黑子)

6#
发表于 Post on 2020-10-26 11:50:35 | 只看该作者 Only view this author
本帖最后由 CHC 于 2020-10-26 14:48 编辑

楼主,你好,我也遇到了相同的问题,你说修改.str中残基名,但是str和mol2文件中的都是一样的,请问是怎么修改?
Error in atomgroup.py: read_mol2_coor_only: no. of atoms in mol2 (86) and top (0) are unequal
Usually this means the specified residue name does not match between str and mol2 files

FAD.mol2

8.92 KB, 下载次数 Times of downloads: 31

FAD.str

5.61 KB, 下载次数 Times of downloads: 29

15

帖子

0

威望

797

eV
积分
812

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2020-10-29 10:30:40 | 只看该作者 Only view this author
CHC 发表于 2020-10-26 11:50
楼主,你好,我也遇到了相同的问题,你说修改.str中残基名,但是str和mol2文件中的都是一样的,请问是怎么 ...

文件里面的名称对应起来

4

帖子

0

威望

145

eV
积分
149

Level 2 能力者

8#
发表于 Post on 2021-6-27 16:14:40 | 只看该作者 Only view this author
wlx542245 发表于 2020-10-29 10:30
文件里面的名称对应起来

楼主,您好,原始文件里的残基名不是对应的嚒,文件里面的名称对应起来的具体操作是什么呀

11

帖子

0

威望

83

eV
积分
94

Level 2 能力者

9#
发表于 Post on 2022-10-9 19:10:46 | 只看该作者 Only view this author
您好!请问这个问题最后解决了吗?我卡在这里好几天了,一直没有查到解决办法

2

帖子

0

威望

11

eV
积分
13

Level 1 能力者

10#
发表于 Post on 2023-8-6 13:29:38 | 只看该作者 Only view this author
请问解决了吗,我也卡在这里

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 21:11 , Processed in 0.197998 second(s), 29 queries , Gzip On.

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