【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は薬物、基質等の体内動態解析方法に関するもので、特に代謝経路の律速段階となる代謝プールの大きさはもとより、代謝容量、代謝速度に関する定量的情報を得ることができる、薬物、基質等の動態解析方法に関する。

【０００２】

【従来の技術】炭素、水素、硫黄、リン等の適当な放射性同位体で標識した薬物、毒物、あるいは基質を静脈注射、腹腔内注射、経口投与等の方法で動物に投与し、血液中の放射能濃度の時間的変化を測定することにより、
動物体内での代謝を解析することができる。 あるいは、
窒素等の適当な安定（非放射性）同位体で標識した薬物等を静脈注射等の方法で動物に投与し、血液中の標識濃度の時間的変化を測定することにより、動物体内での代謝を解析することができる。

【０００３】放射性標識体を静脈注射したとき、血液中の放射能濃度の時間的変化は多くの場合、投与後の時間の真数をグラフの横軸に、血液中放射能濃度の対数を縦軸にプロットしたとき、近似的に、直線または２個もしくは３個以上の直線部分を含む曲線で示されることを利用して、解析されていた。

【０００４】前者、すなわち血液中放射能濃度の時間的変化が片対数グラフ上で直線に近似される場合には、血液中放射能濃度Cpの時間的変化は Cp＝Co exp(-λt) で表され（t は投与後の時間、Co，λは定数）、標識体の代謝に関与する代謝プールは存在しないと推定される。 これに対して、グラフが２つ以上の直線部分から成る場合には、血液中放射能濃度Cpの時間的変化は Cp＝C ₁ exp(-λ ₁ t) ＋ C ₂ exp(-λ ₂ t) ＋ C ₃ exp(-λ ₃ t)＋・・・ (1) で表され（C ₁ , C ₂ , C ₃ ；λ ₁ ,λ ₂ ,λ ₃等は定数）、何らかの代謝プールの存在が推定される。

【０００５】対応するCp，t の値の組合せから、C ₁ ,
C ₂ , C ₃等およびλ ₁ ,λ ₂ ,λ ₃等のパラメータの最適値を求める演算プログラムも市販されている。

【０００６】しかし、グラフ上でｔが充分小さい領域（C ₁ , λ ₁に関して）およびｔが充分大きい領域（C ₂ ,
λ ₂等に関して）からC ₁ , C ₂ , C ₃等およびλ ₁ ,λ ₂ ,λ ₃
等の定数（パラメータ）が精度よく求められるのは、λ
₁がλ ₂等より充分大きい場合、あるいはλ ₂がλ ₃等より充分大きい場合に限られ、λ ₁とλ ₂等との差が比較的小さい場合には直線部分の占める割合が小さくなり、グラフから求められるC ₁ , C ₂およびλ ₂等の精度が悪く、得られたグラフの実測データとの適合が極めて悪い。 このような場合には、グラフを用いないで各パラメータを算出する演算プログラムによる計算の精度も極めて悪くなる。

【０００７】C ₁ , C ₂等およびλ ₁ ,λ ₂等のパラメータの精度が比較的良好な場合でも、それらのパラメータの値を、推定される代謝モデル中のプールの容量や血液の流速等と直接結びつけることは困難であった。

【０００８】また、従来の関数(1) すなわち C ₁ exp(-λ
₁ t), C ₂ exp(-λ ₂ t) 等の単純な指数関数の和で表す方法が適用できるのは、血液中標識濃度Cpが常に減少する場合に事実上限定されており、減少の勾配が途中で増加する場合、特に濃度Cpの時間的変化が顕著な極大を示す場合には、その適用が困難であった。

【０００９】

