AMBER:处理不规则残基中遇到的问题

amber对于蛋白中的残基可以进行处理,但是Amber的标准氨基酸残基库不包含含非标准残基的原子类型和电荷. 但它们是任何分子力学模拟都必需的.所以需要单独构建残基文件,推导非标准残基的原子电荷并判定其原子类型

至于如何判断蛋白中是否存在非标准残基,除开依据经验对pdb文件进行审查之外.,就是查看导入蛋白之后报的错误了.

Loading PDB file: ./5ds3_pre.pdb/public1/soft/amber/18-parallel-new/amber18/bin/teLeap: Warning!
Unknown residue: DSN   number: 20   type: Nonterminal
Creating new UNIT for residue: DSN sequence: 21/public1/soft/amber/18-parallel-new/amber18/bin/teLeap: Warning!
One sided connection. Residue () missing connect0 atom.
Created a new atom named: N within residue: .R<DSN 21>
Created a new atom named: CA within residue: .R<DSN 21>
Created a new atom named: C within residue: .R<DSN 21>
Created a new atom named: O within residue: .R<DSN 21>
Created a new atom named: CB within residue: .R<DSN 21>
Created a new atom named: OG within residue: .R<DSN 21>/public1/soft/amber/18-parallel-new/amber18/bin/teLeap: Warning!
One sided connection. Residue (default_name) missing connect1 atom.total atoms in file: 1934Leap added 1951 missing atoms according to residue templates:1951 H / lone pairsThe file contained 6 atoms not in residue templates
Loading Mol2 file: ./olaparib.mol2
Reading MOLECULE named 09L
Checking Unit./public1/soft/amber/18-parallel-new/amber18/bin/teLeap: Warning!
The unperturbed charge of the unit (4.000000) is not zero.
FATAL:  Atom .R<DSN 21>.A<N 1> does not have a type.
FATAL:  Atom .R<DSN 21>.A<CA 2> does not have a type.
FATAL:  Atom .R<DSN 21>.A<C 3> does not have a type.
FATAL:  Atom .R<DSN 21>.A<O 4> does not have a type.
FATAL:  Atom .R<DSN 21>.A<CB 5> does not have a type.
FATAL:  Atom .R<DSN 21>.A<OG 6> does not have a type.

如果发生上述错误,则表示蛋白中存在非标准残基.可以参照此链接:http://ambermd.org/tutorials/basic/tutorial5/index.htm进行残基的处理.
但现在从pdb数据库中下载的残基,使用如下命令处理后

antechamber -fi ccif -i DSN.cif -bk DSN -fo ac -o dsn.ac -c bcc -at amber <
报了如下错误:

Residue CRO has a type of LINKING. Quitting/

删除下载的cif文件中的

_chem_comp.type “D-PEPTIDE LINKING” <行
又报了另外的错误

processor.acdoctor mode is on: check and diagnosis problems in the input
file.antechamber: mmcif.c:70: rmmcif: Assertion `col2 &gt;= 0' failed.Aborted
(core dumped)*

需要定位到**$AMBERHOME/AmberTools/src/antechamber** 下的mmcif.c中的第68行那一段注释掉,也就是忽略掉蛋白分子中的LINKING部分
然后定位到当下文件夹中,重新安装antechamber

make intall

解决问题,恭喜发财!!!

AMBER:处理不规则残基中遇到的问题相关推荐

amber中生成小分子模板
第一步:生成小分子模板蛋白质中各氨基酸残基的力参数是预先存在的,但是很多模拟过程会涉及配体分子,这些有机小分子有很高的多样性,他们的力参数和静电信息不可能预存在库文件中,需要根据需要自己计算生成模板 ...
insert 语句的选择列表包含的项多于插入列表中的项_如何定义和使用Python列表（Lists）
Python中最简单的数据集合是一个列表(list).列表是方括号内用逗号分隔的任何数据项列表.通常,就像使用变量一样,使用=符号为Python列表分配名称. 如果列表中包含数字,则不要在其周围使用引 ...
六面体体积求解(规则不规则)
规则和不规则六面体求解方法1:立方体映射法(不规则六面体与立方体之间的映射) 这种方法的主要思想是把一个不规则六面体映射到一个立方体当中,通过计算不规则六面体与立方体之间的体积比例(分解为三维空间中 ...
向量几何在游戏编程中的使用
<1>简单的2-D追踪 -Twinsen编写 -本人水平有限,疏忽错误在所难免,还请各位数学高手.编程高手不吝赐教 -我的Email-address: popyy@netease.com ...
向量几何在游戏编程中的使用3
<3>2-D边界碰撞检测 -Twinsen编写 -本人水平有限,疏忽错误在所难免,还请各位数学高手.编程高手不吝赐教 -我的Email-address: popyy@netease.com ...
GraphBind: 通过层次图神经网络学习蛋白质局部结构嵌入规则以用于识别核酸结合残基
<GraphBind: protein structural context embedded rules learned by hierarchical graph neural networ ...
EXCEL 批量生成sheet表+生成超链接目录+某列按多条件去重罗列+提取单元格中的字母、数字
文章目录前言一.EXCEL 公式实现多个条件值匹配二.EXCEL 工作薄下SHEET表太多时你可以这样做 1.生成多个相同表头的SHEET表 1.1.批量生成SHEET表 1.2.批量更改多个S ...
Linux中autoduck批量对接,利用AutoDOCK进行分子对接（一）
8种机械键盘轴体对比本人程序员,要买一个写代码的键盘,请问红轴和茶轴怎么选? 分子对接(Moleculer-docking)理论分子对接就是两个或多个分子之间通过几何匹配和能量匹配相互识别找到最佳 ...
python求不规则图形面积_使用蒙特卡洛方法求解不规则图形的面积
使用蒙特卡洛方法求解不规则图形的面积周绪达 [摘要] [摘要]蒙特卡洛方法是一种以计算机为工具,通过抽样统计作为手段的一种解决问题的方法,在许多领域都有广泛应用.论文基于蒙特卡洛方法, 利用 ...

AMBER:处理不规则残基中遇到的问题

AMBER:处理不规则残基中遇到的问题相关推荐

最新文章

热门文章