计算化学公社

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

[GROMACS] 已解决:em后肽键断裂

[复制链接 Copy URL]

58

帖子

0

威望

297

eV
积分
355

Level 3 能力者

本帖最后由 wyfgg 于 2022-11-30 16:57 编辑

各位老师好,我目前在进行一个人工合成多肽的MD模拟,多肽中包含两个非标准氨基酸。我运行em.mdp后,发现em.gro其中一个非标准氨基酸DLE6与LEU7的肽键断裂,请问这个可能是哪方面出现的问题呢?和我建立的残基rtp有关吗?各位老师有没有什么好的解决办法?
谢谢
更新:我单独取出断裂肽键前后的残基进行独立的em,发现他们并不会断裂,这是不是说明,这个多肽本身不稳定(实际上应该是很稳定的。。。)

202211301156254847..png (10 KB, 下载次数 Times of downloads: 9)

202211301156254847..png

202211301200071277..png (104.52 KB, 下载次数 Times of downloads: 7)

202211301200071277..png

202211301414349001..png (45.01 KB, 下载次数 Times of downloads: 9)

202211301414349001..png

202211301418215401..png (109.78 KB, 下载次数 Times of downloads: 10)

202211301418215401..png

310

帖子

0

威望

1503

eV
积分
1813

Level 5 (御坂)

2#
发表于 Post on 2022-11-30 14:58:47 | 只看该作者 Only view this author
可能模拟的系统中没有这个肽键,可以在top文件中检查,在[ atoms ]部分找到成键的N和C原子的编号,在[ bonds ]部分查找是否有这个键。
出现上述情况可能是由于rtp文件中没有肽键的信息,肽键的写法可以参考力场目录下aminoacids.rtp的写法。

58

帖子

0

威望

297

eV
积分
355

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-11-30 16:56:42 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2022-11-30 14:58
可能模拟的系统中没有这个肽键,可以在top文件中检查,在[ atoms ]部分找到成键的N和C原子的编号,在[ bond ...

谢谢您!已经解决

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 10:23 , Processed in 0.257407 second(s), 23 queries , Gzip On.

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