使用gromacs模拟蛋白配体体系时，关于配体限制势文件在topol.top的位置问题？

12313 · 发表于 Post on 2024-7-30 12:40:49

本帖最后由 12313 于 2024-7-30 14:38 编辑

各位老师好，我要模拟的体系是蛋白配体复合物体系，配体是多肽小分子，我使用pdb2gmx命令，通过输入配体和受体蛋白的pdb文件，分别产生了各自的gro、itp、top文件。之后我将配体gro文件的内容复制到受体的gro文件并将原子数进行了相加，在受体的topol.top文件中引入限制势文件，如图所示。但是在使用gmx grompp进行添加离子的操作时，出现如图的错误，显示没有Protein_chain_B这个分子类型，请问该如何解决？（图中LPPR.itp是配体的限制文件，posre.itp是受体蛋白的限制文件，chainA是受体蛋白，chainB是多肽配体）

Loading0760 · 发表于 Post on 2024-7-30 13:23:06

图呢?

KanbeKotori · 发表于 Post on 2024-7-30 13:56:27

其他的先不管，既然显示没有这个分子类型，那就去蛋白拓扑文件里检查。每种分子类型都必须有先前定义的其对应部分，所以需要在蛋白拓扑文件中引用配体的itp文件以及到[molecules]添加分子

12313 · 发表于 Post on 2024-7-30 14:38:47

Loading0760 发表于 2024-7-30 13:23
图呢?

抱歉抱歉，是我的疏忽，已经添加了

Seyilaxa · 发表于 Post on 2024-7-30 16:04:26

本帖最后由 Seyilaxa 于 2024-7-30 16:09 编辑

图里没看到是否引用了配体的itp文件

12313 · 发表于 Post on 2024-7-30 19:26:14

Seyilaxa 发表于 2024-7-30 16:04
图里没看到是否引用了配体的itp文件

有添加，老师。图中LPPR.itp就是我的配体itp文件

student0618 · 发表于 Post on 2024-7-30 20:04:54

那是位置限制文件还是配体拓朴？一楼指是限制文件，但5楼问的是拓朴吧。有没有include它的拓朴？
如果配体的力场文件不应放在 #ifdef POSRES 部分哦，检查一下它是什么。

建议多了解一下Gromacs topology文件结构，每个报错指那部份；问问题时最好给完整topol.top (有机密的话可涂掉一部份吧)不然其他人很难帮忙的。

12313 · 发表于 Post on 2024-7-30 21:27:49

本帖最后由 12313 于 2024-7-30 21:33 编辑

student0618 发表于 2024-7-30 20:04
那是位置限制文件还是配体拓朴？一楼指是限制文件，但5楼问的是拓朴吧。有没有include它的拓朴？
如果配体 ...

谢谢老师，这个LPPR.itp是通过gmx pdb2gmx指令输入配体的pdb文件后产生的，所以它是位置限制文件吗？因为我总是把它和gmx genrestr产生的itp文件搞混

student0618 · 发表于 Post on 2024-7-30 21:47:42

12313 发表于 2024-7-30 21:27
谢谢老师，这个LPPR.itp是通过gmx pdb2gmx指令输入配体的pdb文件后产生的，所以它是位置限制文件吗？因为 ...

不清楚那是什麼可打开那itp來看看。

12313 · 发表于 Post on 2024-7-30 22:16:51

student0618 发表于 2024-7-30 21:47
不清楚那是什麼可打开那itp來看看。

老师，我已上传

Seyilaxa · 发表于 Post on 2024-7-30 23:17:11

本帖最后由 Seyilaxa 于 2024-7-30 23:18 编辑

12313 发表于 2024-7-30 22:16
老师，我已上传

很明显这是限制文件，你在1L也说了“图中LPPR.itp是配体的限制文件”
gmx pdb2gmx会同时产生两个文件，一个是配体的拓扑文件，一个是限制文件，两个东西不能混为一谈

12313 · 发表于 Post on 2024-7-31 11:39:32

Seyilaxa 发表于 2024-7-30 23:17
很明显这是限制文件，你在1L也说了“图中LPPR.itp是配体的限制文件”
gmx pdb2gmx会同时产生两个文件， ...

谢谢老师，我明白了。那请问我该如何进行修改呢？

Seyilaxa · 发表于 Post on 2024-7-31 12:02:14

12313 发表于 2024-7-31 11:39
谢谢老师，我明白了。那请问我该如何进行修改呢？

把配体的拓扑文件include上就行

12313 · 发表于 Post on 2024-7-31 15:08:25

Seyilaxa 发表于 2024-7-31 12:02
把配体的拓扑文件include上就行

老师，我已经按照正确的顺序在topol.top文件中进行了引用（在LPPR.top的结尾我已include位置限制文件LPPR.itp），但是还是出现了Invalid order for directive defaults，如图所示，请问我哪里出错了呢？

Seyilaxa · 发表于 Post on 2024-7-31 16:04:09

12313 发表于 2024-7-31 15:08
老师，我已经按照正确的顺序在topol.top文件中进行了引用（在LPPR.top的结尾我已include位置限制文件LPPR ...

你看看是不是LPPR.top文件里有[defaults]字段没有删掉，就用这个帖子里的格式来检查http://bbs.keinsci.com/thread-45783-1-1.html
另外如果受体和配体都是多肽和蛋白，没必要在不理解拓扑文件各字段的情况下给配体和受体蛋白分别创建拓扑文件再合并，gmx pdb2gmx可以给蛋白复合物直接创建拓扑文件

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

[GROMACS] 使用gromacs模拟蛋白配体体系时，关于配体限制势文件在topol.top的位置问题？