Next Article in Journal
Correction: Zhang, Y., et al. Clones of FeSOD, MDHAR, DHAR Genes from White Clover and Gene Expression Analysis of ROS-Scavenging Enzymes during Abiotic Stress and Hormone Treatments. Molecules 2015, 20, 20939–20954
Previous Article in Journal
Targeted and Untargeted Metabolomics to Explore the Bioavailability of the Secoiridoids from a Seed/Fruit Extract (Fraxinus angustifolia Vahl) in Human Healthy Volunteers: A Preliminary Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Optimization of Ultrasonic-Assisted Enzymatic Extraction of Atratylodes macrocephala Polysaccharides and Antioxidants Using Response Surface Methodology and Desirability Function Approach

1
Development and Research Center of Official Silkworm Resources, Zhejiang Academy of Traditional Chinese Medicine, Hangzhou 310007, China
2
Key Laboratory of Research and Development of Chinese Medicine of Zhejiang Province, Hangzhou 310007, China
3
College of Pharmacy, Hunan Chinese Medical University, Changsha 410208, China
4
College of Life Science, Zhejiang Chinese Medical University, Hangzhou 310015, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2015, 20(12), 22220-22235; https://doi.org/10.3390/molecules201219837
Submission received: 14 October 2015 / Revised: 30 November 2015 / Accepted: 4 December 2015 / Published: 11 December 2015
(This article belongs to the Section Medicinal Chemistry)

Abstract

:
Rhizoma Atractylodes macrocephala polysaccharides (RAMP) have been reported to have a variety of important biological activities. In this study, an ultrasonic-assisted enzymatic extraction (UAEE) was employed to obtain the highest extraction yield and strongest antioxidant activity of RAMP and optimized by a multi-response optimization process. A three-level four-factor Box-Behnken design (BBD) was performed as response surface methodology (RSM) with desirability function (DF) to attain the optimal extraction parameters. The DPPH scavenging percentage was used to represent the antioxidant ability of RAMP. The maximum D value (0.328), along with the maximum yield (59.92%) and DPPH scavenging percentage (13.28%) were achieved at 90.54 min, 57.99 °C, 1.95% cellulase and 225.29 W. These values were further validated and found to be in good agreement with the predicted values. Compared to the other extraction methods, both the yield and scavenging percentage of RAMP obtained by UAEE was favorable and the method appeared to be time-saving and of high efficiency. These results demostrated that UAEE is an appropriate and effective extraction technique. Moreover, RSM with DF approach has been proved to be adequate for the design and optimization of the extraction parameters for RAMP. This work has a wide range of implications for the design and operation of polysaccharide extraction processes.

1. Introduction

Rhizoma Atractylodes macrocephala (RAM) is the dried rhizome of Atractylodes macrocephala Koidz, which belong to the family Asteraceae (Compositae), found mainly distributed in China, Japan, and Korea [1]. RAM is widely used for the treatment of splenic asthenia, anorexia, oedema, excessive perspiration and abnormal fetal movement [2]. The rhizomes of Atractylodes macrocephala are rich in sesquiterpenes, polyacetylenes, phenylpropanoids, coumarins, glycosides and polysaccharides [1,3,4,5,6]. Recently, interest in Rhizoma Atractylodes macrocephala polysaccharides (RAMP) has been growing, because of the various bioactivities of RAMP, including immune response stimulating [7,8,9], aging defying [10], neuroprotective effects [11], hyperglycemic activity [12] and protective effects against liver ischemia reperfusion injury [13]. Most of the aforementioned activities are related with antioxidant properties and polysaccharides are widely considered to be a strong antioxidants.
However, little attention has been paid to the effective extraction of RAMP. Conventionally, polysaccharides are extracted by refluxing in hot water, which is often time-consuming, gives low yields of polysaccharides and even results in the loss of some of the pharmacological activity [14]. In recent years, various methods have been developed for the extraction of polysaccharides, such as ultrasonic-assisted extraction (UAE) and enzyme-assisted extraction (EAE). Among them, UAE is a rapid, energy efficient method that gives high extraction yields and has a minimal impact on the bioactivities [15]. Also, with the additional of specific enzymes such as cellulase and proteases, EAE can promote the release of intracellular contents by breaking the cell wall and lipid bodies [16]. Thus, it is considered as a mild, efficient and environmentally friendly extraction method, which has been used recently in the extraction of various kinds of compounds [17]. In this study, to achieve a high yield extraction and maximize the antioxidant ability of RAMP we combined the advantages of the UAE and EAE methods, in what is defined as ultrasonic-assisted enzymatic extraction (UAEE).
Response surface methodology (RSM) is an effective tool for optimizing experimental processes when many factors and interactions may affect the response variables [18]. Most of the work on RSM has been focused on the case where there is only one response of interest [15,17]. As a matter of fact, in many situations there are several response variables of interest or the relationship between the response variables and design variables is too complex to be efficiently estimated using traditional surface fitting approaches [19]. In this case, determination of the optimum settings on design variables would require simultaneous consideration of all the response [20], which is called a multi-response optimization (MRO).
Many creative methods have been developed for MRO, and the general purpose of these methods is to convert the MRO problem into a single aggregated objective function and then construct an efficient algorithm to find the optimum solution [21]. Presently, desirability function (DF) is a popular and established technique for the simultaneous determination of optimum settings of variables for multiple responses. It was first developed by Harrington and modified by Derringer and Suich [22,23]. It has been successfully used to optimize several multi-response cases, and its functions and advantages has been discussed in different articles [24,25,26].
However, the feasibility of using UAEE for RAMP and the multi-response optimization of the UAEE procedure by RSM mated with DF have not yet been explored in the literature. There is no research to emphasize the wide range of approaches that can be employed and the great impact on the polysaccharide extraction industries. Hence, this work aim to apply the UAEE method for the high yield and antioxidant activities of RAMP, and to find the optimal extraction condition using the RSM with DF approach in order to ensure rational utilization.

2. Results and Discussion

2.1. Single Factor Experimental Analysis of UAEE

All the factors that may affect the process must be carefully determined and examined. The experimental domain must be defined for each factor and also a control and measurement method must be established [27]. Thus, it is necessary to carry out single factor experiments to determine the experimental variables and interactions that have a significant influence on one or several responses. In this study, according to the priminary experiments, four major influencing factors, including extraction time, extraction temperature, cellulase concentration and ultrasonic power were selected for the next experiments.

