JP4710021B2 - 生体パラメータ出力装置およびプログラム - Google Patents

Description

本発明は、細胞膜の活動電位波形の変化を入力として与えることにより、活動電位シミュレーションを用いて、薬物投入等が効果等を及ぼす生体パラメータを推定する装置である生体パラメータ出力装置等に関するものである。

現在、新規薬物の試験として心電図ＱＴ延長等の評価を始めとして様々な評価を行う必要がある。製薬企業等は培養細胞試験、動物実験、ヒトを対象にした試験等を行い、それら新規薬物が生体の様々なパラメータに対して及ぼす変化を試験し、効果及び安全性等を確認している。

また、薬物の生体組織内での拡散現象を的確に解析することができる薬物の生体組織内拡散のシミュレーション方法及び装置がある（特許文献１参照）。かかるシミュレーション方法は、特定物質の生体内拡散を有限要素法を用いてシミュレーションする方法であって、生体から分離した組織の一部を用いて、有限要素法を使用することなく該生体内の基準拡散特性定数を決定し、該基準拡散特性定数に基づいて基準拡散特性を設定し、解析の対象となる生体構造を有限要素法に基づいて決定し、前記基準拡散特性定数を用いて前記決定した生体構造における拡散を有限要素法に基づいて演算し、該生体内における拡散の有限要素法に基づく演算結果にかかる演算拡散特性と有限要素法によらないで決定した前記基準拡散特性とを比較し、前記演算拡散特性と前記基準拡散特性との偏差が最小となるように基準拡散特性定数を補正して有限要素法に基づく最適拡散特性係数を算出することを特徴とする生体内拡散のシミュレーション方法である。

また、動脈硬化の段階に応じた複数の動脈圧波形の基準情報を記憶しておき、生体の動脈から得られた動脈圧波形情報と比較し、類似している波形の年齢を出力する動脈硬化度測定方法がある（特許文献２、特許文献３参照）。
さらに、関連する技術として、生体パラメータセットを入力として受け取り、細胞をシミュレーションし、活動電位波形情報を取得するシミュレーション装置がある（非特許文献１参照）。
特開平０８−０１６５５１（第１頁、第１図等） 特開昭６１−１１９２５２（第１頁、第１図等） 特開平０３−０１５４３９（第１頁、第１図等） 皿井伸明、野間昭典「ｓｉｍＢｉｏ：生物学的ダイナミックモデル開発基盤」日本エム・イー学会雑誌ＢＭＥ，ｖｏｌ．１８，Ｎｏ．２，ｐ．３−１１，２００４（２００４年２月発行）

しかしながら、上記の試験は、非常に膨大な評価試験であって、時間とコストが膨大にかかってしまう。
また、特許文献１におけるシミュレーション方法は、薬物の生体組織内での拡散現象を解析する方法であり、薬物投入等の効果等を示す生体パラメータを推定する方法ではない。

さらに、特許文献２、３における技術は、動脈硬化の度合いを測定する技術であり、投入薬物が生体パラメータに与える影響を評価すること等はできない。また、特許文献２及び３における年齢決定手段は予め基準パターンを該年齢に応じて複数記憶しておかねばならない。

これらの膨大な評価試験において、より早期の段階のデータを元に、対象薬物の各種チャネル機能への関与、心電図変化の推定を計算機上でシミュレートし、人に対する効果を推測出来れば、創薬プロセスが短縮されることが期待される。そこで、本発明は、培養細胞、もしくは動物実験等により、薬物投与時の活動電位波形変化を測定し、これを用いて、細胞の個々のチャネルに対する投入薬物の作用を評価し、投入薬物が生体パラメータに与える影響を評価する手法を提供することを目的とする。

本第一の発明は、生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を１以上含む生体パラメータセットとそれに対応する活動電位波形を示す情報である活動電位波形情報を有する波形情報を１以上格納している波形情報格納部と、活動電位波形情報の入力を受け付ける波形情報受付部と、前記波形情報受付部が受け付けた活動電位波形情報に最も近い波形情報を再現しうる１以上の生体パラメータ情報を含む生体パラメータセットを識別し、取得する生体パラメータ情報取得部と、前記生体パラメータ情報取得部が取得した生体パラメータセットを出力する生体パラメータ情報出力部を具備する生体パラメータ出力装置である。

また、第二の発明は、第一の発明に対して、前記波形情報格納部に格納されている波形情報のうちの生体パラメータセットは入力された情報であり、前記波形情報格納部に格納されている波形情報のうちの活動電位波形情報は、前記入力された生体パラメータセットに生体シミュレーションを用いて取得した情報である。

