预测蛋白质之间信号流走向的方法
【专利摘要】本发明公开了一种预测蛋白质之间信号流走向的方法。该方法包括以下步骤：S1，建立蛋白质之间信号流走向的标准数据集；S2，分别将阳性数据集和阴性数据集中的上游蛋白质与下游蛋白质的第一氨基酸序列信息中的特征值进行数值运算，得到相互作用对应的第二氨基酸序列信息；S3，利用假设检验方法，计算阳性数据集和阴性数据集中的第二氨基酸序列信息的差异显著性P值；S4，将差异显著性P值小于预定阈值的序列信息作为分类特征，建立分类模型；S5，对分类模型的分类效果进行评估，如果分类模型评估满足预设标准，结束建模过程；否则，返回步骤S4，对分类模型进行校正。应用本发明的技术方案，可以不受蛋白质是否包含结构域的限制。
【专利说明】预测蛋白质之间信号流走向的方法
【技术领域】
[0001]本发明涉及生物信息【技术领域】，具体而言，涉及一种预测蛋白质之间信号流走向的方法。
【背景技术】
[0002]高通量实验技术的发展已经产出了大规模、多物种的蛋白质相互作用数据。在蛋白质相互作用网络中，通常假定相互作用是没有方向的。实际上，在信号转导、转录调控、细胞循环或者代谢途径等多种生物网络中，发生相互作用的蛋白质之间广泛存在着调控和上下游关系。然而，仅有少数蛋白质相互作用被深入研究，大部分相互作用的细节还是未知的，尤其是它们之间的信号流走向还未被揭示。
[0003]蛋白质相互作用的方向性是构成信号通路的前提，有助于从大规模的蛋白质相互作用网络中挖掘新的潜在信号通路。目前，已有一些生物信息学方法基于结构域或功能注释信息来预测蛋白质间信号流的走向。如文章Liu W，et al.Proteome-wide predictionof signal flow direction in protein interaction network based on interactingdomains.Molecular&Cellular Proteomics, 2009, 8 (9): 2063-2070 提出 了一种基于结构域相互作用发现蛋白质之间信号流走向的方法，详细步骤如下:1)以人、小鼠、大鼠、果蝇和酵母中已知方向的蛋白质相互作用作为标准数据集，根据结构域相互作用在数据集中的相对富集程度定义打分函数，发现能够提示信号流走向的结构域相互作用；2)以第一步得到的有向结构域相互作用为基础，统计蛋白质相互作用中结构域相互作用流向的综合打分，用于预测蛋白质相互作用中信号流的方向；3)采用五倍交叉验证和物种间交叉验证对该预测方法的性能进行评估。
[0004]但是，由于现有技术中的预测方法主要是基于结构域等信息预测蛋白质对之间的信号流走向，由于部分蛋白质不包含结构域或仅包含单个的结构域，使得该方法的应用范围非常有限，如大约一半的人蛋白质相互作用无法通过该方法来预测信号流走向。
【发明内容】
[0005]本发明旨在提供一种预测蛋白质之间信号流走向的方法，以解决现有技术中不包含结构域或仅包含单个的结构域的蛋白质预测信号流走向困难的技术问题。
[0006]为了实现上述目的，根据本发明的一个方面，提供了一种预测蛋白质之间信号流走向的方法，包括以下步骤:Si，从信号转导数据库中收集有向的调控关系，建立蛋白质之间信号流走向的标准数据集，从标准数据集中随机选择一半信号流走向的数据作为阳性数据集，剩余一半的信号流走向取反得到的数据作为阴性数据集；S2，分别提取阳性数据集和阴性数据集中蛋白质所包含的第一氨基酸序列信息，并分别将阳性数据集和阴性数据集中的上游蛋白质与下游蛋白质的第一氨基酸序列信息中的特征值进行数值运算，得到阳性数据集和阴性数据集中的相互作用对应的第二氨基酸序列信息；S3，利用假设检验方法，计算阳性数据集和阴性数据集中的第二氨基酸序列信息的差异显著性P值；S4，将差异显著性P值小于预定阈值的序列信息作为分类特征，建立分类模型；S5，采用交叉验证方法或独立测试集，对分类模型的分类效果进行评估，如果分类模型评估满足预设要求，则分类模型建立成功，结束建模过程；否则，返回步骤S4，对分类模型进行校正。
[0007]进一步地，有向的调控关系包括激活、抑制、磷酸化、以及非磷酸化。
[0008]进一步地，蛋白质所包含的第一氨基酸序列信息包括物理化学属性和氨基酸排列信息。
[0009]进一步地，数值运算为选自加、减、乘、除、乘方、开方和取模中的一种或多种。
[0010]进一步地，假设检验方法为选自T检验、U检验、卡方检验、F检验,秩和检验中的一种或多种。
[0011]进一步地，预定阈值Ρ〈0.05。
[0012]进一步地，分类模型为选自支持向量机、决策树、贝叶斯方法中的一种或多种。
[0013]进一步地，上游蛋白质与下游蛋白质之间发生相互作用。
[0014]进一步地，步骤SI包括:建立多个物种中蛋白质之间信号流走向的标准数据集。
[0015]应用本发明的技术方案，通过对蛋白质中氨基酸序列信息的分析，建立合适的预测模型推断蛋白质相互作用之间的信号流走向，从而可以不受蛋白质是否包含结构域的限制，不仅能够用于含有多个结构域的蛋白质相互作用的预测，还能够用于不包含结构域或仅包含单个的结构域的蛋白质相互作用的预测。由于绝大部分蛋白质的氨基酸序列是已知的，该方法的应用范围非常广泛，可以覆盖90%以上的蛋白质相互作用；而且，还可以用于标注大规模蛋白质相互作用的信号流走向，帮助阐释信号分子作用机制，辅助实验设计，节省大量的人力物力；另外，因为信号转导数据库数据丰富，可以考虑多种，甚至上百种序列信息，可靠性更高，且方法简单实用。
【专利附图】
【附图说明】
[0016]构成本申请的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中:
[0017]图1示出了本发明的预测蛋白质之间信号流走向的流程图；
[0018]图2示出了根据本发明一种实施方式的蛋白质序列三元组划分方法示意图；以及
[0019]图3示出了根据本发明一种实施方式的支持向量机分类超平面示意图。
【具体实施方式】
[0020]需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图来详细说明本发明。
[0021]针对现有技术中不包含结构域或仅包含单个的结构域的蛋白质预测信号流走向困难的技术问题，本发明提供了一种基于氨基酸序列信息的蛋白质之间信号流走向的预测方法。
[0022]根据本发明一种典型的实施方式，提供一种预测蛋白质之间信号流走向的方法。如图1所示，该方法包括以下步骤:SI，从信号转导数据库中收集有向的调控关系，建立蛋白质之间信号流走向的标准数据集，从标准数据集中随机选择一半信号流走向的数据作为阳性数据集，剩余一半的信号流走向取反作的数据为阴性数据集；S2，分别提取阳性数据集和阴性数据集中蛋白质所包含的第一氨基酸序列信息，并分别将阳性数据集和阴性数据集中的上游蛋白质与下游蛋白质的第一氨基酸序列信息中的特征值进行数值运算，得到将阳性数据集和阴性数据集中的相互作用对应的第二氨基酸序列信息；S3，利用假设检验方法，计算阳性数据集和阴性数据集中的第二氨基酸序列信息的差异显著性P值；S4，将差异显著性P值小于预定阈值的序列信息作为分类特征，建立分类模型；S5，采用交叉验证方法或独立测试集，对分类模型的分类效果进行评估，如果分类模型评估满足预设要求，则分类模型建立成功，结束建模过程；否则，返回步骤S4，对分类模型进行校正。
[0023]本发明的方法也是一种基于支持向量机模型的蛋白质之间信号流走向的预测方法，其中，信号转导数据库是指KEGG、BioCarta, STEK、NCI。取反是指调换发生相互作用的两个蛋白质的次序，即把上游的蛋白质改至下游，下游的改至上游。步骤S2具体为提取阳性数据集中蛋白质所包含的第一氨基酸序列信息，将其中上游蛋白质与下游蛋白质的第一氨基酸序列信息中的特征值进行数值运算，得到阳性数据集中的相互作用对应的第二氨基酸序列信息；提取阴性数据集中蛋白质所包含的第一氨基酸序列信息，将其中上游蛋白质与下游蛋白质的第一氨基酸序列信息中的特征值进行数值运算，得到阴性数据集中的相互作用对应的第二氨基酸序列信息。步骤S5中的预设要求可以是预测准确率达到预设目标(如80%)以上。
[0024]应用本发明的技术方案，通过对蛋白质中氨基酸序列信息的分析，建立合适的预测模型推断蛋白质相互作用之间的信号流走向，从而可以不受蛋白质是否包含结构域的限制，不仅能够用于含有多个结构域的蛋白质相互作用的预测，还能够用于不包含结构域或仅包含单个的结构域的蛋白质相互作用的预测。由于绝大部分蛋白质的氨基酸序列是已知的，该方法的应用范围非常广泛，可以覆盖90%以上的蛋白质相互作用；而且，还可以用于标注大规模蛋白质相互作用的信号流走向，帮助阐释信号分子作用机制，辅助实验设计，节省大量的人力物力；另外，因为信号转导数据库数据丰富，可以考虑多种，甚至上百种序列信息，可靠性更高，且方法简单实用。
[0025]关于步骤SI中的有向的调控关系，本领域技术人员可以根据实际需要选择，优选的，有向的调控关系包括激活、抑制、磷酸化、以及非磷酸化等，因为这些相互作用在信号网络中是有明显的方向关系的，如上游蛋白质可以激活下游蛋白质。建立蛋白质之间信号流走向的标准数据集优选为建立多个物种的蛋白质之间信号流走向，可以使预测结果更准确。
[0026]根据本发明一种典型的实施方式，蛋白质所包含的第一氨基酸序列信息包括物理化学属性和氨基酸排列信息，其中，物理化学属性包括单个氨基酸的物理化学属性。由于本发明中最优先考虑的是两两相互作用，上游蛋白质与下游蛋白质之间优选是发生相互作用的，每个相互作用中都包含两个蛋白质，因此可获得两组对应的氨基酸序列信息。
[0027]其中，根据蛋白质的一级序列，可提取其物理化学特征，如分子量、疏水性、极性等。疏水性可通过不同氨基酸对应的疏水性值进行加合平均得到。进一步，计算每个蛋白质中氨基酸的出现频率，定义为氨基酸的数目除以蛋白质的长度。根据氨基酸的属性，可以分为微小、小、非极性、极性、带电的和基本的几种。利用Pepstats程序(http://emboss,bioinformatics, nl/cg1-bin/emboss/pep stats)可一次性提取蛋白质的多个序列特征，包括疏水性、分子量、残基数和氨基酸使用频率等。[0028]氨基酸排列信息是指根据蛋白质序列中氨基酸的排列顺序和组成情况来提取序列信息。
[0029]首先，为了减少向量空间的维数，并且考虑到氨基酸同义替换的问题，对氨基酸进行分类。蛋白质相互作用主要包括静电(氢原子)相互作用和疏水的相互作用，这两类相互作用可以通过氨基酸边链的极性和体积规模体现。相应的，这两个参数可以通过密度泛函理论方法B3LYP/6-31G*和分子建模方法得到。根据边链的极性和体积规模，将20个氨基酸分成7类(如表1所示)。在相同分类中，氨基酸具有类似的特性，参与同义替换。
[0030]表1
[0031]
【权利要求】
1.一种预测蛋白质之间信号流走向的方法，其特征在于，包括以下步骤:Si，从信号转导数据库中收集有向的调控关系，建立蛋白质之间信号流走向的标准数据集，从所述标准数据集中随机选择一半信号流走向的数据作为阳性数据集，剩余一半的信号流走向取反得到的数据作为阴性数据集；S2，分别提取所述阳性数据集和所述阴性数据集中蛋白质所包含的第一氨基酸序列信息，并分别将所述阳性数据集和所述阴性数据集中的上游蛋白质与下游蛋白质的所述第一氨基酸序列信息中的特征值进行数值运算，得到所述阳性数据集和所述阴性数据集中的相互作用对应的第二氨基酸序列信息；S3，利用假设检验方法，计算所述阳性数据集和所述阴性数据集中的所述第二氨基酸序列信息的差异显著性P值；S4，将所述差异显著性P值小于预定阈值的序列信息作为分类特征，建立分类模型；S5，采用交叉验证方法或独立测试集，对所述分类模型的分类效果进行评估，如果所述分类模型评估满足预设要求，则所述分类模型建立成功，结束建模过程；否则，返回所述步骤S4，对所述分类模型进行校正。
2.根据权利要求1所述的方法，其特征在于，所述有向的调控关系包括激活、抑制、磷酸化、以及非磷酸化。
3.根据权利要求1所述的方法，其特征在于，所述蛋白质所包含的第一氨基酸序列信息包括物理化学属性和氨基酸排列信息。
4.根据权利要求1所述的方法，其特征在于，所述数值运算为选自加、减、乘、除、乘方、开方和取模中的一种或多种。
5.根据权利要求1所述的方法，其特征在于，所述假设检验方法为选自T检验、U检验、卡方检验、F检验，秩和检验中的一种或多种。
6.根据权利要求1所述的方法，其特征在于，所述预定阈值P〈0.05。
7.根据权利要求1所述的方法，其特征在于，所述分类模型为选自支持向量机、决策树、贝叶斯方法中的一种或多种。
8.根据权利要求1所述的方法，其特征在于，所述上游蛋白质与所述下游蛋白质之间发生相互作用。
9.根据权利要求1所述的方法，其特征在于，所述步骤SI进一步包括:建立多个物种中蛋白质之间信号流走向的标准数据集。
【文档编号】G06F19/18GK103577722SQ201310552227
【公开日】2014年2月12日 申请日期:2013年11月8日 优先权日:2013年11月8日
【发明者】刘伟, 谢红卫 申请人:中国人民解放军国防科学技术大学