24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4070  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:30 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:15 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 3/150 2026-05-16 23:57 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 xx7gd5zq4e 2026-05-15 5/250 2026-05-16 23:57 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 7/350 2026-05-16 23:17 by ue3ir18jc3
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 23:10 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 22:55 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-16 22:25 by lindoudou
[考博] 现在不知道怎么办,感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +6 水和泥不是水泥 2026-05-13 8/400 2026-05-16 19:50 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 2026博士还有哪些学校有名额 +5 小王求读研 2026-05-15 6/300 2026-05-16 10:44 by a4742549
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[基金申请] 青C资助名额大幅增加! +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
信息提示
请填处理意见