(57)【特許請求の範囲】
【請求項１】
多数の既知化合物のＭＳｎマススペクトルを予め記憶するデータベース記憶領域と、
未知化合物のＭＳｎマススペクトルを取得する取得部と、
未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を示すスコアをそれぞれ求めるスコア算出部とを備えるデータ処理装置であって、
未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度の大きさによって数段階の強度スコアにそれぞれ分類する強度スコア算出部と、
未知化合物のＭＳｎマススペクトル中のニュートラルロスのピークを求めるニュートラルロス算出部と、
既知化合物のＭＳｎマススペクトル中の前記ニュートラルロスのピークを含むピークの質量電荷比と、当該ピークの質量電荷比の精度範囲又は当該ピークの質量電荷比に最も近くに存在する未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との差に基づいて、質量電荷比の誤差を示す位置スコアを各ピークについてそれぞれ求める位置スコア算出部とを備え、
前記スコア算出部は、各ピークについて、既知化合物のＭＳｎマススペクトル中のピークの強度スコアと、未知化合物のＭＳｎマススペクトル中のピークの強度スコアと、位置スコアとを積算し、全ピークについての積算値の総和をスコアとして算出することを特徴とするデータ処理装置。
【請求項２】
既知化合物のＭＳｎマススペクトル中のピークの質量電荷比と、未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との内で、一致したペア数から求められた一致度に基づいて、一致度を示すピークスコアを求めるピークスコア算出部とを備え、
前記スコア算出部は、前記スコアに前記ピークスコアを加算することを特徴とする請求項１に記載のデータ処理装置。
【請求項３】
前記強度スコア算出部は、未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度又は当該イオン強度のｌｏｇの変換値の大きさによって３段階、もしくは数段階の強度スコアにそれぞれ分類することを特徴とする請求項１又は請求項２に記載のデータ処理装置。
【請求項４】
前記強度スコア算出部は、未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度の大きさとピークの質量電荷比の大きさとによって数段階の強度スコアにそれぞれ分類することを特徴とする請求項１又は請求項２に記載のデータ処理装置。
【請求項５】
前記強度スコア算出部は、未知化合物のＭＳｎマススペクトルにおいて、イオン強度の大きさが閾値以下であるピークを除去することを特徴とする請求項１又は請求項２に記載のデータ処理装置。
【請求項６】
前記位置スコア算出部は、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比を基準とし、当該基準から離れるにしたがって低くなる位置スコアを、未知化合物のＭＳｎマススペクトル中の各ピークについてそれぞれ求めることを特徴とする請求項１又は請求項２に記載のデータ処理装置。
【請求項７】
前記位置スコア算出部は、未知化合物のＭＳｎマススペクトル中のピークの質量電荷比を、アダクトイオンの質量を用いて補正することを特徴とする請求項１又は請求項２に記載のデータ処理装置。
【請求項８】
前記ピークスコア算出部は、一致したペア数にしたがって増加するピークスコアを求めることを特徴とする請求項２に記載のデータ処理装置。
【請求項９】
前記スコア算出部は、未知化合物のＭＳｎ＋ｍマススペクトルと、多数の既知化合物のＭＳｎ＋ｍマススペクトルとの類似性を示すＭＳｎ＋ｍマススペクトルのスコアをそれぞれ求め、
ＭＳｎマススペクトルのスコアに、ＭＳｎ＋ｍマススペクトルのスコアを加算することを特徴とする請求項１又は請求項２に記載のデータ処理装置。
【請求項１０】
多数の既知化合物のＭＳｎマススペクトルを予め記憶するデータベース記憶領域と、
未知化合物のＭＳｎマススペクトルを取得する取得部と、
未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を示すスコアをそれぞれ求めるスコア算出部とを備えるデータ処理装置を用いたデータ処理方法であって、
未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度の大きさによって数段階の強度スコアにそれぞれ分類する強度スコア算出ステップと、
未知化合物のＭＳｎマススペクトル中のニュートラルロスのピークを求めるニュートラルロス算出ステップと、
既知化合物のＭＳｎマススペクトル中の前記ニュートラルロスのピークを含むピークの質量電荷比と、当該ピークの質量電荷比の精度範囲又は当該ピークの質量電荷比に最も近くに存在する未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との差に基づいて、質量電荷比の誤差を示す位置スコアを各ピークについてそれぞれ求める位置スコア算出ステップと、
各ピークについて、既知化合物のＭＳｎマススペクトル中のピークの強度スコアと、未知化合物のＭＳｎマススペクトル中のピークの強度スコアと、位置スコアとを積算し、全ピークについての積算値の総和をスコアとして算出するスコア算出ステップとを含むことを特徴とするデータ処理方法。
【請求項１１】
既知化合物のＭＳｎマススペクトル中のピークの質量電荷比と、未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との内で、一致したペア数から求められた一致度に基づいて、一致度を示すピークスコアを求めるピークスコア算出ステップと、
前記スコアに前記ピークスコアを加算する加算ステップとを含むことを特徴とする請求項１０に記載のデータ処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、質量分析装置によって得られた未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を示すスコアをそれぞれ求め、そのスコアに基づいて未知化合物を同定するデータ処理方法及びそれを用いたデータ処理装置に関する。
【背景技術】
【０００２】
イオントラップ型質量分析装置等を用いた質量分析においては、ＭＳ／ＭＳ分析（ＭＳ２分析）という手法が知られている。一般的なＭＳ２分析では、まず分析対象の未知化合物から目的とする特定の質量電荷比（ｍ／ｚ）を有するイオンをプリカーサイオン（親イオン）として選別し、その選別したプリカーサイオンをＣＩＤ（Collision Induced Dissociation：衝突誘起分解）によって開裂させ、開裂イオンを生成する。開裂様式は化学構造に依存するため、開裂によって生成したフラグメントイオンを質量分析することにより、未知化合物のＭＳ２マススペクトル（化学構造の情報）を取得することができる。【０００３】
近年、医薬品や農薬やタンパク質等の未知化合物を同定する場合には、未知化合物の分子量が大きくなり、化学構造も複雑になっている。そのため、未知化合物の種類によっては、一段階の開裂操作だけでは充分に小さな質量までイオンが開裂しない場合がある。そうした場合には、開裂操作を複数回繰り返し、生成したフラグメントイオンを質量分析するＭＳｎ分析が行われている。【０００４】
このようなＭＳｎ分析により得られたＭＳｎマススペクトルから未知化合物の化学構造を推定する解析処理として、データベースを用いたパターンマッチング処理が行われている（例えば、特許文献１や非特許文献１参照）。データベースには、多数の化合物名及び化学構造（既知化合物）と、その既知化合物のＭＳｎマススペクトルが含まれている。そして、パターンマッチング処理では、既知化合物のＭＳｎマススペクトルと、未知化合物のＭＳｎマススペクトルとを比較することにより、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を示すスコアを、重み付き内積やピークの出現確率等を用いてそれぞれ求めている。その結果、スコアが高い順番に既知化合物を配列することで、測定者は未知化合物を同定している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００５−２０１８３５号公報
【非特許文献】
【０００６】
【非特許文献１】Stephen E. Stein, and Donald R. Scott "Optimization and Testing of Mas Spectral Library Search Algorithms for Compound Identification" 1994, American Society for Mass Spectrometry, 859-866.
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、同一の化合物をＭＳｎ分析した場合でも、得られたＭＳｎマススペクトルが異なることがある。すなわち、ＭＳｎマススペクトルでは、ピークのイオン強度Ｉやピークの位置（ｍ／ｚ）やピークの有無等が変化しやすい。そのため、上述したようなパターンマッチング処理では、イオン強度Ｉが大きく変化したり、ピークの位置がずれたり、あるべきピークが現れなかったりすることがあり、スコアが低下し、その結果、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を正しく評価していないことがあった。
【課題を解決するための手段】
【０００８】
本件発明者は、上記課題を解決するために、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を正しく評価する検索方法について検討を行った。上述したようなパターンマッチング処理では、ＭＳｎマススペクトル中のピークのイオン強度Ｉやピークの位置（ｍ／ｚ）やピークの有無等が変化しやすいことが考慮されておらず、スコアが低下することがわかった。そこで、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を評価する際に、ＭＳｎマススペクトル中のピークのイオン強度Ｉやピークの位置（ｍ／ｚ）やピークの有無等が変化することを考慮した検索方法を見出した。【０００９】
すなわち、本発明のデータ処理装置は、多数の既知化合物のＭＳｎマススペクトルを予め記憶するデータベース記憶領域と、未知化合物のＭＳｎマススペクトルを取得する取得部と、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を示すスコアをそれぞれ求めるスコア算出部とを備えるデータ処理装置であって、未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度の大きさによって数段階の強度スコアにそれぞれ分類する強度スコア算出部と、未知化合物のＭＳｎマススペクトル中のニュートラルロスのピークを求めるニュートラルロス算出部と、既知化合物のＭＳｎマススペクトル中の前記ニュートラルロスのピークを含むピークの質量電荷比と、当該ピークの質量電荷比の精度範囲又は当該ピークの質量電荷比に最も近くに存在する未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との差に基づいて、質量電荷比の誤差を示す位置スコアを各ピークについてそれぞれ求める位置スコア算出部とを備え、前記スコア算出部は、各ピークについて、既知化合物のＭＳｎマススペクトル中のピークの強度スコアと、未知化合物のＭＳｎマススペクトル中のピークの強度スコアと、位置スコアとを積算し、全ピークについての積算値の総和をスコアとして算出するようにしている。
【００１０】
ここで、「未知化合物のＭＳｎマススペクトル」とは、未知化合物をＭＳｎ分析した際に得られたものであり、「既知化合物のＭＳｎマススペクトル」とは、未知化合物と同様の条件（イオン化法や開裂方法等）で既知化合物をＭＳｎ分析した際に得られたものであってもよく、また、未知化合物と異なる条件で既知化合物をＭＳｎ分析した際に得られたものや、既知化合物の化学構造から計算されたものであってもよい。【００１１】
本発明のデータ処理装置によれば、まず、＜１＞ピークのイオン強度Ｉについては次のように処理する。ピークのイオン強度Ｉは変化しやすいが、大きなイオン強度Ｉのピークは、対象化合物から得られやすいフラグメントイオンのピークであると考えられる。そこで、未知化合物のＭＳｎマススペクトルについて、イオン強度Ｉの大きさによって数段階（例えば、上位、中位、下位）に各ピークをそれぞれ分類するとともに、データベースの既知化合物のＭＳｎマススペクトルについても、イオン強度Ｉの大きさによって数段階（例えば、上位、中位、下位）に各ピークをそれぞれ分類する。これにより、イオン強度Ｉが異なっても、類似したピークを見つけるようにする。
また、ピークの有無は変化しやすいので、未知化合物のＭＳｎマススペクトルと既知化合物のＭＳｎマススペクトルとで同一の位置（質量電荷比）にピークが存在しなくても、互いにニュートラルロスの関係にあるピークが存在する場合がある。そこで、ニュートラルロス算出部は、未知化合物のＭＳｎマススペクトルについて、ピークの位置とプリカーサイオンの質量とからニュートラルロスのピークを求める。これにより、未知化合物のＭＳｎマススペクトル中であるべきピークが現れなくても、ピークを対応付けることができる。なお、ＭＳｎマススペクトル中に追加するニュートラルロスのピークのイオン強度Ｉについては、例えば、元としたピークと同じイオン強度Ｉとする。【００１２】
次に、＜２＞ピークの位置（ｍ／ｚ）については次のように処理する。未知化合物のＭＳｎマススペクトル中のピークの位置はある程度のずれが生じる。そこで、既知化合物のＭＳｎマススペクトル中のピークの位置を基準とし、質量分析装置のピークの位置の精度範囲（例えば、質量分析装置で−１．０〜１．０や、ＴＯＦで−０．００５〜０．００５や、タンデム四重極で−０．５〜０．５）に位置スコア「１．０」を対応付け、その精度範囲外で位置スコア「０.０」に漸近する位置スコアを対応付けることで、ピークの位置が異なっても、類似したピークを見つけるようにする。【００１３】
そして、上述したように、対応付けられた各ピークについて、既知化合物のＭＳｎマススペクトル中のピークの強度スコアと、未知化合物のＭＳｎマススペクトル中のピークの強度スコアと、位置スコアとを積算し、全ピークについての積算値の総和をスコアとして算出する。
【発明の効果】
【００１４】
以上のように、本発明のデータ処理装置によれば、未知化合物のＭＳｎマススペクトル中のピークのイオン強度Ｉやピークの位置（ｍ／ｚ）やピークの有無等が変化しても、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を正しく評価することができる。【００１５】
（その他の課題を解決するための手段及び効果）
また、本発明のデータ処理装置においては、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比と、未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との内で、一致したペア数から求められた一致度に基づいて、一致度を示すピークスコアを求めるピークスコア算出部とを備え、前記スコア算出部は、前記スコアに前記ピークスコアを加算するようにしてもよい。【００１６】
本発明のデータ処理装置によれば、＜３＞ピーク・ニュートラルロスペアの一致度については次のように処理する。既知化合物のＭＳｎマススペクトル中のピークの質量電荷比と、未知化合物のＭＳｎマススペクトル中のピーク（ニュートラルロスのピークも含む）の質量電荷比との内で、一致したペア数から一致度を示すピークスコア（例えば、（一致したピーク数／既知化合物のＭＳｎマススペクトル中のピーク数）×１００）を求める。【００１７】
また、本発明のデータ処理装置においては、前記強度スコア算出部は、未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度又は当該イオン強度のｌｏｇの変換値の大きさによって３段階、もしくは数段階の強度スコアにそれぞれ分類するようにしてもよい。【００１８】
また、本発明のデータ処理装置においては、前記強度スコア算出部は、未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度の大きさとピークの質量電荷比の大きさとによって数段階の強度スコアにそれぞれ分類するようにしてもよい。
本発明のデータ処理装置によれば、ＭＳｎマススペクトルでは、質量電荷比に比例してピークのユニーク度が大きくなるが、質量電荷比が大きなピークは、イオン強度が小さくなる場合があるため、イオン強度を質量電荷比に応じて重み付ける（例えばイオン強度をＩ、質量荷電比をＭとして、あるピークの強度ＩをＩ＝Ｍ×（Ｉ）１／２のように重み付ける）ことで、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を正しく評価することができる。【００１９】
また、本発明のデータ処理装置においては、前記強度スコア算出部は、未知化合物のＭＳｎマススペクトルにおいて、イオン強度の大きさが閾値以下であるピークを除去するようにしてもよい。
本発明のデータ処理装置によれば、適切なピークだけに注目したスコアを得ることができ、信頼性の高い検索ができるようになる。【００２０】
また、本発明のデータ処理装置においては、前記位置スコア算出部は、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比を基準とし、当該基準から離れるにしたがって低くなる位置スコアを、未知化合物のＭＳｎマススペクトル中の各ピークについてそれぞれ求めるようにしてもよい。
本発明のデータ処理装置によれば、例えば、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比を基準とし、その基準から離れるにしたがって低くなる適当な分布関数や、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比を基準とし、適当な傾きで位置スコア「０」に減衰していく近似関数や、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比を基準とし、一定範囲まである位置スコアとし、その一定範囲外は位置スコア「０」とするステップ関数や、これらの組み合わせにすることで、各ピークについて位置スコアをそれぞれ求める。なお、基準とする既知化合物のＭＳｎマススペクトル中のピークの質量電荷比も、一定値とするだけでなく、存在確率や、構造式の分子軌道計算から求められる結合力等に基づいた可変値としてもよい。また、分布関数や近似関数やステップ関数中に、既知化合物のＭＳｎマススペクトル中のピークが複数存在する場合には、最も基準に近いピークや、強度スコアと位置スコアとの積算値が高いピークを選択することとしてもよい。
以上のように、本発明のデータ処理装置によれば、ピークの位置の関数を用いることで、信頼性の高い検索ができるようになる。また、ピークの位置の関数を直線的な分布にした場合には、簡素化され効率的に検索ができるようになる。【００２１】
また、本発明のデータ処理装置においては、前記位置スコア算出部は、未知化合物のＭＳｎマススペクトル中のピークの質量電荷比を、アダクトイオンの質量を用いて補正するようにしてもよい。
本発明のデータ処理装置によれば、アダクトイオン（−Ｈ、−ＯＨ、−ＣＯ２等の付加体）によってピークの位置が１以上ずれることも考えられる。そこで、精度範囲外のピークに対してもアダクトに基づいて対応付けることで、ピークの位置が１以上ずれても、類似したピークを見つけることができる。【００２２】
また、本発明のデータ処理装置においては、前記ピークスコア算出部は、一致したペア数にしたがって増加するピークスコアを求めるようにしてもよい。
本発明のデータ処理装置によれば、一致したペア数にしたがってスコアへの影響度を変えることができ、多くの部分構造が一致する同一の化合物の代謝物等の変化体を正しく検索することができる。【００２３】
また、本発明のデータ処理装置においては、前記スコア算出部は、未知化合物のＭＳｎ＋ｍマススペクトルと、多数の既知化合物のＭＳｎ＋ｍマススペクトルとの類似性を示すＭＳｎ＋ｍマススペクトルのスコアをそれぞれ求め、ＭＳｎマススペクトルのスコアに、ＭＳｎ＋ｍマススペクトルのスコアを加算するようにしてもよい。
本発明のデータ処理装置によれば、関連付けられた子や親のマススペクトルもスコアに反映でき、信頼性の高い検索ができるようになる。【００２４】
そして、本発明のデータ処理方法は、多数の既知化合物のＭＳｎマススペクトルを予め記憶するデータベース記憶領域と、未知化合物のＭＳｎマススペクトルを取得する取得部と、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を示すスコアをそれぞれ求めるスコア算出部とを備えるデータ処理装置を用いたデータ処理方法であって、未知化合物のＭＳｎマススペクトル中の各ピーク及び既知化合物のＭＳｎマススペクトル中の各ピークを、イオン強度の大きさによって数段階の強度スコアにそれぞれ分類する強度スコア算出ステップと、未知化合物のＭＳｎマススペクトル中のニュートラルロスのピークを求めるニュートラルロス算出ステップと、既知化合物のＭＳｎマススペクトル中の前記ニュートラルロスのピークを含むピークの質量電荷比と、当該ピークの質量電荷比の精度範囲又は当該ピークの質量電荷比に最も近くに存在する未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との差に基づいて、質量電荷比の誤差を示す位置スコアを各ピークについてそれぞれ求める位置スコア算出ステップと、各ピークについて、既知化合物のＭＳｎマススペクトル中のピークの強度スコアと、未知化合物のＭＳｎマススペクトル中のピークの強度スコアと、位置スコアとを積算し、全ピークについての積算値の総和をスコアとして算出するスコア算出ステップとを含むようにしている。
【００２５】
また、本発明のデータ処理方法においては、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比と、未知化合物のＭＳｎマススペクトル中のピークの質量電荷比との内で、一致したペア数から求められた一致度に基づいて、一致度を示すピークスコアを求めるピークスコア算出ステップと、前記スコアに前記ピークスコアを加算する加算ステップとを含むようにしてもよい。
【図面の簡単な説明】
【００２６】
【図１】本発明の一実施形態である質量分析装置の概略構成を示すブロック図。
【図２】分布関数記憶領域に記憶された分布関数の一例を示す図。
【図３ａ】データ処理方法の一例について説明するためのフローチャート。
【図３ｂ】データ処理方法の一例について説明するためのフローチャート。
【発明を実施するための形態】
【００２７】
以下、本発明の実施形態について図面を用いて説明する。なお、本発明は、以下に説明するような実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々の態様が含まれる。
【００２８】
図１は、本発明の一実施形態である質量分析装置の概略構成を示すブロック図である。
液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１は、液体クロマトグラフ装置（ＬＣ）１０と、ＭＳｎ分析が可能な質量分析装置（ＭＳ）２０と、コンピュータ（データ処理装置）３０とを備える。
なお、ここでは、ＬＣ／ＭＳ１で得られたマススペクトルを利用する場合を例にとるが、ガスクロマトグラフ質量分析装置（ＧＣ／ＭＳ）等の他のクロマトグラフ質量分析装置や、試料の直接導入を行う質量分析装置を利用する場合も同様である。
【００２９】
ＬＣ１０は、カラムオーブン１１と、カラムオーブン１１に内装されるカラム１２と、カラム１２の入口端に接続されるサンプル注入部１３とを備える。
このようなＬＣ１０によれば、試料は、キャリアガスに押されてサンプル注入部１３からカラム１２内に導入されることになる。これにより、試料に含まれる各成分（未知化合物）は、カラム１２内を通過する間に時間軸方向に分離されて、カラム１２の出口端に到達することになる。なお、カラム１２の出口端はＭＳ２０に接続されている。
【００３０】
ＭＳ２０は、未知化合物をイオン化するイオン源２１と、生成したイオンを質量分離しかつＭＳｎ分析が可能な質量分離部２２と、質量分離されたイオンを検出する検出器２３とを備える。イオン源２１には、エレクトロスプレイイオン源の他、ソニックスプレイイオン源や、イオンスプレイや、マトリックス支援レーザー脱離イオン源等が使用可能である。質量分離部２２としては、例えば、三連四重極型の質量分析装置やイオントラップ型の質量分析装置等を用いることができる。なお、各試料成分由来のイオンから適当な質量電荷比（ｍ／ｚ）を有するイオンがプリカーサイオンとして自動的に選択されるか、事前に指定した質量荷電比（ｍ／ｚ）のピークをプリカーサイオンとして選定し、このプリカーサイオンを開裂させて発生したフラグメントイオンの質量分離・検出が行われる。
【００３１】
このようなＭＳ２０によれば、イオン源２１でイオン化されたイオンは、質量分離部２２へ導入される。質量分離部２２では、イオンを質量分離する。また、測定者の設定にしたがってＭＳｎ分析（ｎ＝２，３，４，・・・）を順次行う。質量分離されたイオンは検出器２３へ送られ、ＭＳｎマススペクトル（ｎ＝２，３，４，・・・）として検出され、ＭＳｎマススペクトル（ｎ＝２，３，４，・・・）のデータがコンピュータ３０へ送られる。
【００３２】
コンピュータ３０は、ＣＰＵ３１とメモリ３２と入力装置３３と表示装置３４とを備える。ＣＰＵ３１が処理する機能をブロック化して説明すると、未知化合物のＭＳｎマススペクトルのデータを取得する取得部３１ａと、強度スコアＩｄ，Ｉｑを算出する強度スコア算出部３１ｂと、位置スコアＳ（ｅｒｒ）を算出する位置スコア算出部３１ｃと、ニュートラルロス算出部３１ｅと、ピークスコアＳ（Ｎ）を算出するピークスコア算出部３１ｆと、スコアＳｃｒを算出するスコア算出部３１ｄとを有する。
【００３３】
メモリ３２は、イオン強度Ｉを蓄積するためのイオン強度記憶領域３２ａと、未知化合物のＭＳｎマススペクトルを記憶するための未知化合物データ記憶領域３２ｂと、多数の既知化合物のＭＳｎマススペクトルを予め記憶するデータベース記憶領域３２ｃと、位置スコアＳ（ｅｒｒ）を算出するための分布関数を予め記憶する分布関数記憶領域３２ｄと有する。
【００３４】
ここで、データベース記憶領域３２ｃに記憶された既知化合物のＭＳｎマススペクトルは、ＭＳ２０で実際に既知化合物をＭＳｎ分析した際に得られたものであり、多数の既知化合物をＭＳｎ分析することにより、データベース記憶領域３２ｃに記憶されている。
そして、データベース記憶領域３２ｃに記憶された既知化合物のＭＳｎマススペクトルは、イオン強度Ｉｄについては強度スコア算出部３１ｂによりイオン強度Ｉｄの大きさによって３段階に各ピークがそれぞれ分類されており、具体的には、イオン強度ＩｄがＩ１以上であるピークを上位の強度スコア「４」に分類し、イオン強度ＩｄがＩ２以上Ｉ１未満であるピークを中位の強度スコア「３」に分類し、イオン強度ＩｄがＩ２未満であるピークを下位の強度スコア「２」に分類している。
また、既知化合物のＭＳｎマススペクトルは、質量電荷比（ｍ／ｚ）については位置スコア算出部３１ｃにより各ピークの質量電荷比（ｍ／ｚ）がそれぞれ求められている。
これにより、例えば、一の既知化合物のＭＳｎマススペクトルでは、ｉ番目のピークについては強度スコア「４」、質量電荷比（ｍｉ／ｚｉ）が割り当てられ、（ｉ＋１）番目のピークについては強度スコア「２」、質量電荷比（ｍ（ｉ＋１）／ｚ（ｉ＋１））が割り当てられるように、全てのピークについて強度スコアと質量電荷比（ｍ／ｚ）とが割り当てられている。そして、全ての既知化合物のＭＳｎマススペクトルについて、全てのピークについて強度スコアと質量電荷比（ｍ／ｚ）とが割り当てられている。
【００３５】
また、図２は、分布関数記憶領域３２ｄに記憶された分布関数の一例を示す図である。横軸は質量電荷比ｍ／ｚの差ｅｒｒとなっており、縦軸は位置スコアＳ（ｅｒｒ）となっている。分布関数は、差ｅｒｒが０であるときには位置スコア「１」を対応付け、差ｅｒｒが大きくなるにしたがって位置スコア「０」に近づくように低くなっている。
【００３６】
取得部３１ａは、未知化合物をＭＳｎ分析（ｎ＝２，３，４，・・・）することで、検出器２３で取得されたイオン強度Ｉｑをメモリ３２に蓄積させた後、イオン強度Ｉｑを縦軸に、ｍ／ｚを横軸にとることにより、ＭＳｎマススペクトル（ｎ＝２，３，４，・・・）を作成して未知化合物データ記憶領域３２ｂに記憶させる制御を行う。
【００３７】
強度スコア算出部３１ｂは、未知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・・）中の各ピークを、イオン強度Ｉｑの大きさによって３段階の強度スコアにそれぞれ分類する制御を行う。具体的には、イオン強度ＩｑがＩ１以上であるピークを上位の強度スコア「４」に分類し、イオン強度ＩｑがＩ２以上Ｉ１未満であるピークを中位の強度スコア「３」に分類し、イオン強度ＩｑがＩ２未満であるピークを下位の強度スコア「２」に分類する。
【００３８】
ニュートラルロス算出部３１ｅは、未知化合物のＭＳｎマススペクトル中のニュートラルロスのピークを求める制御を行う。例えば、各試料成分由来のイオンから適当な質量電荷比（ｍ／ｚ）を有するイオンがプリカーサイオンとして自動的に選択されるか、事前に指定した質量荷電比（ｍ／ｚ）のピークをプリカーサイオンとして選定し、得られた未知化合物のＭＳｎマススペクトルについて、ピークの質量電荷比（ｍ／ｚ）とプリカーサイオンの質量とからニュートラルロスを求める。これにより、未知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・・）中にニュートラルロスのピークを含んだものが作成される。なお、ＭＳｎマススペクトル中に追加するニュートラルロスのピークのイオン強度Ｉについては、例えば、元としたピークと同じイオン強度Ｉとする。
【００３９】
位置スコア算出部３１ｃは、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比（ｍ／ｚ）と、対応する未知化合物のＭＳｎマススペクトル中のピーク（ニュートラルロスのピークも含む）の質量電荷比（ｍ／ｚ）との差ｅｒｒに基づいて、質量電荷比（ｍ／ｚ）の誤差を示す位置スコアＳ（ｅｒｒ）を各ピークについてそれぞれ求める制御を行う。
まず、未知化合物データ記憶領域３２ｂに記憶された未知化合物のＭＳｎマススペクトルについて、各ピーク（ニュートラルロスのピークも含む）の質量電荷比（ｍ／ｚ）をそれぞれ求める。
次に、未知化合物のＭＳｎマススペクトルと、様々な既知化合物のＭＳｎマススペクトルとを次々と比較していくことになるが、Ｘ番目の既知化合物のＭＳｎマススペクトルと、未知化合物のＭＳｎマススペクトルとを比較するときには、Ｘ番目の既知化合物のＭＳｎマススペクトル中のｉ番目のピークの質量電荷比（ｍ／ｚ）に最も近い質量電荷比（ｍ／ｚ）のピーク（ニュートラルロスのピークも含む）を見つけ出す。例えば、Ｘ番目の既知化合物のＭＳｎマススペクトル中のｉ番目のピークの質量電荷比（ｍ／ｚ）に最も近い質量電荷比（ｍ／ｚ）のピークとして、未知化合物のＭＳｎマススペクトル中のｊ番目のピークを見つけ出す。そして、Ｘ番目の既知化合物のＭＳｎマススペクトル中のｉ番目のピークの質量電荷比（ｍ／ｚ）と、未知化合物のＭＳｎマススペクトル中のｊ番目のピークの質量電荷比（ｍ／ｚ）との差ｅｒｒを算出し、差ｅｒｒを図２に示す分布関数に代入することで、位置スコアＳｉｊ（ｅｒｒ）を求める。
【００４０】
さらに、Ｘ番目の既知化合物のＭＳｎマススペクトル中の（ｉ＋１）番目のピークの質量電荷比（ｍ／ｚ）に最も近い質量電荷比（ｍ／ｚ）のピーク（ニュートラルロスのピークも含む）を見つけ出す。例えば、Ｘ番目の既知化合物のＭＳｎマススペクトル中の（ｉ＋１）番目のピークの質量電荷比（ｍ／ｚ）に最も近い質量電荷比（ｍ／ｚ）のピークとして、未知化合物のＭＳｎマススペクトル中の（ｊ＋ｎ）番目のピークを見つけ出す。そして、Ｘ番目の既知化合物のＭＳｎマススペクトル中の（ｉ＋１）番目のピークの質量電荷比（ｍ／ｚ）と、未知化合物のＭＳｎマススペクトル中の（ｊ＋ｎ）番目のピークの質量電荷比（ｍ／ｚ）との差ｅｒｒを算出し、差ｅｒｒを図２に示す分布関数に代入することで、位置スコアＳ（ｉ＋１）（ｊ＋ｎ）（ｅｒｒ）を求める。
このようにして、Ｘ番目の既知化合物のＭＳｎマススペクトル中の全ピークについて対応する未知化合物のＭＳｎマススペクトル中のピーク（ニュートラルロスのピークも含む）を見つけ出し、位置スコアＳ（ｅｒｒ）を求める。
【００４１】
ピークスコア算出部３１ｆは、下記式（１）に基づいてピークスコアＳ（Ｎ）を求める制御を行う。
Ｓ（Ｎ）＝（一致したピーク数／既知化合物のＭＳｎマススペクトル中のピーク数）×１００  ・・・（１）
これにより、既知化合物のＭＳｎマススペクトル中のピークの質量電荷比（ｍ／ｚ）と、未知化合物のＭＳｎマススペクトル中のピーク（ニュートラルロスのピークも含む）の質量電荷比（ｍ／ｚ）との内で、一致したペア数から求められた一致度を示すピークスコアＳ（Ｎ）が求められる。例えば、Ｘ番目の既知化合物のＭＳｎマススペクトル中のピーク数が５個であり、一致したピーク数が３個であれば、一致度を示すピークスコア「６０」となり、（Ｘ＋１）番目の既知化合物のＭＳｎマススペクトル中のピーク数が１０個であり、一致したピーク数が１個であれば、一致度を示すピークスコア「１０」となる。
なお、「一致した」とは、既知化合物のＭＳｎマススペクトルのピークに対して、未知化合物のＭＳｎマススペクトル中において、位置スコアＳ（ｅｒｒ）が０でないピーク或いはニュートラルロスのピークが存在することをいう。
【００４２】
スコア算出部３１ｄは、下記式（２）に基づいてスコアＳｃｒを算出する制御を行う。
Ｓｃｒ＝Ｓｉｊ（ｅｒｒ）×Ｉｄｉ×Ｉｑｊ＋Ｓ（ｉ＋１）（ｊ＋ｎ）（ｅｒｒ）×Ｉｄ（ｉ＋１）×Ｉｑ（ｊ＋ｎ）＋・・・＋Ｓ（Ｎ）  ・・・（２）
なお、「Ｉｄｉ」は、既知化合物のＭＳｎマススペクトル中のｉ番目のピークの強度スコアであり、「Ｉｄ（ｉ＋１）」は、既知化合物のＭＳｎマススペクトル中の（ｉ＋１）番目のピークの強度スコアであり、「Ｉｑｊ」は、未知化合物のＭＳｎマススペクトル中のｊ番目のピークの強度スコアであり、「Ｉｑ（ｊ＋ｎ）」は、未知化合物のＭＳｎマススペクトル中の（ｊ＋ｎ）番目のピークの強度スコアである。
これにより、各ピークについて、既知化合物のＭＳｎマススペクトル中のピークの強度スコアＩｄと、未知化合物のＭＳｎマススペクトル中のピークの強度スコアＩｑと、位置スコアＳ（ｅｒｒ）とを積算し、全ピークについての積算値の総和を算出し、その算出値にピークスコアＳ（Ｎ）を加算することで、未知化合物とＸ番目の既知化合物との類似性を示すスコアＳｃｒが算出される。このようにして、未知化合物のＭＳｎマススペクトルと、様々な既知化合物のＭＳｎマススペクトルとの類似性を示すスコアＳｃｒが算出されていく。
【００４３】
ここで、液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１を用いた未知化合物を同定するデータ処理方法について説明する。図３は、データ処理方法の一例について説明するためのフローチャートである。
まず、ステップＳ１０１の処理において、未知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・・）を作成する。
【００４４】
次に、ステップＳ１０２の処理において、未知化合物のＭＳｎマススペクトルについて、各ピークの質量電荷比（ｍ／ｚ）をそれぞれ求める。
次に、ステップＳ１０３の処理において、未知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・・）中の各ピークを、イオン強度Ｉｑの大きさによって３段階の強度スコアにそれぞれ分類する（強度スコア算出ステップ）。
次に、ステップＳ１０４の処理において、未知化合物のＭＳｎマススペクトル中のニュートラルロスのピークを求める（ニュートラルロス算出ステップ）。
次に、ステップＳ１０５の処理において、既知化合物の順番を示す既知化合物順番パラメータＸ＝１とする。
【００４５】
次に、ステップＳ１０６の処理において、スコアＳｃｒ＝０とする。
次に、ステップＳ１０７の処理において、既知化合物のＭＳｎマススペクトル中のピークの順番を示すピーク順番パラメータｉ＝１とする。
【００４６】
次に、ステップＳ１０８の処理において、既知化合物のＭＳｎマススペクトル中のｉ番目のピークの質量電荷比（ｍ／ｚ）に最も近い質量電荷比（ｍ／ｚ）のピークを見つけ出す。そして、既知化合物のＭＳｎマススペクトル中のｉ番目のピークの質量電荷比（ｍ／ｚ）と、未知化合物のＭＳｎマススペクトル中のｊ番目のピークの質量電荷比（ｍ／ｚ）との差ｅｒｒを算出し、差ｅｒｒを図２に示す分布関数に代入することで、位置スコアＳｉｊ（ｅｒｒ）を求める（位置スコア算出ステップ）。
次に、ステップＳ１０９の処理において、Ｓｃｒ＝Ｓｃｒ＋Ｓｉｊ（ｅｒｒ）×Ｉｄｉ×Ｉｑｊとする。
【００４７】
次に、ステップＳ１１０の処理において、ｉ＝ｉｍａｘ（既知化合物のＭＳｎマススペクトル中で最後の順番のピーク）であるか否かを判定する。ｉ＝ｉｍａｘでないと判定したときには、ステップＳ１１１の処理において、ｉ＝ｉ＋１として、ステップＳ１０８の処理に戻る。
一方、ｉ＝ｉｍａｘであると判定したときには、ステップＳ１１２の処理において、式（１）に基づいてピークスコアＳ（Ｎ）を求める（ピークスコア算出ステップ）。
次に、ステップＳ１１３の処理において、Ｓｃｒ＝Ｓｃｒ＋Ｓ（Ｎ）とする（スコア算出ステップ）。
【００４８】
次に、ステップＳ１１４の処理において、Ｘ＝Ｘｍａｘ（最後の順番の既知化合物）であるか否かを判定する。Ｘ＝Ｘｍａｘでないと判定したときには、ステップＳ１１５の処理において、Ｘ＝Ｘ＋１として、ステップＳ１０６の処理に戻る。
一方、Ｘ＝Ｘｍａｘであると判定したときには、ステップＳ１１６の処理において、スコアＳｃｒが高い順番に既知化合物を配列する。
そして、ステップＳ１１６の処理が終了したときには、本フローチャートを終了させる。
【００４９】
以上のように、液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１によれば、未知化合物のＭＳｎマススペクトル中のピークのイオン強度Ｉｑやピークの位置（ｍ／ｚ）やピークの有無等が変化しても、未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を正しく評価することができる。
【００５０】
＜他の実施形態＞
（１）上述した液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１では、強度スコア算出部３１ｂは、未知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・）中の各ピークと既知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・・）中の各ピークとを、イオン強度Ｉｑ，Ｉｄの大きさによって３段階の強度スコアにそれぞれ分類する構成を示したが、イオン強度Ｉｑ，Ｉｄのｌｏｇの変換値の大きさによって２段階や４段階の強度スコアにそれぞれ分類する構成としてもよい。さらに、強度スコア算出部は、未知化合物のＭＳｎマススペクトルにおいて、イオン強度Ｉｑの大きさが閾値以下であるピークを除去する構成としてもよい。
【００５１】
（２）上述した液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１では、強度スコア算出部３１ｂは、未知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・・）中の各ピークと既知化合物のＭＳｎマススペクトル（ｎ＝２，３，４，・・・）中の各ピークとを、イオン強度Ｉｑ，Ｉｄの大きさによって３段階の強度スコアにそれぞれ分類する構成を示したが、イオン強度Ｉｑ，Ｉｄの大きさとピークの質量電荷比の大きさとによって数段階の強度スコアにそれぞれ分類する構成としてもよい。例えば、イオン強度Ｉｑ，ＩｄがＩ１以上であるピークを上位の強度スコア「４」に分類し、イオン強度Ｉｑ，ＩｄがＩ２以上Ｉ１未満であるピークを中位の強度スコア「３」に分類し、イオン強度Ｉｑ，ＩｄがＩ２未満であるピークを下位の強度スコア「２」に分類するとともに、質量電荷比がｍ１／ｚ１以上であるピークに強度スコア「２」を加算し、質量電荷比がｍ１／ｚ１未満であるピークに強度スコア「１」を加算する。
【００５２】
（３）上述した液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１では、位置スコア算出部は、未知化合物のＭＳｎマススペクトル中のピークの質量電荷比を、アダクトイオンの質量を用いて補正する構成としてもよい。
【００５３】
（４）上述した液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１では、式（１）に基づいてピークスコアＳ（Ｎ）を求める構成を示したが、一致したペア数にしたがって増加するピークスコアＳ（Ｎ）を求める構成としてもよい。例えば、既知化合物のＭＳｎマススペクトル中のピーク数の８割以上が一致した場合にはピークスコア「１００」とし、既知化合物のＭＳｎマススペクトル中のピーク数の５割以上８割未満が一致した場合にはピークスコア「５０」とし、既知化合物のＭＳｎマススペクトル中のピーク数の２割が一致した場合にはピークスコア「２０」とする。
【００５４】
（５）上述した液体クロマトグラフ質量分析装置（ＬＣ／ＭＳ）１では、スコア算出部は、未知化合物のＭＳｎ＋ｍマススペクトルと、多数の既知化合物のＭＳｎ＋ｍマススペクトルとの類似性を示すＭＳｎ＋ｍマススペクトルのスコアをそれぞれ求め、ＭＳｎマススペクトルのスコアに、ＭＳｎ＋ｍマススペクトルのスコアを加算する構成としてもよい。
【産業上の利用可能性】
【００５５】
本発明は、質量分析装置によって得られた未知化合物のＭＳｎマススペクトルと、多数の既知化合物のＭＳｎマススペクトルとの類似性を示すスコアをそれぞれ求め、そのスコアに基づいて未知化合物を同定するデータ処理方法等に利用することができる。
【符号の説明】
【００５６】
１：  液体クロマトグラフ質量分析装置（データ処理装置）
３１ａ：  取得部
３１ｂ：  強度スコア算出部
３１ｃ：  位置スコア算出部
３１ｄ：  スコア算出部
３２ｃ：  データベース記憶領域
【図１】
0006004080.tif 000002
【図２】
0006004080.tif 000003
【図３ａ】
0006004080.tif 000004
【図３ｂ】
0006004080.tif 000005