JP6620894B2 - クロマトグラフ質量分析用データ解析装置 - Google Patents
クロマトグラフ質量分析用データ解析装置 Download PDFInfo
- Publication number
- JP6620894B2 JP6620894B2 JP2018549668A JP2018549668A JP6620894B2 JP 6620894 B2 JP6620894 B2 JP 6620894B2 JP 2018549668 A JP2018549668 A JP 2018549668A JP 2018549668 A JP2018549668 A JP 2018549668A JP 6620894 B2 JP6620894 B2 JP 6620894B2
- Authority
- JP
- Japan
- Prior art keywords
- mass
- data
- chromatograph
- detector
- charge ratio
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004949 mass spectrometry Methods 0.000 title claims description 16
- 238000012545 processing Methods 0.000 claims description 27
- 238000001819 mass spectrum Methods 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 24
- 238000000034 method Methods 0.000 claims description 19
- 238000001228 spectrum Methods 0.000 claims description 19
- 238000004458 analytical method Methods 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 17
- 238000007405 data analysis Methods 0.000 claims description 15
- 239000007788 liquid Substances 0.000 claims description 13
- 230000002123 temporal effect Effects 0.000 claims description 12
- 239000000284 extract Substances 0.000 claims description 9
- 239000000126 substance Substances 0.000 description 19
- 238000005259 measurement Methods 0.000 description 15
- 238000000491 multivariate analysis Methods 0.000 description 12
- 230000014759 maintenance of location Effects 0.000 description 10
- 238000013500 data storage Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 9
- 239000011159 matrix material Substances 0.000 description 9
- 150000001875 compounds Chemical class 0.000 description 7
- 238000000605 extraction Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000002835 absorbance Methods 0.000 description 5
- 238000007619 statistical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 238000000862 absorption spectrum Methods 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 238000001601 dielectric barrier discharge ionisation Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000012628 principal component regression Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
- G01N30/78—Detectors specially adapted therefor using more than one detector
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
- G01N30/72—Mass spectrometers
- G01N30/7206—Mass spectrometers interfaced to gas chromatograph
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
- G01N27/64—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode using wave or particle radiation to ionise a gas, e.g. in an ionisation chamber
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
- G01N30/72—Mass spectrometers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
- G01N30/74—Optical detectors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Toxicology (AREA)
- Library & Information Science (AREA)
- Engineering & Computer Science (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Description
本発明は、液体クロマトグラフ(LC)やガスクロマトグラフ(GC)等のクロマトグラフと質量分析装置とを組み合わせたクロマトグラフ質量分析装置で得られたデータを解析するデータ解析装置に関し、さらに詳しくは、クロマトグラフからの試料中の成分を検出する質量分析装置とは異なる検出器を併設したクロマトグラフ質量分析装置のためのデータ解析装置に関する。
液体クロマトグラフ質量分析装置(LC−MS)は液体クロマトグラフの検出器として質量分析装置を用いたものであるが、これにさらに、紫外可視分光検出器(UV検出器)やフォトダイオードアレイ検出器(PDA検出器)などの別の検出器を併用する場合がある(特許文献1、2等参照)。ガスクロマトグラフの検出器として質量分析装置を用いたガスクロマトグラフ質量分析装置(GC−MS)においても同様に、光イオン化検出器(PID)、水素炎イオン化検出器(FID)などの別の検出器を質量分析装置と併用する場合がある。以下、本明細書では必要に応じて、LC−MSやGC−MSにおいて、主検出器である質量分析装置と併用される検出器をその方式に依らず副検出器ということとする。また、以下の説明では、副検出器を備えたLC−MSやGC−MSを単にLC−MSやGC−MSという。
例えばUV検出器やPDA検出器では波長毎に試料中の成分の吸光度の時間的変化を測定することができる。分子中の特定の化学構造は或る波長を特異的に吸収する場合がある。こうした場合、特定の波長における吸光度の時間的な変化をUV検出器やPDA検出器で測定し、その変化を示すクロマトグラムを作成すると、該クロマトグラムには上記特定の化学構造を有する分子種由来のピークが現れる。そこで、上記LC−MSを利用して、こうしたクロマトグラムに現れるピークを手掛かりとし、特定の化学構造を有する等の着目する分子種に関連する質量情報をマススペクトルから抽出する解析が従来から行われている。
一般に、こうした解析を行う際には、まず、UV検出器等の副検出器による検出信号に基づいて作成されたクロマトグラムにおいてピークを検出する。そして、収集したマススペクトルデータの中で、上記検出したピークに対応する時間範囲に含まれるマススペクトルデータを時間方向に積算する又は平均化することで積算マススペクトルを求める。こうして求めた積算マススペクトルにおいて特徴的な質量電荷比を探索し、保持時間やそうした質量電荷比情報を参照して、試料に含まれる特定の化学構造を有する分子種を同定する。
上述したようにクロマトグラムにおいてピークを検出するには様々なアルゴリズムがあるが、基本的に、ピーク検出に必要な各種のパラメータは操作者(ユーザー)が手動で設定する必要がある。ピークを適切に検出できるパラメータを決める作業はかなり面倒であるし、そうした作業には操作者の技量や経験も要求される。一方で、最近のデータ解析作業においては、そうした作業に関する知識や経験のない操作者であっても適切な解析が自動的に行えることが強く求められている。
本発明はこうした課題に鑑みて成されたものであり、その目的とするところは、副検出器で得られる情報を利用して特定の化学構造を有する等の特徴的な分子種に関連する質量情報を確実に且つ効率良く抽出することができるクロマトグラフ質量分析用データ解析装置を提供することである。
上記課題を解決するために成された本発明は、試料中の成分を時間的に分離するクロマトグラフと、該クロマトグラフで分離された成分を検出する質量分析装置である主検出器と、前記クロマトグラフで分離された成分を検出する質量分析装置とは異なる副検出器と、を具備するクロマトグラフ質量分析装置で収集されたデータを解析するクロマトグラフ質量分析用データ解析装置であって、
a)前記主検出器において繰り返し得られたマススペクトルデータと前記副検出器で得られたクロマトグラムデータとに対して、該クロマトグラムデータにおける信号強度値の時間的変化と前記マスクロマトグラムデータにおける質量電荷比毎の信号強度値の時間的変化との相関性を調べる演算を行う演算処理部と、
b)前記演算処理部での演算により質量電荷比値毎に得られる相関性に関連する指標値に基づいて、相関性が相対的に高い質量電荷比値を抽出する質量電荷比値抽出部と、
を備えることを特徴としている。
a)前記主検出器において繰り返し得られたマススペクトルデータと前記副検出器で得られたクロマトグラムデータとに対して、該クロマトグラムデータにおける信号強度値の時間的変化と前記マスクロマトグラムデータにおける質量電荷比毎の信号強度値の時間的変化との相関性を調べる演算を行う演算処理部と、
b)前記演算処理部での演算により質量電荷比値毎に得られる相関性に関連する指標値に基づいて、相関性が相対的に高い質量電荷比値を抽出する質量電荷比値抽出部と、
を備えることを特徴としている。
本発明に係るクロマトグラフ質量分析用データ解析装置において、クロマトグラフは典型的には液体クロマトグラフ又はガスクロマトグラフである。クロマトグラフが液体クロマトグラフである場合、副検出器は典型的には、フォトダイオードアレイ検出器や紫外可視分光検出器等の、液体クロマトグラフのカラムからの溶出液中に含まれる成分を順次そのまま検出する、つまりは溶出液に含まれる成分の性質や特性に応じた信号の時間的変化を測定する検出器であるが、質量分析装置と同様に、液体試料を微細化して又は気化させて検出する検出器でもよい。
また、クロマトグラフがガスクロマトグラフである場合、副検出器としては、光イオン化検出器、水素炎イオン化検出器、或いは、誘電体バリア放電イオン化検出器などを用いるとよい。
また、クロマトグラフがガスクロマトグラフである場合、副検出器としては、光イオン化検出器、水素炎イオン化検出器、或いは、誘電体バリア放電イオン化検出器などを用いるとよい。
本発明に係るクロマトグラフ質量分析用データ解析装置において、演算処理部は、主検出器である質量分析装置で繰り返し得られた質量電荷比、時間、信号強度という三つの次元を有するマススペクトルデータと、副検出器で得られた時間、信号強度という二つの次元を有するクロマトグラムデータとに対し、所定の統計的解析処理を実施する。より詳しくは、マススペクトルデータとクロマトグラムデータとの間で同じ測定時刻に得られたデータ同士の統計的解析処理を行い、信号強度値の時間的変化の類似性又は相違性を示す指標値を算出する。こうした目的のために利用できる最も単純な統計的解析処理は、クロマトグラムデータにおける信号強度値の時間的変化とマスクロマトグラムデータにおける質量電荷比毎の信号強度値の時間的変化との相関係数を算出する相関分析法である。
また、実質的に相関性を求める手法である多変量解析手法、例えば部分的最小二乗法(PLS=Partial Least Squares)や判別分析などを用いてもよい。即ち、本発明に係るクロマトグラフ質量分析用データ解析装置の一実施態様として、前記演算処理部は部分的最小二乗法による演算を行い、前記指標値として回帰係数を求める構成とすることができる。
上記演算処理部での統計的解析処理によって、質量電荷比値毎に相関性に関連する指標値が得られる。そこで、質量電荷比値抽出部は例えば、その指標値を予め設定された閾値と比較することにより、相関性が相対的に高い質量電荷比値を抽出する。例えば副検出器が特定の化学構造を特異的に検出するものである場合、該副検出器から得られるクロマトグラムデータに基づいてクロマトグラムを作成すると、その化学構造を有する成分(分子)が副検出器に到達する時間付近にピークを有する波形となる。したがって、このクロマトグラムとの相関が高い波形形状が得られる質量電荷比は、上記化学構造を有する成分である可能性が高い。そこで、操作者は質量電荷比値抽出部で抽出された質量電荷比に絞って抽出イオンクロマトグラムの波形を確認したり同じ保持時間に検出される複数の質量電荷比情報を集めたりすることで、上記化学構造を有する成分の同定等を効率良く行うことができる。
また本発明に係るクロマトグラフ質量分析用データ解析装置では、
前記マススペクトルデータに基づいて前記質量電荷比値抽出部により抽出された質量電荷比における抽出イオンクロマトグラムを作成する抽出イオンクロマトグラム作成部と、
前記抽出イオンクロマトグラムと前記クロマトグラムデータに基づくクロマトグラムとを表示部の画面上に表示する表示処理部と、
をさらに備える構成とすることができる。
前記マススペクトルデータに基づいて前記質量電荷比値抽出部により抽出された質量電荷比における抽出イオンクロマトグラムを作成する抽出イオンクロマトグラム作成部と、
前記抽出イオンクロマトグラムと前記クロマトグラムデータに基づくクロマトグラムとを表示部の画面上に表示する表示処理部と、
をさらに備える構成とすることができる。
この構成によれば、操作者が着目する又は確認すべき抽出イオンクロマトグラムが自動的に作成され、副検出器による検出信号に基づくクロマトグラムとともに表示されるので、上記のような特定の成分についての解析作業における操作者の負担が一層軽減される。
また本発明に係るクロマトグラフ質量分析用データ解析装置では、
質量電荷比値毎に得られる前記指標値と質量電荷比との関係を示す指標値スペクトルを作成する指標値スペクトル作成部と、
前記指標値スペクトルを表示部の画面上に表示する表示処理部と、
をさらに備える構成とすることもできる。
質量電荷比値毎に得られる前記指標値と質量電荷比との関係を示す指標値スペクトルを作成する指標値スペクトル作成部と、
前記指標値スペクトルを表示部の画面上に表示する表示処理部と、
をさらに備える構成とすることもできる。
この構成によれば、操作者は表示された指標値スペクトルにより、全ての質量電荷比に対応する指標値を一目で確認することができる。それにより、操作者は相関性が相対的に高い質量電荷比値を自ら判断して抽出することができる。また、上述したように質量電荷比値抽出部により自動的に相関性が相対的に高い質量電荷比値を抽出する場合においては、指標値が閾値に達せずに自動的に抽出されなかったものの、該閾値に近い指標値を示す質量電荷比値を操作者が見つけることができる。
本発明に係るクロマトグラフ質量分析用データ解析装置によれば、紫外可視分光検出器やフォトダイオードアレイ検出器などの質量分析装置とは異なる方式である副検出器で得られる情報を利用して、特定の化学構造を有する等の特徴的な分子種に関連する質量情報を的確に抽出することができる。また、そうした解析に際して操作者による面倒な作業が軽減されるので、解析を効率的に進めることができるとともに、操作者の技量や経験に依存しない信頼性、再現性の高い解析が可能となる。
以下、本発明に係るデータ解析装置を用いたLC−MSシステムの一実施例について、添付図面を参照して説明する。
図1は本実施例のLC−MSシステムの概略構成図である。
図1は本実施例のLC−MSシステムの概略構成図である。
本実施例のLC−MSシステムは、液体クロマトグラフ(LC)部1と、PDA検出器2と、質量分析計3と、データ処理部4と、制御部5と、入力部6と、表示部7と、を含む。
LC部1は、移動相が貯留される移動相容器11と、移動相容器11から移動相を吸引して略一定流速で送給する送液ポンプ12と、移動相中に所定量の試料を注入するインジェクタ13と、試料中の成分を時間的に分離するカラム14と、を含む。カラム14出口から溶出した溶出液はPDA検出器2の試料セル(図示せず)を経て質量分析計3に導入される。PDA検出器2で得られた検出信号及び質量分析計3で得られた検出信号はそれぞれ、図示しないアナログデジタル変換器(ADC)でデジタルデータに変換されてデータ処理部4に入力される。なお、一般に、PDA検出器2の試料セルを通過した溶出液の全量を質量分析計3で測定することは難しいため、試料セルを通過した溶出液の一部をスプリッタ等により分岐して質量分析計3に導入する構成が採られる。
LC部1は、移動相が貯留される移動相容器11と、移動相容器11から移動相を吸引して略一定流速で送給する送液ポンプ12と、移動相中に所定量の試料を注入するインジェクタ13と、試料中の成分を時間的に分離するカラム14と、を含む。カラム14出口から溶出した溶出液はPDA検出器2の試料セル(図示せず)を経て質量分析計3に導入される。PDA検出器2で得られた検出信号及び質量分析計3で得られた検出信号はそれぞれ、図示しないアナログデジタル変換器(ADC)でデジタルデータに変換されてデータ処理部4に入力される。なお、一般に、PDA検出器2の試料セルを通過した溶出液の全量を質量分析計3で測定することは難しいため、試料セルを通過した溶出液の一部をスプリッタ等により分岐して質量分析計3に導入する構成が採られる。
データ処理部4は、後述する特徴的な処理を実施するために、クロマトグラムデータ格納部41、マススペクトルデータ格納部42、多変量解析演算部43、m/z値抽出部44、係数スペクトル作成部45、XIC作成部46、解析結果表示処理部47、などの機能ブロックを備える。
なお、一般に、データ処理部4及び制御部5の実体はパーソナルコンピュータである。そして、このコンピュータにインストールされた専用の制御・処理ソフトウェアを該コンピュータ上で実行することにより、上記各機能ブロックの機能が実現されるものとすることができる。
なお、一般に、データ処理部4及び制御部5の実体はパーソナルコンピュータである。そして、このコンピュータにインストールされた専用の制御・処理ソフトウェアを該コンピュータ上で実行することにより、上記各機能ブロックの機能が実現されるものとすることができる。
本実施例のLC−MSシステムでは、以下のようにして測定が実施されてデータが収集される。
制御部5の制御の下で、送液ポンプ12は一定流速で以て移動相をカラム14に送給する。所定のタイミングでインジェクタ13から多数の成分を含む試料が移動相中に注入されると、該試料は移動相の流れに乗ってカラム14に導入され、カラム14を通過する間に試料中の成分は時間的に分離される。カラム14の出口から溶出する溶出液には、カラム14中で分離された成分が含まれる。PDA検出器2は、図示しないものの、試料セルのほか、光源、分光器、PDA検出素子などを含む。そして、光源から発した光を試料セル中を通過する溶出液に照射し、溶出液を透過した透過光を分光器で波長分散してPDA検出素子で波長毎の光強度を略一斉に検出する。これにより、所定の波長範囲に亘る吸光特性を反映した検出信号が一度に得られる。この検出信号が所定時間間隔でサンプリング及びデジタル化されデータ処理部4に送られ、クロマトグラムデータ格納部41に格納される。
制御部5の制御の下で、送液ポンプ12は一定流速で以て移動相をカラム14に送給する。所定のタイミングでインジェクタ13から多数の成分を含む試料が移動相中に注入されると、該試料は移動相の流れに乗ってカラム14に導入され、カラム14を通過する間に試料中の成分は時間的に分離される。カラム14の出口から溶出する溶出液には、カラム14中で分離された成分が含まれる。PDA検出器2は、図示しないものの、試料セルのほか、光源、分光器、PDA検出素子などを含む。そして、光源から発した光を試料セル中を通過する溶出液に照射し、溶出液を透過した透過光を分光器で波長分散してPDA検出素子で波長毎の光強度を略一斉に検出する。これにより、所定の波長範囲に亘る吸光特性を反映した検出信号が一度に得られる。この検出信号が所定時間間隔でサンプリング及びデジタル化されデータ処理部4に送られ、クロマトグラムデータ格納部41に格納される。
PDA検出器2の試料セルを通過した溶出液は質量分析計3に導入される。質量分析計3では所定の質量電荷比範囲に亘るスキャン測定が繰り返し行われ、1回のスキャン測定毎に上記質量電荷比範囲に亘るマススペクトルデータが得られる。このマススペクトルデータはデータ処理部4に送られ、マススペクトルデータ格納部42に格納される。
本実施例のLC−MSシステムにおいて、特定の化学構造を含む化合物を網羅的に解析する場合の特徴的なデータ解析処理について、図2、図3を参照して説明する。図2はこのデータ解析処理の手順を示すフローチャート、図3はデータ解析処理の概略的な説明図である。
上述したように、特定の部分化学構造を含む化合物は、その部分化学構造に特異的な特定の波長λ1で高い吸光度を示す場合がある。そこで操作者は、着目する部分化学構造を有する化合物に対応する特定の波長λ1の情報を入力部6から指定する(ステップS1)。もちろん、波長λ1の値を入力する代わりに、予め用意された部分化学構造のリストの中から操作者が適宜のものを選択すると、その選択操作を受けて、対応する波長λ1が自動的に設定されるようにしてもよい。
多変量解析演算部43はクロマトグラムデータ格納部41から、ステップS1で指定された又は自動的に設定された波長λ1における信号強度の時間変化を示すクロマトグラムデータを選択的に読み出す(ステップS2)。なお、測定開始時点から測定終了時点までの全ての測定時間範囲におけるクロマトグラムデータを読み出すのではなく、予め(例えばステップS1と同時に)操作者が指定した測定時間範囲のみにおけるクロマトグラムデータを読み出すようにしてもよい。なお、以下の説明では、ここで読み出されたクロマトグラムデータを基準クロマトグラムデータ、該データに基づいて作成される又は形成されるクロマトグラムを基準クロマトグラムという。
また多変量解析演算部43はマススペクトルデータ格納部42から、ステップS2で読み出したクロマトグラムデータに対応する測定時間範囲におけるマススペクトルデータを読み出す(ステップS3)。通常、或る成分を含む溶出液がPDA検出器2において試料セルを通過する時点つまりはPDA検出器2での検出時点から、質量分析計3に導入されて検出される時点までには時間遅れtdがある。そこで、両データにおける保持時間の比較が行えるように時間遅れtdを補正するとよい。
時間遅れtdは溶出液の流路の容積のほか、移動相の流速などのパラメータに依存する。これらパラメータは測定の実行時に確定しているから、時間遅れtdはそれらパラメータに基づいて計算可能である。そこで、そうして時間遅れtdを計算したあと、その時間遅れtdの分だけ基準クロマトグラムデータの時間軸を全体的に遅らせる、或いは、繰り返し得られるマススペクトルデータの時間軸を全体的に進めるように調整することで、時間遅れtdがゼロになるように両データの時間軸を合わせることができる。
なお、データ処理に際して時間遅れtdを補正するのではなく、データを取得する段階で、つまりはデータ格納部41、42にそれぞれデータを格納する段階で時間遅れtdを補正するようにしてもよい。そうした補正が行われている場合には、上述したようなデータ読み出し時における時間遅れtdの補正は不要である。
続いて多変量解析演算部43は、多変量解析の一手法であるPLSを用い、基準クロマトグラムデータにおける信号強度値の時間的変化と、繰り返し得られたマススペクトルデータの中の各質量電荷比における信号強度値の時間的変化、つまりは抽出イオンクロマトグラムにおける信号強度値の時間的変化との相関性を調べる。
具体的には、多変量解析演算部43はまず、基準クロマトグラムデータにおける各測定時刻毎の信号強度値データと、マススペクトルデータにおける質量電荷比毎の信号強度値データとをそれぞれ行列の形式で表す。図3に示すように、基準クロマトグラムデータにおける信号強度値は測定時刻毎の信号強度値y1、y2、…、ynを並べた1次元の行列Yであり、マススペクトルデータにおける信号強度値は測定時刻毎及び質量電荷比値毎の信号強度値x11、x12、…、x1m、x21、x22、…、x2m、…、xn1、xn2、…、xnmを2次元的に並べた2次元行列Xで表すことができる。
そのあと多変量解析演算部43は、マススペクトルデータにおける信号強度値の2次元行列Xを説明変数(入力変数)、基準クロマトグラムデータにおける信号強度値の1次元行列Yを被説明変数(出力変数)とし、周知であるPLSの演算を実行して回帰係数行列を算出する(ステップS4)。この回帰係数行列の要素数はmであり、これは質量電荷比値の回帰係数が並べられた行列である。この回帰係数は、基準クロマトグラムデータに基づく信号強度値データを教師データとした、基準クロマトグラムデータとマススペクトルデータとの質量電荷比毎の相関の高さを示している。したがって、回帰係数が大きい質量電荷比値は、マススペクトルにおける信号強度の時間的変化が基準クロマトグラムデータに基づく信号強度の時間的変化に近い質量電荷比である。
なお、相関性を示す指標値を算出する多変量解析であれば、PLS以外の多変量解析、例えば最小二乗回帰(OLS)、主成分分析(PCA)、主成分回帰(PCR)などを用いてもよい。また、より単純な相関分析法を用いて相関係数を求めるようにしてもよい。
m/z値抽出部44は、上述したように質量電荷比毎に求まった回帰係数を所定の閾値と比較し、回帰係数が閾値以上である質量電荷比を抽出する(ステップS5)。このときの判定基準である閾値は装置内部で定められているデフォルト値でもよいが、操作者が任意に設定できるようにすることが好ましい。また、必要に応じて係数スペクトル作成部45は、回帰係数行列に基づいて、横軸を質量電荷比、縦軸を回帰係数値とした係数スペクトルを作成する。もちろん、m/z値抽出部44は、係数スペクトル作成部45により作成された係数スペクトル上で回帰係数が閾値以上である質量電荷比を抽出してもよい。
XIC作成部46はステップS5で抽出された質量電荷比における抽出イオンクロマトグラム(XIC)を作成する(ステップS6)。そして、解析結果表示処理部47は解析結果として、ステップS2において取得された基準クロマトグラムデータに基づいて作成される基準クロマトグラムとステップS6において作成された抽出イオンクロマトグラムとを共に表示部7の画面上に表示する(ステップS7)。このとき、例えば時間軸を合わせ各クロマトグラムの表示色を異なるものとして重ねて表示するとよい。また、必要に応じて、回帰係数スペクトルも表示部7の画面上に表示するとよい。
ステップS2において取得された基準クロマトグラムデータに基づいて作成される基準クロマトグラムには、操作者が着目する部分化学構造を有する化合物に対応するピークが現れる筈である。したがって、ステップS6で作成された抽出イオンクロマトグラムにも同じ部分化学構造を有する化合物に対応するピークが現れる可能性がある。上述したように、そうした抽出イオンクロマトグラムと基準クロマトグラムとを例えば異なる表示色で以て重ねて表示することで、操作者は両方のクロマトグラムにおいて同じ保持時間に現れるピークを容易に見つけることができ、該ピークの保持時間等に基づいて化合物を特定する等の解析作業を行うことができる。
上記実施例では、保持時間(測定時間範囲)や質量電荷比を指定することなく回帰係数を計算していたが、その場合、図3に示した行列X、Yにおける行数、列数は非常に多くなり、PLS演算の計算量も膨大となる。その結果、計算結果が得られるまでに時間がかかったり高性能なコンピュータが必要になったりするおそれがある。そこで、ステップS4の実行前に、ユーザーが回帰係数を計算する対象の保持時間範囲や質量電荷比値又は質量電荷比範囲を指定し、多変量解析演算部43はその指定された保持時間範囲、及び/又は質量電荷比値若しくは質量電荷比範囲に含まれるデータに限定したPLS演算を行って回帰係数を算出するようにしてもよい。これにより、PLS演算等の計算量を削減することができる。また、ユーザーが興味のある保持時間や質量電荷比に対する結果のみを取得できるため、興味のない結果を確認する必要がなくなり、解析の効率を向上させることができる。
また、全ての保持時間範囲や質量電荷比範囲についてのPLS演算を一度に行うのではなく、保持時間範囲や質量電荷比範囲をそれぞれ複数の範囲に区切り、その区切られた保持時間範囲毎及び/又は質量電荷比範囲毎にPLS演算を実行して回帰係数を求めるようにしてもよい。この場合、区切られた範囲に対する回帰係数が求まる度にその結果を表示すれば、ユーザーが興味のある結果を確認した時点で、それ以降の演算処理を途中で打ち切るようにすることができる。なお、保持時間範囲や質量電荷比範囲を区切る際の区切り(保持時間幅等)はユーザーにより指定されるようにしてもよいし、自動的に決まるようにしてもよい。
上記実施例のLC−MSシステムでは副検出器としてPDA検出器を用いていたが、これはUV検出器でもよい。ただし、一般に、UV検出器の場合、PDA検出器とは異なり分光器で機械的な波長走査を行うため、或る波長範囲に亘る吸光スペクトルを測定すると時間が掛かる。そこで、測定を行う前に、つまりデータ取得前に操作者が着目する波長を一又は複数入力し、その入力された一又は複数の波長のみにおける吸光度を示すデータを取得するようにする、つまりは所定波長範囲全体の吸光度スペクトルを測定するのではなく特定の波長のみの吸光度を測定するとよい。また、副検出器として、PDA検出器、UV検出器以外の検出器、例えば蛍光検出器などを用いてもよいし、さらには核磁気共鳴装置(NMR)を副検出器として用いてもよい。
また本発明に係るデータ解析装置はGC−MSシステムにも適用することができる。その場合、副検出器としては、光イオン化検出器、水素炎イオン化検出器、或いは、誘電体バリア放電イオン化検出器など、特徴的な化学構造や物質を有する化合物を検出可能な検出器を用いるとよい。
また、本発明に係るデータ解析装置による解析に用いられる基準クロマトグラムは、一般的なGCやLCでなく、GC×GC(包括的2次元GC)やLC×LC(包括的2次元LC)等の包括的2次元クロマトグラフ、或いは超臨界流体クロマトグラフ等により取得されたものでもよい。
また、本発明に係るデータ解析装置による解析に用いられる基準クロマトグラムは、一般的なGCやLCでなく、GC×GC(包括的2次元GC)やLC×LC(包括的2次元LC)等の包括的2次元クロマトグラフ、或いは超臨界流体クロマトグラフ等により取得されたものでもよい。
また、質量分析計3で得られるマススペクトルは、一般的なマススペクトル(MS1スペクトル)に限らず、MSnスペクトル(又はMS/MSスペクトル)、電子イオン化(EI)法によるイオン源を備えた質量分析計で得られるようなISD(In Source Decay)スペクトル、米国ウォーターズ(Waters)社によるMSEスペクトルなどでもよい。
さらにまた、上記実施例は本発明の一例にすぎず、上記記載した以外の点において、本発明の趣旨の範囲で適宜変形、修正、追加などを行っても本願特許請求の範囲に包含されることは明らかである。
1…液体クロマトグラフ(LC)部
11…移動相容器
12…送液ポンプ
13…インジェクタ
14…カラム
2…PDA検出器
3…質量分析計
4…データ処理部
41…クロマトグラムデータ格納部
42…マススペクトルデータ格納部
43…多変量解析演算部
44…m/z値抽出部
45…係数スペクトル作成部
46…XIC作成部
47…解析結果表示処理部
5…制御部
6…入力部
7…表示部
11…移動相容器
12…送液ポンプ
13…インジェクタ
14…カラム
2…PDA検出器
3…質量分析計
4…データ処理部
41…クロマトグラムデータ格納部
42…マススペクトルデータ格納部
43…多変量解析演算部
44…m/z値抽出部
45…係数スペクトル作成部
46…XIC作成部
47…解析結果表示処理部
5…制御部
6…入力部
7…表示部
Claims (5)
- 試料中の成分を時間的に分離するクロマトグラフと、該クロマトグラフで分離された成分を含む試料中の成分を検出する質量分析装置である主検出器と、前記クロマトグラフで分離された成分を含む試料中の成分を検出する質量分析装置とは異なる副検出器と、を具備するクロマトグラフ質量分析装置で収集されたデータを解析するクロマトグラフ質量分析用データ解析装置であって、
a)前記主検出器で繰り返し得られたマススペクトルデータと前記副検出器で得られたクロマトグラムデータとに対して、該クロマトグラムデータにおける信号強度値の時間的変化と前記マススペクトルデータにおける質量電荷比毎の信号強度値の時間的変化との相関性を調べる演算を行う演算処理部と、
b)前記演算処理部での演算により質量電荷比値毎に得られる相関性に関連する指標値に基づいて、相関性が相対的に高い質量電荷比値を抽出する質量電荷比値抽出部と、
を備えることを特徴とするクロマトグラフ質量分析用データ解析装置。 - 請求項1に記載のクロマトグラフ質量分析用データ解析装置であって、
前記演算処理部は部分的最小二乗法による演算を行い、前記指標値として回帰係数を求めることを特徴とするクロマトグラフ質量分析用データ解析装置。 - 請求項1に記載のクロマトグラフ質量分析用データ解析装置であって、
前記クロマトグラフは液体クロマトグラフであり、前記副検出器はフォトダイオードアレイ検出器又は紫外可視分光検出器であることを特徴とするクロマトグラフ質量分析用データ解析装置。 - 請求項1に記載のクロマトグラフ質量分析用データ解析装置であって、
前記マススペクトルデータに基づいて前記質量電荷比値抽出部により抽出された質量電荷比における抽出イオンクロマトグラムを作成する抽出イオンクロマトグラム作成部と、
前記抽出イオンクロマトグラムと前記クロマトグラムデータに基づくクロマトグラムとを表示部の画面上に表示する表示処理部と、
をさらに備えることを特徴とするクロマトグラフ質量分析用データ解析装置。 - 請求項1に記載のクロマトグラフ質量分析用データ解析装置であって、
質量電荷比値毎に得られる前記指標値と質量電荷比との関係を示す指標値スペクトルを作成する指標値スペクトル作成部と、
前記指標値スペクトルを表示部の画面上に表示する表示処理部と、
をさらに備えることを特徴とするクロマトグラフ質量分析用データ解析装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2016/083173 WO2018087824A1 (ja) | 2016-11-09 | 2016-11-09 | クロマトグラフ質量分析用データ解析装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2018087824A1 JPWO2018087824A1 (ja) | 2019-04-04 |
JP6620894B2 true JP6620894B2 (ja) | 2019-12-18 |
Family
ID=62109721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018549668A Active JP6620894B2 (ja) | 2016-11-09 | 2016-11-09 | クロマトグラフ質量分析用データ解析装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11215589B2 (ja) |
EP (1) | EP3540425B1 (ja) |
JP (1) | JP6620894B2 (ja) |
CN (1) | CN109983333B (ja) |
WO (1) | WO2018087824A1 (ja) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018071415A1 (en) * | 2016-10-10 | 2018-04-19 | Purdue Research Foundation | Continuous roll-to-roll freeze-drying system and process |
US11959898B2 (en) * | 2018-08-31 | 2024-04-16 | Dh Technologies Development Pte. Ltd. | Identification and scoring of related compounds in complex samples |
KR102214705B1 (ko) * | 2019-02-11 | 2021-02-09 | 한양대학교 에리카산학협력단 | 공정 환경 내에 존재하는 물질의 분석 방법 및 이를 위한 장치 |
US11940426B2 (en) | 2019-03-27 | 2024-03-26 | Shimadzu Corporation | Chromatograph mass spectrometer |
JP7375640B2 (ja) * | 2020-03-23 | 2023-11-08 | 株式会社島津製作所 | イメージング質量分析システム、及び、イメージング質量分析を利用した分析方法 |
JP7460988B2 (ja) * | 2021-01-22 | 2024-04-03 | 株式会社島津製作所 | クロマトグラフ質量分析装置 |
JP2022181723A (ja) * | 2021-05-26 | 2022-12-08 | 株式会社島津製作所 | 分析方法および診断支援方法 |
CN113567605B (zh) * | 2021-08-16 | 2024-09-20 | 苏恺明 | 质量色谱图的自动化解释模型构建方法、装置和电子设备 |
CN114324713B (zh) * | 2022-01-13 | 2023-01-13 | 宁夏医科大学 | Uhplc-hrms数据依赖性采集的信息解析方法 |
CN115963074B (zh) * | 2023-02-23 | 2023-06-02 | 中国人民解放军国防科技大学 | 一种微生物材料孢子菌丝占比的快速检测方法和系统 |
CN117007724B (zh) * | 2023-10-07 | 2023-12-22 | 广州盛安医学检验有限公司 | 一种基于云计算的高效液相串联质谱分析系统 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4353242A (en) * | 1980-12-16 | 1982-10-12 | University Of Utah Research Foundation | Multichannel detection and resolution of chromatographic peaks |
JPS585651A (ja) * | 1981-06-30 | 1983-01-13 | Shimadzu Corp | クロマトグラフデ−タ処理表示装置 |
NL1016034C2 (nl) * | 2000-08-03 | 2002-02-08 | Tno | Werkwijze en systeem voor het identificeren en kwantificeren van chemische componenten van een te onderzoeken mengsel van materialen. |
JP4782278B2 (ja) * | 2000-12-19 | 2011-09-28 | 株式会社島津製作所 | 液体クロマトグラフ質量分析計 |
JP2004045350A (ja) * | 2002-07-16 | 2004-02-12 | Horiba Ltd | ガス分析システムおよびガス分析方法 |
JP2005291715A (ja) | 2004-03-31 | 2005-10-20 | Shimadzu Corp | におい測定装置 |
US9063162B2 (en) * | 2006-05-26 | 2015-06-23 | Waters Technologies Corporation | Apparatus and method for performing mass spectroscopy |
JP5757264B2 (ja) * | 2012-03-16 | 2015-07-29 | 株式会社島津製作所 | クロマトグラフ質量分析データ処理装置 |
JP2014021083A (ja) | 2012-07-24 | 2014-02-03 | Hitachi High-Technologies Corp | 質量分析方法及び質量分析システム |
JP2014219317A (ja) | 2013-05-09 | 2014-11-20 | 株式会社島津製作所 | クロマトグラフ用データ処理装置 |
US10121643B2 (en) * | 2014-07-24 | 2018-11-06 | Shimadzu Corporation | Chromatography/mass spectrometry data processing device |
-
2016
- 2016-11-09 JP JP2018549668A patent/JP6620894B2/ja active Active
- 2016-11-09 EP EP16921395.6A patent/EP3540425B1/en active Active
- 2016-11-09 US US16/347,483 patent/US11215589B2/en active Active
- 2016-11-09 WO PCT/JP2016/083173 patent/WO2018087824A1/ja unknown
- 2016-11-09 CN CN201680090726.6A patent/CN109983333B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
EP3540425A1 (en) | 2019-09-18 |
US20190277817A1 (en) | 2019-09-12 |
JPWO2018087824A1 (ja) | 2019-04-04 |
CN109983333B (zh) | 2022-02-25 |
EP3540425B1 (en) | 2023-12-27 |
US11215589B2 (en) | 2022-01-04 |
CN109983333A (zh) | 2019-07-05 |
EP3540425A4 (en) | 2019-11-27 |
WO2018087824A1 (ja) | 2018-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6620894B2 (ja) | クロマトグラフ質量分析用データ解析装置 | |
CN108508122B (zh) | 色谱仪用数据处理装置 | |
JP6048590B2 (ja) | 包括的2次元クロマトグラフ用データ処理装置 | |
JP6573028B2 (ja) | データ処理装置 | |
JP7409462B2 (ja) | クロマトグラム表示装置 | |
JP5751380B2 (ja) | 多次元クロマトグラフ装置 | |
JPWO2012004855A1 (ja) | 分析データ処理方法及び装置 | |
JP5757264B2 (ja) | クロマトグラフ質量分析データ処理装置 | |
JP6245387B2 (ja) | 3次元スペクトルデータ処理装置及び処理方法 | |
US10317378B2 (en) | Data processing system and data processing method for chromatograph | |
JP6132073B2 (ja) | 包括的2次元クロマトグラフ用データ処理装置 | |
WO2014108992A1 (ja) | クロマトグラフ質量分析用データ処理装置 | |
JP2013068444A (ja) | 包括的2次元クロマトグラフ質量分析用データ処理装置 | |
JP5146344B2 (ja) | クロマトグラフ用データ処理装置 | |
JP6226823B2 (ja) | クロマトグラフ質量分析装置及びその制御方法 | |
JP6583159B2 (ja) | クロマトグラフによる濃度測定方法及びクロマトグラフ | |
WO2018158801A1 (ja) | スペクトルデータの特徴抽出装置および方法 | |
JP6269803B2 (ja) | クロマトグラフ用データ処理装置及び方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20181204 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20191023 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20191105 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6620894 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |