24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 4665  |  回复: 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的回帖

limoujie1987

铜虫 (初入文坛)

引用回帖:
7楼: Originally posted by fanc232 at 2014-09-11 18:21:37
给钙离子加了正电荷,那么附近的受体上的原子就应该加上相应的负电荷吧。...

应该是这样的,我后面做的对接没有考虑到这点,我是把pdbpt文件用记事本打开然后加的电荷。如果跟钙离子连接的原子加了的话,那么紧接着的其他原子电荷也应该要跟着变化,牵一发而动全身。看来应该不是金属离子电荷的问题。
8楼2014-09-12 09:34:30
已阅   回复此楼   关注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,310分求调剂 +3 李潇可 2026-03-26 3/150 2026-03-27 03:39 by wxiongid
[考研] 食品工程专硕求调剂 +3 小张zxy张 2026-03-26 3/150 2026-03-27 01:14 by dgnhs
[考研] 324求调剂 +3 hanamiko 2026-03-26 3/150 2026-03-26 21:00 by sanrepian
[考研] 求调剂 +5 芦lty 2026-03-25 6/300 2026-03-26 20:49 by 不吃魚的貓
[考研] 325求调剂 +3 Aoyijiang 2026-03-23 3/150 2026-03-26 20:46 by 不吃魚的貓
[考研] 调剂 +4 柚柚yoyo 2026-03-26 4/200 2026-03-26 20:43 by fmesaito
[考研] 303求调剂 +7 元夕元 2026-03-20 8/400 2026-03-26 20:38 by 不吃魚的貓
[考研] 281求调剂 +3 亚克西good 2026-03-26 5/250 2026-03-26 19:48 by 不吃魚的貓
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +5 星沉uu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[考研] 086000生物与医药292求调剂 +6 小小陈小小 2026-03-22 9/450 2026-03-26 15:58 by dick_runner
[考研] 资源与环境 调剂申请(333分) +9 holy J 2026-03-21 9/450 2026-03-26 15:47 by 161765490
[考研] 总分293求调剂 +6 加一一九 2026-03-25 8/400 2026-03-26 13:30 by yujianx
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 296求调剂 +4 汪!?! 2026-03-25 7/350 2026-03-25 16:41 by 汪!?!
[考研] 生物学学硕求调剂 +7 小羊睡着了? 2026-03-23 10/500 2026-03-25 02:24 by 清风拂扬。 m
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 300求调剂,材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
信息提示
请填处理意见