さらに、第三の発明は、第一、第二の発明に対して、前記活動電位波形情報は、ＡＰＤ３０、ＡＰＤ６０及びＡＰＤ９０の活動電位波形を特徴づける情報であり、生体パラメータ情報取得部は、前記活動電位波形情報のＡＰＤ３０、ＡＰＤ６０およびＡＰＤ９０の活動電位波形情報を取得し、当該取得した３つの活動電位波形情報と一致するまたは最も近似する１以上の生体パラメータ情報を含む生体パラメータセットを取得する生体パラメータ出力装置である。
かかる構成により、高速かつ少ないＣＰＵパワーで生体パラメータセットを取得することができる。

さらに、第四の発明は、かかる第一、第二の発明の構成により、生体シミュレーションを用いて、薬物投入前後の活動電位波形情報に最も近い波形情報を再現しうる生体パラメータセットを取得でき、薬物投入の効果を生体パラメータとして出力できる。

本発明は、薬物投入等の効果等を示す生体パラメータを出力する装置等を提供できる。

本実施の形態において、薬物投入等の効果を示したり、現在の生体の健康状態等を判断したりするために、生体パラメータを出力する装置等について述べる。なお、生体パラメータには細胞の種々のチャネル（たとえば、Ｎａチャネル、Ｃａチャネル、Ｋ_ＡＴＰチャネル、Ｋｒチャネル、Ｋ１チャネル、Ｋｓチャネルなど）を流れる電流や、各チャネルの開閉速度や、イオン親和性、細胞内外のイオン濃度など、数百にも及ぶパラメータがある。
以下、生体パラメータ出力装置等の実施形態について図面を参照して説明する。なお、実施の形態において同じ符号を付した構成要素は同様の動作を行うので、再度の説明を省略する場合がある。

図１は、生体パラメータ出力装置のブロック図である。生体パラメータ出力装置は、波形情報格納部１１、波形情報受付部１２、生体パラメータ情報取得部１３、生体パラメータ情報出力部１４を具備する。生体パラメータ情報取得部１３は、波形情報比較手段１３１、生体パラメータ情報生成手段１３２を具備する。

波形情報格納部１１は、生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を１以上含む生体パラメータセットと、および活動電位波形を示す情報である活動電位波形情報を有する波形情報を１以上格納している。

なお、波形情報は、通常、図示しない手段により、以下のように取得されたものである。つまり、波形情報を取得する装置等（図示しない）は、生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を１以上含む生体パラメータセットを入力として受け取り、それに対応する活動電位波形情報を取得する。これには、生体を用いて活動電位波形情報を取得しても良い。また、波形情報を取得する装置が、細胞をシミュレーションする細胞モデルに生体パラメータセットを適用し、シミュレーションを用い、活動電位波形情報（以下、適宜、シミュレーションを用いて取得した活動電位波形情報を「再構築活動電位波形情報」という。）を取得しても良い。この細胞モデルは、最も良く実験結果を再現しうるモデルに交換することが出来る。このとき、生体パラメータ出力装置は、波形情報を取得する装置を内蔵しており、シミュレーション結果をその都度生成しても良い。また、生体パラメータ出力装置は、あらかじめ生成された活動電位波形情報をデータベースより読み出しても良い。さらに、活動電位波形情報は、グラフの形式でも良いし、レコードの形式でも良い、つまり、活動電位波形情報のデータ構造は問わない。波形情報格納部１１は、不揮発性の記録媒体が好適であるが、揮発性の記録媒体でも実現可能である。なお、生体パラメータセットを１以上、およびそれに対応する活動電位波形情報を１以上格納している記憶媒体を波形情報データベース（図示しない）とする。なお、生体パラメータセットを入力として受け取り、細胞をシミュレーションし活動電位波形情報を取得するシミュレーション装置（波形情報を取得する装置の一例）については公知技術である（上記、非特許文献１参照）。

波形情報受付部１２は、活動電位波形に関する情報である活動電位波形情報の入力を受け付ける。なお、活動電位波形情報は、通常、個体、臓器、組織、細胞などの生体から測定した情報である。活動電位波形情報は、例えば、薬物投与前後における細胞の活動電位波形の情報や、遺伝子ノックアウト動物や疾患モデル動物の細胞から得られた活動電位波形等の情報である。活動電位波形情報は、膜電位の時系列データ（波形を構成する点の情報の集合）すべてでも良いし、すべてでなくても良い。つまり、活動電位波形情報は、波形のうち、ＡＰＤ３０の値、ＡＰＤ６０値、およびＡＰＤ９０の値等でも良い。ＡＰＤ３０の値、ＡＰＤ６０値、およびＡＰＤ９０の値とは、図２に示す値である。つまり、図２において、波形の最高電位（ａ点）から最低電位（ｂ点）の高さを１００とした場合に、ａ点から３０下の高さにおける、波形の幅がＡＰＤ３０である。同様に、波形の最高電位（ａ点）から最低電位（ｂ点）の高さを１００とした場合に、ａ点から６０下、９０下の高さにおける、波形の幅がＡＰＤ６０、ＡＰＤ９０である。活動電位波形情報の入力手段は、スキャナやキーボードやマウスやメニュー画面によるもの等、何でも良い。なお、活動電位波形情報は、他の装置、例えば、測定装置から出力された情報でも良い。波形情報受付部１２は、スキャナやキーボード等の入力手段のデバイスドライバーや、メニュー画面の制御ソフトウェア等で実現され得る。

