求助：关于如何利用pdb2gmx产生聚合物的top文件的问题（.pdb文件相关）

yihanxu · 发表于 Post on 2019-9-10 05:35:51

本帖最后由 yihanxu 于 2019-9-9 15:37 编辑

tjuptz 发表于 2019-9-9 08:38
请问楼主后来找到更好的方法了吗？我最近也在编rtp，但我的重复单元里分子数有点多，好像从acpype生成的itp ...

还是用的这种方法。如果重复单元的基团比较多，确实改起来很麻烦，我的重复单元没超过三个环，改起来还行。

tjuptz · 发表于 Post on 2019-9-10 12:07:14

本帖最后由 tjuptz 于 2019-9-10 12:54 编辑

sobereva 发表于 2019-3-1 11:55
本来就不会被列出来，只有把xxx.ff放到了top目录下才会在选择力场时出现xxx。得把机制搞清楚
当前用GAFF ...

老师，在GROMACS中，基于acpype获得高分子单元的GAFF参数，并对聚合物用pdb2gmx生成链的top文件，必须新建GAFF.ff而不能基于amber99sb-ildn.ff做吗？因为后面我也想用隐式溶剂。另外，编rtp时需要把itp里的参数都写进去吗？求指点

sobereva · 发表于 Post on 2019-9-11 09:52:55

tjuptz 发表于 2019-9-10 12:07
老师，在GROMACS中，基于acpype获得高分子单元的GAFF参数，并对聚合物用pdb2gmx生成链的top文件，必须新 ...

可以用AMBER的.ff。和GAFF的原子类型并不冲突，力场的default部分也都是兼容的。只要把acpype给出的额外参数补进去就行了

rtp里不是非得直接体现出参数。只要把成键关系在里面定义了就够了，grompp的时候自动就会从力场文件里根据原子类型去找对应的参数。

tjuptz · 发表于 Post on 2019-9-11 13:52:44

本帖最后由 tjuptz 于 2019-9-11 19:37 编辑

sobereva 发表于 2019-9-11 09:52
可以用AMBER的.ff。和GAFF的原子类型并不冲突，力场的default部分也都是兼容的。只要把acpype给出的额外 ...

我明白了，谢谢老师。

tjuptz · 发表于 Post on 2019-9-13 15:56:11

本帖最后由 tjuptz 于 2019-9-15 10:08 编辑

sobereva 发表于 2019-9-11 09:52
可以用AMBER的.ff。和GAFF的原子类型并不冲突，力场的default部分也都是兼容的。只要把acpype给出的额外 ...

再次求助社长。我把低聚体acpype给出的参数补到力场文件里，其中二面角字段的原子名改成原子类型后遇到如下问题：
os  c3  c3  hc 9    0.00 0.00000 0 ;
os  c3  c3  hc 9    0.00 1.04600 1 ;
c3  c3  c3  c3 9    0.00 0.75312 3 ;
c3  c3  c3  c3 9    180.00 0.83680 1 ;
c3  c3  c3  c3 9    180.00 1.04600 2 ;
我查看itp中对应原始原子名字段是如下这样的：
17    18    19    51    9    0.00 0.00000 0 ; O17- C18- C19- H51
17    18    19    51    9    0.00 1.04600 1 ; O17- C18- C19- H51
17    18    19    52    9    0.00 0.00000 0 ; O17- C18- C19- H52
17    18    19    52    9    0.00 1.04600 1 ; O17- C18- C19- H52

30    31    32    76    9    0.00 0.00000 0 ; O30- C31- C32- H76
30    31    32    76    9    0.00 1.04600 1 ; O30- C31- C32- H76
30    31    32    77    9    0.00 0.00000 0 ; O30- C31- C32- H77
30    31    32    77    9    0.00 1.04600 1 ; O30- C31- C32- H77

18    19    20    21    9    0.00 0.75312 3 ; C18- C19- C20- C21
18    19    20    21    9 180.00 0.83680 1 ; C18- C19- C20- C21
18    19    20    21    9 180.00 1.04600 2 ; C18- C19- C20- C21

19    20    21    22    9    0.00 0.75312 3 ; C19- C20- C21- C22
19    20    21    22    9 180.00 0.83680 1 ; C19- C20- C21- C22
19    20    21    22    9 180.00 1.04600 2 ; C19- C20- C21- C22

20    21    22    23    9    0.00 0.75312 3 ; C20- C21- C22- C23
20    21    22    23    9 180.00 0.83680 1 ; C20- C21- C22- C23
20    21    22    23    9 180.00 1.04600 2 ; C20- C21- C22- C23

…………
26    27    28    29    9    0.00 0.75312 3 ; C26- C27- C28- C29
26    27    28    29    9 180.00 0.83680 1 ; C26- C27- C28- C29
26    27    28    29    9 180.00 1.04600 2 ; C26- C27- C28- C29
我的单体里面有长碳链C18-C29，通过acpype获得碳链上的C都是c3原子类型，但这些c3碳原子二面角acpype却给出了多项参数，我想起来，这是多项叠加描述扭转势。这样补到ffbonded.itp中就出现了同样的c3-c3-c3-c3后跟着不同参数的情况；同样的，c3碳原子上的H都是hc原子类型，也遇到了这个问题。我看其他帖子中提到，力常数为零的可以忽略，但对于c3-c3-c3-c3还有三项参数。这时候我该怎么处理呢？把这些参数直接写到rtp里吗，我尝试了一下，参数太多不太现实；而且根据单体的itp参数补还得考虑多聚体成键项改变。

