1、实用生物信息技术课程第 2 次作业1人、小鼠、大鼠血红蛋白及其编码基因序列分析1. 认真阅读 NCBI 书架（Bookshelf） 、PDB 分子月报（Molecule of the Month） 、UniProt蛋白质精选（Protein Spotlight）以及维基百科（Wikipedia）等网站中有关血红蛋白的介绍，了解血红蛋白的生理功能、空间结构、亚基组成等基本知识。血红蛋白是高等生物体内负责运载氧的一种蛋白质，具有四聚体蛋白质四级空间结构。人体内的血红蛋白由四个亚基构成，分别为两个 亚基和两个 亚基，每个亚基由一条肽链和一个血红素分子构成。血红蛋白的每个亚基由一条肽链和一个血红素分子
2、构成，肽链在生理条件下会盘绕折叠成球形，把血红素分子抱在里面，这条肽链盘绕成的球形结构又被称为珠蛋白。血红素分子是一个具有卟啉结构的小分子，在卟啉分子中心，由卟啉中四个吡咯环上的氮原子与一个亚铁离子配位结合，珠蛋白肽链中第 8 位的一个组氨酸残基中的吲哚侧链上的氮原子从卟啉分子平面的上方与亚铁离子配位结合，当血红蛋白不与氧结合的时候，有一个水分子从卟啉环下方与亚铁离子配位结合，而当血红蛋白载氧的时候，就由氧分子顶替水的位置。2. 查阅 ENSEMBL 基因组数据库中已经或正在进行基因组测序的物种树，了解人、小鼠、大鼠三个物种之间演化关系；检索物种分歧时间数据库 TimeTree，了解人和小鼠、
3、小鼠和大鼠之间的分歧时间。3. 从 UniProt 数据库中检索并提取人、小鼠、大鼠血红蛋白 alpha 亚基蛋白质序列，利用WebLab 或 JEMBOSS 软件包中的序列比对程序 Needle，选择默认计分矩阵BLOSUM62 和默认空位罚分进行序列比对，将比对结果填入表 1。分析比对结果，说明得到上述结果的原因和进一步分析思路。表 1 人、小鼠、大鼠血红蛋白 alpha 亚基蛋白质序列比对结果物种Species登录号Accession得分Score相同氨基酸Identity相同和相似氨基酸Similarity空位Gaps人 /小鼠 NP_000549/NP_032244653 123/1
4、42(86.6%) 131/142(92.3%) 0人 /大鼠 NP_000549/ 587 111/142(78.2%) 120/142(84.5%) 0实用生物信息技术课程第 2 次作业2NP_037228小鼠 /大鼠 NP_032244/NP_037228632 120/142(84.5%) 127/142(89.4%) 04. 检索 RefSeq 数据库中人、小鼠、大鼠血红蛋白 alpha 亚基编码基因的 mRNA 序列，下载其编码区序列和所编码的氨基酸序列，进行序列比对，将结果填入下面表 2 中。表 2 人、小鼠、大鼠血红蛋白编码基因序列比对结果物种 基因名 编码区序列登录号编码区序
5、列差异位点氨基酸序列登录号氨基酸序列差异位点人 HBA1/HBA2 NM_000558/NM_0005170 NP_000549/NP_0005080小鼠 Hba-a1/Hba-a2 NM_008218/NM_0010839551 NP_032244/NP_0010774240大鼠 Hba1/Hba2 NM_013096/NM_0010077225 NP_037228/NP_00100772315. 对上述人 HBA1、小鼠 Hba1-a1、大鼠 Hba1 编码区序列进行序列比对（注意选择恰当的计分矩阵和 r） ，将比对结果填入表 3。表 3 人、小鼠、大鼠血红蛋白编因编码区序列比对结果物种
6、基因名 编码区序列登录号编码区序列相同位点氨基酸序列相同位点人/小鼠 HsHBA1/MmHba-a1 NM_000558/NM_008218354/432 127/142人/大鼠 HsHBA1/RnHba1 NM_000558/NM_013096342/430 111/142小鼠/大鼠 MmHba-a1/RnHba1 NM_008218/NM_013096386/430 127/1426. 分析上述人、小鼠、大鼠血红蛋白基因区序列比对结果，阅读血红蛋白基因家族演化论文（Hardison R, 2012） ，说明人、小鼠和大鼠 alpha 血红蛋白基因家族基因结构和演化。7. 将上述数据库检索、
7、序列比对、结果分析的方法、思路、策略用于和你研究课题相关或你熟悉的蛋白质及其编码序列，分析比对结果，说明原因。表 4 人、小鼠、大鼠胃蛋白酶(PGA5)蛋白质序列比对结果得分 相同氨基酸总数（百分比） 相同及相似氨基酸总数（百分比）人 P0DJD9/小鼠 Q9D1061146 219/390 (56.2%) 282/390 (72.3%) 人 P0DJD9/大鼠 Q9JJX21213 222/422 (52.6%) 245/422 (58.1%) 小鼠 Q9D106/大鼠 Q9JJX22742 364/387 (94.1%) 368/387 (95.1%)表 5 人、小鼠、大鼠胃蛋白酶(PGA
8、5)mRNA 编码区序列比对结果得分 相同核苷酸总数（百分比） 空位数（百分比）人 NM_014224.2/小鼠 NM_021453.42190.0 748/1200 (62.3%) 209/1200 (17.4%)人 NM_014224.2 2193.5 748/1179 (63.4%) 167/1179 (14.2%)实用生物信息技术课程第 2 次作业3/大鼠 NM_021753.1小鼠 NM_021453.4/大鼠 NM_021753.14876.0 1028/1094 (94.0%) 0/1094 ( 0.0%) 通过对人、小鼠、大鼠胃蛋白酶(PGA5)蛋白质序列和 mRNA 序列比对，结果如表 3 和表4。显示，小鼠和大鼠的胃蛋白酶蛋白质序列及其编码的 mRNA 序列都具有更高的同源性，可能与二者的膳食结构有关系，而人是杂食性的且结构复杂，所以用于初级消化的蛋白酶差距会与小鼠、大鼠相差很大。进一步分析，通过多三者的 mRNA 序列比对，结果也进一步的验证蛋白质序列比对的结果。参考文献：Hardison RC. Evolution of hemoglobin and its genes. Cold Spring Harb Perspect Med. 2012 Dec 1;2(12):a011627.