生体パラメータ情報取得部１３は、波形情報受付部１２が受け付けた活動電位波形情報および波形情報格納部１１の波形情報の一部または全部に基づいて、１以上の生体パラメータ情報を含む生体パラメータセットを識別し、取得する。具体的な生体パラメータ情報の取得アルゴリズムの例は、後述する。生体パラメータ情報の取得には、実験計画法、遺伝的アルゴリズム、最急降下法、応答局面法等の種々のパラメータ探索アルゴリズムが適用できる。あらゆるシミュレーションの結果と比較し、パラメータを抽出することが出来る。生体パラメータ情報取得部１３は、通常、ＭＰＵやメモリ等から実現され得る。生体パラメータ情報取得部１３の処理手順は、通常、ソフトウェアで実現され、当該ソフトウェアはＲＯＭ等の記録媒体に記録されている。但し、ハードウェア（専用回路）で実現しても良い。

生体パラメータ情報出力部１４は、生体パラメータ情報取得部１３が取得した生体パラメータセットを出力する。出力とは、ディスプレイへの表示、プリンタへの印字、音出力、外部の装置への送信等を含む概念である。生体パラメータ情報出力部１４は、ディスプレイやスピーカー等の出力デバイスを含むと考えても含まないと考えても良い。出力部は、出力デバイスのドライバーソフトまたは、出力デバイスのドライバーソフトと出力デバイス等で実現され得る。

波形情報比較手段１３１は、波形情報格納部１１から取得した活動電位波形情報と波形情報受付部１２が受け付けた活動電位波形情報の近似度を判定する。波形情報比較手段１３１は、例えば、以下のようにＡＰＤ３０等の活動電位特徴情報を取得し、比較しても良い。つまり、波形情報受付部１２が受け付けた活動電位波形情報が波形を形成する電位（ｍＶ）と時間（ｍｓ）の数字列の集合である、とする。その場合、波形情報比較手段１３１は、「最高の電位−（最高の電位−最低の電位）×０．３」の値と対になる２つの時間（Ｔ_１、Ｔ_２）を取得する。次に、波形情報比較手段１３１は、「｜Ｔ_１−Ｔ_２｜」（Ｔ_１とＴ_２の差の絶対値）をＡＰＤ３０として求める。波形情報比較手段１３１は、ＡＰＤ６０、ＡＰＤ９０についても、同様の方法で求める。また、波形情報受付部１２が受け付けた活動電位波形情報は、ＡＰＤ３０、ＡＰＤ６０、およびＡＰＤ９０の値である、とする。そして、波形情報比較手段１３１は、算出したＡＰＤ３０と波形情報受付部１２が受け付けたＡＰＤ３０の差の絶対値、および算出したＡＰＤ６０と波形情報受付部１２が受け付けたＡＰＤ６０の差の絶対値、および算出したＡＰＤ９０と波形情報受付部１２が受け付けたＡＰＤ９０の差の絶対値の合計を算出する。また、Ａｄａｐｔｉｖｅ Ｐｉｅｃｅｗｉｓｅ Ｃｏｎｓｔａｎｔ Ａｐｐｒｏｘｉｍａｔｉｏｎ法などを用いて波形データを離散化しても良いし、比較している波形情報の差分の二乗和等を用い時間に対する電位の数字列の集合を直接比較しても良いし、ユークリッド距離等を計測して比較しても良い。

生体パラメータ情報生成手段１３２は、波形情報比較手段１３１で比較した活動電位波形情報が一致する、または最も近似する前記活動電位情報と対になっている１以上の生体パラメータ情報を含む生体パラメータセットを識別する。生体パラメータ情報生成手段１３２は、波形情報格納部１１から取得できるすべての波形情報を総当り的にチェックして、一致する、または最も近似する波形情報が有する生体パラメータセットを抽出しても良い。例えば、波形情報比較手段１３１が、算出したＡＰＤ３０と波形情報受付部１２が受け付けたＡＰＤ３０の差の絶対値、および算出したＡＰＤ６０と波形情報受付部１２が受け付けたＡＰＤ６０の差の絶対値、および算出したＡＰＤ９０と波形情報受付部１２が受け付けたＡＰＤ９０の差の絶対値の合計を算出する場合、生体パラメータ情報生成手段１３２は、合計値が最小の活動電位情報と対になる1以上の生体パラメータ情報を、波形情報格納部１１から取得する。また、生体パラメータ情報生成手段１３２は、ソートされている波形情報から二分探索等のアルゴリズムにより、高速に１以上の生体パラメータ情報を取得しても良い。また、生体パラメータ情報生成手段１３２は、ハイブリッド・ツリー・インデックス法、実験計画法、遺伝的アルゴリズム、再急降下法、応答局面法等の手法を用い、生体パラメータ情報を抽出しても良い。
以下、本生体パラメータ出力装置の動作について図３のフローチャートを用いて説明する。

