Article information

Abstract

INTRODUCTION

전신마취는 의식소실(amnesia), 진통(analgesia), 근이완(muscle relaxation), 반사소실(abolition of autonomic reflex)의 4가지 요소를 필요로 하며, 이를 충족시키기 위하여 다양한 약제들이 병합 투여되고 있다. 특히 의식소실과 진통을 위하여 투여되는 수면제(hypnotics)와 아편유사제(opioid)의 경우 치료역이 좁은 약제이기 때문에 객관적인 지표에 의하여 정밀하게 적정(titration)할 필요가 있다. 이들 약제는 중추신경계에 작용하는 약제이므로 중추신경계 대리표지자를 이용하여 약제 적정의 지표로 삼을 수 있다. 특히 수면제와 아편유사제를 “적절하게 섞어서” 사용하기 위해서는 개별 약제의 약동, 약력학적 특성을 파악하는 것이 중요하다. 그렇지만 두 약제가 동시에 투여되는 상황에서는 개별 약력학은 더 이상 유효하지 않으며, 두 약제간 약력학적 상호작용의 패턴을 미리 알고 있다면 두 약제의 적정을 보다 효과적으로 할 수 있다(Fig. 1). 여기에서는 두 약제의 약동학적인 상호작용을 제외한 약력학적 상호작용의 관계를 이소볼(isobole)과 반응표면모형(response surface model)을 이용하여 탐색해 보기로 할 것이며, 이를 객관적으로 정량화해 주는 방법 중에서 유연성이 가장 뛰어난 소프트웨어인 NONMEM (nonlinear mixed effects modeling)으로 구현하는 방법까지 소개하기로 한다.

Fig. 1

Necessity for consideration of drug interaction during general anesthesia. PK: pharmacokinetic, ke0: blood-brain equilibration rate constant, PKPD: pharmacokinetic and pharmacodynamic, CNS: central nervous system.

DRUG INTERACTION

약제의 약력학적 상호작용은 크게 3가지로 분류한다. 약제가 병용 투여되는 상황에서 나타나는 효과가 개별 약제가 나타내는 효과의 합과 동일하다면 상가적(additive), 크다면 상승적(synergistic), 작다면 상가이하적(infra-additive)이라고 정의하고, 흔히 상승적을 상가이상적(supra-additive)으로, 상가이하적을 길항적(antagonistic)이라고도 부른다[1]. 이러한 약력학적 상호작용을 설명하기 위하여 상호작용 모형을 이용하는데, 이상적인 상호작용 모형의 특성을 요약하면 다음과 같다.

ISOBOLE

Isobole 또는 isobologram은 동일 효과를 유발하는 두 약제의 조합을 나타내어 주는 동일효과곡선(iso-effect-curve)을 말하며, 이를 Fig. 2에 제시하였다[2]. Isobole은 단지 단일 효과만을 설명하기 때문에 보편적으로 모든 경우에 적용되지는 못하지만 직관적으로 약제의 상호작용을 평가하는 데는 용이하다. 흔히 이용하는 약제의 상가적 상호작용은 Loewe 상가성(Loewe additivity)으로 표시되며 이를 식으로 표시하면 다음과 같다.

Fig. 2

Isobole. (A) additivity, (B) synergy, (C) infra-additivity. D1, and D2, are isoeffective doses of two drugs administered alone. The administration of the two drugs in combination (d1, d2) results in the same effect. If D1, and D2, are the D50 doses, in each case the line represents the 50% isobole.

식(1)에서 C_x는 약제 X의 농도를 의미하고, C_50x는 최대 효과의 반을 나타낼 때의 약제 X의 농도를 의미하며, 이는 약제 Y에 있어서도 동일하게 적용된다. Loewe 상가성 이외에도 Bliss 독립(independence) 또는 중앙-효과 방법(median-effect method)과 같은 방법을 이용할 수는 있지만[3-5], 아직까지 Loewe 상가성이 보편적으로 이용되고 있고[6], 특히 반응표면모형의 Greco 모형에서 이를 받아들여 사용하고 있다[3].

