“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

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

[GROMACS] 怎样处理对共价对接后的复合物在有机溶剂做MD模拟中蛋白质电荷不为0的问题

[复制链接 Copy URL]

5

帖子

0

威望

33

eV
积分
38

Level 2 能力者

本人做的是含有非标准氨基酸的但蛋白质的分子动力学模拟(方法参考的是Gromacs非标准氨基酸残基的蛋白MD建模:sobtop产生rtp文件 - 分子模拟 (Molecular Modeling) - 计算化学公社),所选的溶剂是甲苯溶剂,因为是该蛋白质是酶与酯共价结合形成的酶酰复合物,在网上查找资料所选用共价对接的PH值为7.4,该PH值下复合物的c端脯氨酸是去质子化的,导致在做分子动力学模拟时复合物的电荷值不为零。请问可以直接在maestro中补齐这个氢做局部能量最小化吗?还是在社长计算resp电荷的教程中使用resp2脚本在计算几何优化之后使用几何优化的fchk文件得到pdb文件。还是通过其他方式防止脯氨酸去质子化使复合物的静电荷为0。如果该帖子表达有问题、不清晰会第一时间修改,请各位批评指正。

6万

帖子

99

威望

6万

eV
积分
126401

管理员

公社社长

2#
发表于 Post on 前天 10:57 | 只看该作者 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

5

帖子

0

威望

33

eV
积分
38

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 前天 12:15 | 只看该作者 Only view this author
sobereva 发表于 2026-4-21 10:57
甲苯里不可能去质子化。各氨基酸应当按照电中性质子化状态(体现在结构和拓扑信息层面)做模拟

谢谢老师的回答。是的,应该是做共价对接的时候ph值7.4造成的我查了一下pro的pka小于7.4,请问应该怎样设置电中性质子化状态呢?我现在用maestro补了pro的H,是使用pdb2gro的时候使用 -ter -inter逐个设置他的质子化状态吗?我在做对接的时候ph值设置的为7.4。

6万

帖子

99

威望

6万

eV
积分
126401

管理员

公社社长

4#
发表于 Post on 前天 15:11 | 只看该作者 Only view this author
greatti 发表于 2026-4-21 12:15
谢谢老师的回答。是的,应该是做共价对接的时候ph值7.4造成的我查了一下pro的pka小于7.4,请问应该怎样设 ...

没有叫pdb2gro的东西
水环境下判断质子化态的那些情况对于甲苯环境完全不适用,说pKa那些都没意义
北京科音自然科学研究中心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

5

帖子

0

威望

33

eV
积分
38

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 前天 15:35 | 只看该作者 Only view this author
本帖最后由 greatti 于 2026-4-21 15:48 编辑
sobereva 发表于 2026-4-21 15:11
没有叫pdb2gro的东西
水环境下判断质子化态的那些情况对于甲苯环境完全不适用,说pKa那些都没意义

是的sob老师,我之前的想法有问题,我还有一些问题想请教sob老师
1、请问电中性质子化状态的指的是蛋白质中每一个氨基酸都不会有质子化状态对吗,
2、那这个时候还可以用amber99sb-ildn力场吗?
3、在做分子动力学模拟的时候和普通的水环境的分子动力学模拟除了判断是否加离子还有什么变化吗?
4、请问sob老师我是不是应该使用-ter -inter手动设置氨基酸的质子化状态。

6万

帖子

99

威望

6万

eV
积分
126401

管理员

公社社长

6#
发表于 Post on 前天 18:00 | 只看该作者 Only view this author
greatti 发表于 2026-4-21 15:35
是的sob老师,我之前的想法有问题,我还有一些问题想请教sob老师
1、请问电中性质子化状态的指的是蛋白 ...

1 不解离自己的质子、不额外结合质子
2 AMBER力场是给水溶液环境下模拟设计的。GROMOS 54B7适合模拟这种情况
3 语义不明
4 是
北京科音自然科学研究中心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

5

帖子

0

威望

33

eV
积分
38

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 6 hour ago | 只看该作者 Only view this author
社长好,有三个问题想请教一下您
1、我看了gromacs力场非标准氨基酸的补充文件和amber的不一样,这很正常,我没有找到教程不敢自己写。还可以继续用amber吗?
2、在使用 -inter参数的时候arg精氨酸并没有中性质子化状态,我可以再做一个非标准氨基酸的力场补充文件吗?
3、我是用resp2脚本算出的非标准氨基酸的电荷总量不为零约为0.041,这是正常的吗,我使用的约束方式为约束两端封端分子为0,总体为0.简而言之就是http://bbs.keinsci.com/forum.php ... digest%26digest%3D1这篇教程的约束方式。如果不正常是哪里出了问题应该如何改正呢?
    我本来是截取的非标准氨基酸两端氨基酸残基加氢的封端方式,现在尝试用ace和nme封端,会是这个原因吗?

5

帖子

0

威望

33

eV
积分
38

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 1 hour ago | 只看该作者 Only view this author
sobereva 发表于 2026-4-21 18:00
1 不解离自己的质子、不额外结合质子
2 AMBER力场是给水溶液环境下模拟设计的。GROMOS 54B7适合模拟这种 ...

请问sob老师,我做的是一个在甲苯中的合成酶所以这个复合物是一个中间体,后续对中间体进行采样,然后再做对接,再做模拟,如果使用电中性的质子化状态会对酶活性产生影响吧,怎样最大限度的模拟真实情况呢?全部做电中性处理是可以的吗?

本版积分规则 Credits rule

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

GMT+8, 2026-4-23 23:09 , Processed in 0.172524 second(s), 20 queries , Gzip On.

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