| 5 | 1/1 | 返回列表 |
| 查看: 4930 | 回復(fù): 8 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[交流]
【求助】蛋白質(zhì)的RMSD均方根差指的是什么 已有6人參與
|
||||
amber | gromacs |
金蟲 (文壇精英)
老漢一枚
|
簡(jiǎn)單地說,就是蛋白質(zhì)上的原子們,在各個(gè)時(shí)間上的位置的變化情況。 原子i在t1時(shí)的位置 - 原子i在t2時(shí)的位置,得到原子i的位置變化 計(jì)算其平方 ..... 把多個(gè)原子(i=1...N)的這個(gè)平方,都加起來 除以原子數(shù),得到所謂“均方” 開方 搞成平方,是為了避免出現(xiàn)負(fù)值 再開方,是為了得到長(zhǎng)度的量綱 一種簡(jiǎn)單的小技巧而已 你也可以自己發(fā)明類似的算法 [ Last edited by yalefield on 2010-5-12 at 12:05 ] |
金蟲 (著名寫手)

金蟲 (文壇精英)
老漢一枚
|
http://www.sciencenet.cn/m/user_content.aspx?id=277393 ---> 用來表是蛋白質(zhì)結(jié)構(gòu)之間差異的參數(shù)是兩個(gè)結(jié)構(gòu)之間原子位置的 RMSD。 ---> 計(jì)算RMSD時(shí),可以針對(duì)目標(biāo)蛋白質(zhì)(如: 所有的原子、骨干部份或只考慮 alpha 碳原子等等)。不同的標(biāo)準(zhǔn),計(jì)算RMSD的數(shù)值會(huì)有所差異。 ---> RMSD 距離函數(shù),以一個(gè)結(jié)構(gòu)中的原子與另外一個(gè)結(jié)構(gòu)中對(duì)應(yīng)原子為計(jì)算標(biāo)的,因此,如果兩個(gè)分子在座標(biāo)系統(tǒng)中以不同的位置開始計(jì)算,那么不管其結(jié)構(gòu)是否相似,這兩者之間的 RMSD 必定相當(dāng)大。也因?yàn)檫@樣,我們?yōu)榱艘?jì)算有意義的 RMSD ,兩者的結(jié)構(gòu)要盡可能的重疊。 對(duì)于docking而言,如果有reference ligand,一般不需要額外的重疊,否則會(huì)有偽造數(shù)據(jù)之嫌。 ---> 可以通過計(jì)算 RMSD 來當(dāng)作評(píng)估蛋白質(zhì)結(jié)構(gòu)的可信度: 在模擬過程中,分子會(huì)不斷的發(fā)生變化,而對(duì)于我們而言,必須等到分子結(jié)構(gòu)在穩(wěn)定的狀態(tài)下(fluctuation較小時(shí))再進(jìn)一步進(jìn)行分析,這樣才比較有意義。 對(duì)于序列和長(zhǎng)度不同的蛋白結(jié)構(gòu),比較RMSD似乎意義不大。 對(duì)于存在序列和長(zhǎng)度差異的蛋白,首先我們不知道軟件本身在計(jì)算RMSD的時(shí)候究竟采用了何種計(jì)算方法:所有原子?骨架部分?還是CA?這樣就很難給我們一個(gè)明確的概念。 另外,對(duì)于長(zhǎng)短不一的蛋白,軟件首先會(huì)對(duì)序列進(jìn)行比較,本身比對(duì)的方法不同,顯然會(huì)給后面的結(jié)果帶來影響。 對(duì)于長(zhǎng)短不一的蛋白,軟件究竟取多少殘基作為計(jì)算RMSD的對(duì)象,我們也無從知道。 所以,對(duì)于不是同一個(gè)蛋白的RMSD的計(jì)算,其中的不確定性比較多,計(jì)算出來的結(jié)果只能作為參考,給我們一個(gè)大概的概念。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +3 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +9 | 冬十三 2026-03-15 | 10/500 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑,總分315,考的生物醫(yī)藥,一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|