International Journal of Plant Research

p-ISSN: 2163-2596    e-ISSN: 2163-260X

2014;  4(2): 43-50

doi:10.5923/j.plant.20140402.01

Optimization of the Polyphenolics Extraction from Chamomile Ligulate Flowers Using Response Surface Methodology

Zoran Zeković , Aleksandra Cvetanović , Branimir Pavlić , Jaroslava Švarc-Gajić , Marija Radojković

Department of Biotechnology and Pharmaceutical Engineering, Faculty of Technology, University of Novi Sad, Bulevar Cara Lazara 1,21000 Novi Sad, Serbia

Correspondence to: Aleksandra Cvetanović , Department of Biotechnology and Pharmaceutical Engineering, Faculty of Technology, University of Novi Sad, Bulevar Cara Lazara 1,21000 Novi Sad, Serbia.

Email:

Copyright © 2014 Scientific & Academic Publishing. All Rights Reserved.

Abstract

Different extraction conditions may significantly influence extraction of particular compound groups. In present study response surface methodology (RSM) based on Box-Behnken design was used to define the best combination of extraction temperature (40-80°C), ethanol concentration (50-90%) and extraction time (30-90 min) for maximum yield of antioxidant compounds and maximum antioxidant activity of chamomile ligulate flowers (CLF) extracts. Experimental values of total extraction yields were in the range from 9.36 to 32.02 %. Content of phenolic compounds varied in the range between 35.55 and 57.35 mg/g, whereas the flavonoids were detected in the range from 14.06 to 26.31 mg/g. Antioxidant activity expressed as the 50% inhibition concentration (IC50 value) was in the range from 0.014 to 0.278 mg/mL. The experimental results were fitted to a second order quadratic polynomial model and they have shown a good fit to the proposed model (R2>0.90). Determined optimized conditions for maximizing yield of antioxidant compounds were withing the experimental range. Validation of the model confirmed that experimental values agreed with those predicted, thus indicating suitability of the used model and the RSM approach in optimizing the extraction conditions.

Keywords: Chamomila, Antioxidant activity, Phenols, Flavonoids, Response surface methodology (RSM)

Cite this paper: Zoran Zeković , Aleksandra Cvetanović , Branimir Pavlić , Jaroslava Švarc-Gajić , Marija Radojković , Optimization of the Polyphenolics Extraction from Chamomile Ligulate Flowers Using Response Surface Methodology, International Journal of Plant Research, Vol. 4 No. 2, 2014, pp. 43-50. doi: 10.5923/j.plant.20140402.01.

1. Introduction

In last decades, focus on medicinal plant research increased all over the world. In traditional medicine, plant formulations and combined extracts of plants are used for the treatment of a wide variety of diseases [1]. Plants have been identified as source of various phytochemicals, which possess a range of benefits to human health. Many of these plants have been tested for their antioxidant activities, and results have shown that some crude extracts or isolated pure compounds from them were more effective antioxidants in vitro than BHT or vitamin E [2-4]. Taking this into account, medicinal plants can be regarded as a potential source of natural antioxidants [5]. Their extraction from plant material is very important step in the manufacture of products enriched with phytochemicals. The knowledge on factors influencing the process conditions is necessary to enhance the extraction efficiency for any bioactive compound. Many factors such as solvent composition, time of extraction, temperature, pH, particle size, etc., may significantly influence the solid-liquid extraction [6-13]. Response Surface Methodology (RSM) is a collection of statistical and mathematical methods useful for developing, improving, and optimizing processes [14]. It is effective for optimization of complex processes because it allows efficient and easier interpretation of experiments. Several authors already employed RSM for the optimization of extraction process in order to maximize yield of various polyphenolic compounds from various sources [15-17]. Among other RSM is being used for optimization the extraction of phenolic compounds from black and white mulberry [18], from the extract of hawthorn (Crataeguslaevigata) [19], Pyracantha fortuneana fruit [20]. RSM proved effective due to its ability to analyze effects of independent variables and their mutual interactions on investigated responses.
Matricariarecutita L., commonly known as German chamomile, has a long history of application in herbal medicine. For thousands of years extracts of this plant have been used by folk healers to treat wounds, ulcers, skin irritations, neuralgia, rheumatic pain, gastrointestinal upset, and a vast number of other ailments [21-24]. The plant is used in different commercially available forms such as tea, infusion, and dietary supplements. Several studies have been carried out in order to evaluate the antioxidative potential of this plant. Essential oil, ethanolic extracts, and infusions of chamomile have been tested, and all of them exhibited antioxidative activity [25]. This antioxidant benefit of chamomile is associated with several groups of active components. Among other, polyphenols from this plant possess strong antioxidant activity.
The objective of the present study was to investigate the effects of the extraction temperature, ethanol concentration and extraction time on antioxidative capacity of extracts prepared from chamomile ligulate flowers (CLF), and to optimize the extraction parameters with a consideration of three responses by applying RSM. Additionally, verification of calculated optimal conditions was done.

2. Material and Methods

2.1. Chemicals and Reagents

1,1-Diphenyl-2-picryl-hydrazyl-hydrate(DPPH), Folin – Ciocalteu reagent, chlorogenic acid and rutin were purchased from Sigma-Aldrich (St. Louis, Missouri, USA). Aluminium chloride hexahydrate, sodium carbonate, and sodium acetate trihydrate were produced from Merck (Darmstadt, Germany). All other chemicals and reagents were analytical reagent grade.

2.2. Sample Preparation

In this study white ligulate flowers of chamomile (CLF) were used. Flowers were produced by the Institute of Field and Vegetable Crops, Bački Petrovac, Serbia. Plant material was collected in spring 2011. Extraction process was done with dry plant material. Drying was carried out in solar dryer at temperature of 40°C and was finished when the moisture of the plant material was approximately 12%. Ligulate flowers were separated from other parts by sifting through sieves and, after that, CLF were packed in paper bags and stored in the dark at room temperature until use.
Plant samples were extracted at different temperature (40, 60 or 80°C), using ethanol as a solvent in a different concentration (50, 70 or 90%) and during different extraction time (30, 60 or 90 min). Sample-solvent ratio was 1:50 (w/v) for all experiments. Obtained liquid extracts were stored at 4°C until analysis.

2.3. Determination of Total Extraction Yield

In order to determine the the total extraction yield, certain volume of liquid extracts were evaporated (Devarot, Elektromedicina, Ljubljana, Slovenia) under vacuum. After removing the solvent, drying was performed at 105°C to the constant mass. Further calculating of total extraction yield was done according to procedure described in pharmacopoeia (Ph. Jug. V) [26].

2.4. Determination of Antioxidant Components

Total phenols content was determined spectrophotometrically according to the Folin-Cicolteu colorimetric assay [27, 28]. Chlorogenic acid was used as a standard and the results were expressed as mg of chorogenic acid equivalent (CAE) per g of dry plant material (mg CAE/g). Absorbance was measured at 750 nm, using a JANWAY 6300 VIS-spectrophotometer.
Total flavonoids content was determined by aluminium chloride colorimetric assay [29] using rutin as a standard. The results were expressed as mg of rutin equivalents (RE) per g of dry plant material (mg RE/g). Absorbance was measured at 430 nm, using a JANWAY 6300 VIS-spectrophotometer.

2.5. DPPH Assay

The capacity to scavenge the “stable” free radical 2,2-diphenil-1-picrylhydrazyl (DPPH) was monitored according to the method of Espin [30]. Solution of 90 μM DPPH in methanol (95%) was prepared and 1 ml was added to 3 ml of diluted liquid extract. Absorbance was measured 60 minutes later at 515 nm. Blank was prepared using water instead of CLF liquid extract. Radical scavenging capacity (%RSC) was calculated by the following Eq. (1):
(1)
where A sample is absorbance of the sample solution and Ablank is absorbance of the blank sample. This activity was also expressed as the 50% inhibition concentration (IC50), i.e. the concentration of the test solution required to scavenge 50% of the initial radical.

2.6. Experimental Design

