计算化学公社

 找回密码 Forget password
 注册 Register

使用MCPB.py创建金属蛋白拓扑遇到keyerror等问题

查看数: 1761 | 评论数: 3 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-12-23 19:54

正文摘要:

本帖最后由 zhu2023 于 2023-12-23 20:00 编辑 最近在利用MCPB.py创建金属蛋白拓扑文件,遇到了个问题。 准备:目的蛋白是napa.pdb(体积太大无法上传,未去配体见下图),配体是2个MGD、1个金属MO和一个S原子 ...

回复 Reply

17x7 发表于 Post on 2024-5-25 16:21:44
我已经解决了,怎么说呢,上面说的上面ion id纯属就是胡扯,直接抄原子序号就好了,真正会报这个错误是因为格式转换的错误,就是mcpb是amber下面的程序,如果你是是再外面加氢,H的原子名字就是无法识别,所以就会报key error的错误,可能还会给你标出来是哪一个残基的,但是你优先去看可以形成大小模型的那个氨基酸残基,你直接蛋白dry,之后用amber自带的reduce加氢,之后再pdb4amber转化,就ok,如果你的小配体是正常bcc做的,原子名和H是没有问题的,我习惯用高斯做了生出生成mol2文件,之后形成frcmod文件,所以就是需要用match_atomname来原子名对齐,翻了好多都没有解决,还是我自己来写
17x7 发表于 Post on 2024-5-25 15:17:23
解决了嘛,我也遇到相同的问题,就不是ion_id问题
MilesYYh 发表于 Post on 2024-5-17 18:20:28
是运行MCPB.py -i XXX.in -s 1时出的问题吗?如果是可能是你的XXX.in 文件中的“ion_ids 数字”这个数字写错了,这里该填写的是你金属离子所在的那一行,使得你的文件不是对应的,建议先按照这个流程(http://bbs.keinsci.com/forum.php ... &highlight=MCPB)以例子中运行一遍。

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

GMT+8, 2026-2-19 06:06 , Processed in 0.193097 second(s), 31 queries , Gzip On.

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