24小时热门版块排行榜    

查看: 2639  |  回复: 6
当前主题已经存档。

patent

至尊木虫 (正式写手)

[交流] 如何计算autodock的rmsd值

这个rmsd是拿谁做参照啊
docking前的默认设置是我输入的ligand的pdb为reference,当然出来的值都特别大,有的超过30多。
我想拿结合在蛋白质上的底物做reference,请问怎么操作?
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alwens

铁杆木虫 (正式写手)

老木虫

★ ★
patent(金币+2):thank you
你输入的ligand如果就是蛋白底物的话,你只要用底物的坐标作为对接的初始坐标就ok了
autodock的rmsd只能是同一分子对接前后的坐标标准偏差。因为是逐原子比较,因此如果你的ligand不是底物的话(原子不对应)就比较麻烦。
一个解决的办法是:你要先把你的ligand移动到底物所在的位置,如果ligand比较相象,比如有公共子结构,可以选择这部分把ligand作为source位置,底物作为target位置作superimpose。然后用更新的ligand的位置作为对接的初始位置进行autodock。

[ Last edited by alwens on 2007-1-14 at 19:56 ]
万里夕阳锦背高 翻身犹恨东洋小 太公怎钓?
2楼2007-01-14 17:46:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjmed

木虫 (正式写手)

★ ★ ★
patent(金币+3):非常感谢!!!!!!!!!!!!!!!!!!!!!!
一般做虚拟筛选之前,都是先用晶体结构中的配体进行分子对接,看结果能否重现晶体结构中的结合模式
所以,如果你的对接结果和晶体结构中的结合模式比较接近,也就是对接产生的小分子构象和晶体结构中的参考分子构象比较接近,则认为是成功的,表现就是rmsd值比较小,一般在1.5A以内就可以接受
具体说来可以这样:
我一般在sybyl中操作
先提取晶体结构中的配体分子,另存为mol2分子
然后进行分子对接,autodock得到的结果都是pdb分子,没有关系,一般可能会有10个分子在一个pdb文件中。
在读入sybyl时,会提示读入10个中的拿一个。一个一个的读进去,然后就可以用sybyl中的match计算参考分子和对接结果的rmsd了
当然,这个不是最方便的,建议你使用openeye公司的vida或者AFITT,这个用来观察对接结果非常方便。不足的是不能给出具体的rmsd值
Timeandtidewaitsfornoman!
3楼2007-01-14 19:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

patent

至尊木虫 (正式写手)

引用回帖:
Originally posted by alwens at 2007-1-14 05:46 PM:
一个解决的办法是:你要先把你的ligand移动到底物所在的位置,如果ligand比较相象,比如有公共子结构,可以选择这部分把ligand作为source位置,底物作为target位置作superimpose。然后用更新的ligand的位置作为对接的初始位置进行autodock

这个我不是很明白,老兄,能否详细说明一下,用autodock具体怎样操作呢?
万分感谢!!!!!!!!!!!!!!!!
4楼2007-01-15 00:54:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjmed

木虫 (正式写手)


killl(金币+1):感谢参与应助,欢迎多来求助板块帮助别人。
其实alwens兄弟指的是如何确定受体的活性位点坐标
这个直接可以通过修改.gpf文件中的center grid后的三个坐标文件数值就可以了
可以把晶体结构中配体的坐标作为.gpf文件中的坐标值
Timeandtidewaitsfornoman!
5楼2007-01-15 09:40:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alwens

铁杆木虫 (正式写手)

老木虫


killl(金币+1):感谢参与应助,欢迎多来求助板块帮助别人。
引用回帖:
Originally posted by patent at 2007-1-15 12:54 AM:


这个我不是很明白,老兄,能否详细说明一下,用autodock具体怎样操作呢?
万分感谢!!!!!!!!!!!!!!!!

autodock中无法实现,记得你以前发的帖子好像说过你装了sybyl。
superimpose或者fit atoms。命令用法应该很简单。
万里夕阳锦背高 翻身犹恨东洋小 太公怎钓?
6楼2007-01-15 15:17:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjmed

木虫 (正式写手)


killl(金币+1):感谢参与应助,欢迎多来求助板块帮助别人。
在两个分子空间中,用命令
match m1 m2
Timeandtidewaitsfornoman!
7楼2007-01-15 19:56:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 patent 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见