RSM is a useful technique for developing, improving and optimizing processes. The most extensive application of RSM is in the industry, particularly when several input variables influence some performance measure or quality characteristic of the product or process. This performance measure or quality characteristic is called the response [31].
After determining the preliminary ranges of independent variables, Box-Behnken experimental design was adopted in this optimization study. The extraction temperature, ethanol concentration, and extraction time were applied as independent variables. The range of those variables and their levels are presented in Table 1.
Table 1. Uncoded and Coded Levels of Independent Variables Used in the RSM Design
     
The behavior of the system can be described by the following quadratic equation:
(2)
where Y represents the measured response; β0 is constant, bj, bjj, bij are the linear, quadratic and interactive coefficients of the model, respectively; Xi and Xj are the levels of the independent variables.

3. Results and Discussion

3.1. Response Surface Models

During this study, the target compounds of the extraction process were total phenols (TP) and total flavonoids (TF), and the targeted activity was antioxidant capacity, i.e. IC50 value. Moreover, total extraction yield (EY) was also observed. In order to explore the functional relationship between inputs (extraction temperature, solvent concentration and extraction time) and outputs (extraction yield, target compounds and target activity) RSM was applied. Obtained experimental results of EY, TP and TF content, as well as IC50 values, are shown in Table 2.
Table 2. Experimental Matrix and Values of the Observed Responses
     
The results of EY, TP, TF and IC50 obtained in the experiments were 9.36-32.02%, 35.55-57.35 mg CAE/g, 14.06-26.31 mg RE/g, and 0.014-0.278 mg/ml, respectively. The lowest values of EY, TP and TF were obtained at the same point (temperature of 40°C, 50% of ethanol concentration and 60 min of extraction time).
Correlation between antioxidant activity and polyphenolic content (TP and TF) was confirmed with the highest value of IC50 at the same point where TP and TF were the lowest. The highest EY (32.02%) and TP (57.35 mg CAE/g) were obtained under the experimental conditions of X1 = 80°C, X2 = 50% and X3 = 60 min, whereas the highest TF (26.31 mg RE/g) was obtained under conditions of X1 = 80°C, X2= 70% and X3 = 30 min. The highest antioxidant activity, e.i. the lowest IC50 value was observed for the extract obtained with 70% ethanol at 60°C, after 60 min of extraction.
Experimental results from Table 2 were processed with multiple linear regression using the second-order polynomial model (Eq. (2)). The regression coefficients of the intercept, linear, cross product and quadratic terms are presented in the Table 3. As Joglejar and May suggested [32] good fit of the model (R2) should be at least 0.80, which indicates the adequacy of the applied regression model. The values of R2 and CV for EY, TP and TF were from 0.912 to 0.976 and from 3.075 to 8.404%, respectively. Therefore, it was suggested that quadratic model fitted well with the experimental data.
Table 3. Regression Equation Coefficients for the Selected Responses
     
Suitability of the model was also analyzed through the ANOVA for the model. Calculated statistical parameters are presented in Table 4. According to the p-values of the F-test for suggested model (p<0.05) and the lack of fit (p>0.05), it was suggested that the model is suitable for the investigated extraction system and can relatively predict its variations. Model equations for relationship between EY, TP and TF and independent variables were obtained by applying multiple regression analysis (Table 4). By applying these equations it is possible to predict values of each response. Antioxidant activity depending on independent variables could not be successfully modeled by Eq.2 due to poor R2 (0.743) and high CV (21.44%). F-test applied on the model further confirmed that the quadratic model was not satisfactory for description of this system (p>0.05). Neither linear nor cubic model were suitable in this case. Therefore, relationship between IC50 value (antioxidant activity) and extraction parameters (temperature, ethanol concentration and extraction time) was defined basing on analysis.
Table 4. Analysis of Variance
     

3.2. Effects of Temperature, Solvent Concentration and Extraction Time on Investigated Response Parameters

