[发明专利]一种基于字典学习的氨基酸序列特征提取方法

说明书

本发明公开了一种基于字典学习的氨基酸序列特征提取方法。本发明首先给定一种AAindex，每条氨基酸肽序列数据都可用P＝R₁R₂......R_L表示,通过给定的AAindex，可将上述序列编码成P_E＝E₁E₂......E_L。其中E_i为氨基酸残基的R_i的属性值；接着，对肽序列使用n‑gram模型分解成单词并分割成不同长度的片段，之后对同一长度的片段聚类生成字典；最后将样本基于欧氏距离的词频学习字典的特征构成样本特征矩阵。本发明能够有效的对氨基酸序列进行特征提取，实现降维，得到氨基酸中潜在的规律，特征选择的结果对样本类别判定具有重要意义，有效提高样本类别判定的准确率。

技术领域

本发明属于数据挖掘技术领域，涉及一种基于字典学习的氨基酸序列特征提取方法。

背景技术

随着人类基因组计划的顺利进展，越来越多的蛋白质被测定出来，而通过实验确定其结构与功能的蛋白质则相对较少，且费时、费力、费财，实验中可能还会遇到一些目前无法解决的困难，因此探索利用理论计算方法来研究蛋白质结构和功能具有重要意义。如何从一条氨基酸序列提取它的有用信息，并用适当的数学方法来描述或表示这些信息，使之能正确反映序列与结构或功能之间的关系，对于蛋白质分类研究是至关重要的，也是决定分类质量的关键。目前的氨基酸序列的特征提取方法主要分为两类：一类为仅仅基于氨基酸组成和位置的方法；另一类为基于氨基酸物理化学性质的方法。氨基酸是组成蛋白质的基本单位，一条蛋白质包含的基本信息是20种氨基酸的种类和排列顺序，因此基于氨基酸组成和位置的特征提取算法是最简单、最直观的方法，主要有氨基酸组成(AAC)、熵密度(EDP)、n阶耦联组成(n-OCC)和完全信息集(CIS)。氨基酸的侧链决定了氨基酸的种类、20种氨基酸侧链在形状、大小、负电性、水性以及酸碱性等方面都存在差异，正是这20种氨基酸的差异，使各种不同组合的氨基酸序列形成各种不同的蛋白质结构，并适应各类环境，完成其特定的生理功能。蛋白质的生物学活性和理化性质主要决定其空间结构的完整，因此仅仅知道蛋白质的氨基酸组成和它们的排列顺序并不能完全了解蛋白质的结构，需要考虑氨基酸的性质。目前，基于氨基酸物理化学性质特征提取方法主要有自相关函数、伪氨基酸组成(PseAA)、准序列次序作用和疏水模式组成(HP)。

上述所介绍的方法会丢失氨基酸的顺序信息以及其间的相互作用或计算量非常大等局限性，对于进一步预测蛋白质结构类起到的效果并不是十分明显。

发明内容

本发明的目的是针对现有预测蛋白质结构类的需求，提出一种基于字典学习的氨基酸序列特征提取方法。该方法使用AAindex作为初始特征来分析序列，AAindex是一个包含AAindex1、AAindex2、AAindex3的数字指标数据库，它代表氨基酸和对氨基酸的不同物理化学和生物化学性质。有利于更好地发现理化性质之间的隐藏模式，从而得到潜在的规律，并可以有效的降低特征维数。

本发明的具体步骤如下：

步骤1、输入数据，其中输入数据为蛋白质序列；

步骤2、每条肽序列数据都用P＝R₁R₂......R_L表示,给定一个氨基酸指数AAindex，将上述序列编码成P_E＝E₁E₂......E_L；其中E_i为氨基酸残基的R_i的属性值。

步骤3、构造字典；

在字典构造阶段，分为以下3步：

(1)、肽序列编码：利用n-gram模型将肽序列分解成单词，建立编码序列；

(2)、编码序列分割：将编码序列分成不同长度的片段，如长度为2的片段和长度为3的片段。