[发明专利]基于泊松分布模型的蛋白质二级质谱鉴定方法

权利要求书

1.一种基于泊松分布模型的蛋白质二级质谱鉴定方法，其特征在于，包括如下步骤：

(1)虚拟酶解蛋白质数据库序列，并根据肽段的质量数对酶解后的肽段建立肽段数据库和肽段数据库索引；

(2)根据待分析实验图谱中母离子的核质比在步骤(1)所述的肽段数据库中找出符合要求的候选肽段，将选出符合要求的候选肽段作为理论图谱；

(3)对待分析实验图谱进行去同位素峰和去噪处理；

(4)将步骤(3)中的待分析实验图谱和步骤(2)中每张候选肽段的理论图谱进行匹配打分，选择得分最高的候选肽段作为本次实验图谱的鉴定结果；

所述步骤(4)将待分析实验图谱和理论图谱进行匹配打分包括：基于泊松分布的连续匹配打分，基于泊松分布的b、y离子的匹配打分，具体如下：

1)基于泊松分布的连续匹配打分：

$S_1=-\lg \left(\frac{{\left(0.1789k_1\right)}^{(K_1+\[0.1789k_1\]+1)}\·l^{-0.1789k_1}}{(K_1+\[0.1789k_1\]+1)!}\right)+\lg \left(\frac{{\left(0.1789k_1\right)}^{(\[0.1789k_1\]+1)}\·l^{-0.1789k_1}}{(\[0.1789k_1\]+1)!}\right)$

其中k₁是实验图谱匹配数，K₁是实验图谱连续匹配上的个数，S₁是连续匹配的得分，0.1789是随机匹配的概率值，等于实际连续错误匹配数除以理论连续错误匹配数；

2)基于泊松分布的b、y离子的匹配打分：

$S_2=-\lg \left(\frac{{\left(0.2110k_1\right)}^{(K_2+\[0.2110k_1\]+1)}\·l^{-0.2110k_1}}{(K_2+\[0.2110k_1\]+1)!}\right)+\lg \left(\frac{{\left(0.2110k_1\right)}^{(\[0.2110k_1\]+1)}\·l^{-0.2110k_1}}{(\[0.2110k_1\]+1)!}\right)$

其中K₂是实验图谱和理论图谱中b/y离子的匹配个数，S₂是b/y匹配的得分，0.2110是随机匹配的概率，等于实际错误匹配数除以理论错误匹配数；

3)基于泊松分布模型的总的打分：

$S\_Pep=(S_1+S_2)*\frac{1+\underset{i=1}{\overset{K}{\Σ}}{\mathrm{Pr}}_i}{1+K\×0.155}$

(5)针对所有实验的鉴定结果进行整体的假阳性控制。