24小时热门版块排行榜    

查看: 3194  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

as23p

木虫 (小有名气)

[求助] IndexError: list index out of range,因为蛋白质过大? 已有1人参与

刚刚想用 autodock4.2.6 算算ceruloplasmin,pdb文件请见附件,然后生成pdbqt时显示:
Traceback (most recent call last):
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/AutoDockTools/Utilities24/prepare_receptor4.py", line 179, in <module>
    dict=dictionary)   
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/AutoDockTools/MoleculePreparation.py", line 560, in __init__
    version=version, delete_single_nonstd_residues=delete_single_nonstd_residues)
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/AutoDockTools/MoleculePreparation.py", line 126, in __init__
    self.repairMol(mol, self.repair_type_list)
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/AutoDockTools/MoleculePreparation.py", line 176, in repairMol
    self.newHs = self.addHydrogens(mol)
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/AutoDockTools/MoleculePreparation.py", line 189, in addHydrogens
    HB.addHydrogens(mol)
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/MolKit/hydrogenBuilder.py", line 81, in addHydrogens
    hat = AddHydrogens().addHydrogens(mol.allAtoms, method=self.method)
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/PyBabel/addh.py", line 92, in addHydrogens
    Hatoms = self.place_hydrogens1(atoms, num_H_to_add)
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/PyBabel/addh.py", line 149, in place_hydrogens1
    Hat = Hat + self.add_methyl_hydrogen(a, SP3_O_H_DIST)
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/MGLToolsPckgs/PyBabel/addh.py", line 303, in add_methyl_hydrogen
    c3 = atom2.bonds[1].atom1.coords
  File "/home/zhangss/suotc/x86_64Linux2/mgltools_x86_64Linux2_1.5.6/lib/python2.5/UserList.py", line 28, in __getitem__
    def __getitem__(self, i): return self.data
IndexError: list index out of range
是什么原因呢?
这个蛋白质很大,有一千多residue,是因为太大了吗?我把999以后的基团全删了就没这个错误了,但不想这么干,请问有没有别的办法?
因为我估计是autodock无法读取四位数的index,不知是不是这个原因,以及能不能改?
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : cerulo.pdb
  • 2015-07-24 15:22:27, 652.73 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

as23p

木虫 (小有名气)

★ ★
月只蓝: 金币+2, 感谢分享问题解决方案! 2015-07-28 10:51:59
月只蓝: 回帖置顶 2015-07-28 10:52:02
找到原因了,看起来是分子内结晶水和O2造成的问题。。。
8楼2015-07-28 10:07:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

蓝洋飘雪

木虫 (小有名气)

您好,我想问您一下,用Autodock 进行分子对接,如何选定要对接的区域?恳请您能帮帮忙!非常感谢
人生如茶,不会苦一辈子,但会苦一阵子!
2楼2015-07-24 16:46:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

as23p

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 蓝洋飘雪 at 2015-07-24 16:46:43
您好,我想问您一下,用Autodock 进行分子对接,如何选定要对接的区域?恳请您能帮帮忙!非常感谢

起始位置基本上是用眼看的,然后把盒子尽量做大
3楼2015-07-24 17:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pymol

捐助贵宾 (正式写手)


【答案】应助回帖

感谢参与,应助指数 +1
换一个版本试试,用1.5.6
4楼2015-07-25 11:03:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +8 Kittylucky 2026-08-27 8/400 2026-08-27 18:33 by raiy123
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 17:25 by xxniao123
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见