24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1732  |  回复: 10
当前主题已经存档。

xingyun5182

铜虫 (初入文坛)

[交流] [求助]Autodock对接结果Reference RMSD过大怎么办?

Autodock对接后的结果显示ref rmsd值很大,都在10左右。
我没有sybyl,因此ligand的pdb文件是用下载的复合物pdb文件中的配体部分用写字板复制粘贴做成的。Ligand和MacroMolecule的相对位置与原复合物pdb中相对位置是相同的。
Autodock的FAQs里说在DPF中可以specify the "rmsref" command,但我在用ADT创建DPF文件的时候并没有找到这个选项,请问这个命令应该如何设定呢?

这是FAQs中的原文:
I get very high Reference RMSD values in my DLG (docking log file); what went wrong?
The "Reference RMSD" values that are printed in the "RMSD TABLE" in the DLG are computed from the coordinates of

either the input ligand (PDBQ or PDBQT) file specified by the "move" command in the DPF, if you did not include the "rmsref" command in your DPF;
or the ligand (PDBQ or PDBQT) file you specified in the "rmsref" command in the DPF.
If you do not specify the "rmsref" command, and the ligand input coordinates happen to be translated far from the receptor, you will appear to get high Reference RMSD values. _Don't Panic!_ This is normal, and you don't need to worry about these high values.

You may want to check that the input ligand coordinates are far from the crystallographically observed binding position in active site, by reading in the input ligand and the receptor in ADT.

We usually use the "rmsref" DPF-command to specify the x-ray crystallographic coordinates of a known binding mode taken from a complex PDB structure. This can be a useful way of checking if your redocking is successful, if the Reference RMSD values are less than 2-3 Å from the crystal structure position of the ligand.

谢谢各位指点。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
楼主,我认为你不宜为了减小ref rmsd而去更改这几个参数,首先因为ref rmsd代表不了什么,其次,每次遗传都用相同的初始构象,你想想这会发生什么?这就是削弱了遗传算法威力最大的地方:全局寻优能力。你这样很可能找不到全局最优构象,除非你已经很肯定这个配体就应该在这个地方,只是来优化一下他的侧链位置。所以,我觉得,还是random的好罢!
个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
8楼2008-07-17 09:40:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长


suntao1982(金币+1,VIP+0):xiexie
我认为去掉random是不合适的,应该以控制盒子的位置和大小来控制这种情况
个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
10楼2008-07-17 18:57:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

xxffliu

铜虫 (小有名气)


suntao1982(金币+1,VIP+0):谢谢你!
rmsref参数可以让你指定晶体结构的配体文件作为计算rmsd的参考。具体在Docking->Docking Run Parameters里可以设置。
2楼2008-07-14 18:48:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingyun5182

铜虫 (初入文坛)

能再具体的指一下设定的位置么?
Docking -> Docking Parameters -> ? -> ?
我用的是ADT1.4.6。
谢谢。
3楼2008-07-14 20:12:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weishenme

金虫 (小有名气)


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
复合物中配体去掉了吗?
4楼2008-07-15 06:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingyun5182

铜虫 (初入文坛)

复合物中的配体已经去掉了,我是直接在写字板中删除的。用Viewerliite和ADT打开后均未再发现配体。

[ Last edited by xingyun5182 on 2008-7-15 at 09:03 ]
5楼2008-07-15 09:02:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长


suntao1982(金币+1,VIP+0):辛苦了,欢迎常来!
ref rmsd值太大无关紧要,这主要取决于你以哪个配体结构作为reference,这个是可以在ADT中重新设定和,计算的,也可以重新指定reference ligand,在analysize->conformations->play...中就能。
个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
6楼2008-07-15 09:38:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingyun5182

铜虫 (初入文坛)

嗯,终于找到设置的地方了,ADT -> Docking -> Ligand -> Ligand Parameters,把那几个random点掉就可以固定Ligand的位置了,对接后的ref rmsd就很小了。
谢谢楼上各位的帮忙指点。
7楼2008-07-16 15:12:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingyun5182

铜虫 (初入文坛)

我这次对接就是把从RCSB中下载下来的复合物中的受体和配体拆开后再进行对接的。当我设置random后,发现配体分子到处乱跑,和原来的位置相距甚远,也就是ref rmsd值特别大,这个时候是不是就应该将random点掉了?
9楼2008-07-17 10:39:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xingyun5182 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华南理工大学331分材料求调剂 +5 天下ww 2026-04-09 5/250 2026-04-10 08:36 by 5268321
[考研] 267求调剂 +5 再忙也要吃饭啊 2026-04-09 5/250 2026-04-09 18:47 by stone_128
[考研] 电子信息270求调剂 +11 terminal469 2026-04-07 11/550 2026-04-09 18:28 by hy861222
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +15 努力奋斗112 2026-04-06 17/850 2026-04-08 22:27 by 凯凯要变帅
[考研] 298求调剂 +4 manman511 2026-04-05 4/200 2026-04-08 16:50 by tjzhao
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +20 努力奋斗112 2026-04-07 21/1050 2026-04-08 14:54 by screening
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 080500求调剂 +12 黄宇博 2026-04-06 12/600 2026-04-07 12:41 by upczlm1989
[考研] 一志愿太原理工大学计算机技术专硕348,求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 一志愿上海海洋大学083200食品学硕,求调剂,接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 085400电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-03 6/300 2026-04-04 21:50 by hemengdong
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 一志愿双非085502,267分,过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
信息提示
请填处理意见