计算化学公社

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

[GROMACS] RMSD计算问题

[复制链接 Copy URL]

7

帖子

0

威望

87

eV
积分
94

Level 2 能力者

各位好!我在md完之后的蛋白和配体RMSD计算时遇到了一些问题。
体系是蛋白-配体复合物,为了节约一些时间,盒子选的是十二面体盒子。
但是在计算RMSD时,用vmd计算的结果和用rms计算的结果差的很多,不知道是在哪一步操作上有问题。
希望各位帮助,给出一个详细的计算步骤。
我之前的操作如下:
保持分子完整和去除周期性
gmx trjconv -f ../1kv1_lig_md.xtc -s ../1kv1_lig_md.tpr -pbc mol -o new.xtc
gmx trjconv -f new.xtc -pbc nojump -o nojump.xtc
用vmd->Extensions->Analysis->RMSD Trajectory tool
分别对proteion和resname BMU进行Align->RMSD,得到网盘中1和2两张图
之后用分别用以下两个命令得到蛋白和配体的RMSD
gmx rms -f nojump.xtc -s ../1kv1_lig_md.tpr -o rmsd_protein_GMX.xvg | 4(backbone) | 4(backbone)(网盘中图3)
gmx rms -f nojump.xtc -s ../1kv1_lig_md.tpr -o rmsd_ligand_GMX.xvg | 4(backbone) | 20(BMU))(网盘中图4)
但是和之前用vmd的图差了好多
后来看到之前也有人问这个问题,又将backbone和BMU合并成了一个组24。然后选24|4和24|20分别求RMSD也得到相同的结果。
请高人指点一下哪里出了问题,我对这块一直不太明白。
由于东西较大,都放网盘里了,麻烦大家下载一下,谢谢!
链接:https://pan.baidu.com/s/1c2LIyYG 密码:5dzy

6万

帖子

99

威望

6万

eV
积分
125125

管理员

公社社长

2#
发表于 Post on 2017-12-20 12:50:13 | 只看该作者 Only view this author
最近太忙没时间仔细看,建议你就一律用VMD来计算就完了,VMD里直接给出的曲线和你VMD里看到的完全一样,所见即所得,想指定align哪个组也很方便,RMSD曲线如果有问题则VMD里看到的轨迹肯定也有问题,很容易弄清楚原因。我绘制RMSD基本都用VMD
北京科音自然科学研究中心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

7

帖子

0

威望

87

eV
积分
94

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2017-12-20 14:57:05 | 只看该作者 Only view this author
sobereva 发表于 2017-12-20 12:50
最近太忙没时间仔细看,建议你就一律用VMD来计算就完了,VMD里直接给出的曲线和你VMD里看到的完全一样,所 ...

好的,谢谢!

156

帖子

0

威望

1653

eV
积分
1809

Level 5 (御坂)

4#
发表于 Post on 2017-12-20 17:33:42 | 只看该作者 Only view this author
我估计也是pbc没有处理好吧,你可以在vmd里面看看你的轨迹有没有问题,如果有问题,可能就会出现很大的RMSD值。

我的猜测是,在某些帧,小分子跑出了蛋白,而你是计算的蛋白+小分子的RMSD,你可以试着单独计算蛋白与小分子各自的RMSD看看。

7

帖子

0

威望

87

eV
积分
94

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2017-12-20 22:32:25 | 只看该作者 Only view this author
beyond 发表于 2017-12-20 17:33
我估计也是pbc没有处理好吧,你可以在vmd里面看看你的轨迹有没有问题,如果有问题,可能就会出现很大的RMSD ...

-pbc mol和nojump后的轨迹我从头到尾看过一遍了,轨迹应该是没问题,小分子也没跑出蛋白外,我觉得可能是我后续操作的哪个地方不太对,用VMD算的RMSD看起来正常,和原文献比也基本一致,就是到了用gmx算的时候就是算不对。谢谢啦!

334

帖子

0

威望

2357

eV
积分
2691

Level 5 (御坂)

6#
发表于 Post on 2017-12-20 22:51:10 来自手机 | 只看该作者 Only view this author
你用-pbc mol之后为什么还要再来一次-pbc nojump?

6万

帖子

99

威望

6万

eV
积分
125125

管理员

公社社长

7#
发表于 Post on 2017-12-20 23:57:35 | 只看该作者 Only view this author
霜晨月 发表于 2017-12-20 22:51
你用-pbc mol之后为什么还要再来一次-pbc nojump?

-pbc mol只不过是让轨迹按照PBC记录的时候分子保持完整,nojump是完全去除周期性
北京科音自然科学研究中心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

156

帖子

0

威望

1653

eV
积分
1809

Level 5 (御坂)

8#
发表于 Post on 2017-12-23 01:47:14 | 只看该作者 Only view this author
sobereva 发表于 2017-12-20 23:57
-pbc mol只不过是让轨迹按照PBC记录的时候分子保持完整,nojump是完全去除周期性

我用pbc nujump得到的轨迹,很多时候都是一团乱麻

6万

帖子

99

威望

6万

eV
积分
125125

管理员

公社社长

9#
发表于 Post on 2017-12-23 02:35:27 | 只看该作者 Only view this author
beyond 发表于 2017-12-23 01:47
我用pbc nujump得到的轨迹,很多时候都是一团乱麻


可能是你载入的第一个结构不合适。VMD默认情况下会按照第一个载入的结构中的原子距离来判断连接关系
北京科音自然科学研究中心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

90

帖子

0

威望

330

eV
积分
420

Level 3 能力者

10#
发表于 Post on 2019-7-1 15:35:53 | 只看该作者 Only view this author
请问我想计算MOL的RMSD的话,参考组应该也是选择MOL吗?谢谢

生活需要勇气,一定要坚持下去!

6万

帖子

99

威望

6万

eV
积分
125125

管理员

公社社长

11#
发表于 Post on 2019-7-2 11:13:42 | 只看该作者 Only view this author
青青青 发表于 2019-7-1 15:35
请问我想计算MOL的RMSD的话,参考组应该也是选择MOL吗?谢谢

看你目的
如果你只是想考察这个分子相对于自己的某一帧结构的相对位置波动情况,就两次都选择MOL
北京科音自然科学研究中心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

90

帖子

0

威望

330

eV
积分
420

Level 3 能力者

12#
发表于 Post on 2019-7-5 09:44:36 | 只看该作者 Only view this author
sobereva 发表于 2019-7-2 11:13
看你目的
如果你只是想考察这个分子相对于自己的某一帧结构的相对位置波动情况,就两次都选择MOL

好的,谢谢老师,请问您说的相对于“某一帧结构”,这里的“某一帧"可以自行规定吗?谢谢老师
生活需要勇气,一定要坚持下去!

6万

帖子

99

威望

6万

eV
积分
125125

管理员

公社社长

13#
发表于 Post on 2019-7-5 13:02:50 | 只看该作者 Only view this author
青青青 发表于 2019-7-5 09:44
好的,谢谢老师,请问您说的相对于“某一帧结构”,这里的“某一帧"可以自行规定吗?谢谢老师

可以
北京科音自然科学研究中心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-2-18 06:18 , Processed in 0.185940 second(s), 23 queries , Gzip On.

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