專利名稱：用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種蛋白質(zhì)組分析方法，具體地說，涉及一種預(yù)測肽序列碎裂后產(chǎn)生的碎片離子的分子式的方法。
背景技術(shù)：
在目前利用肽指紋質(zhì)譜及串聯(lián)質(zhì)譜技術(shù)和數(shù)據(jù)庫搜索及直接解序(de novo)方法鑒定肽序列和蛋白質(zhì)的研究中，質(zhì)譜數(shù)據(jù)的預(yù)處理以及鑒定結(jié)果的后處理非常重要。
被鑒定的多肽在質(zhì)譜儀中被碎裂為碎片離子，這些碎片離子的質(zhì)量和豐度被質(zhì)譜儀器測量出來，形成串聯(lián)質(zhì)譜。每一個碎片離子以及其同位素離子都在串聯(lián)質(zhì)譜中形成對應(yīng)的譜峰?？紤]到碎片離子的同位素峰會給肽或蛋白質(zhì)的鑒定過程造成混淆，比如某些氨基酸殘基之間的質(zhì)量差約為0.34，1和1.5da，而同一個碎片離子的一價，二價，三價的同位素峰之間的質(zhì)荷比(m/z)差分別為1、0.5和0.333，這些氨基酸殘基質(zhì)量差值與同位素峰的m/z差值重疊，導(dǎo)致在鑒定過程中需要判斷串聯(lián)質(zhì)譜中的一個譜峰是某個碎片離子峰還是另一個碎片離子的同位素峰；此外，多個氨基酸質(zhì)量求和后與某個碎片離子的同位素峰的重疊現(xiàn)象會更多。因此，傳統(tǒng)的數(shù)據(jù)預(yù)處理任務(wù)之一是識別出一個碎片離子的同位素峰并予以剔除。
然而，事實上，質(zhì)譜中表現(xiàn)出的碎片離子的同位素峰的分布模式與該碎片離子的原子組成(即分子式)是密切相關(guān)的。因此就需要有一種方法能夠利用碎片離子的同位素峰來預(yù)測該碎片離子的分子式，這樣，預(yù)測出的碎片離子的分子式一方面可以為肽鑒定的數(shù)據(jù)庫搜索及de novo方法提供更多更準(zhǔn)確的信息，另一方面，為鑒定結(jié)果進(jìn)行后處理提供更多的依據(jù)。

