24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1737  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 282,电气工程专业,求调剂,不挑专业 +7 jggshjkkm 2026-04-10 7/350 2026-04-10 13:43 by 木子君1218
[基金申请] 有爆料,一个青年教师卖房得400万,然后换了一个四青帽子 +9 babu2015 2026-04-08 9/450 2026-04-10 11:43 by 苏东坡二世
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 一志愿211电子信息347求调剂 +3 554916 2026-04-03 4/200 2026-04-09 22:24 by 554916
[考研] 291分调剂 +3 上岸小莹加油 2026-04-09 4/200 2026-04-09 17:28 by 1753564080
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 086000调剂 +4 十七sa 2026-04-07 4/200 2026-04-09 09:27 by 壹往無前
[考研] 化学0703-一志愿211-338分求调剂 +10 vants 2026-04-05 11/550 2026-04-08 16:02 by screening
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 调剂求助(生物与医药) +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 304求调剂(085602,过四级,一志愿985) +25 化工人999 2026-04-04 26/1300 2026-04-07 22:06 by hemengdong
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 296求调剂 +3 汪!?! 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
信息提示
请填处理意见