计算化学公社

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

[GROMACS] 丝氨酸的左旋SER和右旋DSN两种手性分子的拓扑问题

[复制链接 Copy URL]

107

帖子

0

威望

467

eV
积分
574

Level 4 (黑子)

老师您好,我想分别创建丝氨酸的左旋SER和右旋DSN两种手性分子的拓扑文件。直接用pdb2gmx创建时,发现charmm36力场里本身就有对左旋SER的定义可以直接创建。
但是创建右旋时报错了Reside ‘DSN’ not found in residue topology database.但是我把pdb里的DSN字符替换成SER时发现可以创建,但是得到的结构文件里原子的排列顺序发生了变化,变成了和SER相同的顺序,请问老师我可以通过这种方式获得DSN的拓扑文件吗?谢谢老师

1625486632(1).png (91.67 KB, 下载次数 Times of downloads: 16)

merged.rtp里对SER的定义

merged.rtp里对SER的定义

1625487037(1).png (123.9 KB, 下载次数 Times of downloads: 14)

SER.pdb

SER.pdb

1625487037(1).png (123.9 KB, 下载次数 Times of downloads: 19)

DSN.pdb

DSN.pdb

1625487204(1).png (160.55 KB, 下载次数 Times of downloads: 20)

DSN.gro

DSN.gro

561

帖子

0

威望

3410

eV
积分
3971

Level 5 (御坂)

2#
发表于 Post on 2021-7-5 22:24:40 | 只看该作者 Only view this author
按理说左旋右旋只和坐标有关。你看看坐标的手性对不对就好了

513

帖子

0

威望

4334

eV
积分
4847

Level 6 (一方通行)

truffle

3#
发表于 Post on 2021-7-6 00:29:03 | 只看该作者 Only view this author
根据GROMACS二面角的定义,左旋和右旋的拓扑应该是一样的,感觉可以直接用,就是这个拓扑文件和你的坐标文件的原子顺序能对上就行了
No problem is insoluble in all conceivable circumstances.

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

4#
发表于 Post on 2021-7-6 08:51:06 | 只看该作者 Only view this author
右旋相当于非标准残基,自然在标准的rtp里面没有
你这种做法得到的拓扑文件能不能行实际跑起来看看DSN部分的结构就知道
北京科音自然科学研究中心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

107

帖子

0

威望

467

eV
积分
574

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-7-22 18:03:27 | 只看该作者 Only view this author
sobereva 发表于 2021-7-6 08:51
右旋相当于非标准残基,自然在标准的rtp里面没有
你这种做法得到的拓扑文件能不能行实际跑起来看看DSN部分 ...

谢谢老师,我用pdb2gmx对这两个pdb文件创建得到的结构确实是有差异的,是手性对称的,模拟跑了也没有出错,就是右旋和左旋的拓扑是一样的,用的都是标准的rtp里的,请问老师我可以这样操作吗?谢谢

107

帖子

0

威望

467

eV
积分
574

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2021-7-22 18:04:48 | 只看该作者 Only view this author
naoki 发表于 2021-7-6 00:29
根据GROMACS二面角的定义,左旋和右旋的拓扑应该是一样的,感觉可以直接用,就是这个拓扑文件和你的坐标文 ...

生成的结构文件里原子的排列顺序已经自动变成了拓扑文件里的顺序,运行也没有报错,请问可以这样创建吗

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

7#
发表于 Post on 2021-7-23 03:32:28 | 只看该作者 Only view this author
961106 发表于 2021-7-22 18:03
谢谢老师,我用pdb2gmx对这两个pdb文件创建得到的结构确实是有差异的,是手性对称的,模拟跑了也没有出错 ...

实际跑出来的轨迹看起来没异常就可以
北京科音自然科学研究中心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

107

帖子

0

威望

467

eV
积分
574

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-7-23 14:29:25 | 只看该作者 Only view this author
sobereva 发表于 2021-7-23 03:32
实际跑出来的轨迹看起来没异常就可以

再次感谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 06:04 , Processed in 0.171174 second(s), 23 queries , Gzip On.

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