（ステップＳ３０１）波形情報受付部１２は、活動電位波形情報の入力を受け付けたか否かを判断する。活動電位波形情報の入力を受け付ければステップＳ３０２に行き、活動電位波形情報の入力を受け付けなければステップＳ３０１に戻る。
（ステップＳ３０２）生体パラメータ情報生成手段１３２は、カウンタｉに１を代入する。

（ステップＳ３０３）生体パラメータ情報生成手段１３２は、新たな生体パラメータセット（１以上の生体パラメータ情報）を生成し、ｉ番目の再構築活動電位波形情報を波形情報格納部１１から取得する。なお、生体パラメータセットを生成する手法は様々である。なお、波形情報格納部１１は、図示しない生体シミュレート部におけるシミュレーションの結果、1以上の再構築活動電位波形情報を格納していても良い。
（ステップＳ３０４）波形情報比較手段１３１は、ｉ番目の再構築活動電位波形情報が存在するか否かを判断する。存在すればステップＳ３０５に行き、存在しなければステップＳ３０９に飛ぶ。

（ステップＳ３０５）波形情報比較手段１３１は、ステップＳ３０１で取得した活動電位情報とｉ番目の再構築活動電位波形情報の近似する度合いである近似度を算出する。なお、近似度の算出方法は、種々ある。近似度の算出方法の具体例は後述する。
（ステップＳ３０６）波形情報比較手段１３１は、ｉ番目の再構築活動電位波形情報に対応付けて、ステップＳ３０６で算出した近似度を一時蓄積する。
（ステップＳ３０７）波形情報比較手段１３１は、カウンタｉを１、インクリメントする。

（ステップＳ３０８）波形情報比較手段１３１は、あらかじめ指定された回数に達しておらず、近似度が所定の近似度より低いとき（条件に合致するとき）はステップＳ３０３に戻る。条件に合致しないときはステップＳ３０９に行く。なお、回数とは、あらゆる生体パラメータセットの組み合わせを網羅する回数である。
（ステップＳ３０９）生体パラメータ情報生成手段１３２は、近似度が最も大きい生体パラメータセット（１以上の生体パラメータ情報）を取得する。
（ステップＳ３１０）生体パラメータ情報出力部１４は、ステップＳ３０９で取得した１以上の生体パラメータ情報を出力する。ステップＳ３０１に戻る。
なお、図３のフローチャートにおいて、電源オフや処理終了の割り込みにより処理は終了する。

以下、本実施の形態における生体パラメータ出力装置の具体的な動作について説明する。今、波形情報格納部１１は、図４に示す再構築活動電位波形情報を保持している。本再構築活動電位波形情報は、３つの生体パラメータ識別子「ＩＫｒ」、「ＩＫ１」、および「ＩＫｓ」で識別される生体パラメータと、各生体パラメータの値と活動電位情報の関係を示した情報である。図４（ａ）は、生体パラメータ識別子「ＩＫｒ」の値を「０」から「５．０」まで採り、かつ「ＩＫ１」の値を「０．２」から「２．００」まで採った場合のＡＰＤ３０の値を示す図である。ここでは各生体パラメータは１が正常の状態を表す数値である。なお、本再構築活動電位波形情報は、通常、生体シミュレーションを用い、取得した活動電位波形情報であり、波形情報格納部１１に格納された活動電位情報である。生体シミュレーションの代わりに実際に生体材料から計測した情報を用いても良い。

