24小时热门版块排行榜    

查看: 3825  |  回复: 12

zpliu09

金虫 (小有名气)

[求助] 关于两个蛋白质结构叠合的原理(或者相关的程序)

课题在做蛋白质结构预测这块,经常用到两个蛋白质结构的叠合,然后算出RMSD,用来体现两个蛋白质之间结构的相似程度。但是这个叠合过程究竟是如何完成的?什么原理?
我目前只了解到大概涉及到旋转矩阵方面,坐标变换方面的知识,哪位了解的朋友可以指点下究竟是怎么做的?
有什么参考资料吗?
先感谢各位
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

zpliu09

金虫 (小有名气)

6楼2011-06-28 22:33:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)


我觉得算法过程是这样的:

对于任意给的两个多肽链,在空间任意初始位置,我们总可以计算对应骨架上的原子坐标的偏差平方和MSD,然后不断改变两肽链的空间位置,直到此MSD最小。MSD=sum[(x1i-x2i)^2]   (1,2指代两肽链,i原子序号)

问题有2个:

(1) 如何改变多肽链的空间位置。这个简单,固定一个肽链1,对另外一个肽链2定义一个平移变量,一个旋转变量。肽链2的空间位置变换无非是当前坐标矩阵乘以这2个变量;多肽链的移动就是不断改变这两个变量而已。而总MSD就是这两个变量的函数。

(2) 如何找对应原子。对于两个有差别的肽链,第一步必须确定对应残基。这个有Alignment算法,其实,就是把两肽链的残基字母顺序按某种对应方式摆在一起,然后如果双发字母相同的,加1分,不同的记为0分,反复几次,即可找到最佳的对应方式。确定了两肽链的残基一一对应方式以后,就可以对骨架C-C-N-原子坐标间的MSD计算了。

还有其它问题,最好的alignment(摆放方式)未必对应全局最小的MSD,那么怎么办呢?。。。我想你已经想到了,:-)
9楼2011-07-01 11:55:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yinhui98

金虫 (小有名气)

【答案】应助回帖


youzhizhe(金币+1): 谢谢交流。 2011-05-10 12:13:46
发表愚见仅供参考,呵呵。
1、蛋白质折叠过程及原理:主要由于键与键的相互作用、范德华力、氨基酸的疏水性等等。
2、可以用复杂网络去预测。
3、你可以在ISI knowledge(http://pcs.isiknowledge.com/Sess ... ver.sessionNotFound)网页上,输入protein folding,就会发现很多相关文献。
2楼2011-05-10 09:12:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpliu09

金虫 (小有名气)

yinhui98,很感谢你的关注。但是我的问题不是蛋白质折叠,是两个结构进行叠合,找出相似和不同的区域。叠合主要用于评估两个结构的相似程度。和你说的折叠过程不是一回事。。。所以不好意思,就不给金币了,呵呵
3楼2011-05-10 20:56:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cedarice

木虫 (著名写手)

【答案】应助回帖


youzhizhe(金币+1): 谢谢提示。 2011-05-12 18:23:13
zpliu09(金币+2): 谢谢,但是貌似我的问题还是没有解决。呵呵。这个结构比对的范围有点大了 2011-06-26 22:23:38
搜搜Protein Structure Alignment试试看
生活就像忐忑,没有准确的歌词,却惊心动魄。
4楼2011-05-12 17:11:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cedarice

木虫 (著名写手)

这个有相应的软件可以实现,如chimera,但是不清楚能否指定相应位置进行alignment。另外可以在参考Wiki上的这个介绍:http://en.wikipedia.org/wiki/Structural_alignment
生活就像忐忑,没有准确的歌词,却惊心动魄。
5楼2011-06-26 22:27:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)


贺仪(金币+1): 多谢你的热心帮助! 2011-06-29 21:23:20
引用回帖:
Originally posted by zpliu09 at 2011-06-28 22:33:13:
谢谢你一直的关注阿!
如何使用软件算出RMSD这方面我还是有一点了解。我现在做的蛋白质结构预测,需要自己编程对蛋白数据进行叠合处理。所以没法直接用像你说的这类软件。
我是希望最好能找到源代码(最好是 ...

VMD就可以干。只需一个命令即可。
在RMSD calculator中,alignment。

具体细节涉及到平移和旋转,可以看VMD里相关部分的code。
这一块可以去他们的mailing list上边提问。

» 本帖已获得的红花(最新10朵)

7楼2011-06-29 12:23:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpliu09

金虫 (小有名气)

送鲜花一朵
引用回帖:
Originally posted by yongleli at 2011-06-29 12:23:08:
VMD就可以干。只需一个命令即可。
在RMSD calculator中,alignment。

具体细节涉及到平移和旋转,可以看VMD里相关部分的code。
这一块可以去他们的mailing list上边提问。

十分感谢!!!
你提的方法应该比较可行。我去看看。
再次感谢!!
8楼2011-06-30 20:46:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smutao

禁虫 (著名写手)

本帖内容被屏蔽

10楼2015-04-03 03:23:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zpliu09 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +5 妮妮ninicgb 2026-03-15 8/400 2026-03-16 14:44 by 可淡不可忘
[考研] 材料专硕306英一数二 +4 z1z2z3879 2026-03-16 4/200 2026-03-16 13:53 by laoshidan
[考研] 本人考085602 化学工程 专硕 +9 不知道叫什么! 2026-03-15 10/500 2026-03-16 10:04 by houyaoxu
[考研] 321求调剂 +4 大米饭! 2026-03-15 4/200 2026-03-16 08:41 by Linda Hu
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 285化工学硕求调剂(081700) +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见