Needleman-Wunsch 算法是一種全局聯(lián)配算法,它從整體上分析兩個序列的關(guān)系,即考慮序列總長的整體比較,用類似于使整體相似(global similarity)最大化的方式,對序列進行聯(lián)配。兩個不等長序列的聯(lián)配分析必須考慮在一個序列中一些堿基的刪...[繼續(xù)閱讀]

    Smith-Waterman 算法是在Needleman-Wunsch 算法基礎(chǔ)上發(fā)展而來的,它是一種局部聯(lián)配算法。由于親緣關(guān)系較遠的蛋白質(zhì)序列可能只有一些相互獨立的保守片段,所以進行局部相似性分析有時可能比整體相似性分析更合理。Smith 和Waterman(1981)提...[繼續(xù)閱讀]

    BLAST 算法同樣是利用動態(tài)規(guī)劃算法,與Smith-Waterman 算法類似,其不同之處是引入了所謂“字”或“字符串”(word 或K-tuple,K-mer 等)的檢索技術(shù)。所有序列其實都是由若干字符串組成,例如我們以3個堿基長度的字符串為例,下列DNA 序列包括...[繼續(xù)閱讀]

    采用BLAST 的基本算法目前形成了若干不同的工具,分別用于特定序列數(shù)據(jù)庫和特定目的的序列搜索。以NCBI 提供的在線序列數(shù)據(jù)庫搜索工具BLAST(2014年12月)為例(圖1-3.4),BLASTN 是對核苷酸遞交序列庫搜索核苷酸序列數(shù)據(jù)庫,BLASTP 是在蛋白...[繼續(xù)閱讀]

    BLAST 搜索返回的結(jié)果中,提供了遞交序列與數(shù)據(jù)庫中序列比對結(jié)果的得分(score)和一個統(tǒng)計測驗結(jié)果(E-value)。到目前為止,對局部聯(lián)配的統(tǒng)計學問題已基本搞清楚,特別是那些不含空位的局部聯(lián)配更是如此。我們不妨首先考慮不含空位的...[繼續(xù)閱讀]

    1.目前計分矩陣主要有哪些? 比較它們的異同。2.請利用動態(tài)規(guī)劃Needleman-Wunsch 算法對下列兩條蛋白質(zhì)序列進行全局聯(lián)配,獲得最優(yōu)聯(lián)配結(jié)果:　　P1=AGWGAHEA　　P2=PAWHEAEAG計分系統(tǒng):計分矩陣BLOSUM50,空位罰8分。表1-3.7　BLOSUM50(部分)3.數(shù)據(jù)庫...[繼續(xù)閱讀]

    許多生物學研究都涉及多條序列甚至幾十上百條序列的比較,因此多序列聯(lián)配是生物信息學的一個重要課題。通過多序列聯(lián)配結(jié)果,我們可以確定這些序列的親緣關(guān)系,通過序列保守性判斷功能域或功能位點等。多序列聯(lián)配同樣包括全...[繼續(xù)閱讀]

    多序列全局聯(lián)配算法目前主要是以Clustal 算法為代表的漸進式全局聯(lián)配方法。漸進式全局聯(lián)配算法是20世紀80年代發(fā)展起來的(Waterman 和Perlwitz,1984; Feng 和Doolittle,1987),其中以軟件工具Clustal 算法最為成功。下面就重點介紹Clustal 算法。...[繼續(xù)閱讀]

    具有相同功能的基因往往在序列上存在局部相似性或保守性,這些保守性與相應(yīng)功能和選擇壓等有關(guān)(詳見下節(jié))。圖1-4.3列舉了一個多序列局部保守性。生物信息學的一個重要任務(wù)是找到這些保守序列。圖1-4.3　多序列局部保守性舉例...[繼續(xù)閱讀]

    蛋白質(zhì)功能域(domain)的概念最早由Wetlaufer(1973)在研究蛋白質(zhì)結(jié)構(gòu)時提出。蛋白質(zhì)功能域一般是指一條蛋白質(zhì)序列中一段保守的區(qū)域,該區(qū)域能夠獨立行使功能、進化等。在蛋白質(zhì)結(jié)構(gòu)中,功能域是指一個蛋白質(zhì)結(jié)構(gòu)的一部分,它能形成一...[繼續(xù)閱讀]