相關申請案的交叉參考

本申請案主張2015年2月5日申請的序列號為62/112,212的美國臨時專利申請案的權益，所述臨時專利申請案的全部內容以引用的方式并入本文中。

背景技術：

質譜儀通常與色譜或其它分離系統(tǒng)耦合，以便識別及特征化來自樣本的所關注的洗脫已知化合物。在此耦合系統(tǒng)中，離子化洗脫溶劑，并以規(guī)定的時間間隔從洗脫溶劑獲得一系列質譜。這些時間間隔從例如1秒到100分鐘或更長。所述一系列質譜圖形成色譜或提取離子色譜(xic)。

在xic中發(fā)現(xiàn)的峰值用于識別或特征化樣本中的已知化合物。然而，在復雜混合物中，具有相同質荷比(m/z)的其它峰值的干擾可能使得難以確定表示已知化合物的峰值。在某些情況下，沒有關于已知化合物的預期保留時間的信息。在其它情況下，可能知曉已知化合物的近似保留時間。然而，即使在后一種情況下，如果樣本是復雜的，或如果樣本之間存在不止是少量的保留時間變化，那么已知化合物的確切峰值可能是不明確的。因此，在這些情況下，通常難以識別或特征化已知化合物。

在傳統(tǒng)的分離耦合質譜系統(tǒng)中，選擇已知化合物的碎片或產(chǎn)物離子進行分析。然后在針對包含產(chǎn)物離子的質量范圍的分離的每一間隔執(zhí)行質譜/質譜(ms/ms)掃描。舉例來說，隨著時間推移，收集在每一ms/ms掃描中發(fā)現(xiàn)的產(chǎn)物離子的強度，并將其作為光譜集合或xic來分析。

舉例來說，對于簡單的樣本混合物，通常在已知化合物的預期保留時間在xic中發(fā)現(xiàn)代表產(chǎn)物離子的單峰值。然而，對于更復雜的混合物，除了已知化合物的預期保留時間之外，代表產(chǎn)物離子的兩個或更多個峰值位于光譜集合中的一或多個額外時間間隔處。換句話說，針對產(chǎn)物離子的xic可具有兩個或更多個峰值。

在更復雜的混合物中識別所關注的化合物的一種傳統(tǒng)方法是對已知化合物的兩種或更多種產(chǎn)物離子在其處具有峰值的時間間隔進行定位。舉例來說，此方法在已知序列的肽被定量時用于蛋白質組學。

在典型的多反應監(jiān)測(mrm)方法中，監(jiān)測兩個或更多個mrm轉換，每一者對應于肽的不同產(chǎn)物離子轉換。如果先前的發(fā)現(xiàn)數(shù)據(jù)可用，那么這些轉換基于數(shù)據(jù)中觀察到的最大產(chǎn)物離子。否則，舉例來說，這些轉換基于預測的y離子。針對這兩個或更多個mrm轉換來分析xic。對于所有轉換，在其間都存在產(chǎn)物離子峰值的時間用于特征化已知化合物。

對于復雜樣本，特別是如果未準確知道預期的保留時間，那么產(chǎn)物離子譜的收集可能存在模糊性。舉例來說，可存在多于一個的保留時間或時間間隔，兩個或更多個mrm轉換中的每一者在其間存在產(chǎn)物離子峰值。

幾乎沒有額外信息可用于解決復雜樣本引入的模糊性。在傳統(tǒng)的分離耦合質譜系統(tǒng)中，通常使用窄的前體離子質量窗口寬度來執(zhí)行在每一時間間隔的每一產(chǎn)物離子的每一ms/ms掃描。因此，在數(shù)據(jù)采集后可用的針對每一碎片離子的特定時間間隔的產(chǎn)物離子質譜幾乎不能提供額外的理解。