横軸が「ＩＫｒ」の値、縦軸が「ＩＫ１」の値であり、交点の矩形の色がＡＰＤ３０の値を示す。本図４においては、モノクロであるので、概ね、矩形の色の濃さがＡＰＤ３０の値を示すとことなる。横軸、縦軸および矩形の色の濃さの意味は、図４中の他の図においても、同様に考える。
図４（ｂ）は、生体パラメータ識別子「ＩＫｓ」の値を「０」から「５．０」まで採り、かつ「ＩＫ１」の値を「０．２」から「２．００」まで採った場合のＡＰＤ３０の値を示す図である。
図４（ｃ）は、生体パラメータ識別子「ＩＫｓ」の値を「０」から「５．０」まで採り、かつ「ＩＫｒ」の値を「０」から「５．０」まで採った場合のＡＰＤ３０の値を示す図である。
図４（ｄ）は、生体パラメータ識別子「ＩＫｒ」の値を「０」から「５．０」まで採り、かつ「ＩＫ１」の値を「０．２」から「２．００」まで採った場合のＡＰＤ６０の値を示す図である。
図４（ｅ）は、生体パラメータ識別子「ＩＫｓ」の値を「０」から「５．０」まで採り、かつ「ＩＫ１」の値を「０．２」から「２．００」まで採った場合のＡＰＤ６０の値を示す図である。
図４（ｆ）は、生体パラメータ識別子「ＩＫｓ」の値を「０」から「５．０」まで採り、かつ「ＩＫｒ」の値を「０」から「５．０」まで採った場合のＡＰＤ６０の値を示す図である。
図４（ｇ）は、生体パラメータ識別子「ＩＫｒ」の値を「０」から「５．０」まで採り、かつ「ＩＫ１」の値を「０．２」から「２．００」まで採った場合のＡＰＤ９０の値を示す図である。
図４（ｈ）は、生体パラメータ識別子「ＩＫｓ」の値を「０」から「５．０」まで採り、かつ「ＩＫ１」の値を「０．２」から「２．００」まで採った場合のＡＰＤ９０の値を示す図である。
図４（ｉ）は、生体パラメータ識別子「ＩＫｓ」の値を「０」から「５．０」まで採り、かつ「ＩＫｒ」の値を「０」から「５．０」まで採った場合のＡＰＤ９０の値を示す図である。

なお、図４は、例えば、図５のような再構築活動電位波形情報管理表から取得され得る。なお、図５（ａ）は、「ＩＤ」、「ＩＫｓ」、「ＩＫ１」、「ＡＰＤ３０」、「ＡＰＤ６０」および「ＡＰＤ９０」を有するレコードを複数保持している表である。「ＩＤ」は、レコードを識別する情報であり、表におけるレコード管理のために存在する。「ＩＫｓ」は生体パラメータ識別子「ＩＫｓ」の値、「ＩＫ１」は生体パラメータ識別子「ＩＫ１」の値、「ＡＰＤ３０」、「ＡＰＤ６０」および「ＡＰＤ９０」は、それぞれ、「ＩＫｓ」の値および「ＩＫ１」の値が特定された場合の「ＡＰＤ３０」、「ＡＰＤ６０」および「ＡＰＤ９０」の値である。図５（ｂ）は、生体パラメータ識別子「ＩＫｓ」と「ＩＫ１」の値と「ＡＰＤ３０」、「ＡＰＤ６０」および「ＡＰＤ９０」の関係を示す表である。図５（ｃ）は、生体パラメータ識別子「ＩＫｓ」と「ＩＫｒ」の値と「ＡＰＤ３０」、「ＡＰＤ６０」および「ＡＰＤ９０」の関係を示す表である。
図４は、図５の情報を取得することによって、表示されるビューである、とも言える。
なお、波形からＡＰＤ３０、ＡＰＤ６０、ＡＯＤ９０の値を得ることは、公知技術により可能である。

かかる状況で、図６に示す活動電位波形情報が入力された、とする。すると、波形情報受付部１２は、図６に示す活動電位波形情報の入力を付ける。次に、波形情報比較手段１３１は、波形情報受付部１２が受け付けた活動電位波形情報のＡＰＤ３０、ＡＰＤ６０およびＡＰＤ９０の活動電位を示す情報である活動電位情報を取得する。ここで、波形情報比較手段１３１は、ＡＰＤ３０を７８、ＡＰＤ６０を１１９、ＡＰＤ９０を１２３と取得した、とする。

次に、波形情報比較手段１３１は、１番目から順にすべての再構築活動電位波形情報と、取得した「ＡＰＤ３０：７８」、「ＡＰＤ６０：１１９」、「ＡＰＤ９０：１２３」との近似度を求める。ここでは、再構築活動電位波形情報のＡＰＤの値と取得したＡＰＤの値の差の絶対値の総和の逆数を近似度とする。具体的な近似度の算出方法を、以下に示す。

まず、波形情報比較手段１３１は、１番目の再構築活動電位波形情報「ＩＫｒ＝０」「ＩＫ１＝０．２」「ＩＫｓ＝１」の場合の再構築活動電位波形情報を取得する。具体的には、図５（ａ）より、再構築活動電位波形情報は、「ＡＰＤ３０＝９３」「ＡＰＤ６０＝１６１」「ＡＰＤ９０＝１１１」である。

