生物信息学 第七章:序列比对和数据库搜索
生物信息学中的序列比对技术分析

生物信息学中的序列比对技术分析随着生物技术的不断进步,自动化测序技术的快速发展,大量生物学数据呈爆炸式增长。
同时,对生物信息学分析的需求日益增大,序列比对则成为生物信息学最常见的分析手段之一。
序列比对技术可以对已知序列与未知序列进行匹配、比对,以找出其中的异同点,分析其功能和演化关系,是生物科学、基因组学等分支的核心技术之一。
1. 序列比对的基本概念序列比对是指将两个或多个序列进行对比,找出它们的相似和不同之处的过程。
从基本原理上讲,序列比对是将一条DNA或RNA序列与另一条同源序列进行匹配的过程,而通过比较相同和不同之处来推断它们可能存在的共同祖先。
所谓同源序列,指的是两个或多个序列具有较高的序列相似度,可能来自相同种属的生物体或同一基因家族中的不同基因成员。
同源序列对于了解分子进化、基因结构与功能以及物种关系具有重要的意义。
2. 序列比对的类型在生物信息学领域,基本可以将序列比对分为全局比对和局部比对两种。
(1)全局比对全局比对是指将整个序列与另一条序列进行比对,寻找全长匹配区域。
全局比对适用于已知的高度同源性序列分析。
最常用的全局比对算法包括 Needleman-Wunsch 和 Smith-Waterman 算法。
其中,Needleman-Wunsch 算法较为严谨,适用于匹配全长序列;而 Smith-Waterman 算法则更为灵活,可以匹配任意长度的序列片段,并且可以找到更为相似的匹配序列。
(2)局部比对局部比对是指只比对序列中一部分序列,而不需要考虑整个序列,寻找相似或同源的序列区间。
相较于全局比对,局部比对更适合用于寻找序列中比较短且高度相似的区域。
常用的局部比对算法有 BLAST (Basic Local Alignment Search Tool) 和 FASTA (Fast Alignment Search Tool) 算法。
这些算法适用于较长的未知序列与基因或蛋白质序列数据库进行比对。
生物信息学中的序列比对算法使用方法解析

生物信息学中的序列比对算法使用方法解析序列比对在生物信息学中是一项重要的技术,用于寻找DNA、RNA或蛋白质序列之间的相似性和差异性。
它是理解生物学结构和功能的基石之一。
在本文中,我们将解析生物信息学中常用的序列比对算法的使用方法。
序列比对算法主要分为全局比对和局部比对。
全局比对用于比较完整的序列,而局部比对则更适用于在序列中查找相似区域。
在这两个主要类别中,有几种经典的序列比对算法,包括Pairwise Sequence Alignment、BLAST、Smith-Waterman算法和Needleman-Wunsch算法等。
首先,我们来看Pairwise Sequence Alignment(两两序列比对)算法。
这个算法是基本的序列比对方法,通过比较两个序列中的每一个碱基、氨基酸或核苷酸,并根据其相似性和差异性对它们进行排列。
Pairwise Sequence Alignment算法使用动态规划的思想,通过计算匹配、替代和插入/删除的分数,来确定两个序列的最佳匹配方案。
在生物信息学中,常用的实现包括Needleman-Wunsch算法和Smith-Waterman算法。
Needleman-Wunsch算法是一种全局比对算法,用于比较两个序列的整个长度。
它是通过填充一个二维矩阵来计算最佳匹配路径的。
算法的核心思想是,通过评估每个格子的分数,根据路径选择的最佳分数进行全局比对。
这个算法不仅可以计算序列的相似性,还可以计算每个位置的分数,从而获得两个序列的对应二面的对应关系。
Smith-Waterman算法是一种局部比对算法,用于寻找两个序列中的最佳匹配片段(子序列)。
它与Needleman-Wunsch算法的计算思路相同,但不同之处在于允许负分数,这使得算法能够确定具有高分数的局部匹配片段。
通过动态规划计算,Smith-Waterman算法可以寻找到两个序列中的相似片段,并生成比对的结果。
另一种常用的序列比对算法是基本本地搜索工具(BLAST)。
生物信息学中的序列比对算法综述

