计算化学公社

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

[新手求助] 如何进行大体系下的量化计算?

[复制链接 Copy URL]

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

楼主想计算5000个原子的酶分子与若干有机化合物之间的结合,然后使用multiwfn研究在氨基、苯环等各主要位点上的(sobereva.com)若相互作用、结合能、结合常数。研究方法的选择上遇到了困惑,传统的Gromacs分子动力学可以满足计算量要求,但是其主要针对氢键有关的相互作用比较好,似乎很难描述色散作用和计算结合常数。如果用ORCA计算,更符合楼主习惯和倾向性,但是此情况需要描述的是酶分子和几十个小有机化合物之间的相互作用,似乎初始优化的平衡状态都难以进行。请问各位大神,有无更好的建议?或者已有的研究模式可供参考?

谢谢!

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

2#
发表于 Post on 2020-10-21 15:37:30 | 只看该作者 Only view this author
用QM/MM,活性位点附近100个原子左右用QM,其余用MM
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

339

帖子

0

威望

5053

eV
积分
5392

Level 6 (一方通行)

3#
发表于 Post on 2020-10-21 19:06:46 来自手机 | 只看该作者 Only view this author
只能用GFN-FF来搞了,据大佬宣称这个力场精度可以达到QM的精度。

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

4#
发表于 Post on 2020-10-21 19:43:19 | 只看该作者 Only view this author
本帖最后由 喵星大佬 于 2020-10-21 19:44 编辑

其实2000个原子的话GFN-xTB也干的动原文里就有优化蛋白结构的
GFN-FF也用GFN2-xTB来对比的,半斤八两

393

帖子

1

威望

5145

eV
积分
5558

Level 6 (一方通行)

5#
发表于 Post on 2020-10-21 20:45:59 | 只看该作者 Only view this author
GFN-FF貌似很吃内存,一个2K结构用GFN-FF跑,192GB内存一致报错,从来没有成功过,就此放弃。

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2020-10-21 21:28:08 | 只看该作者 Only view this author
wzkchem5 发表于 2020-10-21 15:37
用QM/MM,活性位点附近100个原子左右用QM,其余用MM

如果研究酶催化中心介导的化学反应,用QM/MM应该没有问题。但是如果一个蛋白质结构依赖于本身的氨基官能团非共价键结合化合物。显然QM/MM将无从下手,因为该结合位点并不表现出高度选择性,而且该氨基结构遍布蛋白质或者酶体系。

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

7#
发表于 Post on 2020-10-21 21:30:46 | 只看该作者 Only view this author
lao7 发表于 2020-10-21 21:28
如果研究酶催化中心介导的化学反应,用QM/MM应该没有问题。但是如果一个蛋白质结构依赖于本身的氨基官能 ...

那我觉得既然是算非共价作用,那可以看考这个博文

http://sobereva.com/214

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2020-10-21 21:39:54 | 只看该作者 Only view this author
喵星大佬 发表于 2020-10-21 19:43
其实2000个原子的话GFN-xTB也干的动原文里就有优化蛋白结构的
GFN-FF也用GFN2-xTB来对比的,半斤八两

理论上,量化计算需要一个大致合理的初始结构。一个大蛋白质分子如果耦合几十个小离子,计算体系的初始结构如何摆弄?是类似于Gromacs的随机放入吗?谢谢。

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

9#
发表于 Post on 2020-10-21 21:45:11 | 只看该作者 Only view this author
lao7 发表于 2020-10-21 21:39
理论上,量化计算需要一个大致合理的初始结构。一个大蛋白质分子如果耦合几十个小离子,计算体系的初始结 ...

你先用力场优化个差不多的出来啊

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2020-10-21 21:49:57 | 只看该作者 Only view this author
喵星大佬 发表于 2020-10-21 21:30
那我觉得既然是算非共价作用,那可以看考这个博文

http://sobereva.com/214

谢谢!受用。之前阅读过。总觉得初始结构不好弄。即一个大蛋白质分子和小化合物分子,非共价键的作用位点可能好多。到底要摆放多少小化合物分子合适?到底怎么给出初始结构?

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

11#
发表于 Post on 2020-10-21 21:55:14 | 只看该作者 Only view this author
lao7 发表于 2020-10-21 21:49
谢谢!受用。之前阅读过。总觉得初始结构不好弄。即一个大蛋白质分子和小化合物分子,非共价键的作用位点 ...

可以先跑MD看看情况,对于能用量化算的体系跑MD轨迹可以跑很久/很多条轨迹,选择合适的力场还是能提供很多信息的

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

12#
发表于 Post on 2020-10-21 23:13:26 | 只看该作者 Only view this author
lao7 发表于 2020-10-21 21:28
如果研究酶催化中心介导的化学反应,用QM/MM应该没有问题。但是如果一个蛋白质结构依赖于本身的氨基官能 ...

哦,我之前没理解,你说的是酶,我就自动理解成只和催化位点作用了
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

13#
发表于 Post on 2020-10-22 04:25:07 | 只看该作者 Only view this author
“似乎很难描述色散作用和计算结合常数”这是不对的。普通分子力场能合理描述色散作用,本身LJ势里就有专门的1/r^6项描述。

如果有很多位点需要考虑,量化计算时候分别挨个考虑,计算哪个位点时就把哪个位点挖成簇模型计算。
北京科音自然科学研究中心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

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

14#
 楼主 Author| 发表于 Post on 2020-10-22 09:20:49 | 只看该作者 Only view this author
wzkchem5 发表于 2020-10-21 23:13
哦,我之前没理解,你说的是酶,我就自动理解成只和催化位点作用了

谢谢您的回复!对于小分子化合物和酶蛋白的相互作用来说,可能非专一性的结合更具有普适性和代表性。

739

帖子

0

威望

1715

eV
积分
2454

Level 5 (御坂)

15#
 楼主 Author| 发表于 Post on 2020-10-22 09:28:49 | 只看该作者 Only view this author
sobereva 发表于 2020-10-22 04:25
“似乎很难描述色散作用和计算结合常数”这是不对的。普通分子力场能合理描述色散作用,本身LJ势里就有专门 ...

谢谢Sob.!请问关于第一个回答:根据您的经验,描述常规的CHON化合物(蛋白质体系)时, 什么力场对色散作用描述的更好些?

第二个回答:提到的簇模型,应该指的是以官能团(靶点)为中心半径截取的模型片段,请问这个截取的模型片段大小及尺寸有无约定俗成的半径范围?谢谢。

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 05:47 , Processed in 0.167053 second(s), 21 queries , Gzip On.

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