2.1.1. Effect of Extraction Time on the Yield of RAMP

Generally speaking, the longer the extraction time, the higher the yield of polysaccharides [28]. The effect of extraction time on extraction yield of RAMP is shown in Figure 1a. The yield of RAMP continued to increase gradually over the extraction time range of 10–70 min, until it reached its highest point at 70 min, which indicated that the polysaccharides were fully extracted as time goes by. When the extraction time reached 90 min, the increase in the yield of RAMP eased steadily. This phenomenon can be attributed to polysaccharide hydrolysis during the extended extraction [17].
Figure 1. Effects of extraction time (a) extraction temperature (b) cellulase concentration (c) and ultrasonic power (d) on the extraction yield of RAMP (%).
Figure 1. Effects of extraction time (a) extraction temperature (b) cellulase concentration (c) and ultrasonic power (d) on the extraction yield of RAMP (%).
Molecules 20 19837 g001

2.1.2. Effect of Extraction Temperature on the Yield of RAMP

The extraction temperature is an important factor influencing the extraction yield. As shown in Figure 1b, the yield of RAMP was significantly increased when the temperature was increased from 30 to 60 °C, and reached a peak value of 8.07%. Then further increase of extraction temperature resulted in a RAMP.extraction yield decease. One possible explanation is that the polysaccharides may be hydrolyzed at a high temperature [29], and the another reason is the activity of enzyme added to the solution were heavily influenced by the temperature, as enzyme activity is well known to be reduced when the temperature is too low or too high [30].

2.1.3. Effect of Cellulase Concentration on the Yield of RAMP

Much research has shown that different cellulase concentrations will significantly affect the yield of polysaccharides [31,32]. As shown in Figure 1c the cellulase concentrationcurve indicated that the first increase of enzyme concentration lead to an obvious RAMP extraction yield increase till 2.0%, and then began to decrease slightly. This indicated that a 2.0% amount of enzyme was sufficient to obtain a high polysaccharide yield. Taking cellulase consumption and the yield of RAMP into account, the enzyme concentration of 2.0% was high enough for the UAEE.

2.1.4. Effect of Ultrasonic Power on the Yield of RAMP

Ultrasonic power is a critical parameter in the UAEE method. Its effect on the extraction yield of RAMP is shown in Figure 1d. The yield increased as ultrasonic power increased from 100 W to 250 W, and decined when it went up to 300 W. The maximum extraction yield of RAMP was observed when the ultrasonic power was 250 W. A possible explaination is that the polysaccharides’ diffusion coefficient and the solubility of the polysaccharides in the extracting solvent increased with the enhanced ultrasonic power, and this eventually caused an incease in the amount of polysaccharide going out from the rhizome of Atractylodes macrocephala into solution [33]. After reaching the peak, the polysaccharide yield droped because of the loss during ultrasonicated post-processing [29]. To summarize the single-factor studies, the following conditions could used for the response surface methodology experiments: an extraction time of 10–90 min, an extraction temperature of 30–70 °C, a cellulase concentration of 0.5%–2.5% and a ultrasound power of 100–300 W.

2.2. Multi-Response Design and Analysis

2.2.1. Statistical Analysis and Model Fitting

The yield of RAMP in the BBD experiments

