24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1719  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +4 努力奋斗112 2026-04-07 4/200 2026-04-07 10:30 by unicornone
[考研] 071000生物学,一志愿深圳大学296分,求调剂 +5 TIckLw 2026-04-06 5/250 2026-04-07 07:56 by dongzh2009
[考研] 308求调剂 +11 墨墨漠 2026-04-06 11/550 2026-04-06 22:56 by 875465
[考研] 0703求调剂383分 +9 W55j 2026-04-03 9/450 2026-04-06 06:50 by houyaoxu
[考研] 化学357分,考研调剂 +11 .Starry. 2026-04-04 12/600 2026-04-06 06:28 by houyaoxu
[考研] 316求调剂 +5 yyx想调剂 2026-04-05 5/250 2026-04-05 22:22 by 咔咔咔咔9
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 271分求调剂学校 +12 zph158488! 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 专硕310求调剂 +5 捞捞我…. 2026-04-04 6/300 2026-04-04 23:33 by barlinike
[考研] 环境285分,过六级,求调剂 +10 xhr12 2026-04-02 10/500 2026-04-04 21:53 by bn53987
[考研] 085400电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-03 6/300 2026-04-04 21:50 by hemengdong
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 10/500 2026-04-03 11:08 by Kittylucky
[考研] 重庆大学材料与化工085600,初试370+,求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 285求调剂 +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
[考研] 农学考研求调剂 +3 dkdkxm 2026-04-01 3/150 2026-04-02 16:04 by wangjagri
[考研] 求调剂推荐 +3 南山南@ 2026-04-01 3/150 2026-04-02 12:09 by xiaoranmu
信息提示
请填处理意见