RESPONSE SURFACE MODEL

Isobole이 특정 효과(예: 50% 효과)를 나타내는 두 약제의 조합을 나타낸 동일효과곡선이라면 반응표면모형은 두 약제의 용량 혹은 농도 조합에서 나타낼 수 있는 반응의 효과를 전체 효과 범위에서 제시한 모형이다. 통상 반응의 효과는 반응이 나타날 확률개념으로 표시를 하게 된다. 반응표면모형에서 50%의 효과를 나타낼 확률에 해당하는 두 약제의 동일효과곡선이 50% isobole이 된다(Fig. 3). 반응표면모형으로 약제의 약력학적 상호작용을 설명할 때 그림을 이용하면 이해하기가 용이하며, 대표적인 예를 Fig. 4에 제시하였다. Figs. 4A–C는 반응이 나타날 동일한 확률에서의 isobole을 이용하여 약제 상호작용을 구분할 수 있다. Fig. 4A의 경우 확률이 1일 때의 isoble이 직선임을 알 수 있기 때문에 이는 상가적 상호작용(additive interaction)을 보인다. Fig. 4B의 경우에는 확률이 1일 때의 isobole이 직선에서 0에 가까운 쪽으로 오목한(concave) 곡선의 형태를 띠게 되며, 이는 상승적 상호작용(synergistic interaction)을 의미한다. Fig. 4C의 경우에는 확률이 1일 때의 isobole이 직선에서 0에서 멀어지는 쪽으로 볼록한(convex) 곡선의 형태를 띠게 되므로, 이는 상가이하적 상호작용(infra-additive interaction)을 의미한다. Figs. 4D–F는 반응이 나타날 확률의 크기로 구분이 가능하다. Fig. 4D의 경우 한 약제의 용량 혹은 농도가 증가하면 S자 모양의 용량 반응 관계가 나오면서 최대효과까지 나오지만 다른 약제의 경우는 용량 혹은 농도를 아무리 올려도 반응이 나타날 확률이 1까지 나타나지 않는 것을 확인할 수 있다. 이는 작용제(agonist)와 부분작용제(partial agonist)의 상호작용을 설명한다. Fig. 4E의 경우는 다른 약제가 전혀 효과를 내지 않고 있으며, 작용제의 용량 혹은 농도가 증가하여도 효과를 억제하는 것을 관찰할 수 있으며, 이는 작용제와 경쟁적 길항제(competitive antagonist)의 상호작용을 설명한다. Fig. 4F의 경우는 원하는 효과가 아닌 반대의 효과가 나타나는 것을 확인할 수 있으며, 이 경우는 역작용제(inverse agonist)와 작용제의 상호작용을 설명한다. 이처럼 약제의 상호작용을 그림으로 그리면 직관적으로 파악하는 것은 용이할 수 있지만 두 약제의 상호작용의 관계를 정량화하기에는 한계가 있다. 약제의 상호작용을 설명하는 반응표면모형의 다양한 종류들을 이용하면 이들 관계들을 객관적으로 파악할 수 있다. 우선 이런 모형들을 다루는 소프트웨어인 NONMEM (nonlinear mixed effects modeling)을 간단히 소개하기로 한다.

Fig. 3

Response surface models showing synergistic (A), additive (B), and infra-additive (C) interactions. Each panel contains a response surface (upper) and an isobologram showing the 50% probability isobole (lower).

Fig. 4

Response surfaces models. For all figures, the effect with no drug is 0. (A) Additive interaction between two agonists, (B) supra-additive (synergistic) interaction between two agonists, (C) infra-additive interaction between two agonists, (D) partial agonist and full agonist; (E) competitive antagonist and full agonist, (F) inverse agonist and full agonist. The isoboles for 10, 20, 30, 40, 50, 60, 70, 80, and 90% response are shown.

NONMEM (NONLINEAR MIXED EFFECTS MODELING)

