24小时热门版块排行榜    

查看: 2450  |  回复: 8

可米kemi

银虫 (小有名气)


[交流] Autodock对接结果可信度讨论

因为最近有在使用autodock,所以有几个问题想请教和大家交流请教以下。
    我是用Autodock做环糊精和药物的自动对接,本来是想先用Autodock找出药物和环糊精最稳定的结合位置和构象,然后把该构象拿到Gaussian进一步优化。
    一开始,我把格子大小设为60*60*60,可是结果出来的重复性不是很好,最高的那一簇能量总是变化,当然每簇之间相差倒也不大。
    后来,我干脆把格子大小设为40*40*40,然后又提高对接次数,结果出来后最高的那一簇能量到是符合了,而且也是能量最低点。可是,我还不满意的是,两次对接结果中,一簇高为81,另一簇高为53,这个偏差就有点大了吧。
    我想请教的是:
     1。其实我把格子设为40*40*40只是把格子大小为60*60*60的对接结果(能量最低点)拿来然后调整格子大小至40*40*40,观察并没有基团伸到格子外,就拿来用了。但我总觉得这样做好象不是很科学。请问,有没有什么方法可以观察这个对接设置是否可行?
     2。我两次对接结果能量最低点确实是最高簇,可是簇的高度相差有点大,一个81,一个53。这样的话,结果是否可信呢?如果不可信,又该怎么办呢?
     3。在对接结果中,总会提到RMSD,也就是均方根偏差,就是结果构象和目标构象的的偏差统计。我的理解就是每簇中的所有对接构像并不都是一样的,只是刚好RMSD在这个范围内,所以才放到同一簇。所以才会有“结果中选簇最高中的能量最低点这个说法”?我的理解就是即使有能量更低的簇,可是分子在那个构像附近出现的概率较低。不知道这个看法对不对。
回复此楼

» 收录本帖的淘帖专辑推荐

分子模拟

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leimiao_hit

木虫之王 (文学泰斗)



可米kemi(金币+1):谢谢参与
autodock
2楼2012-01-14 15:02:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hndxdjh

银虫 (正式写手)


★ ★ ★ ★ ★ ★
可米kemi(金币+1):谢谢参与
chaizhm(金币+5): 谢谢~ 2012-01-15 14:17:38
可米kemi(金币+20): 哦,那如果想要做二次对接或是三次对接,又该怎么做呢?谢谢。 2012-02-08 22:55:02
1.关于盒子大小的设置,用Autodock作对接你首先得充分了解你的活性部位,然后把它们全都包围着就可以了。或大或小应该问题不大。因为本来这个对接就是个参考,计算过程中的参数好多都是近似的。并且算是半柔性的。

2.如果你想更可靠,你可以作二次对接,甚至三次对接,不过花费的时间多。

3.RMSD是可以自己设置的,自己设置的越小分簇可能越多,如果最低能量和最高簇是同一簇,那么结果就会更可靠。

玩儿了近三年了。自己理解的吧。具体看你看文献自己理解了。
5楼2012-01-14 16:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hndxdjh

银虫 (正式写手)



小木虫(金币+0.5):给个红包,谢谢回帖
御剑江湖(金币+4): 谢谢 2012-02-17 22:03:25
在linux下运用命令把dlg文件转化为pdb,选择最低簇和最低能量的那一步的pdb文件复制出来,新建一个pdb文件。然后加氢加电荷做第二次对接。
具体命令的参考文献:J Med Chem. 2005, 48: 2308-2318;J Med Chem. 2006, 49: 2417-2430
7楼2012-02-12 13:45:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

方蓝雨

新虫 (小有名气)


9楼2012-10-05 23:29:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2012-01-14 15:32   回复  
可米kemi(金币+1):谢谢参与
2012-01-14 16:28   回复  
可米kemi(金币+1):谢谢参与
tjcdp6楼
2012-02-05 13:44   回复  
可米kemi(金币+1):谢谢参与
longytu8楼
2012-04-02 17:28   回复  
chaizhm: 编辑内容 2012-04-02 18:48
1 [ Last edited by chaizhm on 2012-4-2 at 18:48 ]
相关版块跳转 我要订阅楼主 可米kemi 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 5/250 2026-08-04 04:44 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:41 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见