生物信息学中的序列比对算法综述序列比对(sequence alignment)是生物信息学中一项重要的任务,其目的是找出两个或多个生物序列中的相似性和差异性。
在生物信息的研究和应用中,序列比对算法起到了至关重要的作用。
本文将对生物信息学中的序列比对算法进行综述。
1. 引言序列比对是生物信息学中的一个基本问题,它在基因组学、蛋白质学、进化生物学等领域都得到了广泛的应用。
通过比对不同生物序列之间的相似性和差异性,可以进一步研究基因功能、蛋白质结构以及物种进化等重要问题。
因此,序列比对算法的研究具有重要的理论价值和实际意义。
2. 序列比对的基本概念在进行序列比对之前,首先需要了解序列之间的相似性和差异性的度量方法。
常用的序列相似性度量方法包括编辑距离、相似度百分比、贝叶斯统计等。
其中,编辑距离是一种常见的度量方式,它衡量了两个序列之间的差异程度。
3. 序列比对算法分类序列比对算法可以分为全局比对和局部比对两类。
全局比对算法着重于找出整个序列的相似性和差异性,常用的算法包括Needleman-Wunsch算法和Smith-Waterman算法。
而局部比对算法则注重于找出序列中的局部相似性和差异性,常用的算法有BLAST和FASTA。
4. 全局比对算法全局比对算法的核心思想是将两个序列通过插入、删除和替换等操作转化为相同长度的序列,然后计算它们的相似性得分。
Needleman-Wunsch算法是一种经典的全局比对算法,通过动态规划的方式找到序列之间的最佳比对方式。
Smith-Waterman算法是基于Needleman-Wunsch算法的改进,它将负得分和局部比对引入到全局比对中,提高了比对的准确性。
5. 局部比对算法局部比对算法主要用于序列中的片段比对,其核心思想是通过寻找序列中的相似片段来找出序列的结构和功能区域。
BLAST算法是一种常用的局部比对算法,它通过生成字典和索引的方式实现快速比对。
FASTA算法则是一种早期的局部比对算法,其基本原理是通过序列片段之间的kmer匹配来寻找相似性。
生物信息学中的序列比对算法综述

生物信息学中的序列比对算法综述序列比对是生物信息学领域中的一个重要问题,指的是比较两个生物序列(DNA,RNA或蛋白质序列)之间的相似性和差异性。
序列比对是许多研究任务中的第一步,如基因识别、物种分类、进化关系的推断等等。
在本文中,我们将介绍序列比对算法的基本概念、方法和软件,包括全局比对、局部比对、多序列比对等方面。
一、序列比对的基本概念序列比对的目的是找出两个序列之间的相似性和差异性,根据相似性分析序列的结构、功能以及进化关系。
相似性可以被表示成一个比对得分,即正数表示相似性,负数表示差异性。
比对得分的计算取决于匹配分、替换分和缺失分。
匹配分是指在比对中找到相同的位置并且相等的分数。
替换分是指找到不同的位置并且不相等的分数。
缺失分是指在任意序列中找不到匹配的分数。
计算得分的方法有很多种,其中最流行的方法是 Needleman-Wunsch 算法和 Smith-Waterman 算法。
二、全局比对算法全局比对算法是一种比较两个序列的整个长度的算法,使得它们之间的相似性或差异性能够被准确地测量。
全局比对算法通常用于比较高度相似的序列或同一物种中相似的序列。
Needleman-Wunsch 算法与 Smith-Waterman 算法是全局比对中最为经典的算法。
Needleman-Wunsch 算法: Needleman-Wunsch 算法是最经典的全局比对算法之一。
该算法通过构建一个二维矩阵,其中每个元素代表在比对过程中两个序列的一个指定位置。
该算法通过分配一个比对得分并使用动态规划来计算所有可能的比对方式。
通过比对得分的计算,算法确定序列之间的最佳比对方式,使比对得分最大化。
该算法常用于比较高度相似的序列,或者已知序列的情况下以寻找相同物种中潜在基因组之间的相似性信息。
Smith-Waterman 算法: Smith-Waterman 算法是一种类似Needleman-Wunsch 算法的全局比对算法。
生物信息学中的序列比对和分析