반응표면모형을 다룰 수 있는 약력분석 소프트웨어가 여러 가지가 있지만 가장 대표적이면서 최초인 것은 University of San Fransico (UCSF)의 내과의사인 Lewis Sheiner와 통계학자인 Stuart Beal이 창안한 NONMEM (nonlinear mixed effects modeling)이 있다. 사실 집단분석(population analysis) 의 개념도 이들에 의하여 최초로 제안된 것이다(Fig. 5). NONMEM은 비선형혼합효과모형(nonlinear mixed effects modeling)에서 유래된 말이면서, 동시에 집단분석용 표준 소프트웨어의 이름이기도 하다. 약력학은 S자모양 최대효과 모형(sigmoid Emax model)으로 설명하기 때문에 비선형이 대부분이라는 것을 알 수 있다. 선형회귀와 비선형 회귀는 추정 방식이 매우 다르기 때문에 이와 같이 비선형이라는 말을 강조하는 것이다. 혼합효과모형(mixed effects modeling)은 말 그대로 무엇인가가 혼합되어 있어 이들을 함께 모형화한다는 의미인데, 약동, 약력학적 모수를 집단의 평균값(이 값은 사람마다 달라지지 않으며 고정효과 모수라고 한다. NONMEM에서는 θ로 표기한다.)과 inter-individual random effects parameter로 나누는 방식을 말한다. 여기서 inter-individual random effects parameter는 random inter-individual variability를 기술해준다. 이러한 random inter-individual variability는 평균이 0이면서 분산이 ω²인 정규분포를 하는 것으로 가정한다. Random inter-individual variability는 생물학적으로 당연히 있을 수 있는 현상이며, 이 중에서 환자가 가진 특성(NONMEM에서는 이것을 covariate, 공변량이라고 부른다.)으로 설명 가능한 부분을 고정효과 모수 쪽으로 가지고 온다. 이것이 바로 집단분석의 개념이다.

Fig. 5

Concepts of population analysis using nonlinear mixed effects modeling. PK: pharmacokinetic, PD: pharmacodynamic, TCI: target-controlled infusion.

이와 같이 하려면 NONMEM의 제어구문(control code, NONMEM에게 회귀분석을 시키는 일종의 명령문의 집합) 내에 다음과 같이 적는다. 예를 들어 목표 효과처농도 조절주입에 필요한 혈액-뇌평형속도상수(k_e0)의 경우,

여기서 TKE0는 k_e0의 대표값(random inter-individual variability가 0인 전형적 인간, 혹은 집단 내 평균값, 여기서는 THETA(1))이며, ETA(1)은 k_e0의 random inter-individual effects parameter이다. 즉, fixed effects parameter인 k_e0이 random inter-individual effects parameter인 ETA(1)과 지수함수를 통하여 연결되었다. 어떤 방식으로 연결하는가 하는 문제는 모형을 하는 사람이 결정할 문제이다. 곱하기나 더하기로 연결할 수도 있으며, 자료를 가장 잘 설명하는 쪽으로 선택하면 된다.

위에서 ETA(1)의 일부를 환자의 연령으로 일부 설명한다고 하자. 이 경우 NONMEM 제어구문은 다음과 같이 바꿀 수 있다.

위와 같은 코드가 농도, 효과 자료를 더 잘 설명하였다면, 최종 모형은 k_e0에 연령이 다음과 같이 연관된 것이 될 것이다.

