计算化学公社

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

[Amber] MM-PBSA中的熵到底能反应什么问题?

[复制链接 Copy URL]

149

帖子

0

威望

1087

eV
积分
1236

Level 4 (黑子)

请问一个比较迷惑的问题:MM-PBSA计算得到的结果包含了熵,这个熵主要是构象熵。不过我不太明白,计算的构象熵到底是什么概念?具体的讲,我研究蛋白-小分子配体作用,我计算了2个小分子与同一个蛋白的构象熵(TΔS),一个-21 kcal/mol,一个-19kcal/mol,他们的差异到底反应了什么问题?

294

帖子

0

威望

3533

eV
积分
3827

Level 5 (御坂)

2#
发表于 Post on 2015-4-9 18:21:15 | 只看该作者 Only view this author
看看这个帖子,不知道对你可有帮助http://bioms.org/forum.php?mod=viewthread&tid=1046

1149

帖子

6

威望

6635

eV
积分
7904

Level 6 (一方通行)

3#
发表于 Post on 2015-4-14 12:48:17 | 只看该作者 Only view this author
反应不了啥问题,熵的计算本来就不是很准确,加上力场的原因,2kcal/mol完全有可能是误差

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

4#
发表于 Post on 2015-4-18 05:55:48 | 只看该作者 Only view this author
模拟体系的熵可以分为构象熵和振动熵。
构象熵是柔性分子由于构象可变对熵的贡献。而像mmpbsa里面做的都是振动熵,而非构象熵。振动熵是特定构象下分子内振动对熵的贡献。受体配体结合前后它们的构象熵是会有改变的,由于结合后运动行为被限制住,可及构象减少,故构象熵会变小,从对结合自由能的贡献上看,构象熵起到的是去稳定化效应。做mmpbsa严格来说也应该把构象熵变化考虑进去,不过相对于mmpbsa已经引入的近似,以及这个过程构象熵变很小,可以忽略了。
北京科音自然科学研究中心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

149

帖子

0

威望

1087

eV
积分
1236

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2015-4-20 15:32:19 | 只看该作者 Only view this author
sobereva 发表于 2015-4-18 05:55
模拟体系的熵可以分为构象熵和振动熵。
构象熵是柔性分子由于构象可变对熵的贡献。而像mmpbsa里面做的都是 ...

我仔细看了一下结果,确实如您所说,里面有三项,Translational,Rotational ,Vibrational。不同配体-蛋白间的主要差距在于第三项也就是震动熵,其他两项差别很小。对于这个振动熵,假如两个配体与同一个蛋白结合的震动熵有差异,又能反应出来什么问题呢?

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

6#
发表于 Post on 2015-4-20 16:44:18 | 只看该作者 Only view this author
如果相差很小,不必在意,计算误差都能导致那种程度。
如果确实相差大,可能一定程度可以体现形成复合物后振动受阻差异。
北京科音自然科学研究中心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

149

帖子

0

威望

1087

eV
积分
1236

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2015-4-21 09:03:09 | 只看该作者 Only view this author
sobereva 发表于 2015-4-20 16:44
如果相差很小,不必在意,计算误差都能导致那种程度。
如果确实相差大,可能一定程度可以体现形成复合物后 ...

嗯。有道理。现在好多的结果都差别不大,尤其是小分子-蛋白作用。可能总自由能也就相差4,5 kcal/mol,所以假如熵相差2kcal/mol的时候,尽管我个人也感觉可能是误差,不过这个差值已经对最终的结果有一定影响了。假如要同时对比多个体系,这个熵值感觉可能会造成最终自由能排序的改变了。这也算是这个方法的缺陷吧。

112

帖子

0

威望

381

eV
积分
493

Level 3 能力者

8#
发表于 Post on 2020-7-29 13:55:40 | 只看该作者 Only view this author
ChaosChiao 发表于 2015-4-21 09:03
嗯。有道理。现在好多的结果都差别不大,尤其是小分子-蛋白作用。可能总自由能也就相差4,5 kcal/mol,所 ...

请问构象熵如何计算的呢?用什么程序计算?

7

帖子

0

威望

275

eV
积分
282

Level 3 能力者

9#
发表于 Post on 2021-9-23 22:15:34 | 只看该作者 Only view this author
hyr95 发表于 2020-7-29 13:55
请问构象熵如何计算的呢?用什么程序计算?

amber中的mmpbsa.py

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 02:29 , Processed in 0.196517 second(s), 20 queries , Gzip On.

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