활성화에너지 분석을 통한 마이야르 반응 기반 TTI (Time-Temperature Indicator)의 고온 건조에서의 수분변화 가시화 가능성 검토

마이야르 반응 물질은 Lee et al. (2018c)에 따라 높은 반응성, 많은 반응물(멜라노이딘) 생성을 통한 명확한 색상변화, 반응물의 가격 등을 고려하여 환원당으로서는 D-xylose (Samchun Chemical Co., Ltd., Seoul, Korea) 아미노산으로서 glycine (Samchun Chemical Co., Ltd., Seoul, Korea)를 사용하였다.

색변화를 측정하기위한 반응용액은 D-xylose 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 M, glycine은 0.5, 1.0, 1.5, 2.0 M를 제조하여 각 농도별로 체적비 1:1로 혼합하여 총 24종류의 혼합물을 제조하였다. 제조된 반응용액은 96-Well 마이크로 플레이트에 농도별로 투입하고 반응용액 표면에 미네랄 오일 층을 만들어 반응액의 수분 증발을 방지하였다.

또한 반응액은 혼합과 동시에 반응을 시작하므로, 96-Well 마이크로 플레이트에 반응액이 투입되고 10분 이내에 60, 70, 80°C의 실험용 건조기(OF3-15, JEIO TECH, Daejeon, Korea)에 넣어 반응을 시작하였으며, 12시간 동안 측정하였다. Fig. 1은 70°C에서 146분 동안 반응시킨 용액에 대한 마이크로 플레이트의 이미지로서 Fig. 1에 나타낸 바와 같이 A-D행, 1-6열에는 1열부터 D-xylose 1.0-6.0 M 농도별 glycine 0.5 M의 혼합용액이 투입되었으며, A-D행, 7-12열에는 7열부터 D-xylose 1.0-6.0 M의 각 농도에 따라 glycine 1.0 M의 혼합용액, E-H행의 1-6열 및 7-12열에는 각각 D-xylose 농도에 따라 glycine 1.5, 2.0 M의 혼합용액을 투입하여 동일한 온도에서 동시에 반응시켰다. 각 조건별 반응용액은 수분 증발에 따른 색변화의 오차 제거 및 실험반복을 위해 Fig. 1에 나타낸 바와 같이 4개의 웰에 같은 농도의 용액을 넣어 반응시켰다.

마이야르 반응 용액의 반응시간 경과에 따른 색변화의 분석을 위해 반응용액이 담긴 마이크로 플레이트를 30분 간격으로 디지털카메라(iPhone SE 2nd generation, ver. 14.6, Apple Inc., CA, USA)를 통해 촬영하였으며, 동일한 광원 및 위치(마이크로 플레이트 중심 및 중심위치 기준 150 mm 높이)에서 촬영하여 광원에 의한 차이를 최소화하였다. 촬영된 이미지는 RGB 색공간(R (red), G (green), B (blue))을 기반으로 분리 및 각 RGB 채널별로 시간경과에 따른 색상 변화를 나타내었다. 색상값은 각 이미지에서 광원의 반사에 의한 오류부분을 제외한 픽셀에 대해 평균값을 계산하여 표현하였다.

농산물의 고온 건조 중 수분변화의 가시화가 가능한 마이야르 반응기반 TTI의 활용가능성을 검토하기 위해 박층 사과를 고온 건조하여 수분변화를 측정하고 마이야르 반응에 따른 색상변화와 비교 검토하였다. 건조시료로 사용될 사과는 2024년 1월 청주지역의 청과물매장에서 부사 품종을 구입하였으며, 4°C의 일반형 냉장고(R-B425GB, LG, Seoul, Korea)에 저장하여 사용하였다. 사과의 초기 함수율은 상압정온건조법(MFDS, 2024)을 이용하여 측정하였으며, 습량기준 함수율(w.b.) 85±0.5%로 나타났다. 건조시료는 사과를 4등분하여 껍질과 씨를 제거하고, 0.5 mm 두께로 박층 절단하였다. 절단된 시료는 수분분석기(MB95, OHAUS Corp., NJ, USA)를 이용하여 시료중량 20, 40, 60 g, 건조온도 60, 70, 80°C의 총 9개 조건으로 건조하고, 5초 간격으로 중량변 화를 측정하였다. 각 조건에 따라 3회 반복 측정 후 평균값으로 나타냈으며, 반복에 따른 차이는 나타나지 않았다. 함수율 변화는 측정된 중량을 기반으로 식(1)을 통해 함수율로 환산하여 나타내었다.