위의 수식이 TCI system 안에 들어가면, 연령에 따라서 용량을 달리할 수 있게 된다. 이와 유사하게 체중, 성별 등 어떠한 약동, 약력 모수에라도 자료를 잘 설명하기만 한다면 공변량(covariate)으로 사용하여 고정효과 모수에 수학적으로 연결시킬 수 있다. 이 모든 것을 합치면 환자가 가진 모든 특성을 용량 적정에 이용할 수 있게 된다. 사실 우리가 임상에서 약물을 투여할 때 의식적이든 무의식적이든 이와 같은 방식을 취하고 있다. 다만 단편적인 지식에 근거하여 수동으로 하던 것을 정확한 수학적 모형을 바탕으로 좀 더 편리하면서도 정밀하게 컴퓨터를 이용하는 것이 다를 뿐이다. 각종 약동, 약력 모수를 고정효과(covariate 포함)와 random inter-individual variability로 설명하고도 다 설명하지 못하는 농도, 효과 곡선 부분(즉, 관찰효과와 모형 예측의 차이; observations - model predictions)은 residual random effects parameter를 사용하여 residual random variability로 넘긴다. 이 residual random variability는 효과 측정의 부정확성, 시간 표기 오류 등 실제 오차를 반영하며, 역시 평균이 0이고 분산이 σ²이라고 가정한다. 요약하면 NONMEM을 이용한 집단분석은 모수를 고정효과(공변량 포함)와 임의 개체간변이로 나누어서 환자집단의 시간-농도 혹은 농도-효과 곡선을 설명하도록 해 준다.

물론 Matlab (MathWorks, Inc., Natick, MA, USA) 또는 R (R Foundation for Statistical Computing, Vienna, Austria)과 같은 프로그램으로 반응표면모형 모수를 추정할 수 있다. 즉, 이 두 프로그램을 이용하여 아래에 제시한 각 반응표면모형의 구조모형을 이용한 제어구문을 생성할 수 있다. 그렇지만 아직까지 Matlab이나 R에서는 혼합모형에 대한 유연한 library를 제공하지 않기 때문에 모형 모수의 추정값과 그 표준편차 밖에 구하지 못하는 단점이 있다. 즉, 개인간 변이 부분에 대한 정보를 제공해 주지 못하며, 이를 설명하기 위한 공변량도 사용자가 자유롭게 제어구문에 포함시키는데 제약이 있다. 이런 단점 때문에 현재까지 마취통증의학과 영역에서 다룬 거의 대부분의 반응표면모형은 NONMEM을 이용하여 분석한 결과를 제시하고 있다.

반응표면모형 종류

Minto model

Greco model이 일반적으로 많이 사용하고 있는 모형이기는 하지만 기본적으로 Greco model의 최대의 약점은 상호작용 약제들의 최대효과와 통상 Hill 계수라고 일컫는 가파름의 정도인 γ가 동일한 경우에만 적용할 수 있는 모형이라는 것이다. 이런 단점을 극복하기 위하여 Minto 등이[2] 새로운 형태의 약물상호작용 모형을 제시하였다. 먼저 각 약물의 농도를 각각의 역가(potency, C₅₀)로 나누어, 단위가 없는 표준화 농도(normalized concentration)로 변환한다.

식(7)UX= [X]C50X,     UY= [Y]C50  Y​

두 약물의 표준화 농도의 비를 θ로 표기하면 식(8)과 같다

식(8)θ= UYUX+UY

