1. 서 론
칡(kudge)은 다년생 식물로서 굵은 줄기로 이루어져 나무로 분류된 다. 이는 오래전부터 자양강장제 등 건강식품으로 이용되었으며 특히 뿌리부분은 항산화 물질이 다량 함유되어있어 최근 많은 연구가 이뤄 지고 있다[1-2]. 칡에는 flavonoid 계통 화합물의 한 종류인 iso-flavone 이 다량 함유되어있으며 이는 여성호르몬인 에스트로겐과 비슷한 식 물성 에스트로겐(phytoestrogen)이다. Iso-flavone은 일반 세포 내에서 estrogen receptor와 결합할 수 있는 능력이 있으며 체내에서 심장병, 골다공증의 발생 빈도를 낮춰주며 여러 갱년기 증상을 완화시켜 주는 것으로 알려졌다[3-4]. 칡에는 다량의 전분 외에 iso-flavone계 성분인 puerarin, daidzein, daidzin, genistein, genistin, formonetin, coumesterol, 4’-7-dihydroxyisoflavone, 7-xyloside puerain, 3’-hydroxypuerain, 3’-methoxypuerin, β-sitostrol, lupeol, lupenone, allantoin 등이 함유되 어 있다[5-6].
본 연구에서는 칡에 함유된 iso-flavonoids인 puerarin, daidzein, daidzin의 추출수율의 최대화를 위한 공정조건의 최적화를 진행하였 다. 최근 천연물로부터 유효성분을 추출하는 공정방법은 열역학적 특 성을 고려한 열수추출, cavitation을 이용한 초음파추출, 마이크로웨이 브 추출, 다양한 용매 추출 등 매우 다양하다[7-9]. 본 연구는 칡으로부 터 항산화 물질을 추출하는 방법으로 초음파 용매추출공정(ultrasonicassisted solvent extraction process)을 이용하였다. 이 공정은 기존의 열수추출방법에 비해 추출시간이 빠르고 효율이 높다고 알려져 있다 [10-12]. 칡으로부터 이소플라보노이드를 추출하는 용매 중 주정 (ethanol)을 선정할 경우 산화스트레스를 줄일 수 있다는 연구결과가 존재하며[13-14], 칡뿌리로부터 puerarin, daidzein, daidzin을 추출하기 위해 에탄올을 용매로 사용한 연구가 진행되고 있다[15]. 또한 추출물 에 함유된 iso-flavonoids (puerarin, daidzein, daidzin)의 분석은 고성능 액체크로마토그래피를 이용하여 수행하였다[16-17].
추출공정의 최적화를 위해 반응표면분석법(response surface methodology, RSM)을 이용하였다. 반응표면분석법은 통계학적 분석방법 으로 기초실험을 통해 변수에 따른 반응결과 값을 예측하여 원하는 반응치에 따른 최적조건을 도출해주는 방법이다. 따라서 이는 복잡한 process들을 최적화하는데 매우 효과적이고 체계적인 기술이며 실험 에서의 각 independent variables에 따른 효과와 최적조건들을 빠른 시 간에 효율적으로 찾는데 많이 활용된다[18-19]. 본 연구에서는 기초실 험을 바탕으로 반응표면분석법을 활용하여 실험횟수를 줄이고 원하 는 결과 값을 도출하는 최적공정조건을 분석했다. 또한 각 공정조건 이 반응치에 미치는 주효과도와 교호효과도를 관찰하고 결정계수 (R2)와 P-value및 F-value를 해석하여 본 실험에서 사용한 분석방법의 적합성을 판단했다[20-21]. 독립변수로 초음파 조사시간, 주정/초순수 부피비, 초음파 조사세기를 설정하였고, 기초 실험 data를 바탕으로 추출효율과 iso-flavonoids 함량을 반응치로 설정하여 최적화과정을 수행하였다.
2. 실험방법
2.1. 초음파 용매추출
건조된 칡 12.5 g을 추출용매 250 mL에 넣고 초음파 용매추출과정 을 진행하였다. 이때 추출용매는 주정(Daehan Ethanol life, 95%)과 초 순수(Vivagenen EXL-III)를 부피비로 혼합하여 사용하였다. 초음파 용매추출을 위해 초음파 발생장치는 Sonics & Material Inc.사의 ultrasonic generator (VCX-600, U.S.A)와 probe type (diameter: 0.5 in, threaded end, amplitude: 124 nm)의 horn을 사용하였다. 또한 반응표 면분석법을 이용한 최적화 과정은 기초실험을 통해 독립변수를 초음 파 조사시간(20.0~30.0 min), 주정/초순수 부피비(35.0~45.0 vol%), 초음파 조사세기(550.0~650.0 W)로 설정하였다. 추출이 완료된 추출 물은 여과 후 진공회전증발기(EYELA N-N type, Aspirator A-3S)를 이용하여 농축하고, -10 ℃로 6 h 동안 동결한 후 냉동건조장치를 이 용하여 -40 ℃, 0.1 mbar의 조건에서 24 h 동안 건조하였다.
2.2. Iso-flavonoids 함량 분석
건조된 추출물 10 mg을 증류수에 녹여 10 mg/mL로 만든 후 10,000 rpm에서 5 min간 원심분리하여 상층액을 분리하였다. 이를 syringe filter (Nylon 66, 0.45 m, Whatman, Kent, UK)로 여과한 후 100 μL를 분취하여 내부표준물질인 geniposide (Sigma-Aldrich)를 20 μL와 혼합 하고 증류수를 넣어 총부피 1 mL로 제조하였다. 제조된 시료를 고성 능 액체크로마토그래피(YL 9100 HPLC System, Young Lin Instrument, Korea)를 이용하여 iso-flavonoids 함량 분석을 수행하였 다. 표준물질로는 puerarin, daidzein, daidzin (Sigma-Aldrich Co.)을 이 용하였으며, 농도범위는 2.5~50 μg/mL이다. 분석에 사용한 컬럼은 SunFireⓇ C18 (4.6 × 150 mm, Waters, USA)이고, 이동상으로는 acetonitrile과 acetic acid, 초순수를 사용하였다. 검출기의 자외선 파장은 250 nm로 설정하였으며, 시료주입량은 10 μL이다.
2.3. 반응표면분석법을 이용한 최적화
반응표면분석법을 이용한 최적화를 위해 Minitab 프로그램(Minitab Inc., USA)을 사용하였으며 실험계획, 수학적 모델링, 최적화 등을 수 행하였다. 반응표면분석법 중 중심합성계획모델(central composite design, CCD)은 3개의 독립변수와 3개의 수준(level), 20번의 실험과정 과 중심점(central point)에서의 6번 반복실험으로 구성된다. 3개의 실 험인자는 초음파 조사시간(x1), 주정/초순수 부피비(x2), 초음파 조사세 기(x3)이고, 이 인자들은 반응표면분석법에 의해 최적화된다. 반응치 는 추출수율(Y1)과 iso-flavonoids 함량(Y2)이며, 이는 아래의 식과 같이 2차 회귀방정식으로 표현된다.(1)
여기서 a0는 상수이며 ai, aii, ajj는 각각 1차, 2차, 상호작용에 따른 계수이다. xi와 xj는 독립변수인 초음파 조사시간, 주정/초순수 부피비, 초음파 조사세기이고, 반응치인 Y1과 Y2를 동시에 만족하는 최적조건 을 도출하였다. 또한 다변량 해석(multivariate analysis)은 MATLAB 프로그램(The MathWorks, MA, USA)을 이용하였다.
3. 결과 및 고찰
3.1. 기초실험 및 반응표면분석법 적용
초음파에너지를 이용하여 칡으로부터 유효성분 추출과정의 최적화 를 위해 독립변수에 따른 추출효율과 iso-flavonoids 함량을 분석하였 다. 실제로 반응표면분석법은 실험횟수를 줄이고 원하는 반응치에 따 른 조건을 최적화하기 위해 여러 연구에서 활용된다[22-24]. 본 연구 에서는 반응표면분석법을 사용하기 위하여 초음파 조사시간(10.00~ 40.00 min), 주정/초순수 부피비(30.00~50.00 vol%), 초음파 조사세기 (500.00~700.00 W)의 범위에서 기초실험을 진행하였다. 기초실험 결 과 높은 반응치를 갖는 각 변수의 적정범위를 초음파 조사시간 (20.00~30.00 min), 주정/초순수 부피비(35.00~45.00 vol%), 초음파 조사세기(550.00~650.00 W)로 좁게 설정한 후 Minitab 프로그램을 활용하여 반응표면분석을 진행하였다. 3개의 독립변수는 초음파 조사 시간(x1), 주정/초순수 부피비(x2), 초음파 조사세기(x3)로 설정하고 20 번의 표면설계를 생성했고 반응치인 추출수율(Y1)과 iso-flavonoids 함 량(Y2)을 반응표면분석법에 의해 최적화했다. 중심합성계획모델에 의 해 설정된 계량인자인 초음파 조사시간, 주정/초순수 부피비, 초음파 조사세기에 따른 반응치인 추출수율과 iso-flavonoids 함량을 Table 1 에 나타내었다.
Table 1의 실험결과를 바탕으로 각 독립변수에 따른 추출수율을 최 적화할 때 결정계수(coefficient of determination, R2), P-value, F-value 등을 Table 2에 나타내었다. 여기서 결정계수는 예측값과 실제값 사이 의 관계를 나타내는 지표로써 추출수율을 해석할 경우 R2 = 0.94로 높 게 나타나 실험값과 예측값의 상호관계가 높은 것으로 해석된다. P-value란 실험값의 검정 통계량이 귀무가설을 지지하는 정도를 확률 로 표현한 것이다. 따라서 P-value가 작을수록 귀무가설을 지지하는 정도가 약하므로 귀무가설을 기각하게 된다. P-value (유의확률 값)가 유의수준보다 작으면 귀무가설 하에 일반적인 형상이 아니므로 귀무 가설을 기각하고 대립가설을 채택하게 되며 일반적인 유의수준으로 는 5%를 설정한다. 또한 F-value란 분산 분석의 모형에 대한 적합성 을 나타내며 값이 클수록 적합하다고 볼 수 있다[25-26].
3.2. 추출수율의 최적화과정
추출수율에 대한 F-value는 16.81로 뛰어난 적합성을 나타내었으며, P-value 또한 매우 작게 나타났다. P-value는 각각의 독립변수의 중요 도를 확인할 수 있는 하나의 도구로써 그 값이 작을수록 더 중요한 독립변수로 판단할 수 있다. Table 2로부터 추출수율에 따른 각 변수 의 P-value는 x1 = 0.187, x2 = 0.587, x3 = 1.76으로 x1의 P-value가 가 장 작다. 따라서 본 연구에서의 추출수율의 경우 초음파 조사시간(x1) 이 가장 큰 영향을 미치는 변수임을 알 수 있었다. 이를 바탕으로 반 응표면분석법을 이용한 각 독립변수에 따른 추출수율의 2차 회귀방정 식을 다음과 같이 산출하였다.(2)
여기서 Y1은 추출수율이고, x1, x2, x3는 각각 독립변수인 초음파 조 사시간, 주정/초순수 부피비, 초음파 조사세기이다.
반응표면분석법을 이용하여 얻은 2차 회귀방정식을 이용하여 3가 지 독립변수 중 2가지 독립변수에 따른 추출수율의 관계를 Figure 1에 나타내었다. 계량인자인 독립변수가 2개 이상일 경우 인자수준의 조 합에서 일어나는 효과를 나타내는 교호효과도를 해석할 수 있다. (초 음파 조사시간 × 초음파 조사세기)와 (주정/초순수 부피비 × 초음파 조사세기)에 비해 (초음파 조사시간 × 주정/초순수 부피비)가 추출수 율에 더 큰 영향을 미치는 것을 알 수 있었다. 이는 Table 2에 x1x2의 F-value (0.10)와 P-value (0.757)를 통해 x1x3와 x2x3에 비해 높은 F-value 와 낮은 P-value 나타내었다.
Figure 2는 각 변수들의 주효과도를 나타낸 그래프이다. 주효과도란 각 변수들이 반응치에 영향을 미치는 정도를 뜻하며 이 그래프에서의 반응치는 추출수율이고 각각의 변수는 초음파 조사시간, 주정/초순수 부피비, 초음파 조사세기이다. Figure 2에서 반응치인 추출수율에 초 음파 조사시간이 상대적으로 가장 큰 주효과를 나타내고 있다. Figure 1 에서 초음파 조사시간에 따른 추출수율은 일정구간 증가하다가 최고 점을 지나 감소하는 경향을 보인다. 이는 초음파 조사시간이 증가할 수록 용액의 온도가 증가하게 되고 이에 따라 추출효율이 증가하지만 일정 구간을 넘어가면 오히려 추출효율이 감소하는 것을 보여준다. 칡에 존재하는 puerarin, daidzein, daidzin은 glycoside iso-flavonoids로 초음파 조사시간이 과도하게 증가할 경우 degradation과 hydrolysation 을 초래한다는 연구결과가 존재한다[27-28]. 주정/초순수 부피비의 경 우 최고점까지 증가하는 추세를 보이다 감소하는 경향을 나타내었다. 이는 주정/초순수 부피비의 경우 고체와 용매 간 물질전달에 영향을 미치기 때문으로 사료된다[29]. 초음파 조사세기의 경우에도 추출수 율이 증가하다가 감소하는 경향을 보이는데, 이는 일정 수준까지 추 출 효율을 높이다가 한계점 이상에서 오히려 높은 에너지로 인해 acoustic cavitation을 유발시키고 cell 장벽의 붕괴를 초래하기 때문으 로 사료된다.
3.3. Iso-flavonoids 함량의 최적화과정
Table 2의 iso-flavonoids 함량의 경우 R2 = 0.75로 나타나 실험값과 예측값의 상호관계가 비교적 높았으며, F-value는 3.4, P-value가 0.035 로 5%로 이내의 유의수준을 나타내었다. 따라서 본 연구에서는 각 독 립변수에 따른 iso-flavonoids 함량의 관계를 다음과 같은 2차 회귀방 정식으로 나타내었다.(3)
여기서 Y2은 iso-flavonoids 함량이고, x1, x2, x3는 각각 독립변수인 초음파 조사시간, 주정/초순수 부피비, 초음파 조사세기이다.
Figure 3은 반응표면분석법을 이용하여 얻은 2차 회귀방정식을 이 용하여 독립변수에 따른 iso-flavonoids 함량을 나타낸 그림이다. iso-flavonoids 함량의 경우에는 (초음파 조사시간 × 주정/초순수 부피 비)와 (주정/초순수 부피비 × 초음파 조사세기)가 (초음파 조사시간 × 초음파 조사세기)에 비해 더 많은 영향을 미치는 것을 알 수 있었다. Table 2로부터 x1x2와 x2x3의 F-value와 P-value가 각각 0.00, 0.954와 0.00, 0.978로 x1x3의 0.02와 0.889보다 더 높은 F-value와 더 낮은 P-value를 나타내었다. Table 2로부터 iso-flavonoids 함량에 따른 각 변수의 P-value는 x1 = 0.035, x2 = 0.971, x3 = 0.809로 x1의 P-value가 가장 작다. 따라서 본 연구에서 iso-flavonoids 함량에 가장 큰 영향을 미치는 변수는 초음파 조사시간(x1)임을 알 수 있었다. 또한 반응치인 iso-flavonoids 함량에 따른 주효과도를 나타낸 Figure 4을 통해 초음 파 조사시간 > 주정/초순수 부피비 > 초음파 조사세기의 순서로 반응 치에 가장 큰 영향을 미친다는 것을 확인할 수 있다. 칡의 주요 iso-flavonoids인 puerarin, daidzein, daidzin의 추출함량을 이용해서 실 험을 진행한 결과 초음파 조사시간이 증가함에 따라 iso-flavonoids 함 량이 증가하다가 감소하는 경향을 보였다. 또한 초음파 조사시간이 증가할수록 용액의 온도가 증가하여 iso-flavonoids의 용해량이 증가 하게 된다. 하지만 일정 시간 이상부에서는 iso-flavonoids의 변형이나 분해 등이 발생하여 오히려 iso-flavonoids 함량이 낮아지게 된다. 주 정/초순수 부피비 또한 한계점 이하에서는 iso-flavonoids 함량이 증가 하다가 감소하는 형태를 보이며 초음파 조사세기 또한 필요이상으로 증가할 경우 오히려 사용된 에너지에 비해 iso-flavonoids 함량이 감소 하는 것으로 알 수 있었다.
3.4. 종합만족도 평가
반응표면분석법을 통해 독립변수인 초음파 조사시간, 주정/초순수 부피비, 초음파 출력에 따른 2개의 반응치인 추출수율과 iso-flavonoids 함량의 최적 조건을 Figure 5에 나타냈다. 두 반응치의 최댓 값을 모두 만족하는 최적조건은 초음파 조사시간(24.75 min), 주정/초 순수 부피비(39.75 vol%), 초음파 조사세기(592.36 W)이다. 최적 반응 치는 추출수율의 경우 27.8363 wt%이고 iso-flavonoids 함량은 30.8064 μg/mg이다. 최적화 그래프에서 또한 최적포인트 기준으로 증 가, 감소 기울기 값이 초음파 조사시간에서 가장 크게 나타났으며 이 는 반응치에 가장 영향력 있는 독립변수라는 것을 의미한다. 종합만 족도를 나타내는 D는 0부터 1로 나타내며 반응치의 적합성을 뜻한다. 즉 D 값이 1에 가까울수록 높은 만족도를 나타내며 0에 가까울수록 낮은 만족도를 나타낸다. 본 그래프에 따른 종합 만족도 D는 0.8938 로 높게 나타났다[30]. 또한 5% 이내의 유의성 수준에 인정된다. 반응 표면분석법을 통한 최적조건이 기초실험을 통해 설정한 변수범위 안 에 들었으므로 성공적인 실험이라고 볼 수 있으며 두 개의 반응치를 같이 분석해야 하는 다중분석 반응표면 분석임에도 불구하고 정확성 이 높게 나타났다.
4. 결 론
본 연구에서는 칡에 포함된 항산화 물질인 iso-flavonoids 성분을 추 출하기 위해 초음파에너지를 이용한 추출공정을 수행했다. 기초실험 을 통해 높은 반응치를 갖는 각 변수의 적정범위로 초음파 조사시간 (20.00~30.00 min), 주정/초순수 부피비(35.00~45.00 vol%), 초음파 조사세기(550.00~650.00 W) 범위를 설정하였다. 3개의 실험인자는 초음파 조사시간(x1), 주정/초순수 부피비(x2), 초음파 조사세기(x3)로 설정하고 20번의 표면설계를 생성했고 반응치인 추출수율(Y1)과 iso-flavonoids 함량(Y2)을 반응표면분석법에 의해 최적화한 결과 다음 과 같은 결론을 얻을 수 있었다.
-
칡의 추출수율에 따른 변수들의 교호효과도는 (초음파 조사시간 × 초음파 조사세기)와 (주정/초순수 부피비 × 초음파 조사세기)에 비 해 (초음파 조사시간 × 주정/초순수 부피비)가 추출수율에 더 큰 영향 을 미치는 것을 알 수 있었다. 또한 추출수율에 따른 각 변수의 P-value 중 x1 = 0.187로 가장 작은 값을 가지며 이는 가장 큰 주효과도를 나타 낸다.
-
Iso-flavonoids 함량의 경우에는 (초음파 조사시간 × 주정/초순수 부피비)와 (주정/초순수 부피비 × 초음파 조사세기)가 (초음파 조사시 간 × 초음파 조사세기)에 비해 더 많은 영향을 미치는 것을 알 수 있 었다. 또한 iso-flavonoids 함량에 따른 각 변수의 P-value는 x1 = 0.035, x2 = 0.971, x3 = 0.809로 x1의 P-value가 가장 작게 나타났으며 이는 초음파 조사시간(x1)이 가장 큰 영향을 미치는 것을 의미한다.
-
반응표면분석법을 이용하여 두 개의 반응치인 추출수율과 iso-flavonoids 함량의 최대값을 갖는 최적조건을 평가한 결과 초음파 조사시간(24.75 min), 주정/초순수 부피비(39.75 vol%), 초음파 조사세 기(592.36 W)에 따라 추출수율(27.8363 wt%)과 iso-flavonoids함량 (30.8064 μg/mg)의 결과를 얻을 수 있었다. 이때 종합 만족도 D는 0.8938로 높게 나타났으며 5% 이내의 유의성 수준에 인정된다. 또한 최적화 과정과 주효과도를 분석한 결과 초음파 조사시간이 두 반응치 모두에 가장 많은 영향을 미치는 인자임을 확인하였다.