生物信息学

Chapter 1 DNA Pairwise Sequence Alignment

| 关于   «  2. 点阵图   ::   目录   ::   4. 本地对齐  »

1.3. 全局比对

动态规划算法 :

可以构造一张隐含包含所有可能比对的图表,其形式是一个矩阵,类似于绘制点阵图时所用的矩阵。一条 序列 的残基作为行索引,另一条序列的残基作为列索引。矩阵中从左上到右下的任意一条路径都对应一种比对。

动态规划如何工作 ?

动态规划的工作方式是:先构造一个二维矩阵,其两条轴分别是要比较的两条序列。残基的匹配依据某个特定的计分矩阵。得分逐行计算。首先用一条序列的第一行,去扫描另一条序列的整个长度,然后扫描第二行。匹配得分随之计算出来。第二行的扫描会考虑第一轮中已经得到的分数。这一过程反复进行,直到所有单元格的值都被填满。于是,分数沿着从左上角到右下角的对角线累积。一旦分数在矩阵中累积完毕,下一步就是找出代表最优比对的那条路径。这要 从矩阵的右下角开始,沿相反顺序回溯穿过矩阵,直至矩阵左上角的原点。最佳匹配路径就是总得分最大的那条路径。如果两条或更多路径达到相同的最高得分,就任意选取其中一条作为最佳比对。路径在某个位置也可以水平或垂直移动,这对应于为两条序列之一引入一个空位,即一次插入或删除。

比对图的动态规划递推式 :

S(i,j):从 (0,0) 到 (i,j) 的最长路径长度
S(i,j)= max [ s(i-1,j) + (i-1,j) 与 (i,j) 之间“垂直”边的权重
s(i,j-1) + (i,j-1) 与 (i,j) 之间“水平”边的权重
s(i-1,j-1) + (i-1,j-1) 与 (i,j) 之间“对角”边的权重
]
Complete binary tree node numbering

成对序列比对的总体目标 :

是找出两条序列的最佳配对,使残基之间具有最大的一致性。为了实现这一目标,需要将一条序列相对于另一条平移,找到匹配数最多的位置。常用的两种不同比对策略是:全局比对与局部比对。

全局比对的动态规划 :

全局比对问题 :在给定计分矩阵下,求两个字符串得分最高的比对。

输入 :两个字符串以及计分矩阵 score。

输出 :两个字符串的一种比对,其比对得分在二者的所有比对中最大。

要解决全局比对问题 :

我们仍需在更新边权以反映计分矩阵中的值之后,找出比对 图 中的一条最长 路径 。
回忆一下,“删除”对应垂直边,“插入”对应水平边,“匹配/错配”对应对角边,
我们便得到下列关于 s(i,j)(从 (0,0) 到 (i,j) 的最长路径长度)的递推式:
S(i,j)= max [ s(i-1,j) + (i-1,j) 与 (i,j) 之间“垂直”边的权重
s(i,j-1) + (i,j-1) 与 (i,j) 之间“水平”边的权重
s(i-1,j-1) + (i-1,j-1) 与 (i,j) 之间“对角”边的权重
]

Needleman–Wunsch 算法更适用于 :

NEEDLEMAN-WUNSCH 算法包含三个步骤 :

1. 初始化得分矩阵与回溯矩阵 :

在初始化过程中,得分矩阵和回溯矩阵的第一行与第一列被初始化。

2. 计算得分并填写得分矩阵与回溯矩阵 :

下一步是迭代地求出结果矩阵中所有元素的得分值。

3. 从回溯矩阵推断比对 :

我们有两个二维矩阵:得分矩阵和回溯矩阵。

Settings

Proficient Saving... Error Saving
Server Error
Resubmit

1.3.1. 得分矩阵

Settings

Proficient Saving... Error Saving
Server Error
Resubmit

1.3.2. 回溯

Settings

Proficient Saving... Error Saving
Server Error
Resubmit

1.3.3. 练习

Settings

Proficient Saving... Error Saving
Server Error
Resubmit

   «  2. 点阵图   ::   目录   ::   4. 本地对齐  »

关闭窗口