【発明が解決しようとする課題】それ故、血液等の体液中での薬物等およびその代謝物の濃度の時間的変化を、
比較的簡単な形の関数で精度よく表現することができ、
代謝モデルの推定や、その中の代謝プールの大きさの半定量的な推定を容易にする、薬物、毒物、または基質の動態解析方法の実現が切に望まれている。 特に、濃度の減少の勾配が途中で増加したり、顕著な極大を示すなど、濃度の時間的変化が複雑なプロフィルを有する場合にも、濃度の時間的変化を比較的簡単な形の関数で精度よく表現し、薬物等の体内動態を解析する方法が切に望まれている。

【００１０】本発明の目的は、第一に、血液等の体液中での薬物等およびそれらの代謝物の濃度の時間的変化を、複雑なプロフィルを有する場合でも比較的簡単な形の関数で精度よく表現できる、薬物等の体内動態解析方法を実現することにある。

【００１１】本発明の目的は、第二に、薬物等の生体内における代謝モデルの推定や、その中の代謝プールの大きさの半定量的な推定を容易にする、薬物等の体内動態解析方法を実現することにある。

【００１２】

【課題を解決するための手段】本発明では、上記目的を達成するため、薬物等の標識体等の体液中での濃度Cpの時間的変化を関数 y＝ A[exp(-k ₁ t)＋b ₀ {1−exp(-k ₁ t)}{1−exp(-k ₂ t')}exp(-k ₃ t)＋ Σb _j f _j [ {1−exp(-k _j2 t' _j )} exp(-k _j3 t)] (2) またはこれと等価の関数で、近似的に表す（t は投与後の時間を、 A，b ₀ ，k ₁ ，k ₂ ，k ₃は正の定数を、j は１または２以上の連続した有限の整数を、b _j ，k _j2 ，k _j3はそれぞれのj に対して正の定数を、f _jは{1−exp(-k
_j1 t)} 又は１を、k _j1はそれぞれのj に対して正の定数を、それぞれ表し、t'はt がd より大きいときt−d に等しく、t がd より大きくないとき０である変数を示す。 d は０または正の定数を示す。 t' _jはt がd _jより大きいとき t−d _jに等しく、t がd _jより大きくないとき０
である変数を示す。 d _jは０または正の定数を示す。 b ₀およびb _jは１を超えることはない。 ）。 A, b ₀ , k ₁ , k ₂ ,
k ₃ , b _j , k _j2 , k _j3 , k _j1等の定数、すなわち各パラメータは、次のようにして求められる。

【００１３】投与後の時間t の真数をグラフの横軸に、
体液中の標識体等の濃度Cpの対数を縦軸にプロットしたとき、t=0 における曲線の切片からA を、t ＝0 付近において曲線の接線の勾配からk ₁を、曲線の接線の勾配が実質的に一定となる部分の接線の勾配からk ₃を、それぞれ近似値として求める。 仮に k ₃ ＝k ₂ ，k ₃ ＝2k ₂ ，k ₃ ＝3k ₂ ，k ₃ ＝4k ₂ ，2k ₃ ＝k ₂ ， 3k ₃ ＝
k ₂ ， 4k ₃ ＝k ₂等とした場合について、それぞれ L＝{1−exp(-k ₂ t)} exp(-k ₃ t) が極大となるt またはそれに近い時点t （t _mとする）において {Cp／A −exp(-k ₁ t _m )} ／{1−exp(-k1t _m )} {1−exp(-k ₂
t _m )} exp(-k ₃ t _m ) の値を求め（A ，k ₁ ，k ₃としては前述の近似値を用いる）、それをb ₀の仮の近似値として、 A [exp(-k ₁ t)＋b ₀ {1−exp(-k ₁ t)}{1−exp(-k ₂ t)} exp(-
k ₃ t)] の値がCp（実測値）に最も近い値を与えるようなk ₂の値（k ₃との相対関係）を見出し、このk ₂を用いて Y ₀ ＝ Ay ₀ ＝A [exp(-k ₁ t)＋b ₀ {1−exp(-k ₁ t)}{1−exp(-k ₂ t)}exp(-k ₃ t)] (2) を体液中の標識体等の濃度Cpの時間的変化を表す関数の第一および第二の項とする。 t の代りに、t が正の定数
d より大きいとき t−d に等しく、t がd より大きくないとき０である変数t'を用いると、より精度が上がる場合がある。 この場合には、適当と思われるd の値をいくつか選んで、それらのうちY0の値がCp（実測値）に最も近い値を与えるようなものを見出す。