As seen in Table 1, the results showed that the yield of polysaccharide (Y1) ranged from 7.22% to 13.10%.
Table 1. BBD matrix and response values for extraction yield and DPPH scavenging percentage of RAMP as well as D value.
Table 1. BBD matrix and response values for extraction yield and DPPH scavenging percentage of RAMP as well as D value.
RunX1X2X3X4Y1 a (%)Y2 a (%)D Value
110500.52008.60 ± 0.6339.11 ± 3.070.167
250701.510010.37 ± 0.5250.51 ± 2.420.243
350501.520012.98 ± 0.6463.13 ± 2.560.337
450501.520012.93 ± 0.6563.26 ± 3.160.336
590501.530012.18 ± 0.6254.16 ± 3.790.287
650501.520013.10 ± 0.5963.87 ± 3.270.342
750701.530012.54 ± 0.6057.74 ± 2.460.308
850700.520010.37 ± 0.5746.60 ± 2.030.227
910501.530011.12 ± 0.7048.70 ± 1.940.247
1090701.520012.15 ± 0.5056.10 ± 2.580.295
1110701.520010.96 ± 0.5648.23 ± 2.950.243
1290301.52009.62 ± 0.7146.13 ± 2.690.213
1350301.530010.35 ± 0.6351.13 ± 2.710.245
1450301.51007.54 ± 0.5542.29 ± 2.100.162
1550501.520012.99 ± 0.7263.88 ± 2.080.340
1650501.520013.06 ± 0.7063.51 ± 2.040.339
1710502.520010.84 ± 0.5949.07 ± 3.000.244
1810501.51008.21 ± 0.7441.75 ± 2.990.172
1950302.520010.67 ± 0.5849.64 ± 3.200.244
2050502.530012.75 ± 0.5055.66 ± 2.710.302
2190502.520012.91 ± 0.7255.63 ± 3.890.305
2250300.52007.22 ± 0.5342.47 ± 3.260.157
2350500.51007.45 ± 0.5039.72 ± 2.230.151
2450502.510010.83 ± 0.5151.88 ± 2.910.256
2590501.51009.95 ± 0.7048.11 ± 2.580.226
2650500.530010.43 ± 0.4550.52 ± 3.230.244
2790500.52009.11 ± 0.6745.91 ± 3.400.203
2850702.520012.51 ± 0.5059.84 ± 2.660.316
2910301.52008.29 ± 0.8043.83 ± 2.650.181
a each value is the mean of triplicate measurements.
Multiple regression analysis was performed to build a mathematical model to find the optimum conditions that maximize the extraction yield of RAMP and study the relationship of the response variable and test variable. The second-order polynomial equation was given by the following expression:
Y1 = 13.01 + 0.66X1 + 1.27X2 + 1.44X3 + 1.25X4 − 0.035X1X2 + 0.39X1X3 − 0.17X1X4 − 0.33X2X3 − 0.16X2X4 − 0.27X3X4 − 1.30X12 − 1.47X22 − 1.34X32 − 1.33X42
The analysis of variance (ANOVA), goodness-of-fit and the adequacy of the regression model were summarized in Table 2. The high model F-value and the low P-value indicates the level of confidence of the selected model [20]. The model F-value of 1337.56 and the associated lower p-value (p < 0.0001) implied the model was highly statistically significant. There was only a 0.01% chance that a “Model F-value” this large could occur due to noise. The value of R2 reflects the proportion of variation in the response attributed to the model rather than to random error [34]. The model has shown a good fit with the high R2 value and adjusted determination coefficient (R2adj) of 0.9993 and 0.9985, respectively. This confirmed that only 0.07% of the total variation is not explained and the experimental data were well fitted by the model. What’s more, the coefficient of the variation (C.V.) refer to the ratio of the standard error (SD) of estimated data to the mean value of the observed response [20]. In this model, the C.V. of 0.68% indicates that the simulation can be considered as resonably reproducible, which also means the precision and the experimental values of the model were of reliability. Table 2 also reveals that the interaction of extraction time and extraction temperature was insignificant (p > 0.05), while all other model terms were highly significant (p < 0.01). The results indicated that the linear coefficients (X1, X2, X3, X4), quadratic term coefficients (X12, X22, X32, X42) and cross product coefficients (X1·X3, X1·X4, X2·X3, X2·X4, X3·X4) were all significantly correlated with the RAMP extraction yield (p < 0.01).
Table 2. ANOVA for dependent variable: the yield of RAMP.
Table 2. ANOVA for dependent variable: the yield of RAMP.
SourceSum of SquaresdfMean SquareFSig. (Prob > F)
Corrected model100.64147.191337.56<0.0001
X15.2015.20967.75<0.0001
X219.28119.283587.29<0.0001
X325.03125.034656.98<0.0001
X418.80118.803498.22<0.0001
X1X2<0.00011<0.00010.910.3558
X1X30.6110.61113.21<0.0001
X1X40.1210.1221.510.0004
X2X30.4310.4379.83<0.0001
X2X40.1010.1019.050.0006
X3X40.2810.2852.27<0.0001
X1211.03111.032052.10<0.0001
X2214.07114.072617.62<0.0001
X3211.56111.562151.64<0.0001
X4211.50111.502139.57<0.0001
Residual0.07514<0.0001
Lack of fit0.05710<0.00011.250.4493
Pure error0.0184<0.0001
Cor total100.7128
R20.9993SD0.60
R2adj0.9985C.V.%0.68
Adeq precision110.900
Response surfaces were plotted by the Design Expert software to explain the interactions of the variables for the maximum response. The corresponding three-dimensional response surfaces are shown in Figure 2. Each figure shows the effects of two factors at a time on the polysaccharide yield while all other factors were kept at zero level. Figure 2a illustrates the polysaccharide yield in response to extraction time and extraction temperature at a fixed ultrasound power of 200 W and a fixed celluose cncentration of 1.5%. The yield was very low at low extraction time and temperature, and increased as extraction time and extraction temperature increased until a peak value was reached. Further increasing the extraction time and temperature led to a decreased polysaccharide yield. Similar trends were observed for the effects of extraction time and cellulase concentration (Figure 2b), of extraction time and ultrasonic power (Figure 2c), of cellulase concentration and extraction temperature (Figure 2d), of ultrasonic power and extraction temperature (Figure 2e) and of ultrasonic power and cellulase concentration (Figure 2f). According to the model, the maximum yield of RAMP was 13.93% and the corresponding variables values were obtained after an extraction time of 61.91 min, at an extraction temperature of 56.99 °C, using a cellulase concentration of 2.00% and 237.96 W of ultrasonic power.
Figure 2. Response surface (3D) showing the effect of extraction parameters on extraction yield of RAMP. (a) extraction time and temperature; (b) extraction time and cellulase concentration; (c) extraction time and ultrasonic power; (d) extraction temperature and cellulase concentration; (e) extraction temperature and ultrasonic power; (f) cellulase concentration and ultrasonic power.
Figure 2. Response surface (3D) showing the effect of extraction parameters on extraction yield of RAMP. (a) extraction time and temperature; (b) extraction time and cellulase concentration; (c) extraction time and ultrasonic power; (d) extraction temperature and cellulase concentration; (e) extraction temperature and ultrasonic power; (f) cellulase concentration and ultrasonic power.
Molecules 20 19837 g002

The DPPH Scavenging Percentage of RAMP in the BBD Experiments

As seen in Table 1, the DPPH scavenging percentage of RAMP (Y2) ranged from 39.11% to 63.88%. The data were analyzed by multiple regression analysis to get the following second-order polynomial equation:
Y2 = 63.53 + 2.95X1 + 3.63X2 + 4.78X3 + 3.64X4 + 1.39X1X2 − 0.06X1X3 − 0.23X1X4 + 1.52X2X3 − 0.40X2X4 − 1.76X3X4 − 8.62X12 − 6.40X22 − 7.46X32 − 6.69X42
ANOVA results of the quadratic model presented in Table 3 show a high R2 of 0.9973 and a low C.V. value of 1.08%, demonstrating that the model can adequately describe the response surface of DPPH scavenging percentage. The high model F-value (376.23) and low p-value (p < 0.0001) suggested the results were highly statistically significant and had a good fit of the model. It can be seen from Table 3 that all the linear coefficients (X1, X2, X3, X4), quadratic term coefficients (X12, X22, X32, X42) and cross product coefficients (X1·X2, X2·X3, X3·X4) were significant model terms, with very small p value (p < 0.001). Moreover, the coefficients of X1·X3, X1·X4, X2·X4 were found non-significant (p > 0.05).
Table 3. ANOVA for the dependent variable: the DPPH scavenging percentage of RAMP.
Table 3. ANOVA for the dependent variable: the DPPH scavenging percentage of RAMP.
SourceSum of SquaresdfMean SquareFSig. (Prob > F)
Corrected model1629.0414116.36376.23<0.0001
X1104.141104.14336.70<0.0001
X2157.911157.91510.55<0.0001
X3274.471274.47887.43<0.0001
X4158.781158.78513.37<0.0001
X1X27.7617.7625.080.0002
X1X30.01410.0140.0470.8323
X1X40.2010.200.650.4320
X2X39.2119.2129.78<0.0001
X2X40.6510.652.100.1698
X3X412.32112.3239.83<0.0001
X12482.071482.071558.67<0.0001
X22265.551265.55858.60<0.0001
X32360.581360.581165.86<0.0001
X42290.381290.38938.89<0.0001
Residual4.33140.31
Lack of fit3.86100.393.270.1321
Pure error0.4740.12
Cor total1633.3728
R20.9973SD0.56
R2adj0.9947C.V.%1.08
Adeq precision60.809
Figure 3a shows the DPPH scavenging percentage with varying extraction time and extraction temperature. From the figure, it can be seen that DPPH scavenging percentage increased as the extraction time and extraction temperature increased during the initial stage and then decreased slightly. Similar trends were observed for the other interactions of the model (Figure 3b–f). According to the model, the maximum DPPH scavenging percentage of 61.49% can be obtained after an extraction time of 85.53 min, at an extraction temperature of 58.29 °C, with a cellulase concentration of 1.84% and 220.05 W of ultrasonic power.
Figure 3. Response surface (3D) showing the effect of extraction parameters on DPPH scavenging percentage of RAMP. (a) extraction time and temperature; (b) extraction time and cellulase concentration; (c) extraction time and ultrasonic power; (d) extraction temperature and ultrasonic power; (e) extraction temperature and cellulase concentration; (f) cellulase concentration and ultrasonic power.
Figure 3. Response surface (3D) showing the effect of extraction parameters on DPPH scavenging percentage of RAMP. (a) extraction time and temperature; (b) extraction time and cellulase concentration; (c) extraction time and ultrasonic power; (d) extraction temperature and ultrasonic power; (e) extraction temperature and cellulase concentration; (f) cellulase concentration and ultrasonic power.
Molecules 20 19837 g003

