(57)【特許請求の範囲】
【請求項１】
２つ以上のディテクタｉ（ｉは１〜ｎ、ｎはディテクタの数）を有するシステムにおいて、前記ディテクタ間のバンド拡大現象を補正するために用いられる拡大モデルの最適フィットパラメータを定めるための方法であって、
ａ）調節可能なパラメータαｊ（ｊは０，・・・，ｌ、ここでｌは調節可能な数）のセットを含む、拡大モデルＢ（α０，α１，・・・，αｌ，τ）を選択するステップ、
ｂ）参照用のサンプルを注入するステップ、
ｃ）各前記ディテクタｉからの単一成分のピークに対応する信号Ｄｉ（ｔ）、ここでピークは単一成分が溶出する時間の範囲である、を時間の関数として集めるステップ、
ｄ）参照用の信号Ｄｎ（ｔ）を時間的に最も拡大された信号として規定するステップ、
ｅ）各ディテクタｉについて、拡大モデルを用いた以下の式のχｉ２モデルを形成し、前記χi２モデルは、参照用の信号Ｄｎ（ｔ）を前記ディテクタｉと対比して、前記ピークにわたって最小化され、
【数１】
0004813784.tif 000035
ここで前記最適フィットパラメータはβｉ、αｉｊおよびτｉであり、
（ａ）βｉはディテクタｉのスケールファクタであり、
（ｂ）αｉｊはディテクタｉとパラメータｊについての拡大の程度を表し、そして、
ｆ）前記χｉ２モデルを最小化して前記ディテクタのピークに対する最適フィットパラメータを算出するステップ、を含む方法。
【請求項２】
前記χｉ２モデルの最小化は、非線形最小二乗アルゴリズムを使用することによって達成される、請求項１に記載の方法。
【請求項３】
前記非線形最小二乗アルゴリズムはマルカルトによって開発された種類のものである、請求項２に記載の方法。
【請求項４】
前記ディテクタ間のバンド拡大は希釈によって引き起こされる、請求項１に記載の方法。
【請求項５】
前記ディテクタ間のバンド拡大は混合によって引き起こされる、請求項１に記載の方法。
【請求項６】
前記ディテクタ間のバンド拡大は、ディテクタのセルおよび／またはその接続部内の機械的不良によって引き起こされる、請求項１に記載の方法。
【請求項７】
前記ディテクタ間のバンド拡大は内部計器効果によって引き起こされる、請求項１に記載の方法。
【請求項８】
前記内部計器効果は電子フィルタ処理によって引き起こされる、請求項７に記載の方法。
【請求項９】
前記内部計器効果は各ディテクタによって測定されるサンプル体積の違いによって引き起こされる、請求項７に記載の方法。
【請求項１０】
単一成分の前記ピークは、単分散画分に対応し、クロマトグラフィ分離により他の画分から分離される、請求項１に記載の方法。
【請求項１１】
単一成分の前記ピークは、それぞれ画分に対応し、前記画分は分離装置によって分離されず、感知できる程度の分離が生じない、請求項１に記載の方法。
【請求項１２】
前記拡大モデルは以下の式によって与えられ、
【数２】
0004813784.tif 000036
ここで、α１とα２は前記ディテクタ間のバンド拡大現象を特徴づけ、α１は内部装置平均によるガウス拡大を特徴づけ、α２はディテクタ間混合体積を特徴づけ、ｔ≧τの場合Ｕ（ｔ−τ）＝１であり、ｔ＜τのときＵ（ｔ−τ）＝０である、請求項１に記載の方法。
【請求項１３】
一連のディテクタを連続的に通過するサンプルの選択された物理特性を引出すための方法であって、前記ディテクタの一部がその信号にバンド拡大を示す場合に、前記ディテクタを通る前記サンプルに応答して前記ディテクタによってもたらされる信号の組合せを使用し、前記方法は
ａ）ディテクタｉについて請求項１に記載の方法により最適フィットパラメータαｉｊおよびτｉを求めるステップ、
ｂ）対応する一連のディテクタ信号を引出すために、前記ディテクタに以下の式の拡大関数を与えるステップを含み、
【数３】
0004813784.tif 000037
ここで、Ｄｉはディテクタ信号であり、Ｂはステップａ）で用いた拡大モデルであり、Ｄiｂ（ｔ）は拡大されたディテクタ信号であり、前記ディテクタ信号はすべて類似の拡大を有し、さらに
ｃ）前記拡大関数を与えた後、拡大された前記ディテクタ信号を用いて前記測定されたサンプルの前記選択された物理特性を引出すステップを含む、方法。
【請求項１４】
一連のディテクタを連続的に通過するサンプルの選択された物理特性を引出すための方法であって、前記ディテクタの一部がその信号にバンド拡大を示す場合に、前記ディテクタを通る前記サンプルに応答して前記ディテクタによってもたらされる信号の組合せを使用し、前記方法は
ａ）ディテクタｉについて請求項１２に記載の方法により最適フィットパラメータα１およびα２を求めるステップ、
ｂ）対応する一連のディテクタ信号を引出すために、前記ディテクタに以下の拡大関数を与えるステップを含み、
【数４】
0004813784.tif 000038
ここで、Ｄｉはディテクタ信号であり、Ｂはステップａ）で用いた拡大モデルであり、Ｄiｂ（ｔ）は拡大されたディテクタ信号であり、前記ディテクタ信号はすべて類似の拡大を有し、さらに
ｃ）前記拡大関数を与えた後、拡大された前記ディテクタ信号を用いて前記測定されたサンプルの前記選択された物理特性を引出すステップを含む、方法。
【請求項１５】
Ｒ（θ）＝Ｋ*ＭｗｃＰ（θ）［１−２Ａ2ＭwｃＰ（θ）］＋Ｏ（ｃ3）の関係から定められる前記選択された物理特性は、濃度信号ｃ（ｔ）から引出された、前記サンプルの重み平均モル質量Ｍwおよび二乗平均半径ｒg、ならびに光散乱ディテクタＤi（ｔ）からなる一連のディテクタからのｉ光散乱信号およびｄＲＩディテクタから引出された過剰レイリー比Ｒ（θ，ｔ）であり、前記ｄＲＩディテクタは前記光散乱ディテクタ信号に相関して拡大を示す濃度信号をもたらし、前記光散乱ディテクタ信号は拡大されている、請求項１３または１４に記載の方法。
【請求項１６】
前記ディテクタ信号はＵＶディテクタから出力されて、続いて多角度光散乱ディテクタから出力されており、前記多角度光散乱信号は拡大されている、請求項１３または１４に記載の方法。
【請求項１７】
前記ディテクタ信号は屈折率ディテクタから出力されて、続いて粘度計ディテクタから出力されており、前記屈折率ディテクタ信号は拡大されている、請求項１３または１４に記載の方法。
【請求項１８】
前記ディテクタｉはクロマトグラフィのディテクタである、請求項１に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はディテクタ間バンド拡大の現象を補正するための方法に関する。
【背景技術】
【０００２】
液体クロマトグラフィ技術による溶液における高分子種の分析は一般に適切な溶媒内にサンプルを用意し、そのアリコートをクロマトグラフに注入することによって行なわれる。クロマトグラフは通過するサンプルを分化するさまざまな種類の分化装置を含む。このような手段によって一般に大きさ、塊、またはカラム親和力に基づいて分けられると、サンプルは光散乱、屈折率、紫外線（ＵＶ）、粘性応答などによって分析を受ける。たとえば分離がサイズ排除カラムによって行なわれる、特定のサンプルの塊およびサイズの分布を定めるために、クロマトグラフは多角度光散乱ＭＡＬＳによって、その後で示差屈折計ｄＲＩによって、サンプルを順次測定する。ｄＲＩは濃度を測定するのに対して、ＭＡＬＳユニットは各溶出画分の角度を関数とした過剰レイリー比を測定する。入射光の波長よりもはるかに小さい分子に対しては、単一角度での光散乱測定で十分な場合が多い。
【０００３】
各画分が検出装置を通ると、「バンド」または「ピーク」と呼ばれる信号が出力される。分散および混合効果により、これらのバンドはサンプルが異なる装置を通るたびに幾分広げられる。単分散タンパク質の低い濃度のアリコートからなるサンプルを考えてみる。この場合、レイリー過剰比はモル質量および濃度に正比例する。光散乱信号および濃度信号は同一の形になるべきであり、これら２つの反応が同じ領域を有するよう正規化されたのなら完全に重なる。しかし、サンプルがＭＡＬＳディテクタを通ってｄＲＩディテクタに入ると、サンプルの分散および混合に寄与する中間領域および接続部を通過することになる。バンド拡大の多くの原因およびその救済策の一部は１９７９年にJohn Wiley & Sonsによって出版された、Yau他の書籍「近代のサイズ排除液体クロマトグラフィ」において詳細に説明されている。上述の例では、ｄＲＩ信号はＭＡＬＳ信号に対して幾分広げられて現われる。
【０００４】
拡大の現象は図１に示され、拡大現象に対して補正がなされていない、計算されたモル質量とともに、水性バッファにおけるウシ血清アルブミンＢＳＡのサンプルのクロマトグラムが示される。補正されていないモル質量はトレース１として示される。これは単量体Ｍw0のモル質量によって除算された形で示されるので、縦軸は単量体には１となり、二量体には２、等となる。トレース２は時間を関数とした、９０°の光散乱信号に対応する。トレース３はｄＲＩ信号に対応する。１６分から１８分にある最も右のピーク領域は純粋単量体からのデータであり、１４分から１６分にあるピーク領域は二量体からのデータであるよう、凝集状態は細分化されている。各ピークにおいて、拡大によりピーク中央の濃度は抑制され、「翼」部分の濃度は高められる。各溶出時においてモル質量を定めるためにこの拡大されたｄＲＩ信号をＭＡＬＳ信号と組み合わせると、ピークの中央近くで定められるモル質量は体系的に高すぎる推定値となり、翼部分のものは低すぎる推定値となる。濃度ディテクタがＭＡＬＳディテクタの上流にある場合、この状態は逆となる。これにより、各ピーク領域内のモル質量データは不定なものとなる。バンド拡大がなければ、データは単量体に対しては１であり、二量体に対しては２である一連の平坦域からなる。細分化されたサンプルがより多分散されると、バンド拡大の現象は目立たなくなる（しかしまだ存在する）。個々の凝集状態からなる十分に細分化されたサンプルでは、この問題は明らかである。図１において、単量体のピークは拡大問題が最も明らかであるよう分離している。凝集体は累進的に分離されなくなるので、問題はそれほど顕著とはならない。
【０００５】
このような歪みを補償するために、さまざまな方法が用いられてきた。その大部分は１９６９年の文献に提示されているL.H. Tungの研究に基づいており、出版されたJournal of Applied Polymer Scienceの第１３巻、頁７７５以降にある。前述のYau他の書籍は、バンド拡大を補正するための手段の一部を記載している。クロマトグラフィの文献のうち、多くはこの補正を行なうための手段を扱っている、すなわち拡大源がなかったのならとっていたであろう形に拡大されたバンドを修復する内容を扱っている。これらの技術のほとんどは数値的に不安定であり、負の濃度または負の分散強度などのような物理的に不合理な結果をもたらすので、めったに用いられない。ＭＡＬＳ測定によるタンパク質のキャラクタライゼーションは特に重要な分野であるが、ここでもバンド拡大の補正はめったに見られない。なぜなら、それを行なうのは非常に困難だからである。
【０００６】
バンド拡大現象は液体クロマトグラフィのさまざまなマルチディテクタ運用で起こる。たとえば、サンプルの固有粘度が決定できるようオンライン粘度計でサンプルが測定されるのなら、ｄＲＩも必要である。サンプルがｄＲＩから粘度計に到達するまでに、バンド拡大を受けて、その結果ｄＲＩ測定から引出された濃度曲線に対して比粘度曲線が広がることとなる。ここでも、拡大された信号を拡大がなかったのなら取っていたであろう形になるよう復元することに伴う演算上の困難により、サンプルが非凝集タンパク質のような単量体である場合など特に問題のある結果をもたらす。
【非特許文献１】１９７９年、Yau他著、John Wiley & Sons出版の書籍「近代のサイズ排除液体クロマトグラフィ」
【非特許文献２】１９６９年Journal of Applied Polymer Scienceの第１３巻、頁７７５以降、L.H. Tung著
【発明の開示】
【発明が解決しようとする課題】
【０００７】
本発明の目的は分析方法を提供することであり、この分析方法によってバンド拡大現象はソフトウェアで容易にかつ正確に補正できる。本発明の別の目的は、従来のアプローチとの違いが明確に区別できるような補正に対する新しいアプローチを提供することである。本発明はタンパク質サンプルの分離分析に対して最も有用であると思われる。これを適用することにより、タンパク質共役および凝集状態のキャラクタライゼーションがより正確なものとなる。本発明のさらなる応用は、サンプルが濃度ディテクタと粘度計との間を通過することに伴うバンド拡大現象を補正することにより、固有粘度測定値の判定を向上させることにある。本発明の重要な応用は、たとえばTrainoffおよびWyattによる、米国特許第６，４１１，３８３号および現在発行される段階にある、２００２年７月２４日に出願された特許出願１０／２０５，６３７号に記載されている、非細分化サンプルの第２のビリアル係数を迅速に定めることへの適用である。この方法をうまく実施することは、各溶出量における濃度の二乗を正確に定める機能に依存する。これには、各溶出量における濃度の正確な測定が必要である。後者の決定は一般に非細分化サンプルを使って行なわれ、そのディテクタ応答は集められた各インターバルにおいて質量組成が同じである単一ピークを示す。この単一のピークは一連のディテクタを通ることにより拡大し、その拡大現象を補正しない限り、引出された結果には体系的誤りを含むこととなる。
【課題を解決するための手段】
【０００８】
本発明は液体クロマトグラフのディテクタ間で起こるバンド拡張現象を補正するための方法を含む。バンドの形を第１のディテクタのものに復元することによって第２のディテクタでのバンド拡大を補正するのではなく、本発明は第１のディテクタのバンドを第２のディテクタでもたらされた拡大バンドと適合するよう第１のディテクタのバンドを広げる。これはすべてのディテクタでのバンド拡大現象を補正するための、より簡単で安定した、かつ実現が容易である手段をもたらす。ピークが部分的に重畳する場合には、この方法はピーク分解能の低下をもたらし得るが、本発明の方法は、十分に分けられかつ分離した
ピークを備えた単分散サンプルに対して主な適用範囲を有すると考えられる。本発明は一連の接続されたディテクタに適用可能であり、そのディテクタの数は任意である。
【発明を実施するための最良の形態】
【０００９】
溶液における分子は一般にその重み平均モル質量Ｍw、その平均平方半径＜ｒg2＞および第２のビリアル係数Ａ2によって特徴付けられる。後者は、分子と溶媒との間の相互作用の測定値である。細分化されない溶液について、その特性は、１９４８年度のブルーノ・ジム(Bruno Zimm)の研究書類であって、Journal of Chemical Physicsの第１６巻、頁１０９３−１０９９に記載される方法を用いて、光を散乱させる態様を測定することによって、定めることができる。より最近では、２００２年７月２４日に出願され、現在発行される段階にある、TrainoffおよびWyattの同時係属中の特許出願１０／２０５，６３７号で記載される方法は、従来のジムアプローチに取って代わる進歩した技術を示す。
【００１０】
少量の溶液からの光散乱は、広い範囲の角度および濃度にわたって測定される。光散乱測定値から引出された特性は、ジムが作成した以下の式によって関連付けられる：
【００１１】
【数１】
0004813784.tif 000002
【００１２】
ここでＲ（θ）はＲ（θ）＝［Ｉs（θ）−Ｉsolv（θ）］ｒ2／ＩOＶとして定義される、単位立体角当りの方向θにおける測定された過剰レイリー比であり、Ｉs（θ）は単位立体角当りの溶液によって散乱した光の強度、Ｉsolv（θ）は単位立体角当りの溶媒から散乱した光の強度、ＩOは入射強度、ｒは散乱体積からディテクタまでの距離、Ｖはディテクタで読取られる体積、Ｐ（θ）は以下の式によって定義される分散分子の形状因子であり、
【００１３】
【数２】
0004813784.tif 000003
【００１４】
上記の式において、Ｎaはアボガドロ数であり、ｄｎ／ｄｃは屈折率増分であり、ｎ0は溶媒の屈折率であり、ｃは濃度であり、λ0は真空における入射光の波長である。入射コリメート光ビームは、光散乱ディテクタを含む面に対して垂直偏波されると想定される。形状因子は以下の式によって平均平方半径＜ｒg2＞に関連付けられる：
【００１５】
【数３】
0004813784.tif 000004
【００１６】
多角度光散乱計器は、角度および時間双方を関数とした、流体サンプルからの光散乱量を測定する。入射光の波長に対して非常に小さい分子では、光散乱測定はたとえば９０°のような単一の散乱角に制限することができる。なぜなら、角度に対して散乱した強度に
ついて測定可能な変動はないからである。しかし、この制限において、モル質量は測定できるが、平均平方半径は測定できない。
【００１７】
サイズ排除クロマトグラフィＳＥＣなどのように、分離系の溶離液を分析するのに計器を用いる場合、サンプル組成および濃度は時間とともに変化する。さらに、ＳＥＣカラムは典型的にサンプル濃度を薄めるので式(１)の右辺の第２および第３の項が第１の項と比べて無視できる。形式的には、これは以下の場合当てはまる：
【００１８】
【数４】
0004813784.tif 000005
【００１９】
典型的に、光散乱計測と組合せられるクロマトグラフィ分離では、Ａ2が前の実験から既知でない限り、上記の条件が想定される。Ａ2が既知の場合、値を直接式(１)に使用することができる。さらに、分離が正しく行なわれるのなら、各時点において、サンプルは本質的に大きさにおいて、さらに一般的に質量において、単分散する。サンプルが質量においても単分散するのなら、以下のように簡素化できる：
【００２０】
【数５】
0004813784.tif 000006
【００２１】
したがって、光散乱信号は濃度信号と正比例し、モル質量は以下の比率によって計算することができる：
【００２２】
【数６】
0004813784.tif 000007
【００２３】
ＭＡＬＳシステムの補正操作をテストするには、分離カラムにほぼ単分散であるサンプルを注入して、単分散ピークについて光散乱信号および濃度信号の両方を測定する。正規化されているのなら、これらの信号は式(５)によって表わされる定常性に対応してどの溶出時間でも完全に重なるはずである。しかし、濃度ディテクタがＭＡＬＳディテクタの下流にあるのなら、２つの現象、すなわち混合および分散により、ピークの形が変わる。この場合、ピークは広げられて現われ、ピークの中央は下げられ、「翼」部分は高められる。したがって、式（５）の比率が作成されると、それは一定ではなく、ピーク近くで引出されたモル質量は体系的に高すぎる推定値となり、翼近くのものは低すぎる値となる。これは、濃度ディテクタがＭＡＬＳディテクタの上流にある場合、逆となる。
【００２４】
細分化せずにクロマトグラフに注入されるサンプルは多くの場合フロー注入と呼ばれるが、あるディテクタから他のディテクタに通過する際にバンドが拡大される。各連続する溶出量でのモル質量分布は一定のままであると予測されるが、バンド拡大現象は選択されたディテクタでの下流の拡大から引出された種々の重み平均特性に影響する。したがって、多分散モル質量分布からなる注入サンプルピークは、たとえば光散乱ディテクタおよびｄＲＩ濃度ディテクタ間を通る場合に広がる。細分化されない光散乱信号が拡大された濃度検出信号と組合せられると、ピークの中央近くの引出された重み平均モル質量はより大
きい値として計算され、翼での対応する値はより小さいものとして計算される。この結果は細分化されたサンプルで見られるものと全く等価である。
【００２５】
サンプルの混合は主にシステムのジオメトリに依存する、すなわち管の長さおよび直径、継ぎ手の数、表面粗さ、および機械的不良によってもたらされる含有物の数、ならびに測定セルのジオメトリに依存する。混合の効果は、粘度のような総合流体特性がサンプルによって実質的に変わらない範囲において、サンプル組成と独立している。それに対して、分散はサンプルの特性に明らかに依存する。
【００２６】
２つのディテクタ間に流れるサンプルに対する分散および混合の相対的重要性を見積もることは有益である。一定の球体については、以下に示されるストークス-アインシュタイン関係は、分散定数ＤTを球体の半径ｒに関連付け、ここでｋBはボルツマン定数であり、Ｔは絶対温度であり、ηは溶媒粘度である。
【００２７】
【数７】
0004813784.tif 000008
【００２８】
混合の源はたくさんあるが、スケールを推定するために、十分に混合された槽の単純なモデルを取上げる。光散乱フローセルは体積Ｖの混合チャンバであると仮定する。流量はｆであり、チャンバに入るサンプルの濃度はｃi（ｔ）であると仮定する。フローセルを出るサンプルの濃度は以下の式によって与えられる：
【００２９】
【数８】
0004813784.tif 000009
【００３０】
ここでｔｆ＝Ｖ／ｆは槽を満たすのにかかる時間である。狭いパルスであって、ｃi（ｔ）＝ｍOδ（ｔ）であるよう、フローセルに注入される質量ｍOのサンプルを含む場合を考えてみる。セルから出るサンプルの濃度は以下の式によって表わされる：
【００３１】
【数９】
0004813784.tif 000010
【００３２】
ここで、ｔ＜０ではθ（ｔ）＝０であり、ｔ≧０ではθ（ｔ）＝１である。本文中において、フローセルの体積はＶ＝８０μｌであり、ｆ＝１．０ｍｌ／分の流量で送られると仮定する。拡大についての指数時間定数はｔｆ＝４．８秒である。この混合チャンバが内径０．２５ｍｍの毛細管によって流入流出されるのなら、最初の局所濃度パルスは長さ１．６ｍの管にわたり不鮮明となる。当然、サンプルが光散乱セルを流れる場合、「十分混合」されないが、これは混合の長さスケールでのよい推定値である。
【００３３】
次に、小さな分子の分散の効果を考えてみる。分子は平均で以下の距離で分散する：
【００３４】
【数１０】
0004813784.tif 000011
【００３５】
ここでＤTは分子の並進拡散定数であり、ｔは時間である。再度、これを本文中に置くために、７．７×１０-7ｃｍ2／秒の分散定数を有するウシ血清アルブミンタンパク質ＢＳＡを考えてみる。濃度スパイクは４．８秒でどれくらい遠く分散するであろうか。式（９）はたった２７μｍの距離を管に沿って分散することを示す。したがって、主要な混合機構はもともとは非分散であるとして、それゆえ分散定数に依存するのではなく、したがって分子の大きさに依存するものではないと考えられた。ディテクタ間拡大の主要なメカニズムは混合によるものであり、サンプルとは無関係である。したがって、システムのジオメトリでの拡大が基準サンプルで特徴付けることができるのなら、もたらされる拡大パラメータを用いて以降のすべてのデータ実行を補正することができる。
【００３６】
ディテクタ間拡大用の単純な線形モデルを考えてみる。説明のため、濃度ディテクタはＭＡＬＳディテクタの下流にあると仮定する。さらに、式（３）は満たされ、単分散基準サンプルについてのデータが集められたものとする。したがって、測定された濃度は、パラメータ化された関数であるＢ（α0，α1，…，αn，τ）でのたたみ込みを介して、ＭＡＬＳディテクタの内在する濃度に依存する。以下の式が与えられる：
【００３７】
【数１１】
0004813784.tif 000012
【００３８】
ここでｃm（ｔ）は拡大された後の下流ディテクタでの測定された濃度であり、ｃ（ｔ）は光散乱ディテクタを通ったときのサンプルの濃度である。パラメータαnはモデルに特有であり、拡大関数の幅およびディテクタ間の遅延量を含む。拡大補正の目的は、これらのパラメータの最適値を見出すことである。式（５）に従い、次のように書くことができる：
【００３９】
【数１２】
0004813784.tif 000013
【００４０】
ここでピークは単分散であるので、モル質量Ｍはピーク中変らず、積分の外に持ってくることができる。最適のフィットパラメータを定めるために、下流の測定された濃度信号ｃm（ｔ）と拡大された上流の光散乱信号との間のカイ自乗を行なって、それを単分散サンプルピークで積算する。
【００４１】
【数１３】
0004813784.tif 000014
【００４２】
χ2を最小にする、Ｋ*Ｍおよびαｉの最適のフィットパラメータを見つけるために、標準非線形最小二乗フィッティングパッケージを使用することができる。これらの最適フィットパラメータαｉ′を記録しておく。
【００４３】
実際には、上記のフィット値に対して他のシステムの特定のパラメータを含めることができる。たとえば、式（１２）の過剰レイリー比は、サンプルに対する溶媒のレイリー比と純粋溶媒のレイリー比との差として定義することができる。純粋溶媒のレイリー比は、溶媒のベースラインを設定することによって手作業で定める代わりに、調節可能なフィットパラメータとして含めることができる。同様に、濃度ディテクタに緩やかなドリフトがあるのなら、ドリフトのオフセットおよび勾配は同様にフィットに含めることができる。式（１２）における以下の内部積分は測定されたレイリー比および拡大関数のみを含むよう、モル質量がピーク中一定であることが絶対必要であることに留意しなければならない。
【００４４】
【数１４】
0004813784.tif 000015
【００４５】
質量が一定でなければ、固有の最適フィットパラメータαｉ′を定めることは不可能である。
【００４６】
最適のフィットパラメータが定められると、以降の分析において用いられる方法は二通りある。１つは、式（１０）のデコンヴォルーションを行なうことによって濃度測定値を「狭める」ようにして、それにより光散乱ディテクタと一致していたのなら濃度ディテクタはどんな値を測定したであろうかという問いに答えようとするものである。代替的に、光散乱結果を人工的に拡大して、濃度ディテクタと一致して下流で行なわれたのなら光散乱結果はどんな値をもたらしたであろうかという概念的問いに答えることができる。従来のアプローチを成している前者の方法は数値的に不安定であり、しばしば反物理的な結果をもたらすのに対して、後者の方法は数値的に安定しているが、測定値の分解能を落とさなければならない。拡大は一般に小さな補正であるので、分解能の低下は最小のものであるので、後者の方法が好ましい。
【００４７】
従来のデコンヴォルーション法を考えてみる。次のように、フーリエ変換を書き出すことにより、既知のｃm（ｔ）からｃ（ｔ）を組み立てることができる：
【００４８】
【数１５】
0004813784.tif 000016
【００４９】
本特許の主題である第２の方法は光散乱ピークを拡大することである。したがって、次の式が示される：
【００５０】
【数１６】
0004813784.tif 000017
【００５１】
これはディテクタ間拡大が補正された、時間を関数としたモル質量の測定値をなす。式（１５）の右辺は測定可能な量のみを含み、前に定められたパラメータで評価された拡大関数を伴う。
【００５２】
拡大関数のモデル
上記のアルゴリズムによって、拡大カーネルのパラメータ化されたモデルが与えられると、ディテクタ間のバンド拡大現象に対して光散乱結果を補正することができる。しかし、アルゴリズムを適用するためには、使用するモデルを選択しなければならない。ここでは、特定のモデルについて記載し、これはカリフォルニア、サンタバーバラのワイアットテクノロジーによって作成されたＭＡＬＳハードウェアのＤＡＷＮ計器でうまく動く。しかし、このアルゴリズムは他の的確なモデルにも十分にうまく適用できることに留意するのは重要である。
【００５３】
ワイアットのＭＡＬＳ計測で用いられる光散乱セルは直径１．２５ｍｍの筒であり、９０°で直径０．１ｍｍの毛細管から供給され、直径０．２ｍｍの毛細管から排出される。典型的な流量は１ｍｌ／分のオーダにある。フローが層流または乱流であるかどうかを定めるのは、Ｒｅ＝ｖaｄ／ｖであって、ｖaは平均フロー速度であり、ｄはシステムの特性寸法であり、ｖはサンプルの粘度であるレイノルズ数を計算することによって定めることができる。入力毛細管の流れとして、ｖa＝１２７ｍｍ／秒、ｄ＝０．１ｍｍおよびｖ＝０．８９ｃｍ2／秒なので、Ｒｅ＝１４である。流体フローは数百のレイノルズ数に達しない限り不安定にはならない。したがって、入力毛細管のフローは層流であると考える。次に、フローセルに入る流体の噴流を考えてみる。流量は同じであるが、特性寸法はｄ＝１．２５ｍｍに増える。したがって、レイノルズ数はＲｅ＝１７８に増える。したがって、セルの入口近くに何らかの乱流混合があると推論し、これはセル体積の一部の混合体積を有する。同様に、濃度ディテクタフローセルの入口に何らかの混合がある。
【００５４】
さらに、何らかの計測上の拡大があり、これはシステムが流体の有限体積を測定するのであって、測定システムは何らかの固有フィルタ時定数を有することからくる。したがって、私どもが取上げるモデルは指数的拡大であり、計測上の拡大を表わすためにガウス項でたたみ込まれたフローセルの混合を表わしている。拡大カーネルのモデルは以下の式によって示される：
【００５５】
【数１７】
0004813784.tif 000018
【００５６】
ここでσは計測上のガウス幅であり、ｗはセルの入口での混合による拡大である。たたみ込みは以下のように定義される：
【００５７】
【数１８】
0004813784.tif 000019
【００５８】
このモデルの適切性を判定するため、以下のテストが行なわれた。２００ｋＤのポリスチレンをトルエンに溶解した５μｌの注入ループを満たし、これを直接光散乱フローセルに注入した。流量は０．５ｍｌ／分であり、サンプル獲得インターバルは０．２５秒であった。注入ループを表わすために、拡大モデルでたたみ込まれたデルタ関数としてデータをモデル化する。これを９０°の光散乱信号と比較するために、χ2を以下のようにした：
【００５９】
【数１９】
0004813784.tif 000020
【００６０】
ここでｔ0はサンプルが注入器からフローセルに流れるまでかかる時間であり、ａは無次元の単位からレイリー比に変換するためのスケールファクタである。４つのフィットパラメータがある：  ａ、ｔ0、ｇ、およびｗ。
【００６１】
図２はフィットの結果を示す。光散乱データ点４はピークが単一の振幅を有するよう正規化される。拡大関数に対するもたらされるフィットは５である。これは注入ループとフローセルの中央との間で起こる拡大のみを含み、ここで測定が行なわれる。拡大関数はデータに十分適合することは明らかであり、ディテクタ間拡大用の適切なモデルである。この拡大モデルは以下に記載される例で用いられる。
【００６２】
一連の任意ディテクタへの方法の適用
このセクションでは、本方法は直列接続される２つ以上の任意のディテクタ用に一般化される。サンプルが連続的にＤ1からＤ2に、…、Ｄnに流れる、ディテクタの応答をＤi（ｔ）として定義する。一般に、サンプルピークはシステムを通って溶出されるにつれ累進的に広がる。データ（Ｄ1（ｔ），Ｄ2，（ｔ），…，Ｄn（ｔ））の分析は、この累進的拡大によって損なわれる。上述のように、上流のディテクタを最も広い応答を備えたディテクタと一致させることが数字的により安定する。一般にこれはチェインの最後のディテクタである。しかし、中間のディテクタのどれかが大きな計測上の拡大を有するのなら、結果はそれを基準にすることができる。説明のため、ここでは上流のディテクタの信号はＤｎと一致するよう広げられると仮定する。
【００６３】
プロシージャは２つの段階を有する。第１の段階において、拡大パラメータを定める。これは、ディテクタ信号が拡大がなければ正比例するのに十分低い濃度で単分散サンプルを測定することによって行なう。実際には、ほぼ単分散するサンプルをクロマトグラフに注入し、拡大パラメータを定めるための、単量体サブピークを選択する。１からｎ−１の各ディテクタｉに対して以下を計算する：
【００６４】
【数２０】
0004813784.tif 000021
【００６５】
ここでχi2（βi，τi，αij）はｉ番目のディテクタに対するカイ2乗パラメータである。パラメータβiはディテクタの信号を同じスケールに持ってくる比例係数であり、τiはサンプルがあるディテクタから次のディテクタに進むのにかかる時間に伴うディテクタ間遅延であり、αijはｉ番目のディテクタに対するｊ拡大パラメータである。これらのχ2は標準非線形最小二乗フィットアルゴリズムを用いて最小化されて、最適フィットパラメータβi，τi，およびαijを定める。このようなフィッティングが可能なソフトウェアは市場に出回っている種々のパッケージで容易に利用できる。そのプロシージャおよび多様なソフトウェアによるアプローチの一部はACM Transactions on Mathematical Software、第７巻、頁１−１６にある「非線形最小二乗問題を解く数学的ソフトウェアの評価」と題された１９８１年度のHiebertの文献において詳細に記載されている。
【００６６】
第２の段階では、第１段階で定められた拡大パラメータを以降のデータ実行を補正するために適用する。この段階において、上流ディテクタの信号はすべてチェインにおける最後のディテクタと一致するかのように拡大される。次に、これらは直接最後のディテクタＤnと比べることができる。拡大されたデータは次のように表わされる：
【００６７】
【数２１】
0004813784.tif 000022
【００６８】
次に、元の分析はデータ（Ｄ1b（ｔ），…，Ｄn-1b（ｔ），Ｄn（ｔ））で進められる。
【００６９】
アルゴリズムの実証
このセクションでは、ＭＡＬＳ＋ｄＲＩ、粘度計＋ｄＲＩ、およびＵＶ＋ＭＡＬＳの計器の組合せに対するサンプルデータのアルゴリズムを示す。図１は日本の昭和電工によって製造されたShodex（登録商標）ＯＨパックＫＷたんぱく質分離カラムに１００μｌのＢＳＡを注入することによって得られたデータを示す。流量は０．５ｍｌ／分である。細分化されたサンプルはまず光散乱ディテクタを通り、次にｄＲＩディテクタを通る。ＢＳＡサンプルは主に単量体からなるが、カラムによって分離された凝集体をかなりの量含有する。１６分から１８分の間のピークは主に純粋単量体からなる。１４．４分および１６分間のピークは純粋二量体であり、後ろのピークは分離が十分ではないより高い凝集である。横軸は分単位の時間である。トレース１は単量体モル質量Ｍw0によって除算される演算されたモル質量Ｍw（ｔ）である。したがって、単量体は縦軸において１と読取られ、二量体は２、などと読取られる。トレース２は９０°の光散乱信号であり、トレース３は屈折率信号である。屈折率および光散乱はピークが同じ最大高さを有するよう正規化される。トレース２および３に対する縦のスケールは任意である。
【００７０】
アルゴリズムの第１のステップは、拡大モデルを選択することにある。図１に設定されるデータに対して、式（１６）に記載されるハイブリッドモデルが用いられた。アルゴリズムにおける第２のステップは、フィットパラメータを定めるための単分散ピークを選択することにある。１６分および１８分間の単量体ピークが用いられた。式（１９）のχ2を形成すると、非線形最小二乗フィッティングによって４つのパラメータを定めなければならない。２つのデータセット間のベースラインのずれを考慮して余分なパラメータも加えられた。したがって、最小化されたχ2は以下のとおりである：
【００７１】
【数２２】
0004813784.tif 000023
【００７２】
ここでＬＳ（ｔ）は時間を関数とした９０°光散乱信号であり、ｄＲＩ（ｔ）は時間を関数としたｄＲＩデータである。このモデルは、マルカルト(D. W. Marquardt)のJournal
of the Society of Industrial and Applied Mathematicsの１９６３年の論文「非線形パラメータの最小二乗推定のためのアルゴリズム」第１１巻、頁４３１−４４１に記載されているような民間のマルカルト非線形最小二乗パッケージを用いることによって最小化された。式（２０）によって記載される拡大補正を行なうためには、βまたはｘ0を必要としないので、以下の結果を示す表には表わされないが、非線形最小化は正しく働いたことを確実にするために含まれている。
【００７３】
【表１】
0004813784.tif 000024
【００７４】
これらの結果は、σ□ｗが予想通り、主要拡大現象は混合によるものであり、ガウス分散または計測的拡大ではないと示していることに留意することよって解釈できる。τ0パラメータはサンプルが２つの計器間を進むのに２１秒かかったことを示しているだけである。この点について、式（２０）は次のように書替えることができる：
【００７５】
【数２３】
0004813784.tif 000025
【００７６】
図３は上記の第１ステップで見出されたフィットパラメータで式（２２）を用いて計算された拡大光散乱信号７を示す。光散乱分析は前述のように、元の光散乱データではなく、トレース７を用いて進められる。その結果はトレース６に示される。注意しなければならないいくつかの特徴がある。第１に、トレース７は単量体ピークにおいてトレース２と重畳する。これはモデルが正しく、かつ非線形最小化が正しく働いたことを示す。第２の観測するべきことは、１４．５分および１６分間の二量体ピークは予想されるように平坦である。１３．５分および１４．５分間の三量体ピークも平坦であるが、そのベースラインが隣接するピークと分けられていないので、その結果はそれほどはっきりしたものではない。次に観測するべきことは、平坦領域が互いの整数倍であることである。これは予想されることである。なぜなら二量体のモル質量は単量体の質量のちょうど２倍となるべきだからである。同様に、三量体の平坦域は単量体の３倍である。
【００７７】
図４はShodex（登録商標）ＯＨパックＫＷたんぱく質分離カラムに１００μｌのＢＳＡを注入することによって得られたデータを示す。流量は０．５ｍｌ／分である。細分化されたサンプルはまずオンラインブリッジ粘度計を通り、次にｄＲＩディテクタを通った。粘度計は以下のように定義される比粘度を測定する：
【００７８】
【数２４】
0004813784.tif 000026
【００７９】
ここでη（ｔ）は流体サンプル組合せの粘度であり、η0は純粋流体の粘度である。ｄＲＩディテクタはサンプルの屈折率を測定し、これはサンプルｄｎ／ｄｃが既知であるのなら濃度に変換することができる。次に固有粘度を計算することができ、これは以下のように定義される：
【００８０】
【数２５】
0004813784.tif 000027
【００８１】
濃度が十分に低ければ、固有粘度は次のように近似化することができる：
【００８２】
【数２６】
0004813784.tif 000028
【００８３】
トレース９は比粘度であり、トレース１０はｄＲＩ信号である。トレース９およびトレース１０は同じ最大高さを有するよう正規化されている。これら２つのトレースの縦のスケールは任意である。この測定の結果はトレース８で示され、これはｍｌ／ｇ単位の縦軸の固有粘度を表わす。固有粘度は各サンプルピーク中一定であるべきである。見られる湾曲はディテクタ間拡大によるものであるが、この場合、粘度計信号は実質的にｄＲＩ信号より広い。これは粘度計がｄＲＩと比べてより大きな流体量で粘度を測定するからである。この実施例では、主要拡大は粘度計に内在する計測上の拡大である。実際、ｄＲＩ信号はサンプルが２つのディテクタを通ることによって拡大され、これによりピーク幅の全体の差を減らす。本発明の方法は前述のように適用できるが、この場合、ｄＲＩ信号は粘度信号と一致するよう拡大される。アルゴリズムの第１のステップは拡大モデルを選択することである。ここでも、ハイブリッドガウス混合モデルが用いられた。アルゴリズムの第２のステップは、フィットパラメータを定めるための単分散ピークを選択することである。１５．２分と１６．２分との間の単量体ピークが用いられた。前述のように、式（１９）は次のように書き直すことができる：
【００８４】
【数２７】
0004813784.tif 000029
【００８５】
フィットパラメータは調整可能なパラメータに対してχ2を最小化することによって定められる。次のフィットパラメータが見出された：
【００８６】
【表２】
0004813784.tif 000030
【００８７】
この場合、パラメータｗは９．０４秒であり、これはディテクタ間混合から起こるのではなく、ｄＲＩフローセルの体積に相関する粘度計毛細管の大きな内部体積から起こる。ガウス項は非常に小さく、これは分散が重要ではないことを示す。ｄＲＩ信号は式（２０）を次のように書き直すことによって広げることができる：
【００８８】
【数２８】
0004813784.tif 000031
【００８９】
分析は広げられたｄＲＩ信号で進めることができる。図５では拡大されたｄＲＩ信号がトレース１２で示され、補正された固有粘度信号がトレース１１に示される。結果は各ピークに対する平坦域からなる。
【００９０】
図６に示される最後の実施例は３．０ｍｇ／ｍｌのＢＳＡのサンプル１００μｌからなり、ＵＶディテクタによって測定され、続いてＭＡＬＳディテクタによって測定された。
この場合、１４に示されるＬＳ信号は１５で示されるＵＶ信号より著しく広く、１３に示される湾曲したモル質量トレースをもたらす。この実施例の場合、ＵＶ信号はＬＳ信号と一致するよう拡大される。式（１９）は次のように書き直される：
【００９１】
【数２９】
0004813784.tif 000032
【００９２】
ここでもハイブリッド混合およびガウスモデルが用いられ、１３．５分と１５分との間の単量体ピークがフィットパラメータを定めるために用いられた。フィットパラメータに対してχ2を最小することによって次のようになる：
【００９３】
【表３】
0004813784.tif 000033
【００９４】
ここでも、拡大の主要な原因は内部混合である。ＵＶデータは図７に示されるように、以下のように書き直された式（２０）を計算することによって広げられている：
【００９５】
【数３０】
0004813784.tif 000034
【００９６】
拡大されたＵＶ信号はトレース１７に示され、モル質量はトレース１６で示される。前述の２つの場合と同様に、データは平坦域からなり、ディテクタ間拡大現象は著しく減じられた。
【００９７】
上記の本発明の方法の実施例は、本明細書に開示される発明の方法によって対処および補正することができる可能な種類のバンド拡大補正のほんの数例しか示されていない。クロマトグラフィ分野の当業者にとって自明であるように、本発明以前は、バンド拡大の現象は複数のディテクタの使用に伴う測定値の精度を向上させるための著しい障害となっていた。本発明の方法およびその変形はバンド拡大現象が著しく減じられるという期待を満たす。
【図面の簡単な説明】
【００９８】
【図１】ｄＲＩディテクタがＭＡＬＳディテクタの下流にある、ＢＳＡサンプルのモル質量の補正されていないプロット図である。
【図２】トルエンに２００ｋＤポリスチレンを５μｌ注入した９０°光散乱トレースに適合されたハイブリッド拡大モデル図である。
【図３】本発明の方法によって補正された図１のデータを示す図である。
【図４】ｄＲＩディテクタが粘度計の下流にある、ＢＳＡサンプルに対する固有粘度対溶出量の補正されていないプロット図である。
【図５】トレース１１が補正された固有粘度である、本発明の方法によって補正された図４のデータを示す図である。
【図６】ＵＶディテクタがＭＡＬＳの上流にある、ＢＳＡサンプルのモル質量の補正されていないプロット図である。
【図７】本発明の方法によって補正された、図６のデータを示す図である。
【符号の説明】
【００９９】
１−３，６−１３  トレース、１４  ＬＳ信号、１５  ＵＶ信号。
【図１】
0004813784.tif 000039
【図２】
0004813784.tif 000040
【図３】
0004813784.tif 000041
【図４】
0004813784.tif 000042
【図５】
0004813784.tif 000043
【図６】
0004813784.tif 000044
【図７】
0004813784.tif 000045