여기서,

마이야르 반응 용액의 색변화 및 사과의 건조 함수율의 측정결과로 나타난 곡선의 동역학적 특성을 나타낼 수 있는 모델식을 검토하기 위해 대표적인 동역학적 모델식에 적용하여 적합도를 분석하였다. 분석에 사용된 동역학적 모델식으로서 벼, 옥수수 등의 박층 건조모델로서 많이 활용되는 Page (1949) 건조모델, Newton의 냉각법칙을 기반으로 박층건조의 모델로 사용되며, 식품의 소비기한 측정에 있어 비선형 품질곡선에 대한 모델식으로도 활용되는 Lewis (1921) 건조모델, 개체군의 성장 및 변화율을 S자 형태로 표현하는 로지스틱(Verhulst, 1838) 모델을 각각 적용하여 비교하였으며, 활용된 식은 아래와 같다(Faal et al., 2015;Park et al., 2016;MFDS, 2023).

여기서,

VS (value scaler)는 데이터의 모델링에서 스케일의 차이에 따른 모델링 오차를 방지하기 위해 색변화 및 함수율 데이터에 대해 값의 범위를 0-1 사이의 값으로 조정한 것으로서, 아래의 식을 이용하여 환산하였다.

여기서,

VS (value scaler)는 건조모델에서는 본래 함수율비(MR, moisture ratio)로서 나타내고 일반적으로 아래의 식(6)으로서 나타내나, 평형함수율을 나타내는 M_e 값의 적용에 따른 차이가 적어 식(7)로서 간소화할 수 있다(Faal et al., 2015). 또한 건조모델뿐만 아니라 마이야르 반응용액의 색변화에 대한 모델식으로서 활용되어야 하므로, VS로서 표현하였다.

여기서,

마이야르 반응 용액의 색변화 및 사과의 건조 함수율 변화에 대한 동역학적 모델식의 적합성은 아래의 식(8)을 통해 결정계수(R²)를 산출하고, 결정계수가 높을수록 높은 적합성을 보이므로, 결정계수가 가장 높은 동역학적 모델식을 적합한 모델식으로 선정하였다.

여기서,

활성화에너지(Ea)는 반응속도에 대한 온도의존도를 나타내는 아레니우스 방정식(Arrhenius equation)인 식(9)를 이용하여 기울기(Ea/R) 및 사전지수인자 D를 선형회귀를 통해 산출하고, 산출된 기울기에 기체상수를 곱하여 활성화에너지 값을 산출하였다.

여기서,

색상분석, 동역학적 특성 분석 및 활성화에너지 분석을 위해 Python (ver. 3.12) 및 통합개발환경으로서 PyCharm Community Edition (ver. 2024. 1.1, JetBrains, Czech Republic.)을 기반으로 Python 라이브러리로 NumPy (ver. 1. 26. 4), Pandas (ver. 2.2.2), Scikit-learn (ver. 1.4.2)을 이용하였다.

