应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » DS/Sybyl/Autodock » IndexError: list index out of range，因为蛋白质过大？
91/1返回列表
查看: 3202  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢
1楼 2015-07-24 15:27:42
已阅   回复此楼   关注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的回帖

as23p

木虫 (小有名气)
★ ★
月只蓝: 金币+2, 感谢分享问题解决方案！ 2015-07-28 10:51:59
月只蓝: 回帖置顶 2015-07-28 10:52:02
找到原因了，看起来是分子内结晶水和O2造成的问题。。。
赞一下(1人) 回复此楼
8楼2015-07-28 10:07:18
已阅   回复此楼   关注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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +9 Kittylucky 2026-08-27 9/450 2026-08-28 01:08 by jurkat.1640
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +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
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
关闭
请填处理意见
关闭