24小时热门版块排行榜    

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

limoujie1987

铜虫 (初入文坛)

[求助] Autodock 对接结果已出,最后一步分析show interactions 时出错,求大神赐教 已有2人参与

对接信息:Windows 7系统,下载安装的autodock4.2.6最新版本,MGLTools 1.5.6,受体大分子下载自PDB数据库,小分子通过Chem3D生成,能量并最小化。大分子去水加全氢(查询资料有的地方说是加极性氢,手册里面的是全氢,我选择的全氢),Grid步骤选择大分子时初始化大分子正常,但提示分子中有些CA没有电荷,我直接点确定跳过了,上面两点是我觉得可能出问题的地方,后面对接步骤都正常,得到结果文件后进行 show interactions 分析,出现错误,求大神赐教。

分析结果时弹出的错误提示如下:

Unable to assign HYB type to atom N
Unable to assign HYB type to atom N
ERROR *********************************************
Traceback (most recent call last):
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\ViewerFramework\VF.py", line 898, in tryto
    result = command( *args, **kw )
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\AutoDockTools\autoanalyzeCommands.py", line 3087, in doit
    self.build()
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\AutoDockTools\autoanalyzeCommands.py", line 2912, in build
    self.intDescr = InteractionDescriptor(lig, macro, percentCutoff)
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\MolKit\interactionDescriptor.py", line 49, in __init__
    self.build(detect_pi=detect_pi)
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\MolKit\interactionDescriptor.py", line 60, in build
    self.buildCloseContactAtoms(percentCutoff)              #
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\MolKit\interactionDescriptor.py", line 69, in buildCloseContactAtoms
    self.macro_atoms, percentCutoff=percentCutoff)
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\MolKit\distanceSelector.py", line 106, in select
    bigC = Numeric.resize(c, (lenC, lenC, 3))
  File "C:\Program Files\MGLTools-1.5.6\lib\site-packages\numpy\core\fromnumeric.py", line 615, in resize
    a = concatenate( (a,)*n_copies)
MemoryError
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pymol

捐助贵宾 (正式写手)


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
limoujie1987: 金币+7 2014-09-11 14:57:30
打开dlg文件,打开大分子,然后play by energy,就可以show interaction!
金属离子的事情可以通过删除金属离子,或者增加金属离子参数来完成对接
4楼2014-09-10 20:23:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

limoujie1987

铜虫 (初入文坛)

加全氢还是加极性氢我已经验证过了,得到的结果都出现这个错误。至于钙离子没带电应该也不是原因。极有可能是我准备大分子有些步骤出错了。有没有人能回复一下自己成功案例是怎么准备受体大分子的步骤的?
2楼2014-09-10 13:34:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanc232

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
limoujie1987: 金币+8 2014-09-11 14:57:21
受体大分子,我平时是拿SYBYL的Biopolymer模块来做受体准备的。原因是里面的界面能够对受体进行各种检查,异常的氨基酸残基能够在界面中显示,方便查看它与对接位点的距离。氨基酸残基的质子化状态和互变异构,不会以报错的形式出现,需要你自己检查。
搜索报错,http://autodock.1369657.n2.nabbl ... m-Mg-td1493897.htmlhttp://mgl.scripps.edu/forum/viewtopic.php?f=11&t=2039都和你的挺像。看样子是加电荷的问题。Ca离子附近是不是有两个N原子啊,是不是那两个N原子没有正确地加电荷?
我很喜欢计算机辅助药物设计
3楼2014-09-10 16:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limoujie1987

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by pymol at 2014-09-10 20:23:19
打开dlg文件,打开大分子,然后play by energy,就可以show interaction!
金属离子的事情可以通过删除金属离子,或者增加金属离子参数来完成对接

忘了具体说明了,我使用的是Autodock Vina 进行的对接,所以木有生成dlg文件。苦恼啊。我觉得应该是大分子的问题。还有请问如何增加金属离子参数?我搜了一下,没找到相关方法。
5楼2014-09-11 15:30:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 (金币+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
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见