生物信息学中的序列比对和分析序列比对和分析是生物信息学中非常基础和重要的一项研究内容。
通过比对和分析序列,可以发现序列之间的相似性和差异性,进而研究生物进化、遗传、表达等方面的问题。
本文将从序列比对和分析的意义、比对方法、分析工具和应用实例等几个方面进行讨论。
一、序列比对的意义和方法序列比对是一种比较两个或多个序列相似度的方法,通过比较序列的相同和不同部分,可以获得有关序列功能、结构和进化的信息。
序列比对的主要目的是确定两个序列之间的相似性程度,从而推断它们的共同祖先、结构和功能。
因此,序列比对是研究生物学、医学和生物工程等领域的必要手段。
序列比对的方法主要包括全局比对和局部比对两种。
全局比对是将一整个序列与另一个序列比对,得到两序列的整体相似性程度。
一般要求两序列中的相似部分要尽可能多,而不注重不同部分的对齐。
常用的全局比对算法有Needleman-Wunsch算法和Smith-Waterman算法。
局部比对是寻找两个序列中任意长度的子序列之间的相似性。
与全局比对不同,局部比对更注重相同的局部片段,忽略不同的片段。
局部比对算法有BLAST和FASTA等。
二、序列分析的工具和方法序列分析是通过对比对后的序列进行进一步分析,获得生物信息的过程。
序列分析的主要内容包括序列注释、序列搜索、序列聚类和序列比较等。
序列分析的方法和工具主要包括多序列比对、单序列比对、序列搜索、聚类分析、进化树分析和功能预测等。
多序列比对是将多个序列进行比对,得到这些序列之间的共同特征和差异。
常用的多序列比对工具有Clustal X和MUSCLE等。
单序列比对是将一个序列和已知的库中所有序列进行比对,以查找相似性和相关性。
常用的单序列比对工具有BLAST和PSI-BLAST等。
序列搜索是在一个已知的序列库中搜索相似的序列。
常用的工具有HMMER、PhyloGenie等。
聚类分析是将相似的序列放在一起形成聚类,便于分析相关性。
生物信息学中的序列比对方法

生物信息学中的序列比对方法序列比对是生物信息学中一项非常重要的工具,其主要目的是将两个或更多的DNA、RNA或蛋白质序列进行比较,以找到它们之间的相似性和差异性。
这样的比对可以用来识别基因、预测蛋白质结构、推断进化关系和研究生物系统的复杂性等。
随着DNA测序技术的快速发展,越来越多的生物学家和生物信息学家开始研究序列比对方法。
序列比对是一项复杂而耗时的任务,需要对大量的序列进行计算和分析。
因此,发展高效的序列比对方法对于生物信息学的发展至关重要。
当前,生物信息学界广泛应用的序列比对方法主要包括全局比对、局部比对和多序列比对等。
一、全局比对全局比对是指将整个序列与另一个相似序列进行比对。
它的应用场景通常是在两个相对较短的序列中查找相似片段,以便在进一步的研究中进行详细的分析。
全局比对方法最常用的是Needleman-Wunsch算法和Smith-Waterman算法。
Needleman-Wunsch(NW)算法是第一个被开发出来的全局比对算法。
该算法基于动态编程的思想,通过将整个序列进行比对,计算出最佳匹配的得分和路径。
然而,这种方法的时间复杂度非常高,随着序列长度的增加,其计算成本也会呈指数级增长。
Smith-Waterman(SW)算法是一种优化的全局比对算法,其核心思想与NW算法类似。
不同之处在于SW算法将匹配的得分设置为正数,而将多余的间隔和未匹配的子序列得分设置为负数。
通过这种方式,SW算法可以得到一个全局最佳的比对结果。
然而,该算法的计算成本也比较高,因此其应用场景受到一定的限制。
二、局部比对局部比对是指在比对序列的过程中,只对部分区域进行比对。
与全局比对不同,局部比对更适用于两个序列之间只有一些片段相似的情况。
常用的局部比对方法主要包括BLAST算法和FASTA算法等。
BLAST算法是一种聚集序列算法,它将大量的搜索序列放入一个空间中,通过加速计算找到最匹配的序列。
通过BLAST算法,可以快速搜索数据库中的所有序列,并找到与目标序列相似的匹配。
生物信息学 第七章:序列比对和数据库搜索

