计算化学公社

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

[GROMACS] MD RMSD分析结构稳定问题

[复制链接 Copy URL]

24

帖子

0

威望

870

eV
积分
894

Level 4 (黑子)

本帖最后由 yuanen 于 2021-6-3 15:30 编辑

本人GROMACS入门小白,想请教一下做金属与有机小分子结合(练手)通过GROMACS的rms命令得到RMSD图如下:
1,关于有机分子的结构是否稳定,虽然在一定值上下波动,但感觉波动范围偏大(0.5埃)
2,将有机小分子和金属写在一个组内整体的RMSD,波动范围更大了(6埃)3,有机小分子(33个原子)通过编辑参数添加到CHARMM27中
请问这是正常的吗?是什么因素引起的?
不太清楚,还请帮忙解答






有机小分子RMSD.jpg (87.68 KB, 下载次数 Times of downloads: 44)

有机分子的RMSD

有机分子的RMSD

有机小分子与金属RMSD.jpg (79.53 KB, 下载次数 Times of downloads: 40)

有机分子与金属的RMSD

有机分子与金属的RMSD

71

帖子

0

威望

1148

eV
积分
1219

Level 4 (黑子)

2#
发表于 Post on 2021-6-3 19:20:37 | 只看该作者 Only view this author
你的体系没有说清楚。如果你的体系为多个一种小分子,想研究它在某金属溶液(确定金属的浓度)中与某金属的结合情况。你可以测金属在小分子周围分布的径向函数。在MD中,如果没有对金属做约束,那么这个离子很显然是个自由移动状态。

6万

帖子

99

威望

6万

eV
积分
125117

管理员

公社社长

3#
发表于 Post on 2021-6-4 09:12:11 | 只看该作者 Only view this author
哪怕是刚性有机小分子,由于热运动,结构自然也会有波动。在VMD里对小分子做align,把多帧同时显示出来就可以清楚地看出小分子的波动情况,比看RMSD直观得多

“有机小分子和金属写在一个组内整体的RMSD”这没有意义,本来金属就可以在小分子附近倒处晃悠。除非你是指配合物那种金属和配体位置关系保持固定的情况。
北京科音自然科学研究中心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

24

帖子

0

威望

870

eV
积分
894

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2021-6-4 09:21:54 | 只看该作者 Only view this author
本帖最后由 yuanen 于 2021-6-4 09:41 编辑
panernie 发表于 2021-6-3 19:20
你的体系没有说清楚。如果你的体系为多个一种小分子,想研究它在某金属溶液(确定金属的浓度)中与某金属的 ...

感谢回答,现在只是一个有机分子和一个金属离子在水模型下跑MD,我是通过直接将金属离子摆放在可能的结合位点,通过观察轨迹文件我发现金属离子与有机小分子是同步移动的。
关于对金属做约束,是什么约束?

24

帖子

0

威望

870

eV
积分
894

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-6-4 09:40:43 | 只看该作者 Only view this author
sobereva 发表于 2021-6-4 09:12
哪怕是刚性有机小分子,由于热运动,结构自然也会有波动。在VMD里对小分子做align,把多帧同时显示出来就可 ...

感谢回答 也就是说关于有机小分子的波动可以暂时看做是合理的,我现在就是直接做的金属与有机分子的配合物这种所以才有这样的疑惑,是缺少什么设置吗?

6万

帖子

99

威望

6万

eV
积分
125117

管理员

公社社长

6#
发表于 Post on 2021-6-5 10:15:19 | 只看该作者 Only view this author
yuanen 发表于 2021-6-4 09:40
感谢回答 也就是说关于有机小分子的波动可以暂时看做是合理的,我现在就是直接做的金属与有机分子的配合 ...

我不知道你说的设置是什么
看轨迹,根据基本常识解释计算出的和结构有关的数据
北京科音自然科学研究中心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

24

帖子

0

威望

870

eV
积分
894

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2021-6-6 18:56:19 | 只看该作者 Only view this author
sobereva 发表于 2021-6-5 10:15
我不知道你说的设置是什么
看轨迹,根据基本常识解释计算出的和结构有关的数据

好的,我再查考一下看看

7

帖子

0

威望

151

eV
积分
158

Level 3 能力者

8#
发表于 Post on 2022-6-27 21:03:32 | 只看该作者 Only view this author
请问各位老师,模拟体系的平均RMSD值,是否就是得到的rmsd.xvg文件最后一行的数值呢?

310

帖子

0

威望

1503

eV
积分
1813

Level 5 (御坂)

9#
发表于 Post on 2022-6-27 21:08:45 | 只看该作者 Only view this author
xingbb 发表于 2022-6-27 21:03
请问各位老师,模拟体系的平均RMSD值,是否就是得到的rmsd.xvg文件最后一行的数值呢?

不是,最后一行是最后一帧的rmsd

7

帖子

0

威望

151

eV
积分
158

Level 3 能力者

10#
发表于 Post on 2022-6-28 08:52:16 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2022-6-27 21:08
不是,最后一行是最后一帧的rmsd

好的,谢谢您!请问我想要得到平均rmsd是需要自己去拟合吗

8

帖子

0

威望

33

eV
积分
41

Level 2 能力者

11#
发表于 Post on 2023-3-13 18:54:08 | 只看该作者 Only view this author
大佬您好,我目前是一名gromacs小白,想用二价铁离子和铁蛋白跑gromacs,想和您请教下相关操作。

本版积分规则 Credits rule

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

GMT+8, 2026-2-16 21:38 , Processed in 0.182298 second(s), 23 queries , Gzip On.

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