24小时热门版块排行榜    

查看: 4975  |  回复: 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的回帖

fanc232

新虫 (小有名气)

引用回帖:
8楼: Originally posted by limoujie1987 at 2014-09-12 09:34:30
应该是这样的,我后面做的对接没有考虑到这点,我是把pdbpt文件用记事本打开然后加的电荷。如果跟钙离子连接的原子加了的话,那么紧接着的其他原子电荷也应该要跟着变化,牵一发而动全身。看来应该不是金属离子电荷 ...

你手里要是有schrodinger2012软件,受体准备模块里面好像有处理金属离子的选项。大概意思是切断金属离子和受体之间的连接键并且分别加上合适的电荷。你可以试着用那个模块准备你的受体,然后试试看。
我很喜欢计算机辅助药物设计
9楼2014-09-12 10:32:59
已阅   回复此楼   关注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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 3/150 2026-08-27 00:09 by Ie9AyIAvGbvs
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +3 a089 2026-08-26 3/150 2026-08-26 23:21 by young90
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 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 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +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
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +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
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见