24小时热门版块排行榜    

查看: 1479  |  回复: 3

Jay_LIN

新虫 (初入文坛)

[求助] 求教Gromacs计算蛋白-蛋白相互作用的top文件建立方法

各位大神好。我最近在计算蛋白-蛋白相互作用时碰到了top文件的建立问题。图和一图二分别是我的蛋白-蛋白复合体系的top文件的头部与底部。我的想法就是将配体蛋白的itp(LEA-LEA.itp)include进top中。

图一

图二

图三是LEA-LEA.itp,我是用gromacs的pdb2gmx指令生成的。

图三

gromacs在进行能量最小化时报出如下错误

可以看出,gromacs并没有识别配体蛋白质的力场文件而是在读取了主体蛋白质的力场后又继续读取主体蛋白质的力场。之前做小分子配体的时候这个建力场方法是行得通的,但配体换成蛋白就出问题了。请问我该如何建立蛋白-蛋白复合体系的top呢?(本人新手小白,第一次提问,问题提的肯能有些不太成熟,请见谅。金币也不多,希望各位大大不要嫌弃 谢谢大家了!)

求教Gromacs计算蛋白-蛋白相互作用的top文件建立方法
蛋白-蛋白复合体系top头部.png


求教Gromacs计算蛋白-蛋白相互作用的top文件建立方法-1
蛋白-蛋白top末尾.png


求教Gromacs计算蛋白-蛋白相互作用的top文件建立方法-2
LEA-LEAitp.png


求教Gromacs计算蛋白-蛋白相互作用的top文件建立方法-3
gromacs能量最小化.png
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hclwonder

新虫 (正式写手)

不知你想研究什么相互作用呢?是计算结合能?那需要index文件设置一下以分开受体配体。
但只是想跑跑蛋白-蛋白复合体的分子动力学,完全不需要任何修改top文件。就是说,跟跑单个蛋白一样的处理方式。

发自小木虫Android客户端
2楼2017-04-15 08:44:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jay_LIN

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by hclwonder at 2017-04-15 08:44:52
不知你想研究什么相互作用呢?是计算结合能?那需要index文件设置一下以分开受体配体。
但只是想跑跑蛋白-蛋白复合体的分子动力学,完全不需要任何修改top文件。就是说,跟跑单个蛋白一样的处理方式。
...

谢谢!也就是说把蛋白质主体与蛋白质客体放到同一个pdb文件中,当作一个蛋白用gromacs处理。是这样吧?
3楼2017-04-15 21:14:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hclwonder

新虫 (正式写手)

引用回帖:
3楼: Originally posted by Jay_LIN at 2017-04-15 21:14:39
谢谢!也就是说把蛋白质主体与蛋白质客体放到同一个pdb文件中,当作一个蛋白用gromacs处理。是这样吧?...

是啊,我是蛋白复合体在一个pdb文件的。

发自小木虫Android客户端
4楼2017-04-15 22:06:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Jay_LIN 的主题更新
信息提示
请填处理意见