蛋白質(zhì)折疊類型分類及其Profile HMM識(shí)別.pdf

1、隨著生物大分子數(shù)據(jù)庫(kù)中蛋白質(zhì)序列數(shù)目的增多,發(fā)展有效的方法,從氨基酸序列提取結(jié)構(gòu)信息成為后基因組時(shí)代的重要研究課題。越來(lái)越多的證據(jù)表明,天然蛋白質(zhì)的折疊類型在總數(shù)上的是有限的,一般認(rèn)為只有數(shù)百到數(shù)千種,遠(yuǎn)小于蛋白質(zhì)所具有的自由度數(shù)。Anfinsen原理指出蛋白質(zhì)的結(jié)構(gòu)在很大程度上由其序列決定,當(dāng)?shù)鞍踪|(zhì)結(jié)構(gòu)數(shù)據(jù)庫(kù)趨于完備以后,結(jié)構(gòu)的解析問(wèn)題就可以轉(zhuǎn)化為折疊識(shí)別問(wèn)題,即使用計(jì)算方法,找到與待預(yù)測(cè)蛋白質(zhì)序列在三維結(jié)構(gòu)上最匹配的已知折疊類型。對(duì)

2、自然界存在的數(shù)百到數(shù)千種折疊類型進(jìn)行系統(tǒng)研究,有助于揭示蛋白質(zhì)的折疊規(guī)律,可為大型的蛋白質(zhì)數(shù)據(jù)庫(kù)提供結(jié)構(gòu)注釋,或者為蛋白質(zhì)的精確結(jié)構(gòu)預(yù)測(cè)提供參考。 目前的蛋白質(zhì)折疊類型識(shí)別基本上都是靠專家來(lái)完成的,不同的庫(kù)分類頗不相同。SCOP通過(guò)觀察將蛋白質(zhì)按照同源性進(jìn)行分類,但在SCOP的一些折疊子中,二級(jí)結(jié)構(gòu)及其走向存在很大差異,為折疊識(shí)別的模型構(gòu)建造成困難；CATH以序列、結(jié)構(gòu)比對(duì)的相似性打分為依據(jù)進(jìn)行"Topology"的分類,并沒(méi)

3、有直接反映出蛋白質(zhì)之間二級(jí)結(jié)構(gòu)及其空間排布上的相似性。實(shí)際上,蛋白質(zhì)的折疊類型反映了蛋白核心結(jié)構(gòu)的拓?fù)浣Y(jié)構(gòu)模式,包括蛋白質(zhì)分子空間結(jié)構(gòu)的三個(gè)主要方面：二級(jí)結(jié)構(gòu)單元、二級(jí)結(jié)構(gòu)單元的相對(duì)排布位置以及蛋白質(zhì)多肽鏈的整個(gè)路由關(guān)系(即肽鏈走向),我們?cè)诘鞍踪|(zhì)折疊結(jié)構(gòu)研究的基礎(chǔ)上,以結(jié)構(gòu)域的拓?fù)洳蛔冃詾橐罁?jù),結(jié)合二級(jí)結(jié)構(gòu)片段的空間排列、取向特征和連接關(guān)系,進(jìn)行蛋白質(zhì)折疊類型分類,建立了LIFCA數(shù)據(jù)庫(kù),為蛋白質(zhì)折疊識(shí)別奠定了基礎(chǔ)。 折疊識(shí)

4、別的一個(gè)重要方面是建立折疊識(shí)別算法。目前,折疊識(shí)別的方法大體上可以分為三類：氨基酸序列的兩兩比較,如使用Blast和Fasta判斷序列之間的相似性；多序列建模,如Profile HMM方法；分類器,如神經(jīng)網(wǎng)絡(luò),支持向量機(jī)等。與兩兩比較方法相比,HMM建立了統(tǒng)一的模型,可以抓住一組同源序列的公共核心,因而對(duì)于那些在已知數(shù)據(jù)庫(kù)沒(méi)有高相似度模板的未知序列有更好的識(shí)別效果;與SVM等分類機(jī)器相比,雖然后者可能獲得較高的準(zhǔn)確率,但是Profile

5、 HMM有較為統(tǒng)一的構(gòu)架,保留了位點(diǎn)信息,有詳細(xì)的序列概形,與SVM相比更有助于對(duì)特定折疊類型進(jìn)行進(jìn)一步的分析和研究。 在本文中,我們?cè)谇捌谡郫B分類的基礎(chǔ)上,對(duì)α類、β類及α/β類中樣本數(shù)量較多的74種折疊類型進(jìn)行研究,利用結(jié)構(gòu)比對(duì)得到多序列比對(duì),繼而產(chǎn)生Profile HMM進(jìn)行識(shí)別,研究工作主要包括以下幾個(gè)方面： 1.蛋白質(zhì)折疊類型分類數(shù)據(jù)庫(kù)LIFCA的建立 選取序列同一性低于25％的2406個(gè)蛋白質(zhì),

6、包含了全α、全β,α/β三個(gè)結(jié)構(gòu)類的所有代表性序列,在蛋白質(zhì)折疊結(jié)構(gòu)研究的基礎(chǔ)上,以結(jié)構(gòu)域的拓?fù)洳蛔冃詾橐罁?jù),結(jié)合二級(jí)結(jié)構(gòu)片段的空間排列、取向特征和連接關(guān)系,進(jìn)行蛋白質(zhì)折疊類型分類,對(duì)于結(jié)構(gòu)域的序列、二級(jí)結(jié)構(gòu)等信息,提供了詳細(xì)的注釋。為蛋白質(zhì)折疊類型識(shí)別奠定基礎(chǔ)。 2.折疊類型的結(jié)構(gòu)比對(duì)研究 使用結(jié)構(gòu)比對(duì)算法對(duì)LIFCA進(jìn)行了同種折疊類型在結(jié)構(gòu)上的差異性研究,以結(jié)構(gòu)比對(duì)為基礎(chǔ)得到了折疊類型的多序列比對(duì)結(jié)果,用于建立折

7、疊類型的模型。 3.HMM模型數(shù)據(jù)庫(kù)的建立 對(duì)成員數(shù)目大于等于4并且結(jié)構(gòu)比對(duì)效果較為顯著的74個(gè)折疊類型分別建立HMM模型,組成隱馬爾科夫模型庫(kù)。使用非冗余的Astra11.65序列庫(kù)進(jìn)行識(shí)別檢驗(yàn),識(shí)別精度74.5％,并保持了很低的假陽(yáng)性率,識(shí)別效果比目前報(bào)導(dǎo)的一些方法識(shí)別效果均要好。 本文從數(shù)據(jù)集篩選及算法兩個(gè)方面對(duì)蛋白質(zhì)的折疊類型識(shí)別方法進(jìn)行了改進(jìn),建立的隱馬爾科夫模型庫(kù)覆蓋范圍較廣,識(shí)別準(zhǔn)確率高,為折