WO2006022226A1

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

Info

Publication number: WO2006022226A1
Application number: PCT/JP2005/015209
Authority: WO
Inventors: Akinori Noma; Nobuaki Sarai; Hiroyuki Kawano; Shohei Hido
Original assignee: Kyoto University
Priority date: 2004-08-26
Filing date: 2005-08-22
Publication date: 2006-03-02
Also published as: EP1783490A4; US20070255504A1; JPWO2006022226A1; EP1783490A1; US7516017B2; JP4710021B2

Abstract

【課題】培養細胞、もしくは動物実験等により、薬物投与時の活動電位波形変化を測定し、これを用いて、薬物の投入が生体パラメータに与える影響を評価することにより、細胞の個々のチャネルに対する作用を評価する手法を提供することを目的とする。【解決手段】本発明は、生体パラメータ識別子と生体パラメータ値の組である生体パラメータ情報を１以上、および活動電位情報を有する波形情報を１以上格納しており、活動電位波形情報の入力を受け付ける波形情報受付部と、前記活動電位波形情報に基づいて、１以上の生体パラメータ情報を取得する生体パラメータ情報取得部と、前記取得した１以上の生体パラメータ情報を出力する生体パラメータ情報出力部を具備する生体パラメータ出力装置であって、かかる生体パラメータ出力装置により、細胞の個々のチャネルに対する作用を評価することができる。

Description

明細書

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

技術分野

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

背景技術

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

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

[0004] また、動脈硬化の段階に応じた複数の動脈圧波形の基準情報を記憶しておき、生体の動脈カゝら得られた動脈圧波形情報と比較し、類似してヽる波形の年齢を出力する動脈硬化度測定方法がある (特許文献 2、特許文献 3参照)。

さらに、関連する技術として、生体パラメータセットを入力として受け取り、細胞をシミユレーシヨンし、活動電位波形情報を取得するシミュレーション装置がある (非特許文献 1参照)。