【００１４】{1−exp(-k ₁ t _m )} が１に近似であるとして、すなわち、式 {Cp／A −exp(-k ₁ t _m )} ／{1−exp(-k ₂ tm)} exp(-k ₃ t _m )
または{Cp／A}／{1−exp(-k ₂ t _m )} exp(-k ₃ t _m ) を用いてb ₀の仮の値を求めてもよい。 あるいは、k ₃の近似値を求めるに用いたtの領域においてCp／Aexp(-k ₃ t)
の値を求め、これをb ₀の仮の値として用いてもよい。

【００１５】関数Y ₀は、 Y ₀ ＝Aexp(-k ₁ t)＋B ₀ {1−exp(-k ₁ t)}{1−exp(-k ₂ t)} exp
(-k ₃ t) と表すこともできるが、b ₀ ＝B ₀ /Aとして(2) 式のように表すほうが計算上簡便である。 b ₀は１を超えることはない。

【００１６】k ₂ ＝mk ₃ (m は正の実数) の関係があるとき、L ＝{1−exp(-k ₂ t)} exp(-k ₃ t)が極大となるt _mは t _m ＝ {ln(m+1)}／mk ₃である。 t _mにおいて得られるL の極大値は {1−1／(m+1)} ／(m+1)1/m である。 例えば、k ₂ ＝2k ₃であるとき、t _m ＝ ln3／(2
k ₃ ) であり、L の極大値は(2/3)(3)-1/2、、すなわち約
1／2.6 である。 k ₂ ＝k ₃であるとき、t _m ＝ln2 ／k ₃であり、L の極大は1/4 である（m が１より大きいとき、L
の極大値はこれより常に大きい）。

【００１７】k ₃ ＝nk ₂ （ nは正の実数）の関係があるとき、{1−exp(-k ₂ t)} exp(-k ₃ t) が極大となるt _mは t _m ＝{ln(n+1)−ln(n)}／k ₂ ＝{ln(n+1)−ln(n)}／(k ₃ /n) である。 t _mにおいて得られるL の極大値は [1−{n／(n
+1)}] {n／(n+1)}n である。 例えばk ₃ ＝4k ₂であれば、L の極大値は(1/5)(4/5)4 、すなわち約0.08である（n が１より大きいとき、L の極大値はk ₂ ＝k ₃であるときの極大値1/4 より常に小さい）。

【００１８】このような関係を利用して、例えば、まず
k ₂ ＝k ₃と仮定し、 {Cp/A − exp(-k ₁ t _m )}／[{ 1−exp(-k ₁ t _m )}/4] ＝4{Cp/A−exp(-k ₁ ln2/k ₃ ) }／{1−exp(-k ₁ ln2/k ₃ )} を仮のb ₀の値として、関数Y ₀すなわち y ₀ ＝A [exp(-k ₁ t)＋b{ 1−exp(-k ₁ t)}{1−exp(-k ₂ t)} e
xp(-k ₃ t)] の値を求め、次いでk ₂ ＝2k ₃と仮定して、 2.6 {Cp/A − exp(-k ₁ t _m )}／{1−exp(-k ₁ t _m )} t _m ＝
ln3/(2k ₃ ) を仮のb ₀の値として、関数Y ₀の値を求め、Cpの実測値との差異が減少するか増大するか、判定する。 もし実測値との差異｜Y ₀ −Cp｜が全体として減少するなら、m＝k ₂ /
k ₃の値を３または４、あるいはそれより大として、｜Y
₀ −Cp｜が最も小さくなる m＝k ₂ ／k ₃の値を見つける。