2.2.2. Optimization Analysis of the UAEE Procedure

According to the previous study, the maximum value of RAMP yield by HWE was 3.13% [35]. We take the HWE maximum yield as our minimum polysaccharide yield in this study. The maximum value of the polysaccharide yield was specified as 50%. Besides, the maximum theoretical value of the DPPH scavenging percentage of RAMP is 100%, and the minimum DPPH scavenging percentage was specified as 20%. Thus, a one-sided transform of polysaccharides yield (d1) and DPPH scavenging percentage (d2) is obtained as follows:
d 1 = { 0 Y i 3.13 [ Y i 3.13 50 3.13 ] 3.13 < Y i < 50 1 Y i 50
d 2 = { 0 Y i 20 [ Y i 20 100 20 ] 20 < Y i < 100 1 Y i 100
The overall desirability D is calculated as:
D = d 1 d 2
By using D as the new response, the optimum values of selected variables can be obtained through regression analysis. In this study, the optimal conditions for highest D (with a D value of 0.328) were: extraction time of 90.54 min, extraction temperature of 57.99 °C, cellulase concentration of 1.95% and 225.29 W of ultrasonic power. The corresponding maximum polysaccharide yield and DPPH scavenging percentage were 13.28% and 59.92%, respectively.

2.2.3. Verification of the Predictive Model

To confirm the suitability of the model equation, three confirmation experiments were conducted under the optimized conditions. Considering the operability in actual production, the optimal conditions were slightly modified as follows: extraction time 90 min, extraction temperature 58 °C, cellulase concentration 1.95% and ultrasonic power 225 W. Under these conditions, the experimental yield of RAMP was 13.18% ± 0.56% (n = 3) and the DPPH scavenging percentage was 60.19% ± 2.99% (n = 3), which matched well with the predicted values of 13.28% and 59.92%, respectively. This confirmed that the model was adequate for optimization of the UAEE process. As a result, RSM coupled with DF approach was considered to be an accurate and decisive tool for predicting the maximum extraction yield and highest antioxidant activity of RAMP using the UAEE technique.

2.3. Comparison with Other Extraction Processes

The results of our comparision with other extraction methods are illustrated in Figure 4. As shown in Figure 4a, the yield of RAM obtained by UAE and EAE was lower than that of UAEE for the same ratio of material and water, extraction time, temperature and pH. In addition, the yield of RAMP obtained by UAEE under the optimal conditions in this study was much higher and obtained faster than with the traditional hot water treatment, which obtained a yield of 5.21% under the conditons of 90 min extraction time and 58 °C temperature. When the extraction time was extended to 180 min and extraction temperature was went up to 80 °C, the RAMP yield by HWE rose to 8.35%, but was still lower than that of UAEE (13.18%). Consequently, UAEE has some obvious time-saving and high efficiency advantages when compared with other extraction processes. Figure 4b shows the DPPH sacvenging percentage among all the extraction processes, which obviously display a similar trend as the yield of RAMP. Thus, it was confirmed that UAEE should be an appropriate and effective extraction technique for obtaining high yields and maximum antioxidant activity of RAMP. Also, the costs of technique is so low that it has a very good industrial application prospects.
Figure 4. The yield (a) and DPPH scavenging percentage; (b) of RAMP of different extraction processes.
Figure 4. The yield (a) and DPPH scavenging percentage; (b) of RAMP of different extraction processes.
Molecules 20 19837 g004

2.4. Preliminary Characterization of RAMP

The total neutral carbohydrate content of RAMP obtained by UAEE under the optimal conditions was estimated at 76.34% by the phenol-sulfuric acid method. RAMP had a negative response to the Bradford method. No absorption was observed at 280 nm in the UV spectrum, which also indicates the absence of protein. Its uronic acid content was under the limit of detection. The FT-IR spectrum of RAMP is shown in Figure 5. A strong broad stretching peak at approximately 3377.80 cm−1 for hydroxyl group and a weak band at approximately 2925.55 cm−1 for the C–H stretching vibration were observed; these two absorption bands are typical characteristic absorptions of polysaccharides [36]. The relatively weak absorption peaks at 1633.04 cm−1 and 1407.05 cm−1 also suggest the characteristic IR absorptions of polysaccharides [37]. The strong band in the region of 1000–1200 cm−1 was attributed to the C–O–C and C–OH bands in the IR spectrum, the peak of 1033.51 cm−1 suggests that the characteristic sugar moieties were of pyranose configuration [36]. Furthermore, no absorption peak was observed at 1740 cm−1, which indicates the absence of uronic acids in the polysaccharide structure [37]. This result further suggest that RAMP is a neutral polysaccharide.
Figure 5. The FT-IR spectrum of RAMP.
Figure 5. The FT-IR spectrum of RAMP.
Molecules 20 19837 g005

3. Experimental Section

3.1. Materials

The rhizoma of Atractylodes macrocephala was purchased from the herb market of Pan’an Country, Zhejiang Province, China. It was authenticated by Professor Jinbao Pu (Zhejiang Academy of Traditional Chinese Medicine).

3.2. Chemical and Apparatus

The cellulase from Trichoderma viride (10,000 U/mg, lot#: C10008391)) was obtained from Shanghai Macklin Biochemical Co., Ltd. (Shanghai, China). d-Glucose and d-galacturonic acid (110833–200503, 111646–200301) were supplied by the National Institutes for Food and Drug Control (Beijing, China). Albumin (Bovine Serum) (10021008) was purchased from Shanghai Shenhang Biological Technology Co., Ltd. (Shanghai, China). 2,2-Diphenyl-1-picrylhydrazyl (DPPH; 100100723) was purchased from Sigma-Aldrich Chemical Company (St. Louis, MO, USA). All solvents and chemicals were at least of analytical grade and obtained from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China).
All UAEE experiments were carried out in an ultrasonic device (JK-300DB, Jiangsu, China) with a usable capacity of 10 L working at a frequency of 40 kHz and an ultrasound power of 100 W to 300 W. A centrifuge (Eppendorf 5415D, Hamburg, Germany) was used for the separation of the extract from the matrix residue. A rotary evaporator (Buchi R-210, Flawil, Switzerland) was used for the concentration of extracts. A UV-visible spectrophotometer (Varian Cary-100, Palo Alto, CA, USA) was used for the analysis of the carbohydrate content, protein content and uronic acid content of the RAMP. A microplate spectrophotometer (Biotek Powerwave XS, Winooski, VT, USA) was used for the analysis of DPPH free radical scavenging activity and a FT-IR spectrometer (Thermo Fisher Nicolet iS10, Waltham, MA, USA) was used for FT-IR spectra analysis of RAMP.

3.3. Methods

3.3.1. UAEE Procedure

The herb material was reduced to a coarse powder with the help of a suitable grinder and passed through a 100 mesh screen. About 200 g of the dry rhizoma powder was doubly extracted with petroleum ether (boiling point: 60–90 °C) at 90 °C for 3 h each time and then pretreated in a Soxhlet extractor with 80% ethanol twice to remove lipids, some colored materials, monosaccharides and oligosaccharides [38]. The residue was then filtered and vacuum dried at 60 °C for 12 h for the following UAEE and conventional extraction.
Pretreated powder (1 g) was doubly extracted with cellulase solution (30 mL, using NaAc-HAc buffer system as our extracting solvent to adjust the pH at 4.6) at the given concentration (ratio of water to material 30 mL/g). The extraction time, extraction temperature, cellulase concentration and ultrasonic power was set according to the experimental design. After removing the matrix residue by centrifuge, the extraction solution was incorporated and concentrated to one-fifth of the initial volume using a rotary evaporator at 65 °C under reduces pressure. The concentrate was precipitated with the addition of anhydrous ethanol to a final concentration of 80% (v/v) and kept at 4 °C for 12 h. The precipitate was collected and deproteinated by the Sevage reagent (1-butanol/chloroform, v/v = 1:4) [39]. After mixing the water solution of polysaccharides and Sevage reagent together (v/v = 5:1), vibrated for 15 min and placed for 2 h to separating the denatured protein. The process was repeated three times and the product obtained then dried to get RAMP.

3.3.2. Conventional Procedure

Hot water extraction (HWE), ultrasonic-assisted extraction (UAE) and enzyme-assisted extraction (EAE) of RAMP were also done to serve as the referance extraction methods. The specific conditions of the different extraction processes are listed in Table 4. All the extracts was obtained and treated according to the procedure mentioned in Section 3.3.1, and both of the yield and DPPH scavenging percentage were determinated.
Table 4. The differrent extraction process of RAMP.
Table 4. The differrent extraction process of RAMP.
Extraction ProcessesExtraction Conditions
ABCDE
HWE1--90587.0
HWE2--180807.0
EAE1.95-90584.6
UAE-22590584.6
UAEE1.9522590584.6
A: cellulase concentration (%), B: ultrasonic power (W), C: extraction time (min), D: extraction temperature (°C), E: pH.

3.3.3. Determination of the Yield of RAMP

The content of RAMP was measured by the phenol-sulfuric acid method [40] using d-glucose as a standard. The yield (%) of RAMP was then calculated as follows:
Y 1 ( % ) = C W × 100 %
where Y1 is the yield of RAMP, C is the weight of polysaccharides and W is the weight of raw material.

3.3.4. Determination of Antioxidant Acitivity of RAMP

The antioxidant activity was studied through the evaluation of the free radical-scavenging effect on the DPPH radical, which was measured by the procedure described previously [41] with slight modification. Briefly, ascorbic acid (Vc) was used as the control, and 30 μL of RAMP solutions (0.1% in water) or control were mixed with 170 μL of the ethanolic DPPH (0.025 g/L). The mixture was shaken vigorously and incubated at 25 °C in the dark for 30 min. Absorption of the samples was measured on the microplate spectrophotometer at 517 nm. The scavenging percentage was calculated as follows:
Y 2 ( % ) = A c ( A i A j ) A c × 100 %
where Y2 is the scavenging percentage of RAMP, Ac is the absorbance of DPPH solution without sample, Ai is the absorbance of the sample mixed with DPPH solution and Aj is the absorbance of the sample without DPPH solution.

3.3.5. Experimental Design

Single Factor Experimental Design

The single factor experiment was performed in a designed conditions. During the optimization of experimental factors, one factor was changed while the other factors kept constant in each experiment. In this study, four factors including extraction time, extraction temperature, cellulase concentration and ultrasonic power were selected for the single factor experiments. In detail, the extraction time range was from 10 min to 90 min, the extraction temperature ranged from 30 °C to 70 °C, the cellulase concentration ranged from 0.5% to 2.5% (w/v) and the ultrasonic power range was from 100 W to 300 W. The effect of each factor was evaluated by determining the extraction yield of RAMP. All the experiments were repeated three times.

Box-Behnken Design

On the basis of the single factor experiments, a Box-Behnken design (BBD) with four factors and three levels was used for the further optimization of the UAEE conditions. Four independent variables including extraction time, extraction temperature, cellulase concentration and ultrasonic power were designated as X1, X2, X3 and X4, respectively. The ranges of values were based on the results of the preliminary experiments, as shown in Table 5. The response variables were the extraction yield of RAMP (Y1) and the DPPH scavenging percentage (Y2). As shown in Table 1, there were a total 29 runs based on the BBD with five center points performed in random order with triplicates in each run.
Table 5. Independent variables and their levels in Box-Behnken design
Table 5. Independent variables and their levels in Box-Behnken design
Independent variablesSymbolLevel
−101
Extraction time (min)X1105090
Extraction temperature (°C)X2305070
Cellulase concentration (%)X30.51.52.5
Ultrasonic power (W)X4100200300

3.3.6. Statistical Analysis and Optimization

The parameters of the response equation and analysis of variance (ANOVA) were performed by Design Expert Software (Version 8.0.6). A second order polynomial model used to fit the response to the independent variables is shown below:
Y = β 0 + i = 1 3 β i X i + i = 1 3 β i i X i 2 + i < j 3 β i j X i X j
where Y is the response, β0 is the intercept parameter and βi, βii, βij are the coefficients of the parameters for linear, squared and interaction effects, respectively.
The statistical significance for each term in the polynomial was evaluated by computing the F-value at a probability p of 0.05. The regression coefficients were then used to make statistical calculations and generate contour maps from the regression models. However, for multi-response, a desirability function approach can be used to transformed several response variables into a desirabitity function, which can be optimized by univariate techniques. A modified desirability approach, proposed by Derringer and Suich [22] is defined as:
D = ( d 1 w 1 d 2 w 2 d 3 w 3 d 4 w 4 d n w n ) 1 / w i
where wi is the relative weight of the ith response, D is the overall desirability , and di is an individual response desiralibity. Then, the optimal setting is determined by the following, which is described previously [42,43].
d i = { 0 Y i Y i min [ Y i Y i min Y i max Y i min ] r Y i min < Y i < Y i max 1 Y i Y i max
where Yi is the response value, Yi-min is the minimum acceptable value for response i, Yi-max is the maximum acceptable value for response i, and r is a weight used to determine scale of desirability and equals 1 in this work.

3.3.7. Preliminary Characterization of RAMP

The total neutral carbohydrate content was determinated by the phenol-sulfuric acid colorimetric method described in Section 3.3.3. The protein content was measured by Coomassie brilliant blue reaction [44] and the UV scanning spectrum at 280 nm. Uronic acid content was determinated by the m-hydroxydiphenyl method using d-galacturonic acid as standard [45]. The characteristic absorption of RAMP was identified by the FT-IR spectrum [36]. The RAMP powder was mixed with KBr powder, ground and pressed for FT-IR measurement in the wavenumber range from 400 to 4000 cm−1.

4. Conclusions

An UAEE process has been optimized for effective extraction of RAMP with high antioxidant activity. The maximum D value of 0.375, along with the maximum yield (13.73%) and scavenging percentage (62.18%) were achieved after an extraction time of 82.04 min, using an extraction temperature of 58.67 °C, a cellulase concentration of 1.99% and 230.67 W of ultrasonic power. These values were further validated by confirmatory experiments to see the efficacy of the model predictability and found to be in good agreement with the predicted values. Compared to other extraction methods, both the extraction yield and DPPH scavenging percentage of RAMP obtained by UAEE was favorable and the method appeared to be time-saving and of high efficiency. These results demostrated that UAEE is an appropriate and effective extraction technique for RAMP. Moreover, RSM with DF approach has been proved to be adequate for the design and optimization of the extraction parameters for RAMP. This works offer a wide range of implications for guiding the design and operation of polysaccharide extraction processes and may have a great impact on the polysaccharide extraction industries.

Acknowledgments

We greatly acknowledge the financial support from the Traditional Chinese Medicine of Zhejiang Province Science and Technology Plan Project (No. 2013ZA006). We are also grateful for the Construction of Key Laboratory of Research and Development of Chinese Medicine of Zhejiang Province (No.2008F3036).

Author Contributions

Bo-Hou Xia, Pan Xu and Jin-Bao Pu designed the research and wrote the paper, Yi-Juan Hu and Hong-Jian Zhang analyze the experiment results, Jing Chen, Jie Zhou and Wei-Qing Liang performed the research. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yao, C.M.; Yang, X.W. Bioactivity-guided isolation of polyacetylenes with inhibitory activity against NO production in LPS-activated RAW264.7 macrophages from the rhizomes of Atractylodes macrocephala. J. Ethnopharmacol. 2014, 151, 791–799. [Google Scholar] [CrossRef] [PubMed]
  2. Chinese Pharmacopoeia Commission. Pharmacopeia of People’s Republic of China; Chemical Industry Press: Beijing, China, 2010; Volume 1, p. 95. [Google Scholar]
  3. Chen, Z.L. The Acetylenes from Atractylodes macrocephala. Planta Med. 1987, 53, 493–494. [Google Scholar] [CrossRef] [PubMed]
  4. Kitajima, J.; Kamoshita, A.; Ishikawa, T.; Takano, A.; Fukuda, T.; Isoda, S.; Ida, Y. Glycosides of Atractylodes ovata. Chem. Pharm. Bull. 2003, 51, 1106–1108. [Google Scholar] [CrossRef] [PubMed]
  5. Dong, H.; He, L.; Huang, M.; Dong, Y. Anti-inflammatory components isolated from Atractylodes macrocephala Koidz. Nat. Prod. Res. 2008, 22, 1418–1427. [Google Scholar] [CrossRef] [PubMed]
  6. Shan, G.S.; Zhang, L.X.; Zhao, Q.M.; Xiao, H.B.; Zhuo, R.J.; Xu, G.; Jiang, H.; You, X.M.; Jia, T.Z. Metabolomic study of raw and processed Atractylodes macrocephala Koidz by LC-MS. J. Pharm. Biomed Anal. 2014, 98, 74–84. [Google Scholar] [CrossRef] [PubMed]
  7. Liu, J.; Chen, X.; Yue, C.; Hou, R.; Chen, J.; Lu, Y.; Li, X.; Li, R.; Liu, C.; Gao, Z.; et al. Effect of selenylation modification on immune-enhancing activity of Atractylodes macrocephala polysaccharide. Int. J. Biol. Macromol. 2015, 72, 1435–1440. [Google Scholar] [CrossRef] [PubMed]
  8. Xie, F.; Li, Y.; Su, F.; Hu, S. Adjuvant effect of Atractylodis macrocephalae Koidz. polysaccharides on the immune response to foot-and-mouth disease vaccine. Carbohydr. Polym. 2012, 87, 1713–1719. [Google Scholar] [CrossRef]
  9. Xie, F.; Sakwiwatkul, K.; Zhang, C.; Wang, Y.; Zhai, L.; Hu, S. Atractylodis macrocephalae Koidz. polysaccharides enhance both serum IgG response and gut mucosal immunity. Carbohydr. Polym. 2013, 91, 68–73. [Google Scholar] [CrossRef] [PubMed]
  10. Shi, N.; Su, J.; Yang, Z.B.; Lv, G.Y.; Chen, S.H. Antioxidantive effect of polyasaccharides from Atractylodes macrocephalaon in d-galactose-induced aging mice. Chin. J. New Drugs. 2014, 23, 577–581. [Google Scholar]
  11. Hu, W.X.; Xiang, Q.; Wen, Z.; He, D.; Wu, X.M.; Hu, G.Z. Neuroprotective effect of Atractylodes macrocephala polysaccharides in vitro on neuronal apoptosis induced by hypoxia. Mol. Med. Rep. 2014, 9, 2573–2581. [Google Scholar] [PubMed]
  12. Shan, J.J.; Tian, G.Y. Studies on physico-chemical properties and hypoglycemic activity of complex polysaccharide AMP-B from Atractylodes macrocephala Koidz. Acta Pharm. Sinica 2003, 38, 438–441. [Google Scholar]
  13. Zhang, P.J.; Jin, C.; Lang, J. Effect of Atractylodes macrocephalaon polysaccharides for reducing liver ischemia/reperfusion injury in rats after orthotopic liver autotransplantation. Chin. J. Integr. Tradit. West. Med. 2010, 30, 1193–1196. [Google Scholar]
  14. Zhao, Z.Y.; Zhang, Q.; Li, Y.F.; Dong, L.L.; Liu, S.L. Optimization of ultrasound extraction of Alisma orientalis polysaccharides by response surface methodology and their antioxidant activities. Carbohydr. Polym. 2015, 119, 101–109. [Google Scholar] [CrossRef] [PubMed]
  15. Lai, J.; Wang, H.; Wang, D.; Fang, F.; Wang, F.; Wu, T. Ultrasonic extraction of antioxidants from Chinese sumac (Rhus typhina L.) fruit using response surface methodology and their characterization. Molecules 2014, 19, 9019–9032. [Google Scholar] [CrossRef] [PubMed]
  16. Moura, J.M.L.N.; Campbell, K.; Mahfuz, A.; Jung, S.; Glatz, C.E.; Johnson, L. Enzyme-assisted aqueous extraction of oil and protein from soybeans and cream de-emulsification. J. Am. Oil Chem. Soc. 2008, 85, 985–995. [Google Scholar] [CrossRef]
  17. Liu, Y.; Gong, G.; Zhang, J.; Jia, S.; Li, F.; Wang, Y.; Wu, S. Response surface optimization of ultrasound-assisted enzymatic extraction polysaccharides from Lycium barbarum. Carbohydr. Polym. 2014, 110, 278–284. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, T.; Li, C.Q.; Zhang, H.; Li, J.W. Response surface optimized extraction of 1-deoxynojirimycin from mulberry leaves (Morus alba L.) and preparative separation with resins. Molecules 2014, 19, 7040–7056. [Google Scholar] [CrossRef] [PubMed]
  19. Salmasnia, A.; Baradaran kazemzadeh, R.; Mohajer Tabrizi, M. A novel approach for optimization of correlated multiple responses based on desirability function and fuzzy logics. Neurocomputing 2012, 91, 56–66. [Google Scholar] [CrossRef]
  20. Khuri, A.; Cornell, J. Response Surfaces: Designs and Analyses, 2nd ed.; Dekker: New York, NY, USA, 1996; Volume 12, pp. 72–75. [Google Scholar]
  21. He, Z.; Zhu, P.F.; Park, S.H. A robust desirability function method for multi-response surface optimization considering model uncertainty. Eur. J. Oper. Res. 2012, 221, 241–247. [Google Scholar] [CrossRef]
  22. Derringer, G. A Balancing Act: Optimizing a Product’s Properties; Quality Press: Perth, Australia, 1994; Volume 27, pp. 51–58. [Google Scholar]
  23. Derringer, G.; Suich, R. Simultaneous optimization of several response variables. J. Qual. Technol. 1980, 12, 214–219. [Google Scholar]
  24. Islam, M.A.; Sakkas, V.; Albanis, T.A. Application of statistical design of experiment with desirability function for the removal of organophosphorus pesticide from aqueous solution by low-cost material. J. Hazard. Mater. 2009, 170, 230–238. [Google Scholar] [CrossRef] [PubMed]
  25. Heidari, H.; Razmi, H. Multi-response optimization of magnetic solid phase extraction based on carbon coated Fe3O4 nanoparticles using desirability function approach for the determination of the organophosphorus pesticides in aquatic samples by HPLC-UV. Talanta 2012, 99, 13–21. [Google Scholar] [CrossRef] [PubMed]
  26. Jeong, I.J.; Kim, K.J. An interactive desirability function method to multiresponse optimization. Eur. J. Oper. Res. 2009, 195, 412–426. [Google Scholar] [CrossRef]
  27. Candioti, L.V.; de Zan, M.M.; Camara, M.S.; Goicoechea, H.C. Experimental design and multiple response optimization. Using the desirability function in analytical methods development. Talanta 2014, 124, 123–138. [Google Scholar] [CrossRef] [PubMed]
  28. Chen, X.P.; Tang, Q.C.; Chen, Y.; Wang, W.X.; Li, S.B. Simultaneous extraction of polysaccharides from Poria cocos by ultrasonic technique and its inhibitory activities against oxidative injury in rats with cervical cancer. Carbohydr. Polym. 2009, 79, 409–413. [Google Scholar] [CrossRef]
  29. Chen, F.; Sun, Y.; Zhao, G.; Liao, X.; Hu, X.; Wu, J.; Wang, Z. Optimization of ultrasound-assisted extraction of anthocyanins in red raspberries and identification of anthocyanins in extract using high-performance liquid chromatography-mass spectrometry. Ultrason. Sonochem. 2007, 14, 767–778. [Google Scholar] [CrossRef] [PubMed]
  30. Fu, X.Y.; Xue, C.H.; Miao, B.C.; Li, Z.J.; Gao, X.; Yang, W.G. Characterization of proteases from the digestive tract of sea cucumber (Stichopus japonicus): High alkaline protease activity. Aquaculture 2005, 246, 321–329. [Google Scholar] [CrossRef]
  31. Ticar, B.F.; Rohmah, Z.; Ambut, C.V.; Choi, Y.J.; Mussatto, S.I.; Choi, S.D. Enzyme-assisted extraction of anticoagulant polysaccharide from Liparis tessellates eggs. Int. J. Biol. Macromol. 2015, 74, 601–607. [Google Scholar] [CrossRef] [PubMed]
  32. Rhein-Knudsen, N.; Ale, M.T.; Meyer, A.S. Seaweed hydrocolloid production: An update on enzyme assisted extraction and modification technologies. Mar. Drugs 2015, 13, 3340–3359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Tahmouzi, S. Optimization of polysaccharides from Zagros oak leaf using RSM: Antioxidant and antimicrobial activities. Carbohydr. Polym. 2014, 106, 238–246. [Google Scholar] [CrossRef] [PubMed]
  34. Kuo, C.H.; Chen, B.Y.; Liu, Y.C.; Chang, C.M.; Deng, T.S.; Chen, J.H.; Shieh, C.J. Optimized ultrasound-assisted extraction of phenolic compounds from Polygonum cuspidatum. Molecules 2013, 19, 67–77. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Liu, J.H.; Ren, H.L.; Guan, B.; Zheng, W.C.; Gao, F. Response surface optimization of extraction process of Atractylodes polysaccharide. Chin. Tradit. Pat. Med. 2008, 30, 667–670. [Google Scholar]
  36. Yi, P.; Li, N.; Wan, J.B.; Zhang, D.; Li, M.; Yan, C. Structural characterization and antioxidant activity of a heteropolysaccharide from Ganoderma capense. Carbohydr. Polym. 2015, 121, 183–189. [Google Scholar] [CrossRef] [PubMed]
  37. Pei, J.J.; Wang, Z.B.; Ma, H.L.; Yan, J.K. Structural features and antitumor activity of a novel polysaccharide from alkaline extract of Phellinus linteus mycelia. Carbohydr. Polym. 2015, 115, 472–477. [Google Scholar] [CrossRef] [PubMed]
  38. Zhao, Z.; Xu, X.; Ye, Q.; Dong, L. Ultrasound extraction optimization of Acanthopanax senticosus polysaccharides and its antioxidant activity. Int. J. Biol. Macromol. 2013, 59, 290–294. [Google Scholar] [CrossRef] [PubMed]
  39. Savag, M.G.; Lackman, D.B.; Smolens, J. The isolation of the components of streptococcal nucleoproteins in serologically active from. J. Biol. Chem. 1938, 124, 425–436. [Google Scholar]
  40. Kettlitz, B. Carbohydrate Analysis: A Practical Approach; Chaplin, M.F., Kennedy, J.E., Eds.; IRL Press: Oxford, UK, 1986; p. 3. [Google Scholar]
  41. Vallverdu-Queralt, A.; Regueiro, J.; Martinez-Huelamo, M.; Rinaldi Alvarenga, J.F.; Leal, L.N.; Lamuela-Raventos, R.M. A comprehensive study on the phenolic profile of widely used culinary herbs and spices: Rosemary, thyme, oregano, cinnamon, cumin and bay. Food Chem. 2014, 154, 299–307. [Google Scholar] [CrossRef] [PubMed]
  42. Costa, N.R.; Lourenco, J.; Pereira, Z.J. Desirability function approach: A review and performance evaluation in adverse conditions. Chemometr. Intell. Lab. 2011, 107, 234–244. [Google Scholar] [CrossRef]
  43. Shi, X.Y.; Li, W.W.; Yu, H.Q. Optimization of H2 photo-fermentation from benzoate by Rhodopseudomonas palustris using a desirability function approach. Int. J. Hydrog. Energy 2014, 39, 4244–4251. [Google Scholar] [CrossRef]
  44. Bradford, M.M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248–254. [Google Scholar] [CrossRef]
  45. Blumenkrantz, N.; Asboe-Hansen, G. New method for quantitative determination of uronic acids. Anal. Biochem. 1973, 54, 484–489. [Google Scholar] [CrossRef]
  • Sample Availability: Samples of the compounds are available from the authors.

Share and Cite

MDPI and ACS Style

Pu, J.-B.; Xia, B.-H.; Hu, Y.-J.; Zhang, H.-J.; Chen, J.; Zhou, J.; Liang, W.-Q.; Xu, P. Multi-Optimization of Ultrasonic-Assisted Enzymatic Extraction of Atratylodes macrocephala Polysaccharides and Antioxidants Using Response Surface Methodology and Desirability Function Approach. Molecules 2015, 20, 22220-22235. https://doi.org/10.3390/molecules201219837

AMA Style

Pu J-B, Xia B-H, Hu Y-J, Zhang H-J, Chen J, Zhou J, Liang W-Q, Xu P. Multi-Optimization of Ultrasonic-Assisted Enzymatic Extraction of Atratylodes macrocephala Polysaccharides and Antioxidants Using Response Surface Methodology and Desirability Function Approach. Molecules. 2015; 20(12):22220-22235. https://doi.org/10.3390/molecules201219837

Chicago/Turabian Style

Pu, Jin-Bao, Bo-Hou Xia, Yi-Juan Hu, Hong-Jian Zhang, Jing Chen, Jie Zhou, Wei-Qing Liang, and Pan Xu. 2015. "Multi-Optimization of Ultrasonic-Assisted Enzymatic Extraction of Atratylodes macrocephala Polysaccharides and Antioxidants Using Response Surface Methodology and Desirability Function Approach" Molecules 20, no. 12: 22220-22235. https://doi.org/10.3390/molecules201219837

APA Style

Pu, J. -B., Xia, B. -H., Hu, Y. -J., Zhang, H. -J., Chen, J., Zhou, J., Liang, W. -Q., & Xu, P. (2015). Multi-Optimization of Ultrasonic-Assisted Enzymatic Extraction of Atratylodes macrocephala Polysaccharides and Antioxidants Using Response Surface Methodology and Desirability Function Approach. Molecules, 20(12), 22220-22235. https://doi.org/10.3390/molecules201219837

Article Metrics

Back to TopTop