24小时热门版块排行榜    

查看: 4655  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 305分求调剂(食品工程) +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
信息提示
请填处理意见