计算化学公社

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

[GROMACS] 求助配受体复合物RMSD等值的分析问题

[复制链接 Copy URL]

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

大家好 本人是一名在读小硕 因为自身没有相关知识背景 特来求助各位前辈有关MD结果分析的一些问题。
我想做的课题是:验证一个药物分子和三种不同蛋白是否能结合 如果能则筛选出一个结合最紧密的一个蛋白(单纯地只考虑是否能结合,然后用实验去验证)。
因此我MD后做了RMSD、Rg、RMSF以及平均自由结合能的计算。部分结果如下图所示,结合能依次是-146、-114、-108(kj/mol).
我想问的是:1.从RMSD值看三个复合物的变化范围都不大 且值接近 我要怎么定性地描述这个结果(一般看文献中RMSD值是一开始会变化明显 然后逐渐稳定)。
                   2.仅仅从结合能是否能得出结合得最紧密的就是第一个复合物,还是说要把各项分析结合起来得出结论。
                   3.回到最初的目的,我做的这些分析是否能验证三种复合物的结合,会不会做得不够或者多余的分析。
特此向给位前辈求助,不胜感激。

rmsd.png (100.93 KB, 下载次数 Times of downloads: 574)

Rmsd值

Rmsd值

rg.png (121.55 KB, 下载次数 Times of downloads: 532)

Rg值

Rg值

rmsf.png (74.58 KB, 下载次数 Times of downloads: 523)

Rmsf

Rmsf

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

2#
发表于 Post on 2019-6-21 00:28:36 | 只看该作者 Only view this author
Rg和RMSF分析没太大必要,说明不了什么

结合自由能分析就基本够了。也可以用VMD做一张图,把蛋白质部分消除平动和转动,然后把小分子在模拟过程中多帧(比如取20帧)叠加显示,从图上可以一目了然了解小分子在蛋白质内位置的波动程度。例如下图


你对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

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2019-6-23 11:00:07 | 只看该作者 Only view this author
sobereva 发表于 2019-6-21 00:28
Rg和RMSF分析没太大必要,说明不了什么

结合自由能分析就基本够了。也可以用VMD做一张图,把蛋白质部分 ...

好的  谢谢sob老师第一次发帖不够规范 抱歉。  我这里的RMSD都是复合物整体。另外想问下您利用VMD取多帧用的是xtc轨迹文件么 因为我的轨迹文件都10个g左右  不知道这样用VMD会不会很吃力

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2019-6-23 11:11:23 | 只看该作者 Only view this author
sobereva 发表于 2019-6-21 00:28
Rg和RMSF分析没太大必要,说明不了什么

结合自由能分析就基本够了。也可以用VMD做一张图,把蛋白质部分 ...

另外想问下老师VMD该如何做到叠加多帧显示呢

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

5#
发表于 Post on 2019-6-24 02:32:09 | 只看该作者 Only view this author
lunan815 发表于 2019-6-23 11:00
好的  谢谢sob老师第一次发帖不够规范 抱歉。  我这里的RMSD都是复合物整体。另外想问下您利用VMD取多帧 ...


你可以在VMD载入轨迹的界面里把stride选项设大,每隔一定帧数才载入一次,这样就不会造成载入过程内存不够

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
lunan815 + 4 谢谢

查看全部评分 View all ratings

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

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

6#
发表于 Post on 2019-6-24 02:34:12 | 只看该作者 Only view this author
lunan815 发表于 2019-6-23 11:11
另外想问下老师VMD该如何做到叠加多帧显示呢


评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
lunan815 + 4 谢谢

查看全部评分 View all ratings

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

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2019-6-24 10:29:23 | 只看该作者 Only view this author

谢谢sob老师!很有帮助!

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2019-6-24 19:40:28 | 只看该作者 Only view this author

sob老师打扰了   还有个问题想问下您 就是关于把蛋白质部分消除平动转动  用RMSDTrajectoryTool 这里面参数我需要怎么设置呢

rmsd_trjtool.png (37.82 KB, 下载次数 Times of downloads: 532)

rmsd_trjtool.png

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

9#
发表于 Post on 2019-6-25 04:22:16 | 只看该作者 Only view this author
lunan815 发表于 2019-6-24 19:40
sob老师打扰了   还有个问题想问下您 就是关于把蛋白质部分消除平动转动  用RMSDTrajectoryTool 这里面参 ...

点一下右上角的align按钮就行了
但往往是对于蛋白质骨架来消平动转动,此时把左边的backbone选上即可
北京科音自然科学研究中心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

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2019-6-25 12:18:25 | 只看该作者 Only view this author
sobereva 发表于 2019-6-25 04:22
点一下右上角的align按钮就行了
但往往是对于蛋白质骨架来消平动转动,此时把左边的backbone选上即可

谢谢老师  终于最后做出了这样的效果图 唯一就是发现把蛋白align后 当我用鼠标滑轮放大时发现会变淡直到消失如第二张图一样

1.JPG (19.25 KB, 下载次数 Times of downloads: 543)

1.JPG

2.JPG (10.77 KB, 下载次数 Times of downloads: 541)

2.JPG

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2019-6-25 12:30:55 | 只看该作者 Only view this author
lunan815 发表于 2019-6-25 12:18
谢谢老师  终于最后做出了这样的效果图 唯一就是发现把蛋白align后 当我用鼠标滑轮放大时发现会变淡直到 ...

已经解决了  谢谢Sob老师~

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2019-6-25 14:52:00 | 只看该作者 Only view this author
sobereva 发表于 2019-6-25 04:22
点一下右上角的align按钮就行了
但往往是对于蛋白质骨架来消平动转动,此时把左边的backbone选上即可

老师好  我根据您指导 做出小分子的叠加轨迹 十分感谢。然后想请教一下您 小分子在口袋内活动越平滑是代表结合越强烈么 还是说具体得看结合能

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

13#
发表于 Post on 2019-6-26 03:07:28 | 只看该作者 Only view this author
lunan815 发表于 2019-6-25 14:52
老师好  我根据您指导 做出小分子的叠加轨迹 十分感谢。然后想请教一下您 小分子在口袋内活动越平滑是代 ...

结合能体现绝对结合强度
分子位置波动范围越小体现结合越紧密
北京科音自然科学研究中心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

9

帖子

0

威望

66

eV
积分
75

Level 2 能力者

14#
 楼主 Author| 发表于 Post on 2019-6-26 13:32:13 | 只看该作者 Only view this author
sobereva 发表于 2019-6-26 03:07
结合能体现绝对结合强度
分子位置波动范围越小体现结合越紧密

学到了  谢谢老师!

2

帖子

0

威望

7

eV
积分
9

Level 1 能力者

15#
发表于 Post on 2022-8-23 10:21:09 | 只看该作者 Only view this author
您好,想请问下您做结合自由能分析的时候是用的什么做的分析呢

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 21:42 , Processed in 0.197964 second(s), 30 queries , Gzip On.

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