そして、次に、波形情報比較手段１３１は、再構築活動電位波形情報のＡＰＤ３０の値「７８」と、ＡＰＤ３０の値の差の絶対値、つまり、「｜７８−９３｜」を算出し、「１５」を得る。次に、波形情報比較手段１３１は、再構築活動電位波形情報のＡＰＤ６０の値「１１９」と、ＡＰＤ６０の値の差の絶対値を算出する。具体的には、波形情報比較手段１３１は、「｜１１９−１６１｜＝４２」を得る。次に、波形情報比較手段１３１は、再構築活動電位波形情報のＡＰＤ９０の値「１２３」と、ＡＰＤ９０の値の差の絶対値を算出する。具体的には、波形情報比較手段１３１は、「｜１２３−１１１｜＝１２」を得る。次に、波形情報比較手段１３１は、ＡＰＤ３０の差、ＡＰＤ６０の差、ＡＰＤ９０の差の総和を求める。つまり、波形情報比較手段１３１は、「１５＋４２＋１２＝６９」を得る。次に、波形情報比較手段１３１は、総和の逆数である近似度「１／６９」を得る。そして、波形情報比較手段１３１は、各生体パラメータの値、差の総和「６９」、近似度「１／６９」を有するレコードを一時記憶し、近似度管理表を構成する。かかる近似度管理表を図７に示す。図７において、近似度管理表は、「ＩＤ」「ＩＫｒ」「ＩＫ１」「ＩＫｓ」「差の総和」「近似度」を有するレコードを１以上、保持している。「ＩＤ」は、レコードを識別する情報であり、表におけるレコード管理のために存在する。「ＩＫｒ」「ＩＫ１」「ＩＫｓ」は、それぞれ各生体パラメータの値である。「差の総和」「近似度」は、上記で求めた値である。

同様に、波形情報比較手段１３１は、２番目の再構築活動電位波形情報「ＩＫｒ＝０．０」「ＩＫ１＝０．４」「ＩＫｓ＝１」の場合の再構築活動電位波形情報を取得する。具体的には、図５（ａ）より、再構築活動電位波形情報は、「ＡＰＤ３０＝１３３」「ＡＰＤ６０＝２４０」「ＡＰＤ９０＝２６９」である。そして、波形情報比較手段１３１は、例えば、差の総和「３２２」、近似度「１／３２２」を算出した、とする。
そして、波形情報比較手段１３１は、ここでは、総当り的にすべての「ＩＫｒ」「ＩＫ１」「ＩＫｓ」の値の組み合わせについて近似度を求める。

すべての「ＩＫｒ」「ＩＫ１」「ＩＫｓ」の値の組み合わせについて近似度を求めた後、生体パラメータ情報生成手段１３２は、最も近似度が大きいレコードを取得する。ここでは、最も近似度が大きいレコードを「ＩＤ＝ｎ」のレコードであるとする。そして、生体パラメータ情報生成手段１３２は、「ＩＤ＝ｎ」のレコードが有する「ＩＫｒ」「ＩＫ１」「ＩＫｓ」の値を取得する。
次に、生体パラメータ情報出力部１４は、「ＩＫｒ」「ＩＫ１」「ＩＫｓ」の値を出力する。かかる出力例を図８に示す。
次いで、薬物「ＡＢＣ」を投薬した後得られた活動電位波形情報を入力し、得られた生体パラメータの値を上記と比較することで、薬剤効果を数値で比較できる。

以上、本実施の形態によれば、薬物投入等の効果等を示す生体パラメータを出力する生体パラメータ出力装置を提供できる。かかる生体パラメータ出力装置により、薬剤の評価試験を減少させ、薬効を迅速に評価することが可能となる。したがって、創薬プロセスが短縮され得る。

なお、本実施の形態の具体例によれば、薬物「ＡＢＣ」を投入した場合の生体パラメータの値を出力したが、例えば、生体パラメータ出力装置は、病気を患う患者の活動電位波形情報を受け付け、かかる患者の生体パラメータの値を出力しても良い。かかる処理により、患者の病気が推定できたり、適切な薬剤を投与できたりする効果が生じる。

また、本実施の形態の具体例によれば、説明の簡単化のために３種類の生体パラメータを用いて、生体パラメータ出力装置の動作を説明したが、上述したように、対象となる生体パラメータは、数百にも及ぶことが通常である。

また、本実施の形態の具体例によれば、生体パラメータを決定するために、再構築活動電位波形情報を総当り的に評価したが、例えば、生体パラメータの値がソートされている場合などは、二分探索等のアルゴリズムを用いて、検索対象の再構築活動電位波形情報を早期に絞り込むなどの検索方法を用いることは好適である。また、ハイブリッド・ツリー・インデックス法、実験計画法、遺伝的アルゴリズム、再急降下法、応答局面法等の手法を用いることも好適である。

また、本実施の形態によれば、ＡＰＤ３０、ＡＰＤ６０およびＡＰＤ９０の活動電位情報に関して、入力された活動電位波形情報と再構築活動電位波形情報の近似度を取得したが、活動電位情報中の他の情報を用いて、入力された活動電位波形情報と再構築活動電位波形情報の近似度を取得しても良いことは言うまでもない。ただし、活動電位波形の特性、および処理の高速化の観点から、ＡＰＤ３０、ＡＰＤ６０およびＡＰＤ９０の３つの活動電位情報を用いることは好適である。