In case model equation for EY, linear terms of temperature and ethanol concentration have significant influence (p<0.01). Temperature exhibited positive influence which could be explained with heat influence on softening of the plant tissues and increasing diffusivity of the extraction solvent into cells. Elevated temperature also reduces solvent viscosity and surface tension and hence to promote the extraction of soluble compounds [33]. Negative effect of ethanol concentration was also expected due to decrease in polarity for higher ethanol concentrations. Lower ethanol concentration contribute to higher EY, however these were not in direct correlation with antioxidant activity of obtained extracts.
Since most of the compounds in the plant material are hydrophilic, extracts obtained with lower ethanol concentrations will have high content of concomitant compounds which does not possess high antioxidant capacity. Interaction between temperature and ethanol concentration showed significant influence on the EY. Positive effect could be explained with the moderate polarity and modified physical properties (viscosity and surface tension) of the extraction solvent at elevated temperatures. These changes increase its capacity for extraction. Quadratic term of temperature exhibited negative effects which suggest that EY will reach the saddle point at certain temperature which the yield starts to decrease. Lower EY at very high temperature is probably due to degradation of extractable compounds and possible formation of aggregates.
Efficiency of extraction process is influenced by multiple parameters. The best way of expressing the effects of different extraction parameters on polyphenolic compounds content and antioxidant activity was to generate response surfaces of the model [18]. The 3D response surface plots are the graphical representations of regression equations. They provide a method to visualize the relationship between responses and experimental levels of each variable and the type of interactions between two test variables. The shapes of the plots indicate whether the mutual interactions between the variables are significant or not [34]. Figure 1A shows the 3D surface plots of the interactive effects of the independent variables corresponding to the extraction yields (EY). It can be seen that increase in solvent concentration leads to reduction of EY. On the other hand EY increasing with temperature up to about 75°C. Further increase of temperature leads to reduction of EY. The same situation can be observed with the extraction time. The yields increase with time, but after 65 minutes of extraction, amount of extracted matters decreases.
In general terms, the extraction efficiency of phenolic compounds is a function of several process variables. Many authors reported the influence of different variables on the extraction of phenols. Among others, the most important factors influencing the recovery of phenols from plant material are solvent type, temperature and contact time. The positive or negative effects of each variable on the mass transfer phenomenon is specific and it is not always obvious [35]. In presented research and derived modellinear and cross product term of temperature (p<0.01) had significant influence on phenol extraction.
However, cross product term of ethanol concentration showed negative effect. On the other hand linear term and quadratic term of extraction time showed significant influence (0.01<p< 0.05). Influence of parameters on TP content was visualized by the chart (Figure 1B). From the chart it is obvious that increase in ethanol concentration up to approximately 60% leads to higher yields for TP, while further increase of solvent concentration causes the decrease of the output value. Considering the fact that the extraction efficiency of the phenolic compounds decreases when using 90% ethanol, due to the number of hydroxyl groups in phenols that are rather hydrophilic, and as such generally more soluble in water-ethanol solutions, those results were expected. In terms of temperature influence, figure 1B shows that with extraction temperature the yield of TP increased. This result is supported with a claim that temperature generally shows positive effects on extraction of antioxidant compounds. Also, a linear increase of TP with extraction time could be observed.
All linear terms of process variables had highly significant influence on TF (p<0.01). Temperature effect was once again positive, while other two variables had negative effects. Elevated temperature modifies polarity of the solvent and increases the solubility of flavonoid compounds in the extraction solvent [36]. Although, linear effects of ethanol concentration and extraction time were significant, they were still minor comparing to the temperature. These effects suggested that less polar solvent (90% ethanol) would not be suitable as 50% or 70% ethanol. Prolonged extraction time could lead to degradation of some flavonoid compounds due to prolonged exposure to environment factors such as light and oxygen [33, 37]. Reduction in extraction time could be significant for technological aspect of the process because it can result in time and cost saving [38]. Temperature-ethanol concentration and temperature-time interaction showed significant cross product effects. Positive effect of temperature-ethanol concentration represented interactive influence between variables where temperature improves solvent properties and increases affinity to desirable group of compounds. However, temperature-extraction time interaction suggests that elevated temperature could produce undesirable effects which probably come from chemical modifications of compounds present in crude extracts. Quadratic terms of all three variables had significant influence. The negative quadratic effect of extraction time confirmed the deceleration in the extraction yield, as Fick’s second law of diffusion predicts a final equilibrium between the solute concentrations in the solid matrix and in the bulk solution after a certain time [15].The relationships between independent variables and TF were illustrated in three-dimensional graphs (Figure 1C). As it can be seen from the figure the amount of TF increased linearly with of temperature. With the increase of ethanol concentration up to approximately 75% increase of the TF was observed, while further increase caused the decrease of the output value. This indicated that there is a saddle point near 75% ethanol concentration.
Figure 1. Response surface plots showing effects of investigated parameters on: A-extraction yields, B-phenolics contents, C-flavonoids contents
The polarity of the solvent plays important role in the selective extraction of different flavonoid families. Furthermore, TF content decreased with the extraction time up to about 50 min, after that total flavonoid amount increased. Influence of independent variables on IC50 value is presented in Figure 2.
Figure 2. Linear Influences of Independent Variables on IC50 Value
Increase of the extraction temperature caused the decrease of the IC50 value (increase of the antioxidant activity). This result is in accordance with previous reports which suggested that extraction yield of thermally stable antioxidants at elevated temperatures was high [16]. Antioxidant activity decreased with the increase of ethanol concentration (Fig 2), indicating that 50% ethanol is a good solvent for most of antioxidants from CLF. On the other hand, extraction time did not show any significant influence on the antioxidant activity. Interaction between independent variables on the antioxidant activity is more complex, but the cross product and quadratic terms could not be determined due to inconsistency between Eq. (2) and experimental results.

