(57)【要約】【課題】  昆虫内の化学物質の薬物動態的振舞を計算するためのコンピュータシステムを提供する。【解決手段】コンピュータシステムは、昆虫の区画内の化学物質の濃度／時間プロフィールを予測するための生理学に基づいた薬物動態シミュレーションモデル（１０２）であって、物質に依存するパラメータを少なくとも一つ有するものと、物質の物理化学特性に基づいて少なくとも一つのパラメータを予測するための予測モジュール（１１０）と、を備える。【選択図】図１
2004238390.tif 000002
【特許請求の範囲】
【請求項１】
昆虫内の化学物質の薬物動態的振舞を計算するためのコンピュータシステムにおいて、
昆虫の区画内の化学物質の濃度／時間プロフィールを予測するための生理学に基づいた薬物動態シミュレーションモデル（１０２）であって、物質に依存するパラメータを少なくとも一つ有するものと、
物質の物理化学特性に基づいて上記少なくとも一つのパラメータを予測するための予測モジュール（１１０）と、
を備えたコンピュータシステム。
【請求項２】
上記少なくとも一つのパラメータは、物質に対する透過率と区画の有効表面積の積であることを特徴とする請求項１のコンピュータシステム。
【請求項３】
上記少なくとも一つのパラメータは、区画間の質量輸送に対するレート係数（λｘ＝ＰＡｘ／Ｖｘ）であり、昆虫の少なくとも一つの器官の容積は時間の関数（Ｖｘ＝Ｖｘ（ｔ））であることを特徴とする請求項１のコンピュータシステム。
【請求項４】
上記少なくとも一つのパラメータは、昆虫の器官と血リンパとの間の平衡係数
2004238390.tif 000022
、クチクラの表面とクチクラとの間の平衡係数（Ｋｃ／ｃｓ）、あるいは、内臓壁と内臓中身との間の平衡係数（Ｋｇｃ／ｇｗ）であることを特徴とする請求項１のコンピュータシステム。
【請求項５】
上記物理化学特性は、物質の水と燐脂質膜の間の分配係数、オクタノール／水分配係数、分子量、溶解度、および／またはこれらパラメータの組み合わせであることを特徴とする請求項１のコンピュータシステム。
【請求項６】
物質の化学構造の記述子から物理化学特性を決定するＱＳＡＲモデルまたはニューラルネットワークを備えた請求項１のコンピュータシステム。
【請求項７】
予測モジュールは、テスト物質の物理化学特性およびこれらテスト物質に対し実験的に得られたパラメータを含むデータベース（１０８）に基づくことを特徴とする請求項１のコンピュータシステム。
【請求項８】
予測モジュールは、物質の脂肪親和性および／または分子量から上記少なくとも一つのパラメータの計算を行う計算関数を含むことを特徴とする請求項１のコンピュータシステム。
【請求項９】
計算関数は、実験的に得られたパラメータ値の直線回帰に基づくことを特徴とする請求項８のコンピュータシステム。
【請求項１０】
昆虫の区画内の化学物質の濃度／時間プロフィールの予測を行うための生理学に基づく昆虫の薬物動態シミュレーションモデルであって、物質に依存するパラメータを少なくとも一つ有するものを用いた、昆虫内の化学物質の薬物動態的振舞を計算するための方法において、
物質の物理化学特性を物質に対する上記少なくとも一つのパラメータを予測するための予想モジュールに入力する工程と、
上記予測された少なくとも一つのパラメータに基づいて、化学物質の濃度／時間プロフィールの予測を行うシミュレーションモデルを用いてシミュレーションを実行する工程と、
を含む方法。
【請求項１１】
上記少なくとも一つのパラメータは、物質に対する昆虫の区画の透過率と区画の有効表面積の積であることを特徴とする請求項１０の方法。
【請求項１２】
上記パラメータは、区画間の質量輸送に対するレート係数（λｘ＝ＰｘＡｘ／Ｖｘ）であり、昆虫の少なくとも一つの器官の容積は時間の関数（Ｖｘ＝Ｖｘ（ｔ））であることを特徴とする請求項１０の方法。
【請求項１３】
上記パラメータは、器官と血リンパとの間の平衡係数
2004238390.tif 000023
、クチクラの表面とクチクラとの間の平衡係数（Ｋｃ／ｃｓ）、あるいは、内臓壁と内臓中身との間の平衡係数（Ｋｇｃ／ｇｗ）であることを特徴とする請求項１０の方法。
【請求項１４】
上記物理化学特性は、物質の水と燐脂質膜の間の分配係数、オクタノール／水分配係数、分子量、溶解度、および／またはこれらパラメータの組み合わせであることを特徴とする請求項１０の方法。
【請求項１５】
上記物理化学特性は、ＱＳＡＲモデルまたはニューラルネットワークにより決定されることを特徴とする請求項１０の方法。
【請求項１６】
上記少なくとも一つのパラメータの予測は、テスト物質の物理化学特性およびこれらテスト物質に対し実験的に得られたパラメータを含むデータベース（１０８）に基づくことを特徴とする請求項１０の方法。
【請求項１７】
上記少なくとも一つのパラメータの予測は、物質の脂肪親和性および／または分子量から計算関数を用いて行われることを特徴とする請求項１０の方法。
【請求項１８】
計算関数は、実験的に得られたパラメータ値の直線回帰に基づくことを特徴とする請求項１０の方法。
【請求項１９】
昆虫の区画内の化学物質の濃度／時間プロフィールを予測するための生理学に基づいた薬物動態シミュレーションモデルであって物質に依存するパラメータを少なくとも一つ有するものを用いて、化学物質の薬物動態的振舞の計算を行うためのプログラムを記憶したデジタル記憶媒体において、
上記計算は、
物質に対する上記少なくとも一つのパラメータを予測するための予想モジュールに物質の物理化学特性を入力する工程と、
上記予測された少なくとも一つのパラメータに基づいて、化学物質の濃度／時間プロフィールの予測を行うシミュレーションモデルを用いてシミュレーションを実行する工程と、
を含むことを特徴とするデジタル記憶媒体。
【請求項２０】
上記少なくとも一つのパラメータは、物質に対する区画の透過率と区画の有効表面積の積であることを特徴とする請求項１０のデジタル記憶媒体。
【請求項２１】
上記少なくとも一つのパラメータは、区画間の質量輸送に対するレート係数（λｘ＝ＰｘＡｘ／Ｖｘ）であり、昆虫の少なくとも一つの器官の容積（Ｖｘ）は、時間の関数（Ｖｘ＝Ｖｘ（ｔ））であることを特徴とする請求項１９のデジタル記憶媒体。
【請求項２２】
上記少なくとも一つのパラメータは、器官と血リンパとの間の平衡係数
2004238390.tif 000024
、クチクラの表面とクチクラとの間の平衡係数（Ｋｃ／ｃｓ）、あるいは、内臓壁と内臓中身との間の平衡係数（Ｋｇｃ／ｇｗ）であることを特徴とする請求項１９のデジタル記憶媒体。
【請求項２３】
上記物理化学特性は、物質の水と燐脂質膜の間の分配係数、オクタノール／水分配係数、分子量、溶解度、および／またはこれらパラメータの組み合わせであることを特徴とする請求項１９のデジタル記憶媒体。
【請求項２４】
上記プログラムとは別にあるいは該プログラムの一部として、物質の化学構造の記述子から上記物理化学特性を決定するＱＳＡＲモデルまたはニューラルネットワークをさらに備えた請求項１９のデジタル記憶媒体。
【請求項２５】
上記予測モジュールは、テスト物質の物理化学特性およびこれらテスト物質に対し実験的に得られたパラメータを含むデータベース（１０８）に基づくことを特徴とする請求項１９のデジタル記憶媒体。
【請求項２６】
上記予測モジュールは、物質の脂肪親和性および／または分子量からパラメータの計算を行う計算関数を含むことを特徴とする請求項２９のデジタル記憶媒体。
【請求項２７】
上記計算関数は、実験的に得られたパラメータ値の直線回帰に基づくことを特徴とする請求項２６のデジタル記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生理学的モデルに基づいて昆虫内の化学物質の薬物動態的振舞を計算するためのコンピュータシステム、対応する方法、およびコンピュータプログラムプロダクトに関する。
【背景技術】
【０００２】
従来技術として哺乳類の生理学的に基づいたモデルが多数知られている（非特許文献１参照。）。さらに、従来技術としてイモムシの生理学に基づいたモデルがある（非特許文献２参照。）。この場合、昆虫は、の個々の器官をそれぞれ表す区画を用いて記述される。個々の区画の「相互作用」は、イモムシに関する公知の生理機能から得られる。このモデルの重要なパラメータは、分配速度を決定する区画間の質量輸送に対するレート係数、および、器官と哺乳類の血流に対応する血リンパとの間の、熱力学的平衡における濃度比を与える器官分配係数である。
【０００３】
こうした区画によるモデルを用いて、実験により決定されたある物質の薬物動態プロフィールを、レート係数および分配係数をマッチングさせることで過去にさかのぼって記述することは公知である（例えば、非特許文献３および４参照。）。
【非特許文献１】チャーニック（Charnick）その他、ジェー・ファーマコキン・ビオファーム（J.Pharmacokin. Biopharm.） 23, 217 (1995)
【非特許文献２】グリーンウッド（Greenwood）その他、ペスティック・サイ（Pestic. Sci.） 30, 97 (1990)
【非特許文献３】ラガディック（Lagadic）その他、ペスティック・バイオケム・フィシオル（Pestic. Biochem. Physiol.） 45, 105 (1993)
【非特許文献４】ラガディック（Lagadic）その他、ペスティック・バイオケム・フィシオル（Pestic. Biochem. Physiol.） 48, 173 (1994)
【発明の開示】
【発明が解決しようとする課題】
【０００４】
そこで、本発明は、昆虫内の化学物質の薬物動態的振舞を予測するコンピュータシステム、対応する方法、およびコンピュータプログラムプロダクトを提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明は、昆虫内の化学物質の薬物動態的振舞の予測値を特に効率的に計算することを可能にする。本発明は特に、物理化学パラメータに基づいて昆虫内の化学物質の吸収、分配、および排泄の評価を可能にする。
【０００６】
この目的のため、昆虫の生理学に基づく薬物動態シミュレーションモデルを用いて、昆虫の区画内の化学物質の濃度／時間プロフィールを予測する。シミュレーションモデルは、調査すべき物質に依存するパラメータを少なくとも一つ含む。シミュレーションモデルの一つまたはそれ以上のパラメータは、特定の物質に対して、該物質の一つまたはそれ以上の物理化学特性に基づいて予測される。
【０００７】
こうした物理化学パラメータは、例えば、水と燐脂質膜との間の分配係数またはオクタノール／水分配係数により記述される物質の脂肪親和性、分子量、あるいは溶解度である。物質の関連する物理化学パラメータは、直接実験により決定したり、あるいは、ＱＳＡＲ（Quantitative Structure Activity Relations)やニューラルネットワークなど公知の方法を用いて物質の化学構造の記述子から直接決定できる。
【０００８】
後者の場合、本発明の方法を用いて、仮想物質すなわちまだ合成されていない物質でさえ、昆虫におけるその吸収特性および分配特性を考慮することで評価できる。さらに、物理化学特性と薬物動態特性との間に確立した関係に基づいて、殺虫剤物質の最適化のための一般基準を導くことができる。
【０００９】
本発明に係る好適な実施形態では、区画間の質量輸送のレート係数は、物質に対する区画の透過性と区画の有効表面積の積に比例し、シミュレーションモデルの物質依存パラメータとして用いられる。したがって、シミュレーションモデルの各区画に対し、物質に対する関連区画の透過性および関連区画の有効表面積を含む物質依存パラメータがあるのが好適である。透過係数は、細胞膜を通る物質のフラックスを記述する量である。
【００１０】
特に有利な点は、物質の物理化学特性に基づいてシミュレーションモデル用のパラメータを決定するために、さらなる実験的研究が不要であることであることよりもむしろ、さらに物質の物理化学特性に基づいて一つまたはそれ以上のパラメータを決定できることである。この決定は、予め種々のテスト物質に対し実験で決定されたデータベースに基づいて行われる。データベースには、テスト物質に対し実験で決定されたシミュレーションモデルの物質依存パラメータと、テスト物質の物理化学特性とが含まれる。このデータベースは、調査すべき物質に対して、物質依存の一つまたはそれ以上のパラメータを予測するのに用いられる。
【００１１】
本発明に係る好適な実施形態では、データベースから直線回帰により計算関数が得られる。例えば、計算関数は、脂肪親和性や分子量の関数である。したがって、調査すべき物質に対するシミュレーションモデル用パラメータを予測するために、調査すべき物質の脂肪親和性および分子量を計算関数を用いて評価する必要があるだけで、これによりパラメータを得ることができる。その結果、予測されたパラメータを用いて、昆虫内の物質の吸収および排泄中の濃度／時間プロフィールの具体的なシミュレーションを実行することができる。
【００１２】
代わりに、直線回帰により得られる計算関数を用いる代わりに、公知の予測方法を用いてもよい。
【００１３】
特に有利なことは、所定数のテスト物質に関するデータベースが一旦得られた後は、他の物質に関して薬物動態的振舞のシミュレーションを行うために、追加の実験が不要であることである。これは、潜在的な殺虫剤に用いられる「候補」を、その薬物動態的振舞の観点から高処理量で評価できることを可能にする。したがって、新しい殺虫剤の調査、開発、および最適化を大きな加速度で行うことが可能である。
【発明を実施するための最良の形態】
【００１４】
以下、本発明に係る好適な実施形態を図面を参照して詳細に説明する。
【００１５】
図１はコンピュータシステム１００を示す。これは公知のパーソナルコンピュータ（ＰＣ）、ワークステーション、あるいはクライアント／サーバシステムが可能である。コンピュータシステム１００はシミュレーションモデル１０２を備える。これは、生理学に基づく昆虫の薬物動態モデルである。この目的のために、シミュレーションモデル１０２は、昆虫の区画を表し、したがって昆虫の区画内の物質に関する濃度／時間プロフィールの予測を可能にする。
【００１６】
シミュレーションモデル１０２は、記述すべき昆虫の種類にのみ依存する生理学的パラメータ１０４と、一つまたはそれ以上の物質依存パラメータ１０６とを含む。したがって、濃度／時間プロフィールのシミュレーション計算を実行するために、評価すべき物質に対する物質依存パラメータ１０６のパラメータ値を入力することが必要である。シミュレーションモデル１０２の一実施形態の例は、図３および４を参照して以下で詳述する。
【００１７】
コンピュータシステム１００はさらにデータバンク１０８を備える。データバンク１０８は、昆虫内のテスト物質の薬物動態的振舞の実験的研究に基づいて得られたデータベースを記憶するために用いられる。予め実験により調査された各テスト物質に対し、データバンク１０８には、関連するテスト物質に対し実験により決定された一つまたはそれ以上のパラメータ値と、関連するテスト物質の少なくとも一つの物理化学特性が含まれる。
【００１８】
データバンク１０８に記憶された上記データベースは、予測モジュール１１０において、新たに調査すべき物質に対するパラメータ値の予測のための基礎をなす。例えば、直線回帰法により、データバンク１０８に記憶されたデータベースから計算規則が得られる。この規則は、調査すべき物質の物理化学特性、例えば脂肪親和性や分子量から、シミュレーションモデル１０２を用いてシミュレーションを実行するのに必要な物質依存パラメータ値を獲得することを可能とする。
【００１９】
コンピュータシステム１００はさらに、調査すべき物質の物理化学特性を入力するための入力／出力モジュール１１２を備える。シミュレーションされた濃度／時間プロフィールはさらに、入力／出力モジュール１１２を介して出力される。
【００２０】
入力／出力モジュール１１２は、現実に存在する物質または例えば潜在的な殺虫剤など仮想の物質の記述子とそれらの物理化学特性とを含むデータバンクに接続してもよい。合成されていない仮想の物質の場合、コンピュータシステム１００への入力に必要な物理化学特性を、例えばＱＳＡＲやニューラルネットワークなど公知の方法を用いてテスト物質の化学構造の記述子から直接決定してもよい。仮想物質の場合、物理化学特性の代わりに、記述子を入力／出力モジュール１１２を介して入力してもよい。この場合、物理化学特性は、記述子に基づいて、予測モジュール１１０用の入力値としてコンピュータシステム１００自体において決定される。
【００２１】
シミュレーション結果は、入力／出力モジュール１１２を介して例えばデータバンクに入力され、これによりシミュレーション結果を後で評価できるようにする。
【００２２】
図２は対応するフローチャートを示す。ステップ２００で、調査すべき物質の物理化学特性を入力する。物理化学特性は、例えば、物質の脂肪親和性や分子量である。
【００２３】
ステップ２０２で、物質に対するシミュレーションモデルの物質依存パラメータのパラメータ値を、評価すべき物質の物理化学特性をベースにしたデータベースを基礎として計算する。このデータベースには、種々のテスト物質に対して予め実験により決定したパラメータ値が含まれる。この計算はステップ２０２で実行される。
【００２４】
ステップ２０４で、ステップ２０２で計算されたパラメータをシミュレーションモデルに入力する。続いて、シミュレーションモデルにおいて、昆虫の区画内での調査すべき物質の濃度／時間プロフィールの計算を行う。ステップ２０６で、この濃度／時間プロフィールを出力し、評価が可能となる。
【００２５】
パラメータ値の予測用にできる限り意味のあるデータベースを決定するために、実験的にデータベースを決定するのに用いられるテスト物質は、化学的にできるだけ多様であるのが有利である。パラメータ値を予測すべき調査対象となる新しい物質は、テスト物質により表されるテスト空間内にあるのが好ましい。
【００２６】
図３は、例として、イモムシ用の生理学的モデルを示す。この生理学的モデルは、Greenwoodらにより与えられるモデルに基づいている（上の記載参照）。
【００２７】
モデルはイモムシの次の８つの区画を含む。血リンパ（ｈｌ）３００、脂肪体（ｆｂ）３０２、筋肉（ｍｓ）３０４、クチクラ表面（ｃｓ）３０６、クチクラ（ｃ）３０８、内臓壁（ｇｗ）３１０、内臓中身（ｇｃ）３１２、神経束（ｎｃ）。評価すべき物質は、クチクラ３０８を介して局所的、あるいは、内臓を介した経口または直接注射を介して血リンパ３００内に投与される。
【００２８】
血リンパ３００は、哺乳類の血と同様、種々の器官の間の主要輸送相(phase)として機能する。
【００２９】
代謝は、クチクラ３０８でレートｋｃ、血リンパ３００でレートｋｈｌ、内蔵壁３１０でレートｋｇｗで行われる。排泄はレートｋｇｃで行われる。
【００３０】
血リンパ３００は、区画により記述されるイモムシ器官と接触状態にある自由に循環する液体としてモデル化される。区画間の輸送プロセスは、受動拡散により、レートを制限するステップとしての膜の透過とともに行われる。区画間の質量輸送のレート係数λｘは、透過率と表面積の積ＰｘＡｘおよび器官の容積Ｖｘにより
2004238390.tif 000003
として決定される。ここで、
Ｐｘ    物質に対する器官Ｘの膜（壁）の透過率
Ａｘ    器官Ｘの膜（壁）の表面積
ｃ    クチクラ
ｍｕ  筋肉
ｆｂ  脂肪体
ｎｃ  神経束
ｇｗ  内蔵壁
【００３１】
この場合、平衡状態は、ｔｘ>>Ｖｘ／（ＰｘＡｘ）経過後に到達する。この場合、平衡状態での周辺区画の濃度の比は、Ｋｘ／ｈｌとして定義される血リンパに関連する分配係数Ｋｘにより
2004238390.tif 000004
として決定される。したがって、Ｋｘは、平衡状態での血リンパと器官ｘとの間の物質の分配係数である。
【００３２】
加えて、クチクラ３０６の表面とクチクラの間の分配係数（Ｋｃ／ｃｓ）、および、内蔵壁３１０と内臓中身の間の分配係数（Ｋｇｃ／ｇｗ）は、既知である必要がある。
【００３３】
Greenwoodらによる８つの区画を有する既知モデルは、
種々の投与形態（血リンパに経口、局所、あるいは注射により行う）
食物通路を介した食物の「オープンループ」質量輸送
幼虫の成長による器官容積の時間依存性
を考慮して修正されている。
【００３４】
図３の生物物理モデルに基づいて、イモムシの各器官ｘに対し、質量平衡関係を、微分方程式により記述できる。このような質量平衡関係は、付録において、イモムシの器官に対し与えられている（式Ａ１〜Ａ８）。
【００３５】
イモムシの器官の質量平衡関係を使って薬物動態プロフィールを計算するには、計２３個のパラメータ値が既知である必要がある。それらは、
器官容積Ｖｘ（８つのパラメータ値、および、幼虫成長によるそれらの時間変動を記述するための他のパラメータ値）
血リンパ３００と周辺区画（ｃ，ｍｕ，ｆｂ，ｎｃ，ｇｗ）との間の物質の分配係数Ｋｘ
クチクラ表面とクチクラの間の分配係数、および、内臓壁と内臓中身の間の分配係数
評価すべき物質が代謝している場合におけるクチクラ、血リンパ、内臓壁の代謝レート定数
区画（ｃ，ｍｕ，ｆｂ，ｎｃ，ｇｗ）間の質量輸送に対するレート係数λｘ＝［ＰｘＡｘ］／Ｖｘ
である。
【００３６】
分配係数、および、透過率−表面積の積したがってレート係数は、昆虫の生物学的パラメータ、および、評価すべき物質の物理化学特性に依存する。代謝レート定数も物質に特異的である。
【００３７】
器官容積は公知の方法により実験的に決定できる。血リンパと各器官ｘの間の物質の分配係数Ｋｘは、次のように計算できる。
2004238390.tif 000005
ここで、
Ｋfat＝平衡状態での水と脂肪の間の物質の分配係数（脂肪親和性）
Ｋprotein＝平衡状態での水とタンパク質の間の物質の分配係数
ｆwater,x＝区画ｘでの水の容量分率
ｆfat,x＝区画ｘでの脂肪の容量分率
ｆprotein,x＝区画ｘでのタンパク質の容量分率
【００３８】
物質の膜親和力（ＭＡ）あるいは代わりにオクタノール／水分配係数（Ｋｏ／ｗ）は、分配係数Ｋfatの概算値として用いることができる。分配係数Ｋproteinは、例えば、人間の血清アルブミンの結合定数（ＫｄＨＳＡ［ｍｍｏｌ］）とタンパク質の分子量（６５ｋＤａ）とから
2004238390.tif 000006
と決定してもよい。
【００３９】
式３は、脂肪の分率およびタンパク質の分率を組み合わせることで一つの器官分率で
2004238390.tif 000007
と簡単化できる。
【００４０】
区画間の質量輸送のレート定数は、レート定数を実験的な薬物動態データとマッチングし、テスト物質の物理化学特性と相関させることで、これらテスト物質に対し実験的に決定される。
【００４１】
近似的には、物質に対する透過係数Ｐは全ての器官ｘで同一、すなわち全ての器官ｘに対してＰｘ＝Ｐであると仮定できる。但し、クチクラを除く。さらに、Ｐは、物質の脂肪親和性（Ｋｆａｔ＝ＭＡまたはＫｆａｔ＝Ｋｏ／ｗ、すなわちオクタノール／水分配係数）と膜拡散係数Ｄｍｅｍとに比例している。
2004238390.tif 000008
【００４２】
膜拡散係数の物質の分子量（ＭＷ）に対する依存度は、指数関係により
2004238390.tif 000009
と記述される。
【００４３】
式６，７を結合して
2004238390.tif 000010
とできる。ここで、αは定数である。等式１を用いて、区画間の質量輸送のレート定数は、以下のように表現できる。
2004238390.tif 000011
【００４４】
この対数をとることにより
2004238390.tif 000012
が得られる。
【００４５】
図４は、種々の区画に関して複数のテスト物質の分子量ＭＷに対するλｘ／Ｋｆａｔの２重対数を示す。図４はまた、測定点から直線回帰を用いて得られた最適適合直線を示す。
【００４６】
したがって、直線回帰により、実験的に決定されたデータベースから、傾斜ｓｍｅｍおよび切片αＡｘ／Ｖｘを決定できる。これにより、新しい物質、しかも仮想の物質であっても、脂肪親和性（Ｋｆａｔ）および分子量ＭＷのみに基づいて、量λｘを計算することのできる計算規則が得られる。式８，９に基づいてＰ・Ａｘが得られる。この物質依存パラメータ、並びに、平衡分配係数Ｋｘすなわち式（３）と（５）のＫｃ／ｃｓおよびＫｇｃ／ｇｗに基づいて、付録の等式系を例えば数値的に解くことができ、個々の区画における調査すべき物質の時間／濃度プロフィールが得られる。
【００４７】
（付録）
イモムシに対する薬物動態シミュレーションモデル用等式系
2004238390.tif 000013
2004238390.tif 000014
2004238390.tif 000015
2004238390.tif 000016
2004238390.tif 000017
2004238390.tif 000018
2004238390.tif 000019
2004238390.tif 000020
Ｃｘ：区画ｘ内の物質の濃度
Ｖｘ：区画ｘの容量；容量は、実験的に決定された関数にしたがって時間変化してもよい（Ｖｘ（ｔ））
［ＰＡ］ｘ：区画ｘの透過率と表面積の積
Ｋｘ／ｙ：区画ｘとｙとの間の分配係数
2004238390.tif 000021
Ｋｘ：代謝に対するレート定数（ｘ＝ｃ，ｈｌ，ｇｗ）
【図面の簡単な説明】
【００４８】
【図１】本発明に係るコンピュータシステムの一実施形態のブロック図を示す。
【図２】本発明に係る方法の一実施形態のフローチャートを示す。
【図３】本発明に係るシミュレーションモデル用のベースとしての生理学に基づくイモムシのモデルを、物質依存のパラメータとともに示す。
【図４】計算関数を決定するために、直線回帰を実行するための図３のシミュレーションモデルに対して得られた新たなデータベースを示す。
【符号の説明】
【００４９】
１００  コンピュータシステム
１０２  シミュレーションモデル
１０４  生理学的定数
１０６  物質依存パラメータ
１０８  データベース
１１０  予測モジュール
１１２  入力／出力モジュール
３００  血リンパ（ｈｌ）
３０２  脂肪体（ｆｂ）
３０４  筋肉（ｍｕ）
３０６  クチクラの表面（ｃｓ）
３０８  クチクラ（ｃ）
３１０  内蔵壁（ｇｗ）
３１２  内臓中身（ｇｃ）
３１４  神経束（ｎｃ）
【図１】
2004238390.tif 000025
【図２】
2004238390.tif 000026
【図３】
2004238390.tif 000027
【図４】
2004238390.tif 000028
【外国語明細書】
2004238390.tif 000029
2004238390.tif 000030
2004238390.tif 000031
2004238390.tif 000032
2004238390.tif 000033
2004238390.tif 000034
2004238390.tif 000035
2004238390.tif 000036
2004238390.tif 000037
2004238390.tif 000038
2004238390.tif 000039
2004238390.tif 000040
2004238390.tif 000041
2004238390.tif 000042
2004238390.tif 000043
2004238390.tif 000044
2004238390.tif 000045
2004238390.tif 000046
2004238390.tif 000047
2004238390.tif 000048