特許文献 1 :特開平 08— 016551 (第 1頁、第 1図等）

特許文献 2 :特開昭 61— 119252 (第 1頁、第 1図等）

特許文献 3 :特開平 03— 015439 (第 1頁、第 1図等）

非特許文献 1：皿井伸明、野間昭典「simBio：生物学的ダイナミックモデル開発基盤」日本ェム 'ィ一学会雑誌 BME, vol. 18, No. 2, p. 3— 11, 2004 (2004年 2月発行）

発明の開示

発明が解決しょうとする課題

[0005] し力しながら、上記の試験は、非常に膨大な評価試験であって、時間とコストが膨大にかかってしまう。

また、特許文献 1におけるシミュレーション方法は、薬物の生体組織内での拡散現象を解析する方法であり、薬物投入等の効果等を示す生体パラメータを推定する方法ではない。

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

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

課題を解決するための手段

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

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

[0010] さらに、第三の発明は、第一、第二の発明に対して、前記活動電位波形情報は、 A PD30、 APD60及び APD90の活動電位波形を特徴づける情報であり、生体パラメータ情報取得部は、前記活動電位波形情報の APD30、 APD60および APD90の活動電位波形情報を取得し、当該取得した 3つの活動電位波形情報と一致するまたは最も近似する 1以上の生体パラメータ情報を含む生体パラメータセットを取得する生体パラメータ出力装置である。

力かる構成により、高速かつ少ない CPUパワーで生体パラメータセットを取得することができる。

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

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

発明を実施するための最良の形態 [0013] 本実施の形態において、薬物投入等の効果を示したり、現在の生体の健康状態等を判断したりするために、生体パラメータを出力する装置等について述べる。なお、生体パラメータには細胞の種々のチャネル（たとえば、 Naチャネル、 Caチャネル、 K

A

チャネル、 Krチャネル、 K1チャネル、 Ksチャネルなど）を流れる電流や、各チヤネ

TP

ルの開閉速度や、イオン親和性、細胞内外のイオン濃度など、数百にも及ぶパラメ一タがある。

以下、生体パラメータ出力装置等の実施形態について図面を参照して説明する。なお、実施の形態において同じ符号を付した構成要素は同様の動作を行うので、再度の説明を省略する場合がある。

[0014] 図 1は、生体パラメータ出力装置のブロック図である。生体パラメータ出力装置は、波形情報格納部 11、波形情報受付部 12、生体パラメータ情報取得部 13、生体パラメータ情報出力部 14を具備する。生体パラメータ情報取得部 13は、波形情報比較手段 131、生体パラメータ情報生成手段 132を具備する。

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

[0016] なお、波形情報は、通常、図示しない手段により、以下のように取得されたものである。つまり、波形情報を取得する装置等（図示しない）は、生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のノメータの値である生体パラメータ値の組である生体パラメータ情報を 1以上含む生体パラメータセットを入力として受け取り、それに対応する活動電位波形情報を取得する。これには、生体を用いて活動電位波形情報を取得しても良い。また、波形情報を取得する装置力細胞をシミュレーションする細胞モデルに生体パラメータセットを適用し、シミュレーションを用い、活動電位波形情報（以下、適宜、シミュレーションを用いて取得した活動電位波形情報を「再構築活動電位波形情報」という。）を取得しても良い。この細胞モデルは、最も良く実験結果を再現しうるモデルに交換することが出来る。このとき、生体パラメータ出力装置は、波形情報を取得する装置を内蔵しており、シミュレーション結果をその都度生成しても良い。また、生体パラメータ出力装置は、あらかじめ生成された活動電位波形情報をデータベースより読み出しても良い。さらに、活動電位波形情報は、グラフの形式でも良いし、レコードの形式でも良い、つまり

、活動電位波形情報のデータ構造は問わない。波形情報格納部 11は、不揮発性の記録媒体が好適であるが、揮発性の記録媒体でも実現可能である。なお、生体パラメータセットを 1以上、およびそれに対応する活動電位波形情報を 1以上格納している記憶媒体を波形情報データベース（図示しない）とする。なお、生体パラメータセットを入力として受け取り、細胞をシミュレーションし活動電位波形情報を取得するシミユレーシヨン装置 (波形情報を取得する装置の一例）につ!/、ては公知技術である（上記、非特許文献 1参照)。

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

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

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

[0020] 波形情報比較手段 131は、波形情報格納部 11から取得した活動電位波形情報と波形情報受付部 12が受け付けた活動電位波形情報の近似度を判定する。波形情報比較手段 131は、例えば、以下のように APD30等の活動電位特徴情報を取得し、比較しても良い。つまり、波形情報受付部 12が受け付けた活動電位波形情報が波形を形成する電位 (mV)と時間 (ms)の数字列の集合である、とする。その場合、波形情報比較手段 131は、「最高の電位（最高の電位最低の電位） X O. 3」の値と対になる 2つの時間 (T、 T )を取得する。次に、波形情報比較手段 131は、「

1 2 I T

1

— T I」（Tと Tの差の絶対値）を APD30として求める。波形情報比較手段 131は

2 1 2

、 APD60、 APD90についても、同様の方法で求める。また、波形情報受付部 12が受け付けた活動電位波形情報は、 APD30、 APD60、および APD90の値である、とする。そして、波形情報比較手段 131は、算出した APD30と波形情報受付部 12 が受け付けた APD30の差の絶対値、および算出した APD60と波形情報受付部 12 が受け付けた APD60の差の絶対値、および算出した APD90と波形情報受付部 12 が受け付けた APD90の差の絶対値の合計を算出する。また、 Adaptive Piecewis e Constant Approximation法などを用いて波形データを離散化しても良いし、比較して!/、る波形情報の差分の二乗和等を用い時間に対する電位の数字列の集合を直接比較しても良いし、ユークリッド距離等を計測して比較しても良い。

[0021] 生体パラメータ情報生成手段 132は、波形情報比較手段 131で比較した活動電位波形情報が一致する、または最も近似する前記活動電位情報と対になっている 1以上の生体パラメータ情報を含む生体パラメータセットを識別する。生体パラメータ情報生成手段 132は、波形情報格納部 11から取得できるすべての波形情報を総当り的にチヱックして、一致する、または最も近似する波形情報が有する生体パラメータセットを抽出しても良い。例えば、波形情報比較手段 131が、算出した APD30と波形情報受付部 12が受け付けた APD30の差の絶対値、および算出した APD60と波形情報受付部 12が受け付けた APD60の差の絶対値、および算出した APD90と波形情報受付部 12が受け付けた APD90の差の絶対値の合計を算出する場合、生体ノメータ情報生成手段 132は、合計値が最小の活動電位情報と対になる 1以上の生体パラメータ情報を、波形情報格納部 11から取得する。また、生体パラメータ情報生成手段 132は、ソートされている波形情報力二分探索等のアルゴリズムにより、高速に 1以上の生体パラメータ情報を取得しても良い。また、生体パラメータ情報生成手段 132は、ハイブリッド'ツリー'インデックス法、実験計画法、遺伝的ァルゴリズム、再急降下法、応答局面法等の手法を用い、生体パラメータ情報を抽出しても良い。

以下、本生体パラメータ出力装置の動作について図 3のフローチャートを用いて説明する。

[0022] (ステップ S301)波形情報受付部 12は、活動電位波形情報の入力を受け付けたか否かを判断する。活動電位波形情報の入力を受け付ければステップ S302に行き、活動電位波形情報の入力を受け付けなければステップ S301に戻る。

(ステップ S302)生体パラメータ情報生成手段 132は、カウンタ iに 1を代入する。

[0023] (ステップ S303)生体パラメータ情報生成手段 132は、新たな生体パラメータセット

(1以上の生体パラメータ情報)を生成し、潘目の再構築活動電位波形情報を波形情報格納部 11から取得する。なお、生体パラメータセットを生成する手法は様々である。なお、波形情報格納部 11は、図示しない生体シミュレート部におけるシミュレーシヨンの結果、 1以上の再構築活動電位波形情報を格納して!、ても良、。

(ステップ S304)波形情報比較手段 131は、潘目の再構築活動電位波形情報が存在するか否かを判断する。存在すればステップ S305に行き、存在しなければステップ S 309に飛ぶ。

[0024] (ステップ S305)波形情報比較手段 131は、ステップ S 301で取得した活動電位情報と i番目の再構築活動電位波形情報の近似する度合いである近似度を算出する。なお、近似度の算出方法は、種々ある。近似度の算出方法の具体例は後述する。

(ステップ S306)波形情報比較手段 131は、 i番目の再構築活動電位波形情報に対応付けて、ステップ S306で算出した近似度を一時蓄積する。

(ステップ S307)波形情報比較手段 131は、カウンタ iを 1、インクリメントする。

[0025] (ステップ S308)波形情報比較手段 131は、あら力じめ指定された回数に達しておらず、近似度が所定の近似度より低いとき（条件に合致するとき）はステップ S303に戻る。条件に合致しないときはステップ S309に行く。なお、回数とは、あらゆる生体パラメータセットの組み合わせを網羅する回数である。

(ステップ S309)生体パラメータ情報生成手段 132は、近似度が最も大きい生体パラメータセット（1以上の生体パラメータ情報）を取得する。

(ステップ S310)生体パラメータ情報出力部 14は、ステップ S309で取得した 1以上の生体パラメータ情報を出力する。ステップ S301に戻る。

なお、図 3のフローチャートにおいて、電源オフや処理終了の割り込みにより処理は終了する。

[0026] 以下、本実施の形態における生体パラメータ出力装置の具体的な動作について説明する。今、波形情報格納部 11は、図 4に示す再構築活動電位波形情報を保持している。本再構築活動電位波形情報は、 3つの生体パラメータ識別子「IKr」、「IK1」、および「IKs」で識別される生体パラメータと、各生体パラメータの値と活動電位情報の関係を示した情報である。図 4 (a)は、生体パラメータ識別子「IKr」の値を「0」から「5. 0」まで採り、かつ「IK1」の値を「0. 2」から「2. 00」まで採った場合の APD30 の値を示す図である。ここでは各生体パラメータは 1が正常の状態を表す数値である。なお、本再構築活動電位波形情報は、通常、生体シミュレーションを用い、取得した活動電位波形情報であり、波形情報格納部 11に格納された活動電位情報である

。生体シミュレーションの代わりに実際に生体材料力計測した情報を用いても良い

[0027] 横軸が「IKr」の値、縦軸が「IK1」の値であり、交点の矩形の色が APD30の値を示す。本図 4においては、モノクロであるので、概ね、矩形の色の濃さが APD30の値を示すとことなる。横軸、縦軸および矩形の色の濃さの意味は、図 4中の他の図においても、同様に考える。

図 4 (b)は、生体パラメータ識別子「IKs」の値を「0」から「5. 0」まで採り、かつ「IK1 」の値を「0. 2」から「2. 00」まで採った場合の APD30の値を示す図である。

図 4 (c)は、生体パラメータ識別子「IKs」の値を「0」から「5. 0」まで採り、かつ「IKr」の値を「0」から「5. 0」まで採った場合の APD30の値を示す図である。

図 4 (d)は、生体パラメータ識別子「IKr」の値を「0」から「5. 0」まで採り、かつ「IK1 」の値を「0. 2」から「2. 00」まで採った場合の APD60の値を示す図である。

図 4 (e)は、生体パラメータ識別子「IKs」の値を「0」から「5. 0」まで採り、かつ「IK1 」の値を「0. 2」から「2. 00」まで採った場合の APD60の値を示す図である。

図 4 (f)は、生体パラメータ識別子「IKs」の値を「0」から「5. 0」まで採り、かつ「IKr」の値を「0」から「5. 0」まで採った場合の APD60の値を示す図である。

図 4 (g)は、生体パラメータ識別子「IKr」の値を「0」から「5. 0」まで採り、かつ「IK1 」の値を「0. 2」から「2. 00」まで採った場合の APD90の値を示す図である。

図 4 (h)は、生体パラメータ識別子「IKs」の値を「0」から「5. 0」まで採り、かつ「IK1 」の値を「0. 2」から「2. 00」まで採った場合の APD90の値を示す図である。

図 4 (i)は、生体パラメータ識別子「IKs」の値を「0」から「5. 0」まで採り、かつ「IKr」の値を「0」から「5. 0」まで採った場合の APD90の値を示す図である。

[0028] なお、図 4は、例えば、図 5のような再構築活動電位波形情報管理表力取得され得る。なお、図 5 (a)は、「ID」、「IKs」、「IK1」、「APD30」、「APD60」および「APD 90」を有するレコードを複数保持している表である。「ID」は、レコードを識別する情報であり、表におけるレコード管理のために存在する。「IKs」は生体パラメータ識別子「IKs」の値、「IK1」は生体パラメータ識別子「IK1」の値、「APD30」、「APD60」および「APD90」は、それぞれ、「IKs」の値および「IK1」の値が特定された場合の「 APD30」、「APD60」および「APD90」の値である。図 5 (b)は、生体パラメータ識別子「IKs」と「IK1」の値と「APD30」、「APD60」および「APD90」の関係を示す表である。図 5 (c)は、生体パラメータ識別子「IKs」と「IKr」の値と「APD30」、「APD60」および「APD90」の関係を示す表である。

図 4は、図 5の情報を取得することによって、表示されるビューである、とも言える。なお、波形から APD30、 APD60、 AOD90の値を得ることは、公知技術により可能である。

[0029] かかる状況で、図 6に示す活動電位波形情報が入力された、とする。すると、波形情報受付部 12は、図 6に示す活動電位波形情報の入力を付ける。次に、波形情報比較手段 131は、波形情報受付部 12が受け付けた活動電位波形情報の APD30、 APD60および APD90の活動電位を示す情報である活動電位情報を取得する。ここで、波形情報比較手段 131は、 APD30を 78、 APD60を 119、 APD90を 123と取得した、とする。

[0030] 次に、波形情報比較手段 131は、 1番目から順にすベての再構築活動電位波形情報と、取得した「APD30 : 78」、「APD60 : 119」、「APD90 : 123」との近似度を求める。ここでは、再構築活動電位波形情報の APDの値と取得した APDの値の差の絶対値の総和の逆数を近似度とする。具体的な近似度の算出方法を、以下に示す。

[0031] まず、波形情報比較手段 131は、 1番目の再構築活動電位波形情報「IKr=0」「I K1 = 0. 2」「IKs = l」の場合の再構築活動電位波形情報を取得する。具体的には、図 5 (a)より、再構築活動電位波形情報は、「APD30 = 93」「APD60= 161」「APD 90 = 111」である。

[0032] そして、次に、波形情報比較手段 131は、再構築活動電位波形情報の APD30の値「78」と、 APD30の値の差の絶対値、つまり、「 I 78— 93 |」を算出し、「15」を得る。次に、波形情報比較手段 131は、再構築活動電位波形情報の APD60の値「11 9」と、 APD60の値の差の絶対値を算出する。具体的には、波形情報比較手段 131 は、「 I 119— 161 I =42」を得る。次に、波形情報比較手段 131は、再構築活動電位波形情報の APD90の値「123」と、 APD90の値の差の絶対値を算出する。具体的には、波形情報比較手段 131は、「 I 123— 111 I = 12」を得る。次に、波形情報比較手段 131は、 APD30の差、 APD60の差、 APD90の差の総和を求める。つまり、波形情報比較手段 131は、「15+42+ 12 = 69」を得る。次に、波形情報比較手段 131は、総和の逆数である近似度「1Z69」を得る。そして、波形情報比較手段 131は、各生体パラメータの値、差の総和「69」、近似度「1Z69」を有するレコードを一時記憶し、近似度管理表を構成する。かかる近似度管理表を図 7に示す。図 7 において、近似度管理表は、「ID」「IKr」「IK1」 riKsj「差の総和」「近似度」を有するレコードを 1以上、保持している。「ID」は、レコードを識別する情報であり、表におけるレコード管理のために存在する。 riKrJ Γΐκυ「IKs」は、それぞれ各生体パラメ一タの値である。「差の総和」「近似度」は、上記で求めた値である。

[0033] 同様に、波形情報比較手段 131は、 2番目の再構築活動電位波形情報「IKr=0.

0」 ΓΙΚ1 = 0. 4」「IKs = l」の場合の再構築活動電位波形情報を取得する。具体的には、図 5 (a)より、再構築活動電位波形情報は、「APD30= 133」「APD60 = 240 」「APD90 = 269」である。そして、波形情報比較手段 131は、例えば、差の総和「3 22」、近似度「1Z322」を算出した、とする。

そして、波形情報比較手段 131は、ここでは、総当り的にすべての「IKr」「IK1」「I KsJの値の組み合わせにつ!/、て近似度を求める。

[0034] すべての「IKr」「IK1」「IKs」の値の組み合わせにつ!/、て近似度を求めた後、生体ノメータ情報生成手段 132は、最も近似度が大きいレコードを取得する。ここでは、最も近似度が大きいレコードを「ID=n」のレコードであるとする。そして、生体パラメータ情報生成手段 132は、「ID=n」のレコードが有する「IKr」「IK1」「IKs」の値を取得する。

次に、生体パラメータ情報出力部 14は、「IKr」「IK1」「IKs」の値を出力する。かかる出力例を図 8に示す。

次いで、薬物「ABC」を投薬した後得られた活動電位波形情報を入力し、得られた生体パラメータの値を上記と比較することで、薬剤効果を数値で比較できる。

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

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

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

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

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

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

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

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

[0044] また、上記実施の形態にお!、て、各処理 (各機能）は、単一の装置 (システム）によつて集中処理されることによって実現されてもよぐあるいは、複数の装置によって分散処理されることによって実現されてもょ、。

本発明は、以上の実施の形態に限定されることなぐ種々の変更が可能であり、それらも本発明の範囲内に包含されるものであることは言うまでもない。

産業上の利用可能性

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

図面の簡単な説明

[0046] [図 1]生体パラメータ出力装置のブロック図

[図 2]APD30、 APD60、 APD90を説明する図

[図 3]生体パラメータ出力装置の動作について説明するフローチャート

圆 4]再構築活動電位波形情報の例を示す図

[図 5]再構築活動電位波形情報管理表を示す図

[図 6]入力された活動電位波形情報を示す図 [図 7]近似度管理表を示す図圆 8]生体パラメータの出力例を示す図

Claims

請求の範囲

[1] 生体のパラメータを識別する生体パラメータ識別子と当該生体パラメータ識別子で識別される生体のパラメータの値である生体パラメータ値の組である生体パラメータ情報を 1以上含む生体パラメータセットと、および活動電位波形を示す情報である活動電位波形情報を有する波形情報を 1以上格納している波形情報格納部と、活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付部と、

前記波形情報格納部の波形情報の一部または全部と前記波形情報受付部が受け付けた活動電位波形情報に基づ、て、 1以上の生体パラメータ情報を含む生体パラメータセットを識別し、取得する生体パラメータ情報取得部と、

前記生体パラメータ情報取得部が取得した生体パラメータセットを出力する生体パラメータ情報出力部を具備する生体パラメータ出力装置。

[2] 前記波形情報格納部に格納されている波形情報のうちの生体パラメータセットは入力された情報であり、

前記波形情報格納部に格納されている波形情報のうちの活動電位波形情報は、前記入力された生体パラメータセットに生体シミュレーションを用いて取得した情報である請求項 1記載の生体パラメータ出力装置。

[3] 活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付ステップと、

格納している波形情報の一部または全部と前記波形情報受付ステップで受け付けた活動電位波形情報に基づいて、 1以上の生体パラメータ情報を含む生体パラメータセットを識別し、取得する生体パラメータ情報取得ステップと、

前記生体パラメータ情報取得ステップで取得した生体パラメータセットを出力する生体パラメータ情報出力ステップを具備する生体パラメータ出力方法。

[4] 前記格納されて、る波形情報のうちの生体パラメータセットは入力された情報であり前記格納されてヽる波形情報のうちの活動電位波形情報は、入力された生体パラメータセットに生体シミュレーションを用いて取得した情報である請求項 3記載の生体パラメータ出力方法。

[5] コンピュータに、

活動電位波形に関する情報である活動電位波形情報の入力を受け付ける波形情報受付ステップと、

前記生体パラメータ情報取得ステップで取得した 1以上の生体パラメータセットを出力する生体パラメータ情報出力ステップを実行させるためのプログラム。

[6] 前記格納されて、る波形情報のうちの生体パラメータセットは入力された情報であり前記格納されてヽる波形情報のうちの活動電位波形情報は、入力された生体パラメータセットに生体シミュレーションを用いて取得した情報である請求項 5記載のプログラム。

[7] 前記波形情報受付部は、生体の挙動を示す任意の時系列値データを受け付け、前記生体パラメータ情報取得部は、

前記波形情報受付部が受け付けた情報に基づ!、て、 1以上の生体パラメータ情報を含む生体パラメータセットを取得する請求項 1または請求項 2いずれか記載の生体パラメータ出力装置。

[8] 前記波形情報受付ステップにおいて、生体の挙動を示す任意の時系列値データを受け付け、

前記生体パラメータ情報取得ステップにおいて、

前記波形情報受付ステップで受け付けた情報に基づ!、て、 1以上の生体パラメータ情報を含む生体パラメータセットを取得する請求項 3または請求項 4いずれか記載の生体パラメータ出力方法。

[9] 前記波形情報受付ステップにおいて、生体の挙動を示す任意の時系列値データを受け付け、

前記生体パラメータ情報取得ステップにおいて、前記波形情報受付ステップで受け付けた情報に基づ!、て、 1以上の生体パラメータ情報を含む生体パラメータセットを取得する請求項 5または請求項 6いずれか記載のプログラム。