마이야르 반응용액의 색상변화는 높은 가시성을 나타냈으며, 온도 및 용액의 농도 증가에 따라 반응속도가 증가하였다. Fig. 2는 분석조건에 따른 색상값의 변화를 비교하여 나타낸 것으로 (a) RGB 채널에 따른 비교, (b) 온도에 따른 비교, (c) D-xylose 농도에 따른 변화, (d) glycine 농도에 따른 변화를 나타낸 것이다. Fig. 2(a)에 보이는 바와 같이 마이야르 반응 용액의 색상변화는 명확한 색상 변화를 나타냈으며, 동일한 농도조건에서도 RGB 채널에 따른 색상값 변화의 차이를 보였다. 마이야르 반응용액의 색상변화는 무색에서부터 노란색, 갈색, 흑갈색, 검은색과 같은 색상변화를 보여 육안으로도 명확한 차이를 확인 가능하였다. 이러한 결과는 마이야르 반응의 결과물로 멜라노이딘이라고 불리는 갈색물질이 생성되고, 갈색물질의 생성량에 따라 황색의 농도가 짙어지며 나타난 것으로 판단된다. 따라서 마이야르 반응 기반으로 제작된 TTI에 있어서 육안을 통해 TTI의 색상변화를 확인하는 것으로 농산물 건조물에서 육안으로 알 수 없는 수분변화에 대해 예측 가능할 것으로 판단되었다.

한편, Fig. 2(a)는 RGB Channel에 따른 차이를 나타내기 위해 대표적으로 반응온도 70°C에서 D-xylose 1.0 M과 glycine 2.0 M 혼합용액의 색상변화를 표현한 것으로 같은 반응조건임에도 color channel에 따라 변화속도의 차이를 보였다. Red channel의 변화속도가 가장 느리고, blue channel의 변화가 가장 빠르게 나타났으며, color channel의 값이 20에 도달하는 시간을 고려하면, red channel은 약 10시간, green channel은 6시간, blue channel의 경우는 약 3시간 수준으로 판단되었다. 이러한 결과는 육안을 통해서도 색상의 변화에 대해 확인할 수 있으나, 디지털 이미지의 분석을 통한 RGB channel의 분석은 더욱 다양한 반응속도 및 명확한 반응시간 판단이 가능하다는 것을 시사할 수 있다. 즉, 피건조물의 수분변화에 대한 예측을 육안을 통한 TTI의 관찰로 단순히 노란색, 갈색, 검은색 등으로 판단하는 것보다 각 color channel 값에 대해 측정값 또는 최소값 수준이 나타나는 시간을 기준으로 판단하는 것이 더욱 명확한 분석 및 예측이 가능할 것으로 판단된다. 또한 피건조물의 건조완료에 대한 판단에 TTI를 활용하는 것을 고려하면, 육안을 통해서는 색상변화가 검은색수준이 되는 것을 기준으로 할 수 있으나, color channel 값을 기준으로 판단하면 red, green, blue channel에 따라 특정 임계값을 설정하거나 각 channel 값이 최소값에 도달하는 수준으로 판단하면 더욱 다양한 시간대에 대한 활용이 가능할 것으로 사료된다.

Fig. 2(b)는 반응온도 따른 차이를 나타내기 위해 대표적으로 D-xylose 3.0 M과 glycine 2.0 M 혼합용액에 대한 색상변화에서 red channel의 변화를 표현한 것으로 반응온도가 높을수록 기울기가 증가하고, 반응속도가 빠른 경향을 보였다. 이러한 경향은 마이야르 반응의 반응온도가 낮은 10°C이하의 저온에서도 반응온도가 높을수록 반응속도가 빠른 경향을 보이는 것으로 보고되어 있다(Sakai et al., 2020).

