计算化学公社

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

[GROMACS] 酶和四糖底物结合能分析,四糖底物拓扑结构生成求助

[复制链接 Copy URL]

26

帖子

0

威望

175

eV
积分
201

Level 3 能力者

各位老师们你们好,最近想做一个多糖裂解酶和四糖底物之间的结合能分析,利用pymol从NCBI复合物晶体结构中直接获取了四糖底物(GMGM)的pdb,利用Autodock vina对接,选择了一个对接位置和氢键数目比较理想的构象,这时候建议用什么软件怎么生成四糖底物的拓扑结构文件。

有想过利用charmm-gui直接生成复合物的gromacs输入文件,也试过单独生成配体的输入文件,但是均在生成PDB过程中出现Mismatch in ligand atom order,想请教各位前辈们有什么建议和看法呢?

202405191926552016..png (85.81 KB, 下载次数 Times of downloads: 9)

202405191926552016..png

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

9#
发表于 Post on 2024-5-23 04:29:21 | 只看该作者 Only view this author
艾小野 发表于 2024-5-22 09:58
谢谢社长老师的回复,后续查到资料复合物在跑NVT和npt平衡中,会对配体加限制,正式MD再取消限制。我当时 ...

当前没明显问题
北京科音自然科学研究中心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

26

帖子

0

威望

175

eV
积分
201

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2024-5-22 09:58:35 | 只看该作者 Only view this author
sobereva 发表于 2024-5-22 05:43
GAFF虽然没有glycam力场描述糖理想,但也没明显问题

谢谢社长老师的回复,后续查到资料复合物在跑NVT和npt平衡中,会对配体加限制,正式MD再取消限制。我当时没注意到这个,就没加,续跑45ns(加上预跑的5ns一共跑了50ns),[img]然后利用gmx energy看了50ns模拟中的温度和压力变化,温度稳定在303K,但是压力却是浮动在0~600bar之间,想请教社长老师这个压力浮动的情况正常吗?

之间看过社长老师发的“如何评价模拟平衡”帖子提到生物大分子体系要看rmsd比较准,也分析了protein_lig、protein、ligand的rmsd,前两者都稳定了,但是配体波动较大,[/img]尤其是后20ns过程中,看了vmd发现主是首糖和末糖单元的糖环发生明显的位置波动,想请教社长老师这与npt没加配体限制有关嘛?

我感觉好像关系不大,一是时间都跑这么久了,二是蛋白建模本身质量偏低,加上酶和底物对接整体效果一般,结合文献的活性位点选了个糖结合比较稳定的对接位置,这个位置首糖和末糖单元不可避免和周围残基作用比较一般,所以出现这个配体结合波动大的结果。

最后想请教社长:这种情况考虑继续跑还是重新换一个对接的配体,配体的波动感觉会好一点呢?

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

7#
发表于 Post on 2024-5-22 05:43:24 | 只看该作者 Only view this author
艾小野 发表于 2024-5-21 16:47
谢谢社长老师的回复!!我的配体就是四糖的寡糖,我觉得它分子量不算大,所以直接图省事方便,用了Sobtop ...

GAFF虽然没有glycam力场描述糖理想,但也没明显问题
北京科音自然科学研究中心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

26

帖子

0

威望

175

eV
积分
201

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2024-5-21 16:47:06 | 只看该作者 Only view this author
sobereva 发表于 2024-5-20 02:40
配体和多糖的拓扑文件应当分别产生
配体用sobtop (http://sobereva.com/soft/Sobtop)产生GAFF的拓扑文件, ...

谢谢社长老师的回复!!我的配体就是四糖的寡糖,我觉得它分子量不算大,所以直接图省事方便,用了Sobtop(真的超级好用,学习成本几乎为零,新手大喜)生成底物的itp gro top文件,然后蛋白用了常用组合力场Amber14sb,没对配体位置限制直接先跑了5ns,VMD下看了轨迹配体没有跑出口袋,这里想请教社长老师:Sobtop可以用来生成寡糖的拓扑文件嘛?然后还想问是不是配体模拟的最开始只要没有脱离的话就说明底物结合比较稳定,可以不用加位置限制预跑?

26

帖子

0

威望

175

eV
积分
201

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-5-20 11:20:43 | 只看该作者 Only view this author
sobereva 发表于 2024-5-20 02:38
有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员 ...

抱歉社长,下次一定注意

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

4#
发表于 Post on 2024-5-20 02:40:09 | 只看该作者 Only view this author
配体和多糖的拓扑文件应当分别产生
配体用sobtop (http://sobereva.com/soft/Sobtop)产生GAFF的拓扑文件,多糖部分可以用AmberTools里的leap产生,然后acpype转成gromacs的拓扑文件
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

3#
发表于 Post on 2024-5-20 02:38:21 | 只看该作者 Only view this author
艾小野 发表于 2024-5-19 19:29
这是Autodock对接前后的配体pdb文件,然后试过未对接的配体可以利用charmm-gui成功得到输入文件

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。
北京科音自然科学研究中心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

26

帖子

0

威望

175

eV
积分
201

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2024-5-19 19:29:54 | 只看该作者 Only view this author
这是Autodock对接前后的配体pdb文件,然后试过未对接的配体可以利用charmm-gui成功得到输入文件

gmgm.pdb

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

subsrate(docking).pdb

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

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

GMT+8, 2026-2-22 21:45 , Processed in 0.172976 second(s), 24 queries , Gzip On.

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