マイクロリアクタ、解離定数測定装置、解離定数測定方法、及び解離定数測定プログラム
【課題】 反応の初期段階において解離定数の予測を行えるようにすること。
【解決手段】 質量センサの下流側にサンプルの濃度センサを備える。これにより質量センサで質量が検出された際のサンプル液の濃度を知ることができる。濃度と質量を所定の理論式に代入すると解離定数の予測値を得ることができる。濃度センサは、対向する電極で構成されたコンデンサとする。電極間にアナライトが存在すると誘電率が変化するため、誘電率によりアナライト濃度を検出することができる。また、コンデンサは小型化できるため、マイクロリアクタを小型に保つことができる。なお、質量センサは水晶振動子の共振周波数を出力し、測定装置はこの周波数の遷移から質量を算出することができる。
【解決手段】 質量センサの下流側にサンプルの濃度センサを備える。これにより質量センサで質量が検出された際のサンプル液の濃度を知ることができる。濃度と質量を所定の理論式に代入すると解離定数の予測値を得ることができる。濃度センサは、対向する電極で構成されたコンデンサとする。電極間にアナライトが存在すると誘電率が変化するため、誘電率によりアナライト濃度を検出することができる。また、コンデンサは小型化できるため、マイクロリアクタを小型に保つことができる。なお、質量センサは水晶振動子の共振周波数を出力し、測定装置はこの周波数の遷移から質量を算出することができる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、マイクロリアクタ、解離定数測定装置、解離定数測定方法、及び解離定数測定プログラムに関し、例えば、異常蛋白質と結合する化合物の探索に用いるものに関する。
【背景技術】
【０００２】
近年、ヒトゲノム（ヒトの遺伝子情報）の解析が終了し、異常な遺伝子構造が生成する異常蛋白質と病気との関係解明が進められつつある。
【０００３】
この関係の解明により、新薬の開発手法が、開発者の薬剤や化合物に対する経験と勘を頼りに行う既存手法から、異常蛋白質に直接作用する化合物を探索して新薬とする手法へと変化している。この手法の採用により、２０年近く必要としていた新薬開発の期間が、今後５年程度に短縮すると見込まれている。
【０００４】
新薬候補の化合物の探索には、異常蛋白質と新薬候補の化合物との物理的な結合力を表す定数である解離定数を指標として用いるのが一般的である。
【０００５】
解離定数の測定方法としては、以前は、酵素、発光物質、放射性同位元素などの標識物質を結合させた化合物を用い、この化合物と異常蛋白を結合させた後、標識物質の量を測定することで、結合した化合物を定量し、解離定数の算出を行っていたが、現在では、標識物質を用いずにマイクロリアクタで行う測定が注目されている。
【０００６】
マイクロリアクタとは、半導体やガラス、樹脂などで構成されたチップの中に導入路と廃液路を形成し、その間に反応槽を設けたデバイスである。反応槽には、予め異常蛋白質を固定したセンサが設置される。
【０００７】
このように構成されたマイクロリアクタにおいて、導入路から化合物を含むサンプル液を流し込むと、サンプル液中の化合物が反応槽に予め設置された異常蛋白質と反応し、反応後の廃液が廃液路から排出される。
【０００８】
なお、反応槽で予め固定されている物質はリガンドと呼ばれ、溶液として供給される物質はアナライトと呼ばれるため、以下ではこの用語を用いる。
【０００９】
反応槽においては、サンプル液中のアナライトのうちあるものはリガンドと結合してセンサに固定される。一方、センサに固定されていたアナライトのうち、リガンドから離れてサンプル液中に放出されるものもある。
【００１０】
アナライトとリガンドの反応が平衡状態に達すると（即ち、センサに固定されるアナライトの量と、センサから離れるアナライトの量が等しくなると）、リガンドに固定されているアナライトの量が一定量となる。このとき計測されるアナライトとリガンドとが結合した化合物、結合していないリガンド、結合していないアナライトの個数（モル濃度）から解離定数が算出され、この解離定数が新薬の探索で必要なデータとなる。
【００１１】
このようなマイクロリアクタに関する技術としては、次の「マイクロリアクタ」がある。
【特許文献１】特表２００３−５２６３５９公報 この技術は、微小な混合室、加熱／分散ユニット、反応室、分離ユニットといった微小反応器をチップに並列に形成し、大量の生成物を生産するものである。
【００１２】
また、リガンドとアナライトの結合を計測する技術としては次の「免疫検査方法」がある。
【特許文献２】特許２９６２０３１ この技術は、リガンドとアナライトの結合を音響センサを用いて検出するものである。
【発明の開示】
【発明が解決しようとする課題】
【００１３】
従来は、解離定数を求めるためにはリガンドとアナライトの反応が平衡状態に達する必要があり、平衡状態に達するまでに長時間を要するため、解離定数を得るのに時間がかかるという問題があった。探索に係るリガンドとアナライトの組み合わせは膨大であるため、全ての組み合わせに対して反応が平衡に達するまで実験を継続するのはコスト的に問題があった。
【００１４】
また、長時間実験を行うことにより高価なサンプルを多量に消費するという問題もあった。
【００１５】
反応の初期の段階で、有望な組み合わせであるか否かをある程度予測できるようにしたいとの要望もあるが、正確な予測は困難であった。ここで、その理由について説明する。
【００１６】
図１０（ａ）は、センサに固定されるアナライトの質量の時間変化を表したグラフである。縦軸は、センサに固定されたアナライトの質量を表し、横軸は時間軸を表している。
【００１７】
まず、区間１０１では、緩衝液がマイクロリアクタに供給される。緩衝液にはアナライトが含まれていないため、センサに固定されるアナライトの質量は０である。このように、サンプル液に先立って緩衝液を流すのは、センサにリガンドとして固定されている異常蛋白質を乾燥させないためである。
【００１８】
マイクロリアクタに供給する溶液を緩衝液からサンプル液に切り替えると、区間１０２に示したようにセンサで固定されるアナライトの質量が徐々に増加してくる。そして、リガンドとアナライトの反応が平衡に達すると区間１０３で示したようにセンサで固定されるアナライトの質量は一定値となる。
【００１９】
ここで、例えば、図１０（ｂ）の時点１０５において、それまでの曲線の形からその後のアナライトの質量を曲線１０７のように予測したとする。しかし、実際に計測を続行すると例えば曲線１０８が得られ、予測値から大きく外れたものとなる。
【００２０】
このように予測が困難となる大きな原因としては、センサに供給する溶液を緩衝液からサンプル液への切り替える際に、センサでサンプルの濃度が徐々に変化するため、センサでのサンプル液の濃度が反応初期段階ではわからないことが挙げられる。
【００２１】
図１０（ｃ）はセンサでのサンプル濃度の時間変化を示したグラフである。曲線１１０に示したように、サンプルの濃度は時間と共に徐々に上昇する。これは、供給する溶液を緩衝液からサンプル液に切り替えた後、導入路内において両者が拡散混合したり、あるいは、流路内の流速分布により、重いアナライトが流れの速い流速域に集中したり、あるいは、センサに至るまでの導入路壁にアナライトが吸着するために起こるものである。
【００２２】
破線１１１に示したように、緩衝液からサンプル液にある時点で二者択一的に変化するのであれば、センサでのサンプルの濃度がわかるが、実際はセンサで質量変化が検出されても、その時点に至るまでのサンプルの濃度変化の経過がわからないため、予測が困難となる。
【００２３】
そこで、本発明の目的は、反応の初期段階において解離定数の予測を行えるようにすることである。
【課題を解決するための手段】
【００２４】
本発明は、前記目的を達成するために、試料を含む溶液を流す流路と、前記流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、を具備したことを特徴とするマイクロリアクタを提供する（第１の構成）。
【００２５】
第１の構成において、前記質量取得手段は、前記獲得物質を固定する共振部材の共振周波数を質量情報として取得するように構成することもできる（第２の構成）。
【００２６】
第１の構成において、前記質量取得手段は、前記獲得物質の上を流れる溶液を伝播する弾性波の速度を質量情報として取得するように構成することもできる（第３の構成）。
【００２７】
第１の構成、第２の構成、又は第３の構成において、前記濃度取得手段は、試料が獲得された後の溶液の誘電率情報を濃度情報として取得するように構成することもできる（第４の構成）。
【００２８】
また、本発明は、試料を含む溶液を流す流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、前記取得した質量情報と濃度情報を用いて、前記試料と前記獲得物質の解離定数を算出する解離定数算出手段と、を具備したことを特徴とする解離定数測定装置を提供する（第５の構成）。
【００２９】
また、本発明は、第１の構成から第４の構成までのうちの何れか１の構成のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求める解離定数算出方法であって、結合物質量算出手段と、未獲得量算出手段と、解離定数算出手段と、を備えたコンピュータにおいて、前記結合物質量算出手段によって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出ステップと、前記未獲得量算出手段によって、前記算出した結合物質の量を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出ステップと、前記解離定数算出手段によって、前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出ステップと、から構成されたことを特徴とする解離定数算出方法を提供する（第６の構成）。
【００３０】
また、本発明は、第１の構成から第４の構成までのうちの何れか１の構成のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求めるための解離定数算出プログラムであって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出機能と、前記算出した結合物質の量の時間変化を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出機能と、前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出機能と、をコンピュータで実現する解離定数算出プログラムを提供する（第７の構成）。
【発明の効果】
【００３１】
本発明によると、リガンドとアナライトの反応の初期段階において解離定数を予測することができる。
【発明を実施するための最良の形態】
【００３２】
（実施の形態の概要）
アナライトを獲得する質量センサの下流側近傍にサンプルの濃度センサを備える。これにより質量センサで質量が検出された際のサンプル液の濃度を知ることができる。
【００３３】
濃度と質量を後述の所定の理論式に代入すると解離定数の予測値を得ることができる。
【００３４】
濃度センサは、対向する電極からなるコンデンサで構成する。電極間にアナライトが存在すると誘電率が変化するため、電気容量を計測することによりアナライト濃度を検出することができる。
【００３５】
また、コンデンサは後述のすだれ構造の採用などにより小型化できるため、マイクロリアクタを小型に保つことができる。
（実施の形態の詳細）
図１は、本実施の形態のマイクロリアクタの構造を模式的に表した平面図である。
【００３６】
マイクロリアクタ１は、例えばガラスなどで構成された筐体部材により構成されており、この筐体部材内に各種の構成要素が形成されている。
【００３７】
マイクロリアクタ１の外寸は、長さ（導入路３の方向）が５０〜１００［ｍｍ］程度、幅が１０〜２０［ｍｍ］程度、厚さが数［ｍｍ］程度である。
【００３８】
マイクロリアクタ１には、導入路３、反応槽４、廃液路５などが形成されている。そして、反応槽４には質量センサ７が設置され、質量センサ７の廃液路５側には質量センサ７と近接して濃度センサ９が設けられている。
【００３９】
なお、緩衝液とサンプル液の流れを制御するバルブや、質量センサ７、濃度センサ９を構成する配線などは図示していない。
【００４０】
導入路３は、マイクロリアクタ１の外部から反応槽４に溶液を導く流路であり、緩衝液、及びアナライトを含んだサンプル液はこの管を通って反応槽４に導かれる。
【００４１】
一方、廃液路５は、反応槽４から排出される溶液をマイクロリアクタ１の外部に導く流路であり、反応槽４で反応を終えたサンプル液はこの管を通って排出される。
【００４２】
本実施の形態では、廃液路５の側から溶液を吸引することにより、導入路３から溶液が供給されるように構成した。
【００４３】
導入路３の側にポンプを設置して、溶液をマイクロリアクタ１に供給することも可能であるが、この場合、ポンプなどに前の実験で使用したアナライトが残留し、今回の測定に混入してしまうことで、測定精度に影響する場合があり（コンタミネーション）、本実施の形態では、溶液を吸引してマイクロリアクタ１に供給することにより、このような問題を防止した。
【００４４】
質量センサ７（質量取得手段）は、表面にリガンドが固定してあり、着脱可能に構成されている。以下に図２を用いて質量センサ７の構造について説明する。
【００４５】
図２は反応槽４の断面を表している。質量センサ７は、金薄膜１３、リガンド１１、水晶振動子１４などから構成されている。
【００４６】
金薄膜１３の上には、図示しない自己組織化膜を介してリガンド１１が予め固定されている。また、金薄膜１３の下面には水晶振動子１４が設けられている。水晶振動子１４は、コルピッツ発信回路の一部を成しており、水晶振動子１４を厚みすべり振動モードで共振させ、金薄膜１３をサンプルの流れと平行な面内で振動させるようになっている。ここで、金薄膜１３、リガンド１１、水晶振動子１４は共振部材を構成している。
【００４７】
質量センサ７の上を流れるアナライト１６のうち、あるものはリガンド１１と結合して（例えば、アナライト１６ａ）質量センサ７に固定される。
【００４８】
また、一旦質量センサ７に固定されたアナライトであっても、アナライト１６ｂのように、リガンドとの結合が切れて再びサンプル液中に放出されるものもある。
【００４９】
このように、リガンド１１は、試料（アナライト１６）を含む溶液から試料を獲得する獲得物質を構成している。
【００５０】
質量センサ７にアナライトが固定されると質量センサ７の質量が増すため、質量センサ７の共振周波数がｆｒからｆｒ＋Δｆへシフトする。このシフト量Δｆに次の式（１）を適用することにより質量センサ７に固定されたアナライトの質量Δｍを計算することができる。ここで、Ａは金薄膜１３の面積であり、ρqは水晶の密度であり、μqは水晶の剪断応力である。
【００５１】
なお、共振周波数は質量センサ７に固定されたアナライトの質量を表す質量情報を構成している。
【００５２】
Δｆ＝２ｆｒ2Δｍ／｛Ａ（ρqμq）1/2｝・・・（１）
次に、濃度センサ９（濃度取得手段）について説明する。
【００５３】
濃度センサ９は、図３（ａ）に示したように、電極２１と電極２２からなるコンデンサにより構成されている。図３（ｂ）は、電極２１、２２の断面を示している。
【００５４】
濃度センサ９は、質量センサ７と同一平面上に近接して配置されている。濃度変化は流路の高さ方向で異なるので、質量センサ７と同一平面に濃度センサ９を設けることで、解離定数の高精度予測が可能となる。
【００５５】
電極２１、２２は、廃液路５の底面に凸状に形成された電極部が廃液路５内で対向するように配置されて構成されている。
【００５６】
濃度センサ９の上をサンプル液が通過すると、対向した電極２１、２２の間にサンプル液中のアナライトが出入りし、このアナライトの量で電極２１、２２間の誘電率が変化する。
【００５７】
誘電率の変化と電気容量の変化は対応させることができるため、電極２１、２２間の電気容量によりサンプル液の濃度を計測することができる。ここで、電気容量は誘電率情報、及び濃度情報を構成している。
【００５８】
また、本実施の形態では、電極２１、電極２２の形状を１対のすだれ形状とすることにより、小さい設置領域面積で大きな電気容量の変化を得ることができるようにしてある。
【００５９】
即ち、すだれ状に電極を配置することにより、図３（ｂ）に示したように、電極２１ａと電極２２ａ、電極２２ａと電極２１ｂ、電極２１ｂと電極２２ｂの３つのコンデンサを形成することができる。
【００６０】
このように濃度センサ９を小型化できることから、濃度センサ９の設置によるマイクロリアクタ１の大型化を防ぐことができ、これによって、貴重な試料の使用量を抑えることができる。
【００６１】
以上のようにマイクロリアクタ１は、質量センサ７の下流側近傍に濃度センサ９を設けたため、質量センサ７でΔｍが検出された際のサンプル液の濃度を知ることができる。
【００６２】
なお、マイクロリアクタ１では、水晶振動子１４の共振周波数の遷移からΔｍを求めたが、弾性波の伝播速度によりΔｍを測定することも可能である。
【００６３】
一般に、共振周波数を用いるタイプの質量センサはＱＣＭ（Ｑｕａｒｔｚ Ｃｒｙｓｔａｌ Ｍｉｃｒｏｂａｌａｎｃｅ）と呼ばれ、弾性波を用いるタイプのものはＳＡＷ（Ｓｕｒｆａｃｅ Ａｃｏｕｓｔｉｃ Ｗａｖｅ）センサと呼ばれている。
【００６４】
図４は、ＳＡＷセンサを搭載したマイクロリアクタ１ａの構成を模式的に示した図である。
【００６５】
質量センサ７の上流側近傍には弾性波発生電極２５が設けられている。
【００６６】
弾性波発生電極２５の構成は、濃度センサ９と同じであり、１対の電極がすだれ状に配置されている。弾性波発生電極２５に電圧を印可すると、例えば、図示しない圧電素子が伸縮して弾性波が発生する。
【００６７】
濃度センサ９は、弾性波の検出を行うと共に、弾性波の検出の合間においてサンプル液の濃度の検出も行う。
【００６８】
このように、ＳＡＷセンサを搭載したマイクロリアクタでは、下流側のコンデンサに濃度検出手段と、弾性波検出手段を兼用させることができる。
【００６９】
弾性波発生電極２５と濃度センサ９の間の距離は予めわかっているため、弾性波発生電極２５が弾性波を発生してから濃度センサ９がこの弾性波を検出するまでの時間差により質量センサ７上を通過した弾性波の速度を求めることができる。
【００７０】
質量センサ７でアナライトが捕獲されると弾性波の伝達速度が遅くなることが知られており、所定の理論式により、弾性波の速度から質量センサ７で固定されたアナライトの量を計算することができる。
【００７１】
なお、ＳＡＷセンサを搭載したマイクロリアクタの場合、質量を検出するための領域（質量センサ７）は、対向する弾性波発生電極２５と濃度センサ９の幅で十分であり（これら電極の幅以外の領域はアナライトを結合しても検出できないため）、反応漕を円形にする必要はない。また、流路の幅を一定にすることができるため、この領域での圧力損失がなく、サンプル液の送液制御が容易になるという利点もある。
【００７２】
次に、マイクロリアクタ１を用いた解離定数の予測方法について説明する。
【００７３】
ここでは、まず、図５を用いて予測に用いる理論式について説明する。以下では、リガンドをＬで表し、アナライトをＡで表す。更にリガンドとアナライトが結合した化合物をＬＡで表す。なお、［］は、その内部に記された物質のモル濃度を示す。
【００７４】
まず、ＬＡ、Ｌ、Ａの間には図５の式（２）で表した関係が成り立つ。
【００７５】
即ち、リガンドとアナライトの結合により質量センサ７で固定されているＬＡのモル数［ＬＡ］の微小時間内での変化量は、質量センサ７でアナライトが固定されて生成されるＬＡのモル数（ｋｐ［Ｌ］［Ａ］）から、一旦捕らえられたアナライトが質量センサ７から放たれて減少するモル数（ｋｍ［ＬＡ］）を減算した値となる。平衡状態では、右辺の第１項と第２項は等しくなる。
【００７６】
なお、質量センサ７で、Ｌ＋Ａ→ＬＡの反応式により固定されるＬＡのモル数は、質量センサ７でまだアナライトと結合していないリガンドのモル数［Ｌ］とサンプル液中のアナライトの濃度（モル濃度）［Ａ］をかけたものに比例し、その比例定数をｋｐとしている。
【００７７】
また、質量センサ７で形成されたＬＡのうち、ＬＡ→Ｌ＋Ａの反応式により結合しているアナライトがサンプル液中に離脱するもののモル数は、質量センサ７で形成されているＬＡのモル数に比例し、その比例定数をｋｍとしている。
【００７８】
式（２）をｋｐ、ｋｍについて解くと、それぞれ式（３）、式（４）が得られる。
【００７９】
一方、式（２）を差分で表すと、式（５）が得られる。これは、微小時間ｄＴの間隔で、［ＬＡ］、［Ｌ］、［Ａ］をサンプリングした値を用いて式（２）を書き換えたものである。添え字のｉは、何番目のサンプリング値であるかを示している。更に、式（５）を［ＬＡ］iについて解くと式（６）が得られる。
【００８０】
ところで解離定数Ｋｄはｋｍ／ｋｐで表され、また平衡状態ではｄ［ＬＡ］／ｄｔ＝０であるので、平衡状態ｃでの解離定数Ｋｄは、式（７）で表される。
【００８１】
また、［Ｌ］c＝［Ｌ］0−［ＬＡ］c（［Ｌ］0はＬＡが形成される前のリガンドのモル数）となるので式（８）が導かれ、更に式（８）から式（９）が導かれる。
【００８２】
次に、データのサンプリングについて説明する。
【００８３】
後述する測定装置は、サンプリング周期ｄＴでマイクロリアクタ１から質量センサ７の共振周波数と濃度センサ９の電気容量を取得し、これらを解析して質量センサ７に固定されたアナライトの質量Δｍと、その時のサンプル液中のアナライトのモル濃度［Ａ］を算出する。
【００８４】
これにより、ｄＴ間隔でΔｍと［Ａ］を得ることができる。
【００８５】
測定装置がサンプリングしたこれらの値をプロットしたのが図６（ａ）、（ｂ）である。
【００８６】
図６（ａ）は、質量センサ７で検出されたアナライトの質量を表しており、図６（ｂ）は、濃度センサ９で検出されたサンプル液中のアナライトのモル濃度を表している。両グラフの時間軸は対応させてある。
【００８７】
図６に示したように、測定装置により、時刻Ｔ0にてΔｍ0と［Ａ］0が得られ、時刻Ｔ1にてΔｍ1と［Ａ］1が得られ、以下同様にして、ｄＴごとにΔｍと［Ａ］が検出される。
【００８８】
次に、図７のフローチャートを用いて、測定装置が解離定数を予測する手順について説明する。
【００８９】
まず、測定装置は初期値の設定を行う（ステップ５）。この初期化では、質量センサ７に予め固定されているリガンドのモル数Ｌｅを［Ｌ］0、マイクロリアクタ導入口でのアナライトのモル濃度Ａｃを［Ａ］ｃに設定し、更にアナライトの分子量Ｍａ、リガンドの分子量ＭＬを設定する。これらの設定は、ユーザが測定装置に入力することにより行われる。また、測定装置は自動的に［ＬＡ］を０に設定する。
【００９０】
次に、測定装置は質量センサ７の出力よりΔｍ0を測定し、更に濃度センサ９の出力より［Ａ］0を測定し、これらの値を記憶装置に記憶する（ステップ１０）。
【００９１】
測定装置は、ｄＴ間待った後（ステップ１５）、Δｍ1と［Ａ］1を測定し、記憶する（ステップ２０）。
【００９２】
次に、測定装置は、記憶したこれらの値を用いて［ＬＡ］1と［Ｌ］1を計算する（ステップ２５）。
【００９３】
［ＬＡ］1（質量センサ７で形成されたＬＡのモル数）は、質量センサ７で検出された質量の増分（Δｍ1−Δｍ0）をアナライトの分子量Ｍａで除すれば得られる。
【００９４】
また、［Ｌ］1（質量センサ７でアナライトと結合していないリガンドのモル数）は、リガンドの初期値［Ｌ］0から生成されたＬＡのモル数（［ＬＡ］1−［ＬＡ］0）を減算すれば得られる。
【００９５】
次に、測定装置は、ここで得られている［ＬＡ］1、［Ｌ］1、［Ａ］1を用いてｋｐ1の仮定値を計算して記憶する（ステップ３０）。計算式はフローチャートに示したとおりである。この値は式（３）においてｋｍ［ＬＡ］＝０とおいたものである。リガンドとアナライトの反応が始まった初期段階であるので、質量センサ７から放たれるアナライトの数を０としたのである。
【００９６】
次に、測定装置は、ｄＴ間待った後（ステップ３５）、Δｍ2と［Ａ］2を測定して記憶する（ステップ４０）。
【００９７】
そして、測定装置は、ステップ２５と同様にして、［ＬＡ］2と［Ｌ］2を計算して記憶する（ステップ４５）。
【００９８】
次に、測定装置は、ｋｍ2とｋｐ2の仮定値を計算して記憶する（ステップ５０）。計算式はフローチャートに示したとおりである。ｋｍ2の計算式は、式（５）を用いたものである。ｋｐ2は、ｄＴ間ではｋｐ1とほとんど差がないと考えられるのでｋｐ1と等しいとした。
【００９９】
次に、測定装置は、ｋｐ2、ｋｍ2などを用いて［ＬＡ］3の予測値を計算する（ステップ５５）。計算式はフローチャートに示したとおりである。この計算式は式（６）から得られたものである。以下、図８のフローチャートに続く。
【０１００】
次に、測定装置は、カウンタｉを３にセットし、ｉが所定の整数Ｎに達するまで以下のループ処理を行う（ステップ６０）。なお、以下の処理では、仮定値を用いて解離定数Ｋｄを計算し、その評価を行うが、解離定数Ｋｄが許容誤差の範囲に入った場合は、ループを抜けてこの解離定数Ｋｄを予測値として出力する。また、整数Ｎは、解離定数Ｋｄの予測値が求まるまでループ処理が継続するように十分に大きい値とする。
【０１０１】
次に、測定装置は、ｄＴ待った後（ステップ６５）、Δｍiと［Ａ］iを測定して記憶する（ステップ７０）。
【０１０２】
次に、測定装置は、ステップ２５と同様にして［ＬＡ］iと［Ｌ］iを計算して記憶装置に記憶する（ステップ７５）。
【０１０３】
次に、測定装置は［ＬＡ］iの予測値の評価を行う（ステップ８０）。評価は、ステップ７５で計算された［ＬＡ］iと、［ＬＡ］iの予測値（ｉ＝３の場合はステップ５５で求めた値、その他の場合はステップ１１０で求めた値）の差を［ＬＡ］iで除した値が予め設定されている誤差許容値未満であるか否かを判定することにより行う。
【０１０４】
［ＬＡ］iの計算値と予測値の差が誤差許容値未満であった場合（ステップ８０；Ｙ）、測定装置は、ｉを１だけデクリメントし（ステップ８５）、解離定数Ｋｄと化学平衡時の［ＬＡ］cを計算して出力し（ステップ１０５）、処理を終了する。Ｋｄと［ＬＡ］cの計算式はフローチャートに示したとおりである。なお、［ＬＡ］cの計算式は式（９）による。
【０１０５】
なお、式（９）は、リガンドとアナライトの反応が平衡している場合に成り立つ式（７）から導出した関係式であるが、最終的には、［Ａ］ｃがマイクロリアクタ導入口に流し込むサンプル液濃度Ａｃになると考えたことによる。
【０１０６】
一方、［ＬＡ］iの計算値と予測値の差が誤差許容値未満でなかった場合（ステップ８０；Ｎ）、測定装置は、ｋｍiとｋｐiの仮定値を計算し、記憶装置に記憶する。計算式はフローチャートに示したとおりである。これらの計算式は何れも式（５）を用いたものであり、現在記憶している最新の値であるｋｍi-1、ｋｐi-1を用いて仮定値を更新するようになっている。
【０１０７】
次に、計測装置は、ステップ９０で仮定したｋｍiとｋｐiを用いて、ステップ５５と同様にして［ＬＡ］iの修正予測値ＮＥ［ＬＡ］iを計算して記憶装置に記憶する（ステップ９５）。計算式はフローチャートに示したとおりである。
【０１０８】
次に、計測装置は、ステップ８０と同様にして修正予測値ＮＥ［ＬＡ］iの評価を行う（ステップ１００）。評価式はフローチャートに示したとおりである。
【０１０９】
［ＬＡ］iの計算値と修正予測値ＮＥ［ＬＡ］iの差が予め設定した誤差許容値未満の場合は（ステップ１００；Ｙ）、ステップ１０５にて解離定数Ｋｄと化学平衡時の［ＬＡ］cを算出して出力する（ステップ１０５）。
【０１１０】
［ＬＡ］iの計算値と修正予測値ＮＥ［ＬＡ］iの差が予め設定した誤差許容値未満でなかった場合は（ステップ１００；Ｎ）、ステップ５５と同様にして次の予測値Ｅ［ＬＡ］i+1を計算する（ステップ１１０）。計算式はフローチャートに示したとおりである。
【０１１１】
そして、測定装置は、カウンタｉを１だけインクリメントし（ステップ１１５）、ステップ６０に戻る。そして、同様の処理を繰り返すことになる。
【０１１２】
このように、測定装置は、マイクロリアクタ１より獲得物質（リガンド）によって獲得された溶液中の試料（アナライト）の質量を表す質量情報（共振周波数）を取得する質量取得手段と、前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報（電気容量）を取得する濃度取得手段と、前記取得した質量情報と濃度情報を用いて、前記試料と前記獲得物質の解離定数を算出する解離定数算出手段と、を具備している。
【０１１３】
また、測定装置は、結合物質量算出手段と、未獲得量算出手段と、解離定数算出手段と、を備えたコンピュータで構成することができ、当該測定装置による測定方法は、前記結合物質量算出手段によって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量（［ＬＡ］c）を算出する結合物質量算出ステップと、前記未獲得量算出手段によって、前記算出した結合物質の量を用いて、前記試料を獲得していない獲得物質の量（［Ｌ］c＝［Ｌ］0−［ＬＡ］c）を算出する未獲得量算出ステップと、前記解離定数算出手段によって、前記濃度取得手段で取得した濃度情報（［Ａ］）、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出ステップと、を備えている。
【０１１４】
以上に説明した手順により、測定装置は、マイクロリアクタ１から出力される質量情報と濃度情報をサンプリングすることにより、リガンドとアナライトが反応を始める初期の段階で、解離定数を予測することができる。そのため、ユーザは、リガンドとアナライトの反応の初期の段階で、そのリガンドとアナライトの組み合わせが新薬候補としてどの程度有望か見当を付けることができる。
【０１１５】
これによって、例えば、新薬候補となる見込みがある場合は、リガンドとアナライトの結合が飽和するまで実験を継続し、見込みがない場合は次の組み合わせの実験に移行するといった運用が可能となる。
【０１１６】
図９は、本実施の形態に係る測定装置のハードウェア的な構成の一例を示したブロック図である。
【０１１７】
測定装置１００は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）５１、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）５２、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）５３、入力装置５４、表示装置５５、印刷装置５６、電子回路５７、記憶装置５８、記憶媒体駆動装置５９、入出力Ｉ／Ｆ（インターフェース）６０などの各機能部がバスライン５０で接続されて構成されている。
【０１１８】
ＣＰＵ５１は、所定のプログラムに従って、各種の演算処理、情報処理、及び測定装置１００を構成する各構成要素の制御を行ったりする中央処理装置である。
【０１１９】
ＣＰＵ５１は、電子回路５７を介してマイクロリアクタ１から出力される共振周波数及び電気容量をサンプリングしたり、サンプリングしたこれらの値を用いて質量センサ７で固定されたアナライトの質量やサンプリング液の濃度を算出するほか、これらの値を用いて解離定数Ｋｄの予測値の計算も行う。
【０１２０】
ＲＯＭ５３は、測定装置１００を動作させるための基本的なプログラムやデータなどを記憶した読み出し専用の記憶装置である。
【０１２１】
ＲＡＭ５２は、ＣＰＵ５１が動作するためのワーキングエリアを提供する読み書き可能な記憶装置である。ＲＡＭ５２は、ＣＰＵ５１が図７で示した情報処理を行う際には、Δｍ０、Δｍ１、・・・など、予測に必要な値を記憶する。
【０１２２】
入力装置５４は、例えば、キーボードやマウスなどの入力装置を備えており、ユーザが測定装置１００を使用するのに必要な情報を入力できるようになっている。図７のステップ５で示した値は入力装置５４から入力することができる。
【０１２３】
表示装置５５は、文字情報や画像情報を表示する表示デバイスを備えており、ユーザが測定装置１００を使用するのに必要なメニュー画面などが表示されるようになっている。また、予測した解離定数Ｋｄも表示装置５５に表示することができる。
【０１２４】
表示デバイスは、例えば、液晶ディスプレイ、プラズマディスプレイ、その他のディスプレイ装置で構成することができる。
【０１２５】
電子回路５７は、マイクロリアクタ１に電力を供給して駆動すると共に、マイクロリアクタ１から出力される質量センサ７の共振周波数や濃度センサ９の電気容量を検出する。
【０１２６】
更に、電子回路５７は、マイクロリアクタ１の出力データ（アナログデータ）を時間間隔ｄＴにてサンプリングしてデジタルデータに変換し、ＣＰＵ５１に提供する。
【０１２７】
ＣＰＵ５１は、サンプリングされたデータからΔｍと［Ａ］を計算する。これによって時間間隔ｄＴでΔｍと［Ａ］が得られる。
【０１２８】
印刷装置５６は、例えば、レーザプリンタやインクジェットプリンタなどの印刷装置であり、実験データを印刷できるようになっている。
【０１２９】
記憶媒体駆動装置５９は、装着された着脱可能な記憶媒体を駆動し、データの読み書きを行う機能部である。
【０１３０】
読み書き可能な記憶媒体としては、例えば、フレキシブルディスク、光磁気ディスク、半導体記憶装置、磁気テープ、紙テープがある。
【０１３１】
また、読み取り専用の記憶媒体としては、例えば、ＣＤ−ＲＯＭなどの光ディスク類がある。
【０１３２】
測定装置１００は、記憶媒体駆動装置５９に装着した記憶媒体からプログラムのインストールなどを行うことができるほか、実験データを記憶媒体に書き込むこともできる。
【０１３３】
記憶装置５８は、例えば、ハードディスクなどで構成された大容量で読み書き可能な記憶装置である。
【０１３４】
記憶装置５８には、プログラム類を格納したプログラム格納部６１とデータ類を記憶したデータ格納部６２が形成されている。
【０１３５】
プログラム格納部６１には、ＯＳ（Ｏｐｅｒａｔｉｎｇ Ｓｙｓｔｅｍ）、マイクロリアクタ実験プログラムなどの各種プログラムがＣＰＵ５１で実行可能に記憶されている。
【０１３６】
ＯＳは、ファイル入出力の管理や各種機能部を制御するなど、測定装置１００を運営する基本的な機能をＣＰＵ５１に発揮させるためのプログラムである。
【０１３７】
マイクロリアクタ実験プログラムは、マイクロリアクタ１から出力されるデータを電子回路５７から取得して解離定数を予測したりなどするプログラムである。
【０１３８】
ＣＰＵ５１は、このプログラムを実行することにより図７で示した情報処理を行うことができる。
【０１３９】
データ格納部６２には、マイクロリアクタ１による実験データなどが記憶される。
【０１４０】
以上に説明したように、本実施の形態のマイクロリアクタは、基板に反応槽部と分析用のサンプル液を供給するサンプル液供給路と、緩衝液を供給する緩衝液供給路と、液体を外部へ排出するための廃液路とが設けられているマイクロリアクタにおいて、圧電基板上に設けられた少なくとも２対のすだれ状電極と、圧電基板上かつすだれ状電極間に設けられた特定物質を捕獲する物質捕獲手段とが反応槽部内に設けられている。そして、少なくとも１対のすだれ状電極間の電気容量を測定する測定手段や、２対のすだれ状電極間の電気容量を測定する測定手段や、２ついのすだれ状電極間で弾性波を発生させる電気的手段を備えることもできる。
【０１４１】
以上、本実施の形態について説明したが、これによって次のような効果が得られる。
（１）マイクロリアクタ１内に濃度センサ９を備えることができる。
（２）濃度センサ９をすだれ状の電極で構成することにより、小型で検出感度の高い濃度センサを実現することができる。
（３）質量センサ７でアナライトが検出された際のサンプル液の濃度を濃度センサ９で検出することができる。
（４）質量センサ７と濃度センサ９の出力をサンプリングし、これを所定の理論式に代入することにより、反応の初期段階で解離定数を高い精度で予測することができる。
（５）解離定数の予測値を用いて、研究資源を有望な新薬候補に集中することができる。
（６）サンプル液を廃液路５側から吸引してマイクロリアクタ１に供給し、マイクロリアクタ１を使い捨て方式とすることにより、コンタミネーションを防止することができる。
【０１４２】
なお、本実施の形態では、リガンドを異常蛋白質とし、アナライトを新薬候補の化合物としたが、逆にリガンドを新薬候補とし、アナライトを異常蛋白質としてもよい。
【０１４３】
また、質量センサ７に固定されたアナライトの量を求める方法は、共振周波数や弾性波の速度を用いたものに限定せず、例えば、表面プラズモン共鳴（ＳＰＲ）を用いることもできる。
【図面の簡単な説明】
【０１４４】
【図１】本実施の形態のマイクロリアクタの構造を模式的に表した平面図である。
【図２】質量センサの構造について説明するための図である。
【図３】濃度センサの構造について説明するための図である。
【図４】ＳＡＷタイプのマイクロリアクタの構成を模式的に示した図である。
【図５】予測に用いる理論式について説明するための図である。
【図６】測定装置がサンプリングした値をプロットしたグラフである。
【図７】測定装置が解離定数を予測する手順について説明するためのフローチャートである。
【図８】フローチャートの続きである。
【図９】本実施の形態に係る測定装置のハードウェア的な構成の一例を示したブロック図である。
【図１０】従来例を説明するための図である。
【符号の説明】
【０１４５】
１ マイクロリアクタ
３ 導入路
４ 反応槽
５ 廃液路
７ 質量センサ
８ 濃度センサ
１１ リガンド
１３ 金薄膜
１４ 水晶振動子
１６ アナライト
２１ 電極
２２ 電極
２５ 弾性波発生電極
【技術分野】
【０００１】
本発明は、マイクロリアクタ、解離定数測定装置、解離定数測定方法、及び解離定数測定プログラムに関し、例えば、異常蛋白質と結合する化合物の探索に用いるものに関する。
【背景技術】
【０００２】
近年、ヒトゲノム（ヒトの遺伝子情報）の解析が終了し、異常な遺伝子構造が生成する異常蛋白質と病気との関係解明が進められつつある。
【０００３】
この関係の解明により、新薬の開発手法が、開発者の薬剤や化合物に対する経験と勘を頼りに行う既存手法から、異常蛋白質に直接作用する化合物を探索して新薬とする手法へと変化している。この手法の採用により、２０年近く必要としていた新薬開発の期間が、今後５年程度に短縮すると見込まれている。
【０００４】
新薬候補の化合物の探索には、異常蛋白質と新薬候補の化合物との物理的な結合力を表す定数である解離定数を指標として用いるのが一般的である。
【０００５】
解離定数の測定方法としては、以前は、酵素、発光物質、放射性同位元素などの標識物質を結合させた化合物を用い、この化合物と異常蛋白を結合させた後、標識物質の量を測定することで、結合した化合物を定量し、解離定数の算出を行っていたが、現在では、標識物質を用いずにマイクロリアクタで行う測定が注目されている。
【０００６】
マイクロリアクタとは、半導体やガラス、樹脂などで構成されたチップの中に導入路と廃液路を形成し、その間に反応槽を設けたデバイスである。反応槽には、予め異常蛋白質を固定したセンサが設置される。
【０００７】
このように構成されたマイクロリアクタにおいて、導入路から化合物を含むサンプル液を流し込むと、サンプル液中の化合物が反応槽に予め設置された異常蛋白質と反応し、反応後の廃液が廃液路から排出される。
【０００８】
なお、反応槽で予め固定されている物質はリガンドと呼ばれ、溶液として供給される物質はアナライトと呼ばれるため、以下ではこの用語を用いる。
【０００９】
反応槽においては、サンプル液中のアナライトのうちあるものはリガンドと結合してセンサに固定される。一方、センサに固定されていたアナライトのうち、リガンドから離れてサンプル液中に放出されるものもある。
【００１０】
アナライトとリガンドの反応が平衡状態に達すると（即ち、センサに固定されるアナライトの量と、センサから離れるアナライトの量が等しくなると）、リガンドに固定されているアナライトの量が一定量となる。このとき計測されるアナライトとリガンドとが結合した化合物、結合していないリガンド、結合していないアナライトの個数（モル濃度）から解離定数が算出され、この解離定数が新薬の探索で必要なデータとなる。
【００１１】
このようなマイクロリアクタに関する技術としては、次の「マイクロリアクタ」がある。
【特許文献１】特表２００３−５２６３５９公報 この技術は、微小な混合室、加熱／分散ユニット、反応室、分離ユニットといった微小反応器をチップに並列に形成し、大量の生成物を生産するものである。
【００１２】
また、リガンドとアナライトの結合を計測する技術としては次の「免疫検査方法」がある。
【特許文献２】特許２９６２０３１ この技術は、リガンドとアナライトの結合を音響センサを用いて検出するものである。
【発明の開示】
【発明が解決しようとする課題】
【００１３】
従来は、解離定数を求めるためにはリガンドとアナライトの反応が平衡状態に達する必要があり、平衡状態に達するまでに長時間を要するため、解離定数を得るのに時間がかかるという問題があった。探索に係るリガンドとアナライトの組み合わせは膨大であるため、全ての組み合わせに対して反応が平衡に達するまで実験を継続するのはコスト的に問題があった。
【００１４】
また、長時間実験を行うことにより高価なサンプルを多量に消費するという問題もあった。
【００１５】
反応の初期の段階で、有望な組み合わせであるか否かをある程度予測できるようにしたいとの要望もあるが、正確な予測は困難であった。ここで、その理由について説明する。
【００１６】
図１０（ａ）は、センサに固定されるアナライトの質量の時間変化を表したグラフである。縦軸は、センサに固定されたアナライトの質量を表し、横軸は時間軸を表している。
【００１７】
まず、区間１０１では、緩衝液がマイクロリアクタに供給される。緩衝液にはアナライトが含まれていないため、センサに固定されるアナライトの質量は０である。このように、サンプル液に先立って緩衝液を流すのは、センサにリガンドとして固定されている異常蛋白質を乾燥させないためである。
【００１８】
マイクロリアクタに供給する溶液を緩衝液からサンプル液に切り替えると、区間１０２に示したようにセンサで固定されるアナライトの質量が徐々に増加してくる。そして、リガンドとアナライトの反応が平衡に達すると区間１０３で示したようにセンサで固定されるアナライトの質量は一定値となる。
【００１９】
ここで、例えば、図１０（ｂ）の時点１０５において、それまでの曲線の形からその後のアナライトの質量を曲線１０７のように予測したとする。しかし、実際に計測を続行すると例えば曲線１０８が得られ、予測値から大きく外れたものとなる。
【００２０】
このように予測が困難となる大きな原因としては、センサに供給する溶液を緩衝液からサンプル液への切り替える際に、センサでサンプルの濃度が徐々に変化するため、センサでのサンプル液の濃度が反応初期段階ではわからないことが挙げられる。
【００２１】
図１０（ｃ）はセンサでのサンプル濃度の時間変化を示したグラフである。曲線１１０に示したように、サンプルの濃度は時間と共に徐々に上昇する。これは、供給する溶液を緩衝液からサンプル液に切り替えた後、導入路内において両者が拡散混合したり、あるいは、流路内の流速分布により、重いアナライトが流れの速い流速域に集中したり、あるいは、センサに至るまでの導入路壁にアナライトが吸着するために起こるものである。
【００２２】
破線１１１に示したように、緩衝液からサンプル液にある時点で二者択一的に変化するのであれば、センサでのサンプルの濃度がわかるが、実際はセンサで質量変化が検出されても、その時点に至るまでのサンプルの濃度変化の経過がわからないため、予測が困難となる。
【００２３】
そこで、本発明の目的は、反応の初期段階において解離定数の予測を行えるようにすることである。
【課題を解決するための手段】
【００２４】
本発明は、前記目的を達成するために、試料を含む溶液を流す流路と、前記流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、を具備したことを特徴とするマイクロリアクタを提供する（第１の構成）。
【００２５】
第１の構成において、前記質量取得手段は、前記獲得物質を固定する共振部材の共振周波数を質量情報として取得するように構成することもできる（第２の構成）。
【００２６】
第１の構成において、前記質量取得手段は、前記獲得物質の上を流れる溶液を伝播する弾性波の速度を質量情報として取得するように構成することもできる（第３の構成）。
【００２７】
第１の構成、第２の構成、又は第３の構成において、前記濃度取得手段は、試料が獲得された後の溶液の誘電率情報を濃度情報として取得するように構成することもできる（第４の構成）。
【００２８】
また、本発明は、試料を含む溶液を流す流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、前記取得した質量情報と濃度情報を用いて、前記試料と前記獲得物質の解離定数を算出する解離定数算出手段と、を具備したことを特徴とする解離定数測定装置を提供する（第５の構成）。
【００２９】
また、本発明は、第１の構成から第４の構成までのうちの何れか１の構成のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求める解離定数算出方法であって、結合物質量算出手段と、未獲得量算出手段と、解離定数算出手段と、を備えたコンピュータにおいて、前記結合物質量算出手段によって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出ステップと、前記未獲得量算出手段によって、前記算出した結合物質の量を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出ステップと、前記解離定数算出手段によって、前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出ステップと、から構成されたことを特徴とする解離定数算出方法を提供する（第６の構成）。
【００３０】
また、本発明は、第１の構成から第４の構成までのうちの何れか１の構成のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求めるための解離定数算出プログラムであって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出機能と、前記算出した結合物質の量の時間変化を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出機能と、前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出機能と、をコンピュータで実現する解離定数算出プログラムを提供する（第７の構成）。
【発明の効果】
【００３１】
本発明によると、リガンドとアナライトの反応の初期段階において解離定数を予測することができる。
【発明を実施するための最良の形態】
【００３２】
（実施の形態の概要）
アナライトを獲得する質量センサの下流側近傍にサンプルの濃度センサを備える。これにより質量センサで質量が検出された際のサンプル液の濃度を知ることができる。
【００３３】
濃度と質量を後述の所定の理論式に代入すると解離定数の予測値を得ることができる。
【００３４】
濃度センサは、対向する電極からなるコンデンサで構成する。電極間にアナライトが存在すると誘電率が変化するため、電気容量を計測することによりアナライト濃度を検出することができる。
【００３５】
また、コンデンサは後述のすだれ構造の採用などにより小型化できるため、マイクロリアクタを小型に保つことができる。
（実施の形態の詳細）
図１は、本実施の形態のマイクロリアクタの構造を模式的に表した平面図である。
【００３６】
マイクロリアクタ１は、例えばガラスなどで構成された筐体部材により構成されており、この筐体部材内に各種の構成要素が形成されている。
【００３７】
マイクロリアクタ１の外寸は、長さ（導入路３の方向）が５０〜１００［ｍｍ］程度、幅が１０〜２０［ｍｍ］程度、厚さが数［ｍｍ］程度である。
【００３８】
マイクロリアクタ１には、導入路３、反応槽４、廃液路５などが形成されている。そして、反応槽４には質量センサ７が設置され、質量センサ７の廃液路５側には質量センサ７と近接して濃度センサ９が設けられている。
【００３９】
なお、緩衝液とサンプル液の流れを制御するバルブや、質量センサ７、濃度センサ９を構成する配線などは図示していない。
【００４０】
導入路３は、マイクロリアクタ１の外部から反応槽４に溶液を導く流路であり、緩衝液、及びアナライトを含んだサンプル液はこの管を通って反応槽４に導かれる。
【００４１】
一方、廃液路５は、反応槽４から排出される溶液をマイクロリアクタ１の外部に導く流路であり、反応槽４で反応を終えたサンプル液はこの管を通って排出される。
【００４２】
本実施の形態では、廃液路５の側から溶液を吸引することにより、導入路３から溶液が供給されるように構成した。
【００４３】
導入路３の側にポンプを設置して、溶液をマイクロリアクタ１に供給することも可能であるが、この場合、ポンプなどに前の実験で使用したアナライトが残留し、今回の測定に混入してしまうことで、測定精度に影響する場合があり（コンタミネーション）、本実施の形態では、溶液を吸引してマイクロリアクタ１に供給することにより、このような問題を防止した。
【００４４】
質量センサ７（質量取得手段）は、表面にリガンドが固定してあり、着脱可能に構成されている。以下に図２を用いて質量センサ７の構造について説明する。
【００４５】
図２は反応槽４の断面を表している。質量センサ７は、金薄膜１３、リガンド１１、水晶振動子１４などから構成されている。
【００４６】
金薄膜１３の上には、図示しない自己組織化膜を介してリガンド１１が予め固定されている。また、金薄膜１３の下面には水晶振動子１４が設けられている。水晶振動子１４は、コルピッツ発信回路の一部を成しており、水晶振動子１４を厚みすべり振動モードで共振させ、金薄膜１３をサンプルの流れと平行な面内で振動させるようになっている。ここで、金薄膜１３、リガンド１１、水晶振動子１４は共振部材を構成している。
【００４７】
質量センサ７の上を流れるアナライト１６のうち、あるものはリガンド１１と結合して（例えば、アナライト１６ａ）質量センサ７に固定される。
【００４８】
また、一旦質量センサ７に固定されたアナライトであっても、アナライト１６ｂのように、リガンドとの結合が切れて再びサンプル液中に放出されるものもある。
【００４９】
このように、リガンド１１は、試料（アナライト１６）を含む溶液から試料を獲得する獲得物質を構成している。
【００５０】
質量センサ７にアナライトが固定されると質量センサ７の質量が増すため、質量センサ７の共振周波数がｆｒからｆｒ＋Δｆへシフトする。このシフト量Δｆに次の式（１）を適用することにより質量センサ７に固定されたアナライトの質量Δｍを計算することができる。ここで、Ａは金薄膜１３の面積であり、ρqは水晶の密度であり、μqは水晶の剪断応力である。
【００５１】
なお、共振周波数は質量センサ７に固定されたアナライトの質量を表す質量情報を構成している。
【００５２】
Δｆ＝２ｆｒ2Δｍ／｛Ａ（ρqμq）1/2｝・・・（１）
次に、濃度センサ９（濃度取得手段）について説明する。
【００５３】
濃度センサ９は、図３（ａ）に示したように、電極２１と電極２２からなるコンデンサにより構成されている。図３（ｂ）は、電極２１、２２の断面を示している。
【００５４】
濃度センサ９は、質量センサ７と同一平面上に近接して配置されている。濃度変化は流路の高さ方向で異なるので、質量センサ７と同一平面に濃度センサ９を設けることで、解離定数の高精度予測が可能となる。
【００５５】
電極２１、２２は、廃液路５の底面に凸状に形成された電極部が廃液路５内で対向するように配置されて構成されている。
【００５６】
濃度センサ９の上をサンプル液が通過すると、対向した電極２１、２２の間にサンプル液中のアナライトが出入りし、このアナライトの量で電極２１、２２間の誘電率が変化する。
【００５７】
誘電率の変化と電気容量の変化は対応させることができるため、電極２１、２２間の電気容量によりサンプル液の濃度を計測することができる。ここで、電気容量は誘電率情報、及び濃度情報を構成している。
【００５８】
また、本実施の形態では、電極２１、電極２２の形状を１対のすだれ形状とすることにより、小さい設置領域面積で大きな電気容量の変化を得ることができるようにしてある。
【００５９】
即ち、すだれ状に電極を配置することにより、図３（ｂ）に示したように、電極２１ａと電極２２ａ、電極２２ａと電極２１ｂ、電極２１ｂと電極２２ｂの３つのコンデンサを形成することができる。
【００６０】
このように濃度センサ９を小型化できることから、濃度センサ９の設置によるマイクロリアクタ１の大型化を防ぐことができ、これによって、貴重な試料の使用量を抑えることができる。
【００６１】
以上のようにマイクロリアクタ１は、質量センサ７の下流側近傍に濃度センサ９を設けたため、質量センサ７でΔｍが検出された際のサンプル液の濃度を知ることができる。
【００６２】
なお、マイクロリアクタ１では、水晶振動子１４の共振周波数の遷移からΔｍを求めたが、弾性波の伝播速度によりΔｍを測定することも可能である。
【００６３】
一般に、共振周波数を用いるタイプの質量センサはＱＣＭ（Ｑｕａｒｔｚ Ｃｒｙｓｔａｌ Ｍｉｃｒｏｂａｌａｎｃｅ）と呼ばれ、弾性波を用いるタイプのものはＳＡＷ（Ｓｕｒｆａｃｅ Ａｃｏｕｓｔｉｃ Ｗａｖｅ）センサと呼ばれている。
【００６４】
図４は、ＳＡＷセンサを搭載したマイクロリアクタ１ａの構成を模式的に示した図である。
【００６５】
質量センサ７の上流側近傍には弾性波発生電極２５が設けられている。
【００６６】
弾性波発生電極２５の構成は、濃度センサ９と同じであり、１対の電極がすだれ状に配置されている。弾性波発生電極２５に電圧を印可すると、例えば、図示しない圧電素子が伸縮して弾性波が発生する。
【００６７】
濃度センサ９は、弾性波の検出を行うと共に、弾性波の検出の合間においてサンプル液の濃度の検出も行う。
【００６８】
このように、ＳＡＷセンサを搭載したマイクロリアクタでは、下流側のコンデンサに濃度検出手段と、弾性波検出手段を兼用させることができる。
【００６９】
弾性波発生電極２５と濃度センサ９の間の距離は予めわかっているため、弾性波発生電極２５が弾性波を発生してから濃度センサ９がこの弾性波を検出するまでの時間差により質量センサ７上を通過した弾性波の速度を求めることができる。
【００７０】
質量センサ７でアナライトが捕獲されると弾性波の伝達速度が遅くなることが知られており、所定の理論式により、弾性波の速度から質量センサ７で固定されたアナライトの量を計算することができる。
【００７１】
なお、ＳＡＷセンサを搭載したマイクロリアクタの場合、質量を検出するための領域（質量センサ７）は、対向する弾性波発生電極２５と濃度センサ９の幅で十分であり（これら電極の幅以外の領域はアナライトを結合しても検出できないため）、反応漕を円形にする必要はない。また、流路の幅を一定にすることができるため、この領域での圧力損失がなく、サンプル液の送液制御が容易になるという利点もある。
【００７２】
次に、マイクロリアクタ１を用いた解離定数の予測方法について説明する。
【００７３】
ここでは、まず、図５を用いて予測に用いる理論式について説明する。以下では、リガンドをＬで表し、アナライトをＡで表す。更にリガンドとアナライトが結合した化合物をＬＡで表す。なお、［］は、その内部に記された物質のモル濃度を示す。
【００７４】
まず、ＬＡ、Ｌ、Ａの間には図５の式（２）で表した関係が成り立つ。
【００７５】
即ち、リガンドとアナライトの結合により質量センサ７で固定されているＬＡのモル数［ＬＡ］の微小時間内での変化量は、質量センサ７でアナライトが固定されて生成されるＬＡのモル数（ｋｐ［Ｌ］［Ａ］）から、一旦捕らえられたアナライトが質量センサ７から放たれて減少するモル数（ｋｍ［ＬＡ］）を減算した値となる。平衡状態では、右辺の第１項と第２項は等しくなる。
【００７６】
なお、質量センサ７で、Ｌ＋Ａ→ＬＡの反応式により固定されるＬＡのモル数は、質量センサ７でまだアナライトと結合していないリガンドのモル数［Ｌ］とサンプル液中のアナライトの濃度（モル濃度）［Ａ］をかけたものに比例し、その比例定数をｋｐとしている。
【００７７】
また、質量センサ７で形成されたＬＡのうち、ＬＡ→Ｌ＋Ａの反応式により結合しているアナライトがサンプル液中に離脱するもののモル数は、質量センサ７で形成されているＬＡのモル数に比例し、その比例定数をｋｍとしている。
【００７８】
式（２）をｋｐ、ｋｍについて解くと、それぞれ式（３）、式（４）が得られる。
【００７９】
一方、式（２）を差分で表すと、式（５）が得られる。これは、微小時間ｄＴの間隔で、［ＬＡ］、［Ｌ］、［Ａ］をサンプリングした値を用いて式（２）を書き換えたものである。添え字のｉは、何番目のサンプリング値であるかを示している。更に、式（５）を［ＬＡ］iについて解くと式（６）が得られる。
【００８０】
ところで解離定数Ｋｄはｋｍ／ｋｐで表され、また平衡状態ではｄ［ＬＡ］／ｄｔ＝０であるので、平衡状態ｃでの解離定数Ｋｄは、式（７）で表される。
【００８１】
また、［Ｌ］c＝［Ｌ］0−［ＬＡ］c（［Ｌ］0はＬＡが形成される前のリガンドのモル数）となるので式（８）が導かれ、更に式（８）から式（９）が導かれる。
【００８２】
次に、データのサンプリングについて説明する。
【００８３】
後述する測定装置は、サンプリング周期ｄＴでマイクロリアクタ１から質量センサ７の共振周波数と濃度センサ９の電気容量を取得し、これらを解析して質量センサ７に固定されたアナライトの質量Δｍと、その時のサンプル液中のアナライトのモル濃度［Ａ］を算出する。
【００８４】
これにより、ｄＴ間隔でΔｍと［Ａ］を得ることができる。
【００８５】
測定装置がサンプリングしたこれらの値をプロットしたのが図６（ａ）、（ｂ）である。
【００８６】
図６（ａ）は、質量センサ７で検出されたアナライトの質量を表しており、図６（ｂ）は、濃度センサ９で検出されたサンプル液中のアナライトのモル濃度を表している。両グラフの時間軸は対応させてある。
【００８７】
図６に示したように、測定装置により、時刻Ｔ0にてΔｍ0と［Ａ］0が得られ、時刻Ｔ1にてΔｍ1と［Ａ］1が得られ、以下同様にして、ｄＴごとにΔｍと［Ａ］が検出される。
【００８８】
次に、図７のフローチャートを用いて、測定装置が解離定数を予測する手順について説明する。
【００８９】
まず、測定装置は初期値の設定を行う（ステップ５）。この初期化では、質量センサ７に予め固定されているリガンドのモル数Ｌｅを［Ｌ］0、マイクロリアクタ導入口でのアナライトのモル濃度Ａｃを［Ａ］ｃに設定し、更にアナライトの分子量Ｍａ、リガンドの分子量ＭＬを設定する。これらの設定は、ユーザが測定装置に入力することにより行われる。また、測定装置は自動的に［ＬＡ］を０に設定する。
【００９０】
次に、測定装置は質量センサ７の出力よりΔｍ0を測定し、更に濃度センサ９の出力より［Ａ］0を測定し、これらの値を記憶装置に記憶する（ステップ１０）。
【００９１】
測定装置は、ｄＴ間待った後（ステップ１５）、Δｍ1と［Ａ］1を測定し、記憶する（ステップ２０）。
【００９２】
次に、測定装置は、記憶したこれらの値を用いて［ＬＡ］1と［Ｌ］1を計算する（ステップ２５）。
【００９３】
［ＬＡ］1（質量センサ７で形成されたＬＡのモル数）は、質量センサ７で検出された質量の増分（Δｍ1−Δｍ0）をアナライトの分子量Ｍａで除すれば得られる。
【００９４】
また、［Ｌ］1（質量センサ７でアナライトと結合していないリガンドのモル数）は、リガンドの初期値［Ｌ］0から生成されたＬＡのモル数（［ＬＡ］1−［ＬＡ］0）を減算すれば得られる。
【００９５】
次に、測定装置は、ここで得られている［ＬＡ］1、［Ｌ］1、［Ａ］1を用いてｋｐ1の仮定値を計算して記憶する（ステップ３０）。計算式はフローチャートに示したとおりである。この値は式（３）においてｋｍ［ＬＡ］＝０とおいたものである。リガンドとアナライトの反応が始まった初期段階であるので、質量センサ７から放たれるアナライトの数を０としたのである。
【００９６】
次に、測定装置は、ｄＴ間待った後（ステップ３５）、Δｍ2と［Ａ］2を測定して記憶する（ステップ４０）。
【００９７】
そして、測定装置は、ステップ２５と同様にして、［ＬＡ］2と［Ｌ］2を計算して記憶する（ステップ４５）。
【００９８】
次に、測定装置は、ｋｍ2とｋｐ2の仮定値を計算して記憶する（ステップ５０）。計算式はフローチャートに示したとおりである。ｋｍ2の計算式は、式（５）を用いたものである。ｋｐ2は、ｄＴ間ではｋｐ1とほとんど差がないと考えられるのでｋｐ1と等しいとした。
【００９９】
次に、測定装置は、ｋｐ2、ｋｍ2などを用いて［ＬＡ］3の予測値を計算する（ステップ５５）。計算式はフローチャートに示したとおりである。この計算式は式（６）から得られたものである。以下、図８のフローチャートに続く。
【０１００】
次に、測定装置は、カウンタｉを３にセットし、ｉが所定の整数Ｎに達するまで以下のループ処理を行う（ステップ６０）。なお、以下の処理では、仮定値を用いて解離定数Ｋｄを計算し、その評価を行うが、解離定数Ｋｄが許容誤差の範囲に入った場合は、ループを抜けてこの解離定数Ｋｄを予測値として出力する。また、整数Ｎは、解離定数Ｋｄの予測値が求まるまでループ処理が継続するように十分に大きい値とする。
【０１０１】
次に、測定装置は、ｄＴ待った後（ステップ６５）、Δｍiと［Ａ］iを測定して記憶する（ステップ７０）。
【０１０２】
次に、測定装置は、ステップ２５と同様にして［ＬＡ］iと［Ｌ］iを計算して記憶装置に記憶する（ステップ７５）。
【０１０３】
次に、測定装置は［ＬＡ］iの予測値の評価を行う（ステップ８０）。評価は、ステップ７５で計算された［ＬＡ］iと、［ＬＡ］iの予測値（ｉ＝３の場合はステップ５５で求めた値、その他の場合はステップ１１０で求めた値）の差を［ＬＡ］iで除した値が予め設定されている誤差許容値未満であるか否かを判定することにより行う。
【０１０４】
［ＬＡ］iの計算値と予測値の差が誤差許容値未満であった場合（ステップ８０；Ｙ）、測定装置は、ｉを１だけデクリメントし（ステップ８５）、解離定数Ｋｄと化学平衡時の［ＬＡ］cを計算して出力し（ステップ１０５）、処理を終了する。Ｋｄと［ＬＡ］cの計算式はフローチャートに示したとおりである。なお、［ＬＡ］cの計算式は式（９）による。
【０１０５】
なお、式（９）は、リガンドとアナライトの反応が平衡している場合に成り立つ式（７）から導出した関係式であるが、最終的には、［Ａ］ｃがマイクロリアクタ導入口に流し込むサンプル液濃度Ａｃになると考えたことによる。
【０１０６】
一方、［ＬＡ］iの計算値と予測値の差が誤差許容値未満でなかった場合（ステップ８０；Ｎ）、測定装置は、ｋｍiとｋｐiの仮定値を計算し、記憶装置に記憶する。計算式はフローチャートに示したとおりである。これらの計算式は何れも式（５）を用いたものであり、現在記憶している最新の値であるｋｍi-1、ｋｐi-1を用いて仮定値を更新するようになっている。
【０１０７】
次に、計測装置は、ステップ９０で仮定したｋｍiとｋｐiを用いて、ステップ５５と同様にして［ＬＡ］iの修正予測値ＮＥ［ＬＡ］iを計算して記憶装置に記憶する（ステップ９５）。計算式はフローチャートに示したとおりである。
【０１０８】
次に、計測装置は、ステップ８０と同様にして修正予測値ＮＥ［ＬＡ］iの評価を行う（ステップ１００）。評価式はフローチャートに示したとおりである。
【０１０９】
［ＬＡ］iの計算値と修正予測値ＮＥ［ＬＡ］iの差が予め設定した誤差許容値未満の場合は（ステップ１００；Ｙ）、ステップ１０５にて解離定数Ｋｄと化学平衡時の［ＬＡ］cを算出して出力する（ステップ１０５）。
【０１１０】
［ＬＡ］iの計算値と修正予測値ＮＥ［ＬＡ］iの差が予め設定した誤差許容値未満でなかった場合は（ステップ１００；Ｎ）、ステップ５５と同様にして次の予測値Ｅ［ＬＡ］i+1を計算する（ステップ１１０）。計算式はフローチャートに示したとおりである。
【０１１１】
そして、測定装置は、カウンタｉを１だけインクリメントし（ステップ１１５）、ステップ６０に戻る。そして、同様の処理を繰り返すことになる。
【０１１２】
このように、測定装置は、マイクロリアクタ１より獲得物質（リガンド）によって獲得された溶液中の試料（アナライト）の質量を表す質量情報（共振周波数）を取得する質量取得手段と、前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報（電気容量）を取得する濃度取得手段と、前記取得した質量情報と濃度情報を用いて、前記試料と前記獲得物質の解離定数を算出する解離定数算出手段と、を具備している。
【０１１３】
また、測定装置は、結合物質量算出手段と、未獲得量算出手段と、解離定数算出手段と、を備えたコンピュータで構成することができ、当該測定装置による測定方法は、前記結合物質量算出手段によって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量（［ＬＡ］c）を算出する結合物質量算出ステップと、前記未獲得量算出手段によって、前記算出した結合物質の量を用いて、前記試料を獲得していない獲得物質の量（［Ｌ］c＝［Ｌ］0−［ＬＡ］c）を算出する未獲得量算出ステップと、前記解離定数算出手段によって、前記濃度取得手段で取得した濃度情報（［Ａ］）、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出ステップと、を備えている。
【０１１４】
以上に説明した手順により、測定装置は、マイクロリアクタ１から出力される質量情報と濃度情報をサンプリングすることにより、リガンドとアナライトが反応を始める初期の段階で、解離定数を予測することができる。そのため、ユーザは、リガンドとアナライトの反応の初期の段階で、そのリガンドとアナライトの組み合わせが新薬候補としてどの程度有望か見当を付けることができる。
【０１１５】
これによって、例えば、新薬候補となる見込みがある場合は、リガンドとアナライトの結合が飽和するまで実験を継続し、見込みがない場合は次の組み合わせの実験に移行するといった運用が可能となる。
【０１１６】
図９は、本実施の形態に係る測定装置のハードウェア的な構成の一例を示したブロック図である。
【０１１７】
測定装置１００は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）５１、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）５２、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）５３、入力装置５４、表示装置５５、印刷装置５６、電子回路５７、記憶装置５８、記憶媒体駆動装置５９、入出力Ｉ／Ｆ（インターフェース）６０などの各機能部がバスライン５０で接続されて構成されている。
【０１１８】
ＣＰＵ５１は、所定のプログラムに従って、各種の演算処理、情報処理、及び測定装置１００を構成する各構成要素の制御を行ったりする中央処理装置である。
【０１１９】
ＣＰＵ５１は、電子回路５７を介してマイクロリアクタ１から出力される共振周波数及び電気容量をサンプリングしたり、サンプリングしたこれらの値を用いて質量センサ７で固定されたアナライトの質量やサンプリング液の濃度を算出するほか、これらの値を用いて解離定数Ｋｄの予測値の計算も行う。
【０１２０】
ＲＯＭ５３は、測定装置１００を動作させるための基本的なプログラムやデータなどを記憶した読み出し専用の記憶装置である。
【０１２１】
ＲＡＭ５２は、ＣＰＵ５１が動作するためのワーキングエリアを提供する読み書き可能な記憶装置である。ＲＡＭ５２は、ＣＰＵ５１が図７で示した情報処理を行う際には、Δｍ０、Δｍ１、・・・など、予測に必要な値を記憶する。
【０１２２】
入力装置５４は、例えば、キーボードやマウスなどの入力装置を備えており、ユーザが測定装置１００を使用するのに必要な情報を入力できるようになっている。図７のステップ５で示した値は入力装置５４から入力することができる。
【０１２３】
表示装置５５は、文字情報や画像情報を表示する表示デバイスを備えており、ユーザが測定装置１００を使用するのに必要なメニュー画面などが表示されるようになっている。また、予測した解離定数Ｋｄも表示装置５５に表示することができる。
【０１２４】
表示デバイスは、例えば、液晶ディスプレイ、プラズマディスプレイ、その他のディスプレイ装置で構成することができる。
【０１２５】
電子回路５７は、マイクロリアクタ１に電力を供給して駆動すると共に、マイクロリアクタ１から出力される質量センサ７の共振周波数や濃度センサ９の電気容量を検出する。
【０１２６】
更に、電子回路５７は、マイクロリアクタ１の出力データ（アナログデータ）を時間間隔ｄＴにてサンプリングしてデジタルデータに変換し、ＣＰＵ５１に提供する。
【０１２７】
ＣＰＵ５１は、サンプリングされたデータからΔｍと［Ａ］を計算する。これによって時間間隔ｄＴでΔｍと［Ａ］が得られる。
【０１２８】
印刷装置５６は、例えば、レーザプリンタやインクジェットプリンタなどの印刷装置であり、実験データを印刷できるようになっている。
【０１２９】
記憶媒体駆動装置５９は、装着された着脱可能な記憶媒体を駆動し、データの読み書きを行う機能部である。
【０１３０】
読み書き可能な記憶媒体としては、例えば、フレキシブルディスク、光磁気ディスク、半導体記憶装置、磁気テープ、紙テープがある。
【０１３１】
また、読み取り専用の記憶媒体としては、例えば、ＣＤ−ＲＯＭなどの光ディスク類がある。
【０１３２】
測定装置１００は、記憶媒体駆動装置５９に装着した記憶媒体からプログラムのインストールなどを行うことができるほか、実験データを記憶媒体に書き込むこともできる。
【０１３３】
記憶装置５８は、例えば、ハードディスクなどで構成された大容量で読み書き可能な記憶装置である。
【０１３４】
記憶装置５８には、プログラム類を格納したプログラム格納部６１とデータ類を記憶したデータ格納部６２が形成されている。
【０１３５】
プログラム格納部６１には、ＯＳ（Ｏｐｅｒａｔｉｎｇ Ｓｙｓｔｅｍ）、マイクロリアクタ実験プログラムなどの各種プログラムがＣＰＵ５１で実行可能に記憶されている。
【０１３６】
ＯＳは、ファイル入出力の管理や各種機能部を制御するなど、測定装置１００を運営する基本的な機能をＣＰＵ５１に発揮させるためのプログラムである。
【０１３７】
マイクロリアクタ実験プログラムは、マイクロリアクタ１から出力されるデータを電子回路５７から取得して解離定数を予測したりなどするプログラムである。
【０１３８】
ＣＰＵ５１は、このプログラムを実行することにより図７で示した情報処理を行うことができる。
【０１３９】
データ格納部６２には、マイクロリアクタ１による実験データなどが記憶される。
【０１４０】
以上に説明したように、本実施の形態のマイクロリアクタは、基板に反応槽部と分析用のサンプル液を供給するサンプル液供給路と、緩衝液を供給する緩衝液供給路と、液体を外部へ排出するための廃液路とが設けられているマイクロリアクタにおいて、圧電基板上に設けられた少なくとも２対のすだれ状電極と、圧電基板上かつすだれ状電極間に設けられた特定物質を捕獲する物質捕獲手段とが反応槽部内に設けられている。そして、少なくとも１対のすだれ状電極間の電気容量を測定する測定手段や、２対のすだれ状電極間の電気容量を測定する測定手段や、２ついのすだれ状電極間で弾性波を発生させる電気的手段を備えることもできる。
【０１４１】
以上、本実施の形態について説明したが、これによって次のような効果が得られる。
（１）マイクロリアクタ１内に濃度センサ９を備えることができる。
（２）濃度センサ９をすだれ状の電極で構成することにより、小型で検出感度の高い濃度センサを実現することができる。
（３）質量センサ７でアナライトが検出された際のサンプル液の濃度を濃度センサ９で検出することができる。
（４）質量センサ７と濃度センサ９の出力をサンプリングし、これを所定の理論式に代入することにより、反応の初期段階で解離定数を高い精度で予測することができる。
（５）解離定数の予測値を用いて、研究資源を有望な新薬候補に集中することができる。
（６）サンプル液を廃液路５側から吸引してマイクロリアクタ１に供給し、マイクロリアクタ１を使い捨て方式とすることにより、コンタミネーションを防止することができる。
【０１４２】
なお、本実施の形態では、リガンドを異常蛋白質とし、アナライトを新薬候補の化合物としたが、逆にリガンドを新薬候補とし、アナライトを異常蛋白質としてもよい。
【０１４３】
また、質量センサ７に固定されたアナライトの量を求める方法は、共振周波数や弾性波の速度を用いたものに限定せず、例えば、表面プラズモン共鳴（ＳＰＲ）を用いることもできる。
【図面の簡単な説明】
【０１４４】
【図１】本実施の形態のマイクロリアクタの構造を模式的に表した平面図である。
【図２】質量センサの構造について説明するための図である。
【図３】濃度センサの構造について説明するための図である。
【図４】ＳＡＷタイプのマイクロリアクタの構成を模式的に示した図である。
【図５】予測に用いる理論式について説明するための図である。
【図６】測定装置がサンプリングした値をプロットしたグラフである。
【図７】測定装置が解離定数を予測する手順について説明するためのフローチャートである。
【図８】フローチャートの続きである。
【図９】本実施の形態に係る測定装置のハードウェア的な構成の一例を示したブロック図である。
【図１０】従来例を説明するための図である。
【符号の説明】
【０１４５】
１ マイクロリアクタ
３ 導入路
４ 反応槽
５ 廃液路
７ 質量センサ
８ 濃度センサ
１１ リガンド
１３ 金薄膜
１４ 水晶振動子
１６ アナライト
２１ 電極
２２ 電極
２５ 弾性波発生電極
【特許請求の範囲】
【請求項１】
試料を含む溶液を流す流路と、
前記流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、
前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、
を具備したことを特徴とするマイクロリアクタ。
【請求項２】
前記質量取得手段は、前記獲得物質を固定する共振部材の共振周波数を質量情報として取得することを特徴とする請求項１に記載のマイクロリアクタ。
【請求項３】
前記質量取得手段は、前記獲得物質の上を流れる溶液を伝播する弾性波の速度を質量情報として取得することを特徴とする請求項１に記載のマイクロリアクタ。
【請求項４】
前記濃度取得手段は、試料が獲得された後の溶液の誘電率情報を濃度情報として取得することを特徴とする請求項１、請求項２、又は請求項３に記載のマイクロリアクタ。
【請求項５】
試料を含む溶液を流す流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、
前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、
前記取得した質量情報と濃度情報を用いて、前記試料と前記獲得物質の解離定数を算出する解離定数算出手段と、
を具備したことを特徴とする解離定数測定装置。
【請求項６】
請求項１から請求項４までのうちの何れか１の請求項に記載のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求める解離定数算出方法であって、
結合物質量算出手段と、未獲得量算出手段と、解離定数算出手段と、を備えたコンピュータにおいて、
前記結合物質量算出手段によって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出ステップと、
前記未獲得量算出手段によって、前記算出した結合物質の量を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出ステップと、
前記解離定数算出手段によって、前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出ステップと、
から構成されたことを特徴とする解離定数算出方法。
【請求項７】
請求項１から請求項４までのうちの何れか１の請求項に記載のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求めるための解離定数算出プログラムであって、
前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出機能と、
前記算出した結合物質の量の時間変化を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出機能と、
前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出機能と、
をコンピュータで実現する解離定数算出プログラム。
【請求項１】
試料を含む溶液を流す流路と、
前記流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、
前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、
を具備したことを特徴とするマイクロリアクタ。
【請求項２】
前記質量取得手段は、前記獲得物質を固定する共振部材の共振周波数を質量情報として取得することを特徴とする請求項１に記載のマイクロリアクタ。
【請求項３】
前記質量取得手段は、前記獲得物質の上を流れる溶液を伝播する弾性波の速度を質量情報として取得することを特徴とする請求項１に記載のマイクロリアクタ。
【請求項４】
前記濃度取得手段は、試料が獲得された後の溶液の誘電率情報を濃度情報として取得することを特徴とする請求項１、請求項２、又は請求項３に記載のマイクロリアクタ。
【請求項５】
試料を含む溶液を流す流路上に設けられた獲得物質によって獲得された前記溶液中の試料の質量を表す質量情報を取得する質量取得手段と、
前記獲得物質によって試料が獲得された際における、前記溶液中の前記試料の濃度を表す濃度情報を取得する濃度取得手段と、
前記取得した質量情報と濃度情報を用いて、前記試料と前記獲得物質の解離定数を算出する解離定数算出手段と、
を具備したことを特徴とする解離定数測定装置。
【請求項６】
請求項１から請求項４までのうちの何れか１の請求項に記載のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求める解離定数算出方法であって、
結合物質量算出手段と、未獲得量算出手段と、解離定数算出手段と、を備えたコンピュータにおいて、
前記結合物質量算出手段によって、前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出ステップと、
前記未獲得量算出手段によって、前記算出した結合物質の量を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出ステップと、
前記解離定数算出手段によって、前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出ステップと、
から構成されたことを特徴とする解離定数算出方法。
【請求項７】
請求項１から請求項４までのうちの何れか１の請求項に記載のマイクロリアクタを用いて前記試料の前記獲得物質に対する解離定数を求めるための解離定数算出プログラムであって、
前記質量取得手段で取得された質量情報を用いて、前記試料と前記獲得物質が結合した結合物質の量を算出する結合物質量算出機能と、
前記算出した結合物質の量の時間変化を用いて、前記試料を獲得していない獲得物質の量を算出する未獲得量算出機能と、
前記濃度取得手段で取得した濃度情報、前記算出した結合物質の量、及び前記算出した獲得物質の量と、を用いて前記試料と前記獲得物質の解離定数を算出する解離定数算出機能と、
をコンピュータで実現する解離定数算出プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【公開番号】特開２００６−２１４７９２（Ｐ２００６−２１４７９２Ａ）
【公開日】平成１８年８月１７日（２００６．８．１７）
【国際特許分類】
【出願番号】特願２００５−２６０３０（Ｐ２００５−２６０３０）
【出願日】平成１７年２月２日（２００５．２．２）
【出願人】（０００００２３２５）セイコーインスツル株式会社 (3,629)
【Ｆターム（参考）】
【公開日】平成１８年８月１７日（２００６．８．１７）
【国際特許分類】
【出願日】平成１７年２月２日（２００５．２．２）
【出願人】（０００００２３２５）セイコーインスツル株式会社 (3,629)
【Ｆターム（参考）】
[ Back to top ]