计算化学公社

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

[GROMACS] 求助关于gromacs中的溶剂可及面积的计算

[复制链接 Copy URL]

23

帖子

0

威望

175

eV
积分
198

Level 3 能力者

各位大神好,我输入的命令是
gmx sasa -f md_nopbc.xtc -s md.tpr -o area.xvg -odg dgsolv.xvg -or resarea.xvg -oa atomarea.xvg -surface  -tv volume.xvg
这样的,为什么我的输出文件中只有两行数据,一行是时间,一行是总的可及面积,我想得到亲水的表面积和亲油的表面积,那我应该怎么计算呢,求大佬指导!感激不尽!


6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

2#
发表于 Post on 2018-8-24 18:59:44 | 只看该作者 Only view this author
参考北京科音分子动力学与GROMACS培训班(http://www.keinsci.com/KGMX)幻灯片:
北京科音自然科学研究中心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

13

帖子

0

威望

111

eV
积分
124

Level 2 能力者

3#
发表于 Post on 2020-7-30 12:18:50 | 只看该作者 Only view this author
请教sob老师,为什么可以按照原子电荷区分疏水和亲水呢

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

4#
发表于 Post on 2020-8-9 05:49:24 | 只看该作者 Only view this author
muxchong 发表于 2020-7-30 12:18
请教sob老师,为什么可以按照原子电荷区分疏水和亲水呢

显然原子电荷较大的原子与水静电相互作用强,可以视为亲水区
北京科音自然科学研究中心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

31

帖子

0

威望

1418

eV
积分
1449

Level 4 (黑子)

5#
发表于 Post on 2021-1-10 22:15:13 | 只看该作者 Only view this author

sob老师,我在用gmx sasa 算的时候直接选择溶质分子,计算的area只是溶质的表面积是吗?如果我的溶质在两相界面处,我想计算溶质和水的接触面积以及溶质和甲苯的接触面积,我应该怎么设置呢?非常感谢sob老师

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

6#
发表于 Post on 2021-1-11 02:44:19 | 只看该作者 Only view this author
weiyi8061225 发表于 2021-1-10 22:15
sob老师,我在用gmx sasa 算的时候直接选择溶质分子,计算的area只是溶质的表面积是吗?如果我的溶质在两 ...

gmx sasa算的是SASA,根据SAS的定义,没法说溶质与比如甲苯的接触面积。本来算SASA的计算就需要定义一个溶剂球半径(程序默认用1.4埃,对应水的)
如果随便瞎算一下的数的话,可以尝试:比如算溶质和甲苯的接触,用溶质+甲苯算一次作为1,单独用溶质算一次作为2,单独用甲苯算一次作为3,然后求(2+3-1)/2。不过这种做法意义不清楚
还一个可以尝试的是:算一次溶质+甲苯,并且用-restrict限制只计算溶质那部分原子的SAS,结果作为1。然后溶质算一次,作为2。用2减去1,就得到由于溶质与甲苯接触,其SASA减少的量。这种做法或许相对合理一些,不过我没试过。
以上SASA在计算,原理上应当设成很小的探针半径。

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

68

帖子

0

威望

407

eV
积分
475

Level 3 能力者

7#
发表于 Post on 2023-2-2 23:37:48 | 只看该作者 Only view this author
sobereva 发表于 2020-8-9 05:49
显然原子电荷较大的原子与水静电相互作用强,可以视为亲水区

请问老师,用-0.2 to 0.2的依据是什么?是有这样的定义,还是得根据自己的体系选择合理的区间?
感谢老师指点!

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

8#
发表于 Post on 2023-2-2 23:53:56 | 只看该作者 Only view this author
Jus 发表于 2023-2-2 23:37
请问老师,用-0.2 to 0.2的依据是什么?是有这样的定义,还是得根据自己的体系选择合理的区间?
感谢老 ...

基本理论化学直觉
对原子电荷有足够多的了解、算过各种体系的原子电荷,自然知道这是个合适的区间
北京科音自然科学研究中心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

68

帖子

0

威望

407

eV
积分
475

Level 3 能力者

9#
发表于 Post on 2023-2-3 13:39:59 | 只看该作者 Only view this author
sobereva 发表于 2023-2-2 23:53
基本理论化学直觉
对原子电荷有足够多的了解、算过各种体系的原子电荷,自然知道这是个合适的区间

好的,谢谢老师

18

帖子

0

威望

101

eV
积分
119

Level 2 能力者

10#
发表于 Post on 2024-6-21 11:14:03 | 只看该作者 Only view this author

请问老师,在计算SASA时遇到这个警告,可以怎么解决,需要替换输入文件吗?
Command line:
  gmx sasa -f md.trr -s md.tpr -n sasa.ndx -o sasa.xvg -or resarea.xvg -tv timevol.xvg -b 140000 -e 200000 -probe 0.14

WARNING: Masses and atomic (Van der Waals) radii will be guessed
         based on residue and atom names, since they could not be
         definitively assigned from the information in your input
         files. These guessed numbers might deviate from the mass
         and radius of the atom type. Please check the output
         files if necessary.
诚盼老师的回复!谢谢老师

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

11#
发表于 Post on 2024-6-22 02:17:50 | 只看该作者 Only view this author
肥力丝 发表于 2024-6-21 11:14
请问老师,在计算SASA时遇到这个警告,可以怎么解决,需要替换输入文件吗?
Command line:
  gmx sasa  ...

不用管,通常对结果没什么影响
北京科音自然科学研究中心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

18

帖子

0

威望

101

eV
积分
119

Level 2 能力者

12#
发表于 Post on 2024-6-22 20:26:33 | 只看该作者 Only view this author
sobereva 发表于 2024-6-22 02:17
不用管,通常对结果没什么影响

好的,十分感谢老师!

19

帖子

0

威望

77

eV
积分
96

Level 2 能力者

13#
发表于 Post on 2025-6-7 09:57:49 | 只看该作者 Only view this author
请问SASA的生物学意义是什么,主要是为了评价什么?如果是在蛋白-小分子相互作用的背景下,有什么指导意义吗?
想把计算化学公社装到我的脑子里

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

14#
发表于 Post on 2025-6-8 02:06:35 | 只看该作者 Only view this author
lele123 发表于 2025-6-7 09:57
请问SASA的生物学意义是什么,主要是为了评价什么?如果是在蛋白-小分子相互作用的背景下,有什么指导意义 ...

没有直接的生物学的意义。结合具体问题才可能有实际意义

蛋白-小分子相互作用能经常用MMPBSA或MMGBSA算结合自由能,其中SA部分是指基于SASA计算溶剂的非极性部分对结合能产生的贡献
北京科音自然科学研究中心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

19

帖子

0

威望

77

eV
积分
96

Level 2 能力者

15#
发表于 Post on 2025-6-10 09:06:00 | 只看该作者 Only view this author
sobereva 发表于 2025-6-8 02:06
没有直接的生物学的意义。结合具体问题才可能有实际意义

蛋白-小分子相互作用能经常用MMPBSA或MMGBSA ...

谢谢sob老师解答
想把计算化学公社装到我的脑子里

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 20:26 , Processed in 1.472677 second(s), 23 queries , Gzip On.

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