또한 두 약물의 합을 신약으로 보고 이를 U(θ)라고 하고, 역가를 U₅₀(θ)로 두자. 약물 X만 있을 경우 θ = 0, 약물 Y만 있을 경우 θ = 1이다. 즉, θ의 범위는 0과 1사이이다(0 ≤ θ ≤ ¹⁾. 두 약이 함께 투여될 경우, 약물 농도는 간단하게 U_X + U_Y이므로, 특정 비(ratio) θ에 대하여 농도-반응 관계는 식(9)로 표시가 가능하다.

식(9)Effect=E0+(Emax(θ)−E0)×(UX+UYU50(θ))γ(θ)1+(UX+UYU50(θ))γ(θ)

여기에서, 특정 비 θ에서, γ(θ)는 농도-반응 곡선의 가파름(steepness)의 정도를 의미하고, U₅₀(θ)는 최대효과의 반을 나타낼 때의 단위수(number of units, U)로 역가를 뜻하며, E_max(θ)는 최대효과를 의미한다. 두 약제의 비(ratio), 즉, θ마다 E_max, C₅₀, γ가 다르며, 이는 θ마다 그 자신의 S자형 농도-반응관계(sigmoidal concentration-response relation)가 있음을 의미한다. 약물 X만 C_50X의 농도로 존재할 경우, U_X = 1, U_Y = 0, θ = 0, 약물 농도(U_X + U_Y)는 1이된다. 따라서 U₅₀(0)도 1이다. 약물 Y만 C_50Y의 농도로 존재할 경우, U_X = 0, U_Y = 1, θ = 1, 약물 농도(U_X + U_Y)는 1이다. 따라서 U₅₀(1)도 1이다.

그러므로 X든, Y든 한 약물만 C₅₀의 농도로 존재할 경우 다음과 같이 일반화할 수 있다.

식(10)U(0,1)  =  UX+UY=1

이 때의 효과는 다음과 같다.

식(11)E=E0+(Emax(0,1)−E0)×(1U50(0,1))γ(0,1)1+(1U50(0,1))γ(0,1)

한 약물만 C₅₀의 농도로 존재할 경우 최대효과의 반에 해당하는 효과를 나타내게 된다.

식(12)E=E0+(Emax(0,1)−E0)×12

따라서 식(11)과 식(12)에 의하여 식(13)이 성립한다.

식(13)(1U50(0,1))γ(0,1)1+(1U50(0,1))γ(0,1)=12

따라서, U₅₀(0,1)은 1이 된다.

이번에는 각 약물이 역가 농도의 절반으로 함께 존재할 경우를 가정해 보자. 즉, U_X = 0.5, U_Y = 0.5, θ = 0.5, 약물 농도(U_X + U_Y)는 1이다.

식(14)E=E0+(Emax(0.5)−E0)×(1U50(0.5))γ(0.5)1+(1U50(0.5))γ(0.5)

이 경우, 약효가 최대효과의 반에 해당하는 효과보다 크다면, 상승적(synergistic)이고, 작다면 상가이하적(infra-additive)이며, 같을 경우 상가적(additive)이다. 즉, U₅₀(0.5) < 1일 경우 상승작용, U₅₀(0.5) > 1일 경우 상가이하작용, U₅₀(0.5) = 1일 경우 상가작용을 보인다.

상승작용을 보일 경우에는 최대효과의 반에 해당하는 효과보다 커야 하므로,

식(15)(1U50(0.5))γ(0.5)1+(1U50(0.5))γ(0.5)>12

U₅₀(0.5) < 1이 된다.

상가이하작용을 보일 경우에는 최대효과의 반에 해당하는 효과보다 작아야 하므로,

식(16)(1U50(0.5))γ(0.5)1+(1U50(0.5))γ(0.5)<12

U₅₀(0.5) > 1이 된다.

상가작용을 보일 경우에는 최대효과의 반에 해당하는 효과보다 작아야 하므로,

식(17)(1U50(0.5))γ(0.5)1+(1U50(0.5))γ(0.5)=12

U₅₀(0.5) =1이 된다.

이상을 종합하면, U₅₀(θ)는 θ = 0, 1일 경우(개별 약물이 단독으로 존재할 때) 1이며, 0 < θ < 1에서 U₅₀(θ) < 1일 경우 상승작용, U₅₀(θ) > 1일 경우 상가이하작용, U₅₀(θ) = 1일 경우 상가작용을 나타낸다.

위 사항을 일반화하면 E_max(θ), U₅₀(θ), γ(θ) 함수에 대한 정보가 많지 않기 때문에 실제 상황을 경험적(empirical)으로 잘 근사할 수 있도록 4차다항식(fourth-order polynomials)을 사용해 볼 수 있다.

식(18)f(θ)=β0+β1 ·  θ +β2 ·  θ2+β3 ·  θ3 +β1 ·  θ4

약물 X만 존재할 경우(θ = 0), E_max(0) = E_maxX, U₅₀(0) = U_50X, γ(0) = γ_X이다. 또한 f(0) = β₀이다. 약물 Y만 존재할 경우(θ = 1), E_max(1) = E_maxY, U₅₀(1) = U_50Y, γ(1)= γ_Y이다. 또한 f(1) = β₀ + β₁ + β₂ + β₃ + β₄이다. 따라서, β₁ = f(1) - β₀ - β₂ - β₃ - β₄이다.

이를 이용하여 E_max(θ), γ(θ), U₅₀(θ)에 대하여 정리하면 다음과 같다.

식(19)Emax(1)=EmaxY=β0,Emax+β1,Emax+β2,Emax+β3,Emax+β4,Emax∴β1,Emax =EmaxY−β0,Emax−β2,Emax−β3,Emax−β4,EmaxEmax(θ)=β0,Emax+β1,Emax×θ+β2,Emax×θ2+β3,Emax×θ3+β4,Emax×θ4β0,Emax =Emax(0)=EmaxX∴β1,Emax =EmaxY−EmaxX−β2,Emax−β3,Emax−β4,Emax∴Emax(θ)=(EmaxX+(EmaxY−EmaxX−β2,Emax−β3,Emax−β4,Emax)×θ+β2,Emax×θ2+β3,Emax×θ3+β4,Emax×θ4)

식(19)에서 이차항까지만 남기면 E_max(θ)는 식(20)으로 정리할 수 있다.

식(20)Emax(θ)=EmaxX+(EmaxY−EmaxX−β2,Emax)×θ+β2,Emax×θ2

식(21)γ(1)=γY=β0,γ+β1,γ+β2,γ+β3,γ+β4,γ∴β1,γ=γY−β0,γ−β2,γ−β3,γ−β4,γγ(θ)=β0,γ+β1,γ×θ+β2,γ×θ2+β3,γ×θ3+β4,γ×θ4β0,γ=γ(0)=γX∴γ(θ)=γX+(γY−γX−β2,γ−β3,γ−β4,γ)×θ+β2,γ×θ2+β3,γ×θ3+β4,γ×θ4

식(21)에서 이차항까지만 남기면 γ(θ)는 식(22)로 정리할 수 있다.

식(22)γ(θ)=γX+(γY−γX−β2,γ)×θ+β2,γ×θ2

식(23)U50(1)=U50  Y=β0,U50+β1,U50+β2,U50+β3,U50+β4,U50∴β1,U50=U50  Y−β0,U50−β2,U50−β3,U50−β4,U50U50  Y=U50(1)=1β0,U50=U50X=U50(0)=1∴β1,U50=−(β2,U50−β3,U50−β4,U50)U50(θ)=β0,U50+β1,U50×θ+β2,U50×θ2+β3,U50×θ3+β4,U50×θ4∴U50(θ)=1−(β2,U50×β3,U50×β4,U50)×θ+β2,U50×θ2+β3,U50×θ3+β4,U50×θ4

식(23)에서 이차항까지만 남기면 U₅₀ (θ)는 식(24)로 정리할 수 있다.

식(24)U50(θ)=1−β2,U50×θ+β2,U50×θ2

β_2,U50 이 0일 경우, 모든 θ에 대하여 U₅₀(θ) = 1이며, 이는 두 약물 간의 상호 작용이 상가작용을 의미한다. β_2,U50 > 0일 경우, 0 < θ < 1에서 U₅₀ (θ) < 1이므로 상승작용을 보인다. β_2,U50 < 0일 경우, 0 < θ < 1에서 U₅₀ (θ) > 1이므로 상가이하작용이다. 이는 두 약물의 최대효과가 동일함을 전제로 한다. 두 약물의 E_max와 γ가 동일하면서, 상승작용을 보일 때의 E_max(θ), γ(θ), U₅₀(θ)를 Fig. 6에 제시하였다. 두 약물의 E_max와 γ가 동일하다고 가정할 경우, NONMEM 제어구문은 다음과 같다.

Fig. 6

The functions for the parameters (Emax(θ), U50(θ), and γ(θ)) for the two drug-interaction showing synergic interaction.

$INPUT ID CONX CONY RESP=DV SEX; This example has SEX as a covariate

$PRED

C50X = THETA(1)*EXP(ETA(1))

C50Y = THETA(2)*EXP(ETA(2))

A = THETA(3)*EXP(ETA(3)); Interaction term

GAM = THETA(4)*EXP(ETA(4))

; Apply covariate effect

IF (SEX.EQ.1) THEN

C50A = C50A*THETA(5)

ENDIF

;Equations of Minto model

; Normalize to potency

UX = CONX/C50X

UY = CONY/C50Y

; Compute “theta” (Q) for Minto interaction model

IF(UX.EQ.0.AND.UY.EQ.0) THEN

Q = 0

ELSE

Q = UX/(UX + UY)

ENDIF

; Compute C50 Function In Combined Domain

U50 = 1 – A*Q + A*Q**2

; Compute “concentration” normalized to potency and corrected for interaction

U = (UX + UY)/U50

; Apply steepness term

UN = U**GAM

; Compute probability

P = UN/(1 + UN)

; Compute probability of observation

Y = RESP*P + (1 – RESP) * (1 – P)

$EST LIKELIHOOD METHOD=1 LAPLACE NOABORT

그렇지만 두 약물의 최대효과가 같지 않아도 상호작용 분석은 가능하다. 식(20), (22), (24)를 NONMEM code에 반영하여, E_maxY, E_maxX, β_2,Emax, γ_Y, γ_X, β_2,γ를 각각 추정하면 된다.

반응표면모형 구축 과정에서 고려할 사항

반응표면모형을 구축할 때는 크게 3가지 접근 방법이 있다. 첫 번째로 슬라이스(slices) 방법이 있는데, 한 약제는 다양한 가성평형상태(pseudo-steady state)의 농도로 유지시키고, 다른 약제는 최대 효과의 97% 정도에 이르도록 16분 동안 영차지속정주(zero-order infusion)를 하는 방법이다. 두 번째는 열십자(crisscross) 방법으로, 연구 참여자의 반은 약제 X의 농도를 일정하게 유지한 상태에서 약제 Y의 농도를 다양하게 변화시키고, 나머지 반은 약제 Y의 농도를 일정하게 유지한 상태에서 약제 X의 농도를 다양하게 변화시킨다. 세 번째는 방사상(radial) 방법으로, 두 약제 모두를 영차지속정주한다. 그리고 각 약제를 단독 정주하고, 표면의 세 가지 방사상 슬라이스(radial slices)를 다음과 같이 취한다. 1) 개별 약제의 집단(population) C₅₀의 함수로서 동일한 정주속도를 사용한다. 2) 이 정주속도를 유지하다가 약제 X의 정주속도는 50% 감소시키고, 약제 Y의 정주속도는 50% 증가시킨다. 3) 약제 X의 정주속도는 50% 증가시키고, 약제 Y의 정주속도는 50% 감소시킨다. 이 세가지 방법 중에서 열십자 설계(crisscross design)가 가장 강건한(robust) 방법임이 밝혀졌고[12], 실제로 여러 문헌에서 이 방법으로 반응표면모형을 구축하였다[11,13]. 반응표면모형에 관한 연구는 주로 건강한 자원자를 대상으로 진행되었는데, 환자를 대상으로 하기 어려운 점이 있기 때문이다. 우선 NONMEM으로 분석할 때 개인간의 변이를 고려하려면 한 연구 참여자에게 최소 3개 이상의 농도값이 필요한데, 각 농도 값에서의 평형에 도달하는 시간이 필요하기 때문에 환자를 대상으로 진행하기에는 어려움이 있다. 두 번째로는 통상 수면제와 아편유사제의 상호작용을 탐색하는 연구들이 많은데, 수면제 없이 아편유사제 단독으로 정주해야 하는 경우에 근육강직과 같은 아편유사제의 이상반응이 나타날 수 있다. 이러한 시간적 제약과 여러 이상반응에 대한 위험을 감수하고 원하는 반응을 임상적인 용량범위에서 모두 관찰하려면 건강한 자원자를 대상으로 진행하는 것이 용이하다. 건강한 자원자를 대상으로 열십자 설계로 진행할 때에는 최소 20명 이상의 연구참여자를 모집하여야 반응표면모형을 적절하게 구축할 수 있는 것으로 알려져 있다[12].

