计算化学公社

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

[GROMACS] 如何将聚合物分子共价连接到蛋白酶上进行模拟?

[复制链接 Copy URL]

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

如题所示,已经将聚合物在真空中模拟蜷缩,如何将聚合物(PEG,共1500个原子)共价连接到蛋白酶的N端(pdb文件中残基编号为1的ASP)进行模拟?两个附件分别是蛋白酶的pdb和聚合物蜷缩后的gro文件

peg_trans.gro

107.11 KB, 下载次数 Times of downloads: 2

oph.pdb

333.54 KB, 下载次数 Times of downloads: 0

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2024-4-9 11:54:18 | 只看该作者 Only view this author
sobereva 发表于 2024-4-9 11:15
如果考虑了周期性,把盒子显示出来弄清楚情况。并且观看优化轨迹

要么就只跑聚合物,先不管蛋白

好的老师,谢谢老师!

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

11#
发表于 Post on 2024-4-9 11:15:51 | 只看该作者 Only view this author
YuniJ 发表于 2024-4-9 11:12
谢谢老师,老师不好意思,我做完以上操作把peg+蛋白酶放到真空中能量最小化后聚合物远离了蛋白酶,如图所 ...

如果考虑了周期性,把盒子显示出来弄清楚情况。并且观看优化轨迹

要么就只跑聚合物,先不管蛋白
北京科音自然科学研究中心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

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2024-4-9 11:12:55 | 只看该作者 Only view this author
sobereva 发表于 2024-4-9 10:11
都可以。但如果你有明确的要结合的位点,先放到真空里蜷缩,有可能就缩到里面去了,不好再接上

谢谢老师,老师不好意思,我做完以上操作把peg+蛋白酶放到真空中能量最小化后聚合物远离了蛋白酶,如图所示,在能量最小化过程中我是限制了蛋白酶的位置。这种情况是为什么呢?请问该怎么解决了。谢谢老师

202404091112105154..png (19.5 KB, 下载次数 Times of downloads: 18)

能量最小化前

能量最小化前

202404091111412826..png (14.84 KB, 下载次数 Times of downloads: 21)

能量最小化后

能量最小化后

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

9#
发表于 Post on 2024-4-9 10:11:21 | 只看该作者 Only view this author
YuniJ 发表于 2024-4-9 08:06
好的谢谢老师,那请问老师我是先共价连接然后把聚合物+蛋白酶体系放到真空中模拟然后位置限制蛋白质进行 ...

都可以。但如果你有明确的要结合的位点,先放到真空里蜷缩,有可能就缩到里面去了,不好再接上
北京科音自然科学研究中心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

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2024-4-9 08:06:52 | 只看该作者 Only view this author
sobereva 发表于 2024-4-8 02:36
字段用的序号是全局序号,即和结构文件里的序号对应

簇模型仅仅是为了帮助你确定连接的地方用什么成键 ...

好的谢谢老师,那请问老师我是先共价连接然后把聚合物+蛋白酶体系放到真空中模拟然后位置限制蛋白质进行聚合物的蜷缩呢还是先单独真空蜷缩聚合物呢?

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

7#
发表于 Post on 2024-4-8 02:36:25 | 只看该作者 Only view this author
YuniJ 发表于 2024-4-7 09:48
好的谢谢老师,我按照您的提示成功生成了bond等信息,我想请问一下我需要更改这里面的原子序号吗?就是扣 ...

[intermolecular_interactions]字段用的序号是全局序号,即和结构文件里的序号对应

簇模型仅仅是为了帮助你确定连接的地方用什么成键相关的参数合适,簇模型的原子和你实际要计算的模型序号没有对应关系
北京科音自然科学研究中心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

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2024-4-7 09:48:23 | 只看该作者 Only view this author
sobereva 发表于 2024-4-7 05:13
显然要抠出一部分,边界进行饱和,构成一个像样的分子,让连接区域的结构处在这个模型体系的中间区域,而 ...

好的谢谢老师,我按照您的提示成功生成了bond等信息,我想请问一下我需要更改这里面的原子序号吗?就是扣出连接区域后添加了氢原子,我又把添加的氢原子增加到了原结构文件里。那么利用[intermolecular_interactions]字段时,我是需要把bonds等字段中的原子序号更改为原gro文件里面的原子序号吗?

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

5#
发表于 Post on 2024-4-7 05:13:26 | 只看该作者 Only view this author
YuniJ 发表于 2024-4-6 08:38
老师打扰您了,我手动合并结构文件后,此时peg末端的C原子和oph中1号残基ASP中的N原子(需要共价连接的两 ...

显然要抠出一部分,边界进行饱和,构成一个像样的分子,让连接区域的结构处在这个模型体系的中间区域,而不是光抠出来两个原子,这样根本什么都不是
北京科音自然科学研究中心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

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-4-6 21:45:20 | 只看该作者 Only view this author
本帖最后由 YuniJ 于 2024-4-6 23:26 编辑
sobereva 发表于 2024-4-3 18:59
建模方面,用VMD载入两个结构文件,挪其中一个到差不多的位置,保存,然后手动合并结构文件

拓扑文件方 ...

请问老师我可以从哪里找到关于此类操作的相关文献呢?


139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-4-6 08:38:40 | 只看该作者 Only view this author
本帖最后由 YuniJ 于 2024-4-6 23:25 编辑
sobereva 发表于 2024-4-3 18:59
建模方面,用VMD载入两个结构文件,挪其中一个到差不多的位置,保存,然后手动合并结构文件

拓扑文件方 ...

老师打扰您了,我手动合并结构文件后,此时peg末端的C原子和oph中1号残基ASP中的N原子(需要共价连接的两个原子)还没有结合(即游离状态)。我把N原子和C原子单独抠出来,用sobtop产生基于GAFF的参数,但是生成的拓扑文件没有bonds等相关信息,这种情况应该怎么解决呢?谢谢老师

peg_oph_bonds.mol2

246 Bytes, 下载次数 Times of downloads: 2

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

2#
发表于 Post on 2024-4-3 18:59:56 | 只看该作者 Only view this author
建模方面,用VMD载入两个结构文件,挪其中一个到差不多的位置,保存,然后手动合并结构文件

拓扑文件方面,pdb2gmx构建蛋白质拓扑文件,sobtop构建聚合物拓扑文件,然后利用[intermolecular_interactions]字段给两部分之间增加恰当的bond、angle、dihedral项(如果不清楚用什么参数合适,可以把连接区域抠一部分,用sobtop产生基于GAFF的参数,然后挪用进去)
北京科音自然科学研究中心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

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

GMT+8, 2026-2-22 23:49 , Processed in 0.183355 second(s), 24 queries , Gzip On.

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