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

金蟲(chóng) (文壇精英)
老漢一枚
|
http://www.sciencenet.cn/m/user_content.aspx?id=277393 ---> 用來(lái)表是蛋白質(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)中以不同的位置開(kāi)始計(jì)算,那么不管其結(jié)構(gòu)是否相似,這兩者之間的 RMSD 必定相當(dāng)大。也因?yàn)檫@樣,我們?yōu)榱艘?jì)算有意義的 RMSD ,兩者的結(jié)構(gòu)要盡可能的重疊。 對(duì)于docking而言,如果有reference ligand,一般不需要額外的重疊,否則會(huì)有偽造數(shù)據(jù)之嫌。 ---> 可以通過(guò)計(jì)算 RMSD 來(lái)當(dāng)作評(píng)估蛋白質(zhì)結(jié)構(gòu)的可信度: 在模擬過(guò)程中,分子會(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é)果帶來(lái)影響。 對(duì)于長(zhǎng)短不一的蛋白,軟件究竟取多少殘基作為計(jì)算RMSD的對(duì)象,我們也無(wú)從知道。 所以,對(duì)于不是同一個(gè)蛋白的RMSD的計(jì)算,其中的不確定性比較多,計(jì)算出來(lái)的結(jié)果只能作為參考,給我們一個(gè)大概的概念。 |
金蟲(chóng) (文壇精英)
老漢一枚
金蟲(chóng) (正式寫(xiě)手)

銀蟲(chóng) (初入文壇)

新蟲(chóng) (初入文壇)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|