计算化学公社

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

[GROMACS] 有关利用Sobtop生成过渡金属蛋白拓扑的问题

[复制链接 Copy URL]

17

帖子

0

威望

201

eV
积分
218

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 pupukoNANA 于 2023-9-20 12:20 编辑

前贴见:如何进行金属酶与底物和金属离子(过渡金属)的模拟?
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=39051&fromuid=51838

方法参考:http://sobereva.com/soft/Sobtop#FAQ7

当前已经完成:利用sobtop对锰离子及配位残基拓扑文件的生成
center.gjf (7.61 KB, 下载次数 Times of downloads: 33) center.itp (42.71 KB, 下载次数 Times of downloads: 25)
(对于Gaussian我完全是小白,若有不当请指出,感谢)
恕学生愚昧,对于生成拓扑之后的工作反复看了试了几天还是没搞成,主要有以下几个问题

1、“
然后,用Sobtop对铁硫簇部分产生参数为空的itp文件。”
在我的体系里面,是指单独对锰离子生成拓扑吗?

2、“之后,用pdb2gmx对纯蛋白质部分产生[ moleculetype ],将之和铁硫簇部分的[ moleculetype ]都引入主top里”

其中:“纯蛋白部分”指的是否是簇模型中的蛋白部分?            “主top”是否是 全 金属蛋白的纯蛋白部分经过pdb2gmx得到的拓扑?

3、“然后在top里加入[ intermolecular_interactions ]字段、、、、、、、”
该步前提是否是:在 全 纯 蛋白生成拓扑之后,对其gro文件手动添加金属离子



GROMACS reminds you: "Focus on science, sublimate youself." (Rick Sanchez)


6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

2#
发表于 Post on 2023-9-20 03:29:32 | 只看该作者 Only view this author
没有叫soptop的东西
关键词只能写opt freq,不能写opt-freq
去掉所有(PDBName...)部分

你的当前情况远比FAQ7说的铁硫簇的情况简单,直接通过[ intermolecular_interactions ]字段给蛋白质和Mn离子这两个[moleculetype]之间加上成键相关的项就完了,其中参数来自于你当前的簇模型用sobtop基于Hessian计算出来的
北京科音自然科学研究中心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

17

帖子

0

威望

201

eV
积分
218

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-9-20 12:24:24 | 只看该作者 Only view this author
本帖最后由 pupukoNANA 于 2023-9-20 23:27 编辑
sobereva 发表于 2023-9-20 03:29
没有叫soptop的东西
关键词只能写opt freq,不能写opt-freq
去掉所有(PDBName...)部分

感谢sob老师的回答。
已经改正帖子中错误的拼写。

17

帖子

0

威望

201

eV
积分
218

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-9-20 23:48:27 | 只看该作者 Only view this author
sobereva 发表于 2023-9-20 03:29
没有叫soptop的东西
关键词只能写opt freq,不能写opt-freq
去掉所有(PDBName...)部分

抱歉再次打扰老师


[ intermolecular_interactions ]字段参数已经修改完。 intermolecular_interactions.txt (4.11 KB, 下载次数 Times of downloads: 50)
我经过多次尝试,拓扑文件的设置总是不对,报错:Invalid order for directive intermolecular_interactions. 请问这个错误的原因和有效的解决方案。
此外,学生不清楚“将之和铁硫簇部分的[ moleculetype ]都引入主top”里,“引入”的具体做法为何。请问是直接在主拓扑里写[ moleculetype ],又或还需要其它文件?
希望老师不吝赐教。

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

5#
发表于 Post on 2023-9-21 01:35:54 | 只看该作者 Only view this author
pupukoNANA 发表于 2023-9-20 23:48
抱歉再次打扰老师

.txt里的内容没插入到top里合适的地方。没完整文件没法说

把itp 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

17

帖子

0

威望

201

eV
积分
218

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-9-21 08:45:16 | 只看该作者 Only view this author
本帖最后由 pupukoNANA 于 2023-9-21 20:04 编辑
sobereva 发表于 2023-9-21 01:35
.txt里的内容没插入到top里合适的地方。没完整文件没法说

把itp include进去,或者直接把内容拷到top ...

感谢老师百忙之中抽空回复学生的问题

已经从gmx官网手册了解到,[size=15.008px][ intermolecular_interactions ]字段需放在拓扑文件最后

14

帖子

0

威望

277

eV
积分
291

Level 3 能力者

7#
发表于 Post on 2023-12-26 22:10:14 | 只看该作者 Only view this author
pupukoNANA 发表于 2023-9-21 08:45
感谢老师百忙之中抽空回复学生的问题

已经从gmx官网手册了解到,[ intermolecular_interactions ]字段 ...

