(57)【要約】本発明は、質量分析法を用いた生体高分子同定技術において、測定開始前の校正作業やサンプルへの内部標準の添加を不要とした高精度の自動同定を可能とすることを目的とする。本発明の生体高分子自動同定方法においては、質量分析によって得られる実測質量値Ｘを所定のデータベースと照合させることによって候補分子を検索し、類似順位の高い任意数の候補分子を選び出し、候補分子を内部標準として用いて実測質量値Ｘを校正し、候補分子の校正質量値Ｘｃと理論質量値Ｍの相対誤差Ｅｃを算出し、相対誤差の標準偏差ＳＥｃを求め、この標準偏差ＳＥｃからデータベース検索における許容誤差Ｔｃを求め、この許容誤差Ｔｃに基づいて再度データベース検索を行う。
2004023132.tif 000002
【特許請求の範囲】
【請求項１】
試料中の生体高分子の質量を質量分析方法に基づいて測定する質量測定手順と、
前記質量測定手順によって得られる実測質量値を所定のデータベースと照合させることによって候補分子を検索するデータベース検索手順と、
類似順位スコアの高い任意数の候補分子を選び出す候補分子選出手順と、
候補分子を内部標準として用いて実測質量値を校正する質量値校正手順と、
前記手順により得られた候補分子の校正質量値と理論質量値の相対誤差を算出し、該相対誤差の標準偏差を求める手順と、
該標準偏差から前記データベース検索手順の許容誤差を求める手順と、
前記許容誤差に基づき再度前記データベース検索手順と、を行うことを特徴とする生体高分子自動同定方法。
【請求項２】
前記質量値校正手順は、前記候補分子選出手順で選び出された候補分子の実測質量値と理論質量値の相対誤差を算出し、
理論質量値と相対誤差のプロットに対する最小二乗直線を作成して実測質量値の系統誤差を見積もる手順と、
この系統誤差を全実測値から差し引くことにより実測質量値を校正する手順と、からなることを特徴とする請求の範囲第１項に記載の生体高分子自動同定方法。
【請求項３】
コンピュータシステムを利用することにより、請求の範囲第１項又は第２項に記載の生体高分子自動同定方法を構成する各手順を実行できるプログラム情報が格納されたことを特徴とする情報記録媒体。
【発明の詳細な説明】
【技術分野】    本発明は、質量分析方法を用いた生体高分子同定技術に関する。より詳しくは、質量分析方法によって得られる質量データの精度向上を目的とする生体高分子自動同定方法に関する。
【背景技術】    質量分析方法は、試料分子をイオン化した後、これを質量／電荷の比（ｍ／ｚ）に従って分離し検出を行う機器分析法であって、得られた質量スペクトルから定性を、イオン量から定量を行うことができる。
この分子質量の測定に用いる質量分析計（以下、「ＭＳ」（ｍａｓｓ  ｓｐｅｃｔｒｏｍｅｔｅｒ）と称する。）は、大別すると、試料のイオン化を行うための「イオン化部（イオン源）」と、イオンを質量／電荷の比であるｍ／ｚ（ｍ：質量、ｚ：電荷数）に従って分離するための「アナライザー」と、分離されたイオンの「検出部（検出器）」及び「データ解析部」と、から構成されている。
前記質量分析計を用いて試料分子の質量分析に当たっては、測定開始前に質量分析計の校正（キャリブレーション）を行う必要がある。具体的には、温度変化や電圧の精度、電気回路ノイズ等の原因により、質量分析計の測定に誤差が発生する場合があるので、測定開始前には、クロマトグラフ等を質量分析計から一旦取り外した上で、質量分析計に所定の質量校正用標準物質を導入して実測質量値を得、この実測質量と既知の理論質量値と比較することによって、質量値に系統誤差が発生しないように予め装置を調整する校正作業（外部標準法によるキャリブレーション作業）を行う必要がある。
更に、高精度の質量値を得るためには前記外部標準法によるキャリブレーション作業に加え、既知物質を試料に混合して質量測定し、その質量値にもとづいて実測質量値を調整する校正作業（内部標準法によるキャリブレーション作業）を行う必要がある。
そして、一般に、この質量分析計（タンデム質量分析計を含む。以下同様。）を用いて行うペプチドやタンパク質等の生体高分子同定方法においては、質量分析によって得られる未知試料分子の実測質量値を、１０万種程度の分子の一次構造又は配列が予め格納されているデータベース（ライブラリー）と照合させて検索（サーチ）し、構造から算出した予想リファレンス（標準）スペクトルの中から測定対象の未知試料分子のスペクトルと類似したものを順位（スコア）付けして選び出していく手順、即ちデータベース検索（又はライブラリーサーチ）を行って、候補分子をリストアップして絞り込み、最終的に未知試料分子の同定を行う。
しかしながら、上記した質量分析計の校正作業（キャリブレーション作業）は、非常に作業が面倒であって、調整時間もかかることから、従来の質量測定作業における作業効率を低下させる主原因であった。即ち、従来は、質量分析計の連続運転（校正作業なしの運転）による効率の良い測定作業を実施することができなかった。また、複数台の質量分析計を用いた測定系においては、各装置について外部標準による校正作業（キャリブレーション作業）を行ったとしても、各装置の精度、信頼性を一元化することは極めて困難であるという問題があった。
外部標準キャリブレーションの場合、従来の上記データベース検索の手順では、外部環境の影響によって発生する質量分析計自体の誤測定による影響を測定データから排除することはできなかった。特に測定環境の微妙な温度変化（０．２℃ぐらいの変化）で生じる測定誤差も時として無視できないものとなっていた。
また、従来の内部標準キャリブレーションによって複雑な生体高分子混合物を測定する場合は、内部標準物質と試料由来のイオンシグナルが重なってしまい、そのイオンを分析できないため、内部標準として試料に入れる物質の種類や濃度の選択はとても難しかった。高い質量精度を広い質量範囲で実現するためには、何点もの内部標準物質を導入する必要があった。
更には、従来は、同定の信頼性が低かったので、その結果を一つ一つ人間が確認しなければならかった。ところが、近年の質量分析計の発達により、より複雑な生体高分子混合物の直接分析が可能になってきたため、データが大量化し、一つ一つのデータを人間が目で確認することが困難になってきたことから、複雑な生体高分子混合物を分析対象とする信頼性の高い自動同定手法の開発が要請されていた。
【発明の開示】
そこで、本発明は、測定開始前の質量分析計の校正作業、あるいはサンプルに予め内部標準を添加することを不要とするとともに、データ処理のみに基づいた高精度で信頼性の高い生体高分子自動同定方法を提供することを目的とする。
上記技術課題を解決するために、本発明では、以下の（１）〜（７）の手順を少なくとも備える生体高分子自動同定方法を提供する。
（１）試料中の生体高分子の質量を質量分析方法に基づいて測定する質量測定手順。（２）前記質量測定手順によって得られる実測質量値を所定のデータベースと照合させることによって候補分子を検索するデータベース検索手順。（３）類似順位スコアの高い任意数の候補分子を選び出す候補分子選出手順。（４）候補分子を内部標準として用いて実測質量値を校正する質量値校正手順。（５）前記手順により得られた候補分子の校正質量値と理論質量値の相対誤差を算出し、該相対誤差の標準偏差を求める手順。（６）該標準偏差から前記データベース検索手順の許容誤差を求める手順。（７）前記許容誤差に基づき再度前記データベース検索手順。なお、上記「データベース」は、分子構造あるいは配列データベースを意味する。
ここで、上記（４）の質量値校正手順は、候補分子選出手順で選び出された候補分子の実測質量値と理論質量値の相対誤差を算出し、理論質量値と相対誤差のプロットに対する最小二乗直線（「ｙ＝ａ×Ｍ＋ｂ」の式で示される直線。Ｍは理論質量値。）を作成して実測質量値の系統誤差を見積もる手順と、この系統誤差を全実測値から差し引くことで、実測質量値を校正する手順を採用することができる。
例えば、飛行時間型質量分析計の場合において、候補分子の系統誤差を上記最小二乗直線から求める。この系統誤差を全実測値から差し引く。具体的には、（Ｘｃ−Ｍ）／Ｍ＝（Ｘ−Ｍ）／Ｍ−（ａＭ＋ｂ）［Ｘは実測質量値、Ｘｃは校正質量値、Ｍは理論質量値］からなる式を変形し、式：Ｘｃ＝Ｘ−Ｍ（ａＭ＋ｂ）を得る。
ここで、理論質量値Ｍは、候補分子については与えられているが、全ての実測値について与えられている訳ではない。このため、全実測値を校正するためには、上記式のＭ（ａＭ＋ｂ）の項を実測値で近似させる必要がある。ａ，ｂの値は、一般にＸ、Ｘｃと比較して非常に小さいため、Ｍ（ａＭ＋ｂ）≒Ｘｃ（ａＸ＋ｂ）とできる。これを上記式に代入して、Ｘｃ＝Ｘ−Ｘｃ（ａＸ＋ｂ）を得る。これを変形し、Ｘｃ＝Ｘ／（１＋（ａＸ＋ｂ））なる式を得て、この式を用いて、全ての実測値を質量校正する。
上記した本発明に係る生体高分子自動同定方法によれば、複雑な生体高分子混合物を対象として、データ処理のみにより、非常に高精度な質量値を得ることができる。得られる質量値の精度が高いと、より一義的に生体高分子を特定、同定することが可能となる。即ち、本発明は、複雑な生体高分子混合物を分析対象とする信頼性の高い自動同定手法を提供できる。
次に本発明では、コンピュータシステムを利用することにより前記生体高分子自動同定方法を構成する各手順を実行できるプログラム情報が格納されているＣＤ−ＲＯＭその他の情報記録媒体を提供する。
上記した手段によれば、測定開始前の質量分析計の校正作業、あるいはサンプルに予め内部標準を添加することを不要とすることができる。また、データ処理のみに基づいた高精度で信頼性の高い生体高分子自動同定方法を実施することができる。
【図面の簡単な説明】    第１図は、実施例１において同定した質量値（ｍ／ｚ）と誤差の関係を示す図である。
第２図は、実施例２において質量校正を行う前の同定結果を示す図である。
第３図は、実施例２において質量校正を行った後の同定結果を示す図である。
第４図は、実施例２において同定した質量値（ｍ／ｚ）と誤差の関係を示す図である。
【発明を実施するための最良の形態】    本発明に係る生体高分子自動同定方法の好適な一実施形態について説明する。なお、本発明は、以下の実施形態に限定されることはない。
まず、試料中の未知生体高分子の質量を、目的に応じた慣用の質量分析方法に基づいて測定し、実測質量値Ｘを得る。質量分析方法は、例えば、タンデム質量計を用いることができる。タンデム質量分析計は、アナライザーをタンデムに複数台結合した構成を備える質量分析計であって、詳しくは、最初のアナライザーで混合物中の特定のイオン（親イオン）を選択し、次のアナライザーで選択したイオンと不活性気体との衝突解離を行い、最後のアナライザーで解離した内部構造情報を示すイオン（生成イオン）を質量分析する構成を備える。
前記質量測定手順によって得られた実測質量値Ｘを、慣用のデータベース検索エンジンが読み込める形式（２値ファイル。質量値と強度。）に変換した上で、その質量値既知の分子が多数記録されたデータベースと照合させて、前記未知生体高分子に該当する可能性のある候補分子の検索を行う。
なお、上記する実測質量値Ｘの形式変換は、質量分析計メーカーから一般に提供されている慣用のＭａｓｓｌｙｎｘ（Ｍｉｃｒｏｍａｓｓ社）等のソフトウェアを適宜用いることによって行うことができ、データベース検索は、市販のＭａｓｃｏｔ（Ｍａｔｒｉｘ  Ｓｃｉｅｎｃｅ社）等のデータベースソフトウェアを用いて好適に行うことができる。
前記データベース検索手順の結果から、類似順位スコアの高い任意数の候補分子（のセット）を選び出す。なお、セットの大きさｎは、統計的処理が出来る程度の任意数である。
続いて、前記した候補分子選出手順によって選び出されてきた各候補分子の実測質量値Ｘと理論質量値Ｍの相対誤差Ｅを、次式（１）に従って算出する。
Ｅ＝（Ｘ−Ｍ）／Ｍ・・・・・（１）
続いて、前記手順によって得られた相対誤差Ｅの平均値ｍＥを次式（２）に基づいて算出する。
ｍＥ＝Σ（Ｅ）／ｎ・・・・・（２）
また、前記相対誤差Ｅの標準偏差ｓＥを次式（３）に基づいて算出する。この標準偏差により、候補分子を内部標準として用いることが妥当かどうかを判定する。なお、ｓＥ＜ｍＥであれば、校正は有効である。
ｓＥ＝｛Σ（Ｅ−ｍＥ）２／（ｎ−１）｝（１／２）・・・・・（３）
次に系統誤差の大きさを見積もり、これを実測質量値Ｘから差し引くことにより、校正質量値Ｘｃを得る。例えば、飛行時間型質量分析計の場合において、候補分子の相対系統誤差は、以下の手順で理論質量値と相対誤差のプロットに対する「最小二乗直線ｙ＝ａｘ＋ｂ」から求めることができる。候補分子の校正後の相対誤差Ｅｃ＝（Ｘｃ−Ｍ）／Ｍとすると、Ｅｃ＝Ｅ−（ａＭ＋ｂ）。したがって、
（Ｘｃ−Ｍ）／Ｍ＝（Ｘ−Ｍ）／Ｍ−（ａＭ＋ｂ）・・・・・（４）
［Ｘは実測質量値、Ｘｃは校正質量値、Ｍは理論質量値］
具体的には、上記（４）式を変形して、次式（５）を得る。
Ｘｃ＝Ｘ−Ｍ（ａＭ＋ｂ）・・・・・（５）
ここで、理論質量値は、候補分子については与えられているが、全ての実測値について与えられている訳ではない。このため、全実測値を校正するためには、上記式（５）の「Ｍ（ａＭ＋ｂ）」の項を実測値で近似させる必要がある。ａ，ｂの値は、一般にＸ、Ｘｃと比較して非常に小さいので、Ｍ（ａＭ＋ｂ）≒Ｘｃ（ａＸ＋ｂ）とできるから、これを上記式（６）に代入して、次式（６）を得る。
Ｘｃ＝Ｘ−Ｘｃ（ａＸ＋ｂ）・・・・・（６）
この式（６）を変形式である次式（７）に基づいて、全ての実測値を質量校正する。
Ｘｃ＝Ｘ／（１＋（ａＸ＋ｂ））・・・・・（７）
なお、前記最小二乗直線における「ｂ」と「ａ」は、それぞれ次式（８）、（９）によって求めることができる。
ｂ＝Σ｛（Ｍ−ｍＭ）×（Ｅ−ｍＥ）｝／Σ｛（Ｍ−ｍＭ）＾２｝・・（８）
ａ＝ｍＥ−ｂ×ｍＭ・・・・・（９）
さらに、ｍＭは、候補分子の理論質量値Ｍの平均値であって、次式（１０）によって求めることができる。
ｍＭ＝Σ（Ｍ）／ｎ・・・・・（１０）
質量校正後の質量値Ｘｃと理論質量値Ｍの相対誤差Ｅｃは、次式（１１）によって求めることができる。
Ｅｃ＝Ｅ−（ａＭ＋ｂ）・・・・・（１１）
続いて、候補分子について得られた相対誤差Ｅｃ＝（Ｘｃ−Ｍ）／Ｍの平均値ｍＥｃ及び標準偏差ＳＥｃを、それぞれ次式（１２）、（１３）に基づいて求める。
ｍＥｃ＝Σ（Ｅｃ）／ｎ・・・・・（１２）
ＳＥｃ＝｛Σ（Ｅ−ｍＥｃ）２／（ｎ−１）｝（１／２）・・・・・（１３）
求めた平均値ｍＥｃから校正を評価する。理想的にはｍＥｃ＝０となる。求めた標準偏差ＳＥｃからデータベース検索に用いる許容誤差Ｔｃを次式（１４）に基づいて算出することによって、一連の校正（キャリブレーション）手順を完了する。
Ｔｃ＝Ｋ×ＳＥｃ・・・・・（１４）
［Ｋ＝１．５〜３．０］
なお、Ｋは、質量値の信頼区間を指定するための経験的な定数を示す。このＫ値は、データベース検索に用いるソフトウェアの精度に応じて適宜決定できる。データベース検索ソフトウェアの同定性能が高いほど、９９．７％の信頼区間であるＫ＝３に近づけることが出来る。なお、Ｍａｓｃｏｔ（Ｍａｔｒｉｘ  Ｓｃｉｅｎｃｅ社）のデータベースソフトウェアの場合では、経験的にＫ＝１．５を採用できる。
得られた前記許容誤差Ｔｃ（Ｔｃ１）に基づいて、再度同様のデータベース検索を行う。必要に応じて、上記した一連の校正及びデータベース検索を複数回繰り返すことによって、許容誤差Ｔｃの範囲を徐々に狭めていき（Ｔ→Ｔｃ１→Ｔｃ２→・・・）、候補分子の選出精度を高める。なお、前記Ｔｃ１は一回目の校正作業によって得られた許容誤差を示し、Ｔｃ２は二回目の校正作業によって得られた許容誤差を示す。
これにより、候補分子同定の確度を高めていくことができる。即ち、未知試料分子の同定精度を向上させることができる。
以上説明した手順を所望のコンピュータプログラム情報に加工し、このプログラム情報をＣＤ−ＲＯＭ、フロッピー（登録商標）ディスクなどの各種情報記録媒体、コンピュータハードウエア、サーバ等に格納し、所望のコンピュータシステムやコンピュータネットワーク（情報通信技術）を介して、該プログラムを実行可能に工夫することができる。
【実施例】    飛行時間型質量分析計は、一定距離Ｌをイオンが飛行する時間を測定し、次の式（１５）で表される質量ｍと飛行時間Ｔの関係から、質量を測定する装置である。
Ｔ＝Ｌ・（２ｅＶ）＾（−１／２）・（ｍ／ｚ）＾（１／２）・・・・・（１５）
（ここで、ｅは電気素量、ｚは電荷数である。）
この装置の測定質量精度は、Ｌと加速電圧Ｖに依存する。Ｌは装置に固有の値であるが、主に温度により膨張収縮することにより変動し、Ｖは電源電圧のドリフトにより変動する。測定条件によっては、これらの変動により１００ｐｐｍ以上の系統的な質量誤差が生じることがある。しかし、一方、質量誤差同士のばらつき（質量分析計の性能を反映）は、系統的な誤差の平均値と比較して小さい。これを利用して、系統的な誤差だけを取り去ることができる。
以下に、実際に本発明の方法により同定精度が向上した例を示す。
【実施例１】    人血清アルブミンのトリプシン消化物１００ｆｍｏｌをＨＰＬＣ−ＭＳ／ＭＳで測定し、市販のデータベース検索ソフトＭａｓｃｏｔを用いてＭＳ／ＭＳ  ｉｏｎｓ  ｓｅａｒｃｈによりデータベース検索を行った。（検索パラメータ、Ｐｅｐｔｉｄｅ  Ｔｏｌｅｒａｎｃｅ  ２５０ｐｐｍ，ＭＳ／ＭＳ  ｔｏｌｅｒａｎｃｅ  ０．５Ｄａ）
検索結果の中から最もスコアが高い２０個のイオンについて同定された理論ｍ／ｚとの相対誤差Ｅ（（Ｘ−Ｍ）／Ｍ、単位ｐｐｍ）を求め、これを理論ｍ／ｚに対してプロットし、第１図に示した。第１図に見られるように、元の相対誤差Ｅ（第１図◆印）の平均値は約１７０ｐｐｍであるが、Ｅのばらつきは１５０−１７５ｐｐｍの範囲に収まり、Ｅ自体の値と比較すると小さかった。
このイオン群に対する最小自乗直線を求め、これを各イオンの誤差から差し引くことで質量を校正した。校正後の相対誤差Ｅｃ（第１図■印）も同様にプロットし、第１図に示した。このＥｃのばらつき（標準偏差で代表）から求めたデータベース検索パラメータは、Ｐｅｐｔｉｄｅ  Ｔｏｌｅｒａｎｃｅ  １８ｐｐｍ，ＭＳ／ＭＳ  ｔｏｌｅｒａｎｃｅ  ０．０８０Ｄａであった。この質量校正により、検索時の許容誤差が２５０→１８ｐｐｍ、０．５→０．０８０Ｄａとそれぞれ約１４倍、６倍狭めることができ、同定の信頼性が向上した。
【実施例２】    次に、本発明の質量校正法により、実際に誤同定を訂正できることを以下に示す。
質量データを用いたデータベース検索により、誤同定しやすいことが知られているペプチドＳＲＬＤＱＥＬＫを定法により合成した。このペプチド１００ｆｍｏｌを上記の人血清アルブミンのトリプシン消化物１００ｆｍｏｌと混合し、同様に実験を行った。通常の検索条件（検索パラメータ、Ｐｅｐｔｉｄｅ  Ｔｏｌｅｒａｎｃｅ  ２５０ｐｐｍ，ＭＳ／ＭＳ  ｔｏｌｅｒａｎｃｅ  ０．５Ｄａ）では、合成ペプチドは第２図に示すように誤同定した。
次に、上記の通り質量校正をしたところ、第３図に示すように、正しいペプチドを同定することができた。
このペプチドのＭＳ／ＭＳスペクトル中の各イオンを同定したそれぞれのペプチド（ＥＫＬＴＱＥＬＫとＳＲＬＤＱＥＬＫ）の理論的な生成イオン（ｂ、ｙイオン系列）にアサインし、その系統誤差をｍ／ｚに対してプロットし、第４図に示した。ＳＲＬＤＱＥＬＫ（第４図◆印）ではすべてのイオンの相対誤差が狭い範囲に収まるのに対し、ＥＫＬＴＱＥＬＫ（第４図■印）では２つの異なった分布を示した。このように、データ処理により質量精度を向上することで、よく似た質量を持ち、ｃ末端部分の配列が同一のペプチドを区別し、正しく同定することが可能となった。
【産業上の利用可能性】    本発明によれば、測定開始前の質量分析計の校正作業、あるいはサンプルに予め内部標準を添加することを不要とすることができるので、質量分析計の連続運転（校正作業による中断のない運転）が可能となる。その結果、作業者は煩わしい装置調整作業から開放され、分子同定作業の効率を向上させることができる。
また、質量分析計自体の誤差の影響を排除し、データ処理のみに基づいた高精度で信頼性の高い生体高分子自動同定方法を実施することができ、そして、複数台の質量分析計を用いた測定系では、各質量分析計から得られるデータ精度の一元化を達成できるので、未知試料分子の誤同定を確実に防止することができる。
【図１】
2004023132.tif 000003
【図２】
2004023132.tif 000004
【図３】
2004023132.tif 000005
【図４】
2004023132.tif 000006
【国際調査報告】
2004023132.tif 000007
2004023132.tif 000008