计算化学公社

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

[GROMACS] 求助含有锌离子的蛋白质的分子动力学模拟后蛋白质结构完全散掉

[复制链接 Copy URL]

152

帖子

0

威望

888

eV
积分
1040

Level 4 (黑子)

本帖最后由 Loading0760 于 2024-9-11 11:50 编辑

老师们好,我最近在计算含有金属离子(锌, +2价)的蛋白质的分子动力学,使用的软件的GROMACS2024,利用CHARMM-GUI的Solution Builder构建的模型,没有进行多余的处理,就直接跑production部分了.
下图是100 ns模拟后的结构对比,去除水和中和离子, 右图为100 ns后的结构,可以看到蛋白质结构已经完全散掉了.
想请问一下各位老师, 这种情况应该怎么解决, 我个人初步分析是由于局部电荷初始的时候太集中了,简而言之,4个锌离子加起来的总电荷高达+8,导致后续模拟蛋白质结构破坏.



图中,
A 初始结构
B step4.0_minimization后结构
C step4.1_equilibration后结构
D 100 ns的md后的结构



蛋白质5WLM :绿色为氧与锌的连线,距离2.3 Å ;黄色为ASP的羧基氧与锌的连线,距离1.9 Å
线为本人自己画的,仅供可视化时参考,实际并不成键.

849

帖子

4

威望

1918

eV
积分
2847

Level 5 (御坂)

A Student

2#
发表于 Post on 2024-9-10 17:35:13 | 只看该作者 Only view this author
先检查一下和Zn离子作用的残基,会不会是加氢加到coordination的地方令em时跑开了?
敬仰一针见血的指责,厌倦别有用心的赞美。

152

帖子

0

威望

888

eV
积分
1040

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-9-10 17:54:27 | 只看该作者 Only view this author
student0618 发表于 2024-9-10 17:35
先检查一下和Zn离子作用的残基,会不会是加氢加到coordination的地方令em时跑开了?

您好,补充了不同时间段的结构

与Zn离子作用的氨基酸应该是HIS,将其质子化状态考虑为HSE,也就是δ位置的氮原子被质子化,与Zn离子更近的氮原子未质子化

96

帖子

0

威望

2571

eV
积分
2667

Level 5 (御坂)

4#
发表于 Post on 2024-9-10 18:18:53 | 只看该作者 Only view this author
考虑使用约束模型或者使用sobtop/MCPB构建金属离子成键模型

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

5#
发表于 Post on 2024-9-11 00:00:59 | 只看该作者 Only view this author
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并清楚、准确反映出帖子具体内容,避免有任何歧义和含糊性,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “含有金属离子的蛋白质的分子动力学模拟” 改了,以后务必注意

弄清楚初始的模型是否靠谱、Zn2+周围残基的质子化状态、溶剂环境这些问题。如果没有负电荷产生抗衡作用,四个Zn2+挤在一起必然容易分离开。必要时可以先加一些限制/约束跑一阵子然后再放开
北京科音自然科学研究中心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

152

帖子

0

威望

888

eV
积分
1040

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2024-9-11 11:30:39 | 只看该作者 Only view this author
fffff 发表于 2024-9-10 18:18
考虑使用约束模型或者使用sobtop/MCPB构建金属离子成键模型

考虑使用约束模型吧,成键模型不太合适,因为没有成键.

152

帖子

0

威望

888

eV
积分
1040

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-9-11 11:42:10 | 只看该作者 Only view this author
sobereva 发表于 2024-9-11 00:00
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此 ...

好的老师,第一次发帖,下次一定注意.
我看了一下相关的文献,PDBID:5WLM. 他是四个锌周围有ASP,还有游离氧,示意图已在提问中补充.这两个的作用使得局部电荷为0.如果想构建这样的体系,游离氧在分子动力学里怎么描述 ?是否需要QM/MM呢?

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

8#
发表于 Post on 2024-9-11 11:50:56 | 只看该作者 Only view this author
Loading0760 发表于 2024-9-11 11:42
好的老师,第一次发帖,下次一定注意.
我看了一下相关的文献,PDBID:5WLM. 他是四个锌周围有ASP,还有游离氧 ...

游离氧可能主要是氢氧根形式存在,可以再查查相关资料,这种情况用sobtop做一个OH-的拓扑文件就行
北京科音自然科学研究中心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

152

帖子

0

威望

888

eV
积分
1040

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2024-9-13 10:44:37 | 只看该作者 Only view this author
sobereva 发表于 2024-9-11 11:50
游离氧可能主要是氢氧根形式存在,可以再查查相关资料,这种情况用sobtop做一个OH-的拓扑文件就行

不太行啊,社长。这个真的就是二价氧,那个区域,两个锌+8价的电荷。
四个谷氨酸,其以阴离子形式存在(-4),加上两个二价氧(-4),刚好是中和了,蛋白质就稳定。感觉是不是得上QM来描述这个区域了。

6万

帖子

99

威望

5万

eV
积分
124691

管理员

公社社长

10#
发表于 Post on 2024-9-13 15:48:19 | 只看该作者 Only view this author
Loading0760 发表于 2024-9-13 10:44
不太行啊,社长。这个真的就是二价氧,那个区域,两个锌+8价的电荷。
四个谷氨酸,其以阴离子形式存在( ...

如果是二价氧,那就不算游离了,而是桥联。这需要结合bond项把Zn和O部分当做一个准刚性的簇对待
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 00:38 , Processed in 0.172926 second(s), 23 queries , Gzip On.

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