第五章计算机在分子生物学中的应用课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《第五章计算机在分子生物学中的应用课件.ppt》由用户(三亚风情)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第五 计算机 分子生物学 中的 应用 课件
- 资源描述:
-
1、第五章第五章 计算机在分子生物学中的应用计算机在分子生物学中的应用DNA双链模型51 计算机在分子生物学中应用的计算机在分子生物学中应用的简介简介 分子生物学研究的对象往往是大规模的实验数据,利用手工计算来处理这些数据显然是力不从心.例如越来越多的物种的基因组将基本上完全地测定。那种倾毕生精力研究一个基因、一条代谢途径、一种生理周期的时代已经过去.genbank数据增长示意图 那种倾毕生精力研究一个基因、一条代谢途径、一种生理周期的时代已经过去。人们正在阐明细胞内的全部互相耦合的调控网络和代谢网络,细胞间的全部信号传导过程,从受精卵到成体的全部生理和病理的基因表达的变化等等。这一切都超出手工分
2、析的可能性,数据的产生、搜集和分析,都必须依靠计算机和网络,都必须发展数据库、算法和程序。计算机科学的发展及其在生物学领域的应用,已经成为生物学发展和进步过程中不可替代的重要力量。计算机在分子生物学发展中的作用是无可替代的。在分子生物学中,DNA、RNA和蛋白质都是表现为特定的序列。不同生物的DNA或蛋白质的相似性是多方面的,可能是核酸或氨基酸序列的相似性,也有可能是结构的相似性。生物功能分子的序列测序与功能预测是从序列中发现基因的两个层次。测序的大致步骤如下:v取DNA目标序列;v查找开放阅读框架(ORF)并将目标序列翻译成蛋白质序列;v据库中进行序列搜索;v进行目标序列与搜索得到的相似序列
3、的整体列线(global alignment);v查找基因家族v查找目标序列中的特定模序 预测目标序列结构 获取相关蛋白质的功能信息把目标序列输入“提醒”服务器 521序列比较中的计算机技术序列比较中的计算机技术 从生物学的角度而言,一个普遍的规律是:从生物学的角度而言,一个普遍的规律是:序列决定结构,结构决定功能。序列的比序列决定结构,结构决定功能。序列的比较一般不考虑空间结构或功能的相似性。较一般不考虑空间结构或功能的相似性。研究序列的相似性的另一个目的是通过序研究序列的相似性的另一个目的是通过序列的相似性,判别序列间的同源性,推测列的相似性,判别序列间的同源性,推测序列间的进化关系。序列
4、间的进化关系。序列比较的作用是:发现生物序列中的功序列比较的作用是:发现生物序列中的功能、结构和进化的信息,从而发现其中的能、结构和进化的信息,从而发现其中的相似性,找出序列间的共同区域,同时辨相似性,找出序列间的共同区域,同时辨别序列之间的差异别序列之间的差异 5211、序列的相似性台戏在计算机内部,不管是DNA、RNA还是蛋白质,都是用特定的字符集来表示的。对于一种未知功能的生物分子,则可以通过将它的序列与已知功能的分子的序列进行比较来推断。序列的相似性可以用定性的方法来描述,也可以用定量的方法表示。在讨论到序列相似性的关系时,经常会遇到同源(homology)和相似(similarity
5、)两个概念。所谓同源序列,简单地说,是指从某一共同祖先经趋异进化而形成的不同序列。相似性(similarity)和同源性(homology)是两个完全不同的概念。相似性概念的含义比较广泛,除了上面提到的两个序列之间相同碱基或残基所占比例外,在蛋白质序列比对中,有时也指两个残基是否具有相似的特性,如侧链基团的大小、电荷性、亲疏水性等。序列比较的基本操作是比对(align),它是一种关于序列相似性的定性描述,反映的主要是在什么部位两条序列相似或差异。如果一个比对方法能够揭示两条序列的最大相似程度或根本差异,就称这个比对是最优比对。1.字符表和序列:在计算机中处理生物功能分子的序列比对时,将其序列抽
6、象为字符串,这些字符串从一个特定的字符集合中抽取,这个字符集合称为:字符表。如教材中的表5.1和表5.2v在分子生物学研究的一些场合,常常要用到子序列,如:分析功能基因或是保守序列,重复序列。生物序列中的子序列在形式上看起来同计算机数据结构中的子串的概念很相近,但实际上子序列和子串还是有些不同的:子序列的范围包含了子串,所有的子串都是子序列,但子序列不一定是子串。子序列可以通过对序列进行选择,删除等操作或取。例如:基因片段1的序列为:ATTTTGCCCTTA,基因片段2的序列是:AGCT,基因片段3的序列是:TTGC。则片段2是片段1的子序列,但2不是1的子串,片段3是片段1的子串。如果有两个
7、生物分子序列分别为t和s,则当t为s的子串时,称s是t 的超串。如果t是s 的子串,也称t是s的连续子序列。生物功能分子中的序列比对根据比较的范围不同可以分成全局比较和局部比较两种。全局比较指的是比较两条完整的序列,而局部比较指的是找出最大相似的子序列。对于两条序列的比对,根据不同的应用场合,常常将序列比较分成以下几种基本操作:(1)判断一个序列是不是另一个序列的子序列;(2)寻找两个序列中的最大相似子序列;(3)寻找两个相似序列中的细微差别;(4)判断一个序列的特定部份(如前缀或后缀)与另一个序列的特定部份是否相同。其中,(1)和(3)是全局比较,(2)和(4)是局部比较。2编辑距离对于两条
8、DNA序列,有时很难看出它们有相似的地方,但是只要对其中的一条序列进行了一些简单的操作,就会发现它们之间有很多的相似之处。例如,有以下两个英文单词“tomorrow”和“sorrow”,我们可以很清楚的看到,只要将sorry错移3个位置,并对起来,就可以发现它们的相似性。tomorrow tomorrowsorrow -sorrow移位前 移位后对于生物序列,有两种方法可以用来定量的表示两条序列的相似程度:一种方法是利用相似度函数来说明,相似度越大,说明两条序列相似的程度越大;另一种方法就是利用两个序列间的距离来说明,距离越大,说明两个序列的相似程序越小。一般说来,相似度较为灵活,所以应用的较
9、多两个序列间的距离,可以用海明距离表示。但对于不同长度的序列用海明距离表示起来不是很精确。而且在实际的实验中,一些生物功能分子如DNA往往会发生像删除或插入一个碱基这样的错误,这时如果用海明距离来表示时,就会产生较大的误差。为了克服海明距离的缺陷,引入了编辑距离的概念,所谓编辑距离(edit distance),指的是:一个字符串变到另一个字符串时插入、删除和替换的最少的字符个数。利用编辑距离来表示两个序列的比对时,一般说来有如下的字符编辑操作:设有两个序列s和t,用-代表空位(或空缺,space)则有如下的操作:vMatch(a,a)-字符匹配;vDelete(a,-)-从s序列中删除一个字
10、符或在t序列中插入一个空位;vReplace(a,b)-以t中的字符b替换s中的字符a,ab;vInsert(-,b)-在s序列中插入空位字符,或在t序列中删除一个字符b。进行序列比较最简单的方法就是利用点标法(Fitch,1969)来实现。这种比较方法的原理是:将两条待比较的序列分别放在二维作标的X轴上(序列的方向是自左向右)和Y轴上(序列的方向是自下而上)。当对应的行与列的字符匹配时,则在作标轴上给出相应的记号,逐个比较所有的字符对,最终形成若干个匹配子串。如下所示:如有两个序列s,t,序列分别为:s:ATCG t:ATGC4.序列比对的数学方法(1)打分矩阵打分矩阵被广泛的用于评价序列比
11、对的质量,通常采用得分(+)、无分(0)和罚分(-)来进行综合的评价。可以定义一个打分函数,用它来表示在序列比对中不同类型的编辑操作所需要的代价。假定有一字符表,字符a,b满足:a,b;则有如下定义:分别对应于得分、无分和失分的情况。1b)(-,(a,-)b)(a 0b)(a,1a)(a,在两条序列s和 t进行比对时的得分等于将s转化为t所用的编辑操作的得分总和;它们间的最优比对是可能的比对中得分最高的一个比对;s和t的真实的编辑距离应当是在打分函数值最大时的距离。这样,进行序列比对的目的就是寻找一个打分函数值最大的比对。(2)核酸打分矩阵与蛋白质打分矩阵:核酸与蛋白质都是常见的生物功能分子,
12、在分子生物学研究中,经常遇到要对它们的序列进行比对的场合。前面所说的打分矩阵方法过于简单,不能考虑到字符替换后实际的生物意义。特别对于蛋白质序列,有些氨基酸的取代是很容易产生而且不会对蛋白质的特性造成太大的影响。也就是说,不同情况下的替代是不等效的。所以,为了区分不同情况下替代对生物功能分子所起的作用,人们提出了核酸与蛋白质的打分矩阵。核酸打分矩阵(i)等价矩阵给出了一种最简单的核酸打分矩阵(等价矩阵),它的设计的原理是,只有相同核苷酸匹配的情况下打分为“1”,其它的情况下,打分均为“0”。这种矩阵过于简单,在实际的应用中很少用到。ATGCA1000T0100G0010C0001 核酸的等价矩
13、阵(ii)转换-颠换矩阵 众所周知,核酸的碱基可以分成两大类:一类是嘌呤,一类是嘧啶。嘌呤的碱基有两个环状结构,而嘧啶的碱基只有一个环。根据这个特点,如果DNA碱基的变化保持环数不变,则称为转换(transition),如G变成A,如果环数发生变化,则称为颠换(transversion),如A转成C。根据这个特性,当两个碱基的替换发生颠换时,它的打分是-5分;当发生转换时,它的打分是-1分;发生匹配时为1分。从而,也可以得到一个矩阵,通常称它为转换-颠换矩阵。ATGCA1-5-5-1T-51-1-5G-5-11-5C-1-5-51转移-颠换矩阵(iii)BLAST矩阵 BLAST(basic
14、local alignment search tool)是一种基本的局部对位排列搜索工具,这里也提供了一个相似性记分矩阵。这个矩阵也相对简单,如果等比较的两个核酸序列是相同的,则打5分,反之,得分为-4分。ATGCA5-4-4-4T-45-4-4G-4-45-4C-4-4-45 BLAST矩阵 2)蛋白质打分矩阵(i)等价矩阵:假设蛋白质的字符表如教材上表5.1所示,则可以构建如下的等价矩阵(如教材上表5-6所示)。它的规则是当组成蛋白质的两种氨基酸相匹配时,打分为“1”,反之,均为“0”。(ii)疏水矩阵v蛋白质由于它的氨基酸残基上的电荷不同,可以分成极性氨基酸、带电氨基酸和疏水氨基酸三大类
15、。所谓的疏水指的是氨基酸与水的亲和力的很小,这主要是因为疏水性强的氨基酸中的原子间仅靠非极性共价键相连,所以,这类氨基酸分子缺少与水分子共同作用的基础。而与疏水性氨基酸相对应的是亲水性氨基酸,这些氨基酸中的原子存在极性的共价键,从而可以与水互相溶解。根据氨基酸的亲水或疏水,也可以生成一个矩阵,称为疏水矩阵,它的设计思想是:如果一个氨基酸残基取代另一个氨基酸残基后,疏水性没有发生太大的变化,就打分高些;反之,如果替换后,疏水性发生了较大的变化,打分就低些。如下图所示:蛋白质疏水矩阵示意图(iii)GCM矩阵生命是不断进化的,在研究分子水平的进化时,常常用到GCM矩阵,它可以方便地描述分子的进化距
16、离,并可以用来绘制进化树。但在蛋白质比对中较少直接用到。GCM矩阵的设计思想是:计算一个氨基酸残基转变成另一个氨基酸残基所需的密码子变化的次数,将变化的次数作为对应矩阵的元素的值。如果一个氨基酸的残基只要有一个碱基发生变化,那么这两个氨基酸的替换代价即为1;如果是发生了两个碱基的变化,则为2,其它依此类推。iv)Dayhoff突变数据矩阵(PAM矩阵)一个PAM的进化距离定义为每100个氨基酸中发生一个点突变的概率。在这个矩阵中,大于0的值表明发生的突变的可能性较大,等于0是中性的(随机突变),小于0的则表示发生突变的可能性较小。一个PAM就是一个进化的变异单位,即1%的氨基酸发生改变,但实际
17、上并不可能说经过100次变化,每个氨基酸都会发生变化。PAM有一系列的的替换矩阵,每个矩阵用于特定的进化距离的序列。但是一般说来,只有当置换速率通过至少具有85%一致性的序列对位排列才能获取。PAM250矩阵V)模块替换矩阵(BLOSUM矩阵)Henikoff(1992)首先从BLOCKS数据库的对位排序列块中导出了一级置换矩阵,称为BLOSUM矩阵。它是从蛋白质序列块(短序列)比对而推导出来的,它用关系较远的序列来获取矩阵元素;而低阶BLOSUM矩阵更多是用来比较亲缘较远的序列。BLOSUM62矩阵图 小结:(I)基于“等价矩阵”的记分 这种记分方法,只考虑序列是否匹配,匹配的位点记正分(通
18、常为1),非匹配的位点记0分。这种方法的优点是:简单明了,适用于高度相似性序列;这种方法的缺点是:没有考虑非匹配位点间的不等价问题,在对相似性较低的序列进行对位排列时,效果尤差。(II)基于“化学相似性”的记分方式 该方法是对一致性记分方法的局部改进。例如,Mclachlan(1972)和Feng et al。(1984)结合氨基酸的性质(如极性、电荷、大小和结构特征),对不同的氨基酸进行了加权。这种方法的优点是考虑了氨基酸和蛋白质的结构与性质;缺点是并非所有的蛋白质的结构与功能的改变都可以用简单的记分描述。(III)基于“遗传密码”的记分 该方法考虑到当一个氨基酸转换成另一个氨基酸时,在基因
19、组水平上碱基变化的最小数目。这种方法的优点是它充分考虑到了在分子水平上的变化,具有一定的分子生物学基础。但是,它忽略了随机因素,例如:碱基变化的数目并不是氨基酸序列间相似性的惟一决定因素。(IV)基于“实验观察”的记分 这种方法考虑了对位排序中所实际观察的频率,从而更有助于解释序列间的进化关系。Dayhoff和BLOSUM矩阵就属于这样的矩阵。Dayhoff矩阵基于进化的突变模型基于蛋白质家族进化过程中从一个共同祖先分化的蛋白质的首次变化的。而BLOSUM矩阵忽略近端和远端的关系,这称为蛋白质进化的星状模型。Dayhoff对相关序列中所有氨基酸位置进行计分,而BLOSUM矩阵则是基于区块中置换
20、和保守位置。因而,Dayhoff模型可用于寻找蛋白质的进化起源,而BLOSUM模型用于发现蛋白质的保守域。计算机在生物序列比对处理中起到的作用计算机在生物序列比对中起到的作用是显著的:1.比对算法是比效率高低的重要基础全局比对和局部比对各有其相应的算法2.数据存储的形式和数据压缩三角形矩阵,稀疏矩阵还有序列的压缩算法可以节省空间,降低大量数据存放时要占用的大量空间5212 序列的两两比对序列的两两比对在生物学中,对各种生物功能分子的序列进行分析是在生物学中,对各种生物功能分子的序列进行分析是一件非常基本的工作。在遗传物质长期的演化过程一件非常基本的工作。在遗传物质长期的演化过程中,一些序列在进
21、化的过程中不免发生一些变化。中,一些序列在进化的过程中不免发生一些变化。在进行比对时,这些序列就不能进行精确的匹配,在进行比对时,这些序列就不能进行精确的匹配,但是他们具有一定的相似性。我们应该如何判定序但是他们具有一定的相似性。我们应该如何判定序列之间的这种相似程度?对于这种情况,生物学家列之间的这种相似程度?对于这种情况,生物学家提出了一种用来评定序列相似性的方法,称为记分提出了一种用来评定序列相似性的方法,称为记分函数的方法函数的方法。1、两两比对的基本算法 进行序列的两两比对最直接的方法就是先生成两条待比较序列的所有可能比对,然后分别计算得分函数的值,在这些结果中寻找一个值最大的比对(
22、也就是代价最小的比对)。生物序列比对算法实际上常常用到的算法是著名的N-W算法与S-W算法,它们都是动态规划算法。其中,N-W算法常用于序列的全局比对,S-W算法常用于序列的局部比对。(1)N-W算法算法 1970年,Needleman和Wunsch提出了著名的Needleman-Wunsch算法,简称为:N-W算法。Needleman-Wunsch算法是一种整体联配(global alignment)算法,最佳联配(两条最佳联配(两条蛋白质序列具有最多匹配残基)蛋白质序列具有最多匹配残基)中包括了全部的最短匹配序列。这一算法是为氨基酸序为氨基酸序列发展列发展的,算法最初寻求的是使两条序列间的
23、距离最小。尽管这类距离的元素是以一种特定的方式定义的,但该算法的良好特性在于它确定了最短距离。这是一个动态规划(dynamic programming)的方法。该算法可以用代数形式加以描述。设有两个序列S和T,Si和Tj(0iLength(S),0jlength(j),length表示求序列的长度)都属于某个字符集,这两个序列间的距离可以用记分函数(S,T)表示。通过评价序列S中的前i个位置和序列T中的前j位置的距离(Si,Tj),递归得到距离(S,T)。由于S和T的长度为m=Length(S)和n=Length(T),所以它的期望距离是(Sm,Tn)。在单元(i,j)内,到达该单元距离增加的
24、三种可能事件是:v从单元(i-1,j)向(i,j)方向垂直移动,相当于在T序列中插入一个空位使相似序列延伸,即:T序列由S序列中的缺失产生,这一事件的权重记作W_(Si);v从单元(i,j-1)向(i,j)方向水平移动,相当于在序列T中增加一个空位使得序列延伸,即:T序列由Tj插入到S序列产生,这一事件权重记作W+(Tj);v从单元(i-1,j-1)向(i,j)对角线移动,相当于增加Si与Tj使得相似序列延伸,即:S序列的Si由T序列的Tj取代所得。这一事件的权重记为W_(Si,Tj);所以,单元(i,j)的距离可以看作是三个相邻单元的距离和相应的权重的和的最小者。)(TW)T,(S)T,W_
25、(S)T,(S)W_(S)T,(Sm)T,(S1-j1-ji1-j1-i1-j1-i1-ij1-ijiax初始条件为:(S0,T0)=0jkjkTWTS10)(),(ikikSWTS10)(),(S-W例题 将待比较的两条序列放在矩阵的两个维上,并按照公式对矩阵进行初始化打分。第一行分别表示S序列的前缀空位与T序列的前面连续j个字符组成的前缀的比对得分;第一列则表示T序列的前缀空位与S序列的前面连续i个字符组成的前缀的比对得分,如下图示:t s A C A C A C T A 0-1-2-3-4-5-6-7-8A-1 G-2 C-3 A-4 C-5 A-6 C-7 A-8 待比较序列在这里规定
展开阅读全文