(57)【要約】多価種の低分解能質量電荷比スペクトルから質量情報を抽出する装置および方法が開示されている。曲線適合は、雑音によって誘発される外乱を除外する質量電荷データの対数を使用して実行される。
2019521329.tif 000002
【特許請求の範囲】
【請求項１】
多価種の内の多原子親分子の質量を同定するために低分解能質量電荷比スペクトルから質量情報を抽出する方法であって、
多価イオンの集団を示すデータセットを質量分析計から受信するステップであり、各イオン上の電荷の数はそのイオンの荷電状態を規定し、各荷電状態は前記イオンの集団内のイオンのサブ集団からなる、ステップと、
前記イオンの荷電状態を規定するデータセットを使用して入力質量電荷比スペクトルを生成するステップであり、前記各荷電状態のサブ集団は前記質量電荷比スペクトルにおける強度によって表される、ステップと、
前記入力質量電荷比スペクトルを処理して信号増強化質量電荷比スペクトル（１０３）を提供するステップであり、前記信号増強化質量電荷比スペクトルは前記質量電荷比スペクトルおよび前記質量電荷比スペクトルの平滑化された表現の商の対数から生成される、ステップと、
規定された荷電質量スペクトル（１０４）を使用して、前記多原子親分子の質量をその中で探索するための規定された荷電質量値の範囲を同定するステップと、
前記規定された荷電質量値の範囲内の各質量に対し、最大荷電状態までの連続した荷電状態で前記質量に対応する前記信号増強化質量電荷比スペクトルにおける値の加算に等しい総和を関数
【数１】
2019521329.tif 000016
を使用して生成するステップであり、ここで、Ｍは規定された荷電質量値の前記範囲内の任意の規定された質量であり、ｍaは電荷付加物の質量であり、ｚは規定された荷電質量スペクトルの次数であり、Ｘは信号増強化質量電荷比スペクトルの分布関数である、ステップと、
前記総和の値を用いて前記多原子親分子の質量を決定するステップと
を含む方法。
【請求項２】
前記総和値を使用して多原子親分子の質量を決定するステップは、前記親分子の前記質量を決定するために、前記総和の範囲にわたって前記総和からの値を正規化するステップを含む、請求項１に記載の方法。
【請求項３】
前記電荷輸送付加物の質量、ｍa、は、プロトン質量を表す１に等しくなるように設定される、請求項１または２に記載の方法。
【請求項４】
前記信号増強化質量電荷比スペクトル（１０３）は関数
【数２】
2019521329.tif 000017
から生成され、
ここで、Ｘは前記信号増強化質量電荷比スペクトルを表し、Ｉは入力質量電荷比スペクトル（２０１）内の強度に対応し、Ｓは前記入力質量電荷比スペクトルの平滑化された表現である、請求項１乃至３のいずれか一項に記載の方法。
【請求項５】
前記規定された荷電質量値の前記規定された範囲内の各質量について、前記方法は、
最大荷電状態の異なる値までの合計を計算するステップと、
前記最大荷電状態の異なる値までの前記合計からの前記値を使用して前記親分子の分子量を決定するステップと
をさらに含む、請求項１乃至４のいずれか一項に記載の方法。
【請求項６】
前記信号増強化質量電荷比スペクトルを形成する前に、前記入力質量電荷比スペクトルのスペクトル前処理をさらに含む、請求項１乃至５のいずれか一項に記載の方法。
【請求項７】
前記スペクトル前処理が、平滑化および／またはベースライン減算のうちの少なくとも１つから選択される、請求項６に記載の方法。
【請求項８】
前記多価種がエレクトロスプレイイオン化によって生成される、請求項１乃至７のいずれか一項に記載の方法。
【請求項９】
前記入力質量電荷比スペクトルを生成するために、質量分析計源内で多原子親分子をイオン化するステップを含む、請求項１乃至８のいずれか一項に記載の方法。
【請求項１０】
前記規定された電荷対質量スペクトルは、全ての同定された種がゼロ荷電値を有するゼロ荷電質量スペクトルである、請求項１乃至９のいずれか一項に記載の方法。
【請求項１１】
前記規定された荷電質量スペクトルが、一価質量値の範囲を同定するために使用される一価質量スペクトルである、請求項１から９のいずれか一項に記載の方法。
【請求項１２】
エレクトロスプレイイオン化ソースおよび検出器を備える質量分析器であって、前記検出器は前記イオン化ソースによって生成された多価イオンの集団を示す出力を生成するように構成され、各イオン上の電荷数はそのイオンの荷電状態を規定し、各荷電状態は前記イオンの前記集団内のイオンのサブ集団からなり、前記分析器は請求項１乃至１１のいずれか一項に記載の方法を実施するように構成されたプロセッサをさらに備える、質量分析器。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、質量分析に関し、特に、多価種の低分解能質量電荷比スペクトルから質量情報を抽出する装置および方法に関する。
【背景技術】
【０００２】
質量分析におけるエレクトロスプレイイオン化の使用は一般的になっており、特に、大きな生体分子の研究におけるその使用が普及している。反応活性種を完全な形で研究することができるので有利ではあるが、多重荷電は、複雑な質量電荷比スペクトルにつながる可能性がある。これらのスペクトルは、異なる質量電荷比でのイオンカウントの観測であるため、この多重荷電により、親分子が複数のピークを示す結果となり、親質量を決定することが困難になる可能性がある。その結果、これらの複素スペクトルから親分子の質量に関する情報を得るために、多くの方法が開発されている。この情報を抽出するためには、観測された親分子イオンの荷電状態に関する知識が必要である。高分解能スペクトル（同位体分解能を有するスペクトル）を用いることにより、これらの荷電状態の割り当てが、隣接する同位体ピーク間の分離から容易に推測されることができる。低分解能の質量電荷比スペクトル（同位体分解能を有するスペクトル）では、親分子質量の抽出はそれほど簡単ではない。
【０００３】
生体分子のような大きな分子の研究のための低分解能質量分析法の使用は、いくつかの利点を提示する。これらの利点は、より低い分解能の質量分析計を低コストで製造する性能に由来し、そのような装置の小型化もより容易である。これらの利点により、オンラインおよびアットライン分析のための低分解能の機器の使用を非常に魅力的になる。サンプリングがサンプル調製がほとんどまたはまったくない環境の範囲で行われることが多いため、この質量分析の使用法にはさらなる課題がある。この結果、質量電荷比スペクトルは雑音が多いように見えることがあり、サンプル媒体に由来する溶媒クラスタおよびイオン付加物によって引き起こされる異質なピークを含む可能性がある。したがって、これらのシナリオでの親分子の抽出は非常に困難である。特に、強い異質なピークおよび劣悪な信号対雑音が質量電荷比スペクトルの特性である場合、ピークピッキングに依存する技術は、容易に失敗する可能性がある。さらに、その性質により、ピークピッキングは、ピークを単一の質量電荷比に減少させる結果となる。したがって、ピークピッキング技法は、肩を有するピークまたは任意の形態の非対称性の誤った表現につながる可能性がある。他のデコンボリューション法は、バックグラウンド雑音の影響を強く受ける可能性があり、親分子の真の質量を表さないアーチファクトピークの生成をもたらす。
【０００４】
オンラインおよびアットラインモニタリングのようなアプリケーションに安価で小型化された低分解能計測器を使用することは、追加の要件を提示する。これらの環境での使用は、質量分析において、必然的に、熟練した実践者による使用から技術および経験がより浅いユーザへの移行をもたらす。結果として、多重荷電質量電荷比スペクトルから親分子質量を抽出するために使用されるいかなる技法も使用が容易であり、最小のユーザ入力を必要となることが望ましい。
【０００５】
同位体分解能のない質量電荷比スペクトルから質量情報を抽出するための多くの方法が実施されている。これらについては以下で説明する。
【０００６】
ＺｈａｎｇとＭａｒｓｈａｌｌは、米国質量分析学会誌，第９巻，２２５−２３３（１９９８）において「ＺＳｃｏｒｅ」という方法を記載している。この方法は、信号対閾値比の対数に基づいてピークピッキングおよびスコアリングシステムを利用する。比率は、バックグラウンド雑音とユーザ定義の信号対雑音比から計算される。この方法はピークピッキングに依存するため、以前に概説したピークピッキングの欠点がある。バックグラウンド雑音を計算するために、ユーザが信号対雑音比および質量電荷比の範囲を定義する必要があることは、このシステムの別の欠点である。特に、雑音レベルが不適切なレベルに設定されていると、有用な情報が失われ得る。
【０００７】
ＭｏｒｇｎｅｒおよびＲｏｂｉｎｓｏｎは、分析化学、第８４巻，２９３９−２９４８（２０１２）において「Ｍａｓｓｉｇｎ」という方法を記載している。上述のＺＳｃｏｒｅ法と同様に、この方法は、正確なピークピッキングおよび２つの閾値レベルの形でのユーザ入力に依存する。さらに、分離が不十分なピークは、計算に含めるために手動で同定する必要がある。したがって、この方法は、低質量の質量電荷比スペクトルから親質量を抽出するのには適していない。
【０００８】
質量分析の国際学術誌、第２９０巻、１−８（２０１０）および分析化学、第７７巻、１１１１−１９（２００５）において、Ｍａｌｅｋｎｉａ他は、荷電状態からの整数の固有の比を利用して、低分解能の質量電荷比スペクトルにおけるピークの荷電状態を計算する「ｅＣＲＡＭ」法を記載している。この方法はまた、分析されるスペクトルのピークをはっきりと同定することにも依存する。
【０００９】
Ｗｉｎｋｌｅｒは、Ｒａｐｉｄ  Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ  ｉｎ  Ｍａｓｓ  Ｓｐｅｃｔｏｍｅｔｒｙ、第２４巻、２８５−２９４（２０１０）において、質量電荷比スペクトルにおいて観察されるピークを使用して、異なる荷電状態で物質種の質量を計算する「ＥＳＩｐｒｏｔ」法を記載している。正しい荷電状態は、ベッセルの補正を用いて最低標準偏差を生じる荷電状態の集合を計算することによって同定される。この方法の性能はまた、質量スペクトルにおける関心のあるピークを明確に同定する能力にも依存する。この方法で同定され使用されるピークが連続した一連の多価イオンからのものであるという仮定はまた、この方法の重大な制限である。
【００１０】
Ｍａｎｎ他は、分析化学、第６１巻、１７０２−１７０８（１９８９）において、米国特許第５，１３０，５３８号にも記載されている２つの方法を記載している。第１の方法は、相対ピーク位置を用いて荷電状態を推定し、親分子の質量を計算する平均化アルゴリズムである。この方法はピークピッキングに依存しており、雑音および異質のピークの影響を強く受ける可能性がある。これらの２つの出典に記載された第２の方法は、デコンボリューションアルゴリズムである。この方法は、変換関数を使用して親分子質量の試行値を評価する。この変換関数は、質量電荷比スペクトルにおけるイオンカウントの分布関数から計算される。このデコンボリューションの方法は、ピークピッキングに依存しないため恩恵を受けるが、どちらの引例でも示すように、バックグラウンド雑音の影響を強く受ける。これは、デコンボリューションされたスペクトル内に出現する、親分子の倍数または分数の形での誤った結果の生成がデコンボリューションされたスペクトルに現れる可能性がある。デコンボリューションされたスペクトルにおける質量に伴うバックグラウンド雑音の増加は、この方法の付加的な望ましくないアーチファクトである。
【００１１】
デコンボリューションアルゴリズムを利用する方法は、分析化学、第６６巻、１８７７−１８８３（１９９４）および米国特許第５，３５２，８９１号にも記載されている。上述のＭａｎｎアルゴリズムとは異なり、イオンカウントは、一緒に加算されるのではなく、乗算相関アルゴリズムを形成するために一緒に乗算される。この方法論の変更により、バックグラウンド雑音の影響とアーチファクトのピークの導入が減少する。しかしながら、この乗法技術は、ゼロに等しい単一のイオンカウントまたは異常に高いイオンカウントの存在によって圧倒的に不利な条件に置かれる可能性がある。このゼロイオンカウントを計算に含めることは、デコンボリューションされたスペクトルにおける信号の完全な抑制につながる。同様に、単一の大きなイオンカウント（例えば、付加イオンまたは雑音から生じる）の包含は、デコンボリューションされたスペクトルにおける信号に過度に影響を及ぼす可能性がある。
【００１２】
最初にＲｅｉｎｈｏｌｄ他により米国質量分析学会誌、第３巻、２０７−２１５（１９９２）に記載された、フォワードワーキング最大エントロピー法は、アーチファクトピークに対する判別を加えたので有利である。しかしながら、この方法は計算コスト高になる可能性があり、親分子質量の効果的な抽出には、質量対強度比スペクトルにおける相対強度およびピーク形状の先験的な知識が必要である。
【００１３】
必要とされるのは、低品質の低分解能の質量電荷比スペクトル、すなわち同位体分解能、低信号対雑音および異質な付加ピークを含むスペクトルから多価親分子の質量を抽出することができる堅牢な技法である。さらに、この技法は、ピークピッキング技法の使用に依存すべきではなく、バックグラウンド雑音の影響を強く受け、アーチファクトピークの生成を受けやすいものであるべきである。この技法はまた、最小限のユーザ入力と事前知識を必要とすべきであり、そのため、質量分析の専門家でないユーザがこの方法を正常に実施することが可能になる。
【先行技術文献】
【特許文献】
【００１４】
【特許文献１】米国特許第５，３５２，８９１号公報
【特許文献２】米国特許第５，１３０，５３８号公報
【非特許文献】
【００１５】
【非特許文献１】ZhangとMarshall、Journal of the American Society for Mass Spectrometry、１９９８年、第９巻、p.２２５−２３３
【非特許文献２】MorgnerとRobinson、Analytical Chemistry、２０１２年、第８４巻、p.２９３９−２９４８
【非特許文献３】Maleknia  他、the International Journal of Mass Spectrometry、２０１０年、第２９０巻、p.１−８
【非特許文献４】Maleknia  他、Analytical Chemistry、２００５年、第７７巻、p.１１１−１１９
【非特許文献５】Winkler、Rapid Communications in Mass Spectrometry、２０１０年、第２４巻、p.２８５−２９４
【非特許文献６】Mann  他、Analytical Chemistry、１９８９年、第６１巻、p.１７０２−１７０８
【非特許文献７】Reinhoid  他、the Journal of the American Society for Mass Spectrometry、１９９２年、第３巻、p.２０７−２１５
【発明の概要】
【００１６】
これらおよびその他の課題は、本教示により、多価種の低分解能質量電荷比スペクトルから質量情報を抽出できる装置および方法適応型アパーチャチューナブルアンテナのためのシステムおよび方法によって対処される。
【００１７】
本教示の第１の態様によると、親分子の多価イオンの集団を生成するために多原子親分子の初期イオン化を含む方法が提供される。各イオンについて、存在する電荷の数は、そのイオンの荷電状態を規定し、各荷電状態は、イオンの集団内のサブ集団からなる。これらのサブ集団の分析は、質量電荷比スペクトルをもたらし、各質量電荷比の強度は、各荷電状態の母集団の直接的な表現である。
【００１８】
質量電荷比スペクトルが生成された後、本教示は、雑音を低減し、得られた質量電荷比スペクトルのベースラインをゼロに減少させるための前処理ステップを使用し得る。前処理の後、質量電荷比スペクトルは、以下の関数を使用してスペクトルの新しい表現に変換される。
【００１９】
【数１】
2019521329.tif 000003
【００２０】
ここで、Ｘは質量電荷比スペクトルの新しい表現であり、Ｉは前処理された入力質量電荷比スペクトルの強度に対応し、ＳはＩの平均中心表示である。質量電荷比スペクトルのこの新しい表現は、ゼロ荷電で全ての物質種を有する質量スペクトルを生成するように変換することができる。同様に、すべての一価の物質種を有する質量電荷比スペクトルが同じ趣旨で生成され得ることが理解されるであろう。
【００２１】
本教示をゼロ荷電の例に限定することは意図していないが、本教示の理解を助けるために、ゼロ荷電状態の表現を詳細な説明で使用し、以下の関数を使用して計算する。
【００２２】
【数２】
2019521329.tif 000004
【００２３】
所定の範囲内のＭの異なる試行値に対して、信号増強化質量電荷比スペクトル（Ｘ）の値は、最大荷電状態ｎ’maxに合計される。最大荷電状態ｎ’maxの異なる値に対して反復することにより、Ｆ（Ｍ、ｎ’max）の値は、オーバーフィッティングおよびゼロ荷電質量スペクトルにおけるアーチファクトの出現を防止するように最適化される。
【００２４】
一価質量スペクトルを評価するための関数は、上記のゼロ荷電関数から導出され、一価質量値の範囲を特定するために使用し得る。一価状態の表現は、以下の関数を使用して計算し得る。
【００２５】
【数３】
2019521329.tif 000005
【００２６】
上記の式は、以下の関数の使用によってさらに一般化し得る。
【００２７】
【数４】
2019521329.tif 000006
【００２８】
上記のゼロ荷電関数は、一般関数においてｚを０と等しくすることによって得ることができ、一方、一価状態表現についての上記関数は、一般関数においてｚを１と等しくすることによって得ることができる。
【００２９】
ゼロ荷電状態表現が最も一般的に使用されており、明確化のために、詳細な説明を通して使用される。一価表示では、付加物イオンの質量も考慮される。この表現は、ＥＳＩ−ＭＳ実験からの一価種の未処理の質量電荷スペクトルとの直接比較が必要な場合に有用である。
【００３０】
このような方法は、完全な質量電荷比スペクトルを有利に分析し、このようにして全ての成分が等しく処理されるので、ピークピッキング法および雑音レベルの先験的な知識に依存しない。変換された質量電荷比スペクトルの使用は、バックグラウンド雑音が出力に過度に影響することを防止する、従来技術の他の方法で見られる現象である。その結果、この方法は、最小限のユーザ入力と先験的な知識で、低品質で雑音の多い質量電荷比スペクトルからの親分子質量の効果的な抽出を可能にする。したがって、本教示の第１の態様は、多価種の低原子質量電荷比スペクトルから質量情報を抽出して、多価種内の多原子親分子の質量を同定する方法を提供し、本方法は、
多価イオンの集団を示すデータセットを質量分析計から受信するステップであり、各イオン上の電荷の数はそのイオンの荷電状態を規定し、各荷電状態は前記イオンの集団内のイオンのサブ集団からなる、ステップと、
イオンの荷電状態を規定するデータセットを使用して質量電荷比スペクトルを生成するステップであり、各荷電状態のサブ集団は質量電荷比スペクトルにおける強度によって表される、ステップと、
入力質量電荷比スペクトルを処理して信号増強型質量電荷比スペクトル（１０３）を提供するステップであり、信号増強化質量電荷比スペクトルは質量電荷比スペクトルおよび質量電荷比スペクトルの平滑化された表現の商の対数から生成される、ステップと、
規定された荷電質量スペクトル（１０４）を使用して、多原子親分子の質量をその中で探索するための規定された荷電質量値の範囲を同定するステップと、
規定された荷電質量値の範囲内の各質量に対し、最大荷電状態までの連続した荷電状態で質量に対応する信号増強化質量電荷比スペクトルにおける値の加算に等しい総和を関数
【００３１】
【数５】
2019521329.tif 000007
【００３２】
を使用して生成するステップであり、ここで、Ｍは規定された荷電質量値の範囲内の任意の規定された質量であり、ｍaは電荷付加物の質量であり、ｚは規定された荷電質量スペクトルの次数であり、Ｘは信号増強化質量電荷比スペクトルの分布関数である、ステップと、
前記総和値を用いて多原子親分子の質量を決定するステップとを含む。
【００３３】
本教示のさらなる発展において、総和値を使用して多原子親分子の質量を決定するステップは、親分子の質量を決定するために、総和の範囲にわたって総和からの値を正規化するステップを含む。
【００３４】
本教示のさらなる発展において、電荷輸送付加物の質量、ｍa、は、プロトン質量を表す１に等しくなるように設定される。
【００３５】
本教示のさらなる発展において、信号増強化質量電荷比スペクトル（１０３）は関数
【００３６】
【数６】
2019521329.tif 000008
【００３７】
から生成され、
ここで、Ｘは信号増強化質量電荷比スペクトルを表し、Ｉは入力質量電荷比スペクトル（２０１）内の強度に対応し、Ｓは入力質量電荷比スペクトルの平滑化された表現である。
【００３８】
本教示のさらなる発展において、規定された荷電質量値の規定された範囲内の各質量について、方法は、最大荷電状態の異なる値までの合計を計算するステップと、
最前記大荷電状態の異なる値までの前記合計からの前記値を使用して前記親分子の分子量を決定するステップとを含む。
【００３９】
本教示のさらなる発展において、方法は、信号増強化質量電荷比スペクトルを形成する前に、入力質量電荷比スペクトルのスペクトル前処理をさらに含む。
【００４０】
本教示のさらなる発展において、スペクトル前処理が、平滑化および／またはベースライン減算のうちの少なくとも１つから選択される。
【００４１】
本教示のさらなる発展において、多価種がエレクトロスプレイイオン化によって生成される。
【００４２】
本教示のさらなる発展において、方法は、力質量電荷比スペクトルを生成するために、質量分析計源内で多原子親分子をイオン化するステップをさらに含む。
【００４３】
本教示のさらなる発展において、規定された電荷対質量スペクトルは、全ての同定された物質種がゼロ荷電値を有するゼロ荷電質量スペクトルである。
【００４４】
本教示のさらなる発展において、規定された荷電質量スペクトルが、一価質量値の範囲を同定するために使用される一価質量スペクトルである。
【００４５】
本発明の第２の態様により、エレクトロスプレイイオン化ソースおよび検出器を備える質量分析器であって、検出器は前記イオン化ソースによって生成された多価イオンの集団を示す出力を生成するように構成され、各イオン上の電荷数はそのイオンの荷電状態を規定し、各荷電状態はイオンの集団内のイオンのサブ集団からなり、分析器は先行する上記の方法を実施するように構成されたプロセッサをさらに備える、質量分析器が提供される。
【００４６】
本教示のこれらおよび他の態様は、本教示の理解を助けるために提供される以下の図を参照して説明されるが、決して限定的なものとして解釈されるべきではない。
【図面の簡単な説明】
【００４７】
【図１】図１は、本教示の第１の態様で使用される処理ステップのフローチャートである。
【図２】図２は、ウシ血清アルブミンの質量電荷比スペクトルにおける図１のステップ１０２の実施の例である。
【図３】図３は、信号増強化質量電荷比スペクトルを生成するための図１のステップ１０３の実施の例である。
【図４】図４は、図２に示す質量電荷比スペクトルについて図１を参照して一般的に説明した技術を用いて生成されたゼロ荷電スペクトルの例である。
【図５】図５は、本教示の第２の態様による、図１の方法の変更例で使用されるステップのフローチャートである。
【図６】図６は、図２のデータに対して図５の方法を用いて生成されたゼロ荷電スペクトルの例である。
【図７】図７は、図１の方法の別の変更例で使用されるステップのフローチャートである。
【図８】図８は、信号対雑音が非常に低いウシ血清アルブミンの質量電荷比スペクトルである。
【図９】図９は、図８のデータに対するステップ７０２の実施の例である。
【図１０】図１０は、図９の平滑化されたデータに対するステップ７０３の実施を示す。
【図１１】図１１は、図１０のデータに対するステップ７０４の実施の例である。
【図１２】図１２は、図１１のデータに対するステップ７０５の実施を示す。
【図１３】図１３は、図８に示すデータに対して図７の方法を使用して生成されたゼロ荷電スペクトルの例である。
【図１４】図１４は、図１、図５または図７のいずれか１つの方法を実施するために、本教示の文脈の中で使用され得るコンピュータ処理装置の例である。
【図１５】図１５は、３つの電荷ピークを同定する理想化された質量スペクトルの例である。
【図１６Ａ】図１６Ａおよび図１６Ｂは、図１５のデータに対する荷電質量値の範囲の計算における加算技術の効果を示す。
【図１６Ｂ】図１６Ａおよび図１６Ｂは、図１５のデータに対する荷電質量値の範囲の計算における加算技術の効果を示す。
【図１７】図１７は、図１５のデータの平均中心データセットの対数を使用して信号増強化質量電荷比スペクトルを作成する適用の効果の例である。
【図１８Ａ】図１８Ａおよび図１８Ｂは、荷電質量値の範囲の計算において図１７の信号増強化質量電荷比スペクトルへの変換を用いる効果を示す。
【図１８Ｂ】図１８Ａおよび図１８Ｂは、荷電質量値の範囲の計算において図１７の信号増強化質量電荷比スペクトルへの変換を用いる効果を示す。
【図１９】図１９は、単一の不純物を添加した図１５と同様のスペクトルを示す。
【図２０Ａ】図２０Ａは、単純な雑音の減算を用いて図１９のデータを処理する効果を示す。
【図２０Ｂ】図２０Ｂは、図１９のデータを商およびログ関数で処理する効果を示す。
【発明を実施するための形態】
【００４８】
図１は、本教示による低分解能の質量電荷比スペクトルから質量情報を抽出する方法を表す流れ図を示す。
【００４９】
ステップ１０１において、質量スペクトルは、従来の質量分析技術を用いて測定され、質量電荷比スペクトルの形態で保存される。当業者には理解されるように、このような技術は、親分子の異なる荷電状態のサブ集団を異なる質量電荷比でイオンカウントとして表すことを可能にする。
【００５０】
ステップ１０２において、入力質量電荷比スペクトルからの傾向を表すためにデータに対して平滑化技術が用いられる。平滑化技術のタイプは様々であり、例えば、平滑化は、中心化移動平均を使用して効果的にデータを中心化し得る。代替の構成では、データを平滑化または中心化を行うために、平均値を含むがこれに限定されない他の移動平均を利用することができる。指数平滑法など他の平滑化方法も実施することができる。図２は、ウシ血清アルブミンの入力質量電荷比スペクトルに対するステップ１０２の実施を示す。この例示的なデータセットでは、ウシ血清アルブミンの質量電荷比スペクトルは、「雑音のある」シグナル２０１によって表され、２０２は、平均中心化データを用いたこの態様で、信号２０１のデータに対する平滑化アルゴリズムの効果を示す。
【００５１】
ステップ１０３において、ピークの存在を強化し、雑音の影響を低減する質量電荷比スペクトル（以下、「信号増強化質量電荷比スペクトル」と呼ぶ）が生成される。このスペクトルの生成は、（ステップ１０１からの）入力スペクトルからのデータの商−データ２０１によって表されるデータセット−およびステップ１０２において生成された平均中心化スペクトルの対数を取ることによって達成される。
【００５２】
【数７】
2019521329.tif 000009
【００５３】
この式において、Ｘは信号増強化質量電荷比スペクトルを表し、Ｉは入力質量電荷比スペクトル（２０１）における強度に対応し、ＳはＩ（２０２）の平均中心化表示である。典型的には、このステップは自然対数の使用を含む。代替的な構成では、任意の塩基の対数を用いて同様の効果を得ることができる。このステップの処理を図２に示すデータに適用することのデモンストレーションを図３に示す。図３から、ＩおよびＳの商を使用することにより、各ピークからの寄与がより均等になり、単一の付加ピークまたは雑音アーチファクトに起因する可能性のある少数の非常に強いピークによって不当に影響されないことが保証されることが明らかである。この対数変換によって達成される技術的効果は、未解決ピーク間の雑音および極小値の合計からの正味の寄与が無効にされることである。この処理ステップは、事前スペクトルピークピッキングに依存せずに、入力スペクトル内のすべてのデータを使用することを可能にする。Ｘの出力は、データ配列内のデータ点のセットによって等しく表されることができる。
【００５４】
ステップ１０４では、ステップ１０８で生成されるゼロ荷電スペクトルの最大値と最小値が定義される。これらの値は、ステップ１０１からの入力スペクトルの最大値および最小値から推測することができる。このステップでは、ステップ１０８で評価される最大充電状態の値もまた、ｎmaxとして定義される。本教示のこの態様によれば、ｎmaxは、ステップ１０４で定義されたゼロ荷電スペクトルの最大値および最小値、およびステップ１０１で使用される質量電荷比スペクトルの範囲から計算される。代替構成では、出力の特定のユーザー最適化が必要な場合、ｎmaxはユーザーにより定義できる。式３を参照して以下の議論から明らかなように、そのような最適化は自動的に達成できる。他の構成では、ｎmaxは、処理ルーチン内で一定値として予め定義されてもよい。ステップ１０４は、図１の構成においてステップ１０１から１０３に続くものとして示されている。本教示の文脈内で使用される代替の処理ルーチンは、ステップ１０４がステップ１０５以前の任意の時点で行うことを可能にする。
【００５５】
ステップ１０５において、関数Ｆ（Ｍ）は、ステップ１０４で定義されたゼロ荷電質量値の範囲について、下記を使って評価される。
【００５６】
【数８】
2019521329.tif 000010
【００５７】
式２において、Ｍは、ステップ１０４で定義された範囲内の任意のゼロ荷電質量である。変数ｍaは、電荷運搬付加の質量であり、Ｘは、ステップ１０３で評価される信号増強化質量電荷比スペクトルの分布関数である。
【００５８】
一連の値をＩの適切な値まで近似すると、計算時間が短縮される。典型的には、付加イオン（ｍa）の質量は、プロトンを表す１に等しく設定されるが、代替的な配置は、最も一般的な付加イオンに対して適切な任意の値にｍaを設定することができ、ｍaの値の実際の選択必要な方法論と精度の詳細に依存する。このように、異なる方法論では、異なる値に設定する必要があり、これは調査されるサンプルの先験的知識を必要とする。
【００５９】
対数関数を含む信号増強化質量電荷比スペクトルを用いることにより、Ｍが親分子の質量に対応する場合にはＦ（Ｍ）の総和が増強され、Ｍが対応しない場合には減少する。Ｍが親分子の質量に対応する場合、（Ｍ／ｎ＋ｍa）の値は元の質量電荷比スペクトルのピークと一致し、Ｘは正の値を生じる。Ｍが親分子の質量に対応していない場合、（Ｍ／ｎ＋ｍa）の値は元の質量電荷比スペクトルの雑音または最小値と一致し、合計に対する正味の寄与は負となる。Ｍは実際にＦ（Ｍ）の値に負の影響を与えるため、Ｆ（Ｍ）の評価のための信号増強化質量電荷比スペクトルＸの使用はまた、親分子質量の複数の値におけるアーチファクピークの減少が存在することを意味する。これは、バックグラウンドの追加サンプリングにより、親分子の質量の倍数が強調され、Ｍの増加に伴ってベースラインが増加する強度（Ｉ）の単純な合計とは対照的である。
【００６０】
Ｆ（Ｍ）が一旦評価されると、ステップ１０６でゼロ荷電スペクトルを生成することができる。次いで、ステップ１０７でＦ（Ｍ）の値を正規化して、専門知識のレベルよってエンドユーザにとってより管理しやすい出力を生成することができる。図１の方法を使用して処理されたときの図２のデータのゼロ荷電スペクトルが図４に示されている。このデータ分析の提示は、採用可能な様々な異なる技術のうちの１つに過ぎないことが理解されるであろう。例えば、代替的な構成では、出力は、質量値の表の形式であってもよいし、単なる質量の値であってもよい。
【００６１】
この方法で使用される対数の総和のために、この方法によって生成されたピークの高さは、存在する物質種の定量的尺度として使用することができない。しかしＦ（Ｍ）の最大値に対してデータを正規化して解釈が容易な出力を得ることができる。この方法によって調査される親分子の質量が得られると、元の質量電荷比スペクトルから他の定量データを容易に抽出することができる。例えば、親分子の質量を用いて適切な質量電荷比を計算することができ、元の質量電荷比スペクトルにおけるこれらの値の振幅を合計することができる。したがって、図１の技法によれば、調査されている分子の性質に関する事前情報をユーザに有利に提供し、この情報をより詳細な処理ステップで使用することができることが理解されよう。
【００６２】
本教示による低分解能の質量電荷比スペクトルから質量情報を抽出するための上記の方法は、ステップ１０８において一価スペクトルを出力するように実行され得ることが理解されよう。
【００６３】
例えば、ステップ１０４において、単電荷スペクトルの最大値と最小値が定義される。次いで、関数Ｆ（Ｍ）は、ステップ１０５で、ステップ１０４で定義された一価質量値の範囲について、下記式を使って評価される。
【００６４】
【数９】
2019521329.tif 000011
【００６５】
式２ａにおいて、Ｍは、ステップ１０４で定義された範囲内の任意の一価質量である。変数ｍａは、電荷輸送付加物の質量であり、Ｘは、ステップ１０３で評価された信号増強化質量電荷比スペクトルの分布関数である。
【００６６】
式２ａによるＦ（Ｍ）の評価に続いて、図１のステップ１０６において、一価スペクトルを生成し得る。ステップ１０７でＦ（Ｍ）の値を正規化して、ステップ１０８で一価質量スペクトル出力を生成し得る。
【００６７】
ゼロおよび一価の表現は、さらに一般化し得る。この一般化された例では、ゼロまたは一価スペクトルの最大値および最小値がステップ１０４で定義される。関数Ｆ（Ｍ）は、ステップ１０４で定義された零および一価質量値の両方について、下記式を使ってステップ１０５で評価され得る。
【００６８】
【数１０】
2019521329.tif 000012
【００６９】
ゼロ荷電表現（式２）は、ｚを０に等しくすることによって得られ、一価荷電表現（式２ａ）は、ｚを１に等しくすることによって得ることができる。
【００７０】
式（２ｂ）において、Ｍは、ステップ１０４で定義された範囲内の任意のゼロまたは単電荷の質量である。変数ｍaは電荷輸送付加物の質量であり、ｚはスペクトルの次数である（ゼロ荷電表現に対してはｚ＝０であり、一価表現に対してはｚ＝１）、Ｘはステップ１０３で評価された信号増強質量電荷比スペクトルの分布関数である。
【００７１】
式２ｂに従ってＦ（Ｍ）を評価した後、図１のステップ１０６においてゼロまたは単一電荷スペクトルを生成し得る。Ｆ（Ｍ）の値は、ステップ１０８においてｚ−荷電質量スペクトル出力を生成するようにステップ１０７において正規化し得る。
【００７２】
図１の構成の変形例では、図１のステップ１０１から１０４は保持されるが、そのデータセットの追加の処理が使用される。この変形例は、図５のフローチャートを参照して説明され、ステップ５０１から５０４は、先に説明したステップ１０１から１０４と同一である。ステップ５０５から５０９は、親分子の決定された質量の精度を改善し、ｎmaxを割り当てることの困難を取り除くために導入され得る追加の最適化技術を表す。
【００７３】
ステップ５０５において、一時的なパラメータｎ’maxが定義され、最初は１に等しく設定される。
【００７４】
図１を参照して概説したように、信号増強化質量電荷比スペクトルの使用は、バックグラウンドまたは最小値をサンプリングすることによる正味の結果が総和に対する負の寄与であるので有益である。これは、ｎmaxが増加し、バックグラウンドがオーバーサンプリングされると、Ｆ（Ｍ）の値が減少することを意味する。これを入力質量電荷比スペクトル（Ｉ）からの強度の簡単な合計と比較すると、そのような技術は、Ｆ（Ｍ）の値が増加させるｎmaxの高い値でバックグラウンドのオーバーサンプリングを使用することが理解されるだろう。この洞察を使用して、図5を参照して例示される修正において、これは、ステップ５０６において関数F（Ｍ、ｎ’max）を評価することによって有利に使用され得る。
【００７５】
【数１１】
2019521329.tif 000013
【００７６】
式３において、Ｍは、ステップ５０４において定義された範囲内の任意のゼロ荷電質量である。変数ｍaは、付加イオンの質量であり、Ｘは、ステップ５０３で評価された信号増強質量電荷比スペクトルの分布関数関数F（Ｍ、ｎ’max）は、最大荷電状態の値を最適化するために使用される。ステップ５０６から５０８でｎ’maxを循環し、ｎmaxまでの異なる最大荷電状態でＦ（Ｍ、ｎ’max）を評価することにより、関数F（Ｍ、ｎ’max）の最大値を決定することができる。この最適化の利点は２倍である。第１に、それは、特に、入力質量電荷比スペクトルの信号対雑音が悪い場合に、この方法が親分子の正確な質量を見つける最良の機会を有することを保証する。さらに、方法はｎmaxの固定値だが高い値を使用することができるため、ｎmaxを実験的またはユーザー定義の変数ではなく方法において定数に減少させることができる。
【００７７】
この反復的アプローチは、親分子のための最も適切な質量を有利に同定し、荷電質量スペクトルにおけるアーチファクトの存在を低減することができることが理解されよう。高分解能質量分析において、当業者は、荷電状態が同位体分離から容易に推測できることを理解するであろう。しかしながら、低分解能質量分析データを使用する場合、親質量を効果的に抽出するために、最大荷電状態（ｎmax）に関する仮定を行う必要がある。これは、推定が高すぎる場合には問題となることがあり、これは図４および図５を参照すると容易に理解することができる。質量Ｍ＝９００を有する親分子の多価質量スペクトルの非常に単純な例を示す図１５から１８を参照されたい。
【００７８】
この例では、最大電荷（ｎmax）は３に設定され、単純な強度の加算に基づいて従来技術ごとにＦ（Ｍ）が評価された場合、Ｆ（Ｍ）は１５の値に対応する（即ち、それぞれ５の強度で３つのピーク）。
【００７９】
ｎmaxがあまりにも高く設定された場合、例えばｎmax＝６である場合（図１６Ａ）、これは、ゼロ荷電質量に出現するアーチファクトピーク（Ｍ＝９００での実際のピークと同様の大きさを有する）を導入する効果を有する。強度の簡単な合計の場合、Ｆ（Ｍ）≒１５．４の出力は、Ｍの９００と１８００の２つの値（図１６Ｂ）で生成することができる。さらに、ｎmaxの値が高すぎると、ベースラインからの雑音の余計なサンプリングにつながる可能性がある。
【００８０】
したがって、本教示がなく、Ｆ（Ｍ）を計算するために強度の単純な加算を使用する技術を使用すると、結果の計算は、ｎmaxが増加するにつれてベースライン雑音のサンプリングを増加させるだけであることが理解されよう。これはゼロ荷電質量スペクトルにおける偽／アーチファクトピークの大きさの増加をもたらす。元の質量電荷比スペクトルが特に雑音が多い低分解能質量分析データに本発明の技術を適用した場合、これはこれらのスペクトルを分析する問題のみを複合化するだけである。したがって、ｎmaxを高く推定しすぎるとアーチファクトのピークにつながり、ｎmaxを低くと推定しすぎると親分子（Ｍ）の正確な質量値が失われる可能性があることは理解されよう。
【００８１】
値をｎmaxに間違って帰属させることによって生じるこれらの潜在的な問題を認識して、本発明者は、平均中心スペクトルのログを用いて、上記の式1に従って、雑音またはピーク間の領域を負にし、Ｆ（Ｍ）の最大値を小さくする。図１７に示すように、図１５に示した同じ単純化されたスペクトルの例を考慮するが、その代わりにｎmaxx＝３の式１の変換を使用すると、等価Ｆ（Ｍ）値は１５．４とは対照的に約９である。これらのピークの高さは、図１５の比較可能な結果と比較して低減されるが、依然として容易に同定可能である。
【００８２】
本発明の技法の意義は、Ｍ＝９００の場合にｎmax＝６でｎmaxを高く設定する例では、Ｍ＝１８００での結果として得られるアーチファクトで、より容易に明らかである。図１８Ａおよび図１８Ｂに示すように、本教示によるＦ（Ｍ）の結果は、Ｆ（Ｍ）を約１５．４とした対応する従来技術とは対照的に、約６である。以上のように、信号増強化質量電荷比スペクトル（式１）への変換を使用し、ｎmaxが高すぎる場合、これは、単純な合計の場合に見られるものとは対照的に、Ｆ（Ｍ）の最大値の減少をもたらす。したがって、ステップ５０６、５０７および５０８においてｎmaxを循環し、ｎmaxがＦ（Ｍ）の最高値に対応する電荷スペクトルを出力することによって、Ｆ（Ｍ）の最大値を探すことが可能である。これは、ｎmaxの初期推定およびユーザによるさらなる最適化を行う必要がない点で、従来技術の手法よりも有益である。ｎmaxの値は、最初は定数として任意に高い値に固定することがでる。ｎmaxの値は、親分子の正しい質量を得るために最適化することができます。図１５および図１６の例を参照して示されるように、この最適化は、生の質量電荷比スペクトルの単純な加算を用いては不可能である。ステップ５０６、５０７および５０８におけるこの最適化は、多くの場合、数千の要素を有する２Ｄアレイを作成するために、複数の充電状態にわたるすべてのデータ点を合計することを含むので、これは、以下に参照するような計算リソースを必要とする複雑な計算を表す図１４に示す。
【００８３】
Ｆ（Ｍ、ｎ’max）が最大であるｎ’maxの決定は、平均スペクトルの対数分析ー上記の式１の計算を利用することが理解されよう。商とログの使用は、いくつかの非常に強いピークが最終結果に不当に影響しないことを確実にする。実験質量スペクトルは、付加物または不純物から生じる強いピークを含むことが多いことが理解されよう。この商およびログ関数を使用することにより、多重電荷分子で起こるようなピークの集合が、本発明者が実現した単一の強烈なピークよりも影響が大きいことが確実になる。平均化または正規化されたスペクトルデータセットを生データサンプルから差し引いた。これは、図１５の単純化された質量電荷比スペクトルのデータの拡張から明らかである。この延長では、図１９に示すように、Ｍ＝９００（ｍ／ｚ＝９００、４５０および３００）の親の質量に無関係であるｍ／ｚ＝１５００の単一の強い不純物ピークが、典型的な質量スペクトルに存在する可能性のある不純物データの種類を例示する。この例では、図１９のデータは、ｍ／ｚ＝１５００で不純物を添加した生の質量電荷比スペクトルを表している。ここで、この単一のピークは、Ｍ＝９００（ｍ／ｚ＝９００、４５０および３００）での親分子からの３つのピークの合計強度（１５．５１）よりもはるかに強い（７４．９）ことが分かる。雑音技術の簡単な減算を使用すると、ｍ／ｚ＝１５００でのピークがｍ／ｚ＝９００、４５０および３００でのピークの複合強度よりも強いままである図２０Ａに図形的に表されるようなデータセットが得られる。対照的に、商およびｌｏｇ関数に基づく変換を有利に用いる式１の技法を用いると、ｍ／ｚ＝１５００での不純物ピークの強度は、我々の親分子（Ｍ＝９００）からのｍ／ｚ＝９００、４５０および３００でのピークの複合強度よりも小さくなる。
【００８４】
一旦Ｆ（Ｍ、ｎ’max）が最大であるｎ’maxの値が決定されると、ステップ５１０および５１１において、ゼロ荷電スペクトルを生成することができる。これらのステップの図１を参照して説明したものとの比較それらがステップ１０７および１０８と同一であることを確認する。これらの技術を使用して、図６に明らかな出力を提供することができる。図６の例は、図２のデータセットに対する図５の方法の適用を示しており、当然のことながら、図４と同じピークが特定されている。図４を参照して上述したように、この表現は、データが様々な異なる形態、例えば質量値の表の形で、または単に単一の質量値の出力として出力され得るため、限定的ではない。
【００８５】
関数Ｆ（Ｍ、ｎ’max）は、一価の場合のゼロ荷電式３から得られることが理解されよう。以下の式３ａに従うと、ステップ５１０および５１１において、一価スペクトルを出力するように関数Ｆ（Ｍ、ｎ’max）を評価し得る。
【００８６】
【数１２】
2019521329.tif 000014
【００８７】
より一般的には、関数Ｆ（Ｍ、ｎ’max）は、ゼロおよび一価スペクトルの両方について導出し得ることが理解されよう。以下の式３ｂのように、関数Ｆ（Ｍ、ｎ’max）は、ステップ５１０および５１１においてゼロまたは一価スペクトルを出力するように評価され得る。
【００８８】
【数１３】
2019521329.tif 000015
【００８９】
ｚが０に等しいとき、式３ｂの関数Ｆ（Ｍ、ｎ’max）は、式３のゼロ荷電関数に戻ることは明らかである。式３ａの一価の関数は、式３ｂの一般的な関数から得られ、ｚは１に等しい。
【００９０】
図１および図５の技術が図２のデータセットに適用された場合、類似の出力を提供するが、本発明者は、データセットの追加の前処理の形でこれらの技術をさらに変更することが採用され得ることが分かった。これまでに説明したものに対するこの変更例は、図７のフローチャートのステップのシーケンスを参照して提示される。図５を参照して上述した処理技術を共有するこの技法において、２つの追加のステップ−ステップ７０２および７０３−が使用される。これらのステップを図５を参照して上述した技術に含めることにより、元の質量電荷比スペクトルから質量情報を抽出する方法の強化がなされる。これらの前処理工程を使用するこの技術は、信号対雑音が低い質量電荷比スペクトルから親分子質量を抽出する際に特に有益である。例えば、図２に示すウシ血清アルブミンの質量電荷比スペクトルこのデータセットと図２のデータとの比較は、図２のデータセットよりも低いイオンカウントレベルが問題であり、したがって、固有の雑音レベルによるエラーの影響をより受けやすいことを示している。
【００９１】
図８に示されるようなデータセットが提供される場合、本教示は、ピーク構造を維持するが、スペクトル内の雑音を除去するための前処理平滑化ステップ（ステップ７０２）を用いる。好ましい構成では、Ｓａｖｉｔｚｋｙ−Ｇｏｌａｙフィルタが使用される。このフィルタは、質量電荷比スペクトルにおける元の特徴の形状を維持する上で非常に良好であるからである。
【００９２】
他の移動平均平滑法または指数平滑法を利用することもできる。図９は、図８のスペクトルに対するこの前処理平滑化の使用を示しており、これから「よりクリーンな」データセットを生成できることが明らかであろう。
【００９３】
図９のこのフィルタリングされたデータセットを使用し、図７の方法は、ステップ７０３で追加の処理ステップ−バックグラウンド減算−を提供する。凸包、ウェーブレット、およびメジアンフィルタのような様々なベースライン減算技術を利用することができるが、好ましい態様では、Ｎｕｃｌ．  Ｉｎｓｔｒｕｍ．  Ｍｅｔｈ．  Ｂ，第３４巻、３９６−４０２（１９９８）に示されるような統計的に感受性のある非線形反復ピーククリッピング（ＳＮＩＰ）ベースライン減算が使用される。ＳＮＩＰは、多種多様な背景形状に対処する能力が実証されているので、質量分析において広く使用されていることは、質量分析の当業者には理解されるであろう。図１０は、図９に提示された平滑化されたデータに対するこのベースライン減算の使用を示す。ベースライン減算の使用は、図９の元のウシ血清アルブミンスペクトル１００１に対して平滑化された質量電荷比スペクトル１００２を生成することが理解されよう。
【００９４】
ベースライン減算ステップ（７０３）は平滑化ステップに従う必要はなく、ステップ７０２および他の手法は平滑化ステップ（７０２）の前処理ステップとしてベースライン減算を使用することができることが理解されよう。代替技術は、ベースライン減算または平滑化を単独で利用することもできる。
【００９５】
ステップ７０１から７０３の後で、前処理されたスペクトルは、図１および図５のステップ１０２およびステップ５０２にそれぞれ概説されたものと同じ方法で、ステップ７０４において平均中心化される。図１１は、図１０の前処理されたデータに対するこのステップの効果を示している。図１０に関して前述したウシ血清アルブミンの前処理された質量電荷比スペクトル１００１は、図１１において１１０１で表され、１１０２は、ステップ７０４の平均中心化データ技術を使用するそのデータセットの処理を表す。
【００９６】
ステップ７０４に続いて、ステップ７０５において、信号増強化質量電荷比スペクトルが、上述のように、ステップ１０３および５０３と同一の方法で式１を用いて評価される。図１２は、図１１に示すデータに適用した場合のこの出力を示す。
【００９７】
ステップ７０６から７１２では、信号増強化質量電荷比スペクトルが、上で概説したステップ５０４から５１０と同じ方法で式３の関数を使用して評価される。
【００９８】
ステップ７１３において、ゼロ荷電スペクトルが生成される。図１３は、図７の技術が、図８に示すような雑音の多いデータセットに使用されたときのこの出力を示している。図４および図６の表現と同様に、グラフ形式のこのデータ表現は、これに限定されるものではなく、代替の構成では、この出力は、質量値の表または単一の質量値として提示することができる。
【００９９】
ステップ７０６において、信号増強化質量電荷比スペクトルは、ステップ５０４から５１０に関して上に概説したのと同じ方法で、式３ａまたは式３ｂの関数のいずれかを使用して評価し得ることが理解されよう。特定の方程式を選択すると、ステップ７１３においてゼロまたは一価スペクトルがそれぞれ生成されることが理解されよう。
【０１００】
上述のように、本教示による方法は、質量分析計によって生成されるデータセットを使用し、本教示の機能性は、そのような質量分析計の既存の機能性と統合され得る。本教示内で有用に使用され得る既知の質量分析計の機能の例には、欧州特許第１８６５５３３号または欧州特許第２３７２７４５号のような当初の出願に記載されているものが含まれる。これらのタイプの質量分析計を使用することにより、本教示によれば、追加の処理を統合することによって、または結合することによって、それらの機能性を同じものに拡張することが可能である。そのような統合の一例が図１４に示されており、質量分析計１４００は、エレクトロスプレイイオン化ソースであり得るイオン化ソース１４３０を含む。イオン化ソースは、検出器１４２０によって検出されたときに、上述の入力質量電荷比スペクトルを生成する多原子親分子をイオン化するように構成される。次いで、入力質量電荷比スペクトルは、コンピュータシステムまたは他の処理装置６００に中継される。一実装では、処理装置６００は、通常、少なくとも１つの処理部６０２およびメモリ６０４を含む。処理装置６００、メモリ６０４は、揮発性（例えば、ＲＡＭ）、不揮発性（例えば、ＲＯＭおよびフラッシュメモリ）、またはその両方の組み合わせであってもよい。処理装置６００の最も基本的な構成は、破線６０６によって示されるように、処理部６０２およびメモリ６０４のみを含む必要がある。主または基本的オペレーティングシステムは、不揮発性メモリ６０４において処理装置６００の基本機能を制御するように構成される。
【０１０１】
処理装置６００は、メモリ格納または検索のための追加の装置をさらに含み得る。これらのデバイスは、磁気および光学媒体上のメモリ格納および検索のための、例えば、メモリカード、磁気ディスクドライブ、磁気テープドライブ、および光学ドライブなどのリムーバブルストレージデバイス６０８または固定型ストレージデバイス６１０であってもよい。ストレージ媒体は、取り外し可能型および固定型の両方である、揮発性および不揮発性媒体を含み得、数の構成のいずれかで、例えば、ＲＡＭ、ＲＯＭ、ＥＥＰＲＯＭ、フラッシュメモリ、ＣＤ−ＲＯＭ、ＤＶＤ、またはその他の光学ストレージ媒体、または磁気カセット、磁気テープ、磁気ディスク、または他の磁気記憶装置、またはデータを格納するために使用することができ、処理部６０２によってアクセスすることができる任意の他のメモリ技術または媒体で提供され得る。特定用途向けの処理装置６００を作成するために基本オペレーティングシステムと相互作用する、例えば、ソフトウェアの形態の、追加の命令、この実装では、質量分析計検出器１４２０から受け取った質量分析計データを処理するためのデータアレイ形式の命令は、例えば、コンピュータ可読命令、データ構造、およびプログラムモジュールなどのデータの記憶のための任意の方法または技術を使用して、メモリ６０４または記憶装置６１０に格納し得る。
【０１０２】
処理装置６００は、処理装置６００が、例えば、質量分析計などの他の装置と通信することを可能にする１つ以上の通信インターフェース６１２を有してもよい。通信インターフェース６１２は、ネットワークに接続されてもよい。ネットワークは、ローカルエリアネットワーク（ＬＡＮ）、ワイドエリアネットワーク（ＷＡＮ）、電話ネットワーク、ケーブルネットワーク、光ネットワーク、インターネット、直接有線接続、無線ネットワーク、例えば、無線周波数、赤外線、マイクロ波、または音響、またはデバイス間のデータ転送を可能にする他のネットワークであってもよい。データは、一般に、変調されたデータ信号、例えば、搬送波または他の搬送媒体を介して、ネットワークを介して通信インターフェース６１２に、およびそこから送信される。変調されたデータ信号は、信号内のデータを符号化するように設定または変更することができる特性を有する電磁信号であることが理解されよう。このようにして、図１４の構成は、破線の輪郭１４００内の完全な統合を示しているが、処理６００の機能は、実際の検出器１４２０およびイオン源１４３０から遠隔で行い得る。
【０１０３】
処理装置６００はさらに、様々な入力デバイス６１４および出力デバイス６１６を有してもよい。例示的な入力デバイス６１４は、ビデオカメラ、レコーダまたは再生ユニット、キーボード、マウス、タブレット、および／またはタッチスクリーンデバイスを含み得る。例示的な出力デバイス６１６は、ビデオディスプレイ、オーディオスピーカ、および／またはプリンタを含み得る。そのような入力デバイス６１４および出力デバイス６１６は、コンピュータシステム６００と一体化されてもよく、またはそれらは有線または無線で、例えば、ＩＥＥＥ８０２．１１またはＢｌｕｅｔｏｏｔｈプロトコルを介してコンピュータシステム６００に接続されてもよい。これらの入力デバイスおよび出力デバイスは、質量分析計用のユーザインターフェース１４１０と通信するか、またはユーザインターフェース１４１０と一体化し得る。これらの一体化されたまたは周辺入出力デバイスは、一般に周知であり、本明細書ではこれ以上説明しない。ネットワーク通信トランザクションを処理するなどの他の機能は、処理装置６００の不揮発性メモリ６０４内のオペレーティングシステムによって実行されてもよい。
【０１０４】
本明細書で使用される場合、含む／含んでいるという用語は、記載された特徴、整数、ステップまたは構成要素の存在を特定するものであるが、１つまたは複数の他の特徴、整数、ステップ、構成要素またはそれらのグループの存在または追加を排除しない。
【０１０５】
本明細書に記載の方法は、質量分析計または質量分析計の構成要素に結合されているか、またはそれらと電子的に通信している１つまたは複数のシステムにおける論理演算および／またはモジュールとして実施し得る。論理演算は、１つまたは複数のコンピュータシステムで実行される一連のプロセッサ実装ステップとして、および１つまたは複数のコンピュータシステム内の相互接続された機械または回路モジュールとして実装されてもよい。
【０１０６】
同様に、様々な構成要素モジュールの説明は、モジュールによって実行または達成される動作に関して提供されてもよい。結果として生じる実装は、記述された技術を実装する基礎となるシステムの性能要件に依存する、選択の問題である。したがって、本明細書で説明される技術の実施形態を構成する論理演算は、動作、ステップ、オブジェクト、またはモジュールとして様々に参照される。さらに、別段明示的に請求されない限り、または特定の順序がクレーム言語によって本質的に必要とされない限り、論理演算は任意の順序で実行されてもよいことを理解されたい。
【０１０７】
いくつかの実装形態では、製品は、本発明を実施するためにコンピュータシステム上で動作の具体化を実行するコンピュータプログラム製品として提供される。コンピュータプログラム製品の１つの実装形態は、コンピュータプログラムを符号化するコンピュータシステムによって読み取り可能な非一時的なコンピュータプログラム記憶媒体を提供する。記載された技術は、パーソナルコンピュータとは独立した特定用途向けの装置で使用され得ることをさらに理解されたい。上記の明細書、実施例およびデータは、特許請求の範囲で画定される本発明の例示的な実施形態の構造および使用の完全な説明を提供する。特許請求の範囲に記載された発明の様々な実施形態および態様は、ある程度の特殊性をもって、または１つまたは複数の個別の実施形態を参照して上述されているが、当業者は、特許請求されている本発明の精神または技術範囲を逸脱することなく開示された実施形態に多数の変更を行うことができる。上記の説明に含まれ、添付の図面に示された全ての事項は、特定の実施形態の単なる例示として解釈されるものであり、限定的なものではないものとする。以下の特許請求の範囲で画定される本発明の基本的な要素から逸脱することなく、詳細なまたは構造の変形を実施し得る。
【図１】
2019521329.tif 000018
【図２】
2019521329.tif 000019
【図３】
2019521329.tif 000020
【図４】
2019521329.tif 000021
【図５】
2019521329.tif 000022
【図６】
2019521329.tif 000023
【図７】
2019521329.tif 000024
【図８】
2019521329.tif 000025
【図９】
2019521329.tif 000026
【図１０】
2019521329.tif 000027
【図１１】
2019521329.tif 000028
【図１２】
2019521329.tif 000029
【図１３】
2019521329.tif 000030
【図１４】
2019521329.tif 000031
【図１５】
2019521329.tif 000032
【図１６Ａ】
2019521329.tif 000033
【図１６Ｂ】
2019521329.tif 000034
【図１７】
2019521329.tif 000035
【図１８Ａ】
2019521329.tif 000036
【図１８Ｂ】
2019521329.tif 000037
【図１９】
2019521329.tif 000038
【図２０Ａ】
2019521329.tif 000039
【図２０Ｂ】
2019521329.tif 000040
【手続補正書】
【提出日】平成30年3月26日(2018.3.26)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
多価種の内の多原子親分子の質量を同定するために低分解能質量電荷比スペクトルから質量情報を抽出する方法であって、
多価イオンの集団を示すデータセットを質量分析計から受信するステップであり、各イオン上の電荷の数はそのイオンの荷電状態を規定し、前記イオンの集団内のイオンのサブ集団は同じ電荷状態を有する前記イオンの集団のイオンを含む、ステップと、
前記イオンの荷電状態を規定するデータセットを使用して入力質量電荷比スペクトルを生成するステップであり、前記各荷電状態のサブ集団は前記質量電荷比スペクトルにおける強度によって表される、ステップと、
前記入力質量電荷比スペクトルを処理して信号増強化質量電荷比スペクトル（１０３）を提供するステップであり、前記信号増強化質量電荷比スペクトルは前記質量電荷比スペクトルおよび前記質量電荷比スペクトルの平滑化された表現の商の対数から生成される、ステップと、
前記多原子親分子の質量をその中で探索するための規定された荷電質量スペクトルの規定された荷電質量値の範囲を同定する（１０４）ステップと、
前記規定された荷電質量値の範囲内の各質量に対し、最大荷電状態までの連続した荷電状態で前記質量に対応する前記信号増強化質量電荷比スペクトルにおける値の加算に等しい総和を関数
【数１】
2019521329.tif 000041
を使用して生成するステップであり、ここで、Ｍは規定された荷電質量値の前記範囲内の任意の規定された質量であり、ｍaは電荷付加物の質量であり、ｚは規定された荷電質量スペクトルの次数であり、Ｘは信号増強化質量電荷比スペクトルの分布関数である、ステップと、
前記総和の値を用いて前記多原子親分子の質量を決定するステップと
を含む方法。
【請求項２】
前記総和値を使用して多原子親分子の質量を決定するステップは、前記親分子の前記質量を決定するために、前記総和の範囲にわたって前記総和からの値を正規化するステップを含む、請求項１に記載の方法。
【請求項３】
前記電荷輸送付加物の質量、ｍa、は、プロトン質量を表す１に等しくなるように設定される、請求項１または２に記載の方法。
【請求項４】
前記信号増強化質量電荷比スペクトル（１０３）は関数
【数２】
2019521329.tif 000042
から生成され、
ここで、Ｘは前記信号増強化質量電荷比スペクトルを表し、Iは入力質量電荷比スペクトル（２０１）内の強度に対応し、Ｓは前記入力質量電荷比スペクトルの平滑化された表現である、請求項１乃至３のいずれか一項に記載の方法。
【請求項５】
前記規定された荷電質量値の前記規定された範囲内の各質量について、前記方法は、
最大荷電状態の異なる値までの合計を計算するステップと、
前記最大荷電状態の異なる値までの前記合計からの前記値を使用して前記親分子の分子量を決定するステップと
をさらに含む、請求項１乃至４のいずれか一項に記載の方法。
【請求項６】
前記信号増強化質量電荷比スペクトルを形成する前に、前記入力質量電荷比スペクトルのスペクトル前処理をさらに含む、請求項１乃至５のいずれか一項に記載の方法。
【請求項７】
前記スペクトル前処理が、平滑化および／またはベースライン減算のうちの少なくとも１つから選択される、請求項６に記載の方法。
【請求項８】
前記多価種がエレクトロスプレイイオン化によって生成される、請求項１乃至７のいずれか一項に記載の方法。
【請求項９】
前記入力質量電荷比スペクトルを生成するために、質量分析計源内で多原子親分子をイオン化するステップを含む、請求項１乃至８のいずれか一項に記載の方法。
【請求項１０】
前記規定された電荷対質量スペクトルは、全ての同定された種がゼロ荷電値を有するゼロ荷電質量スペクトルである、請求項１乃至９のいずれか一項に記載の方法。
【請求項１１】
前記規定された荷電質量スペクトルが、一価質量値の範囲を同定するために使用される一価質量スペクトルである、請求項１から９のいずれか一項に記載の方法。
【請求項１２】
エレクトロスプレイイオン化ソースおよび検出器を備える質量分析器であって、前記検出器は前記イオン化ソースによって生成された多価イオンの集団を示す出力を生成するように構成され、各イオン上の電荷数はそのイオンの荷電状態を規定し、各荷電状態は前記イオンの前記集団内のイオンのサブ集団からなり、前記分析器は請求項１乃至１１のいずれか一項に記載の方法を実施するように構成されたプロセッサをさらに備える、質量分析器。
【国際調査報告】
2019521329.tif 000043
2019521329.tif 000044