| 研究課題/領域番号 |
02807143
|
|---|
| 研究種目 |
一般研究(C)
|
| 配分区分 | 補助金 |
|---|
| 研究分野 |
麻酔学
|
|---|
| 研究機関 | 東京大学 |
|---|
研究代表者 |
矢島 直 東京大学, 医学部, 助手 (20191108)
|
|---|
| 研究分担者 |
花岡 一雄 東京大学, 医学部, 教授 (80010403)
菅井 直介 東京大学, 医学部, 講師 (10010563)
|
|---|
| 研究期間 (年度) |
1990 – 1991
|
| キーワード | 薬物動態学 / 薬力学 / 非脱分極性筋弛緩薬 / 最小2乗法 / 筋弛緩モニター / 単収縮高 / ディジタルフィルタ / Sheinerのモデル |
|---|
| 研究概要 |
常用量の非脱分極性弛緩薬を投与後、筋弛緩モニターのデータをオンラインでコンピュータに取り込み、患者の薬物動態学的および薬力学的パラメータ(pharmacokinetic-dynamic parameters)を推定するのに、精度がよくかつ簡便な方法を開発した。以下に示すように(1)作用発現時、(2)作用回復時に分けて各パラメータを推定する理論式を作成した。
筋弛緩薬の薬物動態学は次の2コンパートメントモデルに従い、
C(t)=Ae^<-αt>+Be^<-βt>(α>β) (1)
薬力学はSheinerのモデルに従うと仮定した。すると筋弛緩薬投与後、時刻tでの単収縮高を推定する式として次式が導かれる。
log[(l_O-l)/l]=γlogσ+γlog{(ξ(e^<-αt>-e^<-Keo t>)/(K_<eo>-α)+(e^<-βt>-e^<-Keo t>)/(K_<eo>-β))/(ξ+1)} (2)
ここで、l_Oはcontrol twitchの高さ、lは投与後t秒でのtwitchの高さ、γとk_<eo>とC^<55>_P50はpharmacodynamicsの定数である。また、ξ=A/B,δ=(k_<eo>・D)/(C<55PKP>50・Vd)である。
(1)作用発現時:Xが十分小さいときにはe^<-X>≒1-Xが成り立つので、式(1)より
log[(l_O-l)/l]=γlogδ+γlog t (3)
が成り立つ。この式より、初回投与量Dと作用発現時間tとは反比例することが示唆される。実際、たとえばベクロニウムの95%作用発現時間は、初回投与量に反比例することが、ヒトの臨床データで確認された。
(2)作用回復時間:時間tが十分に大きいときには
(e^<-αt>-e^<keo t>)/(k_<eo>-α)<<(e^<-βt>-e^<keo t>)/(k_<eo>-β) (∵α>keo>β)
が成り立つので、式(1)は次式のように変形される。
log[(l_O-l)/l]=γlogδ′-γβt+γlog(1-e^<-c・t>) (4)
ただしδ′=(k_<eo>・D)/(C^<55>_P50・Vd^<ext>・C)、Vd^<ext>=(1+ξ)Vd、c=k_<eo>-βである。
最小2乗法を用いて、これらの式に筋弛緩モニターのデータをあてはめて求めたパラメータγ、β、k_<eo>などの値は、実際に血中濃度を測定して求めた値と良く一致していた。また、初回投与時のデータから、追加投与後の筋弛緩薬の薬理学的作用をよく推定することができた。
さらに、筋弛緩薬投与後の単収縮高の経時変化のデータから、雑音成分を除去するためのF1Rディジタルフィルタを作成し、電気メスやその他の物理的外乱を取り除くプログラムを開発した。
|