3.3. Verification of Experiments

The final step in this research was the verification of suitability of the defined model equation. By using RSM total extraction yield (EY) and antioxidant compounds (TP, TF) were optimized. Optimized system was developed for each of response Results of optimized parameters, predicted and observed values are reported in Table 5. All experiments were done at triplicate.
Table 5. Comparison of Predicted and Experimental Values for the Response Variables
     
The model has shown good prediction for two of three output variables (EY, TP), however, observed values of TF were bellow expected. The good correlation between predicted and observed values confirmed that the response model was suitable for the intended optimization.

4. Conclusions

High yield of extraction process is very important goal, especially for large-scale industrial applications. Considering the great role of antioxidants, there is increasing number of studies oriented towards the production of antioxidant-rich products. Polyphenolic compounds are recognized as one of the most important antioxidants. In order to produce extracts with high level of total phenols, total flavonoids and antioxidant activity, extraction of chamomile ligulate flowers (CLF) was successfully optimized by response surface methodology (RSM). Optimal extraction parameters defined by this approach for all target compounds were validated. In two of three cases, experimental results were close to the predicted calculated, confirming the validity and adequacy of the proposed models. Considering the fact that this is the first report on optimizing extraction technology of antioxidant compounds of CLF, these results can be useful for developing new CLF-based products.

ACKNOWLEDGMENTS

Financial support of this work by the Serbian Ministry of Education, Science and Technological Development, Project No. TR 31013 is greatfully acknowledged. We are greateful to Dr Dušan Adamović, Institute of Field and Vegetable Crops, BačkiPetrovac, Serbia, for his support in supplying plant material.

References

