(57)【特許請求の範囲】
【請求項１】
少なくとも、工程（１）ないし（５）
（１）評価すべき多成分系薬剤および基準群を構成する他の同種の多成分系薬剤の三次元高速液体クロマトグラフィーのフィンガープリントデータを得る工程
（２）上記（１）で得られたフィンガープリントデータから、等間隔ごとのデータポイントを採取して圧縮したフィンガープリントデータを得る工程
（３）前記（２）の圧縮したフィンガープリントデータについて、多成分系薬剤の番号と、溶出時間又は検出波長とに対してＭＴ法における変数軸を割り振り、シグナル強度をＭＴ法における特徴量とする工程
（４）前記（３）の特徴量から、ＭＴ法を用いて単位空間を得る工程
（５）前記（４）で得た単位空間から、評価すべき多成分系薬剤のマハラノビス距離をＭＴ法を用いて得る工程
により得られたマハラノビス距離を用い、前記評価すべき多成分系薬剤の、基準群として選ばれた複数の多成分系薬剤からの相違の程度を判定することを特徴とする多成分系薬剤の評価方法。
【請求項２】
少なくとも、工程（１）ないし（８）
（１）評価すべき多成分系薬剤および基準群を構成する他の同種の多成分系薬剤の三次元高速液体クロマトグラフィーのフィンガープリントデータを得る工程
（２）上記（１）で得られたフィンガープリントデータから、等間隔ごとのデータを採取して圧縮したフィンガープリントデータを得る工程
（３）前記（２）の圧縮したフィンガープリントデータについて、多成分系薬剤の番号と、溶出時間又は検出波長の何れか一方とに対してＭＴ法における変数軸を割り振り、シグナル強度をＭＴ法における特徴量とする工程
（４）前記（３）の特徴量から、ＭＴ法を用いて単位空間を得る工程
（５）前記（４）で得た単位空間から、検出波長又は溶出時間ごとに、全ての多成分系薬剤のマハラノビス距離をＭＴ法を用いて得る工程
（６）多成分系薬剤の番号と、工程（３）で変数軸を割り振らなかった他方に対してＭＴ法における変数軸を割り振り、工程（５）で得られたマハラノビス距離をＭＴ法における特徴量とする工程
（７）前記（６）の特徴量から、ＭＴ法を用いて第２の単位空間を得る工程
（８）前記（７）で得た第２の単位空間から、評価すべき多成分系薬剤のマハラノビス距離をＭＴ法を用いて得る工程
により得られたマハラノビス距離を用い、前記評価すべき多成分系薬剤の、基準群として選ばれた複数の多成分系薬剤からの相違の程度を判定することを特徴とする多成分系薬剤の評価方法。
【請求項３】
工程（２）において、ＭＴ法における項目軸（ｘ軸）に溶出時間を割り振り、ＭＴ法における番号列軸（ｙ軸）に多成分系薬剤の番号を割り振り、特徴量をシグナル強度とする請求項１記載の多成分系薬剤の評価方法。
【請求項４】
工程（６）において、ＭＴ法における項目軸（ｘ軸）に検出波長を割り振り、ＭＴ法における番号列軸（ｙ軸）に多成分系薬剤の番号を割り振り、特徴量を、工程（５）で得られた各検出波長ごとのマハラノビス距離とする請求項２記載の多成分系薬剤の評価方法。
【請求項５】
工程（２）において、基準群が、ある一の会社が製造した多成分系薬剤の群である請求項１ないし請求項４の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項６】
工程（２）において、基準群が、フィンガープリントデータが得られている全ての多成分系薬剤の群である請求項１ないし請求項４の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項７】
データの圧縮を、多成分系薬剤の有効成分がカラムから溶出する時間範囲に限定し、溶出時間３秒〜３０秒から選ばれる一定の秒数ごとのデータポイントを選択採取することにより行う請求項１ないし請求項６の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項８】
データの圧縮を、溶出時間として、１００〜１０００個のデータポイントを選択採取することにより行う請求項１ないし請求項７の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項９】
データの圧縮を、２００ｎｍ〜３００ｎｍの範囲から選ばれた複数の検出波長におけるデータに限定することにより行う請求項１ないし請求項８の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項１０】
データの圧縮を、検出波長２ｎｍ〜２５ｎｍから選ばれる等間隔の波長ごとのデータポイントを選択採取することにより行う請求項１ないし請求項９の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項１１】
データの圧縮を、検出波長として、１０〜１００個のデータポイントを選択採取することにより行う請求項１ないし請求項１０の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項１２】
多成分系薬剤が、漢方製剤である請求項１ないし請求項１１の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項１３】
多成分系薬剤の三次元高速液体クロマトグラフィーのフィンガープリントデータを、多成分系薬剤の各化学成分ごとのシグナル強度に分離しないで行う請求項１ないし請求項１２の何れかの請求項記載の多成分系薬剤の評価方法。
【請求項１４】
請求項１ないし請求項１３の何れかの請求項記載の多成分系薬剤の評価方法を用い、一の多成分系薬剤の単位空間からのマハラノビス距離が、所定の値より大きい場合に、その一の多成分系薬剤を不合格にする多成分系薬剤の品質管理方法。
【請求項１５】
請求項１ないし請求項１３の何れかの請求項記載の多成分系薬剤の評価方法を用い、一の多成分系薬剤の単位空間からのマハラノビス距離が、所定の値以下となるように製造する多成分系薬剤の製造方法。
【請求項１６】
コンピューターに、少なくとも、手段（Ａ）ないし（Ｄ）
（Ａ）基準群を形成する複数の多成分系薬剤の三次元高速液体クロマトグラフィーから得られるフィンガープリントデータから、等間隔ごとのデータポイントを採取して圧縮したフィンガープリントデータを保存する手段
（Ｂ）評価すべき多成分系薬剤の三次元高速液体クロマトグラフィーから得られるフィンガープリントデータから、等間隔ごとのデータポイントを採取して圧縮したフィンガープリントデータを入力し、保存された上記（Ａ）の圧縮したフィンガープリントデータと組み合わせ一つのデータ群とする手段
（Ｃ）上記（Ｂ）のデータ群について、多成分系薬剤の番号と、溶出時間又は検出波長とに対してＭＴ法における変数軸を割り振り、シグナル強度をＭＴ法における特徴量として、ＭＴ法を用いて単位空間を得る手段
（Ｄ）ＭＴ法を用いて、評価すべき多成分系薬剤の、前記単位空間からのマハラノビス距離を得る手段
を実行させるための、一の多成分系薬剤の複数の多成分系薬剤からの相違の程度を一元的に判定するプログラム。
【請求項１７】
請求項１６記載のプログラムが記録されている情報記録媒体。
【請求項１８】
請求項１６記載のプログラムを伝送する情報伝送媒体。
【請求項１９】
（I）及び（II）
（I）コンピューターに、手段（Ａ）ないし（Ｄ）
（Ａ）基準群を形成する複数の多成分系薬剤の三次元高速液体クロマトグラフィーから得られるフィンガープリントデータから、等間隔ごとのデータポイントを採取して圧縮したフィンガープリントデータを保存する手段
（Ｂ）評価すべき多成分系薬剤の三次元高速液体クロマトグラフィーから得られるフィンガープリントデータから、等間隔ごとのデータポイントを採取して圧縮したフィンガープリントデータを入力し、保存された上記（Ａ）の圧縮したフィンガープリントデータと組み合わせ一つのデータ群とする手段
（Ｃ）上記（Ｂ）のデータ群について、多成分系薬剤の番号と、溶出時間又は検出波長とに対してＭＴ法における変数軸を割り振り、シグナル強度をＭＴ法における特徴量として、ＭＴ法を用いて単位空間を得る手段
（Ｄ）ＭＴ法を用いて、評価すべき多成分系薬剤の、前記単位空間からのマハラノビス距離を得る手段
を実行させるプログラム
（II）基準群を構成する複数の同種の多成分系薬剤の三次元高速液体クロマトグラフィーから得られるフィンガープリントデータ
が記録された情報記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は多成分系薬剤の評価方法に関し、更に詳しくは、多成分系薬剤を構成する各成分の相違の程度を、三次元高速液体クロマトグラフィー（以下、「３Ｄ−ＨＰＬＣ」と略記することがある）によって測定されたデータに、特定の解析手段を加えることによって得られた数値を用いる多成分系薬剤の評価方法に関するものである。
【背景技術】
【０００２】
多成分系薬剤、特に漢方製剤等の天然物由来の製剤の定量的、定性的プロフィールは、使用する原料生薬についての、地質学的要因、生態学的要因、収集時期、収集場所、収集年代、生育期の天候等が原因で変化する。そのため、漢方製剤等の多成分系薬剤については、その品質、安全性、効力等の評価のために一定の判定基準を規定し、その基準に基づいて、国家の監督機関、化学的組織、製造業者等が評価を行っている。【０００３】
しかしながら、多成分系薬剤の品質等の判定基準は、多成分系薬剤中のある特徴的な成分の一つないしは複数を適当に選択して、その含有量等に基づいて作成することが一般的であった。【０００４】
例えば、月刊薬事  vol.28,  No.3,  67-71  (1986)  には、多成分系薬剤において本質的成分の同定ができていない場合は、定量分析が可能、水に溶けやすい、熱水中で分解しない、他の成分と化学反応をしない等の物性を持った複数の成分を選択して、それら成分の含有量を化学分析し、得られた数値を評価の基準とすることが記載されている。【０００５】
一方、多成分系薬剤のクロマトグラフィーを測定し、各溶出時間ごとに紫外可視吸収スペクトル（以下、「フィンガープリントデータ」と略記する）を得て、その成分情報から評価の基準を作成することも知られている。例えば、特開２００２−２１４２１５号公報には、フィンガープリントデータの中からいくつかのピークを選択し、バーコード化することによって多成分系薬剤を評価することが記載されている。【０００６】
しかしながら、上記した方法には、「特定成分の含有量」又は「特定成分のクロマトピーク」という概念があって、その定量を行っているため、特定成分のクロマトピークをクロマトチャート上で、他の成分のクロマトピークから分離したり、特定成分のピーク面積や高さを求めるために、コンピューターによる波形処理操作が必要になり、このことがデータの正確さを低下させる一因になっていた。すなわち、ピークの溶出時間には厳密な意味での再現性がなく、そのためピークの波形処理にはバラツキ（dispersion）が生じていた。特に、このバラツキは、小さくブロードなピークや連続したピークを波形処理するときに著しく、評価方法として信頼性に欠けるものであった。更には、コンピューターによる波形処理は、一般に膨大な時間が必要であった。【０００７】
また上記した方法では、情報量（データポイント数）が特定成分のピーク数に限定されるため、情報量の増減が自由にはできず、データ処理時間の調整もできないため、最適な評価方法が確立できない場合があった。【０００８】
更に上記した方法では、複数の特定成分の含有量が複数の数値として得られるため、それらの数値を総合して判断する必要があり、多成分系薬剤を一の数値で評価することはできなかった。すなわち、評価すべき一の多成分系薬剤が、複数の多成分系薬剤群から、どの程度相違しているかの程度を、一の数値で表すことができなかった。【０００９】
そこで、バラツキが少なく信頼性が高く、データ処理時間も改善され、また、一の数値で判定することができる多成分系薬剤の評価方法が望まれていた。
【発明の開示】
【００１０】
本発明者らは、上記実情に鑑み鋭意検討を行った結果、３Ｄ−ＨＰＬＣのフィンガープリントデータを用い、このフィンガープリントデータの各項目を特定の方法で処理し、数理解析を行うことで、多成分系薬剤の品質、安全性、効力等を一の数値で評価できる方法を見出し、本発明に至った。
【００１１】
すなわち、本発明は、少なくとも、工程（１）ないし（５）
（１）評価すべき多成分系薬剤の３Ｄ−ＨＰＬＣのフィンガープリントデータを  
得る工程
（２）上記（１）で得られたフィンガープリントデータを、基準群を構成する他
の同種の多成分系薬剤の３Ｄ−ＨＰＬＣのフィンガープリントデータと組み  
合わせる工程
（３）前記（２）のフィンガープリントデータについて、多成分系薬剤の番号と、
溶出時間又は検出波長とに対してＭＴ法における変数軸を割り振り、シグナ
ル強度をＭＴ法における特徴量とする工程
（４）前記（３）の特徴量から、ＭＴ法を用いて単位空間を得る工程
（５）前記（４）で得た単位空間から、評価すべき多成分系薬剤のマハラノビス
距離をＭＴ法を用いて得る工程
により得られたマハラノビス距離を用い、前記評価すべき多成分系薬剤の、基準群として選ばれた複数の多成分系薬剤からの相違の程度を判定することを特徴とする多成分系薬剤の評価方法を提供するものである。
【００１２】
また本発明は、少なくとも、工程（１）ないし（８）
（１）評価すべき多成分系薬剤の３Ｄ−ＨＰＬＣのフィンガープリントデータを  
得る工程
（２）上記（１）で得られたフィンガープリントデータを、基準群を構成する他
の同種の多成分系薬剤の３Ｄ−ＨＰＬＣのフィンガープリントデータと組み  
合わせる工程
（３）前記（２）のフィンガープリントデータについて、多成分系薬剤の番号と、
溶出時間又は検出波長の何れか一方とに対してＭＴ法における変数軸を割り  
振り、シグナル強度をＭＴ法における特徴量とする工程
（４）前記（３）の特徴量から、ＭＴ法を用いて単位空間を得る工程
（５）前記（４）で得た単位空間から、検出波長又は溶出時間ごとに、全ての多  
成分系薬剤のマハラノビス距離をＭＴ法を用いて得る工程
（６）多成分系薬剤の番号と、工程（３）で変数軸を割り振らなかった他方に対    
してＭＴ法における変数軸を割り振り、工程（５）で得られたマハラノビス  
距離をＭＴ法における特徴量とする工程
（７）前記（６）の特徴量から、ＭＴ法を用いて第２の単位空間を得る工程
（８）前記（７）で得た第２の単位空間から、評価すべき多成分系薬剤のマハラ
ノビス距離をＭＴ法を用いて得る工程
により得られたマハラノビス距離を用い、前記評価すべき多成分系薬剤の、基準群として選ばれた複数の多成分系薬剤からの相違の程度を判定することを特徴とする多成分系薬剤の評価方法を提供するものである。
【００１３】
また更に本発明は、上記評価方法を実施するためのプログラム及び／又は単位空間を得るための基準群のフィンガープリントデータが記録された記録媒体を提供するものである。
【発明を実施するための最良の形態】
【００１４】
本発明は、フィンガープリントデータを、ＭＴ法を用いて数理解析して得られたマハラノビス距離によって、多成分系薬剤を評価する方法である。【００１５】
本発明方法には、工程（１）ないし（５）を含み、変数軸割り振りを１回行う方法（以下、「本発明１方法」という）と、工程（１）ないし（８）を含み、変数軸割り振りを２回行う方法（以下、「本発明２方法」という）がある。【００１６】
本発明１方法を実施するにはまず、工程（１）として、多成分系薬剤の３Ｄ−ＨＰＬＣのフィンガープリントデータを得る必要がある。【００１７】
ここで、多成分系薬剤とは、複数の有効化学成分を含有する薬剤と定義され、特に限定されるものではないが、好ましくは、植物抽出物、生薬や、これらを使用する漢方製剤等が例示される。また剤形も特に限定されず、例えば、湯剤、漢方エキス顆粒、漢方エキス液、漢方錠剤、漢方カプセル剤等が含まれる。
本発明において用いられる３Ｄ−ＨＰＬＣとは、所定の複数の溶出時間ごとに得られる溶出成分のスペクトルを意味する。【００１８】
この３Ｄ−ＨＰＬＣに使用されるカラムとしては、通常、高速液体クロマトグラフィーに用いられるものが使用される。【００１９】
また、３Ｄ−ＨＰＬＣでの検出波長は、特に限定はないが、好ましくは１５０ｎｍ〜９００ｎｍの範囲であり、特に好ましくは２００ｎｍ〜４００ｎｍの紫外可視吸領域、更に好ましくは２００ｎｍ〜３００ｎｍから選ばれた複数の波長である。【００２０】
スペクトルのシグナル強度については、透過率（transmittance）であっても、吸光度（absorbance）であってもよいが、吸光度が好ましい。【００２１】
工程（１）で得られるフィンガープリントとは、上記３Ｄ−ＨＰＬＣを用い、例えば、溶出時間と検出波長に対して、上記スペクトルのシグナル強度が３次元的に表示されたものであり、例えば図１（ａ）に示したようなものである。従ってフィンガープリントデータとは、少なくとも、多成分系薬剤の番号（ロット番号）、溶出時間、検出波長及びシグナル強度をデータとして有するものである。
３Ｄ−ＨＰＬＣを用いたフィンガープリントデータは、市販の装置によっても得ることができ、かかる市販の装置としては、島津製作所社製ＬＣ−ＶＰシステム等が挙げられる。【００２２】
本発明１方法における工程（２）は、上記（１）で得られたフィンガープリントデータを、基準群を構成する他の同種の多成分系薬剤の３Ｄ−ＨＰＬＣのフィンガープリントデータと組み合わせる工程である。【００２３】
基準群を構成する他の同種の多成分系薬剤の３Ｄ−ＨＰＬＣのフィンガープリントデータも上記した方法で得ることができる。測定条件等の差には特に限定はないが、評価すべき多成分系薬剤も基準群を構成する他の同種の多成分系薬剤も、同一カラムを用いて同条件で測定することが好ましい。【００２４】
基準群を構成する他の同種の多成分系薬剤とは、例えば漢方製剤の場合、同じ名称で呼ばれているものであっても、収集時期、収集場所、収集年代、生育期の天候等が異なるため、地質学的要因、生態学的要因、気候学的要因等によって、その成分組成に変動が生じるので、それらの要因が異なっていてもよい複数の多成分系薬剤をいう。【００２５】
基準群は、複数の多成分系薬剤から選ばれるが、ある一の会社が製造した多成分系薬剤の群であることが好ましい。また、フィンガープリントデータが得られている全ての多成分系薬剤の群であることも好ましい。【００２６】
更に、標準となる地質学的要因、生態学的要因、気候学的要因等が分かっている場合には、その要因の範囲に入っている多成分系薬剤を基準群として選択することも好ましい。例えば採取年を限定したり、採取場所を限定することもできる。また、一定期間内に特定の会社から販売された多成分系薬剤を選択することも好ましい。さらに特に一定の薬効があったり、基準としたい複数の多成分系薬剤がある場合には、それらを選択することも好ましい。また、複数の会社から同一名で販売等されている全ての多成分系薬剤を選択することも好ましい。また、評価したい一の多成分系薬剤自身も、基準群を構成するために用いられていてもよい。
基準群を構成する多成分系薬剤の数は特に限定はないが、５以上が好ましく、１００以上が特に好ましい。【００２７】
本発明１方法における工程（３）は、前記フィンガープリントデータのうち、多成分系薬剤の番号と、溶出時間又は検出波長の一方とに対して、ＭＴ法における変数軸を割り振り、シグナル強度をＭＴ法における特徴量とする工程である。【００２８】
ここで本発明のＭＴ法とは、現在品質工学でＭＴ法として一般に知られている計算手法を意味するものであるが、例えば、「品質工学の数理」日本規格協会発行（２０００）第１３６−１３８頁、品質工学応用講座「化学・薬学・生物学の技術開発」日本規格協会編（１９９９）第４５４−４５６頁及び品質工学，１１（５），７８−８４（２００３）に記載の方法が挙げられる。【００２９】
また一般に市販されているＭＴ法プログラムソフトも使用できる。市販のＭＴ法プログラムソフトとしては、（株）プローブ社のＰＲＡＴ  ｆｏｒ  Ｒｅｓｅａｒｃｈ  Ｖ１．０、ＰＲＡＴ  ｆｏｒ  Ｄｅｖｅｌｏｐｍｅｎｔ  Ｖ１．０；（財）日本規格協会のＴＭ−ＡＮＯＶＡ；（株）オーケンのＭＴＳ  Ｖｅｒ．２．０  ｆｏｒ  Ｅｘｃｅｌ、ＭＴ  ｆｏｒ  ｗｉｎｄｏｗｓ等が挙げられる。【００３０】
本発明１方法では、少なくとも、多成分系薬剤の番号と、溶出時間又は検出波長の何れか一方とに対してＭＴ法における変数軸を割り振り、シグナル強度をＭＴ法における特徴量とすることが必要である。
変数軸の割り振りには特に限定はないが、ＭＴ法におけるいわゆる項目軸（ｘ軸）に溶出時間を割り振り、いわゆる番号列軸（ｙ軸）に多成分系薬剤の番号を割り振り、ＭＴ法におけるいわゆる特徴量にシグナル強度を割り振ることが好ましい。【００３１】
ここで、上記項目軸（ｘ軸）と番号列軸（ｙ軸）は、以下のように定義される。すなわちＭＴ法においては、下記表１中のデータセットＸijについて、平均値mjと標準偏差σjを求め、Ｘijを規準化した値であるｘij＝（Ｘij−ｍj）／σjから、iとjの相関係数rを求めて、単位空間やマハラノビスの距離を得るが、この時、「平均値mjと標準偏差σjは、項目軸（ｘ軸）の値ごとに、番号列軸（ｙ軸）の値を変化させて求める」というように項目軸（ｘ軸）と番号列軸（ｙ軸）は定義される。【００３２】
【表１】
0004735539.tif 000002
【００３３】
本発明１方法における工程（４）は、工程（３）で軸が割り振られたデータと特徴量から、ＭＴ法を用いて、基準点と単位量（以下、「単位空間」と略記することがある）を得る工程である。ここで、基準点、単位量及び単位空間は、上記ＭＴ法の文献に記載に従い定義される。【００３４】
本発明１方法における工程（５）は、評価すべき多成分系薬剤の、前記単位空間からのマハラノビス距離（ＭＴ法においてＤ２と記載される）をＭＴ法を用いて得る工程である。ここでマハラノビス距離（Ｄ２）は、上記ＭＴ法の文献の説明と同様に定義され、またマハラノビス距離は、上記文献に記載の方法で求められる。【００３５】
変数軸として多成分系薬剤の番号と溶出時間を選択したときは、各検出波長ごとに、各多成分系薬剤のマハラノビス距離が得られ、変数軸として多成分系薬剤の番号と検出波長を選択したときは、各溶出時間ごとに各多成分系薬剤のマハラノビス距離が得られる。例えば、項目軸（ｘ軸）に溶出時間を割り振り、番号列軸（ｙ軸）に多成分系薬剤の番号を割り振り、ＭＴ法における特徴量にシグナル強度を割り振った場合には、各検出波長ごとにマハラノビス距離（Ｄ２）が得られる。【００３６】
このようにして得られたマハラノビス距離を用いて、評価すべき多成分系薬剤の、基準群として選ばれた複数の多成分系薬剤からの相違の程度を好適に判定することができる。【００３７】
更に本発明２方法は、一般に分割合成法、マルチＭＴ法又は多階ＭＴ法と呼ばれている方法の適用である。
本発明２方法を実施するには、工程（１）ないし工程（５）は、本発明１方法と同様に行われる。【００３８】
更に本発明２方法における工程（６）は、多成分系薬剤の番号と、工程（３）で変数軸を割り振らなかった他方に対してＭＴ法における変数軸を割り振り、マハラノビス距離をＭＴ法における特徴量とする工程である。【００３９】
例えば、工程（３）において、ＭＴ法におけるいわゆる項目軸（ｘ軸）に溶出時間を割り振り、いわゆる番号列軸（ｙ軸）に多成分系薬剤の番号を割り振り、ＭＴ法におけるいわゆる特徴量にシグナル強度を割り振った場合、各検出波長ごとに、各多成分系薬剤の番号に対するマハラノビス距離が得られるが、次いで新たに、項目軸（ｘ軸）に検出波長を割り振り、いわゆる番号列軸（ｙ軸）に多成分系薬剤の番号を割り振り、ＭＴ法におけるいわゆる特徴量に、上記工程（１）ないし工程（５）で得られたマハラノビス距離を割り振る方法が、好ましい方法として挙げられる。【００４０】
本発明２方法における工程（７）は、前記工程（４）と同様に、新たに軸を割り振られたデータ群に基づき、第２の単位空間を得る工程である。【００４１】
本発明２方法における工程（８）は、前記（７）で得た第２の単位空間から、評価すべき多成分系薬剤のマハラノビス距離をＭＴ法を用いて得る工程である。
工程（８）で得られたマハラノビス距離（以下「合成マハラノビス距離」ということがある）を用い、前記一の評価すべき多成分系薬剤の、基準群として選ばれた複数の多成分系薬剤からの相違の程度を好適に判定することができる。
合成マハラノビス距離を用いると、溶出時間、検出波長及び多成分系薬剤の番号のすべてを考慮した評価が可能になるので、３Ｄ−ＨＰＬＣの特長が生かせるので特に好ましい。【００４２】
本発明の好ましい態様としては、特定の会社から同一名称で製造又は販売された多成分系薬剤で基準群を構成させ、３Ｄ−ＨＰＬＣフィンガープリントデータを得ておき、別会社の販売した評価すべき一の多成分系薬剤の単位空間からのマハラノビス距離を求め、その別会社の多成分系薬剤が、基準群からどれだけ相違しているかの程度を判断する態様である。また、一定の基準群を選択し、その基準群から得られた単位空間を全国的に又は国際的に一定に定めておけば、漢方製剤等の多成分系薬剤を、全国的に又は国際的に、一律に客観的に評価できるので好ましい。【００４３】
本発明においては、前記フィンガープリントデータを、適当に圧縮することも好ましい。すなわち、フィンガープリントデータの情報量を任意に減量することも好ましい。本発明においては、シグナル強度、多成分系薬剤の番号（ロット番号）、溶出時間、検出波長等の各データ（測定点）の数を圧縮することが好ましい。特に溶出時間については、有効成分がカラムから溶出する時間範囲が限定されている場合があるので、その溶出時間帯のデータのみを選択する等により圧縮することが好ましい。検出波長についても有効成分の吸収スペクトルが限定された波長域にある場合があるので、その限定された波長域のみのデータを選択すること等によって圧縮することが好ましい。また、データの圧縮は、データの選択範囲を限定することのほか、多くのデータポイントから適当な数のデータポイントを、等間隔に抜き取る方法等によってもなされる。【００４４】
溶出時間については、有効成分がカラムから溶出し始める時間から溶出し終わる時間までのデータに限定することが好ましい。
具体的には、溶出時間については、有効成分がカラムから溶出する時間範囲のデータのみを選択したり、また通常の市販の３Ｄ−ＨＰＬＣ装置は、溶出時間約０．３秒〜１秒ごとにシグナル強度等のデータが得られるが、そのデータから、好ましくは、溶出時間約３秒〜３０秒から選ばれる等間隔の秒数ごとのデータだけを採取する方法、特に好ましくは５秒〜２０秒から選ばれる等間隔の秒数ごとのデータを採取する方法等が挙げられる。【００４５】
更に、通常得られる３Ｄ−ＨＰＬＣフィンガープリントデータのデータポイントの数を、１／４〜１／１００に減じることが好ましい。また、データポイントの数としては、１００〜１０００個とすることが好ましい。特に好ましくは、２００〜８００個である。【００４６】
検出波長については、重要な有効成分が特徴的な吸収を有する波長域を選択することが好ましい。例えば２００ｎｍ〜４００ｎｍのデータのみに限定することが好ましい。特に好ましくは２００ｎｍ〜３００ｎｍに限定することである。【００４７】
また、通常得られる３Ｄ−ＨＰＬＣフィンガープリントデータのデータポイントの数を、１／２〜１／５０に減じることが好ましい。更に、データポイントの数としては、１０〜１００個とすることが好ましい。特に好ましくは、２０〜８０個である。具体的には、通常、３Ｄ−ＨＰＬＣ装置で得られる波長１ｎｍオーダーごとのシグナル強度等のデータから、好ましくは波長２ｎｍ〜２５ｎｍから選ばれる一定の波長間隔ごと、特に好ましくは５ｎｍ〜２０ｎｍから選ばれる一定の波長間隔ごとのデータのみを採取する方法等が挙げられる。【００４８】
上記多成分系薬剤の評価方法は、品質管理に好ましく用いることができる。この場合、品質管理の対象となる一の多成分系薬剤の単位空間からのマハラノビス距離が、一定の範囲から選ばれたある値より大きい場合に、その一の多成分系薬剤を不合格にする品質管理方法が特に好ましい。また上記一定の範囲とは、多成分系薬剤の種類によって異なるが、例えば２〜１０００が好ましく、特に好ましくは１０〜１００の範囲が挙げられる。【００４９】
また、上記多成分系薬剤の評価方法は、多成分系薬剤の製造に好ましく用いることができる。この場合、単位空間からのマハラノビス距離が、一定の範囲から選ばれたある値以下となるように製造することが好ましい。また上記一定の範囲とは、多成分系薬剤の種類によって異なるが、例えば２〜１０００が好ましく、特に好ましくは１０〜１００の範囲が挙げられる。【００５０】
具体的には、評価すべき一の多成分系薬剤の単位空間からのマハラノビス距離がある値以下の場合のみ、その一の多成分系薬剤を製品として選択する多成分系薬剤の製造方法や、ある値以下となるように、複数の多成分系薬剤又はそれらの原料を混合する多成分系薬剤の製造方法等が挙げられる。【００５１】
本発明の評価方法は、その手段をプログラムとしてコンピューターに読み込ませて計算させることができる。
すなわち、コンピューターに、少なくとも、手段（Ａ）ないし（Ｄ）
（Ａ）基準群を形成する複数の多成分系薬剤の３Ｄ−ＨＰＬＣから得られるフィ  
ンガープリントデータを保存する手段
（Ｂ）評価すべき多成分系薬剤の３Ｄ−ＨＰＬＣから得られるフィンガープリン  
トデータを入力し、保存された上記（Ａ）のフィンガープリントデータと組
み合わせ一つのデータ群とする手段
（Ｃ）上記（Ｂ）のデータ群について、多成分系薬剤の番号と、溶出時間又は検  
出波長とに対してＭＴ法における変数軸を割り振り、シグナル強度をＭＴ法
における特徴量として、ＭＴ法を用いて単位空間を得る手段
（Ｄ）ＭＴ法を用いて、評価すべき多成分系薬剤の、前記単位空間からのマハラ  
ノビス距離を得る手段
を実行させるための、一の多成分系薬剤の複数の多成分系薬剤からの相違の程度を一元的に判定するプログラムが好適に使用できる。【００５２】
このプログラムが記録された情報記録媒体又は伝送される情報伝送媒体としても活用できる。ここで、情報記録媒体とは、汎用コンピュータが読み書き可能な情報格納手段（半導体メモリ、フレキシブルディスク、ハードディスク等）、光読出手段（ＣＤ−ＲＯＭ、ＤＶＤ等）等を言い、情報伝送媒体とは、プログラム情報を搬送波として伝搬させて供給するためのコンピュータネットワーク（ＬＡＮ、インターネット等のＷＡＮ、無線通信ネットワーク等）、システムにおける通信媒体（光ファイバや無線回線等）等を言う。【００５３】
更に、上記記録媒体には上記プログラムの他に、基準群のデータが記録されていてもよい。【００５４】
すなわち、（I）及び（II）
（I）上記手段（Ａ）及び手段（Ｂ）を実行させるプログラム
（II）基準群を構成する複数の同種の多成分系薬剤の３Ｄ−ＨＰＬＣから得られるフィンガープリントデータ
が記録された情報記録媒体である。
また、各多成分系薬剤ごとに、例えば各漢方製剤ごとに、上記（II）のフィンガープリントデータのみが記録された記録媒体も評価に有用である。【００５５】
本発明によると、ＨＰＬＣピークを波形処理しないので、データにバラツキが少なく信頼性が高く、情報量（データポイント数）が特定成分のピーク数に限定されないため、情報量の増減が自由にでき、また、複数の成分の含有量の数値を組み合わせて判定する必要がなく、一の数値で判定することができるので、簡単に評価すべき一の多成分系薬剤の基準群からの相違の程度を、正確に、客観的に、簡便に、一元的に判定できる。
【実施例】
【００５６】
以下に実施例及び試験例を挙げて本発明をさらに詳細に説明するが、本発明はこれらに何ら制約されるものではない。【００５７】
実  施  例  １
［工程（１）］及び［工程（２）］
株式会社ツムラ社が１９９８年から２００３年の期間に製造した桂枝茯苓丸漢方製剤（エキス顆粒）（以下、「ＴＪ−２５」と略記する）３４４ロット及びＡ社製の桂枝茯苓丸漢方製剤について、下記の方法で測定試料を調製し、下記の条件で３Ｄ−ＨＰＬＣフィンガープリントデータを取得した。【００５８】
＜測定試料の調製方法＞
各ロットのエキス顆粒３ｇとメタノール１００ｍＬを混合し、ホモジナイザーで１００秒間撹拌した。２分間放置後、液面と底面の中間位置から抽出液を採取し、プレフィルター付きメンブランフィルターで濾過して測定試料を調製した。【００５９】
＜３Ｄ−ＨＰＬＣフィンガープリントデータ取得条件＞
測定装置  ：ＬＣ−ＶＰシステム（島津製作所社製）
展開カラム：ＴＳＫ−ＧＥＬ  ８０ＴＳ（東ソー社製）
移動相    ：Ａ液  ５０ｍＭ  酢酸−酢酸アンモニウム  バッファー
Ｂ液  アセトニトリル
グラジエント条件：リニアグラジエント
カラム温度      ：４０℃
流速            ：１．０ｍＬ／ｍｉｎ
注入量          ：３０μＬ
検出器          ：ＳＰＤ−Ｍ１０Ａｖｐ（島津製作所社製）
測定波長        ：２００ｎｍ〜４００ｎｍ
解析ソフトウェア：ＣＬＡＳＳ−ＶＰ  (島津製作所社製)【００６０】
取得された３Ｄ−ＨＰＬＣフィンガープリントデータについて、表２に示すようにデータの圧縮を行った。
【００６１】
【表２】
0004735539.tif 000003
【００６２】
これにより、データポイント数は以下のようになった。
多成分系薬剤の番号・・・・３４４個
検出波長・・・・・・・・・・・・・・・１１個
溶出時間・・・・・・・・・・・・・・１０８個
データ圧縮前後のフィンガープリントの１例を図１に示す。これにより、問題なくＭＴ法で単位空間、マハラノビス距離が計算できるようになった。【００６３】
［工程（３）］
取得された３Ｄ−ＨＰＬＣフィンガープリントデータについて、表３の右欄に示すように、各検出波長ごとにＭＴ法における軸の割り振りを行った。
「各検出波長ごと」として、２００ｎｍ〜３００ｎｍの範囲を、１０ｎｍごとに等分割して得られる１１波長ごととした。
【００６４】
【表３】
0004735539.tif 000004
【００６５】
これによって、通常取得されるフィンガープリントデータが、ＭＴ法による数値解析に適した形式になった。【００６６】
［工程（４）］
「品質工学の数理」日本規格協会発行（２０００）第１３６−１３８頁、品質工学応用講座「化学・薬学・生物学の技術開発」日本規格協会編（１９９９）第４５４−４５６頁及び品質工学，１１（５），７８−８４（２００３）に記載されたように、ＭＴ法を用いて、単位空間、すなわち基準点と単位量を得た。
計算には、（株）オーケン社のＭＴ  ｆｏｒ  ｗｉｎｄｏｗｓ  を用いた。【００６７】
［工程（５）］
Ａ社製の桂枝茯苓丸漢方製剤の、工程（４）で得た基準点からのマハラノビス距離を、上記文献の記載に従ってＭＴ法により求めた。【００６８】
次いで工程（１）において、Ａ社製の桂枝茯苓丸漢方製剤に代えて、表４に記載した他の１２個を「評価すべき一の多成分系薬剤」としたほかは、上記工程（１）ないし工程（５）と同様にして、それぞれの桂枝茯苓丸漢方製剤のマハラノビス距離を計算した。
マハラノビス距離は、２００ｎｍないし３００ｎｍの範囲を１０ｎｍごとに等分割して得られる１１波長ごとに得られたが、一例として、２３０ｎｍにおけるマハラノビス距離を表４に示す。【００６９】
【表４】
0004735539.tif 000005
【００７０】
工程（１）ないし工程（５）で得られた各検出波長ごとのマハラノビス距離を用いることによって、ある一の桂枝茯苓丸の基準群からの相違の程度が判定できた。これによりこの評価方法が品質管理に有効であることが示された。【００７１】
実  施  例  ２
実施例１では、各検出波長ごとに、各多成分系薬剤の番号に対するマハラノビス距離が得られたが、実施例１で得られたマハラノビス距離をＭＴ法における特徴量とすることにより、以下のように多階ＭＴ法により合成マハラノビス距離を求めた。【００７２】
［工程（１）］ないし［工程（５）］
実施例１と同様に行った。【００７３】
［工程（６）］
表５に示すように、ＭＴ法における軸の割り振りを行った。【００７４】
【表５】
0004735539.tif 000006
【００７５】
すなわち、データセットしては、表６のようになる。【００７６】
【表６】
0004735539.tif 000007
【００７７】
表中、(Di,j)2は、多成分系薬剤の番号ｉ、検出波長ｊにおける、実施例１で得られたマハラノビス距離である。【００７８】
［工程（７）］及び［工程（８）］
実施例１と同様にして、単位空間を求め、合成マハラノビス距離を算出した。【００７９】
図２に、単位空間を得るために共通して使用したツムラ社製桂枝茯苓丸（ＴＪ−２５）の３４４ロットの合成マハラノビス距離を横軸に、そのマハラノビス距離を有する試料の個数を縦軸にして作成した図を示す。なお、単位空間を得るために使用したＴＪ−２５の３４４ロットのマハラノビス距離の平均値は、その定義から常に１である。【００８０】
評価すべき一の多成分系薬剤に該当する、ツムラ社製桂枝茯苓丸２個とツムラ社製以外の桂枝茯苓丸漢方製剤１１個の、基準点からの合成マハラノビス距離を表７に示す。
【００８１】
【表７】
0004735539.tif 000008
【００８２】
得られた合成マハラノビス距離を用いることによって、３４４個のツムラ社製桂枝茯苓丸漢方製剤（ＴＪ−２５）からの相違の程度が、一の数値で判定できた。そして例えば、Ｇ社製やＣ社製−IIは、３４４個のＴＪ−２５からなる基準群から最も相違していること等も明らかになった。
また、合成マハラノビス距離を用いることで、広い検出波長範囲による網羅的評価が可能であった。
さらにこの評価方法を用い、品質の安定した桂枝茯苓丸漢方製剤が製造可能であることが判った。
さらにこのツムラ社製桂枝茯苓丸漢方製剤（ＴＪ−２５）について、ＭＴ法の単位空間のデータが記録された記録媒体は、評価に有用であることが判った。【００８３】
実  施  例  ３
株式会社ツムラ社が１９９８年から２００３年の期間に製造した補中益気湯（以下、「ＴＪ−４１」と略記する）１００ロットについて、実施例１、２と同様の方法で測定試料を調製し、同様の条件で３Ｄ−ＨＰＬＣフィンガープリントデータを取得し、ＭＴ法を用いて、単位空間、すなわち基準点と単位量を得た。
上記単位空間を得るために使用した全ＴＪ−４１とツムラ社製以外の補中益気湯６試料の単位空間からのマハラノビス距離を求めた。マハラノビス距離を用いると、一の補中益気湯の、基準点からの相違の程度が簡単に正確に客観的に判断できることが判った。【００８４】
図３（ａ）にマハラノビス距離が１．０の補中益気湯（ＴＪ−４１）のフィンガープリントを示す。また、図３（ｂ）にマハラノビス距離が２７０４の補中益気湯のフィンガープリントを示す。フィンガープリントを見ただけでは、それらの相違の程度は不明であったが、本発明により得られたマハラノビス距離では、１．０と２７０４というように数値が大きく相違していることから、マハラノビス距離が、多成分系薬剤同士の相違を定量的に判断するために有力であることが判った。【００８５】
実  施  例  ４
［工程（１）］及び[工程（２）]
株式会社ツムラ社が１９９９年から２００３年の期間に製造した桃核承気湯漢方製剤（エキス顆粒）（以下、「ＴＪ−６１」と略記する）２５１ロットおよびＡ社製の桃核承気湯漢方製剤について、下記の方法で測定試料を調製し、下記の条件で３Ｄ−ＨＰＬＣフィンガープリントデータを取得した。【００８６】
＜測定試料の調製方法＞
各ロットのエキス顆粒２ｇとメタノール１００ｍＬを混合し、ホモジナイザーで１００秒撹拌した。２分間放置後、液面と底面の中間位置から抽出液を採取し、プレフィルター付きメンブランフィルターで濾過して測定試料を調製した。【００８７】
＜３Ｄ−ＨＰＬＣフィンガープリントデータ取得条件＞
測定装置  ：ＬＣ−ＶＰシステム（島津製作所製）
展開カラム：ＴＳＫ−ＧＥＬ  ８０ＴＳ（東ソー社製）
移動相    ：Ａ液  ５０ｍＭ  酢酸−酢酸アンモニウム  バッファー
Ｂ液  アセトニトリル
グラジエント条件：リニアグラジエント
カラム温度      ：４０℃
流速            ：１．０ｍＬ／ｍｉｎ
注入量          ：３０μＬ
検出器          ：ＳＰＤ−Ｍ１０Ａｖｐ（島津製作所製）
測定波長        ：２００ｎｍ〜４００ｎｍ
解析ソフトウェア：ＣＬＡＳＳ−ＶＰ（島津製作所製）【００８８】
取得された３Ｄ−ＨＰＬＣフィンガープリントデータについて、表８に示すようにデータの圧縮を行った。【００８９】
【表８】
0004735539.tif 000009
【００９０】
これにより、データポイント数は以下のようになった。
多成分系薬剤の番号・・・・・・・・２５１個
検出波長・・・・・・・・・・・・・・・・・・・２１個
溶出時間・・・・・・・・・・・・・・・・・１４５個【００９１】
[工程（３）]
取得された３Ｄ−ＨＰＬＣフィンガープリントデータについて、表９の右欄に示すように、各検出波長ごとにＭＴ法における軸の割り振りを行った。
「各検出波長ごと」として、２００ｎｍ〜４００ｎｍの範囲を、１０ｎｍごとに等分割して得られる２１波長ごととした。【００９２】
【表９】
0004735539.tif 000010
【００９３】
これによって、通常取得されるフィンガープリントデータが、ＭＴ法による数値解析に適した形式になった。【００９４】
[工程（４）]
ＭＴ法を用いて、単位空間、すなわち基準点と単位量を得た。
計算には（株）オーケン社の  ＭＴ ｆｏｒ ｗｉｎｄｏｗｓを用いた。【００９５】
[工程（５）]
Ａ社製の桃核承気湯漢方製剤の、工程（４）で得た基準点からのマハラノビス距離を、ＭＴ法により求めた。【００９６】
次いで工程（１）において、Ａ社製の桃核承気湯漢方製剤に代えて、表１０に記載した他の９個を「評価すべき一の多成分系薬剤」としたほかは、上記工程（１）ないし工程（５）と同様にして、それぞれの桃核承気湯漢方製剤のマハラノビス距離を計算した。
マハラノビス距離は、２００ｎｍないし４００ｎｍの範囲を１０ｎｍごとに等分割して得られる２１波長ごとに得られたが、一例として、２８０ｎｍにおけるマハラノビス距離を表１０に示す。【００９７】
【表１０】
0004735539.tif 000011
【００９８】
工程（１）ないし工程（５）で得られた各検出波長ごとのマハラノビス距離を用いることによって、ある一の桃核承気湯の基準群からの相違の程度が判定できた。これによりこの評価方法が品質管理に有効であることが示された。【００９９】
［工程（６）］
各検出波長ごとに、各多成分系薬剤の番号に対するマハラノビス距離が得られたが、工程（１）ないし工程（５）で得られたマハラノビス距離をＭＴ法における特徴量とすることにより、以下のように多階ＭＴ法により合成マハラノビス距離を求めた。
表１１に示すように、ＭＴ法における軸の割り振りを行った。【０１００】
【表１１】
0004735539.tif 000012
【０１０１】
すなわち、データセットとしては、表１２のようになる。【０１０２】
【表１２】
0004735539.tif 000013
【０１０３】
表中、(Di,j)2は、多成分系薬剤の番号ｉ、検出波長ｊにおける、工程（５）で得られたマハラノビス距離である。【０１０４】
［工程（７）］および［工程（８）］
工程（１）ないし工程（５）と同様にして、単位空間を求め、合成マハラノビス距離を算出した。【０１０５】
図４に、単位空間を得るために共通して使用したツムラ社桃核承気湯（ＴＪ−６１）の２５１ロットの合成マハラノビス距離を横軸に、そのマハラノビス距離を有する試料の個数を縦軸にして作成した図を示す。なお、単位空間を得るために使用したＴＪ−６１の２５１ロットのマハラノビス距離の平均値は、その定義から常に１である。
評価すべき一の多成分系薬剤に該当する、ツムラ社製桃核承気湯３個とツムラ社製以外の桃核承気湯漢方製剤７個の、基準点からの合成マハラノビス距離を表１３に示す。【０１０６】
【表１３】
0004735539.tif 000014
【０１０７】
得られた合成マハラノビス距離を用いることによって、２５１個のツムラ社製桃核承気湯漢方製剤（ＴＪ−６１）からの相違の程度が、一の数値で判定できた。そして例えば、Ｆ社製−IやＤ社製は、２５２個のＴＪ−６１からなる基準群から最も相違していること等も明らかとなった。【０１０８】
図５（ａ）にマハラノビス距離が１．０の桃核承気湯（ＴＪ−６１）のフィンガープリントを示す。また、図５（ｂ）にマハラノビス距離が３５１７８３３０の桃核承気湯のフィンガープリントを示す。
【産業上の利用可能性】
【０１０９】
本発明によれば、信頼性が高く、情報量の増減が自由にできるためデータ処理時間を適当に短縮でき、一の数値で簡単に評価すべき一の多成分系薬剤の基準群からの相違を判定できる。従って、多成分系薬剤の相違の程度を容易に評価できるため、社内での品質管理に使用できるものである。また、これを用いて好適に多成分系薬剤が製造できるものである。
また更には、全国的又は国際的統一基準を本発明の単位空間で定義しておくことによって、ある多成分系薬剤のその基準からの相違を判定できるものであり、安定品質の漢方製剤等の多成分系薬剤を供給することができる。
【図面の簡単な説明】
【０１１０】
【図１】データ圧縮前後の桂枝茯苓丸漢方製剤フィンガープリントを示す図である。図中、（ａ）はデータ圧縮前、（ｂ）はデータ圧縮後のものを示す。
【図２】３４４個の桂枝茯苓丸漢方製剤（ツムラ社製ＴＪ−２５）のマハラノビス距離と個数との関係を示す図である。
【図３】マハラノビス距離が１．０の補中益気湯（ツムラ社製ＴＪ−４１）とマハラノビス距離が２７０４の補中益気湯のフィンガープリントを示す図である。図中、（ａ）はマハラノビス距離が１．０、（ｂ）はマハラノビス距離が２７０４のものを示す。
【図４】２５１個の桃核承気湯漢方製剤（ツムラ社製ＴＪ−６１）のマハラノビス距離と個数との関係を示す図である。
【図５】マハラノビス距離が１．０の桃核承気湯（ツムラ社製ＴＪ−６１）とマハラノビス距離が３５１７８３３０の桃核承気湯のフィンガープリントを示す図である。
【図１】
0004735539.tif 000015
【図２】
0004735539.tif 000016
【図３】
0004735539.tif 000017
【図４】
0004735539.tif 000018
【図５】
0004735539.tif 000019