Fig. 2(c) 및 (d)는 반응농도에 따른 차이를 비교하기위해 각각 반응온도 60°C, glycine 1.5 M 농도에서 D-xylose 농도에 따른 green channel의 변화 및 반응온도 80°C, D-xylose 2.0 M 농도에서 glycine 농도에 따른 blue channel의 변화를 나타낸 것이다. Fig. 2(c) 및 (d)에 보이는 바와 같이 농도가 높을수록 반응속도가 증가하는 경향을 보였으며, 반응속도의 증가율은 농도가 높아짐에 따라 감소하는 것으로 나타났다. 주된 기울기를 보이는 직선구간에 대해 단순선형회귀를 통한 기울기를 분석한 결과, Fig. 2(c)의 반응온도 60°C, glycine 1.5 M 농도에서 D-xylose 농도가 증가할 수록 기울기가 각각 8.9, 13.4, 17.9, 26.1, 30.8, 39.2로 나타났다. 또한 반응온도 80°C, D-xylose 2.0 M 농도에서 glycine 농도에 따른 기울기는 22.2, 53.2, 107.9, 106.7을 보였다. 즉, 높은 농도에서 높은 기울기가 나타나 최소값이 나타나는 시간이 단축되었으며, 농도의 증가에 따라 기울기의 증가율도 증가하는 경향을 보이다가 D-xylose 농도 3.0 M와 4.0 M 및 glycine 농도 1.0 M와 1.5 M에서 최대 증가량을 보인 후 다시 감소하였다. 특히, Fig. 2(d)에서 glycine 1.5 M과 2.0 M의 변화곡선이 매우 유사하여 큰 차이가 없었다. 이러한 결과에 따라 동일한 반응온도에 있어서도 농도의 조절에 따라 다양한 반응속도를 얻을 수 있으므로, 마이야르 반응 기반의 TTI 제작에 있어 농도의 적절한 조절을 통해 다양한 건조속도를 가지는 농산물의 건조에 있어서도 적용 가능할 것으로 사료되었다.

Table 1은 마이야르 반응용액의 반응조건(온도, 농도)에 따른 반응 완료시간을 나타낸 것으로 각 RGB channel 값이 20 이하에 도달한 시간을 기준으로 반응 완료시간을 판단하였다. Table 1에 보이는 바와 같이 반응 완료시간은 1시간에서 12시간 사이의 다양한 시간을 보였으며, (-)로서 표시된 부분은 RGB channel 값이 반응 12시간 후에도 20이하의 값이 나타나지 않은 조건으로 실제로는 12시간 이상 반응시간이 소요될 것으로 예상되므로 농도의 조절에 따라 1시간에서 12시간 이상 수분변화의 표현이 가능한 TTI의 제작이 가능할 것으로 판단되었다. 단, 가장 낮은 온도인 60°C에서는 가장 높은 농도에도 불구하고 red, green, blue channel에 따라 각각 6.1, 4.1, 2.6시간을 보여 color channel의 선택에 따라 측정가능한 최소시간에 대한 차이가 존재하였으며, 최소시간 이내의 예측 및 분석을 위해서는 검은색이 아닌 짙은 갈색 수준으로 종료시간을 판단하거나 임계치의 값을 20보다 높이는 방안이 필요할 것으로 사료되었다.

마이야르 반응 기반의 TTI의 실질적인 활용가능성을 확인하기 위해 마이야르 반응용액 색상변화에 대한 결과를 실제 농산물의 건조 시 수분변화 특성과 비교하였다. Fig. 3은 건조온도 및 피건조물의 중량에 따른 함수율의 변화곡선을 나타낸 것으로 각 건조온도 (a) 60°C, (b) 70°C, (c) 80°C에 따른 비교 및 (d) 건조중량 40 g의 피건조물에 대해 건조온도에 따른 차이를 표현하였다. Fig. 3에 보이는 바와 같이 건조중량이 가볍고, 건조온도가 높을 수록 건조속도가 빠르고 건조함수율 최저 수준까지 도달하는 시간도 단축되는 것으로 나타났다. 박층사과의 건조에 있어 건조종료시간을 함수율 0.2 수준으로 판단한 경우, Fig. 3에 보이는 바와 같이 건조시간은 2시간에서 16시간 정도의 범위로 나 타났으며, 이러한 범위는 마이야르 반응 용액의 색상변화에 대한 반응종료시간과 비교하여도 비슷한 범위를 보였다. 또한 동일한 온도에 대해 비교한 경우에도, 측정된 건조종료 시간에 대한 값을 마이야르 반응 용액에 대한 농도 조절 및 색상변화에 대한 color channel의 선택에 따라 충분히 예측 가능할 것으로 판단되었다. 다만, 건조물의 중량을 고려하였을 때 중량의 증가에 따라 건조시간도 증가하므로 60 g 이상의 피건조물에 대해서는 마이야르 반응기반의 TTI를 활용하지 못할 것으로 판단될 수 있다. 그러나 Park et al. (2016)은 함수율 28%의 콩을 45°C로 열풍건조한 경우 건조시간 300분 후 10% 함수율을 나타낸 것으로 보고하였으며, Lee et al. (2017)은 초기함수율 68%의 증자된 고구마 4 kg를 25%까지 건조하는데, 가장 느린 건조시간이 나타난 조건인 건조두께 10 mm, 건조온도 60°C, 송풍속도 0.6 m/s의 조건에서도 417분이 나타났다고 발표하였다. 따라서 이와 같은 연구와 비교하면 건조방식 및 피건조물의 전처리 상태에 따라 중량 또는 조건이 변화하더라도 60, 70, 80°C의 건조온도 범위에서 건조시간이 마이야르 반응 기반 TTI의 반응시간 범위 이내이거나 유사한 결과를 보여 고온의 온도범위에서의 활용가능성이 높을 것으로 사료될 수 있다.