さらに、本実施の形態によれば、再構築活動電位波形情報は、１以上含む生体パラメータセットを入力として受け取り、シミュレーションを用いて取得した情報であったが、例えば、人手で入力した情報であっても良い。また、再構築活動電位波形情報は、生体パラメータ出力装置が有しないシミュレータもしくは実験計測装置により取得された情報であっても良い。また、上述した生体パラメータ出力装置は、活動電位波形の代わりに任意の生体の挙動を示す時系列値を実波形データとして入力し、生体シミュレーションとして入力された波形を再現しうる任意のシミュレーションを用いて、生体パラメータを出力しても良い。つまり、生体パラメータを取得する元になる情報は、活動電位波形情報に限らない。具体的には、上述した生体パラメータ出力装置において、前記波形情報格納部は、生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を１以上、および任意の生体の挙動を示す時系列値である実波形データを有し、前記波形情報受付部は、生体から取得した情報を受け付け、前記生体パラメータ情報取得部は、前記波形情報受付部が受け付けた情報に基づいて、１以上の生体パラメータ情報を取得する。例えば、任意の生体の挙動を示す時系列値を実波形データとは、心電図波形、血圧の時系列データ、心拍出量の時系列データ、血糖値や血中酸素濃度などの時系列データであり、それらの波形に影響を及ぼす生体パラメータとしては心筋収縮力などが挙げられる。

なお、かかる場合の生体パラメータ出力装置は、上記の生体パラメータ出力装置に対して、前記波形情報受付部は、生体の挙動を示す任意の時系列値データを受け付け、前記生体パラメータ情報取得部は、前記波形情報受付部が受け付けた情報に基づいて、１以上の生体パラメータ情報を取得する生体パラメータ出力装置である。また、方法およびプログラムは、上記の方法およびプログラムに対して、前記波形情報受付ステップにおいて、生体の挙動を示す任意の時系列値データを受け付け、前記生体パラメータ情報取得ステップにおいて、前記波形情報受付ステップで受け付けた情報に基づいて、１以上の生体パラメータ情報を取得する方法およびプログラムである。

上記のことから、本実施の形態において、以下の生体パラメータ出力方法を実現すれば良い。つまり、本方法は、活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付ステップと、前記活動電位波形情報に基づいて、１以上の生体パラメータ情報を取得する生体パラメータ情報取得ステップと、前記生体パラメータ情報取得ステップで取得した１以上の生体パラメータ情報を出力する生体パラメータ情報出力ステップを具備する生体パラメータ出力方法である。

さらに、本実施の形態における処理は、ソフトウェアで実現しても良い。そして、このソフトウェアをソフトウェアダウンロード等により配布しても良い。また、このソフトウェアをＣＤ−ＲＯＭなどの記録媒体に記録して流布しても良い。なお、本実施の形態における生体パラメータ出力装置を実現するソフトウェアは、以下のようなプログラムである。つまり、このプログラムは、コンピュータに、活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付ステップと、前記活動電位波形情報と実験情報に基づいて、１以上の生体パラメータ情報を取得する生体パラメータ情報取得ステップと、前記生体パラメータ情報取得ステップで取得した１以上の生体パラメータ情報を出力する生体パラメータ情報出力ステップを実行させるためのプログラム、である。なお、このプログラムを実行するコンピュータは、単数であってもよく、複数であってもよい。すなわち、集中処理を行ってもよく、あるいは分散処理を行ってもよい。

また、上記実施の形態において、各処理（各機能）は、単一の装置（システム）によって集中処理されることによって実現されてもよく、あるいは、複数の装置によって分散処理されることによって実現されてもよい。
本発明は、以上の実施の形態に限定されることなく、種々の変更が可能であり、それらも本発明の範囲内に包含されるものであることは言うまでもない。

以上のように、本発明にかかる生体パラメータ出力装置は、薬物投入等の効果等を示す生体パラメータを推定することができるという効果を有し、創薬のために利用するシミュレーション装置等として有用である。

生体パラメータ出力装置のブロック図 ＡＰＤ３０、ＡＰＤ６０、ＡＰＤ９０を説明する図 生体パラメータ出力装置の動作について説明するフローチャート 再構築活動電位波形情報の例を示す図 再構築活動電位波形情報管理表を示す図 入力された活動電位波形情報を示す図 近似度管理表を示す図 生体パラメータの出力例を示す図

Claims (9)