第七章:序列比对和数据库搜索Gregory D.SchulerNational Center for Biotechnology InformationNational Library of Medicine. National Institutes of HealthBethesda. Maryland引言在生物学的研究中,有一个常用的方法,就是通过比较分析获取有用的信息和知识。
达尔文正是研究比较了galapagos finches同其它一些物种的形态学特征,从而提出了自然选择学说。
今天,我们对基因和蛋白质序列进行比较,从本质上来讲是同达尔文一样,进行同样的分析,只不过更加精细,更加详尽。
在这个意义上,我们从核酸以及氨基酸的层次去分析序列的相同点和不同点,以期能够推测它们的结构、功能以及进化上的联系。
最常用的比较方法是序列比对,它为两个或更多个序列的残基之间的相互关系提供了一个非常明确的图谱。
在这一章,我们只讨论一下双重比对,即只比较两个序列,至于较多的序列即多序列比对,将在第八章介绍。
七十年代以来,DNA测序方法的飞速发展,极大地引发了序列信息量的扩增,从而使可供比较的序列数量呈现爆炸式增长。
分子生物学家应该意识到,将未知序列同整个数据库中的已知序列进行比较分析已经成为他们手中一个强有力的研究手段。
在过去的三十年里,即使不提及计算机的应用,序列比较的各种算法也已经发展得越来越迅速,也越来越成熟,已经能够跟上序列数据库增长的步伐。
今天,我们已经拥有一些小的模式物种的基因组的全序列,还拥有人类基因序列的一些较大的样品,我们已经进入比较基因组时代,也就是说,对两个物种进行全基因组序列比较已经不再是一个梦想。
序列比对的进化基础进行序列比对的目的之一是让人们能够判断两个序列之间是否具有足够的相似性,从而判定二者之间是否具有同源性。
值得注意的是,相似性和同源性虽然在某种程度上具有一致性,但它们是完全不同的两个概念。
生物信息学中的序列分析与比对