아레니우스 방정식 활성화에너지는 반응이나 변화에 대한 온도의존성을 검토할 수 있으며, 온도 및 시간의 적산값에 대해 나타내는 TTI에 있어 예측 정확도의 판단을 위해 주로 사용되는 방법이다(Yang & Xu, 2021). 즉, 농산물 품질 변화에 있어서의 활성화에너지와 TTI의 색상변화에 대한 활성화에너지가 비슷한 경우 TTI의 색상변화를 통해 농산물의 품질변화에 있어서도 예측 가능한 것으로 판단한다. 따라서 마이야르 반응용액의 색상변화와 박층 사과의 함수율변화에 대한 활성화에너지 분석 및 비교가 필요하나, 활성화에너지의 산출을 위해서는 비례상수 K 값이 필요하므로 곡선의 동역학적 특성에 대한 검토가 필요하다. 마이야르 반응용액의 색변화에 대한 동역학적 특성은 적합한 비례상수 K의 값을 얻을 수 있도록 충분히 색상이 변화된 조건 즉, VS 값이 0.4 이하까지 감소한 조건에 대해서만 분석하였다.

Fig. 4는 각 모델식에 따른 실제 값과의 적합도를 확인하기 위해 측정값과 각 모델에 대한 곡선을 비교하여 나타낸 것이다. Fig. 4(a)는 마이야르 반응용액의 색변화에 대한 것으로 대표적으로 건조온도 70°C, D-xylose 2.0 M, glycine 1.0 M 조건에 대한 green channel의 색상변화에 대한 측정값 곡선과 선형회귀분석을 통해 산출된 값으로 표현한 각 모델의 예측곡선을 비교하여 나타낸 것이다. 또한 Fig. 4(b)는 박층 사과의 수분변화에 대한 것으로 대표적으로 건조온도 70°C, 건조중량 60 g의 조건으로 건조된 박층사과의 함수율 변화에 대한 측정값의 곡선과 각 모델의 예측곡선을 비교하여 표현하였다. Fig. 4에 보이는 바와 같이 측정값에 대한 곡선과 비교하여 Page model 및 Logistic model의 경우는 유사한 곡선을 나타내었으나 Lewis model은 곡선의 유사도가 비교적 낮은 경향을 보였다. 결정계수(R²)의 결과에서도 Page model 및 Logistic model의 결정계수는 0.97 이상의 매우 높은 값을 나타내었으나. Lewis model은 0.80 수준의 낮은 값이 나타나 적합성이 낮았다. Table 2는 마이야르 반응 용액의 색변화와 박층사과의 함수율 변화에 대해 각 동역학적 모델식을 적용한 경우 나타난 결정계수(R²) 값을 나타낸 것으로 Table 2에 보이는 바와 같이 Lewis model의 결정계수가 가장 낮은 값을 보였고, Page model 및 Logistic model의 결정계수는 유사한 수준의 높은 값이 나타났다. Lewis model은 평균 0.7974 였으며, Page model 및 Logistic model의 평균값은 각각 0.9834, 0.9859의 값을 보였다. 이러한 경향은 실제 측정값의 변화속도 및 변화 곡선의 형태에 의한 것으로 판단되며, 반응 및 건조 초기에는 측정값의 변화가 느리게 진행되다가 점차 변화속도가 증가하여 일정한 속도로 변화가 진행된 후 변화속도가 다시 느려지는 역 S 형태의 곡선이 나타나고 있으나 Lewis model 식은 변화 초기의 변화에 대해 명확하게 표현하지 못하고 측정값과의 차이가 발생하여 전체적으로 적합성이 낮아진 것으로 사료된다. Page model와 Logistic model을 비교할 경우 전반적으로 높은 적합도 및 결정계수를 보였으나 Page 모델의 경우 Fig. 4(b)에 나타낸 것처럼 낮은 적합성을 보인 일부 조건이 존재하였으며 Table 2의 건조온도 70°C, 시료 중량 60 g의 조건에서 Logistic model은 0.9892의 값을 보였으나 Page model에서는 0.9672의 다소 낮은 값을 보여 전체적인 조건에 대해 높은 적합성을 보인 Logistic model이 더 적합한 것으로 판단되었다.

