(57)【要約】クロマトグラフィ／分光測定法、特に２次元多変量統計解析を用いる液体クロマトグラフィ（ＬＣ）／質量分析（ＭＳ）から得た３次元データを解析する方法がここに開示されている。ＬＣ部はサンプル中の分析物の分離を可能にする。このような方法から得られた情報は一般に保持時間（Ｒ）に依存する。この分析物はシステムを通って進行するとともに、ＬＣ／ＭＳシステムのＭＳ領域に入る。この分析物はそこでイオン化され、その後質量検出器がこれらイオン化種を検出する。この方法で入手した情報は一般に、対応するｍ／ｚ値の強度として報告される。したがって、ＬＣ／ＭＳシステムは少なくとも３個の情報を供給する。主成分分析（ＰＣＡ）は、異なったサンプル間でこのタイプのデータの多変量解析を行う頑強たる方法である。しかしながら、一般に、ＰＣＡ分析は２次元データだけを使用して行われる。したがって、ＬＣ／ＭＳ操作から得られるデータの少なくとも１つのパラメータが失われることが予期されるが、ここで開示した方法では、３つの数学的次元はすべて保存される。
2005528606.tif 000002
【特許請求の範囲】
【請求項１】
（ａ）サンプルについてデータを得るステップ（前記データは少なくとも３つのパラメーターによって特徴づけられ、前記パラメーターの少なくとも１つは時間と相関関係にある。）、
（ｂ）ステップ（ａ）からの前記データをビンニングして複数のビンニングされたデータ集合を形成するステップ（ビンサイズは時間パラメーターに基づいて選択され、前記データ集合におけるデータ点は残りのパラメーターの２つを用いて特徴づけられる。）、
（ｃ）前記ビンデータ集合のそれぞれを整列されたビンデータ集合に変換するステップ（すべてのビンデータ集合は同じ数のデータ点を有し、および前記整列されたデータ集合における各第１パラメーターについて少なくとも１つのビンデータ集合がゼロでない第２パラメーターを有するように、ゼロの第２パラメータを有する整列データ点が前記ビンデータ集合に加えられる。）
（ｄ）前記整列されたビンデータ集合を、１つの軸上の前記第１パラメータおよび第２の軸上のビン番号および交差部分に記録されたデータ点の前記第２パラメータにより揃えることにより２次元配列を形成するステップ、および
（ｅ）ステップ（ｄ）から得た前記データを２次元多変量統計解析にかけるステップ
を含む、２次元解析法を使用してサンプルから得た３次元データを解析する方法。
【請求項２】
前記各ビンが共通の要素をもたない請求項１に記載の方法。
【請求項３】
前記各ビンが部分的に共通するところがある請求項１に記載の方法。
【請求項４】
前記ビンサイズが前記時間パラメータの全範囲を包含するビンから前記時間パラメータの２つ以下の同定し得る最小の増分を包含するビンにまでに及ぶ請求項１に記載の方法。
【請求項５】
前記２次元多変量統計解析が、古典的最小二乗分析（ＣＳＬ）、主成分分析（ＰＣＡ）、部分最小二乗分析（ＰＬＳ）および部分最小二乗判別式分析（ＰＬＳ−ＤＡ）から成る集合から選択される請求項１に記載の方法。
【請求項６】
前記２次元多変量統計解析が主成分分析（ＰＣＡ）である請求項５に記載の方法。
【請求項７】
ステップｃの変換が、各データ点が関連していることを前記第１のパラメータの各値が示す場合は前記データ点を併合することをさらに含み、前記併合は前記第２のパラメータの各値を加えることによって達成される請求項１に記載の方法。
【請求項８】
前記３次元データが分光測定法と組み合わされたクロマトグラフィ法から得られる請求項１に記載の方法。
【請求項９】
前記クロマトグラフィ法が液体クロマトグラフィ法である請求項８に記載の方法。
【請求項１０】
前記分光測定法が、紫外線分光法、赤外線分光法、質量分析法および核磁気共鳴から成る集合から選ばれる請求項８に記載の方法。
【請求項１１】
前記分光測定法が質量分析法である請求項１０に記載の方法。
【請求項１２】
非時間ベースのパラメータが、前記サンプルの物理的／化学的性質に対応する物理的にアクセスできる変数および各ビンに含まれる各種の変数に対する前記変数の強度であり、前記時間パラメータがＲｔの継続時間である請求項８に記載の方法。
【請求項１３】
前記サンプルの物理的／化学的性質に対応する前記物理的にアクセスできる変数がｍ／ｚである請求項１２に記載の方法。
【請求項１４】
ステップｃの変換が、各データ点が同じイオンの同位元素であること示すｍ／ｚパラメータを有する前記データ点を併合することをさらに含み、前記併合が各強度値を加えることによる請求項１３に記載の方法。
【請求項１５】
前記３次元データが液体クロマトグラフィ／質量分析法から得られる請求項８に記載の方法。
【請求項１６】
（ａ）サンプルについてクロマトグラフィ／質量分析データを得る段階（前記データは少なくとも３つのパラメーターによって特徴づけられ、前記パラメーターの少なくとも１つはＲｔ、ｍ／ｚおよび強度からなる集合から選択される。）
（ｂ）ステップ（ａ）からのデータをビンニングして複数のビンニングされたデータ集合を形成するステップ（ビンサイズはＲｔに基づいて選択され、前記ビンデータ集合におけるデータ点は、各ビン内に含まれる各イオンについての１つまたはそれ以上のｍ／ｚおよび強度の値を用いて特徴づけられる。）、
（ｃ）前記ビンデータ集合のそれぞれを整列されたビンデータ集合に変換するステップ（すべてのビンデータ集合は同じ数のデータ点を有し、および前記整列されたデータ集合における各ｍ／ｚ値について少なくとも１つの整列されたビンデータ集合が前記ｍ／ｚ値に関連するゼロでない強度値を有するように、ゼロの強度パラメータを有する整列データ点が前記ビンデータ集合に加えられる。）
（ｄ）前記整列されたビンデータ集合を、１つの軸上の前記ｍ／ｚ、第２の軸上のビン番号および交差部分に記録された強度値により揃えることにより２次元配列を形成するステップ、および
（ｅ）ステップ（ｄ）から得た前記データをＰＣＡ分析にかけるステップ
を含む、主成分分析（ＰＣＡ）を使用して３次元データを解析する方法。
【請求項１７】
前記クロマトグラフィ／質量分析データが液体クロマトグラフィ／質量分析法から得られる請求項１６に記載の方法。
【請求項１８】
ステップｃの変換が、各データ点が同じイオンの同位元素であることをｍ／ｚパラメータが示す場合は前記データ点を併合することをさらに含み、前記併合は各強度値を加えることによる請求項１６に記載の方法。
【請求項１９】
前記２次元配列がスプレッドシートに形成される請求項１６に記載の方法。
【請求項２０】
前記各ビンが共通の要素をもたない請求項１６に記載の方法。
【請求項２１】
前記各ビンが部分的に共通するところがある請求項１６記載の方法。
【請求項２２】
前記ビンサイズが前記時間パラメータの全範囲を包含するビンから前記時間パラメータの２つ以下の同定し得る最小の増分を包含するビンにまでに及ぶ請求項１６に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、データビンニングのデータへの応用によってより少ない次元に適用できる解析法により多次元データを解析する方法に関する。特に、本発明は、データビンニングを用いるクロマトグラフィ／分光測定データの解析方法に関する。本発明は、主成分分析と組み合わされたデータビンニングを用いる液体クロマトグラフィ／質量分析（ＬＣ／ＭＳ）データ解析を行う方法に特に使用される。
【背景技術】
【０００２】
本出願は、参照によりその内容全体が本明細書に特に組み込まれている２００２年５月３１日に出願した米国仮出願番号６０／３８４，７１２（整理番号ＷＡＡ−３０６）の優先権を主張するものである。
【０００３】
メタボノミクス（Ｍｅｔａｂｏｎｏｍｉｃｓ）は急速に発展している科学的研究の分野である。それは生体代謝特性を研究するためのシステムズアプローチであり、特定の疾患、毒性、および遺伝子機能に関する情報を提供することができる。メタボノミクスにおいて、動物または生物全体に対する薬剤候補の効果は、化合物投与に続くある時間の間に起こる代謝の変化を研究することにより調べられる。これらの研究で生じた分析データは、主成分分析（以後「ＰＣＡ」と称す）などの多変量数学的手法によって解析される。この数学的手法は調べているサンプルの微妙な差異と全体の差異の両方を強調するために用いられる。
【０００４】
今日まで、この分野の研究の大部分が、選択した分析方法として陽子核磁気共鳴（ＮＭＲ）を利用してきた。ＮＭＲは非常に有効であるが、いくつかの欠点、すなわち、感度が劣る、分析に時間が掛かる、およびいくつかの化学物質類、たとえば、硫酸塩類の検出ができないなどの欠点を有している。ＮＭＲのさらなる欠点は、すべての信号が１つのスペクトルに含まれるので、１つの大きな化合物が簡単に低濃度分析物の微妙なしかし重要な変化を隠してしまうことである。さらに、生体異物に関連する化合物をＮＭＲスペクトルから除去する必要があるが、この除去が関心のある内生化合物から得られる信号も除去し、そのため後に続くＰＣＡ分析に使用するデータ集合を低減してしまう。
【０００５】
紫外線分光法、赤外線分光法、核磁気共鳴または質量分析法などの分光測定法と組み合わせられたガスおよび液体の両クロマトグラフィは、強力な手法に進化した。液体クロマトグラフィと結合されたエレクトロスプレー質量分析法（ＬＣ／ＭＳ）は、定量的および定性的両方のバイオ分析に対して選択される手法になっている。この手法は、確固たるものであり、感度が高く、かつ選択性があり、ｐｇ／ｍＬまでの範囲の感度が容易に達成される。短いカラムおよび急勾配の使用により、ＬＣ／ＭＳが良好なクロマトグラフィの分解能と感度を備えた１分当たり１サンプル領域内の分析サイクル回数を達成することが可能になる。それゆえにＬＣ／ＭＳのメタボノミクスへの応用は論理的な手順である。質量分析計が高感度の定量化と構造上の情報の両方を与える一方で、クロマトグラフィのステップは、分析物を分離することにより時間分解能を与え、重なり合う信号の問題に対処する。ＬＣ／ＭＳの使用は、クロマトグラフィのステップが任意のある時間における質量分析計に入る競争イオンの数を低減することから、イオン抑制を低減するので、質量分析法の導入以上に好ましい。ＮＭＲに対するＬＣ／ＭＳの別の利点はデータ解析において観測される。生体異物に関連する物質を除去する場合、１つまたは２つの特定質量におけるわずか一切れの時間だけがデータ集合から取り除かれることになる。したがって、残りのＬＣ／ＭＳデータは変更されないで残り、数理解析、たとえばＰＣＡなどに利用できる。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
主成分分析はＭＳから得られるデータを解析する非常に有効な数学的デバイスである。しかしながら、ＰＣＡは２次元の手法であり、それに対してＬＣ／ＭＳから得られるデータは３次元である。現在、従来のＰＣＡ分析において通常失われるクロマトグラフィ分離情報をＰＣＡが保持することを可能にする必要性が存在する。
【０００７】
一般に、結果についての情報が、多数のパラメータにより特徴づけられるデータの中に含まれる実験条件が多数ある。説明のツールがこれらのパラメータを取り扱えない場合は情報が失われる。ＰＣＡ法により分析されるＬＣ／ＭＳデータの解析の分解能を向上する方法は、他のデータ集合に対しても適用できる。
【課題を解決するための手段】
【０００８】
本発明は、より少ない次元を取り扱うツールを使用して多次元データを解析する方法に関する。特定の実施形態において、本発明は、主成分分析（ＰＣＡ）などの２次元多変量統計解析手法を用いて、液体クロマトグラフィ／質量分析（ＬＣ／ＭＳ）試験などのクロマトグラフィ／分光測定試験から得た３次元データを解析する方法に関する。一般には、ＰＣＡ分析は２次元データだけを使用して行われる。本発明においては、３次元データを使用して解析が行われる。
【０００９】
クロマトグラフィ操作から得られたデータは、いくつかのパラメータによって特徴づけられる。特に、クロマトグラフィは時間の関数としてサンプルを分離する。保持時間、すなわちサンプルの１成分が固相上に保持される時間の長さは特に情報を与えるものである。このサンプルの種々の成分は、条件により異なった時間に固相から溶離する。特性に依存するこの時間は、サンプルの特定の成分を区別する場合に有用である。
【００１０】
さらなるパラメータは、クロマトグラフィ分離の後に使用される検出方法に特有のものでもよい。たとえば、蛍光の検出において、関心のあるパラメータは光の波長および強度であり、一方、紫外線の検出おいては、関心のあるパラメータは吸収された光の量である。多くの場合、解析方法で取り扱うことができるパラメータの数が多ければ多いほど、その結果の識別力が高くなる。１つのパラメータを完全に無視しなければならない場合、そのパラメータだけが異なる各データ点を識別する機会が失われる。
【００１１】
ＬＣ／ＭＳ操作から得られるデータは３つのパラメータから成る。これらのパラメータの第一番目は、このシステムのＬＣ部に固有のものである。液体クロマトグラフィは分析物を時間の関数として分離する。したがって、ＬＣ部から得られるデータは時間ベース、すなわち保持時間（Ｒｔ）である。種々の分析物が、溶離液などに応じて異なった時間にＬＣシステム（カラム）の固相エレメントから溶離することになる。この識別分離は、各分析物に対してＲｔ特性を考慮している。この情報は、たとえば、特定の分析物の素性を決定する際に役立つことができる。
【００１２】
第２および第３のパラメータは、ＬＣ／ＭＳシステムのＭＳ部に固有のものである。ＬＣを経由してサンプルがいったん分離されると、個々の分析物はこのシステムのＭＳ部に入る。ＭＳ部の条件により、分析物はある一定容積内でイオン化され、その後分子質量検出器に向けて加速される。質量電荷比（ｍ／ｚ）および強度は質量検出に従って算出することができる。この情報は２つ以上のサンプル間の差異を求めるのに使用される。
【００１３】
一般に、ＰＣＡ分析は、このｍ／ｚおよび強度パラメータだけを使用してこの定量試験を行うために使用される。それゆえ、Ｒｔデータは全体的な分析には寄与しない。質量分析ではイオンピーク当たりおおよそ５から１０原子質量単位（ＡＭＵ）点がある豊富なデータを利用できるためこの犠牲は許容されてきた。本発明は、データビンニングを組み込むことにより、この過ちに対処するものである。データビンニングにより、ＰＣＡ分析においてＲｔ情報が失われず、そのため、３次元データ解析が容易になる。サンプル間の比較は現時点でＲｔ、ｍ／ｚおよび強度データを含む。この手法は、大量の複雑なデータの中に単純な疑問に対する答えを見いださなければならない複雑なサンプルの結果を調査するのに特に役に立つ。その結果、従来検出できなかったサンプル間の微妙な変化も現在では解明できる。
【００１４】
１つの実施形態において、個々のＬＣ／ＭＳクロマトグラムは所定の部分（またはビン）に分割される。これらのビンはＲｔパラメータを反映している。このビンのサイズは操作者が決定する。たとえば、ビンサイズは、クロマトグラムを通る時間が１分であることを反映して１分にすることができる。各ビン内で、そのビン時間の間に行われるすべてのスキャンについての強度情報は、対応するＲｔ時に存在するそれぞれの種（ｍ／ｚ）について合計される。ビン、またはもっと正確に言えばその中に含まれるデータは、その後、変換され、市販のスプレッドシート、たとえば、Ｅｘｃｅｌなどのデータ行列にエクスポートされる。したがって、元のグラフによるデータは表形式のデータに変換される。変換された各ビン内の情報はクロマトグラフィの時間、ｍ／ｚおよび強度である。２つ以上のビンをこの変換処理にかけることができる。たとえば、６０分のＬＣ手順の場合、１分につき１個のビンに相当する６０個のビンとすることができる。別法として、重複するやり方で各ビンを構成することができるので、たとえば、４０秒につき１回発生する９１分用ビンが使用されて６０分のＬＣ手順を分析する。変換がいったん完了すると、各データ集合（表内のデータ）はゼロを用いて整列される。この各データ集合の整列は単に隙間を埋める手順を指す。たとえば、ｍ／ｚ＝１００におけるビン３内に強度値１０２があり、ｍ／ｚ＝１００におけるビン４内に強度値が無い場合は、ビン４内のｍ／ｚ＝１００に対応する強度の位置にゼロが置かれる。言い換えると、ゼロはプレースホルダの役割を果たし、その結果適切なデータ整列が可能になる。したがって、各データ集合は、各ｍ／ｚにつき１つの強度値（ゼロを含む）を有する。さらに、同位体存在量は一緒に、たとえば、Ｃ１２およびその同位体Ｃ１３を合わせた値にすることができる。最後に、スプレッドシートのファイルが市販のＰＣＡシステム、たとえばＰＣＡ分析用ＭａｔＬａｂに移される。
【発明を実施するための最良の形態】
【００１５】
（詳細な説明）
本発明は、より少ない次元を取り扱うツールを使用して多次元データを解析する方法に関する。本発明は、その後に分光測定方法が続くクロマトグラフィ方法から得たデータの解析に特に応用できる。これらのケースにおけるデータ集合は時間成分、ならびに一般にサンプルの物理的または化学的性質に対する強度の値を有している。特定の実施形態において、本発明は、主成分分析（ＰＣＡ）を用いて、液体クロマトグラフィ（ＬＣ）／質量分析（ＭＳ）方法から得たデータを解析する方法に関する。一般には、ＰＣＡ分析は２次元データだけを使用して行われる。本発明においては、３次元データを使用して解析が行われる。ＬＣ／ＭＳおよびＰＣＡ分析の詳細説明では、一般的ケースならびにこの特定の実施形態の方法論について示す。
【００１６】
このシステムのクロマトグラフィ部はサンプルマトリックス中の分析物の分離を可能にする。このような方法から得られた情報は一般に保持時間（Ｒｔ）として報告される。この分析物がシステムを通って進行するとともに、ＬＣ／ＭＳシステムのＭＳ領域に入る。それらはそこでイオン化され、その後質量検出器が、ｍ／ｚおよび強度の情報をもたらすこれらイオン化種を検出する。したがって、ＬＣ／ＭＳシステムは少なくとも３個の情報を供給することができる。多くの多変量統計解析手法が使用できると同時に、主成分分析（ＰＣＡ）はデータを解析する確固たる方法であり、異なったサンプルを比較する場合、かなり有効である。しかしながら、ＰＣＡ分析は一般に２次元データ集合だけを使用して行われる。したがって、ＬＣ／ＭＳ操作から得られるデータの少なくとも１つのパラメータが失われることが予期される。本発明においては、ＰＣＡ分析に加えて３次元データ集合を使用して解析が行われる。
【００１７】
ガスクロマトグラフィ（ＧＣ）はガス状のサンプルに対して同じ様なデータを与えるが、サンプルマトリックス中の分析物を分離するのに一般に用いられる方法はクロマトグラフィ、特に、液体クロマトグラフィ（ＬＣ）である。高性能液体クロマトグラフィ（以下において「ＨＰＬＣ」と称す）により、あるサンプル内の分析物の迅速で効率的な分離および特性記述が可能になる。ＨＰＬＣシステムの構成部品には、システムを通る水相および／または有機相の移動を容易にする高圧ポンプが含まれる。この水相（すなわち移動相）は、ＨＰＬＣシステムを最初に平衡させるのに使用される溶剤を含む。この溶剤は、ＨＰＬＣシステム全体を通って移動する分析物に水性環境も与える。最後に、移動相はＨＰＬＣカラムから分析物を溶離する溶剤を含む。
【００１８】
ＨＰＬＣシステムの別の構成部品はクロマトグラフィ分離カラムである。このカラムは固相を含む。この固相は移動相と組み合わさって、サンプルマトリックスに含まれる分析物の時間差のある分離を行う。この固相は一般に、ある種の分析物と相互に作用する化学重合物から成る。たとえば、逆相ＨＰＬＣカラムは、疎水性の力によって分析物と相互に作用する固相の化学的性質（たとえば、シリカのビーズに付いている１８個の炭素から成る炭化水素鎖）を有している。
【００１９】
分析物がいったんカラムから溶離されると、検出器内に移動しそこを通りぬける。ＨＰＬＣシステムで用いることができるものとして多様な検出器システムがある。たとえば、紫外線（「ＵＶ」）域内の分析物を検出するＵＶ検出器、赤外線（ＩＲ）域内の分析物を検出するＩＲ検出器、電解電量法を使用して分析物を検出する電気化学検出器、磁場に対する反応に基づいて分析物を検出する核磁気共鳴検出器および質量電荷比に基づいて分析物を検出する質量分析検出器がある。
【００２０】
ＬＣ／ＭＳシステムにおいて、サンプルマトリックスは最初ＬＣ分離にかけられ、これによってサンプル内の分析物の分離を行う。多くの場合分析物のＲｔによって伝えられるこの情報は、分析物の特性を記述し、分析物を特定するのに使用することができる。しかしながら、特に異なるサンプル間で比較を行うために、分析物をＭＳ検出にかけることによってさらなる情報を得ることができる。一般に質量分析計は、ある容積内の分子をイオン化し、その後イオン化された分子を分子質量検出器に向けて加速する。分子のイオン化は、電子イオン化、化学イオン化、エレクトロスプレーイオン化、または光イオン化により生じ得る。このイオン化プロセスは、大気圧を含む種々の圧力条件下で発生することができる。ＭＳから得られる情報には特定のｍ／ｚ値に対する応答の強度が含まれる。２つ以上のサンプルを比較する場合に重要なのはこの強度（ならびにｍ／ｚ値）情報である。
【００２１】
ＭＳ情報は比較分析に利用することができる。たとえば、人間の細胞系について実験しているとする。この研究の目標は、もしあるとすれば、化合物Ｘは特定の代謝産物に対してどのような作用を及ぼすか判定することである。話を続けて、操作者は、対照サンプルすなわち化合物Ｘがなんら加えられていない細胞培養からの細胞抽出であるサンプル１と、化合物Ｘが加えらた細胞培養系から取った細胞抽出であるサンプル２の２つのサンプルを得たとする。これら２つのサンプルはＬＣ／ＭＳにかけることができる。さらに、関心のある代謝産物が、逆相カラムを用いるクロマトグラフィ単離に基づいて分析できるとする。ＬＣを通っていったん分離されると、サンプル１とサンプル２のサンプルは別々にＭＳにかけることができる。２つのサンプルから発生したＭＳ信号は、その後、特に関心のある代謝産物について比較を行って、化合物Ｘの添加によってその強度が増加するか、減少するか、または影響を受けないか決定することができる。
【００２２】
図１は、ＬＣ／ＭＳ実験に適用した場合の本発明の様々なステップを記した流れ図を示す。図１ａは、仮想サンプルに関してＬＣ／ＭＳ手順から得られたデータをグラフ形式で示す。図示されるように、このグラフは３つのパラメータを規定している。Ｘ軸はｍ／ｚパラメータを規定しており、たとえば、Ｘ軸は１００から８００ｍ／ｚまで０．１ｍ／ｚ刻みで延びることができる。図１ａのＹ軸は強度または応答パラメータを表す。強度は、あるイオンのパーセンテージ、たとえば基底のピークに設定することができる。最後に、時間パラメータを表すＺ軸がある。図１ａに示されているのは、３つのパラメータすべてではないにしても少なくとも１つが異なるイオン種である。それでも依然として各イオン種は図１ａのパラメータによって一意的に定められる。
【００２３】
現存のＰＣＡ分析では、図１ａに示されたデータを分析前に２次元データ集合に変換することが要求される。たとえば、従来のＰＣＡ分析はｍ／ｚおよび強度の各パラメータを使用して実行することができる。しかしながら、このことにより時間パラメータが削除される。より完璧な分析は、時間次元を含め入手できるすべての関連情報を含む。図１ｂは、３つのパラメータすべてを保存し分析にかける方法を表している。データビンニングは、多数の測定点に相当する１つのデータ集合を、測定点を集めた多数のデータ集合に分離するプロセスである。ビンは１つのパラメータの範囲、すなわちビンパラメータに対応する。各ビン内のビンパラメータの変動は、ビンの値におけるビン内のすべてのデータ点を累積することなどにより除かれる。データ集合全体を包含するビンはビンパラメータを無視するのと同じ効果をもつことに注目すること。共通の要素をもたないビンを用いることは最も一般的であるが、ある分析に対しては一部共通するところがあるビンを使用してもよい。図１ａから図１ｂへの移行において、所定のビンサイズが選択されるとデータビンニングのプロセスが生じる。ビンサイズは時間次元の関数である。たとえば、ビンサイズは２分にすることができる。多数のビンが図１ｂに示されている。各ビンは、スキャン全体で見いだされ、そのビンの範囲内に入るイオン種についてのｍ／ｚおよび強度の両方の情報を含んでいる。ビン（または領域）サイズは、経験および試験されるサンプルに基づいて従業者によって選択される。選択方法は当業者によく知られている。たとえば、ソフトウェアのＭａｓｓＬｙｎｘを用いることができる。このソフトウェアを使用して、ビンが以下のように集められる。（１）全体のイオンクロマトグラムを開く、（２）マウス右クリックでビン領域を選択しビンエリア上にドラッグすれば、すべてを含めた質量スペクトル（ビン）が現れる、（３）すべてを含めた質量スペクトル（ビン）のウィンドウ内に「リストスペクトル（Ｌｉｓｔ  Ｓｐｅｃｔｒｕｍ）」のオプションが現れ、始動する、（４）リストスペクトル（Ｌｉｓｔ  Ｓｐｅｃｔｒｕｍ）が選ばれると、ウィンドウが対応する存在量を備えたｍ／ｚを２行形式で表示する、（５）２つのカラム内に含まれるデータを複写する、および市販の（６）Ｅｘｃｅｌのようなスプレッドシートに貼り付ける。
【００２４】
データビンニングプロセスの後に、変換ステップがある。（図１ｂから図１ｃ参照）。この変換は、得られた分離リストをグラフによるデータ形式から行列データ表現に変換する。従業者は、単に「切り取りおよび貼り付け」手順を用い、２行のデータ集合情報をＥｘｃｅｌのような市販スプレッドシートに与えることでこの変換を始めることができる。変換のこの部分がいったん完了すると、データ集合はゼロを用いて整列させる必要がある。各ビンは、分析されるすべてのイオンの値を有する必要がある。分析の忠実度を確保するために、そのビンについてイオン信号が観測されないところではゼロが強度値として挿入され、各データ集合が各ｍ／ｚパラメータに対して１つの強度値を有する結果になる。その後、ｍ／ｚパラメータの各行は、行列の１つの軸として作用する１つの行に融合させることができる。別法として、この変換を自動化するためにソフトウェアプログラムを使用することができる。表の形式で表したものが図１ｃに示されている。図１ｃにはｍ／ｚの行とビンの列がある。ｍ／ｚの行は分析されるｍ／ｚイオンを含んでいる。ビンの列は分析される望ましいビンを含んでいる。行列の交差部分は、適切なｍ／ｚおよびビンパラメータに対応する強度値を含んでいる。さらに、データを解析にかける前に、試験によって引き起こされたｍ／ｚおよび強度における知られた変換に基づいて、操作者がデータ操作を行ってもよい。たとえば、同位体存在量は一緒にすることができ、その結果、たとえばＣ１２およびその同位体Ｃ１３に対する各値が一緒にされる。
【００２５】
この行列またはスプレッドシートは次に、各種サンプル間の比較を行うことができる市販のＰＣＡプログラムに移され、解析される。このような市販のＰＣＡプログラムの１つが、マサチューセッツ州ネーティックのＭａｔｈＷｏｒｋｓ社が販売しているＭａｔＬａｂである。別のＰＣＡプログラムパッケージは、ワシントン州ウディンビル（Ｗｏｏｄｉｎｖｉｌｌｅ）にあるＩｎｆｏＭｅｔｒｉｘ社から入手でき、Ｐｉｒｏｕｅｔｔｅと呼ばれている。
【００２６】
複数のスキャンの結果として生じたデータをビンニングを通じて合同することにより質量分析法特有の高分解能のいくつかをぼかすことによって、２次元解析は得られた３次元データを取り扱うことができる。これにより、試料間の変化が予期されるが種々のスキャンにおける変化の位置は分からない大量のデータを解析することが可能になる。変動の位置が分かった後で、後半の分析は、生データ本来の高分解能を使用してさらに差異の特徴を調べることができる。このビンニング手法は、データ集合を多くのタイプの多変量統計解析に適合させるように使用することができる。
【００２７】
本発明の特徴およびその他の詳細について以下の実施例でより詳細に説明し指摘する。本発明の特定の実施形態を図を用いて示すが、本発明を限定するものとして示すものではないことが理解されるだろう。本発明の主な特徴は発明の精神と範囲から逸脱することなく様々な実施形態で使用することができる。
【実施例】
【００２８】
ラットの尿サンプル
合計２０匹のラットのサンプルがこの実験に使用された。これらのラットは４つのグループに分けられた。グループＩは投与マトリックスだけを受け入れ、グループＩＩ、ＩＩＩおよびＩＶは化合物Ａ、ＢおよびＣをそれぞれ経口投与された。ラットの尿のサンプルは２つの時点、投与後０〜８および８〜２４時間に集められた。ラットの尿は分析前に−２０℃に冷凍保存された。
【００２９】
クロマトグラフィ
クロマトグラフィは、カラムオーブンを備えたＷａｔｅｒｓ社のＡｌｌｉａｎｃｅ  ２７９５ＨＴおよびＷａｔｅｒｓ社の２９９６ＰＤＡ検出器を使用して実施した。ラットの尿（蒸留水で１：４に希釈）を１０ｃｍ×２．１ｍｍのＷａｔｅｒｓ社Ｓｙｍｍｅｔｒｙ  Ｃ１８３．５μｍカラムに２０μＬ／分で注入した。カラムは４０℃に保持し、流量６００μＬ／分において勾配条件下、溶離液「Ａ」は０．１％蟻酸水溶液、溶離液「Ｂ」はアセトニトリルとして溶離した。カラムは、最初の１分間１００％Ａに保持し、その後、移動相を次の７分間にわたって３０％Ｂまで急激に増加させた。その後、移動相を０．１分間にわたって９５％Ｂまで急激に増加させた。この移動相は、カラムを洗浄するために０．５分間保持され、その後０．１分間にわたって最初の状態に戻された。最初の注入の１０分後に次の注入が行われた。
【００３０】
質量分析
質量分析の解析は、エレクトロスプレーインターフェース（「ＥＳＴ」）を備えたＭｉｃｒｏＭａｓｓ社Ｑｕａｔｔｒｏ  Ｍｉｃｒｏで行われた。この計測器は陰イオンモードでキャピラリ電圧３ｋＶの状態で操作され、コーン電圧は２５Ｖに設定、ネブライザガスは脱溶媒温度１５０℃およびガス供給温度７０℃の状態で６００Ｌ／ｈｒに設定された。カラム溶離液は５：１の比率に分割され、その結果１００Ｌ／ｈｒが質量分析供給源に入った。この計測器はスキャン時間２００ミリ秒およびスキャン間遅れ時間５０ミリ秒で１００から８００ｍ／ｚまでスキャンするフルスキャンモードで操作された。データは０分から１０分まで集められた。クロマトグラムの最初の部分は導入されず廃棄した。
【００３１】
データ解析
段階Ｉおよび段階ＩＩの代謝産物などの任意の薬物関連生成物の存在について、ＬＣ／ＭＳクロマトグラムを検査した。個々のＬＣ／ＭＳクロマトグラムはその後それぞれ１０画分（１画分１００スキャン）に分割された。１００スキャンの個別クロマトグラム画分のそれぞれに対してＭＳ信号が組み合わされた。ピークリストおよびイオン存在量はその後Ｅｘｃｅｌファイルにエクスポートされた。次に、イオン信号が観測されない箇所にゼロが強度値として挿入されてデータ集合が整列され、その結果各データ集合が、１００から８００の各ｍ／ｚ整数に対して１つの強度値を有した。その後、社内のプログラムを使用してＥｘｃｅｌファイルがＰＣＡ分析用ＭａｔＬａｂに移された。
【００３２】
検討
ラットの尿は数千の化合物を含む可能性があり、その大部分は不明のままである。これら化合物の大多数は中程度の極性かまたは高い極性であり、したがって、尿サンプルの全成分の完全な溶離を行うには有機含有率が０から３０％の勾配を用いるだけでよい。陰イオンエレクトロスプレーイオン化（ＥＳＩ）ＭＳは、陽イオンＥＳＩよりも情報が多いデータ集合をもたらすので、ＭＳ検出モードとして陰イオンを選んだ。この研究で用いたクロマトグラフィのカラム長は、勾配時間、および溶離液流量は、サンプル処理量とクロマトグラフィ分解能の間で最良のバランスを得るように選んだ。その結果、サンプル処理量は１時間当たり６サンプル、ピークキャパシティは７８となった。このピークキャパシティは非常に大きくはないが、１００から８００ｍ／ｚをスキャンする質量分析法の使用ははるかに大きな有効ピークキャパシティをもたらす。本研究においては、単純な逆相クロマトグラフィシステムを用いた。アミノ酸や糖などの極性の高い多くの化合物は保持されないことが十分に了解されるであろう。このシステムを変更することによってこれら化合物の解明が可能になる。
【００３３】
対照ラットの尿と時点１で投与された動物の尿サンプルのベースピーク強度（「ＢＰＩ」）クロマトグラムの予備的な比較によって定性的な差異が明らかになっている（図２参照）。これらのクロマトグラムから確認できるように、投与されたサンプルグループのクロマトグラムには、対照グループのそれと比較した場合に高いレベルを示すいくつかのピークがある。これらのピークは、クロマトグラムの５〜９分の時間域で溶離する。これらピークの徹底した調査によって、ピークはいずれも投与した化合物に起因する代謝産物ではないことが確認された。したがってこれらのピークは、動物の代謝状態の変化によるものでなければならない。
【００３４】
化合物「Ａ」の場合の時点２のサンプルから得られたＢＰＩクロマトグラムは、時点１のサンプルと比較した場合、いくつかのピーク強度におけるいくつかの微妙な変化を示している（図３参照）。時点２のサンプルにおいて、５から７分の間に溶離する分析物は、時点１の同一分析物と比較した相対濃度を低下させているようである。化合物「Ｂ」および「Ｃ」のサンプルは、時点１と時点２の間の変動が少ないことを示している（図４および５参照）。これらの時点１と時点２のサンプルは昼と夜のサンプルなので、生体異物の投与が無くても２つのサンプルタイプの間には著しい変化が見られることが予期される。このことは、昼と夜の間の動物代謝の変化に関連した化合物は保持されなかったかまたは検出されなかったことを示唆している。
【００３５】
対照サンプルと比較して投与サンプルで変化したと識別されるピークを表１に挙げる。（ｍ／ｚ＝３３８におけるこれらイオンの１つの抽出されたイオンクロマトグラムおよびこのピークのＭＳスペクトルが図６に示されている。）
【００３６】
【表１】
2005528606.tif 000003
【００３７】
表１のデータは、対照サンプルおよび投与サンプルのこれらのピークについて観察されたピーク強度の変化を示している。このデータから観察できるように、ピーク濃度の変化は非常に著しい。これらの化合物は対照サンプル中にも観測されるという事実は、これらのピークが投与化合物の代謝生成物ではないことを立証している。
【００３８】
データ集合全体のＰＣＡ分析は図７に呈示されている。この図で、対照サンプルに関係するデータは丸で囲んだ部分内に含まれる。この対照サンプルは投与動物のサンプルから明らかに分けられている。このＰＣＡ分析は、保持時間およびｍ／ｚ値のデータ集合全体を使用して実施される。このクロマトグラムを１０個の時間域に細分化することによって、同様にクロマトグラムの特定時間域でＰＣＡ分析を行うことが可能になった。
【００３９】
５分から６分のデータのＰＣＡ分析が図８に呈示されている。この図では、丸で囲んだ部分内部に含まれている対照サンプルと投与サンプルとを分離することが可能である。１つに合わせたＭＳスペクトルは、投与サンプルおよび対照サンプルをＰＣＡプロットの異なる領域に位置づけることができるＭＳピークはｍ／ｚ＝２８３およびｍ／ｚ＝２６１であったことを示している。ｍ／ｚ＝２８３のイオンは投与サンプルの強度が増加したのに対して、ｍ／ｚ＝２６１のイオンは投与サンプルの強度が低下した。この結果は、ＬＣ／ＭＳデータを解析することによって得られた観測値を確認している。
【００４０】
図９は３分から４分のＬＣ／ＭＳデータのＰＣＡ分析を示している。このデータでは対照サンプルと投与サンプルは分離されていない。このデータは、クロマトグラムの領域内では、投与サンプルと対照サンプルの間に有意な差が無かったことを示している。
【００４１】
ＰＣＡデータ分析に結合されたＬＣ／ＭＳの応用は、３つの候補薬剤の投与に続くラットの尿の選別に使用することに成功した。この方法を使用して、対照サンプルと投与サンプルを区別することができた。ＰＣＡ分析を使用して、分離を担うＭＳスペクトルの各成分を特定することもできる。このデータは、薬の発見および開発におけるメタボノミックな応用のために、ＬＣ／ＭＳが陽子ＮＭＲに取って代わる実行可能な方法であるかまたはこれを補完さえするものであることを明確に実証している。
【００４２】
本発明をその実施形態に関連して特に示し、説明したが、本明細書に添付した請求項で定義したように本発明の精神と範囲から逸脱することなく形態と詳細の数多の変更が可能であることを当業者は了解するであろう。
【図面の簡単な説明】
【００４３】
【図１Ａ】本発明の実施形態に伴う各ステップを示す概略図である。
【図１Ｂ】本発明の実施形態に伴う各ステップを示す概略図である。
【図１Ｃ】本発明の実施形態に伴う各ステップを示す概略図である。
【図２Ａ】対照ラットおよび投与ラットに対する陰イオンＢＰＩＬＣ／ＭＳクロマトグラムの比較図である。
【図２Ｂ】対照ラットおよび投与ラットに対する陰イオンＢＰＩＬＣ／ＭＳクロマトグラムの比較図である。
【図２Ｃ】対照ラットおよび投与ラットに対する陰イオンＢＰＩＬＣ／ＭＳクロマトグラムの比較図である。
【図２Ｄ】対照ラットおよび投与ラットに対する陰イオンＢＰＩＬＣ／ＭＳクロマトグラムの比較図である。
【図３Ａ】化合物Ａを投与したラットの時点１および時点２における陰イオンＢＰＩＬＣ／ＭＳ尿クロマトグラムの比較図である。
【図３Ｂ】化合物Ａを投与したラットの時点１および時点２における陰イオンＢＰＩＬＣ／ＭＳ尿クロマトグラムの比較図である。
【図４Ａ】化合物Ｂを投与したラットの時点１および時点２における陰イオンＢＰＩＬＣ／ＭＳ尿クロマトグラムの比較図である。
【図４Ｂ】化合物Ｂを投与したラットの時点１および時点２における陰イオンＢＰＩＬＣ／ＭＳ尿クロマトグラムの比較図である。
【図５Ａ】化合物Ｃを投与したラットの時点１および時点２における陰イオンＢＰＩＬＣ／ＭＳ尿クロマトグラムの比較図である。
【図５Ｂ】化合物Ｃを投与したラットの時点１および時点２における陰イオンＢＰＩＬＣ／ＭＳ尿クロマトグラムの比較図である。
【図６Ａ】抜出したイオンクロマトグラムの図である。
【図６Ｂ】ｍ／ｚ３３８のピークについてのＭＳスペクトル図である。
【図７】ラットの尿サンプル１〜２４から成る全データ集合ＰＣＡプロット図である。
【図８Ａ】統合して重ね合わせた質量スペクトルの図である。
【図８Ｂ】５分から６分の部分のＰＣＡプロット図である。
【図９Ａ】統合して重ね合わせた質量スペクトルの図である。
【図９Ｂ】３分から４分の部分のＰＣＡプロット図である。
【図１Ａ】
2005528606.tif 000004
【図１Ｂ】
2005528606.tif 000005
【図１Ｃ】
2005528606.tif 000006
【図２Ａ】
2005528606.tif 000007
【図２Ｂ】
2005528606.tif 000008
【図２Ｃ】
2005528606.tif 000009
【図２Ｄ】
2005528606.tif 000010
【図３Ａ】
2005528606.tif 000011
【図３Ｂ】
2005528606.tif 000012
【図４Ａ】
2005528606.tif 000013
【図４Ｂ】
2005528606.tif 000014
【図５Ａ】
2005528606.tif 000015
【図５Ｂ】
2005528606.tif 000016
【図６Ａ】
2005528606.tif 000017
【図６Ｂ】
2005528606.tif 000018
【図７】
2005528606.tif 000019
【図８Ａ】
2005528606.tif 000020
【図８Ｂ】
2005528606.tif 000021
【図９Ａ】
2005528606.tif 000022
【図９Ｂ】
2005528606.tif 000023
【国際調査報告】
2005528606.tif 000024