计算化学公社

 找回密码 Forget password
 注册 Register

AMBER非标准氨基酸时出现could not find bond parameter

查看数: 2208 | 评论数: 15 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-10-20 21:47

正文摘要:

在体系准备时处理含非标准氨基酸的蛋白,生成非标准氨基酸的参数后进行tleap的时候check出现问题:

回复 Reply

zxczxc 发表于 Post on 2025-5-15 15:56:46
你好,我的好像和你的类似,缺少结晶水的角度参数,请问该如何解决呢?
/opt/amber22/bin/teLeap: Error!
Could not find angle parameter for atom types: HW - Y4 - HW
        for atom H1 at position 8.762411, -11.683037, -11.870854,
            atom O at position 7.805411, -11.683037, -11.870854,
        and atom H2 at position 7.565411, -10.756037, -11.870854.

/opt/amber22/bin/teLeap: Error!
Could not find angle parameter for atom types: HW - Y5 - HW
        for atom H1 at position 10.797411, -9.612037, -12.714854,
            atom O at position 9.840411, -9.612037, -12.714854,
        and atom H2 at position 9.600411, -8.685037, -12.714854.

/opt/amber22/bin/teLeap: Error!
Could not find angle parameter for atom types: HW - Y6 - HW
        for atom H1 at position 8.124411, -8.884037, -11.312854,
            atom O at position 7.167411, -8.884037, -11.312854,
        and atom H2 at position 6.927411, -7.957037, -11.312854.
Building proper torsion parameters.
Building improper torsion parameters.
total 2760 improper torsions applied
Building H-Bond parameters.
Incorporating Non-Bonded adjustments.

/opt/amber22/bin/teLeap: Warning!
Parameter file was not saved.
ovolcano 发表于 Post on 2024-7-22 13:36:06
kongmenghua 发表于 2024-5-11 17:56
你好,这个问题解决了嘛,我也出现了这样的问题

非标准氨基酸的两端不要只用氢原子补全,用ACE和NME封端后再进行参数生成,就不会有参数缺少的问题了。
kongmenghua 发表于 Post on 2024-5-11 17:56:27
你好,这个问题解决了嘛,我也出现了这样的问题
Frozen-Penguin 发表于 Post on 2023-10-23 18:16:56
ovolcano 发表于 2023-10-23 16:39
我的非标准氨基酸残基主链应该是NT-CT-C(NT是与上一个残基的C相连的原子,CT是αC原子,C是与下一个残基 ...

直接写 source leaprc.protein.ff19SB 就可以了,19SB的参数里面有C-N*。用其他的AMBER力场也可以,根据实际的需求决定。
ovolcano 发表于 Post on 2023-10-23 16:42:59
八月的雨季 发表于 2023-10-22 15:25
就是你漏掉了一些键、角参数,补上就行了

是需要自己手动补上吗?
ovolcano 发表于 Post on 2023-10-23 16:39:01
Frozen-Penguin 发表于 2023-10-21 20:54
C-NT这个键应该是肽键里面的吧,所以不会出现在残基内,如果你是按照antechamber的流程做的,过程中应该 ...

我的非标准氨基酸残基主链应该是NT-CT-C(NT是与上一个残基的C相连的原子,CT是αC原子,C是与下一个残基的N相连的原子),C和NT没有相连,要说C-NT键存在那应该是非标准残基与相邻残基形成的键,那我生成的非标准残基的参数确实是没有这个键的参数的。关于您说的这个“肽键的参数之间用amber力场的”应该怎么做?
八月的雨季 发表于 Post on 2023-10-22 15:25:34
就是你漏掉了一些键、角参数,补上就行了
rpestana94 发表于 Post on 2023-10-22 00:35:38
ovolcano 发表于 2023-10-21 03:42
The parameter information of the angle that was not found also should not exist in my protein. Bec ...

Sometimes the error is due to the name of the residue but usually when you load the prep and frcmod file that you create this error is solve
Frozen-Penguin 发表于 Post on 2023-10-21 20:54:02
ovolcano 发表于 2023-10-21 09:52
这个是我生成的prepin文件,从文件看C和NT应该没有成键,所以我觉得生成的frcmod文件应该没有缺少信息, ...

C-NT这个键应该是肽键里面的吧,所以不会出现在残基内,如果你是按照antechamber的流程做的,过程中应该写过一个主链信息文件,这个文件里写了这个键的信息。如果非标准残基的参数用了gaff2,肽键的参数之间用amber力场的就可以了。
ovolcano 发表于 Post on 2023-10-21 16:42:40
rpestana94 发表于 2023-10-21 10:54
Which is the non-standard residue and what about the angle parameters which are the other giving you ...

The parameter information of the angle that was not found also should not exist in my protein. Because there is no C-NT bond in the protein, there should naturally be no information about the angles formed by C, NT  and other atoms. And after checking, the BOND section in the obtained frcmod file includes all the bonds of the non-standard residue. I don't see any problems with the preprin and frcmod files. I want to know why it still looks for information about this C-NT bond which should not already exist?
rpestana94 发表于 Post on 2023-10-21 10:54:52
Which is the non-standard residue and what about the angle parameters which are the other giving you errors
rpestana94 发表于 Post on 2023-10-21 00:37:13
What you show is just the parameters for bond and C-T bond dosen't exist, you should create the parameters for this residue and the problem will be solve, this can be done using gaussian to calculate the charges and obtain the prepi and frcmod file

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

GMT+8, 2025-8-16 22:22 , Processed in 0.186115 second(s), 31 queries , Gzip On.

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