生物信息学中的序列分析与比对生物信息学是一门综合运用计算机科学、数学、物理学、化学、生物学等学科方法来研究生物信息学的交叉学科。
生物信息学的研究方法多样,其中序列分析与比对是其中重要的一个分支。
序列是指生物学中的一种信息记录方式,表示基因、蛋白质、RNA等分子的特定序列。
序列分析和比对是研究序列信息的一种重要手段,也是对生物大分子结构、功能、进化等各方面研究的基础。
本文将从序列概述、序列分析、序列比对三方面来详细讨论序列分析与比对的内容。
一、序列概述DNA序列、RNA序列、蛋白质序列是生物学中最基础的三种序列类型。
DNA是一种双螺旋结构,由四种碱基(腺嘌呤、鸟嘌呤、胸腺嘧啶和鳞状细胞嘧啶)组成, RNA是单链结构,基本上由DNA中的碱基组成,而蛋白质则是由20种不同的氨基酸通过共价键结合在一起形成的链状分子。
DNA序列和RNA序列是由碱基组成的线性序列,蛋白质序列是由氨基酸组成的线性序列,序列是描述生物大分子结构和功能的最基本的语言。
二、序列分析1. 序列比较序列比对是对两条或两条以上生物序列之间相似性或差异性进行分析的方法。
序列比对是将两条或多条序列中相同或类似的部分找出来,同时也可以找出它们中不同或不同的部分,并计算它们之间的相似性和差异性,也是比较序列进化关系的一种方法。
2. 基因预测基因预测是确定核酸序列中哪些区域包含已知的基因,并预测这些基因与蛋白质的长链和功能的相关性质。
序列分析的目标之一就是预测基因和蛋白质的序列。
基因预测估计原核生物中基因组序列的含义比真核生物要容易得多。
3. 重复序列分析许多基因组中都存在着许多多次出现的特定、比较长的DNA 序列,这种序列被称为重复序列。
重复序列的分析是基因组学的必要内容,使人们理解基因组进化和功能等方面的信息。
通过对重复序列的分析,人们可以了解基因组结构的重要细节,解决许多生物学问题,比如基因家族的起源,基因组的演化,基因转座子活动和某些疾病和种群的进化关系等。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第七章:序列比对和数据库搜索Gregory D.SchulerNational Center for Biotechnology InformationNational Library of Medicine. National Institutes of HealthBethesda. Maryland引言在生物学的研究中,有一个常用的方法,就是通过比较分析获取有用的信息和知识。
达尔文正是研究比较了galapagos finches同其它一些物种的形态学特征,从而提出了自然选择学说。
今天,我们对基因和蛋白质序列进行比较,从本质上来讲是同达尔文一样,进行同样的分析,只不过更加精细,更加详尽。
在这个意义上,我们从核酸以及氨基酸的层次去分析序列的相同点和不同点,以期能够推测它们的结构、功能以及进化上的联系。
最常用的比较方法是序列比对,它为两个或更多个序列的残基之间的相互关系提供了一个非常明确的图谱。
在这一章,我们只讨论一下双重比对,即只比较两个序列,至于较多的序列即多序列比对,将在第八章介绍。
七十年代以来,DNA测序方法的飞速发展,极大地引发了序列信息量的扩增,从而使可供比较的序列数量呈现爆炸式增长。
分子生物学家应该意识到,将未知序列同整个数据库中的已知序列进行比较分析已经成为他们手中一个强有力的研究手段。
在过去的三十年里,即使不提及计算机的应用,序列比较的各种算法也已经发展得越来越迅速,也越来越成熟,已经能够跟上序列数据库增长的步伐。
今天,我们已经拥有一些小的模式物种的基因组的全序列,还拥有人类基因序列的一些较大的样品,我们已经进入比较基因组时代,也就是说,对两个物种进行全基因组序列比较已经不再是一个梦想。
序列比对的进化基础进行序列比对的目的之一是让人们能够判断两个序列之间是否具有足够的相似性,从而判定二者之间是否具有同源性。
值得注意的是,相似性和同源性虽然在某种程度上具有一致性,但它们是完全不同的两个概念。
相似性是指一种很直接的数量关系,比如部分相同或相似的百分比或其它一些合适的度量,而同源性是指从一些数据中推断出的两个基因在进化上曾具有共同祖先的结论,它是质的判断。
基因之间要么同源,要么不同源,绝不象相似性那样具有多或少的数量关系。
如图7.1所示,比较家鼠和小龙虾的同源的胰蛋白酶序列,发现它们具有41%的相似性。
由于受到研究进化关系这一目的的影响,大多数比对方法很自然地都希望能够在某种程度上建立起分子进化的模型。
我们通常都假定同源序列是从某一共同祖先不断变化而来,但事实上,我们无法得知这个祖先序列到底是什么样子,除非能够从化石中获得它的DNA,我们所能够做到的只是从现存物种中,探求真相。
从祖先序列以来所发生的变化包括取代、插入以及缺失。
在理想情况下,同源基因或蛋白质序列在相互比较时,残基之间相互对应,从而使取代的情况很明显地表现出来。
在某些位置,一个序列中拥有某些残基而另一个序Bioinformatics: A Practical Guide to the Analysis of Genes and ProteinsEdited by A.D.Baxevanis and B.F.F.OuelletteISBN 0-471-19196-5.pages 145-171. Copyright© 1998 Wiley-Liss. Inc.列中缺少这种残基,表明这些残基是插入到前者或是从后者中丢失的。
这些空位在序列比对时用连续的短线填补。
如图7.1,在序列比对中,发现了5个空位。
|------ S-S-------*|Mouse IVGGYNCEENSVPYQVSLNS-----GYHFCGGSLINEQWVVSAGHCYK-------SRIQVCrayfish IVGGTDAVLGEFPYQLSFQETFLGFSFHFCGASIYNENYAITAGHCVYGDDYENPSGLQI*Mouse RLGEHNIEVLEGNEQFINAAKIIRHPQYDRKTLNNDIMLIKLSSRAVINARVSTISLPTACrayfish VAGELDMSVNEGSEQTITVSKIILHENFDYDLLDNDISLLKLSGSLTFNNNVAPIALPAQ |---- S-S--------|Mouse PPATGTKCLISGWGNTASSGADYPDELQCLDAPVLSQAKCEASYPG-KITSNMFCVGFLE Crayfish GHTATGNVIVTGWG-TTSEGGNTPDVLQKVTVPLVSDAECRDDYGADEIFDSMICAGVPE ◇ *|-------------S-S------------------|Mouse GGKDSCQGDSGGPVVCNG----QLQGVVSWGDGCAQKNKPGVYTKVYNYVKWIKNTIAAN Crayfish GGKDSCQGDSGGPLAASDTGSTYLAGIVSWGYGCARPGYPGVYTEVSYHVDWIKANAV-- 图7.1、保守位点通常在功能上极为重要。
对老鼠的胰蛋白酶(Swiss-ProtP07146)和小龙虾的胰蛋白酶(Swiss-Prot P00765)作比对,相同的残基用下标线标出,在比对上方标出的是三个二硫键(-S-S),这些二硫键中的半胱氨酸残基极为保守,打星号的残基的侧链参与电荷传递系统,打菱形符号的活性位点的残基负责底物的特异性。
在残基-残基比对中,很明显,某些位置的氨基酸残基相对于其它位置的残基具有较高的保守性,这个信息揭示了某些残基对于一个蛋白质的结构和功能是极为重要的。
如图7.1所示,处于活性位点的残基都是极为保守的,比如形成二硫键的半胱氨酸,参与电子传递的氨基酸残基以及决定底物特异性的氨基酸残基。
这些保守的残基对于保持蛋白的结构与功能非常重要,另一方面,由于历史原因,某些保守位置对蛋白功能并无太大的重要性。
当我们处理非常相近的物种时必须十分小心,因为相似性在某些情况下更多地是历史的反映而不是功能的反映,比如,mouse和rat的某些序列具有高度的相似性,可能仅仅是因为没有足够的时间进行分化而已。
尽管如此,系列比对仍然是从已知获得未知的一个十分有用的方法,比如通过比较一个新的蛋白同其它已经经过深入研究的蛋白,可以推断这个未知蛋白的结构与功能的某些性质。
必须指出的是,不能够仅仅是通过比较分析这一判据来断定结论是否正确,结论还必须经过实验验证。
当我们发现两个基因或蛋白质具有惊人的相似性时,我们会认为他们之间具有一段共同的进化历程,从而我们判断他们会具有相似的生物学功能,但是,这个推断在成为结论之前必须经过实验的验证。
例如,ζ-晶状物是脊椎动物眼睛里晶状体基质的组成部分,根据序列相似性的基础,它在E.coli中的同源物是代谢酶苯醌氧化还原酶(如图7.2),不管二者的共同祖先如何,它们的功能在进化中已经改变了(Gonzalez et al.,1994)。
这就好象火车变成了铁路餐车,虽然对二者的外部结构的观察揭示了它们结构的历史,但是仅仅根据这一信息往往会得出有关其功能的错误结论。
当一个基因适应了一个新的功能时,保守位置通常也会发生一些形式上的变化,比如,当蛋白具有催化功能时,活性为点的残基相当保守,而当蛋白功能改变时,这些残基将会发生漂移。
Human-ZCr MATGQKLMRAVRVFEFGGPEVLKLRSDIAVPIPKDHQVLIKVHACGVNPVETYIRSGTYSEcoli-QOR------MATRIEFHKHGGPEVLQA-VEFTPADPAENEIQVENKAIGINFIDTYIRSGLYP. . ******. . . * …. . . * *.* ..****** *Human-ZCr RKPLLPYTPGSDVAGVIEAVGDNASAFKKGDRVFTSSTISGGYAEYALAADHTVYKLPEKEcoli-QOR-PPSLPSGLGTEAAGIVSKVGSGVKHIKAGDRVVYAQSALGAYSSVHNIIADKAAILPAA* ** *.. **.. ** . * **** . . * *. **Human-ZCr LDFKQGAAIGIPYFTAYRALIHSACVKAGESVLVHGASGGVGLAACQIARAYGLKILGTAEcoli-QOR ISFEQAAASFLKGLTVYYLLRKTYEIKPDEQFLFHAAAGGVGLIACQWAKALGAKLIGTV. * * ** . * * * .. .* * * * *.***** *** *.* * *..**Human-ZCr GTEEGQKIVLQNGAHEVFNHREVNYIDKIKKYVGEKGIDIIIEMLANVNLSKDLSLLSHGEcoli-QOR GTAQKAQSALKAGAWQVINYREEDLVERLKEITGGKKVRVVYDSVGRDTWERSLDCLQRR** . . *. ** .* * ** …. * * * . .. . . . . * * .Human-ZCrGRVIVVG-SRGTIEINPROTMAKES----SIIGVTLFSSTKEEFQQYAAALQAGMEIGWLEcoli-QOR GLMVSFGNSSGAVTGVNLGILNQKGSLYVTRPSLQGYITTREELTEASNELFSLIASGVI* .. * * *.. . . . . . .*.** . . * . . * .Human-ZCr KPVIGSQ--YPLEKVAEAHENIIHGSGATGKMILLLEcoli-QOR KVDVAEQQKYPLKDAQRAHE-ILESRATQGSSLLIP* . * *** *** *. . * .*.图7.2、最佳全局比对:对人类ζ-晶状物(Swiss-Prot Q08257)和E.coli苯醌氧化还原酶(Swiss-Prot P28304)的氨基酸序列进行比对。
这是一个由CLUSTAL W程序(Higgins et al., 1996)得到的最佳全局比对结果。
在比对下方,星号表示残基相同,打点表示这个残基是保守的。
早期的序列比对方法只应用于那些在全长范围内具有简单相似性的一些序列。
全序列比对就是对序列进行全程扫描,进行比较。
以上讨论的胰蛋白酶和ζ-晶状物之间的比较就属于全序列比对。
具有简单的球形结构域的蛋白一般可以使用全序列比对的策略,以为所有的同源序列尚未经过实质上的变化蛋白质的模块性质许多蛋白质在全程范围内并不具有相似性,但却似乎是由众多的模块结构域搭建而成。
图7.3描述了这样的一个例子,如图所示的是在血凝过程中的两种蛋白的组成结构,它们是凝血因子XII(F12)和组织型血纤蛋白溶酶原活化因子(PLAT),除了具有丝氨酸蛋白酶活性的催化结构域,这两种蛋白还具有不同数量的其它结构域单元,包括两种纤连蛋白重复,一个类似于上皮生长因子的结构域以及一个成为“kringle”域的单元。