標題為“將窗口化質譜數(shù)據(jù)用于保留時間確定或確認(useofwindowedmassspectrometrydataforretentiontimedeterminationorconfirmation)”的第14/368,874號美國專利申請案(以下稱為“'874申請案”)描述收集額外ms/ms數(shù)據(jù)并使用此數(shù)據(jù)來解決由復雜樣本引入的模糊性的方法。在'874申請案中，使用分離耦合質譜系統(tǒng)，所述系統(tǒng)使用一或多個循序質量窗口寬度在每一時間間隔執(zhí)行ms/ms掃描，以便橫跨整個質量范圍。換句話說，在分離過程中，在每一時間間隔獲得整個質量范圍的光譜信息。一種用于使用一或多個循序質量窗口寬度來執(zhí)行ms/ms掃描以便橫跨整個質量的方法是absciex的swathtm技術。

高分辨率及高通量儀器允許使用具有相鄰或重疊質量窗口寬度的多次掃描在時間間隔內準確掃描質量范圍。將多次掃描的結果拼接在一起，以在每一時間間隔產(chǎn)生整個質量范圍的光譜。然后可針對分離過程的任何時間間隔上的質量范圍中的任何質量計算xic。

在'874申請案中，使用窗口化獲取方法收集的整個質量范圍的光譜信息用于解決復雜混合物中的保留時間模糊性。換句話說，當在分離過程中的兩個或更多個不同的時間間隔處發(fā)現(xiàn)產(chǎn)物離子在光譜集合中具有兩個或更多個峰值時，對不同時間間隔中的每一者處的整個質量范圍的產(chǎn)物離子質譜進行分析以確定實際保留時間。使用各種準則來分析整個質量范圍的質譜，所述準則包括電荷狀態(tài)、同位素狀態(tài)，質量準確度以及與已知化合物的已知碎片化分布相關聯(lián)的一或多個質量差異?；谶@些準則，對兩個或更多個時間間隔處的產(chǎn)物離子的每一峰值進行評分。已知化合物的保留時間經(jīng)識別處于具有最高得分的峰值處。

'874申請案的方法也可使用已知化合物的多于一種的產(chǎn)物離子來解決保留時間模糊性。兩種或更多種產(chǎn)物離子中的每一者的峰值在兩個或更多個時間間隔處被獨立評分，并且兩種或更多種產(chǎn)物離子的峰值的得分在兩個或更多個時間間隔中的每一者處組合。然后根據(jù)兩個或更多個時間間隔中的每一者處的組合得分來確定保留時間。換句話說，'874申請案的方法在每一時間間隔分組已知化合物的兩種或更多種不同產(chǎn)物離子的峰值，并在每一時間間隔比較群組的組合得分。

然而，'874申請案的方法沒有描述如何選擇已知化合物的兩種或更多種不同產(chǎn)物離子的xic峰值。如上文關于mrm轉換所描述，如果先前的發(fā)現(xiàn)數(shù)據(jù)可用，那么所選擇的mrm轉換是基于在數(shù)據(jù)中觀察到的最大產(chǎn)物離子。換句話說，選擇具有最強產(chǎn)物離子峰值的mrm轉換。

舉例來說，假設'874申請案的方法類似地通過首先選擇第一產(chǎn)物離子的最強峰值并且接著選擇在時間上最接近于第一產(chǎn)物離子的最強峰值的頂點的其它產(chǎn)物離子的峰值來對峰值分組?，F(xiàn)在假設第一產(chǎn)物離子的最強峰值來自并非是已知化合物的前體離子。然后對其它產(chǎn)物離子的峰值進行分組并用錯誤峰值進行評分。因此，'874申請案的方法的結果高度依賴于峰值群組的適當峰值選擇。

目前，質譜行業(yè)缺乏從大的質量范圍上所收集的多種產(chǎn)物離子的xic中選擇已知化合物的峰值群組的可靠且準確的方法。

技術實現(xiàn)要素：

揭示一種用于分組已知化合物的兩種或更多種產(chǎn)物離子的提取離子色譜(xic)峰值的系統(tǒng)。所述系統(tǒng)包含分離裝置、質譜儀及處理器。

所述分離裝置將已知化合物從樣本混合物分離。所述質譜儀在多個保留時間的每一保留時間使用一或多個循序質量窗口寬度對分離樣本混合物執(zhí)行一或多個質譜/質譜(ms/ms)掃描，以便橫跨整個質量范圍，從而產(chǎn)生針對所述多個保留時間的所述整個質量范圍的產(chǎn)物離子光譜集合。

所述處理器接收針對所述多個保留時間的所述整個質量范圍的所述產(chǎn)品離子光譜集合，并選擇所述已知化合物的m種產(chǎn)物離子。所述處理器進一步計算來自所述產(chǎn)物離子光譜集合的所述m種產(chǎn)物離子中的每一者的xic，從而產(chǎn)生m個xic。所述處理器進一步將所述m個xic中的每一xic從其它m個xic中的每一者減去，從而產(chǎn)生個減法曲線。根據(jù)第一xic及第二xic計算每一減法曲線。在每一保留時間，使所述第一xic在所述每一保留時間的強度歸一化，且所述第一xic在兩個或更多個相鄰保留時間的兩個或更多個強度被歸一化。所述第二xic在所述每一保留時間的強度及所述第二xic在所述相鄰保留時間的兩個或更多個強度被歸一化。從所述第一xic的對應歸一化強度減去所述第二xic的歸一化強度。計算差值強度的統(tǒng)計測量。

所述處理器進一步識別一或多個保留時間的一或多個區(qū)域，其中所述個減法曲線的一或多個減法曲線具有在零的閾值內的值。針對所述一或多個區(qū)域中的每一區(qū)域，所述處理器獲得用以計算識別所述區(qū)域的所述一或多個減法曲線的兩個或更多個xic，并將所述區(qū)域中的所述兩個或更多個xic的每一峰值添加到峰值群組。

揭示一種用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的方法。在多個保留時間獲得整個質量范圍的產(chǎn)物離子光譜集合。使用分離裝置將已知化合物從樣本混合物分離。使用質譜儀在多個保留時間的每一保留時間使用一或多個循序前體離子質量窗口寬度對所述分離樣本混合物執(zhí)行一或多次ms/ms掃描，以便橫跨整個質量范圍，從而產(chǎn)生針對所述多個保留時間的所述整個質量范圍的所述產(chǎn)物離子光譜集合。使用處理器選擇所述已知化合物的m種產(chǎn)物離子。使用所述處理器計算來自所述產(chǎn)物離子光譜集合的所述m種產(chǎn)物離子中的每一者的xic，從而產(chǎn)生m個xic。

使用所述處理器將所述m個xic中的每一xic從其它m個xic中的每一者減去，從而產(chǎn)生個減法曲線。根據(jù)第一xic及第二xic計算每一減法曲線。在每一保留時間，所述第一xic在所述每一保留時間的強度及所述第一xic在兩個或更多個相鄰保留時間的兩個或更多個強度被歸一化。所述第二xic在所述每一保留時間的強度及所述第二xic在所述相鄰保留時間的兩個或更多個強度被歸一化。從所述第一xic的對應歸一化強度減去所述第二xic的歸一化強度。計算差值強度的統(tǒng)計測量。

使用所述處理器識別一或多個保留時間的一或多個區(qū)域，其中所述個減法曲線的一或多個減法曲線具有在零的閾值內的值。針對所述一或多個區(qū)域中的每一區(qū)域，使用所述處理器獲得用以計算識別所述區(qū)域的所述一或多個減法曲線的兩個或更多個xic，并將所述區(qū)域中的所述兩個或更多個xic的每一峰值添加到峰值群組。

揭示一種計算機程序產(chǎn)品，其包含非暫時性及有形計算機可讀存儲媒體，所述存儲媒體的內容包含具有指令的程序，所述指令在處理器上執(zhí)行以便執(zhí)行一種用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的方法。所述方法包含提供系統(tǒng)，其中所述系統(tǒng)包括一或多個相異軟件模塊，且其中所述相異軟件模塊包括測量模塊及分析模塊。

所述測量模塊在多個保留時間獲得整個質量范圍的產(chǎn)物離子光譜集合。使用分離裝置將已知化合物從樣本混合物分離。使用質譜儀在多個保留時間的每一保留時間使用一或多個循序前體離子質量窗口寬度對所述分離樣本混合物執(zhí)行一或多次ms/ms掃描，以便橫跨整個質量范圍，從而產(chǎn)生所述多個保留時間的所述整個質量范圍的產(chǎn)物離子光譜集合。

所述分析模塊選擇所述已知化合物的m種產(chǎn)物離子。所述分析模塊計算來自所述產(chǎn)物離子光譜集合的所述m種產(chǎn)物離子中的每一者的xic，從而產(chǎn)生m個xic。所述分析模塊將所述m個xic中的每一xic從其它m個xic中的每一者減去，從而產(chǎn)生個減法曲線。根據(jù)第一xic及第二xic計算每一減法曲線。在每一保留時間，所述第一xic在所述每一保留時間的強度及所述第一xic在兩個或更多個相鄰保留時間的兩個或更多個強度被歸一化。類似地，所述第二xic在所述每一保留時間的強度及所述第二xic在所述相鄰保留時間的兩個或更多個強度被歸一化。從所述第一xic的對應歸一化強度減去所述第二xic的歸一化強度。計算差值強度的統(tǒng)計測量。

所述分析模塊識別一或多個保留時間的一或多個區(qū)域，其中所述個減法曲線的一或多個減法曲線具有在零的閾值內的值。針對所述一或多個區(qū)域中的每一區(qū)域，所述分析模塊獲得用以計算識別所述區(qū)域的所述一或多個減法曲線的兩個或更多個xic，并將所述區(qū)域中的所述兩個或更多個xic的每一峰值添加到峰值群組。

本文闡述申請者的教示的這些及其它特征。

附圖說明

所屬領域的技術人員將理解，下文描述的圖式僅用于說明目的。圖式并不希望以任何方式限制本教示的范圍。

圖1是說明根據(jù)各種實施例的計算機系統(tǒng)的框圖。

圖2是根據(jù)各種實施例的已知化合物的五種產(chǎn)物離子的五個提取離子色譜(xic)的示范性繪圖。

圖3是根據(jù)各種實施例的在圖2中所示的保留時間50與65之間的五個xic的詳細部分的示范性繪圖。

圖4是展示根據(jù)各種實施例的從圖2中所示的xic中的兩個xic的強度的局部減法計算的平均值的減法曲線的示范性繪圖。

圖5是展示根據(jù)各種實施例的從圖2中所示的xic中的兩個xic的強度的局部減法計算的標準偏差值的減法曲線的示范性繪圖。

圖6是展示根據(jù)各種實施例的從圖2的五個xic計算的標準偏差值的個減法曲線的示范性繪圖。

圖7是展示根據(jù)各種實施例的在圖6中所示的保留時間50與65之間的標準偏差值的個減法曲線的詳細部分的示范性繪圖。

圖8是展示根據(jù)各種實施例的用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的系統(tǒng)的示意圖。

圖9是展示根據(jù)各種實施例的用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的方法的示范性流程圖。

圖10是根據(jù)各種實施例的包含執(zhí)行用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的方法的一或多個相異軟件模塊的系統(tǒng)的示意圖。

在詳細描述本教示的一或多個實施例之前，所屬領域的技術人員應了解，本教示在其應用中不限于下文詳細描述中闡述的或圖式中說明的構造的細節(jié)、組件的布置及步驟的布置。此外，應理解，本文所用的措辭及術語是為了描述的目的，而不應被視為限制。

具體實施方式

計算機實施的系統(tǒng)

圖1是說明可在其上實施本教示的實施例的計算機系統(tǒng)100的框圖。計算機系統(tǒng)100包含用于傳達信息的總線102或其它通信機構，以及與總線102耦合以用于處理信息的處理器104。計算機系統(tǒng)100還包含存儲器106，其可為隨機存取存儲器(ram)或其它動態(tài)存儲裝置，其耦合到總線102以用于存儲待由處理器104執(zhí)行的指令。存儲器106還可用于在執(zhí)行待由處理器104執(zhí)行的指令期間存儲臨時變量或其它中間信息。計算機系統(tǒng)100進一步包含只讀存儲器(rom)108或耦合到總線102以用于存儲用于處理器104的靜態(tài)信息及指令的其它靜態(tài)存儲裝置。存儲裝置110(例如磁盤或光盤)被提供并耦合到總線102以用于存儲信息及指令。

計算機系統(tǒng)100可經(jīng)由總線102耦合到顯示器112(例如陰極射線管(crt)或液晶顯示器(lcd))以用于向計算機用戶顯示信息。包含字母數(shù)字及其它鍵的輸入裝置114耦合到總線102以用于將信息及命令選擇傳達到處理器104。另一類型的用戶輸入裝置是光標控制件116(例如鼠標、軌跡球或光標方向鍵)，其用于將方向信息及命令選擇傳達到處理器104并用于控制顯示器112上的光標移動。此輸入裝置通常具有兩個軸(第一軸(即，x)及第二軸(即，y))中的兩個自由度，其允許裝置指定平面中的位置。

計算機系統(tǒng)100可執(zhí)行本教示。與本教示的某些實施方案一致，計算機系統(tǒng)100響應于處理器104執(zhí)行存儲器106中含有的一或多個指令的一或多個序列來提供結果。此類指令可從另一計算機可讀媒體(例如存儲裝置110)讀入存儲器106中。對存儲器106中含有的指令序列的執(zhí)行致使處理器104執(zhí)行本文描述的處理。替代地，可使用硬連線電路代替軟件指令或與軟件指令組合以實施本教示。因此，本教示的實施方案不限于硬件電路及軟件的任何特定組合。

本文所使用的術語“計算機可讀媒體”是指參與向處理器104提供指令以用于執(zhí)行的任何媒體。此類媒體可呈現(xiàn)許多形式，包括(但不限于)非易失性媒體、易失性媒體及傳輸媒體。非易失性媒體包含(例如)光盤或磁盤，例如存儲裝置110。易失性媒體包含動態(tài)存儲器，例如存儲器106。傳輸媒體包含同軸纜線、銅線及光纖，其包含包括總線102的線。

計算機可讀媒體的常見形式包含(例如)軟磁盤、軟性磁盤、硬盤、磁帶或任何其它磁性媒體、cd-rom、數(shù)字視頻盤(dvd)、藍光盤、任何其它光學媒體、拇指驅動器，存儲卡、ram，prom及eprom、flash-eprom、任何其它存儲器芯片或盒式磁盤，或計算機可從其讀取的任何其它有形媒體。

各種形式的計算機可讀媒體可涉及將一或多個指令的一或多個序列攜載到處理器104以用于執(zhí)行。舉例來說，指令最初可攜載于遠程計算機的磁盤上。遠程計算機可將指令加載到其動態(tài)存儲器中，并使用調制解調器通過電話線發(fā)送指令。在計算機系統(tǒng)100本地的調制解調器可接收電話線上的數(shù)據(jù)，并使用紅外傳輸器將數(shù)據(jù)轉換為紅外信號。耦合到總線102的紅外檢測器可接收紅外信號中攜載的數(shù)據(jù)并將數(shù)據(jù)放置在總線102上。總線102將數(shù)據(jù)攜載到存儲器106，處理器104從存儲器106檢索并執(zhí)行指令。由存儲器106接收的指令可任選地在由處理器104執(zhí)行之前或之后存儲在存儲裝置110上。

根據(jù)各種實施例，經(jīng)配置以由處理器執(zhí)行以執(zhí)行方法的指令被存儲在計算機可讀媒體上。計算機可讀媒體可為存儲數(shù)字信息的裝置。舉例來說，計算機可讀媒體包含所屬領域已知的用于存儲軟件的光盤只讀存儲器(cd-rom)。計算機可讀媒體由適于執(zhí)行經(jīng)配置以被執(zhí)行的指令的處理器存取。

為了說明及描述的目的，已呈現(xiàn)本教示的各種實施方案的以下描述。所述描述不是窮盡的，并不將本教示限制于所揭示的精確形式。根據(jù)上述教示，修改及變化是可能的，或可從實踐本教示中獲取。另外，所描述的實施方案包含軟件，但本教示可實施為硬件與軟件的組合或僅在硬件中實施。本教示可由面向對象及非面向對象的編程系統(tǒng)來實施。

通過曲線減法發(fā)現(xiàn)的類似區(qū)域

如上文所描述，在復雜樣本中，已知化合物的產(chǎn)物離子的提取離子色譜(xic)可具有兩個或更多個的峰值。因此，可從xic發(fā)現(xiàn)的已知化合物的實際保留時間是模糊的?？赏ㄟ^在每一時間間隔在質量范圍上收集額外產(chǎn)物離子數(shù)據(jù)并使用額外產(chǎn)物離子數(shù)據(jù)來確定正確的保留時間來解決此模糊性?？墒褂妙~外產(chǎn)物離子數(shù)據(jù)來通過已知化合物的其它產(chǎn)物離子的xic峰值來分組產(chǎn)物離子的xic峰值。被發(fā)現(xiàn)與已知化合物的其它產(chǎn)物離子的xic峰值相對應的產(chǎn)物離子的xic峰值較可能是已知化合物的峰值，并且較有可能具有正確的保留時間。遺憾的是，到目前為止，質譜行業(yè)尚無法從在大的質量范圍上收集的多種產(chǎn)物離子的xic獲得選擇已知化合物的峰值群組的系統(tǒng)的、可靠的及準確的方法。

在各種實施例中，通過局部比較xic對中的相鄰強度的群組來可靠且準確地識別來自已知化合物的不同產(chǎn)物離子的xic峰值群組。與'874申請案的方法一樣，此方法依賴于在分離過程的每一時間間隔收集一定質量范圍的額外產(chǎn)物離子數(shù)據(jù)。舉例來說，通過使用一或多個循序質量窗口寬度在每一時間間隔執(zhí)行ms/ms掃描以便橫跨整個質量范圍來收集此額外產(chǎn)物離子數(shù)據(jù)。將來自多次掃描的結果拼接在一起，以在每一時間間隔產(chǎn)生整個質量范圍的光譜。然后可在分離過程的任何時間間隔針對質量范圍中的任何質量計算xic。absciex的swath^tm技術是收集此數(shù)據(jù)的示范性方法。

圖2是根據(jù)各種實施例的已知化合物的五種產(chǎn)物離子的五個xic的示范性繪圖200。從通過使用一或多個循序質量窗口寬度在分離過程的每一時間間隔執(zhí)行ms/ms掃描以便橫跨整個質量范圍而收集的數(shù)據(jù)來計算繪圖200的五個xic。繪圖200的五個xic都表現(xiàn)為在保留時間56周圍具有類似區(qū)域。

在各種實施例中，類似區(qū)域是通過局部比較或減去xic對中的相鄰保留時間處的強度群組來發(fā)現(xiàn)的。舉例來說，并非簡單地在相同的保留時間下將xic220的強度從xic210的強度中減去。而是，在每一保留時間，xic210的強度及在兩個或更多個相鄰保留時間的xic210的強度除以在所述保留時間的xic210的強度，從而有效地歸一化xic210的第一群組強度。在同一保留時間，xic220的強度及在兩個或更多個相鄰保留時間的xic220的強度除以在所述保留時間的xic220的強度，從而有效地歸一化xic220的第二群組強度。接著將第二群組的每一強度從第一群組的對應強度減去，從而產(chǎn)生不同值的集合。通過計算差值集合的統(tǒng)計測量，獲得每一保留時間的單個值。統(tǒng)計測量可為(但不限于)差值集合的平均值、眾數(shù)、中值、方差或標準偏差。

圖3是根據(jù)各種實施例的在圖2中所示的保留時間50與65之間的五個xic的詳細部分的示范性繪圖300。繪圖300展示(例如)如何通過局部比較或減去xic對中的相鄰保留時間的九個強度n＝9的群組來發(fā)現(xiàn)類似區(qū)域。在繪圖300中，在保留時間57處歸一化的xic210的九個強度被展示為a53、a54、a55、a56、a57、a58、a59、a60及a61。這九個值中的每一者例如通過除以a57被歸一化。在保留時間57被歸一化的xic220的對應九個強度為(例如)b53、b54、b55、b56、b57、b58、b59、b60及b61(未展示)。類似地，這九個值中的每一者例如通過除以b57被歸一化。在歸一化之后，從xic210的對應九個強度減去xic220的這九個強度，從而產(chǎn)生不同值的集合。通過計算差值的集合的統(tǒng)計測量，獲得保留時間57的單個值。對比較的xic的每一保留時間計算的統(tǒng)計測量可繪制為比較或減法曲線。

圖4是展示根據(jù)各種實施例的從圖2中所示的xic中的兩個xic的強度的局部減法計算的平均值的減法曲線的示范性繪圖400。例如，根據(jù)等式(1)計算繪圖400中所示的每一平均值μ。

在每一保留時間i+m，第二xic的n個相鄰強度值中的每一強度值b被歸一化，并將其從第一xic的n個相鄰強度值的每一對應歸一化值a減去。n是奇數(shù)，且m是n的中點。

返回圖3，繪圖300展示例如使用九個(n＝9)強度來計算在每一保留時間的平均值。在繪圖300中，展示用于計算在保留時間57的平均值的xic210的九個強度。這九個強度是a53、a54、a55、a56、a57、a58、a59、a60及a61。用于計算在保留時間57的平均值的xic220的對應九個強度為(例如)b53、b54、b55、b56、b57、b58、b59、b60及b61(未展示)。

當從xic210減去xic220時，通過根據(jù)等式(1)從xic210的九個點減去xic220的九個點b53、b54、b55、b56、b57、b58、b59、b60及b61(未展示)來計算保留時間57的平均值μ57。例如，根據(jù)以下等式計算在保留時間57處的平均值μ57

其中九個點的中點m是5。

返回到圖4，可使用展示平均值的減法曲線來識別兩個xic類似的區(qū)域。在兩個xic類似的區(qū)域中，平均值應接近零。繪圖400展示在保留時間56附近的保留時間區(qū)域410，其中平均值接近零。因此，繪圖400表明兩個xic的峰值可分組在保留時間區(qū)域410中。然而，在繪圖400中，平均值經(jīng)?？邕^零值，使得稍微難以區(qū)分類似區(qū)域。

在各種實施例中，通過計算展示標準偏差值的減法曲線，可進一步區(qū)分xic的類似區(qū)域。像平均值μ一樣，從橫跨保留時間的兩個xic中的每一者的區(qū)域中的奇數(shù)n個xic值a及b計算在每一保留時間的標準偏差σ。在數(shù)學上，對于具有數(shù)個點n并且具有m作為n個點的中點的保留時間i+m，由等式(2)給出標準偏差。

返回圖3，例如，當從xic310減去xic320時，通過根據(jù)等式(2)從xic310的九個點減去xic320的九個點b53、b54、b55、b56、b57、b58、b59、b60及b61(未展示)來計算在保留時間57的方差的平方根σ57。舉例來說，根據(jù)以下等式計算方差保留時間57的平方根σ57

其中九個點的中點m是5。

圖5是展示根據(jù)各種實施例的從圖2中所示的xic中的兩個xic的強度的局部減法計算的標準偏差值的減法曲線的示范性繪圖500。繪圖500展示保留時間56附近的保留時間區(qū)域510，其中方差的平方根的值接近零。因此，繪圖500表明兩個xic的峰值可分組在保留時間區(qū)域510中。圖5與圖4的比較展示，相較于從平均值區(qū)分出，兩個xic的類似區(qū)域更容易從方差的平方根的值區(qū)分出。

在各種實施例中，對于已知化合物的m種產(chǎn)物離子的每一集合，執(zhí)行xic的減法，從而產(chǎn)生個減法曲線。然后根據(jù)個減法曲線對m個xic的峰值進行分組。識別具有在零的閾值內的值的個減法曲線中的一或多者的保留時間。換句話說，針對統(tǒng)計比較測量值接近零的位置檢查減法曲線。對于其中統(tǒng)計比較測量值接近零的一或多個保留時間，獲得用于計算識別一或多個保留時間的一或多個減法曲線的兩個或更多個xic。根據(jù)一或多個保留時間內的兩個或更多個xic的峰值創(chuàng)建峰值群組。

圖6是展示根據(jù)各種實施例的從圖2的五個xic計算的標準偏差值的個減法曲線的示范性繪圖600。繪圖600展示在保留時間56附近的保留時間610中，個減法曲線具有接近零的值。這意味著所有五個xic在保留時間區(qū)域610中具有類似的峰值形狀。

圖7是展示根據(jù)各種實施例的在圖6中所示的保留時間50與65之間的標準偏差值的個減法曲線的詳細部分的示范性繪圖700。繪圖700更清楚地展示個減法曲線在保留時間55與57之間的保留時間區(qū)域710中都具有小于0.1的值。因此，保留時間區(qū)域710中的所有五個xic的峰值可被分組。然后可使用類似于'874申請案中使用的準則的準則來對群組的峰值進行評分。然后使用具有最高得分的群組來識別及/或定量已知化合物。

用于分組xic峰值的系統(tǒng)

圖8是展示根據(jù)各種實施例的用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的系統(tǒng)800的示意圖。系統(tǒng)800包含分離裝置810、質譜儀820及處理器830。分離裝置810將已知化合物從樣本混合物分離。分離裝置810可包含(但不限于)電泳裝置、色譜裝置或移動裝置。

舉例來說，質譜儀220是串聯(lián)質譜儀。質譜儀220可包含執(zhí)行兩個或更多個質量分析的一或多個物理質量分析儀。串聯(lián)質譜儀的質量分析儀可包含(但不限于)飛行時間(tof)、四極桿、離子阱、線性離子阱、軌道阱、磁四扇區(qū)質量分析儀、混合四極桿飛行時間(q-tof)質量分析儀或傅里葉(fourier)變換質量分析儀。質譜儀220可分別包含在空間或時間上獨立的質譜分析階段或步驟。

質譜儀220在多個保留時間中的每一保留時間使用一或多個循序質量窗口寬度對分離樣本混合物執(zhí)行一或多次質譜/質譜(ms/ms)掃描，以便橫跨整個質量范圍，從而產(chǎn)生針對多個保留時間的整個質量范圍的產(chǎn)物離子光譜集合。

處理器230與質譜儀220通信。處理器230還可與分離裝置210通信。處理器230可為(但不限于)圖1的系統(tǒng)、計算機、微處理器，或能夠將控制信號及數(shù)據(jù)發(fā)送到串聯(lián)質譜儀220或從所述串聯(lián)質譜儀接收控制信號及數(shù)據(jù)以及處理數(shù)據(jù)的任何裝置。

處理器230從質譜儀220接收多個保留時間的整個質量范圍的產(chǎn)物離子光譜集合。處理器230選擇已知化合物的m種產(chǎn)物離子。處理器220計算來自產(chǎn)物離子光譜集合的m種產(chǎn)物離子中的每一者的xic，從而產(chǎn)生m個xic。處理器220將m個xic中的每一xic從其它m個xic中的每一者減去，從而產(chǎn)生個減法曲線。根據(jù)第一xic及第二xic計算每一減法曲線。在每一保留時間，第一xic在所述每一保留時間的強度及第一xic在兩個或更多個相鄰保留時間的兩個或更多個強度被歸一化。類似地，第二xic在所述保留時間的強度及第二xic在相鄰保留時間的兩個或更多個強度被歸一化。將第二xic的歸一化強度從第一xic的對應歸一化強度減去。計算差值強度的統(tǒng)計測量。統(tǒng)計測量可包含(但不限于)差值強度的平均值、眾數(shù)，中值，方差或標準偏差。

處理器230識別一或多個保留時間的一或多個區(qū)域，其中個減法曲線中的一或多個減法曲線具有在零的閾值內的值。舉例來說，閾值可預先確定或從用戶接收。對于一或多個區(qū)域中的每一區(qū)域，處理器獲得用于計算識別所述區(qū)域的一或多個減法曲線的兩個或更多個xic，并將所述區(qū)域中的兩個或更多個xic的每一峰值添加到峰值群組。

在各種實施例中，在每一保留時間i+m，處理器230根據(jù)以下等式計算差值強度的平均值μi+m

其中n是在每一保留時間計算的差值強度的數(shù)目，n是奇數(shù)，m是n的中點，且對于j＝1到n，是針對保留時間i+m計算的差值強度。

在各種實施例中，在每一保留時間i+m，處理器230根據(jù)以下等式計算差值強度標準偏差σi+m

其中n是在每一保留時間計算的差值強度的數(shù)目，n是奇數(shù)，m是n的中點，且對于j＝1到n，是針對保留時間i+m計算的差值強度，且μi+m是根據(jù)以下等式計算的平均值

在各種實施例中，處理器230進一步對一或多個區(qū)域的每一峰值群組進行評分，并使用最高得分峰值群組來識別或量化已知化合物。

在各種實施例中，處理器230在通過質譜儀220分析樣本混合物之后使用減法曲線分析峰值群組的產(chǎn)物離子光譜集合。換句話說，已知化合物的兩種或更多種產(chǎn)物離子的xic峰值在后處理步驟中進行分組。

分組xic峰值的方法

圖9是展示根據(jù)各種實施例的用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的方法900的示范性流程圖。

在方法900的步驟910中，在多個保留時間獲得整個質量范圍的產(chǎn)物離子光譜集合。使用分離裝置將已知化合物從樣本混合物分離。使用質譜儀在多個保留時間的每一保留時間使用一或多個循序前體離子窗口寬度對分離樣本混合物執(zhí)行一或多次質量ms/ms掃描，以便橫跨整個質量范圍，從而產(chǎn)生針對多個保留時間的整個質量范圍的產(chǎn)物離子光譜集合。

在步驟920中，使用處理器選擇已知化合物的m種產(chǎn)物離子。

在步驟930中，使用處理器計算來自產(chǎn)物離子光譜集合的m個產(chǎn)物離子中的每一者的xic，從而產(chǎn)生m個xic。

在步驟940中，使用處理器將所述m個xic中的每一xic從其它m個xic中的每一者減去，從而產(chǎn)生個減法曲線。根據(jù)第一xic及第二xic計算每一減法曲線。在每一保留時間，第一xic在所述保留時間的強度及第一xic在兩個或更多個相鄰保留時間的兩個或更多個強度被歸一化。類似地，第二xic在所述保留時間的強度及第二xic在相鄰保留時間的兩個或更多個強度被歸一化。將第二xic的歸一化強度從第一xic的對應歸一化強度減去。計算差值強度的統(tǒng)計測量。

在步驟950中，使用處理器識別一或多個保留時間的一或多個區(qū)域，其中個減法曲線的一或多個減法曲線具有在零的閾值內的值。

在步驟960中，針對一或多個區(qū)域中的每一區(qū)域，使用處理器獲得用以計算識別區(qū)域的一或多個減法曲線的兩個或更多個xic，并將區(qū)域中的兩個或更多個xic的每一峰值添加到峰值群組

用于分組xic峰值的計算機程序產(chǎn)品

在各種實施例中，計算機程序產(chǎn)品包含非暫時性及有形計算機可讀存儲媒體，所述存儲媒體的內容包含具有指令的程序，所述指令在處理器上執(zhí)行，以便執(zhí)行用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的方法。此方法由包含一或多個相異軟件模塊的系統(tǒng)執(zhí)行。

圖10是根據(jù)各種實施例的包含執(zhí)行用于分組已知化合物的兩種或更多種產(chǎn)物離子的xic峰值的方法的一或多個相異軟件模塊的系統(tǒng)1000的示意圖。系統(tǒng)1000包含測量模塊1010及分析模塊1020。

測量模塊1010獲得多個保留時間的整個質量范圍的產(chǎn)物離子光集合。使用分離裝置將已知化合物從樣本混合物分離。使用質譜儀在多個保留時間的每一保留時間使用一或多個循序前體離子質量窗口寬度對分離樣本混合物執(zhí)行一或多次質量ms/ms掃描，以便橫跨整個質量范圍，從而產(chǎn)生針對多個保留時間的整個質量范圍的產(chǎn)物離子光譜集合。

分析模塊1020選擇已知化合物的m種產(chǎn)物離子。分析模塊1020從產(chǎn)品離子光譜集合計算m種產(chǎn)物離子中的每一者的xic，從而產(chǎn)生m個xic。

分析模塊1020將m個xic中的每一xic從其它m個xic中的每一者減去，從而產(chǎn)生個減法曲線。根據(jù)第一xic及第二xic計算每一減法曲線。在每一保留時間，第一xic在所述保留時間的強度及第一xic在兩個或更多個相鄰保留時間的兩個或更多個強度被歸一化。類似地，第二xic在所述保留時間的強度及第二xic在相鄰保留時間的兩個或更多個強度被歸一化。將第二xic的歸一化強度從第一xic的對應歸一化強度減去。計算差值強度的統(tǒng)計測量。

分析模塊1020識別一或多個保留時間的一或多個區(qū)域，其中個減法曲線的一或多個減法曲線具有在零的閾值內的值。針對一或多個區(qū)域中的每一區(qū)域，分析模塊1020獲得用以計算識別區(qū)域的一或多個減法曲線的兩個或更多個xic，并將區(qū)域中的兩個或更多個xic的每一峰值添加到峰值群組。

雖然結合各種實施例描述了本教示，但并不希望將本教示限于此類實施例。相反，如所屬領域的技術人員應了解，本教示涵蓋各種替代例，修改及等效物。

此外，在描述各種實施例時，說明書可呈現(xiàn)作為特定步驟序列的方法及/或過程。然而，在所述方法或過程不依賴于本文所闡述的步驟的特定次序的程度上，所述方法或過程不應限于所描述的特定步驟序列。如所屬領域的一般技術人員將了解，其它步驟序列可為可能的。因此，說明書中闡述的步驟的特定次序不應被解釋為對權利要求書的限制。另外，針對方法及/或過程的權利要求不應限于按照書面次序執(zhí)行其步驟，并且所屬領域的技術人員可容易地了解，序列可變化并且仍然保持在各種實施例的精神及范圍內。