你好可以看一下你的top文件怎么设置的吗?

27

帖子

0

威望

347

eV
积分
374

Level 3 能力者

8#
发表于 Post on 2023-12-28 06:52:49 | 只看该作者 Only view this author
pupukoNANA 发表于 2023-9-21 08:45
感谢老师百忙之中抽空回复学生的问题

已经从gmx官网手册了解到,[ intermolecular_interactions ]字段 ...

能把你成功的top文件和itp文件发上来参考一下吗?需要为gro文件添加金属离子吗

39

帖子

0

威望

219

eV
积分
258

Level 3 能力者

9#
发表于 Post on 2025-8-30 02:41:07 | 只看该作者 Only view this author
sobereva 发表于 2023-9-20 03:29
没有叫soptop的东西
关键词只能写opt freq,不能写opt-freq
去掉所有(PDBName...)部分

老师,我有一个小问题,我用簇模型用高斯去计算酶中心(又一个锌离子和三个与之配位的氨基酸),我的初始结构使用的af3预测的蛋白和锌离子的复合物结构,截取锌离子周围5A的几个残基进行了opt freq计算,计算后的结构和初始结构发生了一定的变化,那么请问老师我是应该将这些计算后的氨基酸的构象替换初始预测的结构的构象吗?如果不是这么做的话又如何通过[ intermolecular_interactions ]字段给蛋白质和锌离子这两个[moleculetype]之间加上成键相关的项呢?感谢!

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

10#
发表于 Post on 2025-8-30 05:20:41 | 只看该作者 Only view this author
wang011127 发表于 2025-8-30 02:41
老师,我有一个小问题,我用簇模型用高斯去计算酶中心(又一个锌离子和三个与之配位的氨基酸),我的初始 ...

AF3直接给你预测的当然没有优化之后的准确,自然不应该用“初始预测的结构的构象”
北京科音自然科学研究中心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

39

帖子

0

威望

219

eV
积分
258

Level 3 能力者

11#
发表于 Post on 2025-8-30 22:46:33 | 只看该作者 Only view this author
sobereva 发表于 2025-8-30 05:20
AF3直接给你预测的当然没有优化之后的准确,自然不应该用“初始预测的结构的构象”

明白了,感谢!那么请问如何将优化后的簇和原来的蛋白“拼”在一起呢?需要用sobtop对簇生成拓扑然后替换到原蛋白中的残基的拓扑吗?

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

12#
发表于 Post on 2025-8-31 07:49:08 | 只看该作者 Only view this author
wang011127 发表于 2025-8-30 22:46
明白了,感谢!那么请问如何将优化后的簇和原来的蛋白“拼”在一起呢?需要用sobtop对簇生成拓扑然后替换 ...

根本就不需要拼结构
本来动力学跑起来之后结构该怎么变就怎么变,初始结构并不重要
配体部分的拓扑信息用sobtop产生的
北京科音自然科学研究中心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

39

帖子

0

威望

219

eV
积分
258

Level 3 能力者

13#
发表于 Post on 2025-8-31 12:38:23 | 只看该作者 Only view this author
sobereva 发表于 2025-8-31 07:49
根本就不需要拼结构
本来动力学跑起来之后结构该怎么变就怎么变,初始结构并不重要
配体部分的拓扑信息 ...

感谢回复!我再理解了一下,相当于是我对于和锌离子成键的几个氨基酸在原蛋白里加一些 intermolecular_interactions,参数利用sobtop获得,这个参数包括键长(是否包括键角和二面角以及我是否应该添加键角和二面角的信息?),从而用于限制残基和锌离子的配位形态和距离。

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

14#
发表于 Post on 2025-9-1 07:05:53 | 只看该作者 Only view this author
wang011127 发表于 2025-8-31 12:38
感谢回复!我再理解了一下,相当于是我对于和锌离子成键的几个氨基酸在原蛋白里加一些 intermolecular_in ...

此字段里加入[bonds]、[angles]项,平衡键长、键角用优化完的结构的,力常数随便猜一个足够大的就完了
北京科音自然科学研究中心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

18

帖子

0

威望

129

eV
积分
147

Level 2 能力者

15#
发表于 Post on 2025-12-12 12:35:21 | 只看该作者 Only view this author
wang011127 发表于 2025-8-30 02:41
老师,我有一个小问题,我用簇模型用高斯去计算酶中心(又一个锌离子和三个与之配位的氨基酸),我的初始 ...

大佬,您好,我是一个刚刚接触这方面的小白,也在做锌指蛋白,请问能否求下您做分析时的gif文件、top和itp文件呢?

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 13:15 , Processed in 0.214168 second(s), 22 queries , Gzip On.

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