计算化学公社

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

[GROMACS] 酶蛋白配体复合物里结晶水的处理问题

[复制链接 Copy URL]

188

帖子

0

威望

1498

eV
积分
1686

Level 5 (御坂)

本帖最后由 uenh1998 于 2021-11-23 21:37 编辑

各位老师好,这两天在做酶蛋白催化小分子的问题。其中有一个关键的模拟是看一下小分子与酶蛋白关键位点的结合问题。从蛋白数据库中找到的结合方式是第一个图(已关闭酶蛋白Polymer显示),图中间为配体小分子。采用的是酶蛋白Amber力场+配体GAFF力场。


可以看到有几个涉及氢键作用的结晶水必须保留。请问各位老师我该如何处理这个体系呢?已经预先处理将不涉及图里的结晶水都去掉了。
PDB中结晶水与小分子均在非标准残基中,是应该将结晶水与小分子配体一起单独保存为一个PDB,然后利用acpype获得拓扑文件,include到主top文件里。

还是说要将图里面结晶水,配体,以及与结晶水形成氢键的酶蛋白相关的氨基酸一起单独保存成PDB,然后利用acpype获得拓扑文件,include到主top文件里呢?
还是说用其他的办法。请各位老师指点一下!
感谢各位老师!

111.png (62.54 KB, 下载次数 Times of downloads: 26)

111.png

头.png (45.06 KB, 下载次数 Times of downloads: 23)

头.png

尾.png (24.41 KB, 下载次数 Times of downloads: 24)

尾.png
顶天立地,树德明智,奋斗创新

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2021-11-23 08:31:11 | 只看该作者 Only view this author
结晶水和小分子又不是化学键连在一起的,当然不能合在一起用acpype来获得拓扑文件。和相关的氨基酸一起给acpype搞更是不可思议

水分子的拓扑文件本来就有现成的,诸如spce.itp。蛋白质标准氨基酸的定义直接在rtp里就有,不需要任何折腾。小分子拿acpype搞拓扑文件。
北京科音自然科学研究中心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

188

帖子

0

威望

1498

eV
积分
1686

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2021-11-23 08:42:27 | 只看该作者 Only view this author
sobereva 发表于 2021-11-23 08:31
结晶水和小分子又不是化学键连在一起的,当然不能合在一起用acpype来获得拓扑文件。和相关的氨基酸一起给ac ...

所以只需要按常规操作,把小分子的拓扑文件include到主top文件里,不需要有其他多余的操作了吧
顶天立地,树德明智,奋斗创新

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2021-11-23 09:06:17 | 只看该作者 Only view this author
uenh1998 发表于 2021-11-23 08:42
所以只需要按常规操作,把小分子的拓扑文件include到主top文件里,不需要有其他多余的操作了吧


但这几个水分子需要加氢
而且做正式模拟之前的限制性动力学的时候要把这几个水的氧原子限制住免得跑走
北京科音自然科学研究中心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

188

帖子

0

威望

1498

eV
积分
1686

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2021-11-23 21:43:15 | 只看该作者 Only view this author
sobereva 发表于 2021-11-23 09:06

但这几个水分子需要加氢
而且做正式模拟之前的限制性动力学的时候要把这几个水的氧原子限制住免得跑 ...

社长,我已经完成了酶蛋白带上那几个关键结晶水的拓扑文件的获得,配体拓扑文件已获得,配体的限制势也已经获得并include到配体itp文件里,准备把配体itp文件include到主top文件里。1楼的第2、3张图片是主top文件的开头和末尾部分。

有一个问题想再请教一下社长,之前说要把这几个氧原子也要限制住。请问在top文件开头结尾这两张图片中需要怎么进行操作才能达到目的呢?麻烦社长详细说说了
     准备下个月分子动力学Gromacs模拟课上仔细学学,但是这两天研究就要用到,所以提前请教一下社长了。
     谢谢社长!
顶天立地,树德明智,奋斗创新

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2021-11-24 08:05:43 | 只看该作者 Only view this author
不得24小时内在思想家公社QQ群和计算化学公社论坛里同一问题,注意看群文件里群规,以及论坛置顶的论坛新社员必读贴


北京科音自然科学研究中心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

188

帖子

0

威望

1498

