本发明涉及用于预测蛋白质的至少一个适应度值(fitnessvalue)的方法和相关的电子系统，所述蛋白质包含氨基酸序列。本发明还涉及包括软件指令的计算机程序产品，所述程序产品被计算机执行，当被计算机执行时，所述程序产品执行这种方法。发明背景蛋白质是由至少一个氨基酸序列链组成的生物分子。蛋白质彼此相异主要在于它们的氨基酸序列，序列之间的差异被称为“突变”。蛋白质工程的最终目标之一是具有所需性质(统称为“适应度”)的肽、酶、蛋白质或氨基酸序列的设计和构建。用工程化的氨基酸置换、氨基酸的缺失或插入或氨基酸嵌段(嵌合蛋白质)修饰的氨基酸序列(即“突变体”)的构建允许评估任何特定氨基酸在适应度中的作用以及了解蛋白质结构与其适应度之间的关系。定量结构-功能/适应度关系分析的主要目的是调查和数学地描述蛋白质结构变化对其适应度的影响。突变的影响与不同氨基酸的物理化学和其它分子性质相关，并且可通过统计学分析来逼近。探索适应度景观，调查n个单点置换的所有可能的组合(排列)是一项非常困难的任务。确实，突变体的数量增加非常快(表1)。单点突变的数目突变体的数目24416664825610102412409614163841665536401.1x1012表1.针对n个突变的可能的突变体数目探索所有可能的突变体在实验上是困难的，特别是当n增加时。在实践中，在湿实验室中用单点置换产生突变体是相当容易和便宜的。对于它们中的每一个来说，适应度都可以很容易地表征。但在湿实验室中组合单点置换并不容易。产生被靶向的n个单点置换的所有可能的(2n)组合可以是非常荷求和昂贵的。大规模评估适应度是有问题的。已经开发了混合的体外和电脑模拟方法来帮助蛋白质的定向进化过程。它们需要从湿实验室(通过定点、随机或组合诱变)构建突变体文库，以从文库中检索有限数量样本的序列和/或结构(称为“学习数据集”)并评估每个取样突变体的适应度。它们还要求从电脑模拟提取每个突变体的描述符，以使用多变量统计方法来建立描述符和适应度(学习阶段)之间的关系，并建立模型来对未经实验测试的突变体进行预测。已经提出了基于3d结构的称为定量结构-功能关系(qfsr)的方法(damborskyj,prot.eng.(1998)jan；11(1):21-30)。提出了仅基于序列而不基于3d结构并使用统计建模进行电脑模拟合理筛选的其它方法(foxr.等，proteineng.(2003)16(8):589-97；foxr.,journaloftheoreticalbiology(2005),234:187–199；minshullj.等，curropinchembiol.2005apr；9(2):202-9；foxr.等，naturebiotechnology(2007),25(3):338-344；foxr.andhuismangwtrendsbiotechnol.2008mar；26(3):132-8)。最著名的是prosar(foxr.,journaloftheoreticalbiology(2005),234:187–199；foxr.等，naturebiotechnology(2007),25(3):338-344)，其基于二进制编码(0或1)。qsfr方法是有效的，并且考虑到与非变体残基可能相互作用的信息。然而，qsfr需要目前仍然有限的关于3d蛋白质结构的信息，并且该方法还慢。相比之下，prosar不需要知道3d结构，因为它仅基于一级序列计算，并且可以使用线性和非线性模型。然而，prosar仍然有缺点，且其筛选能力有限。特别地，只有那些正在变化的残基被包括在建模中，因此，缺少关于突变残基和其他非变体残基之间可能的相互作用的信息。prosar依赖于不考虑氨基酸的物理化学或其他分子性质的突变的二进制编码(0或1)。另外，(i)可以测试的新序列只是在用于构建模型的学习集中使用的位置上具有突变或突变组合的序列；(ii)待筛选的新序列中突变的位置数目不能与训练集中的突变数目不同；和(iii)当在超级计算机上引入非线性项以建立模型时计算时间是非常长的(对于100个非线性项，直至2周)。因此，仍然需要通用且快速的电脑模拟方法来帮助蛋白质定向进化过程。本发明提供了满足这些要求并且基于数字信号处理(dsp)的方法。数字信号处理技术是分解和处理信号以揭示其中内含的信息的分析程序。信号可能是连续的(无休止的)，或离散的诸如蛋白质残基。在蛋白质中，傅里叶变换方法已被用于生物序列(dna和蛋白质)比较，蛋白质家族的表征和模式识别、分类和其它基于结构的研究(诸如对称性和重复结构单元或模式的分析，二级/三级结构的预测，疏水核心、基序、保守结构域的预测，膜蛋白的预测，保守区域的预测，蛋白质亚细胞定位的预测)，用于氨基酸序列中二级结构含量的研究和用于蛋白质中的周期性的检测。最近提出了用于检测蛋白质结构中螺线管结构域的新方法。数字信号处理技术帮助了蛋白质相互作用的分析(cosici.,ieeetransbiomedeng.(1994)41(12):1101-14)，并使生物学功能性可计算。这些研究已在nwankwon.和sekerh.(jproteomicsbioinform(2011)4(12):260-268)中进行了详细的综述。在这些方法中，首先使用来自数据库aaindex(kawashima，s.和kanehisa，m.nucleicacidsres.(2000)，28(1)：374；kawashima，s等人，nucleicacidsres.jan2008；36)的可用aaindex之一将蛋白质残基转化成数值序列，该数值序列代表每个氨基酸的生物化学性质或物理化学参数。然后通过离散傅立叶变换(dft)处理这些数值序列，以信息谱的形式呈现蛋白质的生物学特征。该程序称为信息频谱法(ism)(veljkovicv.等，ieeetransbiomedeng.1985may；32(5):337-41)。ism程序已被用于研究钙结合蛋白(viaria.等，computapplbiosci.1990apr；6(2):71-80)和流感病毒(veljkovicv.，等bmcstructbiol.2009apr7；9:21,veljkovicv.等bmcstructbiol.2009sep28；9:62)中的主要排列(principalarrangement)。采用称为电子-离子相互作用电位(eiip)的氨基酸参数的ism的变体被称为共振识别模型(rrm)。在该程序中，生物功能性被呈现为频谱特征。这种物理数学过程基于以下事实：具有相同生物学特性的生物分子在其价电子振荡并随后在电磁场中反射(reverberate)时识别它们自己并生物关联于(bio-attach)它们自己(cosici.,ieeetransbiomedeng.(1994)41(12):1101-14；cosici.,theresonantrecognitionmodelofmacromolecularbioactivitybirkhauserverlag,1997)。共振识别模型涉及四个步骤(参见nwankwon.和sekerh.,jproteomicsbioinform(2011)4(12):260-268)：-步骤1：将蛋白质残基转化成电子-离子相互作用电位(eiip)参数的数值。-步骤2：补零/上采样。该过程使用补零来在任何位置上填充待分析的蛋白质的序列中的缺口，因为信号处理要求所有蛋白质的窗口长度相同。-步骤3：使用快速傅里叶变换(fft)处理数值序列以产生频谱特征(sc)，并在步骤4中将其逐点相乘以产生交叉频谱(cs)特征。-步骤4：交叉频谱分析：交叉频谱(cs)分析表示频谱特征(sc)的逐点乘法。因此，所述cs分析已经被定性地用于基于配体与受体频谱之间的共同频率(共振)来预测例如配体-受体结合。另一个实例是通过将rrm应用于ha-rasp21蛋白质序列来预测ras样活性存在与否，即转化细胞的能力存在与否。这些现有技术方法提供的信息是有用的，但又不足以鉴定通过定向进化产生的最有价值的蛋白质突变体。发明概述因此，本发明涉及用于预测蛋白质的至少一个适应度值的方法，该方法在计算机上执行，并且包括以下步骤：-根据蛋白质数据库将蛋白质的氨基酸序列编码成数值序列，数值序列包含序列中每个氨基酸的值；-根据数值序列计算蛋白质谱；和对于每个适应度：-将计算的蛋白质谱与预定数据库的蛋白质谱值进行比较，所述数据库包含针对所述适应度的不同值的蛋白质谱值，-根据比较步骤预测所述适应度的值。因此，本发明人开发的方法涉及对蛋白质谱的定量分析，这使得可能预测蛋白质的适应度值，而且还可能预测给定活性的存在与否。根据本发明的其它有利方面，该方法包括单独地或根据所有技术上可能的组合采用的以下特征中的一个或多个特征：-计算的蛋白质谱包括至少一个频率值，并将计算的蛋白质谱就每个频率值与所述蛋白质谱值进行比较；-在蛋白质谱计算步骤期间，将傅里叶变换(诸如快速傅里叶变换)应用于编码步骤后获得的数值序列；-每个蛋白质谱验证以下方程式：其中j为蛋白质谱|fj|的指数(index-number)；数值序列包括指示xk的n数值，其中0≤k≤n-1且n≥1；以及i定义虚数，使i2＝-1；-在编码步骤期间，蛋白质数据库包括生物化学或物理化学特性值的至少一个指数，每个特性值被赋予相应的氨基酸；并且对于每个氨基酸，数值序列中的值等于给定的指数中所述氨基酸的特性值；-在编码步骤期间，蛋白质数据库包括特性值的几个指数；并且所述方法还包括基于样本蛋白质的测量的适应度值与根据每个指标对所述样本蛋白质预先获得的预测的适应度值的比较来选择最佳指标的步骤；然后使用所选择的指数进行编码步骤；-在选择步骤期间，所选择的指数是具有最小均方根误差的指数，其中每个指数的均方根误差验证以下方程式：其中yi为第i个样本蛋白质的测量的适应度，为第i个样本蛋白质用第j个指数预测的适应度，以及s为样本蛋白质的数量；-在选择步骤期间，所选择的指数是决定系数最接近1的指数，其中每个指数的决定系数验证以下方程式：其中yi为第i个样本蛋白质的测量的适应度，为第i个样本蛋白质用第j个指数预测的适应度，s为样本蛋白质的数量，为s样本蛋白质的测量的适应度的平均值，以及为s样本蛋白质的预测的适应度的平均值；-该方法还在编码步骤之后并且在蛋白质谱计算步骤之前包括以下步骤：+通过将数值序列的每个值减去数值序列值的平均值对通过编码步骤获得的数值序列进行归一化；然后对归一化的数值序列进行蛋白质谱计算步骤；-该方法还在编码步骤之后并且在蛋白质谱计算步骤之前包括以下步骤：+通过在所述数值序列的一端添加m个零对通过编码步骤获得的数值序列进行补零，其中m等于(n-p)，其中n是预定整数，p是所述数值序列中的值的数目；然后对补零步骤后获得的数值序列进行蛋白质谱计算步骤；-比较步骤包括在针对所述适应度的不同值的蛋白质谱值的预定数据库中，根据预定标准确定与计算的蛋白质谱最接近的蛋白质谱值，所述适应度的预测值则等于在所述数据库中与确定的蛋白质谱值相关联的适应度值；-在蛋白质谱计算步骤期间，根据几个频率范围计算所述蛋白质的几个蛋白质谱，且其中，在预测步骤期间，根据比较步骤针对每个蛋白质谱估计适应度的中间值，然后使用中间适应度值来计算适应度的预测值，优选利用对中间适应度值进行回归，诸如偏最小二乘回归来计算；以及该方法包括以下步骤：-根据计算的蛋白质谱分析蛋白质，以筛选突变体文库，该分析优选使用多因辨别分析(factorialdiscriminantanalysis)或主成分分析进行。本发明还涉及包括软件指令的计算机程序产品，当所述指令通过计算机执行时，执行如上定义的方法。本发明还涉及用于预测蛋白质的至少一个适应度值的电子预测系统，该预测系统包括：-编码模块，其被配置用于根据蛋白质数据库将氨基酸序列编码为数值序列，所述数值序列包含序列中每个氨基酸的值；-计算模块，其被配置用于根据所述数值序列计算蛋白质谱；和-预测模块，对于每个适应度，其被配置用于：+将计算的蛋白质谱与预定数据库的蛋白质谱值进行比较，所述数据库包含针对所述适应度的不同值的蛋白质谱值，以及+根据所述比较预测所述适应度的值。附图简述在阅读以下仅通过实例并参照附图给出的描述后，本发明将被更好地理解，其中：-图1是用于预测蛋白质的至少一个适应度值的电子预测系统的示意图，该预测系统包括被配置用于将氨基酸序列编码为数值序列的编码模块，被配置用于根据数值序列计算蛋白质谱的计算模块；以及被配置用于预测每个适应度的至少一个值的预测模块；-图2是根据本发明的用于预测蛋白质的至少一个适应度值的方法的示意性流程图；-图3表示针对人glp1蛋白的天然和突变形式获得的蛋白质谱的曲线；-图4是说明细胞色素p450家族的一组蛋白的热稳定性的预测值和测定值的一组点，每个点与相应的蛋白质相关，其中纵坐标对应于预测值，横坐标对应于测量值，使用蛋白质谱中包括的所有频率；-图5和图6是分别针对来自细胞色素p450家族的蛋白质组的训练和验证子集获得的类似图4的图示，训练子集用于计算包含针对热稳定性的不同值的蛋白质谱值的数据库，验证子集与训练子集不同，并用于测试与相应测量值相比较的预测值的相关性；-图7是关于一组glp1突变体的结合亲和力的预测值和测量值的类似于图4的图示；-图8是关于一组glp1突变体的效力的预测值和测量值的类似于图4的图示；-图9和图10是关于分别针对一组肠毒素see和sea的训练和验证子集获得的热稳定性的预测值和测量值的类似于图4的图示，训练子集用于计算含有针对所述热稳定性的不同值的蛋白质谱值的数据库，所述验证子集与训练子集不同，并用于测试预测值的相关性；-图11和12是关于分别针对一组tnf突变体的训练和验证子集获得的结合亲和力的预测值和测量值的类似于图4的图示，训练子集用于计算含有针对所述结合亲和力的不同值的蛋白质谱值的数据库，所述验证子集与所述训练子集不同，并用于测试所述预测值的相关性；-图13是使用来自蛋白质谱的频率值的选择的类似于图4的图示；-图14是关于一组环氧化物水解酶家族的蛋白的对映选择性的预测值和测量值的类似于图4的图示；-图15表示512个环氧化物水解酶突变体的文库的筛选；-图16表示使用多变量分析(主成分分析)进行蛋白质筛选来进行的10个环氧化物水解酶突变体的蛋白质谱的分类。-图17是关于布鲁顿酪氨酸激酶变体的蛋白质表达水平的预测值和测量值的类似于图4的图示；-图18是关于k562细胞系中rna的mrna表达水平的预测值和测量值的类似于图4的图示；-图19是关于心脏细胞中蛋白质的蛋白质表达水平的预测值和测量值的类似于图4的图示；和-图20是关于肾细胞中蛋白质的蛋白质表达水平的预测值和测量值的类似于图4的图示。优选实施方案的详述如本文中所用，“蛋白质”意指通过肽键连接在一起的至少2个氨基酸。术语“蛋白质”包括蛋白质、寡肽、多肽和肽。肽基可包含天然存在的氨基酸和肽键，或合成的拟肽结构，即“类似物”，诸如类肽。氨基酸可以是天然存在的或非天然存在的。在优选实施方案中，蛋白质包含至少10个氨基酸，但较少的氨基酸也可处理。蛋白质的“适应度”是指其对标准诸如催化效力，催化活性，动力学常数，km，keq，结合亲和力，热稳定性，溶解度，聚集，效力，毒性，变应原性，免疫原性，热力学稳定性，柔性的适应，根据本发明，“适应度”也称为“活性”，并且在下面的描述中认为适应度和活性是指相同的特征。催化效力通常以s-1.m-1表示，是指kcat/km的比率。催化活性通常以mol.s-1表示，是指酶催化中的酶活性水平。动力学常数kcat通常以s-1表示，是指定量反应速度的数值参数。km通常以m表示，是指反应速度为其最大值的一半时的底物浓度。keq通常以(m、m-1或无单位)表示，为表示化学反应中的化学平衡的量，结合亲和力通常以m表示，是指一种或多种蛋白质与配体(肽或小化学分子)之间的相互作用的强度。热稳定性通常以℃表示，通常是指测量的活性t50，其定义为在10分钟的孵育时间后50％的蛋白质被不可逆变性时的温度。溶解度通常以mol/l表示，是指在溶液变饱和之前每升溶液可溶解的物质(溶质)的摩尔数。聚集通常使用聚合指数(来自280nm和340nm处的简单吸收测量)来表示，并且是指其中错误折叠的蛋白质在细胞内或细胞外聚集(即，堆集和簇集在一起)的生物学现象。效力通常以m表示，是指以产生给定强度的效果所需的量表示的药物活性的量度。毒性通常以m表示，是指物质(毒素或毒药)可对人或动物造成伤害的程度。变应原性通常以生物等效变态反应单位/ml(bau/ml)表示，是指抗原性物质产生立即超敏反应(变态反应)的能力。免疫原性通常表示为样品中抗体量的单位，是指特定物质(诸如抗原或表位)在人或动物体内引发免疫应答的能力稳定性通常表示为δδg(kcal/mol-1)，是指蛋白质的快速、可逆和协同地解折叠和重折叠的热力学稳定性。柔性通常以a°表示，是指蛋白质紊乱和构象变化。在图1中，用于预测蛋白质的至少一个适应度值的电子预测系统20包括数据处理单元30、显示屏32和用于将数据输入数据处理单元30的输入装置34。数据处理单元30由例如存储器40和与存储器40联接的处理器42组成。显示屏32和输入装置34本身是已知的。存储器40适用于存储编码计算机程序50(其被配置用于根据蛋白质数据库51将氨基酸序列编码为数值序列)和计算计算机程序52(其被配置用于根据数值序列计算下文所示的蛋白质谱|fj|，其中j为蛋白质谱的指数号)。存储器40还适用于存储建模计算机程序54，所述程序54被配置用于预先确定包含针对所述适应度的不同值的蛋白质谱值的蛋白质谱数据库55。存储器40适用于存储预测计算机程序56，所述程序56被配置用于对于每个适应度，将计算的蛋白质谱与所述预定数据库的蛋白质谱值进行比较，以及用于根据所述比较预测所述适应度的值；以及任选地进一步用于筛选突变体文库。可选地，存储器40适用于存储筛选计算机程序58，所述程序58被配置用于根据计算的蛋白质谱分析蛋白质，从而筛选突变体文库，该分析优选是多因辨别分析或主成分分析。处理器42被配置用于执行编码、计算、建模、预测和筛选计算机程序50、52、54、56、58中的每一个。编码、计算、建模、预测和筛选计算机程序50、52、54、56、58，当由处理器42来执行它们时，分别形成用于根据蛋白质数据库将氨基酸序列编码成数值序列的编码模块；用于根据数值序列计算蛋白质谱的计算模块；用于预先确定含有蛋白质谱值的数据库的建模模块；用于将计算的蛋白质谱与所述预定数据库的蛋白质谱值进行比较，并根据所述比较预测所述适应度的值以及进行筛选的预测模块；以及用于根据计算的蛋白质谱分析蛋白质的筛选模块。或者，编码模块50、计算模块52、建模模块54、预测模块56和筛选模块58以可编程逻辑组件的形式，或以专用集成电路的形式存在。编码模块50适用于根据蛋白质数据库51将氨基酸序列编码成数值序列，数值序列包含序列的每个氨基酸的值xk。数值序列由p值xk构成，其中0≤k≤p-1，p≥1，k和p为整数。将蛋白质数据库51例如存储在存储器40中。或者，将蛋白质数据库51存储在与存储器40不同的远程存储器(未显示的)中。蛋白质数据库51优选为氨基酸指数数据库，也称为aaindex。氨基酸指数数据库可从http://www.genome.jp/dbget-bin/www_bfind？aaindex(发布版9.1，8月6日)获得。蛋白质数据库51包括生物化学或物理化学特性值的至少一个指数，每个特性值被赋予相应的氨基酸。蛋白质数据库51优选地包括生物化学或物理化学特性值的几个指数。每个指数对应于例如aaindex代码，如以下将根据各个实例所示的。用于编码氨基酸序列的选择的aaindex代码是例如：扩展结构的d归一化频率、d电子-离子相互作用电位值、总蛋白质的aa组成的dsd、来自ifh量表(scale)的dpk-c或d权重。为了编码氨基酸序列，则使编码模块50适用于确定每个氨基酸在给定的指数中所述氨基酸的特性值，数值序列中的每个编码值xk则等于相应的特性值。另外，以任选的方式，当蛋白质数据库51包括特性值的几个指数时；编码模块50被进一步配置用于基于样本蛋白质的测量的适应度值与先前根据每个指数针对所述样本蛋白质获得的预测的适应度值的比较，选择最佳指标；然后用于使用所选择的指数编码氨基酸序列。所选择的指数是例如具有最小均方根误差的指数，其中每个指数的均方根误差验证以下方程式：其中yi是第i个样本蛋白质的测量的适应度，是第i个样本蛋白质用第j个指数预测的适应度s样品蛋白质的数量。或者，所选择的指数是具有最接近于1的决定系数的指数，其中每个指数的决定系数验证以下方程式：其中yi第i个样本蛋白质的测量的适应度，是第i个样本蛋白质用第j个指数预测的适应度，s是样品蛋白的数量，是s样本蛋白质的测量的适应度的平均值，以及是s样本蛋白质的预测的适应度的平均值。另外，以任选的方式，编码模块50被进一步配置用于对所获得的数值序列进行归一化，例如通过将数值序列的每个值xk减去数值序列值的平均值换句话说，表示为的每个归一化值验证以下方程式：平均值为例如算术平均值并且满足：或者，平均值是几何平均值、调和平均值或均方平均值。另外，以任选的方式，编码模块50被进一步配置用于通过在所述数值序列的一端添加m个零来对所获得的数值序列进行补零，其中m等于(n-p)，其中n是预定整数并且p是所述数值序列中的值的初始数量。因此，n是补零后的数值序列中的值的总数。计算模块52被配置用于根据数值序列计算蛋白质谱。计算的蛋白质谱包括至少一个频率值。计算模块52被配置用于计算蛋白质谱|fj|，优选地通过对所获得的数值序列应用傅里叶变换，诸如快速傅里叶变换。每个蛋白质谱|fj|因此验证例如以下方程：其中j是蛋白质谱|fj|的指数指数号；以及i定义虚数，使得i2＝-1。另外，当通过编码模块50对数值序列进行归一化时，计算模块52被进一步配置用于对归一化的数值序列进行蛋白质谱计算。换句话说，在该情况下，每个蛋白质谱|fj|因此验证例如以下方程式：另外，当通过编码模块50对数值序列进行补零时，计算模块52被进一步配置用于对进一步获得以用于补零的数值序列计算蛋白质谱|fj|。换句话说，在该情况下，每个蛋白质谱|fj|因此验证例如以下方程式：另外，当通过编码模块50对数值序列进行归一化和补零时，计算模块52被进一步配置用于对进一步获得以用于补零的归一化数值序列计算蛋白质谱|fj|。换句话说，在该情况下，每个蛋白质谱|fj|因此验证例如以下方程式：建模模块54适用于根据从编码模块50下发的学习数据和从计算模块52下发的学习蛋白质谱，预先确定蛋白质谱数据库55(也称为模型)。学习蛋白质谱对应于学习数据并且学习数据各自与给定的适应度相关，并且优选地针对所述适应度的不同值。蛋白质谱数据库55包含针对每个适应度的不同值的蛋白质谱值。优选地，使用至少10个蛋白质谱和10个不同的适应度来构建蛋白质谱数据库55。当然，蛋白质谱数量和相关蛋白质适应度越高；在适应度预测方面的结果越好。在下面的实例中，用作学习数据的蛋白质谱和适应度的数量范围为8至242个(242个蛋白质谱和242个蛋白质适应度；8个蛋白质谱和8个蛋白质适应度)。对于每个适应度，预测模块56适用于将计算的蛋白质谱与蛋白质谱数据库55的蛋白质谱值进行比较，并根据所述比较来预测所述适应度的值。预测模块56被进一步配置用于在蛋白质谱数据库55中根据预定标准确定最接近计算的蛋白质谱的蛋白质谱值，然后所述适应度的预测值等于在蛋白质谱数据库55中与所确定的蛋白质谱值相关联的适应度值。预定标准是例如计算的蛋白质谱与蛋白质谱数据库55中包含的蛋白质频谱值之间的最小差异。或者，预定标准是计算的蛋白质谱与蛋白质谱数据库55中包含的蛋白质谱值之间的相关系数r或决定系数r2。当蛋白质谱|fj|包含几个频率值时，对于每个频率值将计算的蛋白质谱|fj|与蛋白质谱值进行比较。或者，仅将一些频率值考虑用于计算的蛋白质谱|fj|与所述蛋白质谱值的比较。在该情况下，例如根据频率值与适应度的相关性来对频率值进行排序，并且仅将最佳频率值考虑用于计算的蛋白质谱的比较。另外，以任选的方式，预测模块56被进一步配置用于在根据几个频率范围计算所述蛋白质的几个蛋白质谱时估计每个蛋白质谱的适应度的中间值。然后，预测模块56被进一步配置用于通过对所述中间适应度值进行回归(诸如偏最小二乘回归，也称为plsr的)来计算适应度的预测值。或者，预测模块56被配置用于使用人工神经网络(ann)来计算适应度的预测值，其中输入变量是所述中间适应度值，输出变量是适应度的预测值。另外，以任选的方式，预测模块56允许获得突变体文库的筛选，因为将在下面在图15的图示中以对映选择性作为适应度来描述其。另外，以任选的方式，筛选模块58适用于根据计算的蛋白质谱分析蛋白质，以及使用数学处理(诸如多因辨别分析(factorialdiscriminantanalysis)或主成分分析，然后例如进行k-means)根据其各自的蛋白质谱对蛋白质序列进行分类。例如，可以进行分类，以确定在蛋白质谱家族中是否存在不同的组：具有高、中和低适应度的组；作为例子，具有适应度的表达的组和没有适应度的表达的组。在下文中，将根据图16进一步说明该筛选。现在将参照图2描述根据本发明的电子预测系统20的操作，所述图2表示用于预测蛋白质的至少一个适应度值的方法的流程图。在初始步骤100中，编码模块50根据蛋白质数据库51将蛋白质的氨基酸序列编码成数值序列。编码步骤100可使用氨基酸指数数据库(也称为aaindex)来进行。在编码步骤100期间，编码模块50针对每个氨基酸确定给定的指数中(例如在给定的aaindex代码中)的所述氨基酸的特性值，然后下发等于所述特性值的编码值xk。另外，当蛋白质数据库51任选地包括几个特性值的指数时；编码模块50还基于样本蛋白质的测量的适应度值与先前根据每个指数针对所述样本蛋白质获得的预测的适应度值的比较，选择最佳指数；然后使用所选择的指数编码氨基酸序列。例如，使用方程式(1)或方程式(2)选择最佳指数。此外，编码模块50任选择地对所获得的数值序列进行归一化，例如通过根据方程式(3)将数值序列的每个值xk减去数值序列值的平均值。另外，编码模块50任选地通过在所述数值序列的一端添加m个零来对所获得的数值序列进行补零。在编码步骤100结束时，编码模块50将学习数值序列和验证数值序列传送至计算模块52，并将学习数据传送至建模模块54。两个蛋白质频谱的实例示于图3中，第一曲线102表示人glp1蛋白的天然形式的蛋白质谱，第二曲线104表示人glp1蛋白的突变形式(单突变)的蛋白质谱。对于每个曲线102、104，蛋白质谱的连续离散值彼此相连。在下一步骤110中，计算模块52计算从编码模块50下发的每个数值序列的蛋白质谱|fj|。对应于学习数值序列的蛋白质谱也称为学习谱，对应于验证数值序列的蛋白质谱也称为验证谱。步骤110也称为谱变换步骤。蛋白质谱|fj|优选地通过使用傅立叶变换诸如快速傅立叶变换，例如取决于任选的归一化和/或补零，根据方程式(5)至(8)中的方程式来计算。然后，建模模块54在步骤120中根据在编码步骤100期间获得的学习数据和在谱变换步骤110期间获得的学习蛋白质谱，确定蛋白质谱数据库55。在步骤130中，对于每个适应度，预测模块56将计算的蛋白质谱与从蛋白质谱数据库55中下发的蛋白质谱值进行比较，然后根据所述比较预测的适应度值。更精确地，预测模块56在蛋白质谱数据库55中根据预定标准确定最接近计算的蛋白质谱的蛋白质谱值，然后预测的适应度值等于与蛋白质谱数据库55中确定的蛋白质谱值相关联的适应度值。任选地，只将一些频率值考虑用于计算的蛋白质谱|fj|与所述蛋白质谱值的比较。另外，当根据几个频率范围任选地计算所述蛋白质的几个蛋白质谱时，预测模块56估计每个蛋白质谱的中间适应度值。然后，预测模块56通过对所述中间适应度值进行回归(诸如plsr)来计算预测的适应度值。或者，人工神经网络(ann)被预测模块56用于基于所述中间适应度值计算适应度的预测值。然后，预测模块56允许通过根据预测的适应度对蛋白质谱进行分级来进行蛋白质筛选。最后和任选地，筛选模块58在步骤140中使用数学处理(诸如多因辨别分析或主成分分析)根据其各自的蛋白质谱对蛋白质序列进行分析、分类。或者，例如通过与预定值进行比较，直接对计算的蛋白质谱进行分析以筛选突变体文库。因此，其允许获得更好的突变体文库筛选。此步骤也称为多变量分析步骤。应当注意，分析步骤140紧接谱变换步骤120，并且另外可在分析步骤140之后执行预测步骤130，以用于预测一些或所有分类的蛋白质的适应度值。将潜在分量计算为原始变量的线性组合计算；选择潜在分量的数量以使rmse(均方根误差)最小化。将潜在分量计算为原始变量(频率值)的线性组合；通过逐个添加分量来选择潜在分量的数量以最小化rmse(均方根误差)。实施例将根据以下实施例来进一步说明本发明。实施例1：细胞色素p450(图4至6)在本实施例中，使用以下aaindex代码将细胞色素p450的氨基酸序列编码成数值序列：扩展结构的d归一化频率(maxfield和scheraga,biochemistry.1976；15(23):5138-53)。第一数据集(来自li等人，2007:natbiotechnol25(9):1051–1056.；romero等，pnas.2013:january15，第11卷，n°3:e193-e201)来自围绕细胞色素p450家族，特别是细胞色素p450bm3a1、a2和a3的序列/稳定性-功能关系的研究，该研究旨在提高细胞色素的热稳定性。含有血红素的氧化还原酶的多功能细胞色素p450家族羟化广泛的底物，以产生具有重大医学和工业重要性的产品。利用从这三种不同亲本中的任一亲本继承的八个连续片段构建新的嵌合蛋白。测量的活性是t50，其被定义为在10分钟的孵育时间后50％的蛋白质被不可逆地变性时的温度。输出数据集(out-comingdataset)由242个变体序列组成，其中t50实验值范围为39.2至64.48℃。cyp102a1及其同源物cyp102a2(a2)和cyp102a3(a3)的血红素结构域的重组允许产生由八个片段组成的242个嵌合p450序列，每个片段选自三个亲本之一。根据片段组成书写嵌合体：例如，23121321代表继承亲本a2的第一片段、a3的第二片段、a1的第三个片段(依此类推)的蛋白质。表2：cytp450学习集图4显示使用留一交叉验证(loocv)r2＝0.96和rmse＝1.21对蛋白质序列的整个集合进行模型之后获得的结果。这表明可以使用这种方法捕获与蛋白质的适应度相关的信息。嵌合体t50嵌合体t50嵌合体t501133221247.83131323251.92221322350.83233223149.42333222146.421331332522331311156.92211133250.91131333353.82333331145.72233222250.332311323523133133147.321131121532313231144.52123123350.62123223249.52133321155.92111212250.331212232513231233357.82211321151.12321321147.42231233259.123333233513223213143.92231233363.51333321145.72213321247.21232233347.92321331149.52131331156.92131223162.83233232348.521332231602231133360.12221321250.521113312532131123163.22213221246.62231222356.22131221159.32111133362.42223212149.72221233358.23211323247.93133223349.9表3：cytp450测试组图5和图6给出了模型预测细胞色素p450的突变的组合的能力。此处，数据集被在196个序列中被分为学习序列，在46个序列中被分为验证序列。实施例2：人胰高血糖素样肽-1(glp1)预测的类似物(图7和8)在本实施例中，使用以下aaindex代码将glp1的氨基酸序列编码为数值序列：d电子-离子相互作用电位值(cosic,ieeetransbiomedeng.1994dec；41(12):1101-14.)。他司鲁肽和extendin-4是glp1类似物，其作为胰高血糖素样肽(glp)受体的肽激动剂，并且正处于下用于治疗ii型糖尿病的临床开发(他司鲁肽)中。人glp1haegtftsdvssylegqaakefiawlvkgr(seqidno:1)他司鲁肽haegtftsdvssylegqaakefiawlvkar(seqidno:2)已实施本发明的方法以提供相对于天然人glp1和他司鲁肽提高结合亲和力(与受体的相互作用)和/或提高效力(受体-腺苷酸环化酶活性的激活)的glp1受体的候选激动剂。已经始于人glp1的序列，通过进行单点位点饱和诱变：用19个其它天然氨基酸置换氨基酸序列的每个位置来通过电脑模拟设计突变体文库。因此，如果蛋白质序列由n＝30个氨基酸组成，则生成的文库将由30×19＝570个单点变体组成。已经运行单点突变的组合。adelhorstk等(jbiolchem.1994mar4；269(9):6275-8)先前描述了通过ala扫描(即通过用l-丙氨酸连续替换每个氨基酸以鉴定与glp-1受体相互作用所需的侧链官能团)制备的一系列glp-1类似物。在l-丙氨酸是亲本氨基酸的情况下，已用在胰高血糖素相应位置上发现的氨基酸进行了置换。已在针对大鼠glp-1受体的结合测定(ic50)中测定了这些类似物，并已进一步监测了效力(通过检测腺苷酸环化酶活性ec50测量的受体激活)。将这些类似物(30个单突变体)及其报告的活性(分别与野生型人glp1的ic50或ec50相比标准化的log(ic50)和log(ec50))用作学习数据集来建立预测模型(参见图7和图8)。表4:glp1学习集测试肽logic50glp1t5a0.54407glp1l14a0.23045glp1q17a0.04139glp1f22a2.54531表5:glp1测试序列(结合)测试肽logec50glp1v10a0.8451glp1f22a3.41497glp1v27a0.30103野生型glp10.41497表6:glp1测试序列(效力)对于结合亲和力，它们的活性范围为-0.62至2.55(logic50)，对于效力，它们的活性范围为-0.30至4.00(logec50)。结果表明，对于结合亲和力(图7)，r2和rmse分别为0.93和0.19，对于效力(图8)分别为0.94和0.28，因此表明可以以非常有效的方式捕获与两个适应度相关的信息。针对人glp1、他司鲁肽和最佳电脑模拟的类似物(基于预测模型)评估的结合和效力如表7所示：结合(ic50)nm竞争力(ec50)nm人glp10.272.6他司鲁肽0.790.39最佳电脑模拟的类似物0.0020.021表7：对人glp1和类似物评估的结合和效力对于glp1的肽配体类似物对其受体的结合亲和力实现了135倍的提高。获得124倍的效力提高。这说明本发明的方法可以用于同时提高多个参数。实施例3：环氧化物水解酶的对映选择性的演化(图14和15)在本实施例中，使用以下aaindex代码：总蛋白质的aa组成的dsd(nakashima等，proteins.1990；8(2):173-8)将环氧化物水解酶的氨基酸序列编码成数值序列。对映选择性是在化学反应中一种立体异构体相对于另一种的优先形成。对映选择性对于许多与工业相关的化学物质的合成很重要，并且难以实现。绿色化学利用重组酶(因为酶具有很高的特异性)来合成目标化学产品。因此，在绿色化学中特别寻求具有提高的效率的酶。reetz，等(ang2006feb13；45(8):1236-41)描述了来自黑曲霉(aspergillusniger)的环氧化物水解酶的对映选择性突变体作为缩水甘油醚1的水解动力学拆分(其中形成二醇(r)-2和(s)-2)中的催化剂的定向进化。该模型建立在reetz等(同上)描述的一组10个学习序列上。环氧化物水解酶δδg(kcal/mol)wt-0.85l215f-1.50a217n-1.17r219s-0.85l249y-0.85t317w-1.50t318v-0.85m329p-1.08l330y-0.85c350v-0.97表8：学习集将湿实验室中产生的32个突变体的结果与使用我们的方法预测的结果进行了比较。定量值显示在图14的右侧：表示实验值和预测值。所获得的预测值非常接近实验值，平均偏差为-0.011kcal/mol。这表明即使对少量学习序列和学习数据，也可以获得具有提高的参数的良好突变体。在图15中，构建并筛选了512个突变体的文库。在湿实验室中鉴定的最佳突变体确实是一个好的突变体(箭头150)，但不是最好的。最好的突变体由图15中的椭圆160标识。箭头170指向野生型蛋白质。表9：测试序列实施例4：肠毒素sea和see(图9和10)的热稳定性(tm)的预测在本实施例中，使用以下aaindex代码：dpk-c(fasman，1976)将肠毒素的氨基酸序列编码成数值序列。第四数据集(来自cavallina.等，2000:biolchem.jan21；275(3):1665-72.)与肠毒素see和sea的热稳定性相关。超抗原(sag)，诸如葡萄球菌肠毒素(se)，是已知引起食物中毒或中毒性休克的非常强效的t细胞活化蛋白。通过将这些肠毒素与肿瘤反应性抗体融合，已经探索了将由它们诱导的强细胞毒性用于癌症治疗。tm定义为变性温度ec50值，对于由12个蛋白质序列组成的数据集(wtsae+wtsee+10个突变体，包括形成1个单突变至21个多重突变)，范围为55.1至73.3℃。表10.sea和see的突变区的细节。see/aa、-f、-h和-ah是分别具有来自sea的区域a、f、a和a+h的see，而sea/e-bdeg是具有来自see的区域b+d+e+g的sea。肠毒素tmsea_d227a55.1sea_h187a57.5sea_233aa(野生型)61.4sea/e-bdeg68.4see/a-h69see/a-a_d227a69.3see_233aa(野生型)71.3see/a-a75.3表11:学习集肠毒素tmsee_a-f70see_a-ah69.1see_d227a67.4sea_d227a_f47a55.4表12：测试序列将我们的预测与湿实验室结果进行了比较(cavallina.2000)。此处再次地，使用小的学习序列(8个学习序列)和学习数据，有可能捕获与热稳定性相关的信息，并预测新突变体的该参数。应当注意，在对应于图10的验证集的蛋白质序列(4个蛋白质序列)中，2个在对应于图9的训练集中未被采样的位置处包含突变(具有7个新突变的1个序列，和具有超过2的1个新突变的1个序列(1sequenceavec1newmutationover2))。因此，这些结果证实，有可能鉴定新的突变体，包括在训练集中未被采样的突变的位置。结果表明，对于训练集(图9)，r2和rmse分别为0.97和1.16，对于验证集(图10)，r2和rmse分别为0.96和1.46，因此表明在该情况下可以高效地预测与热稳定性相关的信息。实施例5：具有改变的受体选择性的突变型tnf(图11和12)在本实施例中，使用以下aaindex代码：来自ifh量表的d权重(jacobsandwhite,biochemistry.1989；28(8):3421-37)将tnf的氨基酸序列编码为数值序列。肿瘤坏死因子(tnf)是抑制致癌作用并排除感染性病原体以维持内环境稳定的重要细胞因子。tnf激活其两种受体，tnf受体tnfr1和tnfr2。mukaiy等(jmolbiol.2009jan30；385(4):1221-9)产生了仅激活一种tnfr的受体选择性tnf突变体。mukai等(同上)公开的21个突变体的受体选择性已使用数据突变体(包括1个单突变至6个多重突变的wt+20个突变体)和在本文中公开为学习数据集进行了预测。tnf多肽受体选择性wt0k11m,k65s,k90p,k98r,k112n,k128p0.079l29i0.079a84t,v85h,s86k,q88p,t89q0.544a84s,v85k,s86t,q88s,t89h0.663l29q,r32w0.826l29k,r31a,r32g,e146s,s147t0.924a84s,v85t,s86n,q88n,t89g0.869a84s,v85s,s86h,q88r,t89f1.079a84s,v85p,s86l,q88p,t89k1.217a84t,v85s,s86a,q88g,t89p1.230a84t,v85t,s86a,q88s,t89g1.310a145r,e146t,s147d1.301a145k,e146d,s147t2.870a145r,e146e,s147t2.228a145a,e146d,s147d1.949a145a,e146n,s147d2.462表13：tnf学习集基于elisa测量预测tnf与tnfr1(r1)和tnfr2(r2)的竞争性结合，如mukaiy等在论文中所述的。使用r1和r2的相对亲和力(％kd)来计算logr1/r2的比率。相对亲和力log10(r1/r2)的范围为0至2.87。在第一步中，该方法已被应用于整个数据集。对于tnf的结合亲和力，r2和rmse分别等于0.97和0.11。这再次证明，该方法能够捕获与适应度相关的信息。在第二步中，将17个突变体用作学习序列，4个用作验证序列。tnf多肽受体选择性l29t_r31g_r32y0.380l29t_r31k_r32y1.127l29t_r32f_e146t2.026a84s_v85k_s86t_q88t_t89h0.924表14：tnf测试序列结果表明，对于训练集(图11)r2和rmse分别为0.93和0.21，对于验证集(图12)r2和rmse分别为0.99和0.17，因此表明有可能使用该方法对tnf突变体优先(比率r1/r2)与一种类型的受体结合的能力建模。在所有上述实施例1至5中，为了进行预测，使用全蛋白质谱。在下面的实施例6中，我们证明了根据本发明的方法以非常高效的方式仅使用部分蛋白质谱进行工作。实施例6：使用来自蛋白质谱的频率值的选择来预测细胞色素p450的热稳定性(图13)在本实施例中，使用以下aaindex代码：扩展结构的d归一化频率(maxfield和scheraga,biochemistry.1976；15(23):5138-53)将细胞色素p450的氨基酸序列编码成数值序列。此处，使用来自蛋白质谱的最相关频率的选择来进行预测。根据其与适应度的相关性对频率值进行排序，并且仅考虑最佳频率值。数据集与实施例1中的相同。结果表明，r2和rmse分别为0.91和1.75，从而表明，仅使用部分(选择)来自蛋白质谱的频率也可高效地预测适应度(此处为热稳定性)。这说明可使用全蛋白质谱或来自蛋白质谱的部分(选择)频率来使用本发明的方法。实施例7：使用多变量分析对蛋白质谱进行分类以进行蛋白质筛选(图16)使用包括10个具有低适应度值和高适应度值(对映选择性)的蛋白质谱的环氧化物水解酶亚组(如在实施例3中)。进行了pca(主成分分析)。适应度的低值和高值分别在小椭圆180和大椭圆190中，因此表明应用于蛋白质谱的多变量分析有助于蛋白质筛选。x、y和z轴是从pca产生的3个主要成分，并且考虑到58.28％(以惯用的x、y和z轴，分别为21.51％、19.72％和16.05％)的与蛋白质谱收集相关的全局信息。因此，在上述实施例中获得的几个适应度的预测值与测量值之间的r2和rmse表明，根据本发明的预测系统20和方法允许对不同蛋白质的不同适应度值进行高效预测。另外，根据本发明的方法允许测试在其它位置上具有突变或突变的组合的新序列(验证/测试序列)，这些新序列用于学习序列集以用于建立模型。该方法还允许测试新序列(验证/测试序列)，所述新序列具有与学习序列集中使用的突变位置数相比不同的突变位置数。该方法还允许测试新序列，所述序列包括在训练集中未被采样的突变的位置。在该情况下，给出肠毒素作为实施该方法的实例。另外，该方法还允许测试新序列(验证/测试序列)，所述序列具有与用于构建模型的学习序列集的长度相比不同的长度(在氨基酸数量方面)。该方法使得能够使用相同的学习序列和一个或多个不同的编码aaindex和不同的适应度/活性值作为学习数据来预测用于学习序列或验证序列的适应度(验证/测试数据)：即使用该新方法预测蛋白质序列的2个或更多个活性/适应度的能力。glp1在本文档中用作实例：作为实例，使用相同的aaindex进行对glp1受体的结合亲和力的预测和效力的预测。使用该方法，有可能使用非常小的学习序列和学习数据来实现非常好的预测以及获得具有提高的适应度的突变体。作为例子，给出其中仅使用10个蛋白质序列的环氧化物水解酶(epoxydehydrolase)。该方法还允许使用嵌合蛋白而不是具有单点突变或单点突变的组合的蛋白质序列。在本文档中将细胞色素p450作为实例给出。使用不同p450的片段组合。本发明使得考虑到氨基酸序列中不同位置处的不同氨基酸之间的相互作用的作用成为可能。图3显示单点突变在每个频率下影响整个蛋白质谱。另外，该方法非常高效，因此在用于预测适应度的编码步骤之后需要不超过10分钟，同时使用50个蛋白质序列用于学习序列以及20个蛋白质序列用于验证序列。另外，蛋白质的“适应度”还指其对诸如蛋白质表达水平或mrna表达水平的标准的适应。因此，蛋白质的“适应度”是指其对标准诸如催化功效、催化活性、动力学常数、km、keq、结合亲和力、热稳定性、溶解度、聚集、效力、毒性、变应原性、免疫原性、热力学稳定性、柔性、蛋白表达水平和mrna表达水平的适应。如上所述，“适应度”也称为“活性”，并且在本说明书中认为适应度和活性是指相同的特征。根据以下实施例，将进一步说明适应度诸如蛋白质表达水平或mrna表达水平。实施例8：预测布鲁顿酪氨酸激酶变体的蛋白质表达水平(图17)在本实施例中，布鲁顿酪氨酸激酶(btk)是参与b细胞发育和成熟的关键蛋白。实际上，btk诱导成熟b细胞产生抗体，有助于消除感染。此外，该蛋白的功能障碍可能引起疾病如x连锁无丙种球蛋白血症或布鲁顿无丙种球蛋白血症(b细胞未能成熟)。将18种蛋白质变体(futatanit.等1998,《deficientexpressionofbruton'styrosinekinaseinmonocytesfromx-linkedagammaglobulinemiaasevaluatedbyaflowcytometricanalysisanditsclinicalapplicationtocarrierdetection.》,blood.1998jan15；91(2):595-602；kaneganeh.等2000,《detectionofbruton'styrosinekinasemutationsinhypogammaglobulinaemicmalesregisteredascommonvariableimmunodeficiency(cvid)inthejapaneseimmunodeficiencyregistry》,clinexpimmunol.2000jun；120(3):512-7)、野生型btk用于本实施例，如下表15所示。表15-btk变体的序列和蛋白质表达水平值在图17中，测量的活性对应于btk的蛋白质表达水平的体外测量，预测的活性对应于通过根据本发明的方法对btk的蛋白质表达水平预测的值。该值以蛋白质表达水平的百分比给出，其中100％对应于野生型的蛋白质表达水平。留一交叉验证(loocv)用于构建模型并预测蛋白质表达值。结果表明，r2和rmse分别为0.98和1.5，从而表明适应度(此处为蛋白质表达水平)也可被高效预测。蛋白质序列使用优化的相对分配能量(optimizedrelativepartitionenergies)-方法b(miyazawa-jernigan,1999self‐consistentestimationofinter‐residueproteincontactenergiesbasedonanequilibriummixtureapproximationofresidues.proteins:structure,function,andbioinformatics,34(1),49-68)编码。来自embl-ebi(http://www.ebi.ac.uk/gxa)的表达图谱提供了关于不同细胞类型、生物体部位、发育阶段、疾病和其它病况的动物和植物样品中的基因和蛋白质表达水平的信息。关于哪些基因产物在“正常”条件(例如组织、细胞类型)下以多少丰度存在的信息，技术人员可参考petryszak等,2016《expressionatlasupdate—anintegrateddatabaseofgeneandproteinexpressioninhumans,animalsandplants.》,nucl.acidsres.(04january2016)44(d1):d746-d752.doi:10.1093/nar/gkv1045。实施例9：预测k562细胞系中的mrna表达水平(图18)根据本发明的方法还适用于预测k562细胞系中的mrna表达水平值(fonsecana等2014rna-seqgeneprofiling-asystematicempiricalcomparison.plosone9(9):e107026.doi:10.1371/journal.pone.0107026)。由于rna序列与蛋白质序列之间存在共线性，因此使用与每个基因相关的蛋白质序列来建立模型。蛋白质相异在于反映rna序列和长度的氨基酸组成和长度。在下表16中提供了97种rna的数据集(序列和蛋白表达水平)。图18显示使用留一交叉验证(r2：0.81，rmse：10.3)获得的结果，从而说明根据本发明的方法也适用于通过与rna相关的蛋白质序列预测mrna表达水平。在双态模型中使用基于自我信息值(25％可达性)的亲水量表(hydropathyscale)编码蛋白质序列(naderi-manesh等，2001predictionofproteinsurfaceaccessibilitywithinformationtheory.proteins:structure,function,andbioinformatics,42(4),452-459)。实施例10：预测心脏细胞中不同蛋白质的蛋白质表达水平(图19)根据本发明的方法还用于预测心脏细胞中不同蛋白质的蛋白表达水平值。蛋白质相异在于氨基酸组成和长度。下表17提供了85个蛋白质的数据集(序列和蛋白质表达水平)。图19显示使用留一交叉验证(loocv,r2:0.87,rmse:20.22)获得的结果。在图19中，将值乘以10000。因此，根据本发明的方法也适用于预测心脏细胞中不同蛋白质的蛋白质表达水平值。使用暴露残基的百分比编码蛋白质序列(janin等，1978conformationofaminoacidside-chainsinproteins.journalofmolecularbiology,125(3),357-386)。实施例11：预测肾细胞中不同蛋白质的蛋白质表达水平(图20)在本实施例中，根据本发明的方法还用于预测肾细胞中不同蛋白质的蛋白质表达水平值。蛋白质相异在于氨基酸组成和长度。在下表18中提供了数据集(序列和蛋白质表达水平)。图20显示对于130个蛋白质序列使用留一交叉验证(loocv,r2:0.83,rmse:1.75)获得的结果。因此，根据本发明的方法也适用于预测蛋白质(特别是对于肾脏细胞中的不同蛋白质)的表达水平值。使用处于中间的相对偏好值(richardson-richardson,1988aminoacidpreferencesforspecificlocationsattheendsofalphahelices.science,240(4859),1648-1652)编码蛋白质序列。因此，在上述实施例中获得的几个适应度诸如蛋白质表达水平或mrna表达水平的预测值与测量值之间的r2和rmse表明，根据本发明的预测系统20和方法允许高效预测不同蛋白质或蛋白质变体的不同适应度值，也适用于预测蛋白质表达水平和mrna表达水平。sequencelisting<110>cadet,frédéric<120>用于预测蛋白质的至少一个适应度值的方法和电子系统、相关计算机程序产品<130>bet16p0498<150>ep15305552.0<151>2015-04-14<160>3<170>patentinversion3.5<210>1<211>30<212>prt<213>人工序列<220><223>合成肽<400>1hisalagluglythrphethrseraspvalsersertyrleuglugly151015glnalaalalysglupheilealatrpleuvallysglyarg202530<210>2<211>30<212>prt<213>人工序列<220><223>合成肽<400>2hisalagluglythrphethrseraspvalsersertyrleuglugly151015glnalaalalysglupheilealatrpleuvallysalaarg202530<210>3<211>659<212>prt<213>智人<400>3metalaalavalileleugluserilepheleulysargserglngln151015lyslyslysthrserproleuasnphelyslysargleupheleuleu202530thrvalhislysleusertyrtyrglutyraspphegluargglyarg354045argglyserlyslysglyserileaspvalglulysilethrcysval505560gluthrvalvalproglulysasnproproprogluargglnilepro65707580argargglyglugluserserglumetgluglnileserileileglu859095argpheprotyrpropheglnvalvaltyraspgluglyproleutyr100105110valpheserprothrglugluleuarglysargtrpilehisglnleu115120125lysasnvalileargtyrasnseraspleuvalglnlystyrhispro130135140cysphetrpileaspglyglntyrleucyscysserglnthralalys145150155160asnalametglycysglnileleugluasnargasnglyserleulys165170175proglyserserhisarglysthrlyslysproleuproprothrpro180185190glugluaspglnileleulyslysproleuproprogluproalaala195200205alaprovalserthrsergluleulyslysvalvalalaleutyrasp210215220tyrmetprometasnalaasnaspleuglnleuarglysglyaspglu225230235240tyrpheileleuglugluserasnleuprotrptrpargalaargasp245250255lysasnglyglngluglytyrileproserasntyrvalthrgluala260265270gluaspserileglumettyrglutrptyrserlyshismetthrarg275280285serglnalagluglnleuleulysglngluglylysgluglyglyphe290295300ilevalargaspserserlysalaglylystyrthrvalservalphe305310315320alalysserthrglyaspproglnglyvalilearghistyrvalval325330335cysserthrproglnserglntyrtyrleualaglulyshisleuphe340345350serthrileprogluleuileasntyrhisglnhisasnseralagly355360365leuileserargleulystyrprovalserglnglnasnlysasnala370375380proserthralaglyleuglytyrglysertrpgluileaspprolys385390395400aspleuthrpheleulysgluleuglythrglyglnpheglyvalval405410415lystyrglylystrpargglyglntyraspvalalailelysmetile420425430lysgluglysermetsergluaspglupheilegluglualalysval435440445metmetasnleuserhisglulysleuvalglnleutyrglyvalcys450455460thrlysglnargproilepheileilethrglutyrmetalaasngly465470475480cysleuleuasntyrleuargglumetarghisargpheglnthrgln485490495glnleuleuglumetcyslysaspvalcysglualametglutyrleu500505510gluserlysglnpheleuhisargaspleualaalaargasncysleu515520525valasnaspglnglyvalvallysvalserasppheglyleuserarg530535540tyrvalleuaspaspglutyrthrserservalglyserlysphepro545550555560valargtrpserproprogluvalleumettyrserlyspheserser565570575lysseraspiletrpalapheglyvalleumettrpgluiletyrser580585590leuglylysmetprotyrgluargphethrasnsergluthralaglu595600605hisilealaglnglyleuargleutyrargprohisleualaserglu610615620lysvaltyrthrilemettyrsercystrphisglulysalaaspglu625630635640argprothrphelysileleuleuserasnileleuaspvalmetasp645650655glugluser当前第1页12