1. Ｎａチャネル、Ｃａチャネル、Ｋ_ＡＴＰチャネル、Ｋｒチャネル、Ｋ１チャネル、Ｋｓチャネルのうちの１以上のチャネルを含む細胞のチャネルを流れる電流、前記チャネルの開閉速度、イオン親和性、細胞内外のイオン濃度のうちの一の生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を２以上含む生体パラメータセットと、および生体の細胞から得られた活動電位波形の情報または膜電位の波形を構成する点の情報の集合である活動電位波形情報を有する波形情報を１以上格納している波形情報格納部と、
   活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付部と、
   前記波形情報格納部の波形情報の一部または全部と前記波形情報受付部が受け付けた活動電位波形情報に基づいて、１以上の生体パラメータ情報を含む生体パラメータセットを識別し、取得する生体パラメータ情報取得部と、
   前記生体パラメータ情報取得部が取得した生体パラメータセットを出力する生体パラメータ情報出力部を具備する生体パラメータ出力装置。
2. 前記波形情報格納部に格納されている波形情報のうちの生体パラメータセットは入力された情報であり、
   前記波形情報格納部に格納されている波形情報のうちの活動電位波形情報は、前記入力された生体パラメータセットに生体シミュレーションを用いて取得した情報である請求項１記載の生体パラメータ出力装置。
3. 記憶媒体に、
   Ｎａチャネル、Ｃａチャネル、Ｋ_ＡＴＰチャネル、Ｋｒチャネル、Ｋ１チャネル、Ｋｓチャネルのうちの１以上のチャネルを含む細胞のチャネルを流れる電流、前記チャネルの開閉速度、イオン親和性、細胞内外のイオン濃度のうちの一の生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を２以上含む生体パラメータセットと、および生体の細胞から得られた活動電位波形の情報または膜電位の波形を構成する点の情報の集合である活動電位波形情報を有する波形情報を１以上格納しており、
   活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付ステップと、
   前記記憶媒体に格納されている波形情報の一部または全部と前記波形情報受付ステップで受け付けた活動電位波形情報に基づいて、１以上の生体パラメータ情報を含む生体パラメータセットを識別し、取得する生体パラメータ情報取得ステップと、
   前記生体パラメータ情報取得ステップで取得した生体パラメータセットを出力する生体パラメータ情報出力ステップを具備する生体パラメータ出力方法。
4. 前記格納されている波形情報のうちの生体パラメータセットは入力された情報であり、
   前記格納されている波形情報のうちの活動電位波形情報は、入力された生体パラメータセットに生体シミュレーションを用いて取得した情報である請求項３記載の生体パラメータ出力方法。
5. 記憶媒体に、
   Ｎａチャネル、Ｃａチャネル、Ｋ_ＡＴＰチャネル、Ｋｒチャネル、Ｋ１チャネル、Ｋｓチャネルのうちの１以上のチャネルを含む細胞のチャネルを流れる電流、前記チャネルの開閉速度、イオン親和性、細胞内外のイオン濃度のうちの一の生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を２以上含む生体パラメータセットと、および生体の細胞から得られた活動電位波形の情報または膜電位の波形を構成する点の情報の集合である活動電位波形情報を有する波形情報を１以上格納しており、
   コンピュータに、
   活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付ステップと、
   前記記憶媒体に格納されている波形情報の一部または全部と前記波形情報受付ステップで受け付けた活動電位波形情報に基づいて、１以上の生体パラメータ情報を含む生体パラメータセットを識別し、取得する生体パラメータ情報取得ステップと、
   前記生体パラメータ情報取得ステップで取得した1以上の生体パラメータセットを出力する生体パラメータ情報出力ステップを実行させるためのプログラム。
6. 前記格納されている波形情報のうちの生体パラメータセットは入力された情報であり、
   前記格納されている波形情報のうちの活動電位波形情報は、入力された生体パラメータセットに生体シミュレーションを用いて取得した情報である請求項５記載のプログラム。
7. 前記波形情報受付部は、生体の挙動を示す任意の時系列値データを受け付け、
   前記生体パラメータ情報取得部は、
   前記波形情報受付部が受け付けた情報に基づいて、１以上の生体パラメータ情報を含む生体パラメータセットを取得する請求項１または請求項２いずれか記載の生体パラメータ出力装置。
8. 前記波形情報受付ステップにおいて、生体の挙動を示す任意の時系列値データを受け付け、
   前記生体パラメータ情報取得ステップにおいて、
   前記波形情報受付ステップで受け付けた情報に基づいて、１以上の生体パラメータ情報を含む生体パラメータセットを取得する請求項３または請求項４いずれか記載の生体パラメータ出力方法。
9. 前記波形情報受付ステップにおいて、生体の挙動を示す任意の時系列値データを受け付け、
   前記生体パラメータ情報取得ステップにおいて、
   前記波形情報受付ステップで受け付けた情報に基づいて、１以上の生体パラメータ情報を含む生体パラメータセットを取得する請求項５または請求項６いずれか記載のプログラム。