使用最长公共子序列算法进行序列比对

介绍

在分子生物学中，DNA 和蛋白质可以表示为字母序列。

DNA 序列由 A、T、G、C 组成，代表核苷碱基(nucleobases) 腺嘌呤、胸腺嘧啶、鸟嘌呤和胞嘧啶。

蛋白质由 20 个不同的字母组成，表示 20 种不同的氨基酸。

比较来自同一生物体或来自不同生物体的两个序列，称为 序列比较 (Sequence comparison)，是分子生物学中的一项重要任务。

它有助于为许多生物学问题提供解决方案，例如：

1. 预测蛋白质的结构和功能
2. 推断物种的进化历史和相关性
3. 定位基因/蛋白质中的常见子序列以识别常见基序，
4. 作为 DNA 测序的基因组组装中的一个子问题

为了进行序列比较，我们首先进行 序列比对(sequence alignment)。

什么是序列比对？

序列比对是将一个序列置于另一个序列之上以识别相似字符或子串之间的对应关系的一种方式。

它可以在脱氧核糖核酸 (DNA)、核糖核酸 (RNA) 或蛋白质序列上进行。

来自不同生物体的序列可能具有不同的大小。对齐需要在序列的任意位置插入空格，以便两者具有相同的大小。

在序列的开头或结尾插入“空格”或“间隙”。

让我们考虑一个例子。 假设我们有以下两个 DNA 序列 GACGGATTAG 和 GATCGGAATAG。

这两个序列之间的一种可能比对如下所示

GA-CGGATTAG
|| |||| |||
GATCGGAATAG

请注意根据此比对的两个序列的差异。 第二个序列中的额外 T，以及从 T 到 A 的替换。在第一个序列的合适位置添加一个“空格”以对齐序列。

比对类型

根据我们正在寻找的内容，有两种类型的比对方式：

全局比对( Global alignmen )：

沿序列的整个长度将查询序列与目标序列对齐。 全局比对是跨越两个查询序列全长的全局优化过程。

局部比对( Local alignment )：

将查询序列的子串与目标序列的子串对齐，即找到两个序列之间相似度高的局部区域。

根据被比对的序列数量，也有两种比对方式：

双序列比对（ pairwise alignment ）：
涉及两个序列，查询和与之对齐的目标。

多序列比对（ multiple alignment ）：
涉及两个以上序列的比对。

序列比对的软件工具

定位查询序列与数据库序列的相似区域是生物信息学中的一项具有挑战性的任务。

BLAST 和 FASTA 等工具有助于检测生物体之间的相似区域。它涉及使用局部成对比对、详尽的启发式算法和动态规划方法（如 Smith-Waterman 算法）来检测查询序列与数据库序列的相似区域。

最长公共子序列 (LCS) 问题

检测两个或多个序列相似性的一种方法是找到它们的最长公共子序列(Longest common subsequence)。

两个序列的最长公共子序列 (LCS) 是两个序列共有的具有最大可能长度的子序列。

举例

让我们考虑一个例子。 假设我们有以下两个 DNA 序列：TAGTCACG 和 AGACTGTC。

两个序列的LCS是AGACG，可以从下面的比对中得到。

TAGTCAC-G--
 ||  || |
-AG--ACTGTC

还有其他可能的较短长度的公共/共有序列 ( common sequences )，例如 AGCG 或 AGTC。

尽管通常多个 LCS 是可能的，但对于这个特定示例只有一个 LCS，即这两个序列没有其他长度为 5 的公共子序列。

LCS 还有助于计算两个序列的相似程度：LCS 越长，相似性越高。

LCS算法

假设我们有两个长度分别为 m 和 n 的序列 S1 和 S2，其中
S1 = a₁ a₂ … a_m
S2 = b₁ b₂ … b_n。

我们将构造一个矩阵 A，其中 A_i,j 表示 a₁ a₂ … a_i 和 b₁ b₂ … b_j 的最长公共子序列的长度。

序列 S1 = TAGTCACG 和 S2 = AGACTGTC 的示例矩阵。

序列 S1 垂直写入，S2 水平写入。 将代表行的每个字母与代表列的每个字母进行比较。 我们将逐行、逐列进行。

1. 如果 a_i = b_j，那么我们找到了匹配项。 当前匹配的得分为 1，而 A_i-1,B_j-1 来自我们已经从子串 a₁ … a_i-1 和 b₁ … b