eV
积分
1686

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2021-11-24 15:00:38 | 只看该作者 Only view this author
sobereva 发表于 2021-11-24 08:05
不得24小时内在思想家公社QQ群和计算化学公社论坛里同一问题,注意看群文件里群规,以及论坛置顶的论坛新社 ...

好的社长,一定注意!
顶天立地,树德明智,奋斗创新

188

帖子

0

威望

1498

eV
积分
1686

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2021-11-26 09:29:04 | 只看该作者 Only view this author
本帖最后由 uenh1998 于 2021-11-26 11:03 编辑
sobereva 发表于 2021-11-24 08:05
不得24小时内在思想家公社QQ群和计算化学公社论坛里同一问题,注意看群文件里群规,以及论坛置顶的论坛新社 ...

   社长您好,我将力场关于水分子的spc.itp文件复制一份改名为FIXWATER作为结晶水的itp文件include到主top文件里,且itp文件里所有的SOL改名为FIXWATER。同时紧接着他include了限制势文件water.itp。          FIXWATER.itp文件的settle部分已去掉,能量最小化能正常进行。但是用VMD看能量最小化结束后的gro文件,结构有问题。11个水分子只有一个是正常的。其余的氢氧原子都各自跑到别的地方去了。如图。

    我觉得问题可能出在FIXWATER,water这两个itp文件上。结晶水为11个,33个原子。而从力场文件拿过来的文件FIXWATER.itp里只有一个水分子。water.itp里是对11个结晶水的pdb文件生成的限制势文件。我需不需要将FIXWATER.itp里也弄成11个水分子,数量对应一下呢。还是说要改限制势文件为对单个水分子的限制势。之后我试着将FIXWATER.itp文件里的原子部分改成第二个图所示,但还是出错。

    相关文件已上传,请社长指点一下,谢谢社长。

ERROR.png (14.02 KB, 下载次数 Times of downloads: 39)

ERROR.png

改.png (68.53 KB, 下载次数 Times of downloads: 30)

改.png

FIXWATER.itp

702 Bytes, 下载次数 Times of downloads: 8

water.itp

1.55 KB, 下载次数 Times of downloads: 3

顶天立地,树德明智,奋斗创新

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

9#
发表于 Post on 2021-11-26 11:02:42 | 只看该作者 Only view this author
uenh1998 发表于 2021-11-26 09:29
社长您好,我将力场关于水分子的spc.itp文件复制一份改名为FIXWATER作为结晶水的itp文件include到主t ...

带有限制势字段[position_restraint]的文件是在被限制的分子的[moleculetype]之后引入的。[moleculetype]里的[atoms]里有多少个原子,你若想限制所有这类分子的原子,[position_restraint]里就要有多少原子。
北京科音自然科学研究中心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

188

帖子

0

威望

1498

eV
积分
1686

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2021-11-26 11:09:10 | 只看该作者 Only view this author
本帖最后由 uenh1998 于 2021-11-27 17:03 编辑
sobereva 发表于 2021-11-26 11:02
带有限制势字段的文件是在被限制的分子的[moleculetype]之后引入的。[moleculetype]里的[atoms]里有多少 ...

社长,water.itp里是我单独把11个水分子拿出来生成的限制文件。FIXWATER.itp是力场文件的单个水分子的itp文件。我把FIXWATER.itp原子部分改成楼上最下方那个图了。。我是直接复制粘贴并改了一下序号,发现还是出现氢氧原子各自跑到一边去的错误。具体我也没有查到这个itp的文件格式该怎么写。请问社长我这个图里的写法该怎么改呢。


顶天立地,树德明智,奋斗创新

3

帖子

0

威望

67

eV
积分
70

Level 2 能力者

11#
发表于 Post on 2024-8-23 09:29:46 | 只看该作者 Only view this author
uenh1998 发表于 2021-11-26 11:09
社长,water.itp里是我单独把11个水分子拿出来生成的限制文件。FIXWATER.itp是力场文件的单个水分子的itp ...

前辈您好,我最近也在做酶-配体-底物-结晶水的gromacs模拟,结晶水也需要保留。请问如何对结晶水进行位置约束?可以求您之前的水的限制势文件和用于做模拟的complex.gro和topol.top文件吗?我想对着修改一下。非常感谢您。

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 05:43 , Processed in 0.169505 second(s), 23 queries , Gzip On.

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