计算化学公社

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

[GROMACS] 求助:关于如何利用pdb2gmx产生聚合物的top文件的问题(.pdb文件相关)

[复制链接 Copy URL]

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 yihanxu 于 2019-2-26 02:19 编辑

老师好,大家好,我试了一下这个PE的例子: https://mailman-1.sys.kth.se/pip ... 9-March/040125.html ,但是链上有12个C的时候就不行了,请问是我的.pdb格式错了(比如哪列没有对齐)吗?而且用gview产生的.pdb在这种情况下感觉不太好用,比如会有多余的列,要手动处理一下,若是文件长的时候就很麻烦,请问怎么样能直接生成好用的.pdb格式文件呢?
谢谢指教。

PE.pdb

695 Bytes, 下载次数 Times of downloads: 61

ffoplsaa.hdb

131 Bytes, 下载次数 Times of downloads: 17

ffoplsaa.rtp

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

live each day as if it was your last

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

2#
发表于 Post on 2019-2-26 16:27:18 | 只看该作者 Only view this author
pdb格式对列有严格的要求,看http://www.wwpdb.org/documentati ... /format33/v3.3.html。某些列可以没有(比如元素),但不能错位。你的pdb明显存在错位。利用ultraedit等程序的列模式进行修改很容易。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yihanxu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2019-2-27 02:28:31 | 只看该作者 Only view this author
sobereva 发表于 2019-2-26 02:27
pdb格式对列有严格的要求,看http://www.wwpdb.org/documentation/file-format-content/format33/v3.3.html ...

谢谢老师。用gview画碳链,碳原子被标记为HETATM,是为什么呢?要手动改成ATOM,感觉麻烦呢,能直接生成高分子pdb的时候就是ATOM吗?
live each day as if it was your last

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

4#
发表于 Post on 2019-2-27 02:44:47 | 只看该作者 Only view this author
yihanxu 发表于 2019-2-27 02:28
谢谢老师。用gview画碳链,碳原子被标记为HETATM,是为什么呢?要手动改成ATOM,感觉麻烦呢,能直接生成 ...

没什么麻烦的,用ultraedit列模式五秒钟就搞定了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2019-2-27 02:58:38 | 只看该作者 Only view this author
sobereva 发表于 2019-2-26 12:44
没什么麻烦的,用ultraedit列模式五秒钟就搞定了

谢谢老师
live each day as if it was your last

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2019-3-1 09:37:29 | 只看该作者 Only view this author
本帖最后由 yihanxu 于 2019-2-28 20:48 编辑

老师好,用gmx pdb2gmx命令后,会让选择力场,用的GAFF并没被列出来,请问这时候应该怎么办呢?选择 6: AMBER99SB-ILDN protein, nucleic AMBER94是好方法吗?谢谢。
live each day as if it was your last

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

7#
发表于 Post on 2019-3-1 11:55:28 | 只看该作者 Only view this author
yihanxu 发表于 2019-3-1 09:37
老师好,用gmx pdb2gmx命令后,会让选择力场,用的GAFF并没被列出来,请问这时候应该怎么办呢?选择 6: AMB ...

本来就不会被列出来,只有把xxx.ff放到了top目录下才会在选择力场时出现xxx。得把机制搞清楚
当前用GAFF的参数,和这里选什么毫无关系
如果你基于GAFF自己写了rtp,应当新建个力场目录,比如叫GAFF.ff
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2019-3-1 13:32:09 | 只看该作者 Only view this author
本帖最后由 yihanxu 于 2019-2-28 23:57 编辑
sobereva 发表于 2019-2-28 21:55
本来就不会被列出来,只有把xxx.ff放到了top目录下才会在选择力场时出现xxx。得把机制搞清楚
当前用GAFF ...

谢谢老师。我把gaff.ff文件夹放在top目录里了,gaff.ff文件夹里有两个文件:gaff.hdb和gaff.rtp。但是使用pdb2gmx命令时备选力场里没有gaff,是怎么回事呢?

gaff.hdb

634 Bytes, 下载次数 Times of downloads: 14

gaff.rtp

5.06 KB, 下载次数 Times of downloads: 34

live each day as if it was your last

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2019-3-1 14:42:14 | 只看该作者 Only view this author
老师好,GROMACS官网给的聚乙烯的例子中.rtp文件中没有[bondedtypes]字段,但是看GROMACS手册说这是强制的,这是怎么回事呢?哪里可以查到下面各项的定义呢?谢谢。
[ bondedtypes ]
; bonds  angles  dihedrals  impropers all_dihedrals nrexcl HH14 RemoveDih
     1       1          9          4        1         3      1     0
live each day as if it was your last

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

10#
发表于 Post on 2019-3-1 18:17:52 | 只看该作者 Only view this author
yihanxu 发表于 2019-3-1 13:32
谢谢老师。我把gaff.ff文件夹放在top目录里了,gaff.ff文件夹里有两个文件:gaff.hdb和gaff.rtp。但是使 ...

照着其它.ff里的内容去写。明显你还没设置forcefield.doc

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yihanxu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

11#
发表于 Post on 2019-3-1 18:18:41 | 只看该作者 Only view this author
yihanxu 发表于 2019-3-1 14:42
老师好,GROMACS官网给的聚乙烯的例子中.rtp文件中没有字段,但是看GROMACS手册说这是强制的,这是怎么回事 ...

如果rtp里直接定义了参数,则bondedtypes虽然有,但内容可以留空,反正也不需要了。直接定义了参数的情况程序就不会再自动从bondedtypes里试图找匹配的参数

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yihanxu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

12#
 楼主 Author| 发表于 Post on 2019-3-2 15:40:50 | 只看该作者 Only view this author
老师好,请问蛋白质和高分子为什么不能用antechamber生成拓扑文件,而要用pdb2gmx命令?
谢谢。
live each day as if it was your last

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

13#
发表于 Post on 2019-3-2 15:49:35 | 只看该作者 Only view this author
yihanxu 发表于 2019-3-2 15:40
老师好,请问蛋白质和高分子为什么不能用antechamber生成拓扑文件,而要用pdb2gmx命令?
谢谢。

本来antechamber就是搞小分子的
巨大的分子本身超过了antechamber程序设计角度能处理的极限,而且其调用SQM进行半经验优化、计算AM1电荷的步骤,对于几千原子体系也根本算不动。
更何况antechamber产生的拓扑文件是基于GAFF力场的,对于蛋白质有明显更适合的蛋白质专用力场。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
yihanxu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

384

帖子

0

威望

3344

eV
积分
3728

Level 5 (御坂)

14#
 楼主 Author| 发表于 Post on 2019-3-4 12:58:32 | 只看该作者 Only view this author
本帖最后由 yihanxu 于 2019-3-5 11:49 编辑

老师好,生成的高分子的.top文件中没有bond、angle、dihedral等成键参数。为了获得GAFF下的这些参数,我是用antechamber+acpype.py生成了低聚物的.itp文件,然后手动修改低聚物的.itp文件得到的高聚物的成键参数。但是这样耗时有点多,尤其是把低聚物.itp文件中原子序号替换成原子类型时很累眼睛,请问有更好、更快捷的办法吗?谢谢。
live each day as if it was your last

516

帖子

1

威望

4765

eV
积分
5301

Level 6 (一方通行)

15#
发表于 Post on 2019-9-9 22:38:01 | 只看该作者 Only view this author
请问楼主后来找到更好的方法了吗?我最近也在编rtp,但我的重复单元里分子数有点多,好像从acpype生成的itp改rtp很不容易……

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 08:30 , Processed in 0.301843 second(s), 28 queries , Gzip On.

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