[1]  Ansarullah, R.N.J., Thounaojam, M.C., Patel, V., Devkar, R.V., and Ramachandran, A.V., 2009, Antihyperlipidemic potential of a polyherbal preparation on triton WR 1339 (Tyloxapol) induced hyperlipidemia: a comparasion with lovastatin, International Journal of Green Pharmcy, 3, 119-124.
[2]  Gordon, M.H., and Weng, X. C., 1992, Antioxidant Properties of extracts from Tanshen (Salvia militiorrhiza Bunge), Food Chemistry, 44,119-122.
[3]  Gu, L.W., and Weng, X.C., 2001, Antioxidantactivity and components of Salvia plebeia R.Br.-a Chinese Herb, Food Chemistry, 73, 299-305.
[4]  Pyo, Y.H., Lee, T.C., Logendrac, L., and Rosen, R. T., 2004, Anto-ioxidant activity and phenolic compounds of Swiss Chard (Beta Vulgaris Subspecies Cycla) Extracts, Food Chemistry, 85, 19-26.
[5]  Cesquini, M.,Torsoni, M.A., Stoppa, G.R., and Ogo, S.H., 2003, T-BOOH-Induced oxidative damage in sickle red blood cells and the role of flavonoids, Biomedicine & Pharmacotherapy, 57, 124-129.
[6]  Pinelo, M., Del Fabbro, P., Manzocco, L., Nunez, M.J., and Nicoli, M.C., 2005, Optimization of continuous phenol extraction from Vitisvinifera byproducts, Food Chemistry, 92, 109-117.
[7]  Spigno, G., and De Faveri, D.M., 2007, Antioxidants from grape stalks and marc: influence of extraction procedure on yield, purity and antioxidant power of the extracts, Journal of Food Engineering, 78, 793-801.
[8]  Spigno, G., Tramelli, L, and De Faveri, D.M., 2007, Effects of extraction time, temperature and solvent on concentration and antioxidant activity of grape marc phenolics, Journal of Food Engineering, 81, 200 -208.
[9]  Pompeu, D.R., Silva E.M., and Rogez, H., 2009, Optimisation of the solvent extraction of phenolic antioxidants from fruits of Euterpeoleracea using Response surface methodology, Bioresurce Technology, 100, 6076-6082.
[10]  Wijangaard, H.H., and Brunton, N., 2010, The optimization of solid-liquid extraction of antioxidants from apple pomace by response surface methodology, Journal of Food Engineering, 96, 134-140.
[11]  Cacace, J.E., and Mazza, G., 2003, Mass transfer process during extraction of phenolic compounds from milled berries, Journal of Food Engineering, 59, 379-389.
[12]  Myers, R.H., and Montgomery, D.C., 1995, Response Surface Methodology: Process and Product in Optimization using Designed Experiments. Wiley, New York,
[13]  Švarc-Gajić, J., 2011, Samples & Sample Preparation in Analytical Chemistry, Novapublishers, New York.
[14]  Adamczyka, J., Hornya, N., Tricoteauxa, A., Jouana, P.Y., and Zadamb, M., 2008, On the use of response surface methodology to predict and interpret the preferred c-axis orientation of sputtered AlN thin films, Applied Surface Science, 254, 1744–1750.
[15]  Silva, E. M., Rogez, H., and Larondelle, Y., 2007, Optimization of extraction of phenolics from Inga edulis leaves using response surface methodology, Separation and Purification Technology, 55, 381–387.
[16]  Liyana-Pathirana, C., and Shahidi, F., 2005, Optimization of extraction of phenolic compounds from wheat using response surface methodology, Food Chemistry, 93, 47-56.
[17]  Fan, G., Han, Y., Gu, Z., & Chen, D. (2008). Optimizing conditions for anthocyanins extraction from purple sweet potato using response surface methodology (RSM). LWT-Food Science and Technology, 41(1), 155-160.
[18]  Radojković, M.., Zeković, Z., Sudar, R., Jokić, S., and Cvetanović, A., 2013, Optimization of solid-liquid extraction of antioxidants and saccharides from mulberry fruit by response surface methodology, Journal of Food and Nutrition Research, 52, 146-15.
[19]  Cui, H.Y., Jia, X.Y., Zhang, X., Zhang, J., Zhang, Z.Q., 2011, Optimization of high-speed counter-current chromatography for separation of polyphenols from the extract of hawthorn (Crataeguslaevigata) with response surface methodology, Separation and Purification Technology 77, 269–274.
[20]  Zhao, C.F., Li, S., Li, S.J., Songa, G.H., Yu, L.J., Zhang, H.,2013, Extaction optimization approach to improve accessibility of functional fraction based on combination of total polyphenol, chromatographic profiling and antioxidant activity evaluation: Pyracanthafortuneana fruit as an example, Journal of functional foods, 5, 715-728.
[21]  Ashnagar, A., Naseri, N., and Alavi, S.Y., 2009, Isolation and identification of the major chemical constituents in the capitula of Matricariachamomilla grown in Khuzestan Province of Iran, Asian Journal of Chemistry 21, 4981–4986.
[22]  McKay, D.L., and Blumberg, J.B., 2006, A Review of the bioactivity and potential health benefits of chamomile tea (Matricariarecutita L.), Phytotherapy Research 20, 519–530.
[23]  Srivastava, J.K, and Gupta, S., 2007, Antiproliferative and apoptotic effects of chamomile extract in various human cancer cells, Journal of Agricultural and Food Chemistry 55, 9470–9478.
[24]  Iverson, C.I., Zahid,S., Li, Y.-, Shoqafi, A.H., Ata, A., and Samarasekera R., 2010, Glutathione S-transferase inhibitory, free radical scavenging, and anti-leishmanial activities of chemical constituents of Artocarpusnobilis and Matricaria Chamomilla, Phytochemistry Letters, 3, 207–211.
[25]  Petronilho, S., Maraschin, M., Coimbra, M.A., and Rocha, S.M., 2012, In vitro and in vivo studies of natural products: A challenge for their valuation. The case study of chamomile (Matricaria recutita L.), Industrial Corps and products, 40, 1-12.
[26]  Pharmacopoea Jugoslavica edito qarta, (Ph. Jug. IV), 1984, Federal Office of Public Health, Beograd.
[27]  Singleton, V. L., and Rossi Jr, J. A., 1965, Colorimetry of total phenolics with phosphomolybdic phosphotungstic acid reagents, American Journal of Enology and Viticulture, 16, 144-158.
[28]  Kähkönen, M.P., Hopia, A.I., Vuorela, H.J., Rauha, J.P., Pihlaja, K., Kujala, T.S., and Heinonen, M., 1999, Antioxidant activity of plant extracts containing phenolic compounds, Journal of Agriculture and food chemistry, 47, 3954–3962.
[29]  Markham, K. R., 1989, Flavones, flavonoids, and their glycosides, In: Harborne, J. B. Dey, P. M. (Ed.): Methods in plant biochemistry. Vol. 1: Plant phenolics, Academic Press, London.
[30]  Espin, J. C., Soler-Rivas,C., and Wichers, H. J., 2000, Characterization of the total free radical scavenger capacity of vegetable oils and oil fractions using 2,2-diphenyl-1- picrylhydrazyl radical, Journal of Agricultural Food Chemistry, 48, 648-656.
[31]  Myers, R.H., and Montgomery, D.C., 2011, Anderson-Cook, C.M., Response Surface Methodology: Process and Product Optimization Using Designed Experiments, 3rd Edition, Wiley, New York
[32]  Joglekar, A. M., and May, A. T., 1987, Product excellence through design of experiments, Cereal Foods World, 32, 857-868.
[33]  Juntachote, T., Berghofer, E., Bauer, F., and Siebenhandl S., 2006, The application of response surface methodology to the production of phenolic extracts of lemon grass, galangal, holy basil and rosemary, International Journal of Food Science and Technology, 41,  121–133.
[34]  Zhong, K., and Wang, Q., 2010, Optimization of ultrasonic extraction of polysaccharides from dried longan pulp using response surface methodology, Carbohydrate Polymers, 80, 19-25.
[35]  Angela, M., and Meireler A., 2008, Extracting Bioactive Compounds for Food Products Theory and Application, Edit by, Published by CRC Press.
[36]  Cacace, E., and Mazza, G., 2003, Mass transfer process during extraction of phenolic compounds from milled berries, Journal of Food Engineering, 59, 379-389.
[37]  Chirinos, R., Rogez, H., Campos, D., Pedreschi, R., and Larondelle, Y., 2007, Optimization of extraction conditions of antioxidant phenolic compounds from mashua (Tropaeolumtuberosum Ruíz & Pavón) tubers, Separation and Purification Technology, 55, 217-225.
[38]  Chew, K.K, Khoo, M.Z., Ng, S. Y., Thoo, Y. Y., Wan Aida, W. M., and Ho, C. W., 2011, Effect of ethanol concentration, extraction time and extraction temperature on the recovery of phenolic compounds and antioxidant capacity of Orthosiphonstamineus extracts, International Journal of Food Research, 18, 1427-1435.