This study was conducted to investigate the membrane potential to apply Doksam-tang to a human heart failure model. The human heart model was built by Luo et al. CellML model, Priebe et al. CellML model, and a human heart mesh file. Doksam-tang gives channel the half maximal inhibitory concentration(IC ₅₀ ), half maximal effective concentration(EC ₅₀ ) values and compounds concentrations. These data load into the laptop with Ubuntu OS, and build the library with the data. While results of the study with the heart failure model shows abnormal membrane potential from the normal heart model, the study with applying Doksam-tang to heart failure model shows restoring membrane potential that is similar to normal heart model. These results of the testings suggest that a conception of novel technique to investigate the effects of Korean herbal medicine. 생물정보학 분야의 지식이 쌓여가면서 생체에 관한 많은 정보가 데이터베이스화 되었다. 연구자들은 그동안 쌓인 데이터베이스에 알려진 처리방식을 결합해 생체시스템에 대한 라이브러리(Library) ^1-3) 를 만들었다. 라이브러리는 데이터의 처리방식과 관련된 파일들의 모임으로 정보에 접근하는 것을 용이하게 해주는데, CellML(Cell XML database) ²⁾ 과 CHASTE(Cancer, Heart And Soft Tissue Environment) ³⁾ 도 그 중 하나이다. CellML은 270여 가지의 전기생리학적 세포모델을 담고 있고, XML(extensible markup language)을 이용하여 만들어졌다. XML은 HTML을 개선하여 만들어진 프로그램 언어이며 HTML과는 달리 구조화된 데이터베이스를 사용자가 직접 조작할 수 있다. 그러므로 CellML의 데이터베이스는 몇 가지 과정을 거친 후 시스템상의 라이브러리에 적용할 수 있다. CHASTE는 다양한 형식의 데이터를 변환하여 CHASTE 라이브러리에 통합할 수 있는 기능과 미분방정식을 푸는 기능을 핵심기능(core functionality)으로 갖고 있다. 그 중 세포의 생리학적 기전을 풀 때 사용되는 ODE(Ordinary differential equation, 상미분방정식) solver와 PDE(Partial differential equation, 편미분방정식) solver는, 컴퓨터의 메모리를 직접 조작할 수 있는 C++로 코딩되어 있기 때문에 처리속도가 상용라이브러리보다 빠르다. 또한 공개운영체제인 리눅스에서 개발되어 사용자가 직접 확장기능을 추가할 수 있으며, 사용과 배포에 있어서 상용라이브러리와는 달리 저작권문제에서 자유로울 수 있다. 지금까지 시뮬레이션 도구들은 주로 컴퓨터 공학자들과 생물학자들이 만들었다. 그렇기 때문에 프로그램 자체의 오류는 없는지 ^4-9) 와 약리학적으로 알려진 방정식을 정확히 묘사하는지 ^10-15) 가 주된 연구내용이었다. 또한 대부분 상용프로그램이었기에 그 적용에 한계가 있었다. 그러나 CHASTE나 CellML 등의 공개된 소스를 이용하면, 생물정보학 분야의 지식을 약물과 인체의 복잡한 상호작용에 관한 이해 ¹⁶⁾ 를 심화시키는데 적용할 수 있을 것이다. 본 연구에서는 정상적인 심장세포를 묘사하는 Luo et al. CellML 모델 ¹⁷⁾ 과 심부전심장세포를 묘사하는 Priebe et al. CellML 모델 ¹⁸⁾ 그리고 심장세포의 간소화된 좌표를 담고 있는 파일 ³⁾ 을 이용하여, 정상심장모델과 심부전심장모델을 만들고 그들 간의 차이를 확인하고자 하였다. 그 후 숨이 차고, 기운이 없고, 맥이 몹시 약한데 처방하는 獨蔘湯 ¹⁹⁾ 의 K ⁺ 채널, Na ⁺ 채널, Ca ²⁺ 채널 등의 IC ₅₀ , EC ₅₀ 값을 변수로 설정했을 때 나타나는 심부전과 관련된 심장의 막전위 개선효과를 예측하고자 하였다. 이러한 결과를 통해서 한약의 심장 전기역학적 생리반응성 평가에 새로운 방법론을 제시하였다. 운영체제는 Linux OS Ubuntu 64bit 12.04 LTS와 Linux OS Ubuntu 64bit 14.04 LTS, 컴파일러는 gcc 4.8.1(GNU Compiler Collection), 빌드는 SCons v2.1.0.r5357, 편집기는 vim editor와 그 옵션인 ctags, cscope 등을 사용하였다. SVN 등을 이용한 프로젝트 버전관리에는 관여하지 않았다. Oxford Univ.의 Chaste executable 3.2 64bit version, Chaste-source code 3.2 64bit version을 사용하였다. 심장세포는 Luo et al. ¹⁷⁾ , Priebe et al. ¹⁸⁾ 의 CellML모델을 사용하였고, 심장의 기하학적 정보를 담고 있는 격자(mesh)파일은 인체심장의 좌표 ³⁾ 를 가공하여 사용하였다. 수치해석을 통해 계산된 값들은 모든 세포의 위치에서 저장되는데, 시각화도구를 사용하면 전문적인 지식이 없이도 생리적인 결과를 쉽게 이해할 수 있도록 Gnuplot 4.4 version ²⁰⁾ , Paraview 3.14.1 version ²¹⁾ , meshylyzer ³⁾ 를 사용하였다. 빌드는 프로그램 소스를 사용할 수 있도록 개인의 컴퓨터 환경에 맞게 적절히 적용하는 작업이다. 본 연구에는 작성된 프로그래밍 코드로 Chaste Heart Executable을 빌드하고, 빌드된 Executable 파일로 시뮬레이션을 진행하였다. ode, pde, printing을 각각 0.005 ms, 0.01 ms, 0.1 ms로 하였다. 계산방식은 bi-domain 방식으로 하였다. bi-domain은 전류가 한 면으로 흐른다고 가정하는 mono-domain과는 달리, 전류가 설정한 격자의 안쪽과 바깥쪽 면을 따라 흐른다고 가정한다. 전기역학법칙에 따르면 bi-domain 방식이 mono-domain 방식보다 실제상황을 더 잘 설명하며, 설정은 소스 내에서 사용자가 직접 해주어야 한다 ²²⁾ . Mono-domain을 푸는 방정식은 수식 (A)와 (B)로 이루어져 있다 ²²⁾ . (A)에서 χ 는 표면적과 체적의 비율인데 공통적으로 1400(1/cm)이라 설정되어 있다. I_ion 은 단위면적당 이온의 흐름을 나타내고 그 값은 특정한 벡터 와 전위 V 에 대한 함수이다. I ^(vol) 은 단위부피당 stimulus current를 나타낸다. 효과전도도 σ 값은, 특정상수 λ가 있을 때, 이다. mono-domain의 (A), (B)에 (D)가 추가된 (C), (D), (E)가 bi-domain을 푸는 수식이다. 경계조건은 (F), (G), (H)로 주어진다 ²²⁾ . (F)는 (D)와 의미가 같다. (G)와 (H)는 안쪽 표면과 바깥쪽 표면을 흐르는 전류를 계산하는 수식이다. mono-domain과 bi-domain의 차이는 가 된다. 심장의 세포는 종류에 따라 이온채널에 관한 많은 수식을 갖고 있고, bi-domain 수식을 통하여 세포끼리 상호 연결된다. 계산은 모든 세포에서 동시에 진행되는데, 만일 세포 중에 하나라도 정상적인 심장의 막전위를 벗어나는 값이 계산되거나, 이온채널의 정상적인 작동범위를 벗어나는 값이 계산되면 라이브러리 및 시스템에서 오류가 보고된다. 심장의 막전위는 ECG를 이해하는데 기본적인 요소가 되므로 ¹⁵⁾ 본 연구에서는 심장의 막전위(mV)를 관찰하였다. 심장의 세포 내, 외의 전기전도도는 모델마다 세포의 상호작용을 보기 위해 2.5(mS/cm)를 넘지 않는 값으로 설정하였다. 500 (ms), -80000.0 (uA/cm^3), 1400 (1/cm), 1.0 (uF/cm^2)을 상수로 하였고, DrugAction은 Ca ²⁺ 채널들, Na ⁺ 채널들, K ⁺ 채널들에 IC ₅₀ 값을 변수로 설정한 후 결과를 관찰하였다. LuoRudy et al. CellML모델과 사람의 심장격자 파일을 사용하였는데, 심장의 기하학적 크기는 약 2 cm이며 심장의 Apex가 원점이다 ^3,21) . 전기자극은 동방결절에서 시작되는데, 본 논문에선 심장세포들의 기하학적 좌표를 고려하여 심장의 Apex(Apex의 위치는 {(x, y, z) (cm)| (x, y)=(0, 0), z=0.04248})까지 도달한 이후의 전위를 관찰하였다. 세포 내와 세포 외의 longi, trans , normal 방향의 전기전도도 모두 0.2 mS/cm로 설정하였다. stimulus는 심장의 Apex 위치에서 2 ms간 – 50000.0 (uA/cm^3)으로 발생시켰고, 세 위치(#1, #2, #3)에서 측정한 막전위를 대표로 도시하였다( Fig. 1 ). Priebe et al. CellML 모델은 심부전심장의 세포모델인데, 정상심장 설정의 좌심실을 지나는 부위에 본 모델을 배치시켰다. 세 위치에서 측정한 막전위를(#1, #2, #3)대표로 도시하였고 정상심장 설정의 모델과 stimulus는 동일하게 적용했다.( Fig. 2 .) {(x, y, z) (cm) | 1.2 < x < 1.4, -100 < y < 100, -100 < z <100}}이 Priebe et al. CellML이 배치된 곳인데, 이 중 원래 심장의 좌표를 벗어나서 설정되는 범위는 프로그램 내에서 자동으로 무시된다(직육면체로 배열 된 Priebe et al. 세포 중 심장의 좌표에 들어가는 것만 시뮬레이션 시 적용된다). Na ⁺ 채널, Ca ²⁺ 채널, K ⁺ 채널에 人蔘의 성분이 복합적으로 작용하므로, 각 채널의 IC ₅₀ 값을 Cisapride와 獨蔘湯의 약물농도를 참조해 설정하였다. Na ⁺ , Ca ⁺ , K ⁺ 의 IC ₅₀ 값은 각각 0, 0, 0.0445 μM 로 하였다. 약물농도는 5 mM 로 하였다 ²³⁾ . 獨蔘湯은 人蔘을 40∼80 g으로 달이며 위급한 상황에서 투여하는 처방이다. 한 첩은 150 ml로 설정했다. 獨蔘湯의 투여는 1)추출, 2)장내 성분별 흡수 및 대사, 3)장기 내 성분별 농도, 4)성분별 작용의 순서로 고찰하였다. Ginsenoside를 비롯한 다양한 물질이 함유되어 있는 人蔘 각 성분의 심장이온채널에 대한 IC ₅₀ , EC ₅₀ 값과 농도를 최종적인 변수로 사용하였다. 人蔘은 추출용매, 온도, 시간에 따라 추출물의 함량이 달라지나, 100 ℃의 물에서 1시간 동안 달일 경우 22.86%의 엑기스가 추출되고 그 중 5.22%는 ginsenoside이다. 따라서 인삼 80 g으로 끓는 물에서 달이는 獨蔘湯의 경우, 22.86%는 엑기스로 추출되는데 그 용량은 18.3 g이다. 그 중 5.22%에 해당하는 0.96 g이 獨蔘湯 안에 들어있는 ginsenoside 용량이다. 전신순환에 도달하는 투여약물의 분획인 생체 이용률을 구하였다. 사람에게 ginsenoside를 경구투여 할 때, 위와 장에서 흡수 및 대사가 이루어진다. 이는 위(Stomach)에서만 흡수 및 대사가 이루어지는 경우, 장에서만 흡수 및 대사가 이루어지는 경우, 위에서 대사된 것이 장에서 다시 대사를 거치고 흡수되는 경우, 장에서 대사된 것이 장에서 한 번의 대사를 더 거치고 흡수되는 경우의 총 4가지 경로로 나뉜다. 각 경로마다 성분의 흡수시간이 상이하다. 그 중 0~3 시간 내에 흡수가 이루어졌다고 추정할 수 있는 물질로는 Rh ₁ , Rg ₁ , Rd, Re, Rb ₂ , Rc 등이 있는데, PPT(protopanaxatriol)는 최대 1.2%, PPD(protopanaxadiol)은 최대 0.2%가 urine 내에서 검출된다. ginsenoside의 경우 일반적으로 urine을 통해 얻은 데이터와 plasma를 통해 얻은 데이터는 일치한다 ²⁷⁾ . Urine 내의 검출량을 plasma 내 잔존량이라 할 때 심장에 작용하는 약물의 농도는 다음과 같이 구하였다( Table 1 ). 엑기스의 여러 성분이 다양한 작용을 일으키는데 ²⁸⁾ , 각 추출방법에 따른 人蔘 추출물의 성분 값들은 알려져 있다 ^24,26,28) . Ginsenoside는 38가지 이상인데 ^27,29,30) , 심장의 이온채널에 작용하는 ginsenoside의 성분에 따른 IC ₅₀ 값들은 서로 상이하다 ^23,31-38) . 그러므로 심장의 이온채널에서 ginsenoside의 각 성분에 대한 IC ₅₀ , EC ₅₀ 값을 장기 내 성분별 농도에 가중치를 주어 평균을 계산한 후 최종채널에 적용되는 값을 정하였다( Table 2 ). #1, #2, #3 세 위치( Fig. 1. )에서의 정상심장에 관한 모델은 정상적인 심실세포의 전위와 같은 모습을 보여주었다( Fig. 3. , Fig 4. ). #1, #2, #3 세 위치( Fig. 2. )에서의 비정상심장에 관한 모델은 정상적인 심실세포의 전위와는 다른 양상을 보여주었다( Fig. 5. ). Priebe et al. CellML 모델을 사용한 비정상 심장은 약물검사를 위해 Cisapride의 변수를 설정했을 때 심장의 막전위 수치는 비정상적인 변화를 보여주었다( Fig. 6. , Fig. 7. ). 비정상심장의 약물반응모델에 獨蔘湯에 해당하는 약물변수를 넣으면(#2) 정상적인 심실의 전위패턴( Fig. 3. )을 회복하는 변화를 보여주었다( Fig. 8. ). 獨蔘湯은 人蔘 한 가지만을 고농도로 달여 먹는 처방으로 氣脫, 氣虛血脫 등에 쓰인다 ¹⁹⁾ . 人蔘은 여러 가지 성분과 효능이 보고된 약물이나, 다양한 성분이 동시에 작용할 때 심장의 전위에 어떤 효과가 발생하는지에 대해서는 아직도 미지수이다. 따라서 본 논문에서는 심근세포에 다양한 채널에 작용하는 ginsenoside를 복합적으로 적용시켜서 심부전을 일으킬 수 있는 심장의 막전위가 어떻게 변화하는지 시뮬레이션을 통해 확인하였다. 변수화한 獨蔘湯에 관한 심장의 반응을 시뮬레이션을 통해 살펴 본 결과, 비정상적인 패턴을 보이던 심실의 전위가 정상적인 패턴으로 돌아오는 것을 확인하였다. 본 연구는 이미 알려져 있는 표준화된 세포의 모델을 통해 진행하므로 다중회기분석 등의 통계적인 방법을 사용할 필요가 없다. 또한 실험 세포 등의 조건은 이미 CellML에서 연구된 자료를 기초로 하므로, 다양한 변이가 있는 생명체에 실험을 진행하는 것과는 달리 실험대상의 조건은 충분히 서술되어 있다(그러나 표준화된 모델이 획일적인 모델을 뜻하는 것은 아니다. 알려진 자료로 조합하여 만들 수 있는 생체모델의 경우의 수는 무한히 많다.). 본 논문에서는 랩탑 계산능력의 한계로 인해 다소 단순한 실험을 진행 하였지만, 클러스팅 컴퓨터로 진행되는 연구는 인체 내 조직을 실제와 유사하게 설정하여 어떠한 변화가 나타나는지를 정밀하게 다룰 수 있다. 그러면 채널에 적용할 수 있는 모든 요소를 고려하면서 더 다양한 세포들의 군집을 실제 심장모델에 적용할 수 있으므로, 정확한 시뮬레이션이 될 수 있다. 다만 본 논문에서는 대사정보가 충분한 인삼을 대상으로 시뮬레이션을 진행하였는데, 대사정보가 충분하지 않은 다른 약물의 경우 본 연구기법을 바로 적용하기에는 무리가 따른다. 그러므로 앞으로는 효능 및 약물 대사정보에 대한 연구가 이루어져야한다. 이를 토대로 한약에 관한 구체적인 정보를 얻을 수 있을 것이고, 이 정보와 함께 생체 시뮬레이션 정보를 조합해서 새로운 양생방법과 치료방법을 얻을 수 있을 것이다. 본 연구는 소스코드를 직접 이용하는 방식을 택했고, SVN 등을 이용한 프로젝트의 버전관리에는 관여하지 않았다. 앞으로의 연구는 프로젝트의 버전관리에 참여하면서 클러스팅 컴퓨터를 이용해 보다 많은 변수를 관찰하는 실험이 수행되어, 시뮬레이션 상의 생체와 약물이 복합적으로 작용하였을 때 어떤 반응이 일어나는지 살펴보아야 할 것이다. 정상심장과 비정상심장의 막전위는 서로 상이했다. 비정상심장의 막전위는 정상적으로 동작하지 않았는데, 정상세포에서 측정한 막전위까지 변화를 보였다. Drug-Action을 적용하기 위해 비정상심장모델의 특정한 채널을 변화시켰을 때 심부전이 일어나는 것을 관찰할 수 있었는데, 獨蔘湯의 알려진 성분(Na ⁺ , Ca ²⁺ , K ⁺ IC ₅₀ 값)들을 모델에 투여했을 때 비정상심장의 막전위는 정상적인 심장의 막전위 패턴을 보여주었다. 이러한 연구는 특정질환의 모델과 그 모델에 유효하며 안전할 것으로 예측되는 약물을 시뮬레이션을 통해 사전에 찾아낼 수 있는 유익한 방법이 될 것이다.