请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 126|回复: 6

[Amber] 求助:leap自动键连如何消除

[复制链接]

56

帖子

0

威望

158

eV
积分
214

Level 3 能力者

发表于 2019-9-28 23:46:38 | 显示全部楼层 |阅读模式
为一个超级复杂的多肽搭建分子动力学立场多肽中有两组桥连形成内环,且都是非标准氨基酸,只能自己搭建力场,我把两组力场都搭建完成后运行tleap

它提示出现了许多非正常距离的键长,但是没有给我是哪两个原子之间的
/home/apps/chem/amber18/bin/teLeap: Warning!
There is a bond of 6.413129 angstroms between:
/home/apps/chem/amber18/bin/teLeap: Warning!
There is a bond of 9.451063 angstroms between:
/home/apps/chem/amber18/bin/teLeap: Warning!
There is a bond of 3.084539 angstroms between:
/home/apps/chem/amber18/bin/teLeap: Warning!
There is a bond of 3.432854 angstroms between:
/home/apps/chem/amber18/bin/teLeap: Warning!
There is a bond of 3.059956 angstroms between:
/home/apps/chem/amber18/bin/teLeap: Warning!
There is a bond of 11.831135 angstroms between:

(部分,更多的就不展示了)

其次出现了许多未知的键角力场参缺失,但是我所有的S都是做在自定义的力场AAA和BBB中,其中一组包含S的缺失显然不是由于力场错误导致的,而可能就是之前那些奇怪的键连导致的
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: N - C - N
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: N - CT - N
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: N - C - N
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: N - CT - N
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: N - C - N
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: N - CT - S
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: O2 - C - N
/home/apps/chem/amber18/bin/teLeap: Error!
Could not find angle parameter: O2 - C - N


求助:
1. tleap报错如何纠正,是否能继续下去
2. 这样的力场搭建方式是否可行(因为一般的共价结合搭建力场都是在一个氨基酸上做修饰,但这个桥连是在两个氨基酸之间连接,如果两个氨基酸一起取出来做力场就需要像“工”字形一样人为bond四个地方,这样就不得不打断原有的氨基酸序列形成间隔,这时如果leap自动把跳开的两个残基键连就会产生一些奇怪的键长)

leap.out

7.21 KB, 下载次数: 2

leap.in

526 Bytes, 下载次数: 3

BBB.mol2

3.64 KB, 下载次数: 2

BBB.frcmod

7.45 KB, 下载次数: 2

AAA.mol2

2.28 KB, 下载次数: 2

AAA.frcmod

3.54 KB, 下载次数: 2

RES.pdb

23.5 KB, 下载次数: 3

上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU

10

帖子

0

威望

52

eV
积分
62

Level 2 能力者

发表于 7 天前 | 显示全部楼层
原理我不太懂,操作的话,我把你的pdb文件放在viewer pro里面测量了精氨酸N-C-N的角度,然后把参数补充到frcmod文件中,再次运行leap.sh,得到的out文件少了3个错误。按照该操作过程,找到其他的参数,就不会有Error了。

10

帖子

0

威望

52

eV
积分
62

Level 2 能力者

发表于 7 天前 | 显示全部楼层
嘟比妹妹 发表于 2019-10-7 15:53
原理我不太懂,操作的话,我把你的pdb文件放在viewer pro里面测量了精氨酸N-C-N的角度,然后把参数补充到fr ...

如图
QQ图片20191007160210.png
QQ图片20191007160221.png

2万

帖子

25

威望

2万

eV
积分
46654

管理员

公社社长

发表于 6 天前 | 显示全部楼层
对这个问题我暂时没时间细看,可以把缺的参数在Google里搜,如果有就拿来用;或者看看GAFF里有没有能借用的;还不行就量化优化一下结构,扫一下键角势能曲线,把平衡键角和力常数补进去
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

56

帖子

0

威望

158

eV
积分
214

Level 3 能力者

 楼主| 发表于 5 天前 | 显示全部楼层

非常感谢您的回答,如果是补充键的信息我也知道,问题在于出现了不少不应该出现的键角和键长。
也就是说,你补的这个键角理论上不应该由我来补充,要么是AMBER没有识别出来要么就是这里错误成键,如果不加以甄别只是把缺失的信息补上是可以成功tleap但是其实拓扑关系是错误的,继续计算没有意义。
我比较在意的是为什么会出现这个错误。。。。。。
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU

77

帖子

0

威望

1053

eV
积分
1130

Level 4 (黑子)

发表于 5 天前 | 显示全部楼层
DoubeeTwT 发表于 2019-10-9 16:27
非常感谢您的回答,如果是补充键的信息我也知道,问题在于出现了不少不应该出现的键角和键长。
也就是说 ...

稍微看了一下,发现你的mol2文件中的bond有一些错误,如BBB中键2 - 4 对应N2 - C2长度是12A多,那些长键的信息,就是来自这里。你再检查检查你的mol2文件吧.

56

帖子

0

威望

158

eV
积分
214

Level 3 能力者

 楼主| 发表于 3 天前 | 显示全部楼层
beyond 发表于 2019-10-9 18:02
稍微看了一下,发现你的mol2文件中的bond有一些错误,如BBB中键2 - 4 对应N2 - C2长度是12A多,那些长键 ...

哦!太感谢了!!!!
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-10-14 22:22 , Processed in 0.169452 second(s), 26 queries .

快速回复 返回顶部 返回列表