24小时热门版块排行榜    

查看: 3223  |  回复: 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的回帖
回帖置顶 ( 共有1个 )

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的回帖
普通回帖

蓝洋飘雪

木虫 (小有名气)

您好,我想问您一下,用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的回帖

as23p

木虫 (小有名气)

引用回帖:
4楼: Originally posted by pymol at 2015-07-25 11:03:47
换一个版本试试,用1.5.6

就是用的1.5.6
5楼2015-07-25 11:31:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pymol

捐助贵宾 (正式写手)


【答案】应助回帖

★ ★
月只蓝: 金币+2, 感谢指导! 2015-07-28 10:51:53
引用回帖:
5楼: Originally posted by as23p at 2015-07-25 11:31:33
就是用的1.5.6...

没有问题的,是不是你的路径中有中文?
IndexError: list index out of range,因为蛋白质过大?
1.png

6楼2015-07-25 19:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

as23p

木虫 (小有名气)

引用回帖:
6楼: Originally posted by pymol at 2015-07-25 19:22:37
没有问题的,是不是你的路径中有中文?

1.png
...

没有中文。。。问题是出在加氢的时候,莫非你的可以顺利生成 pdbqt 文件?
7楼2015-07-26 13:23:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小zxx

铁虫 (小有名气)

引用回帖:
8楼: Originally posted by as23p at 2015-07-28 10:07:18
找到原因了,看起来是分子内结晶水和O2造成的问题。。。

楼主 请问最后这个问题您是怎么解决的 谢谢
9楼2016-12-04 16:11:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 as23p 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:30 by cITaGg3p5edr
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-18 02:28 by cITaGg3p5edr
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-17 23:28 by cuPZTDXS3VOd
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 23:28 by cuPZTDXS3VOd
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 LwdutQ8HoqWP 2026-09-13 8/400 2026-09-17 22:53 by cuPZTDXS3VOd
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 7/350 2026-09-17 16:40 by 北京莱茵编辑
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-17 13:29 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 yaCh4X3wd045 2026-09-14 4/200 2026-09-17 13:17 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 12:29 by ObVzqQQh1rrE
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 04:52 by k06BKNblrGRH
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:26 by k06BKNblrGRH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 04:04 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 17:12 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
信息提示
请填处理意见