1. 서론
세팔로마닌(cephalomannine)은 파클리탁셀(paclitaxel)과 함께 주목(yew tree) 에서 발견된 대표적 디테르페노이드 계 항암물질(Fig. 1)로 난소암, 유방암, 폐암, 백혈병 등의 치료에 효과적이다[1]. 이들 성분들은 주로 주목에서 직접 추출하거나 합성 또는 식물세포배양법으로 생산되는데, 이중 식물세포배양법은 외부인자(기후, 환경)의 영향을 받지 않고 생물반응기에서 안정적인 생산이 가능하기 때문에 일정한 품질의 제품을 대량 생산할 수 있다[2,3]. 하지만 식물세포배양으로부터 생산된 세팔로마닌과 파클리탁셀의 구조와 물리적 특성의 유사성 때문에 이들 두 성분의 분리는 상당히 어렵다[4]. 세팔로마닌과 파클리탁셀의 대표적 분리 방법으로 크로마토그래피법 [1], 액-액 추출법[4], 오존처리법[5]이 있다. 크로마토그래피법은 분리효율은 우수하나 처리량이 적고 조업비용이 많이 발생하고, 액-액 추출법은 조업은 간단하나 한정된 추출용매와 낮은 선택도로 분리 효율이 떨어진다. 또한 오존처리법은 효과적이고 경제적이나 조업이 복잡하고 유용성분을 산화시킬 우려가 있다. 따라서 식물세포배양법에 의해 생산된 이들 성분의 효율적 분리에 대한 연구가 여전히 요구되고 있다. 최근 식물세포배양으로부터 세팔로마닌/파클리탁셀 분리를 위하여 여러 종류의 분리/정제 공정을 도입한 결과, 흡착 공정에서 세팔로마닌과 파클리탁셀이 가장 효과적으로 분리되었다[6]. 흡착 공정에서 안정성이 확보된 상용흡착제로 활성백토, 활성탄, 실로퓨트(Sylopute) 등이 사용될 수 있는데, 이중 실로퓨트가 다른 흡착제들보다 흡착 효율이 우수하고 흡착 후 여과 속도가 빨라 가장 효과적이다[7]. 하지만 기존 문헌들[6,7]은 흡착 공정 변수의 영향 및 최적화에 대한 정성적 연구와 파클리탁셀 흡착 특성 연구에 집중되어 있으며, 세팔로마닌의 흡착 특성에 관한 정량적인 연구는 전무한 실정이다. 따라서 본 연구에서는 상용흡착제 실로퓨트에 의한 세팔로마닌의 흡착 특성을 정량적으로 조사하였다. 대표적 등온흡착식(Langmuir, Freundlich, Temkin, and Dubinin-Radushkevich isotherms)을 이용하여 흡착 평형 데이터에 가장 적합한 모델을 선정하였다. 또한 속도론적 해석을 위하여, 유사 일차 및 이차 반응속도식과 입자 내 확산 모델식을 이용하였다. 동시에 엔탈피 변화, 엔트로피 변화, 깁스 자유에너지 변화, 활성화 에너지, 등량 흡착열을 계산하여 실로퓨트에 의한 세팔로마닌 흡착에 대한 열역학적 해석을 수행하였다.
Fig. 1. The chemical structure of paclitaxel (A) and cephalomannine (B).
2. 재료 및 방법
2-1. 세팔로마닌 시료
본 흡착 실험에 사용된 세팔로마닌(순도: 98%)은 식물세포 Taxus chinensis 배양액으로부터 정제 되었으며[8], ㈜삼양바이오팜으로부터 제공받았다.
2-2. 세팔로마닌 분석
세팔로마닌함량분석을위해 HPLC 시스템(SCL-10AVP, Shimadzu, Japan)과 Capcell Pak C18 column(250 × 4.6 mm, Shiseido, Japan)을 사용하였다. 이동상은 아세토니트릴(acetonitrile)과 증류수(distilled water) 혼합용액(35/65~65/35, v/v, gradient mode)이며 유속은 1.0 mL/min이었다. 시료 주입량은 20 μL이며 227 nm에서 UV에 의해 검출하였다[9].
2-3. 등온흡착평형 실험
메틸렌 클로라이드에 세팔로마닌(농도: 2,000~3,250 mg/L)을 녹이고 흡착제 실로퓨트(Fuji Silysia Chemical Ltd., Japan) 20 g/L를 첨가하여 자력 교반 항온조(EYELA, PS-1000)에서 온도(293, 303, 308, 313, 318 K)를 달리하여 30분 동안 등온 흡착 하였다. 흡착제 실로퓨트의 기공 지름, 기공 부피, 표면적은 각각 40~60 nm, 1.81 cm3 /g, 331 m2 /g이다[10]. 흡착 후 여과하고 여과액을 313 K, 감압상태에서 건조하였으며 세팔로마닌 함량은 HPLC로 분석하였다. 용액에서의 흡착질 평형 농도(Ce)에 따른 평형 흡착 용량(qe)의 변화에 대한 실험 데이터를 네 종류의 등온흡착식 Langmuir, Freundlich, Temkin, Dubinin-Radushkevich에 각각 적용하였다[10-12]. 모델 유효성을 확인하기 위해 통계 방법으로 결정계수(coefficient of determination, r2)와 평균 제곱근 표준편차(root mean square deviation, RMSD)를 이용하였다[13].
2-4. 흡착속도 실험
메틸렌 클로라이드에 세팔로마닌(3,250 mg/L)을 녹이고 실로퓨트 20 g/L를 첨가하여 자력 교반 항온조(EYELA, PS-1000)에서 온도(293, 303, 308, 313, 318 K)를 달리하여 흡착 시간 1, 3, 6, 12, 18, 24, 30분 동안 각각 실험을 수행하였다. 흡착 후 여과로 얻은 여과액 을 313 K, 감압상태에서 건조하였으며, 세팔로마닌 함량은 HPLC로 분석하였다. 각 온도 별 평형 데이터를 유사 일차 및 이차 반응속도식과 입자 내 확산 모델식에 적용하였다[12-14]. 또한 정규화된 표준편차(normalized standard deviation), Δq(%)를 이용하여 모델의 유효성을 확인하였다.
2-5. 열역학적 해석
열역학적 데이터를 vant’s Hoff 식, Arrhenius 식, Clausius-Clapeyron 식에 각각 적용하여 표준 엔탈피 변화(DH0), 표준 엔트로피 변화(S0 ), 표준 깁스 자유에너지 변화(DG0 ), 활성화 에너지(Ea), 등량 흡착열(DHX)을 계산하였다[10,11,13].
3. 결과 및 고찰
3-1. 등온흡착식 해석
본 연구에서는 흡착 온도(293, 298, 303, 308, 313 K)에 따른 등온 흡착평형 실험을 수행하여 실로퓨트에 의한 세팔로마닌의 흡착 양상을 조사하였다. 흡착평형농도가 증가함에 따라 흡착량이 증가하다 점차평형에 도달하였다(Fig. 2). 흡착 특성을 조사하기 위하여, 흡착 평형 데이터를 선형의 Langmuir, Freundlich, Temkin, Dubinin-Radushkevich 등온흡착식에 각각 적용하여 파라미터를 계산하고 이를 Table 1에 정리하였다.
Fig. 2. Adsorption isotherm for cephalomannine onto Sylopute at different temperatures.
Table 1. Langmuir, Freundlich, Temkin, and Dubinin-Radushkevich isotherm constants for the adsorption of cephalomannine onto Sylopute
Langmuir 식으로부터 흡착 온도(293, 298, 303, 308, 313 K) 증가에 따라 최대흡착용량(qmax: 139.86, 138.89, 137.93, 137.17, 136.43 mg/g)은 감소하였다. 이는 온도 증가에 따라 용매에 대한 흡착질의 용해도가 증가하게 되고, 이로 인해 흡착질-흡착제 보다 흡착질-용매간의 상호작용이 더 강해지기 때문이다[10]. 이러한 결과는 동일 흡착제에서 온도 증가에 따라 최대흡착용량이 증가하는 파클리탁셀 흡착(500.0~625.0 mg/g at 298~318 K)과는 다른 양상을 보였으며, 주어진 온도에서 파클리탁셀의 qmax보다 작은 값을 보였다[14]. Langmuir 상수(KL)는 0.148 > 0.137 > 0.124 > 0.114 > 0.103 L/mg으로 온도 증가에 따라 그 값이 감소하였다. 무차원 분리계수 RL은 범위에 따라, 0L<1인 경우 흡착 공정이 적합함을, RL>1인 경우 흡착 공정이 적합하지 않음을 나타낸다[15]. Table 1에서 보는 바와 같이 모든 온도에서 RL값(0.002077~0.004816)이 0과 1사이에 존재하므로 흡착공정이 적합함을 알 수 있었다. Freundlich 식으로부터 흡착 온도(293~313 K)가 증가함에 따라 흡착 용량과 관련된 Freundlich 상수(KF: 74.080~70.307 (mg/g)(L/mg)1/n)는 감소하였다. 흡착 강도와 관련된 무차원 분리계수(1/n)는 흡착 공정의 유효성 판단의 지표로 사용된다. 즉, 1/n 값이 0과 1 사이에 존재할 때 흡착공정이 적합함을 나타내는데, 이는 흡착제 표면부위(surface site)의 분포 또는 표면 밀도의 증가에 따른 흡착제-흡착질간의 상호작용의 감소를 유발하는 요인에 기인한다[12]. 본 연구의 경우, 모든 온도에서 1/n 값이 0.112~ 0.113을 나타내어 흡착 공정이 적합함을 알 수 있었다. Temkin 식으로부터 흡착 온도(293~313 K) 증가에 따라 흡착열 대응 상수(B: 13.541~13.096 J/mol)는 감소하므로 흡착 과정에서 에너지가 방출되는 발열 반응임을 알 수 있었다[16]. 또한 낮은 흡착열(B<20 J/mol)은 물리적 흡착을 의미하며, 흡착제-흡착질간의 약한 결합(van der Waals 결합)이 존재함을 알 수 있었다[12]. 평형결합상수(KT)는 온도 증가에 따라 107.61 > 103.46 > 95.70 > 86.98 > 86.11 mg/L로 감소하여 흡착 온도가 증가함에 따라 흡착질과 흡착제의 결합력이 감소함을 알 수 있었다. Dubinin-Radushkevich 식으로부터 흡착 온도 (293~313 K) 증가에 따라 이론포화흡착용량(qD: 134.14~130.76 mg/g)은 감소하여 흡착 온도가 높을수록 흡착이 원활하지 않음을 의미한다[10]. 흡착에너지(E)는 온도 증가에 따라 0.187 > 0.173 > 0.156 > 0.140 > 0.128 kJ/mol로 E<8 kJ/mol 영역에 존재하므로 실로퓨트에 의한 세팔로마닌의 흡착이 물리적 흡착임을 재확인 하였다[12,15]. 네 종류의 등온흡착식을 비교한 결과, Temkin 등온흡착식의 경우 가장 높은 r2(>0.981)와 가장 낮은 RMSD (<2.009)로 실로퓨트를 이용한 세팔로마닌의 흡착에 가장 적합함을 알 수 있었다.
3-2. 속도론적·열역학적 해석
세팔로마닌 농도를 3,250 mg/L으로 고정하고 실로퓨트 20 g/L를 첨가하여 흡착시간에 따른 흡착량을 조사하였다. 흡착 시간에 따라 흡착량의 변화가 거의 없었는데, 이는 높은 순도(~98%)의 세팔로마닌이 실로퓨트(세공지름: 40~60 nm)의 기공으로 빠르게(1분 이내) 흡착되는 것으로 판단된다[10,17]. 실험값을 선형의 유사 일차 반응 속도식에 적용한 결과, 평형상태에서의 흡착량 qe 와 특정 흡착 시간에서의 흡착량 qt 의 차이가 미미하여 속도론적 해석이 불가하였다. 선형의 유사 이차 반응속도식의 경우, t/qt 와 t를 도식화하고(Fig. 3), 직선의 기울기와 절편으로부터 파라미터 값을 계산하였다(Table 2). 흡착 속도 상수 k2는 흡착 온도(293, 298, 303, 308, 313 K)가 증가할수록 0.276 > 0.255 > 0.242 > 0.239 > 0. 231 g/mg·min로 감소하였다. 모든 온도에서 높은 r2 (0.999)와 낮은 Dq(%)(0.650~0.908)를 가져 실로퓨트에 대한 세팔로마닌의 흡착 공정은 유사 이차 반응속도식을 잘 따름을 알 수 있었다. 흡착 공정에서 율속 단계를 확인하기 위하여, 입자 내 확산 모델에 적용하였다. 일반적으로 입자 내 확산 모델에 의해 흡착은 필름 확산, 입자 내 확산, 평형의 세 단계로 구분된다[12,18,19]. 흡착 속도 실험값을 입자 내 확산 모델에 적용한 결과(data not shown), 필름 확산과 입자 내 확산은 확인되지 않았으며 실로퓨트에 대한 세팔로마닌의 흡착이 매우 빠르게 평형에 도달함을 알 수 있었다. 즉, 필름 확산과 입자내 확산이 율속 단계에 영향을 미치지 않았다. 또한 평형 단계에서의 속도상수는 매우 작은 값(<0.414 mg/g·min1/2)을 보였다. 일반적으로 평형 단계의 속도 상수는 율속 단계에 비해 매우 작은 값(~780배)을 가지므로 반응속도에 미치는 영향은 미미하다[14]. 세팔로마닌과 파클리탁셀의 분자 구조가 매우 유사한데, 공간채움모형(space-filling model)에 근거한 파클리탁셀의 분자 크기(1.0×1.5×2.0 nm)를 고려하면 실로퓨트의 세공(지름: 40~60 nm) 내로 세팔로마닌이 원활히 확산 가능함을 알 수 있다[20].
Fig. 3. Pseudo-second-order plot for the adsorption of cephalomannine onto Sylopute at cephalomannine concentration 3,250 mg/L, Sylopute 20 g/L.
Table 2. Parameters of pseudo-second-order kinetic model for the adsorption of cephalomannine onto Sylopute at different temperatures
Arrhenius 식에 의해 ln k2 대 1/T를 도식화하고(Fig. 4(A)), 직선의 기울기로부터 계산된 Ea는 6.353 kJ/mol이었다(Table 3). Ea 값이 <40 kJ/mol인 경우 물리적 흡착, 40~800 kJ/mol인 경우 화학적 흡착을 나타내는데[21], 실로퓨트에 대한 세팔로마닌의 흡착은 물리적 흡착임을 알 수 있었다. van’t Hoff 식에 의해 ln Ke 대 1/T를 도식화하고 (Fig. 4(B)), 직선의 기울기와 절편으로부터 DH0 와 DS0를 구하고 이로부터 DG0 를 계산하여 Table 3에 나타내었다. DH0 는 -9.155 kJ/ mol로 음의 값을 나타내어 흡착 공정이 발열임을 확인하였다. DS0 는 음의 값(-0.017 kJ/mol·K)을 보여 흡착인 진행됨에 따라 무질서도가 감소함을 알 수 있었다. DG0는 온도 293, 298, 303, 308, 313 K에서 각각 -4.861, -4.710, -4.612, -4.578, -4.522 kJ/mol이었으며 모두 음의 값을 보여 흡착 공정이 자발적으로 수행됨을 알 수 있었다[22]. 일반적으로 DG0 가 -20 kJ/mol 이상이면 물리적 흡착인데[23], 본 흡착은 물리적 흡착임을 재확인할 수 있었다. 등량 흡착열 DHX을 조사하기 위하여, 흡착량 qe값을 100, 105, 110, 115, 120, 125, 130, 135 mg/g으로 설정하고 이에 해당되는 Ce 값을 Clausius-Clapeyron 식에 적용하여 ln Ce와 1/T를 도식화하였으며(Fig. 5(A)), 직선의 기울기로부터 DHX를 구하였다. 또한 흡착량(qe)에 따른 흡착에너지 분포를 확인하기 위하여, DHX와 qe를 도식화하여 Fig. 5(B)에나타내었다. 주어진 qe (100~135 mg/g)에서 DHX값은 -20 kJ/mol 정도로 거의일정하였는데, 이는 흡착량이 많아져도 흡착제-흡착질 상호작용이 존재하지 않아 흡착제 표면의 에너지는 균일함을 의미한다[16,24]. 또한 DHX가 80 kJ/mol 이하이면 물리적 흡착을 의미하여[16], 실로퓨트에 대한 세팔로마닌의 흡착이 물리적 흡착임을 재확인하였다.
Fig. 4. Plots of ln k2 versus 1/T (A) and ln Ke versus 1/T (B) for the adsorption of cephalomannine onto Sylopute at different temperatures.
Fig. 5. Plots of ln Ce versus 1/T (A) and DHX versus qe (B) for the adsorption of cephalomannine onto Sylopute.
Table 3. Thermodynamic parameters for the adsorption of cephalomannine onto Sylopute at different temperatures
4. 결론
본 연구에서는 상용흡착제인 실로퓨트에 대한 세팔로마닌의 흡착 공정에 대한 등온흡착식 및 속도론적·열역학적 해석을 통해 흡착 특성을 정량적으로 조사하였다. 등온 흡착 평형 데이터를 Langmuir, Freundlich, Temkin, Dubinin-Radushkevich 등온흡착식에 적용한 결과, Temkin 식이 가장 높은 r2 (>0.981)와 가장 낮은 RMSD (<2.009)로 실로퓨트를 이용한 세팔로마닌 흡착에 가장 적합함을 알 수 있었다. 흡착 온도가 증가함에 따라 흡착 용량은 오히려 감소하였다. 속도론적 해석에서 유사 이차 반응 속도식에 적합함을 확인하였으며, 입자 내 확산 모델로부터 흡착이 매우 빠르게 평형에 도달하였으며 입자 내 확산과 필름 확산이 율속 단계에 거의 영향을 미치지 않았다. 열역학적 파라미터인 DH0 (-9.155 kJ/mol)와 DS0 (-0.017 kJ/mol·K)는 모두 음의 값을 보여 흡착 공정이 발열이며 흡착인 진행됨에 따라 무질서도가 감소하였다. DG0 (-4.861~-4.522 kJ/mol)는음의 값을 보여 흡착 공정이 자발적으로 수행됨을 알 수 있었다. 또한 Ea (6.353 kJ/ mol)를 통해 물리적 흡착임을 확인하였다. 주어진 qe (100~135 mg/g)에 서 DHX값은 -20 kJ/mol 정도로 거의 일정하였다. 따라서 흡착제-흡착질 상호작용이 존재하지 않아 흡착제의 표면에너지가 균일한 물리적 흡착임을 재확인하였다.
감사
이 논문은 2018년도 정부(교육부)의 재원으로 한국연구재단의 지원을 받아 수행된 기초연구사업의 결과물입니다 (과제번호 : 2018r1da3a03000683).
참고문헌
- Dong, H., Luo, L., Zhou, S., Bi, P., Liu, K. and Zhao, J., "Separation and Purification of Taxol and Cephalomannine from Taxus cuspidada by Normal Phase Chromatography and Twice-reversed-phase Chromatography," Chem. Nat. Compd., 43, 478-480(2007). https://doi.org/10.1007/s10600-007-0169-z
- Zwawiak, J. and Zaprutko, L., "A Brief History of Taxol", J. Med. Sci., 83, 47-52(2014). https://doi.org/10.20883/medical.e43
- Georgiev, M. I., Weber, J. and Maciuk, A., "Bioprocessing of Plant Cell Cultures for Mass Production of Targeted Compounds," Appl. Microbiol. Biotechnol., 83, 809-823(2009). https://doi.org/10.1007/s00253-009-2049-x
- McPartland, J. T., Patil, R. A., Malone, M. F. and Roberts, S. C., "Liquid-liquid Extraction for Recovery of Paclitaxel from Plant Cell Culture: Solvent Evaluation and Use of Extractants for Partitioning and Selectivity," Biotechnol. Prog., 28, 990-997(2012). https://doi.org/10.1002/btpr.1562
- Beckvermit, J. T., Anziano, D. J. and Murray, C. K., "An Improved Method for Separating Paclitaxel and Cephalomannine Using Ozone and Girard Reagents," J. Org. Chem., 61, 9038-9040(1996). https://doi.org/10.1021/jo961402l
- Lee, C. G. and Kim, J. H., "Separation Behavior of Paclitaxel and its Semi-synthetic Precursor 10-deacetylpaclitaxel from Plant Cell Cultures," Korean Chem. Eng. Res., 54, 89-93(2016). https://doi.org/10.9713/kcer.2016.54.1.89
- Lee, C. G. and Kim, J. H., "Optimization of Adsorbent Treatment Process for the Purification of Paclitaxel from Plant Cell Cultures of Taxus chinensis," Korean Chem. Eng. Res., 52, 497-502(2014). https://doi.org/10.9713/kcer.2014.52.4.497
- Lee, S. H. and Kim, J. H., "Kinetic and Thermodynamic Characteristics of Microwave-assisted Extraction for the Recovery of Paclitaxel from Taxus chinensis," Process Biochem., 76, 187-193 (2019). https://doi.org/10.1016/j.procbio.2018.11.010
- Lee, C. G. and Kim, J. H., "Increasing the Efficiency of the Separation and Purification Process for Paclitaxel by Pre-treatment with Water," Process Biochem., 51, 1738-1743(2016). https://doi.org/10.1016/j.procbio.2016.07.025
- Lim, Y. S. and Kim, J. H., "Isotherm, Kinetic, and Thermodynamic Studies on the Adsorption of 13-dehydroxybaccatin III from Taxus chinensis onto Sylopute," J. Chem. Thermodyn., 115, 261-268(2017). https://doi.org/10.1016/j.jct.2017.08.009
- Bang, S. Y. and Kim, J. H., "Isotherm, Kinetic, and Thermodynamic Studies on the Adsorption Behavior of 10-deacetylpacli-Taxel Onto Sylopute," Biotechnol. Bioproc. Eng., 22, 620-630(2017). https://doi.org/10.1007/s12257-017-0247-4
- Shin, H. S. and Kim, J. H., "Isotherm, Kinetic and Thermodynamic Characteristics of Adsorption of Paclitaxel onto Daion HP-20," Process Biochem., 51, 917-924(2016). https://doi.org/10.1016/j.procbio.2016.03.013
- Park, S. H. and Kim, J. H., "Isotherm, Kinetic, and Thermodynamic Characteristics for Adsorption of 2,5-xylenol Onto Activated Carbon," Biotechnol. Bioproc. Eng., 23, 541-549(2018). https://doi.org/10.1007/s12257-018-0259-8
- Kim, Y. S. and Kim, J. H., "Isotherm, Kinetic and Thermodynamic Studies on the Adsorption of Paclitaxel Onto Sylopute," J. Chem. Thermodyn., 130, 104-113(2019). https://doi.org/10.1016/j.jct.2018.10.005
- Lee, J. J., "Isotherm, Kinetic and Thermodynamic Characteristics for Adsorption of Congo Red by Activated Carbon," Korean Chem. Eng. Res., 53, 64-70(2015). https://doi.org/10.9713/kcer.2015.53.1.64
- Chowdhury, S., Mishra, R., Saha, P. and Kushwaha, P., "Adsorption Thermodynamics, Kinetics and Isosteric Heat of Adsorption of Malachite Green Onto Chemically Modified Rice Husk," Desalination, 265, 159-168(2011). https://doi.org/10.1016/j.desal.2010.07.047
- Oh, H. J., Jung, K. Y. and Kim, J. H., "Effect of Pore Size of Mesoporous Spherical Silica for the Purification of Paclitaxel From Plant Cell Cultures," KSBB J., 28, 208-212(2013). https://doi.org/10.7841/ksbbj.2013.28.3.208
-
Grecel, O., Ozcan, A., Ozcan, A. S. and Gercel, H. F., "Preparation of Activated Carbon from a Renewable Bio-plant of Euphorbia Rigida by
$H_2SO_4$ Activation and Its Adsorption Behavior in Aqueous Solutions," Appl. Surf. Sci., 253, 4843-4852(2007). https://doi.org/10.1016/j.apsusc.2006.10.053 - Tan, I. A. W., Ahmad, A. L. and Hameed, B. H., "Adsorption of Basic Dye on High-surface-area Activated Carbon Prepared from Coconut Husk: Equilibrium, Kinetic and Thermodynamic Studies," J. Hazard. Mater., 154, 337-346(2008). https://doi.org/10.1016/j.jhazmat.2007.10.031
- Hata, H., Saeki, S., Kimura, T., Sugahara, Y. and Kuroda, K., "Adsorption of Taxol Into Ordered Mesoporous Silica with Various Pore Diameters," Chem Mater., 11, 1110-1119(1999). https://doi.org/10.1021/cm981061n
- Boparai, H. K., Joseph, M. and O'Carroll, D. M., "Kinetics and Thermodynamics of Cadmium Ion Removal by Adsorption Onto Nano Zerovalent Iron Particles," J. Hazard. Mater., 186, 458-465 (2011). https://doi.org/10.1016/j.jhazmat.2010.11.029
- Zulfikar, M. A., "Effect of Temperature on Adsorption of Humic Acid from Peat Water Onto Pyrophyllite," Int. J. Chem., Environ. Biol. Sci., 1, 88-90(2013).
- Celekli, A., Ilgun, G. and Bozkurt, H., "Sorption Equilibrium, Kinetic, Thermodynamic, and Desorption Studies of Reactive Red 120 on Chara contraria," Chem. Eng. J., 191, 228-235(2012). https://doi.org/10.1016/j.cej.2012.03.007
- Saha, P. and Chowdhury, S., "Insight Into Adsorption Thermodynamics," Prof. Mizutani Tadashi (Ed.), ISBN: 978-953-307-544-0, InTech, Available from: http://www.intechopen.com/books/thermodynamics/insight-into-adsorption-thermodynamics (2011).