多重序列比对的数学模型与算法自美国提出组织的人类基因组计划(Human Genome Proreet)简称为HGP以来,美国每年拔出相当大的经费支持,日本、法国、英国、德国等纷纷响应,它们的工作使新的交叉学科生物信息论得以诞生和发展,生物信息论是用数理和信息科学的观点、理论和方法去研究生命现象,组织和分析呈指数增长的生物学数据。生物信息学是一门综合学科,是计算机科学、数学、物理、生物学的结合。生物信息学的基础是各种数据库的建立和分析工具的发展。目前,生物学数据库已达500个以上,共有四大类:基因组数据库,核酸和蛋白质一级结构数据库、生物大分子三维空间结构数据库及其以她们为基础构建的二级数据库。生物信息学主要研究基因组测序及其信息分析、生物大分子的结构与功能预测及其模拟和药物设计、大规模基因表达数据的分析与基因芯片设计,以及基因与蛋白质相互作用网络等四方面的问题。
8 v4 O/ p$ d) F: ]6 E多重序列比对是计算分子生物学中最重要的运算。多重序列比对的基本问题就是找出适当安排删减与插入尽量少的空格,使得两个序列达到最大程度的一致的方案。比如给出下列三个序列:
$ W% G; {3 q' _( f- ~8 j AC_G 7 j. X* L# A+ R9 A4 ?6 \) _/ d6 ]
AGTCC (1)
2 Z0 q" r1 Z' x* m$ I ACT
+ v" Z; |0 D2 e G. D, d' s/ ~+ P' }# Z$ x* w0 Q# q% X+ o& S e
) e% N0 y0 `, I$ i* L/ @8 L1 Q
我们适当安排删减与插入空格得到:
7 ~" G& Y9 Q) @& o ]" S: k ACG___
. Q# H2 t. E- b. B4 q/ A A_GTCC (2). u2 ~! {* Z) U4 x' u5 ~/ \
AC_T__9 v6 E% ]/ x( A7 p# C+ W
3 i6 ^) c ]+ _2 _1 J(2)就是多重序列的一个比对。
! _6 ^ z7 \; ^局部分段比对是其中更为常见的运算。上世纪80年代,Smith-Waterman提出了两个序列的局部比对的明确的模型。1998—1999年,相继出现利用k-tuple的快速容错分段比对搜索法。2002年开始出现对完整基因组及其异常基因的比较研究以及多重序列比对问题的研究,2003年刘军Mayetri Gupta和刘军得到Motif的搜索算法。
" r- `: \# X+ Q2 I5 ?6 g人类基因组计划后,目前已经进入后基因时代,主要就是对人类基因组计划实施得到的基本数据库进行信息分析、加工和利用,提取有用信息,用来研究生命现象中的重大问题。多重序列比对问题是生物信息学的基本问题,多重序列比对技术也是生物信息学的基本工具,有着十分广泛的应用,比如基因是否为同一个家族,癌症患者的基因与正常时的基因比对分析等等。因此,请您们就基因的多重序列比对,设计合理的衡量比对好坏的定量描述模型,建立多重序列比对的基本问题的数学模型,并设计一种求解的算法。最后就附录一中的12个序列,请您们利用你们得到的模型与算法,给出使序列有最大相似程度的比对。
) t1 q' s3 t* X6 l! e! G8 J附录一:
! I, P0 U, S0 cCATTTCTTTTTAGGGATTTTAAAAGTTGTCTTTTCTT
8 f% A& q- c- j8 p, ]; y& {* l5 DCATTTCTTTTTAAGGTTTTAAAAATTGTCTTTTTT6 X- k0 O3 \3 h
CATTTCTTTTTAAGGGTTTTAAAAATTGTCTTTTCTT8 E2 ]& p" L9 h% ], c; X
CATTTTTTCTTAAGTGTTTTGGTATTTATCTTTTTCTT6 h4 s3 v5 ~' x! \+ y0 I
CATTTTTGCTTATGTATTTATAGTGGGTTGTCTTTTTGACTT& x) F8 ^# M; A# W1 K, B
CATTTCTTTTGAAGTGATTTGAGATTTATCTTTTTCTT
, Q. s3 C! w8 d- x' VCATTTCTTTTTAAGGGTTTTAAAAATTGTCTTTTCTT, D6 `8 b4 O, L0 b, y0 W
CATTTCTTTTTATGTTGAGATATTTGTCTGTTTTCTT/ p g! c9 {; `6 i( a' o D+ z6 h
CATTTTTACTATGTGTTGATTGTGGATTGTCTTTTCTT; ~/ V* |2 a6 Y9 e2 p8 d
CATTTCTTTTATTGAGTGAAGAAGAGATTTTGTCTTGTTTTGAT
8 ?2 |* w# u0 l( |! j# W5 rCATTTTTCTTAGTGTTTTGGTATTTATCTTTTTCTT6 ~ u0 a- a4 n: I( E% d% K
CATTTCTTTTAAGGGTTTTAAAAATTGTCTTTTCTT |