更新：我看了手册，说是二面角函数类型9就可以用多参数，所以我就这样补进ffbonded.itp就可以了吧？

diaok · 发表于 Post on 2019-9-14 12:29:28

tjuptz 发表于 2019-9-13 15:56
再次求助社长。我把acpype给出的二面角字段的原子名改成原子类型后遇到如下问题：
os c3 c3 hc 9 ...

写rtp时有没有涉及到片段连接时应该补全的键角二面角之类？
要是不好解决可以直接做成分开的片段，然后用parmed合并成一个拓扑分子，再根据原子序号补上连接的键角之类参数

tjuptz · 发表于 Post on 2019-9-14 12:39:25

diaok 发表于 2019-9-14 12:29
写rtp时有没有涉及到片段连接时应该补全的键角二面角之类？
要是不好解决可以直接做成分开的片段，然后 ...

涉及到了，我用共聚体重复单元的itp应该能获得这些参数吧

diaok · 发表于 Post on 2019-9-14 16:02:34

tjuptz 发表于 2019-9-14 12:39
涉及到了，我用共聚体重复单元的itp应该能获得这些参数吧

gmx给的那个rtp示例里也没看到键角二面角之类信息，应该是用ffbond里面给补全了
实现gaff应该需要把有关参数都补充完整
It is easier to create .rtp file first, since you only need partial charge and bond information. The dihedral and angle will be automatically generated. Then you run grompp to check for any non-standard bond/angle/dihedral, correct for them (modify the ffbonded.itp).
https://www.researchgate.net/pos ... mulation_in_Gromacs

我一般用的另一种方法，先做好单独的片段，然后直接整合pdb连成完整的聚合物，用antechamber来判断拓扑结构得到各种键连参数，结构需要稍微优化下让自动判断不出错

tjuptz · 发表于 Post on 2019-9-14 16:54:45

diaok 发表于 2019-9-14 16:02
gmx给的那个rtp示例里也没看到键角二面角之类信息，应该是用ffbond里面给补全了
实现gaff应该需要把有关 ...

对于很长的聚合物，不能直接用antechamber吧

diaok · 发表于 Post on 2019-9-15 00:10:32

tjuptz 发表于 2019-9-14 16:54
对于很长的聚合物，不能直接用antechamber吧

antechamber只限制做单个分子，长了也不要紧，起码我做的800原子可以，只有些说键太多的warning，另外check检查缺失参数部分可以直接拿片段的用，不然要检查特别久。

tjuptz · 发表于 Post on 2019-9-15 08:10:21

本帖最后由 tjuptz 于 2019-9-15 10:01 编辑

diaok 发表于 2019-9-15 00:10
antechamber只限制做单个分子，长了也不要紧，起码我做的800原子可以，只有些说键太多的warning，另外che ...

用于800个原子，antechamber做AM1单点不知道我的机子做不做得了；我今天试了下我写的rtp，结果如下：
Now there are 4 residues with 194 atoms
Making bonds...
Number of bonds was 200, now 197
Generating angles, dihedrals and pairs...
Before cleaning: 514 pairs
Before cleaning: 514 dihedrals
Keeping all generated dihedrals
Making cmap torsions...
There are  514 dihedrals, 20 impropers,  364 angles
         502 pairs,    197 bonds and    0 virtual sites
Total mass 1143.680 a.m.u.
Total charge -0.590 e
Writing topology

Writing coordinate file...
            --------- PLEASE NOTE ------------
You have successfully generated a topology from: biru3.pdb.
The Gaff force field is used.

请问Number of bonds was 200, now 197 是提示出有什么问题吗？

diaok · 发表于 Post on 2019-9-15 10:57:29

tjuptz 发表于 2019-9-15 08:10
用于800个原子，antechamber做AM1单点不知道我的机子做不做得了；我今天试了下我写的rtp，结果如下：
No ...

antechamber这时候只用来生成拓扑，电荷用片段已经算过了，键的数目要对比好，不匹配肯定不行

tjuptz · 发表于 Post on 2019-9-30 20:19:59

yihanxu 发表于 2019-2-27 02:28
谢谢老师。用gview画碳链，碳原子被标记为HETATM，是为什么呢？要手动改成ATOM，感觉麻烦呢，能直接生成 ...

感觉用MS建聚合物的模型更好一些，格式更规范，而且另存时原子顺序不会乱

tjuptz · 发表于 Post on 2019-11-30 23:01:36

sobereva 发表于 2019-3-1 11:55
本来就不会被列出来，只有把xxx.ff放到了top目录下才会在选择力场时出现xxx。得把机制搞清楚
当前用GAFF ...

再次借楼请教老师，比如我有很多重复单元要补到ffbond文件中，那对于从不同重复单元获得的同样成键类型（比如 ca h）需要去重吗？对于多次定义的dihedral类型9会叠加，那bond angle呢？是用最后定义的吗？

sobereva · 发表于 Post on 2019-12-2 02:53:46

tjuptz 发表于 2019-11-30 23:01
再次借楼请教老师，比如我有很多重复单元要补到ffbond文件中，那对于从不同重复单元获得的同样成键类型（ ...

会用最后一次定义的

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

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

浏览过的版块