發(fā)明內(nèi)容
本發(fā)明的目的在于提供一種利用串聯(lián)質(zhì)譜中的碎片離子的同位素峰來預(yù)測該碎片離子的分子式的方法。
為了實現(xiàn)上述目的，本發(fā)明提供一種用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，包括步驟1)從串聯(lián)質(zhì)譜中獲取一碎片離子的單同位素及其至少一個同位素的譜峰，計算所述碎片離子的單同位素的質(zhì)量、所述碎片離子的單同位素的譜峰和所述碎片離子的至少一個同位素的譜峰之間的相對豐度；步驟2)提供碎片離子的一通用分子式，所述通用分子式中各元素的原子個數(shù)待定；步驟3)用所述通用分子式得到碎片離子的理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素的相對豐度；所述理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素離子的相對豐度為所述通用分子式中待定的原子個數(shù)的函數(shù)；步驟4)將步驟3)中得到的質(zhì)量和相對豐度與步驟1)中從串聯(lián)質(zhì)譜質(zhì)量和相對豐度做匹配，以獲得所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解，從而得到所述碎片離子的分子式。
在上述技術(shù)方案中，步驟1)和步驟3)中所述的碎片離子的至少一個同位素包括碎片離子的第一同位素和第二同位素。
在上述技術(shù)方案中，將步驟1)中得到的所述碎片離子的單同位素的質(zhì)量、所述碎片離子的單同位素的譜峰和所述碎片離子的至少一個同位素的譜峰之間的相對豐度構(gòu)成一實驗的同位素分布向量；將步驟3)中得到的碎片離子的理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素的相對豐度構(gòu)成一理論的同位素分布向量；步驟4)中的所述匹配是用所述實驗的同位素分布向量與所述的理論的同位素分布向量之間的歐氏距離作為匹配分?jǐn)?shù)。
在上述技術(shù)方案中，還包括用使獲得的分子式符合化學(xué)意義的化學(xué)規(guī)則約束條件約束所述匹配。
在上述技術(shù)方案中，通過所述匹配獲得的所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解包括通過所述匹配獲得所述通用分子式中待定的各元素的原子個數(shù)的實數(shù)解；在所述實數(shù)解的領(lǐng)域內(nèi)搜索得到所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解。
在上述技術(shù)方案中，還包括對步驟4)中得到的所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解進(jìn)行過濾的步驟。所述過濾包括平均同位素分布模式方法，該方法用碎片離子的理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素的相對豐度之間的統(tǒng)計關(guān)系過濾所述非負(fù)整數(shù)解。所述過濾包括用使獲得的分子式符合化學(xué)意義的化學(xué)規(guī)則約束條件過濾所述非負(fù)整數(shù)解。所述過濾包括用兩個碎片離子的非負(fù)整數(shù)解進(jìn)行交叉驗證以過濾所述兩個碎片離子的非負(fù)整數(shù)解。
本發(fā)明的優(yōu)點在于1)本方法是對串聯(lián)質(zhì)譜中碎片離子的同位素譜峰信息的充分利用；2)本方法能通過串聯(lián)質(zhì)譜碎片離子的同位素譜峰的模式，快速準(zhǔn)確地計算此碎片離子對應(yīng)的分子式(準(zhǔn)確程度與質(zhì)譜的精度相關(guān)，精度越高，計算出的分子式越可靠)；3)本方法可以提供碎片離子的準(zhǔn)確的分子式信息，可對鑒定多肽序列的數(shù)據(jù)庫搜索方法提供的候選序列進(jìn)行鑒別；4)本方法計算出的離子分子式可以指導(dǎo)求解多肽序列的de novo方法產(chǎn)生高可靠的候選的序列。
具體實施例方式
下面結(jié)合附圖
和具體實施方式
對本發(fā)明作進(jìn)一步詳細(xì)描述。
將一個碎片離子的單同位素記為P，此碎片離子的第一同位素記為P1，第二同位素碎片離子記為P2，依此類推，第N同位素離子記為PN。在這里，碎片離子單同位素P是指在該離子的各種組成元素均為單同位素(即質(zhì)子數(shù)和中子數(shù)相同)。而碎片離子的同位素是指與單同位素碎片離子具有相同的分子式、但是比單同位素帶有更多額外中子的離子，例如碎片離子的第一同位素P1比碎片離子的單同位素P多帶有一個額外的中子，第二同位素P2比單同位素P多帶有兩個額外的中子，依此類推。在本發(fā)明中，碎片離子的同位素是在整體上比碎片離子的單同位素帶有額外的中子的離子。
肽序列進(jìn)入質(zhì)譜儀被離子化，且在質(zhì)譜儀中，具有特定質(zhì)荷比(m/z)的肽離子(這些肽離子通常也有相同的氨基酸序列)在碰撞-誘導(dǎo)的分離(Collision-InducedDissociation，CID)作用下裂解為多個碎片離子。這些碎片離子的m/z被檢測量出來從而形成串聯(lián)質(zhì)譜，在一個串聯(lián)質(zhì)譜中，其橫坐標(biāo)表示碎片離子的質(zhì)荷比(m/z)，其縱坐標(biāo)為檢測到的碎片離子的豐度。
在串聯(lián)質(zhì)譜中，挑選出一個碎片離子的單同位素P以及其同位素P1～PN中至少一個對應(yīng)的譜峰，本發(fā)明的目標(biāo)則是通過這些同位素峰的分布情況來預(yù)測碎片離子的單同位素P對應(yīng)的分子式。在本發(fā)明的一個實施例中，從串聯(lián)質(zhì)譜中僅挑選出該碎片離子的單同位素P以及其第一同位素P1和第二同位素P2。從后面的描述中本領(lǐng)域的技術(shù)人員很容易理解，在本發(fā)明的其它實施例中，對于同位素碎片離子，也可以僅挑選出碎片離子的一個同位素的譜峰——例如第一同位素碎片離子P1，或者也可以挑選出更多的同位素的譜峰，不同數(shù)目同位素碎片離子的選取都可以實現(xiàn)本發(fā)明的方法，但是會影響到本發(fā)明實施時計算的復(fù)雜度和精度。
從串聯(lián)質(zhì)譜中還可以得到單同位素碎片離子P的離子質(zhì)量Me，這是本領(lǐng)域的技術(shù)人員所熟知的。
為了方便于下面的計算，首先定義一實驗的同位素分布向量eIPV＝(Me，I1，I2)，其中，Me為從串聯(lián)質(zhì)譜中獲得的碎片離子的單同位素P的離子質(zhì)量，I1和I2分別對應(yīng)碎片離子的第一同位素P1和第二同位素P2的譜峰相對于單同位素P的譜峰的相對豐度，這些數(shù)據(jù)均可從串聯(lián)質(zhì)譜中獲得。
然后，再定義一理論的同位素分布向量tIPV＝(M，T1，T2)，該理論同位素分布向量tIPV可從碎片離子的通用分子式獲得。設(shè)碎片離子的通用分子式為Cn1Hn2Nn3On4Sn5，其中該分子式中表示各原子組成個數(shù)的n1～n5為待定參數(shù)。這樣，在理論同位素分布向量tIPV中，M為從通用分子式獲得的碎片離子的質(zhì)量，T1和T2分別為從通用分子式獲得的第一同位素碎片離子和第二同位素碎片離子關(guān)于單同位素碎片離子相對豐度。理論同位素分布向量tIPV可具體可通過公式得到M＝V×X (1)T1＝n1qC+n2qH+n3qN+n4qO1+n5qS1(2)T2=n4qO2+n5qS2+12T12-12(n1qC2+n2qH2+n3qN2+n4qO12+n5qS12)---(3)]]>其中V＝[12，1，14，16，32]，V中的數(shù)字為各元素的原子量，X＝[n1，n2，n3，n4，n5]T；qC、qH和qN分別是自然界中13C相對于12C、D相對于H、14N相對于15N的相對豐度，q01和q02則分別是自然界中17O相對于18O、18O相對于16O的相對豐度，qs1和qs2是自然界中33S相對于32S、34S相對于32S的相對豐度，這些相對豐度均為已知數(shù)值。
可見，對于理論的同位素分布向量tIPV＝(M，T1，T2)，其中的M、T1和T2均為X＝[n1，n2，n3，n4，n5]的函數(shù)。
在本發(fā)明中，將理論的同位素分布向量tIPV＝(M，T1，T2)與實驗同位素分布向量eIPV＝(Me，I1，I2)做匹配，以便獲得與實驗的同位素分布向量最匹配的分子式，也即通用分子式中的原子組成向量X＝[n1，n2，n3，n4，n5]的一個非負(fù)整數(shù)解。
在本發(fā)明的一個實施例中，用理論的同位素分布向量tIPV和實驗的同位素分布向量eIPV之間的歐氏距離E作為tIPV與eIPV的匹配分?jǐn)?shù)
E=δm2+δ12+δ22=(M-Me)2+(T1-I1)2+(T2-I2)2---(4)]]>將公式(1)～(3)代入(4)，得到δm＝n1*12+n2*1+n3*14+n4*16+n5*32-Me1，(5)δ1＝n1*qC+n2*qH+n3*qN+n4*qO1+n5*qs1-I1，(6)δ2=n4*qO2+n5*qS2-12(n1*qC2+n2*qH2+n3*qN2+n4*qO12+n5*qS12)]]>+(n1*qC+n2*qH+n3*qN+n4*qO+n5*qS)*I1-12I12-I2+12δ12.---(7)]]>忽略公式(7)中的 項，則有[δmδ1δ2]＝AX+B，得到Q(X)=E2=δmδ1δ2δmδ1δ2=XTATAX+2BTAX+BTB,---(8)]]>則有E=Q(X)=XTATAX+2BTAX+BTB---(9)]]>這里在公式(9)中，X＝[n1，n2，n3，n4，n5]T是待定的碎片離子的原子組成向量，A和B是由已知量構(gòu)成的常數(shù)矩陣，這里已知量包括從串聯(lián)質(zhì)譜中獲得的Me、I1和I2，和V＝[12，1，14，16，32]以及公式(2)和(3)中的各同位素的相對豐度。
將公式(9)所描述的歐氏距離E最小化，即可得到X的一個解。通常，為了使獲得的分子式符合化學(xué)意義，優(yōu)選還要對公式(9)設(shè)置一些化學(xué)規(guī)則約束條件，例如●用X獲得的分子式對應(yīng)的碎片離子質(zhì)量一定要在范圍[Me-δ，Me+δ]內(nèi)，δ是m/z誤差的最大范圍，δ可由質(zhì)譜儀的測量精度來確定。也就是要滿足|VX-Me|≤δ。
●對于碎片離子分子式中的某種元素，用離子的m/z除以這種元素質(zhì)量最低的同位素的質(zhì)量數(shù)，取所得結(jié)果的整數(shù)部分就是此元素個數(shù)的上限。例如元素O的原子量為16，若離子的質(zhì)荷比為m/z，則碎片離子中O元素的個數(shù)的上限為 即在X中 類似地，對于碎片離子中其它元素也可獲得相似的約束條件。
●在碎片離子中，C的個數(shù)一定小于H的個數(shù)(即在X中n1＜n3)、O和N的個數(shù)一定小于C的個數(shù)(即在n4＜n1和n3＜n1)等等。這些約束條件隱含在氨基酸殘基的分子組成方式和主要離子類型的組成方式中，本領(lǐng)域的技術(shù)人員很容易根據(jù)它們的特點總結(jié)出來。
●在帶一個電子的離子中，H和N的個數(shù)之和為奇數(shù)。原因是如果離子帶有一個電荷，那么就有一個不飽和化學(xué)鍵存在，并且，H和N都有奇數(shù)個化合價而C、O、S都有偶數(shù)個化合價。
應(yīng)當(dāng)理解，本領(lǐng)域的技術(shù)人員也可從使碎片離子的分子式符合化學(xué)意義的目的出發(fā)構(gòu)造出其它的約束條件。
上述約束條件或者其它約束條件中的一部分或者全部可表示為一個線性不等式DX≤G。這樣，結(jié)合公式(9)，可以通過標(biāo)準(zhǔn)的二次規(guī)劃方法來解決歐氏距離E的這個最小化問題，如公式(10)所示 從公式(10)用二次規(guī)劃方法求出的X的最優(yōu)解為一個實數(shù)域內(nèi)的解XR，為了尋找真正的分子式，可以將XR當(dāng)作起始點，然后在它的鄰域內(nèi)局部搜索X的非負(fù)整數(shù)候選解。確切地說，就是對每一個與XR存在一個距離d的范圍內(nèi)的非負(fù)整數(shù)候選候選解分子式進(jìn)行打分，或者說用公式(9)評價這些非負(fù)整數(shù)候選解的匹配度。d的值是與離子質(zhì)量范圍相適應(yīng)的。這樣避免了枚舉所有可能的分子式，能夠在大質(zhì)量范圍內(nèi)預(yù)測離子分子式并且確保較高的可靠性和運(yùn)行效率。
經(jīng)過局部搜索，仍會產(chǎn)生一定數(shù)量的候選分子式，其中包括一些不合法的和與實驗串聯(lián)質(zhì)譜不匹配的分子式(可分別稱為無效的和不可能的分子式)，為了提高預(yù)測的精確度，優(yōu)選需要盡可能多的排除它們。在本發(fā)明中可利用包括平均同位素分布模式、化學(xué)規(guī)則約束和交叉驗證中的一種或者多種方法來過濾候選分子式。這些方法具體描述如下A.平均的同位素分布模式所說的平均的同位素分布模式是理論同位素分布向量tIPV＝(M，T1，T2)中的組成部分M、T1和T2之間的統(tǒng)計關(guān)系。為了尋找碎片離子的理論平均同位素分布模式，發(fā)明人計算了現(xiàn)有的蛋白質(zhì)數(shù)據(jù)庫中所有蛋白質(zhì)的trypsin水解對應(yīng)的多肽的理論碎片離子的同位素的平均分布和標(biāo)準(zhǔn)差，揭示了tIPV的組成部分M、T1和T2之間的關(guān)系。具體地說，發(fā)明人首先將SWISS-PROT中的蛋白質(zhì)進(jìn)行理論酶切計算得到多肽；然后選擇質(zhì)量在(60u～3000u)內(nèi)的多肽，這個范圍對應(yīng)著Q-TOF MS/MS實驗質(zhì)譜的標(biāo)準(zhǔn)范圍。另外，值得注意的是S的同位素+2S在自然界中的含量很高(出現(xiàn)的幾率是0.04210，大約是18O的20倍)，而多數(shù)情況下能夠包含五個以上的S的多肽十分少見。因此，我們可以將上述分子式分成六類S0，S1，S2，S3，S4和S5+，分別對應(yīng)所含S的個數(shù)為0，1，2，3，4和5個及5個以上的肽段。發(fā)明人按這六個類別對做了統(tǒng)計。統(tǒng)計結(jié)果顯示T1與質(zhì)量M呈線性關(guān)系，T2則與M呈二次關(guān)系，而T2隨著T1增加而增加并且與T1成二次函數(shù)關(guān)系。
這樣，通過T1、T2與M的上述分布關(guān)系可以對候選分子式進(jìn)行過濾，以排除那些無效的和/或不可能的分子式。
B.化學(xué)規(guī)則約束這里的化學(xué)規(guī)則約束與公式(10)中的約束條件DX≤G相類似，其區(qū)別在于在公式(10)中，約束條件DX≤G用于約束公式E=Q(X)=XTATAX+2BTAX+BTB]]>以便得到在此約束條件下X的一個實數(shù)域內(nèi)的解XR。而在這里，這些約束條件用于約束在XR的領(lǐng)域內(nèi)搜索得到的非負(fù)整數(shù)解候選分子式，以便對這些候選分子式進(jìn)行過濾。
C.交叉驗證特別地，一個肽段的b系列的碎片離子都是同源的，包括b-，a-，b*-，a*-，b°-，a°型離子，它們共享一個相同的原始氨基酸序列，由此可推測它們的同位素分布模式很相似。y系列離子也是這樣。如果質(zhì)譜中某兩個碎片離子的Me相差28、17或18，并且這兩個碎片離子的I1和I2很接近，就可認(rèn)為兩個碎片離子對應(yīng)的eIPV同源的。而后，我們就可以使用同源的eIPV對預(yù)測結(jié)果進(jìn)行交叉驗證。例如對于同源的兩個碎片離子，在一個碎片離子中的候選分子式列表中有Ca1Ha2Na3Oa4Sa5，如果Ca1-1Ha2Na3Oa4-1Sa5沒有出現(xiàn)在另一個碎片離子的候選分子式列表里，那么就可以認(rèn)為候選分子式Ca1Ha2Na3Oa4Sa5是隨機(jī)匹配上的結(jié)果而將它排除。
權(quán)利要求
1.一種用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，包括步驟1)從串聯(lián)質(zhì)譜中獲取一碎片離子的單同位素及其至少一個同位素的譜峰，計算所述碎片離子的單同位素的質(zhì)量、所述碎片離子的單同位素的譜峰和所述碎片離子的至少一個同位素的譜峰之間的相對豐度；步驟2)提供碎片離子的一通用分子式，所述通用分子式中各元素的原子個數(shù)待定；步驟3)用所述通用分子式得到碎片離子的理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素的相對豐度；所述理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素離子的相對豐度為所述通用分子式中待定的原子個數(shù)的函數(shù)；步驟4)將步驟3)中得到的質(zhì)量和相對豐度與步驟1)中從串聯(lián)質(zhì)譜質(zhì)量和相對豐度做匹配，以獲得所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解，從而得到所述碎片離子的分子式。
2.根據(jù)權(quán)利要求1所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，步驟1)和步驟3)中所述的碎片離子的至少一個同位素包括碎片離子的第一同位素和第二同位素。
3.根據(jù)權(quán)利要求1所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，將步驟1)中得到的所述碎片離子的單同位素的質(zhì)量、所述碎片離子的單同位素的譜峰和所述碎片離子的至少一個同位素的譜峰之間的相對豐度構(gòu)成一實驗的同位素分布向量；將步驟3)中得到的碎片離子的理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素的相對豐度構(gòu)成一理論的同位素分布向量；步驟4)中的所述匹配是用所述實驗的同位素分布向量與所述的理論的同位素分布向量之間的歐氏距離作為匹配分?jǐn)?shù)。
4.根據(jù)權(quán)利要求1所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，還包括用使獲得的分子式符合化學(xué)意義的化學(xué)規(guī)則約束條件約束所述匹配。
5.根據(jù)權(quán)利要求1所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，通過所述匹配獲得的所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解包括通過所述匹配獲得所述通用分子式中待定的各元素的原子個數(shù)的實數(shù)解；在所述實數(shù)解的領(lǐng)域內(nèi)搜索得到所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解。
6.根據(jù)權(quán)利要求1所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，還包括對步驟4)中得到的所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解進(jìn)行過濾的步驟。
7.根據(jù)權(quán)利要求6所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，所述過濾包括平均同位素分布模式方法，該方法用碎片離子的理論上的單同位素的質(zhì)量、碎片離子的單同位素和其至少一個同位素的相對豐度之間的統(tǒng)計關(guān)系過濾所述非負(fù)整數(shù)解。
8.根據(jù)權(quán)利要求6所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，所述過濾包括用使獲得的分子式符合化學(xué)意義的化學(xué)規(guī)則約束條件過濾所述非負(fù)整數(shù)解。
9.根據(jù)權(quán)利要求6所述的用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，其特征在于，所述過濾包括用兩個碎片離子的非負(fù)整數(shù)解進(jìn)行交叉驗證以過濾所述兩個碎片離子的非負(fù)整數(shù)解。
全文摘要
本發(fā)明公開了一種用串聯(lián)質(zhì)譜中碎片離子的同位素峰預(yù)測離子分子式的方法，該方法從串聯(lián)質(zhì)譜和從各元素的原子個數(shù)待定的通用分子式中分別獲取碎片離子的單同位素的質(zhì)量以及各同位素譜峰相對于單同位素的相對豐度；將分別獲取的質(zhì)量和相對豐度做匹配以獲得所述通用分子式中待定的各元素的原子個數(shù)的非負(fù)整數(shù)解，得到碎片離子的分子式。本發(fā)明的方法利用串聯(lián)質(zhì)譜中碎片離子的同位素譜峰信息，通過串聯(lián)質(zhì)譜碎片離子的同位素譜峰的模式計算此碎片離子對應(yīng)的分子式。本發(fā)明的方法可以提供碎片離子準(zhǔn)確的分子式信息，可對鑒定多肽序列的數(shù)據(jù)庫搜索方法提供的候選序列進(jìn)行鑒別；以及為求解多肽序列的de novo方法產(chǎn)生高可靠候選序列提供依據(jù)。
文檔編號G06F19/00GK1773276SQ20041009080
公開日2006年5月17日 申請日期2004年11月12日 優(yōu)先權(quán)日2004年11月12日
發(fā)明者高文, 張京芬, 蔡津津, 賀思敏, 曾嶸, 陳潤生, 王海鵬 申請人:中國科學(xué)院計算技術(shù)研究所