반응표면모형의 임상적 적용

현재까지 문헌으로 소개된 반응표면모형들은 주로 약제의 평형시간(equilibration time)이 짧아서 혈장과 효과처 농도 간의 항정상태에 빨리 도달하는 약제들 위주로 연구가 되었고, sevoflurane과 remifentanil의 약력학적 상호작용[8,14-18], propofol과 remifentanil의 상호작용이 대부분을 차지하고 있다[7,11,13,19-23]. 물론 이 모형들이 실제로 어떻게 작용하는지를 다소 모호하다고 느낄 수도 있을 것이다. 우선, 동등역가비(equipotency ratio)로 약제를 대체하는데 이용할 수 있다. 사람에서 최소폐포농도(minimum alveolar concentration)를 감소시키는 것을 토대로 fentanyl, sufentanil, alfentanil, remifentanil의 역가비(potency ratio)를 계산하면 1 : 12 : 0.0625 : 1.2가 나왔다[24,25]. 이를 토대로 fentanyl의 용량에 1.2를 곱하면 remifentanil을 fentanyl로 대체하여 사용할 수 있다. 약동학 연구를 기반으로 한 다른 연구에 의하면 fentanyl, sufentanil, alfentanil, remifentanil의 역가비가 1 : 9 : 0.014 : 0.435로 계산되기도 하였다[26,27]. 또한, 각성 시간을 시뮬레이션(simulation)하는 데 반응표면모형의 결과가 이용되기도 한다. Ting 등은[18] 척추측만증(scoliosis) 수술을 받는 환자를 대상으로 수술 중 각성 검사(Stagnara wake-up test)를 시행할 때 필요한 적절한 각성 시간을 예상하기 위하여 반응표면모형을 이용하였다. Sevoflurane과 remifentanil의 반응표면모형을 이용하여 임상에서 사용하는 desflurane과 fentanyl 상호작용으로 확장하여 각성의 예측시간을 계산하였는데, 실제 환자가 각성한 시간과 반응표면모형을 이용하여 예측한 각성시간 간의 시간 차이가 평균 2.6분으로 상당히 정확하게 예측하였다[18]. 이 뿐만 아니라 반응표면모형을 이용하면 환자를 가장 빨리 각성시킬 수 있는 적절한 약제의 농도조합을 선택할 수 있는 장점도 있다.

결론

이제까지 이소볼과 반응표면모형(response surface model)의 개념과 반응표면모형의 종류 및 임상적 적용에 대하여 살펴보았다. 이소볼이 특정 효과(예: 50% 효과)를 나타내는 두 약제의 조합을 나타낸 동일효과곡선이라면 반응표면모형은 두 약제의 용량 혹은 농도 조합에서 나타낼 수 있는 반응의 효과를 전체 효과 범위에서 제시한 모형이다. 이를 잘 이해하고 있으면, 현재 임상에서 사용하는 여러 약제의 병용 투여 과정에서 약제의 세밀한 적정을 보다 효과적으로 할 수 있을 것이다.

1성별, 체중, 연령, 유전자 변이, 질환의 정도 등과 같은 환자 특이적인 공변량을 이용하여 약동, 약력 모수의 개인간 변이(inter-individual variability)를 설명하려는 분석방법. 이를 위하여 개인별 시간농도 곡선을 따로 분석하였던 standard two stage method와는 달리 모든 개체를 한꺼번에 분석함으로써 어떤 특성이 약동, 약력 모수와 연관되는지 탐색할 수 있다.

Article information Continued