计算化学公社

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

[GROMACS] Lys和Asp侧链异肽键相关力场参数的添加

[复制链接 Copy URL]

6

帖子

0

威望

365

eV
积分
371

Level 3 能力者

本帖最后由 MELO 于 2024-6-1 13:45 编辑

各位老师好,最近课题组在进行融合蛋白的研究,需要对力场文件添加基于ASP和LYS参数修改后含异肽键的ASQ和LYT。
      我的问题是模拟的融合蛋白中包含异肽键,即Lys侧链中的NZ和Asp侧链中的CG进行连接,如图一所示。在修改amber14sb_parmbsc1.ff力场中aminoacids.hdbaminoacids.rtp,以及top文件夹中residuetypes.dat和specbond.dat后加氢成功,但无法正常生成肽键-NH-CO-,需要手动在topol.top中手动添加这一异肽键。

虽然加氢后的初始构象没问题,但是在经过能量最小化后的构象出现如图二所示的错误结构。     

这几个力场文件还需要调整哪些参数,还需要设置力场中哪些文件,还请各位老师能够指点一下!谢谢!





能量最小化后的构象.png (76.32 KB, 下载次数 Times of downloads: 5)

图二 能量最小化后异肽键的结构

图二 能量最小化后异肽键的结构

加氢后的构象.png (237.62 KB, 下载次数 Times of downloads: 5)

图一 正确加氢后的初始构象

图一 正确加氢后的初始构象

aminoacids.hdb

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

加氢

aminoacids.rtp

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

specbond.dat

479 Bytes, 下载次数 Times of downloads: 0

79

帖子

2

威望

743

eV
积分
862

Level 4 (黑子)

分子模拟晶戈

2#
发表于 Post on 2024-5-31 16:45:10 | 只看该作者 Only view this author
我当初做过一个类似的建模,是把两个MET的S用一个有机基团连接起来,实验组的同学把这个叫做“硫盐”。我当时的实现是准备好为成键的两边,用pdb2gmx生成拓扑文件,然后再手动生成连接部分的各种键连关系。这个方法有点繁琐,但应该是可以work的

6

帖子

0

威望

365

eV
积分
371

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-6-1 13:33:48 | 只看该作者 Only view this author
对抗路达摩 发表于 2024-5-31 16:45
我当初做过一个类似的建模,是把两个MET的S用一个有机基团连接起来,实验组的同学把这个叫做“硫盐”。我当 ...

我也是想从根源上解决这个问题,直接从力场文件中添加相关参数,直接通过pdb2gmx生成这一个键,但是目前没有解决这个问题。

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 21:46 , Processed in 0.725923 second(s), 24 queries , Gzip On.

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