【００１９】k ₂ ＝2k ₃としたとき、k ₂ ＝k ₃である場合より実測値との差異｜Y ₀ −Cp｜が増大するなら、k ₃ ＝nk ₂
(n＝2,3,4,等）として、上と同様の手順で最適のn を見つける（m,n は整数でなくてもよい）。

【００２０】この最適のm またはn の値を用いて得られるb ₀の値に対して、関数Y ₀の値とCpの実測値の差異（標準偏差の和として表す）を求め、b ₀およびk ₂の値、さらに必要ならk ₁ ，k ₃の値を調節して、最適の関数Y ₀を決定する。 関数Y ₀の増加する領域において、k ₂の値を調節しても関数Y ₀の増加率とCpの実測値の増加率の差異が比較的大きい場合には、t の代りに、t が正の定数d より大きいとき t−d に等しく、t がd より大きくないとき０
である変数t'を用い、適当なd の値を選択することにより、近似の精度を改善することができる。

【００２１】y _j ＝b _j f _j [ {1−exp(-k _j2 t' _j )}exp(-k _j3 t)]
(j＝1,2,3,・・・) についても、同様の要領で、 y ₁ ＝b ₁ f ₁ [ {1−exp(-k ₁₂ t' ₁ )}exp(-k ₁₃ t)]， y ₂ ＝b ₂ f ₂ [ {1−exp(-k ₂₂ t' ₂ )}exp(-k ₂₃ t)]， y ₃ ＝b ₃ f ₃ [ {1−exp(-k ₃₂ t' ₃ )}exp(-k ₃₃ t)]， 等の各項について、各々が支配的な（主として寄与する）t の領域（例えば、y _jの極大付近）について、それらとCp/Aとの差を最も小さくする各パラメータを求めることができる。 各項の寄与が重なり合うt の領域についても、各項の和と実測値との差が最も小さくなるように、各パラメータを調節する。 この調節は、適当な隔たりで選ばれた複数のt の値に対して各項の値を計算して、計算表を作成することにより、困難なく行うことができる。

【００２２】このようにして、関数(2) Y＝ A[exp(-k ₁ t)＋b ₀ {1−exp(-k ₁ t)}{1−exp(-k ₂ t')}ex
p(-k ₃ t)＋Σb _j f _j [ {1−exp(-k _j2 t' _j )} exp(-k _j3 t)] の各パラメータを決定し、関数(2) を具体的に決定することができる。

【００２３】上記の計算は電子計算機を用いると容易に行うことができる。 電子計算機は、例えば、A について有効数字１から9.999 まで、10進桁0.01から1000まで、
b ₀ ，b _j ，k ₁ ，k ₂ ，k ₃ ， d，k _j1 , k _j2 , k _j3 , d _jについて、有効数字１から9.9 又は 9.99まで、10進桁0.001から10までの数値信号を演算部へ選択的に出力できるようにしておけば、目的を達する。

【００２４】本発明の代謝解析方法は、静脈注射、腹腔内注射、経口投与等、種々の投与方法に対して適用できるが、特に静脈（門脈を含む）注射により投与する場合に有用である。 静脈注射により放射能標識体を投与し、
血液中の放射能濃度の時間的変化を測定するとき、本発明で用いる関数は血液中の放射能濃度の時間的変化の実測値に極めてよく一致する。 放射能標識体の代りに安定同位体で標識された化合物を用いても同様である。 標識体を用いず、化合物自体およびその代謝物の少なくとも一部を追跡する場合（液体クロマトグラフィ、ガスクロマトグラフィを等の微量分析法を用いて）にも、有効である。

【００２５】本発明の代謝解析方法は、新薬の開発においてその薬物動態を明らかにする上で極めて有用である。

【００２６】

【実施例】以下に実施例を示し、本発明のさらに具体的な説明とする。 ［実施例１］6.03 MBq／mgの¹⁴ Ｃで１位に標識したD-th
reo-[dichloroacetyl-1- ¹⁴ C]CHLORAMPHENICOL 55マイクロク゛ラム（0.33MBq ）を含む溶液0.18マイクロリットルを体重 180g のウィスター系雄ラットの尾静脈に注射し、１分、５分、
15分、30分、45分、１時間、２時間、４時間、６時間、
８時間、24時間、48時間、72時間、96時間後にそれぞれ尾静脈から約 100マイクロリットルずつ採血し、その放射能を液体シンチレーション計数法により測定した。 動物３個体について平均した測定結果を表１に示す。

【００２７】

【００２８】投与後の時間t の真数を横軸に、相対放射能強度Cpの対数を縦軸にプロットしたグラフを作成し、
t＝１分（0.0167時間）、５分（0.0833時間）、15分（0.25時間）に対応する点の最も近傍を通る直線を引き、その勾配および縦軸(t＝0)の切片を求めた。 勾配は−12、切片は 187であった。 またグラフ上で t＝30，4
5，60分の各点の最も近傍を通る直線の勾配を求めたところ、約3.0 であった。

【００２９】そこで仮に前記関数Y におけるA として 1
87を、k ₁として12.0を、k ₃として3.0 を、それぞれ選び、次にk ₂の値を決定する。 k ₂ ＝mk ₃におけるm または
nk ₂ ＝k ₃におけるn を選ぶために、L すなわち{1−exp
(-k ₂ t)} exp(-k ₃ t) の、極大値を与えるt _mをまず推定する。 k ₃ ＝3.0 であるとき、t _mはm との間に表２に示す関係を、n との間に表３に示す関係を、それぞれ有する。

【００３０】表 ２ m {ln(m+1)}／m ＝k ₃ t _m t _m 1 ln2 ＝0.693 0.231 2 ln3／2 ＝0.549 0.183 3 ln4／3 ＝0.462 0.154 4 ln5／4 ＝0.402 0.134

【００３１】表 ３ nn{ln(n+1)−ln(n)}＝k ₃ t _m t _m 2 2{ln3−ln2}＝0.811 0.270 3 3{ln4−ln3}＝0.863 0.288 4 4{ln5−ln4}＝0.893 0.298

【００３２】グラフ上の曲線の形状から、t _mは15分付近にあると推定し、仮にm を１とすればk ₂ ＝mk ₃の値は３、n を２とすればk ₂ ＝k ₃ ／n の値は1.5 となる。 t _m ＝
0.231または0.270 においては勿論、時点30分ないし45
分においてexp(-k ₁ t _m )の値は０に近く、この時間領域において Cp/A −exp(-k ₁ t)の値はCp/A とみなしてよいので、30分および45分における Cp/A の値0.15および0.10
0 と L＝{1−exp(-k ₂ t)} exp(-k ₃ t)＝{1−exp(-k ₂ t)}ex
p(-3.0t) （k ₂ ＝1.5 又は3.0)の値の比較により、b ₀の値を仮に1.0 に選ぶ。 A として 187を、k ₁として12.0
を、b ₀として1.0 、k ₃として3.0 、k ₂として1.5 又は3.
0 をそれぞれ用いて、 Y ₀ ＝ Ay ₀ ＝ A[exp(-k ₁ t)＋b ₀ {1−exp(-k ₁ t)}{1−exp(-k
₂ t)} exp(-k ₃ t) の値を計算すると、時点５分ないし１時間において表４
の通りとなる。

【００３３】

【００３４】５分の計算値から見るとk ₁として12.0は小に過ぎ、また15分ないし30分の計算値から見るとk ₂として3.0 は過大で、1.5 はやや小に過ぎることがわかる。
k ₁およびk ₂の値を調節して、k ₁として13.7、k ₂として1.
7 が計算値と実測値の差を最も小さくする値であった。
k ₁を13.7としたとき、15分から45分の領域で上記の関数
Y ₀の値はCpの実測値に比して低いので、 Aの値を 195に修正した。 これによる１分でのY ₀の値の増加は、実測値の誤差から見て許容範囲にある。

【００３５】２時間から９６時間の時間領域での相対放射能強度Cpは、８ないし24時間の範囲に極大を有すると見られるので、この領域に主として寄与する関数 y ₁ ＝b ₁ {1−exp(-k ₁₁ t)} {1−exp(-k ₁₂ t)} exp(-k ₁₃ t) を想定し、その各パラメータを求めた。 k ₁₁はk ₃に等しいと仮定した。 48時間と72時間の間の放射能強度Cpの減衰から、k ₁₃は0.009 と推定された。 {1−exp(-k ₁₂ t)}
exp(-k ₁₃ t)の極大を与えるt _mは、m との間に表５に示す関係を有する。

【００３６】

【００３７】m＝30を選ぶと、k ₁₂ (＝mk ₁₃ ) は0.27となる。 ２時間以上の時間領域においてCp/A −exp(-k ₁ t)
の値は Cp/A とみなしてよいので、 Cp/A の値と L＝{1−exp(-k ₁₂ t)} exp(-k ₁₃ t)＝{1−exp(-0.27t)} ex
p(-0.009t) の値の比較により、b ₁の値を仮に0.12に選んだ。 しかし
k ₁₃ ＝0.009、k ₁₂ ＝0.27の組合せで与えられる関数 Y＝
A(y ₀ ＋y ₁ ) の値は２時間ないし６時間の領域で誤差がやや大きいので、 k ₁₂ ＝0.28に修正して誤差をさらに減少させた。

【００３８】結局、Cpの時間的変化を近似的に表す関数
Yとして Y＝ 187[exp(-13.7t)＋{1−exp(-13.7t)}{1−exp(-1.7
t)}exp(-3.0t)＋0.12{1−exp(-3.0t)} {1−exp(-0.028
t)}exp(-0.009t)] が得られた。 各測定時間における Yの値とCpの実測値とを、表６に比較して示した。 また、図１にグラフとして示した。 図１には、投与後の時間t の真数が横軸に、相対放射能強度Cpの対数が縦軸に示されている。

【００３９】

【００４０】［実施例２］実施例１と同様に、 ¹⁴ Ｃで１
位に標識したD-threo-[dichloroacetyl-1- ¹⁴ C]CHLORAMP
HENICOL の溶液をウィスター系雄ラットの尾静脈に注射し、１分、５分、15分、30分、45分、１時間、２時間、
４時間、６時間、８時間、24時間、48時間、72時間、96
時間後にそれぞれ尾静脈から採取した血液から、常法により血漿を得て、その放射能を液体シンチレーション計数法により測定した。 動物３個体について平均した測定結果を表７に示す。

【００４１】

【００４２】実施例１と同様の関数を想定し、１分から１時間の領域で実施例１と同様の方法で A，k ₁ ，k ₂ ，
b ₀ ，k ₃の値を選定し、48時間から96時間の領域について
k ₁₃の値を選定した。 Aとして 146を、k ₁として12、k ₂
として 6.8、b ₀として0.25、k ₃として1.12、 k ₁₃として
0.0095を選んだ。 b ₁の値は、48時間以上においてexp(-k
₁₂ t)の値が１に比して無視し得るとして、 Cp/A の値と
exp(-k ₁₃ t)の値の比較から0.047 とした。 k ₁₂は0.2 以下でないと、１時間での Yの値が過大となる。 y ₁における{1−exp(-k ₁₁ t)} の項は１に近似した。

【００４３】k ₁₂の値を0.2 以下とすると、４ないし８
時間の領域で Yの値は４にも満たないので、第４の項として y ₂ ＝b ₂ {1−exp(-k ₂₁ t)} {1−exp(-k ₂₂ t')}exp(-k ₂₃ t) t'＝ t−d を付加することにした。 Cp/A の値とy ₀ ＋y ₁の値の差を求め、この差を関数y ₂によって表すものとして、最適の各パラメータb ₁ ，k ₂₂ ，k ₂₃ ，d を選定した。 k ₂₁はk ₃
と等しい1.12とした。 まず d＝０として、y ₀ ，y ₁の場合と同様の手順でk ₂₃ ，k ₂₂ ，b ₂を順次求めた。 Cp/A− (y
₀ ＋y ₁ ) の値の６時間と２４時間の間での変化から、
k ₂₃は約 0.1が適当と見られた。 しかし、パラメータを様々に変えても、１ないし２時間の Yの値が相対的に過大となる。 そこでt'＝ t−d において d＝１，1.4 ，２
等とし、t'を用いて{1−exp(-k ₂₂ t')}exp(-k ₂₃ t)が極大となるt _mが６時間付近になるようにk ₂₂ ，k ₂₃の組合せを選び、その結果各パラメータの最適値として d ＝1.4 ，k ₂₂ ＝0.26，k ₂₃ ＝0.12，b ₂ ＝0.24 を得た。

【００４４】前記関数 Yは A (y ₀ ＋y ₁ ＋y ₂ ) に相当するから、血液中の濃度を表す関数として Y＝A(y ₀ ＋y ₁ ＋y ₂ ) ＝ A[exp(-k ₁ t) ＋b ₀ {1−exp(-k ₁ t)}{1−exp(-k ₂ t)} exp(-k ₃ t) ＋b ₁ {1−exp(-k ₁₁ t)} {1−exp(-k ₁₂ t)} exp(-k ₁₃ t) ＋b ₂ {1−exp(-k ₂₁ t)} {1−exp(-k ₂₂ t')}exp(-k ₂₃ t)] t'＝ t−d ＝ 146[exp(-12t) ＋0.25{1−exp(-12t)}{1−exp(-6.8t)} exp(-1.12t) ＋0.047 {1−exp(-0.2t)} exp(-0.0095t) ＋0.24{1−exp(-1.12t)}{1−exp(-0.26t')}exp(-0.12t)] t'＝ t−1.4 が得られた。 各測定時間における Yの値とCpの実測値とを、表８に比較して示した。

【００４５】

【００４６】

【発明の効果】本発明の薬物、毒物、または基質の動態解析方法によれば、血液等の体液中での標識体等の濃度の時間的変化を、比較的簡単な形の関数で精度よく表現することができる。 特に、極大を示す等複雑なプロフィルを有する変化は、従来のコンパートメントモデルでは表現が困難であったけれども、本発明によると、複雑な計算を行うことなく、比較的簡単な形の関数で精度よく表現できる。

【００４７】また本発明によれば、薬物、毒物、基質等の生体内における代謝モデルを推定し、代謝プールの大きさ等を少なくとも半定量的に推定することが容易になる。 さらに、代謝プールの大きさと体液の流量や組織の摂取速度との相対関係も推定でき、一方についての情報があれば他方について推定することができる。

【図面の簡単な説明】

【図１】 放射能強度の時間変化を示す関数及び実測値を示すグラフである。

─────────────────────────────────────────────────────

【手続補正書】

【提出日】平成６年２月７日

【手続補正１】

【補正対象書類名】明細書

【補正対象項目名】図面の簡単な説明

【補正方法】削除

【手続補正２】

【補正対象書類名】明細書

【補正対象項目名】発明の詳細な説明００３８

【補正方法】変更

【補正内容】

【００３８】結局、Cpの時間的変化を近似的に表す関数
Yとして Y＝ 187[exp(-13.7t)＋{1−exp(-13.7t)}{1−exp(-1.7
t)}exp(-3.0t)＋0.12{1−exp(-3.0t)} {1−exp(-0.028
t)}exp(-0.009t)] が得られた。 各測定時間における Yの値とCpの実測値とを、表６に比較して示した。