따라서 활성화에너지 분석을 위한 비례상수 K는 Logistic model을 기준으로 산출되었으며, 기울기에 대한 부분인 실험상수 B 값을 기준으로 활성화에너지를 산출하였다. Table 3은 마이야르 반응용액의 색상변화에 대한 활성화에너지 산출 결과이며, Table 4는 박층 사과건조에서 수분변화에 대한 활성화에너지이다. Table 3에 나타낸 바와 같이, 마이야르 반응용액의 색상변화에 대한 활성화에너지의 범위는 red channel, green channel, blue channel에 따라 각각 43-111, 56-107, 27-113 kJ/mol·K의 범위를 나타냈으며, 박층 사과건조에서 수분변화에 대한 활성화에너지와 유사한 범위를 보였다. 따라서 마이야르 반응 용액은 사과건조의 수분변화와 유사 한 온도의존성을 가지는 것을 알 수 있으며, 이는 마이야르 반응 기반의 TTI의 활용은 박층 사과건조의 수분변화의 가시화 및 예측에 있어 충분히 활용 가능할 것을 시사하였다. 또한 Taoukis et al. (2001)은 활성화에너지의 차이가 ±25 kJ/mol·K 이내인 경우 품질변화의 오차는 15% 이하인 것으로 보고하였으며, 이를 바탕으로 고려할 경우 선택 가능한 농도의 범위는 더욱 증가할 것으로 사료된다.

이와 같은 결과를 바탕으로 박층 사과의 각 건조중량에 따른 적절한 농도조건 및 color channel을 박층 사과의 중량 별 건조 종료시간, 각 RGB color channel의 색상변화 및 활성화에너지 값의 유사도를 기준으로 하여 검토하였다. 박층 사과의 건조종료시간은 함수율 0.2 이하로 판단하였으며, 각 RGB color channel의 색상변화는 20 이하의 값이 나타난 시점을 검은색이 나타난 반응종료 시간으로 고려하였고, 활성화에너지의 유사도는 Taoukis et al. (2001)에 따라 ±25 kJ/mol·K 이내인 수준으로 판단하였다. 따라서 활성화에너지의 차이가 ±25 kJ/mol·K인 조건 중 건조 및 색상변화의 종료시간이 유사한 조건을 기준으로 사과의 각 건조중량에 따른 적절한 농도조건 및 color channel을 검토하였다. 그 결과, 건조중량 20 g에서는 D-xylose 4.0 M, glycine 2.0 M 용액에 대한 색상변화의 blue channel을 이용하는 것이 적절할 것으로 판단되었으며, 건조중량 40 g 및 건조중량 60 g에서는 각각 D-xylose 3.0 M, glycine 1.0 M의 green channel 및 D-xylose 1.0 M, glycine 2.0 M의 red channel을 활용할 경우 수분변화에 대해 높은 정확도로 예측이 가능할 것으로 사료되었다.