Next Article in Journal
Phytochemical Characterization and Heavy Metal and Thermal Analyses of Saussurea hypoleuca Root and Evaluation of Its Anthelmintic and Antioxidant Activity In Vitro and In Silico
Next Article in Special Issue
Inhibitory Effect and Mechanism of Dill Seed Essential Oil on Neofusicoccum parvum in Chinese Chestnut
Previous Article in Journal
Isotopic Characterization of Italian Industrial Hemp (Cannabis sativa L.) Intended for Food Use: A First Exploratory Study
Previous Article in Special Issue
Enantiomeric Separation and Molecular Modelling of Bioactive 4-Aryl-3,4-dihydropyrimidin-2(1H)-one Ester Derivatives on Teicoplanin-Based Chiral Stationary Phase
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimization of Steam Distillation Process and Chemical Constituents of Volatile Oil from Angelicaesinensis Radix

1
Key Laboratory of Modern Preparation of TCM, Ministry of Education, Jiangxi University of Chinese Medicine, No. 1688 Meiling Avenue, Xinjian District, Nanchang 330004, China
2
Pharmaceutical Informatics Institute, College of Pharmaceutical Science, Zhejiang University, Hangzhou 310058, China
3
Jinhua Institute of Zhejiang University, Jinhua 321016, China
*
Authors to whom correspondence should be addressed.
Submission received: 22 April 2022 / Revised: 18 May 2022 / Accepted: 24 May 2022 / Published: 30 May 2022

Abstract

:
In this study, the steam distillation process of volatile oil from Angelicae sinensis Radix was optimized according to the concept of quality-by-design. A homemade glass volatile oil extractor was used to achieve better cooling of the volatile oil. First, the soaking time, distillation time, and liquid–material ratio were identified as potential critical process parameters by consulting the literature. Then, the three parameters were investigated by single factor experiments. The volatile oil yield increased with the extension in the distillation time, and first increased and then decreased with the increase in soaking time and liquid–material ratio. The results confirmed that soaking time, distillation time, and liquid–material ratio were all critical process parameters. The kinetics models of volatile oil distillation from Angelicae sinensis Radix were established. The diffusion model of spherical particle was found to be the best model and indicated that the major resistance of mass transfer was the diffusion of volatile oil from the inside to the surface of the medicinal herb. Furthermore, the Box–Behnken experimental design was used to study the relationship between the three parameters and volatile oil yield. A second-order polynomial model was established, with R2 exceeding 0.99. The design space of the volatile oil yield was calculated by a probability-based method. In the verification experiments, the average volatile oil yield reached 0.711%. The results showed that the model was accurate and the design space was reliable. In this study, 21 chemical constituents of volatile oil from Angelicae sinensis Radix were identified by gas chromatograph-mass spectrometer(GC-MS), accounting for 99.4% of the total volatile oil. It was found that the content of Z-ligustilide was the highest, accounting for 85.4%.

1. Introduction

Angelicaesinensis Radix is the dry root of Angelica sinensis (Oliv.) Diels [1]. It is widely used as a herbal medicine and has the effects of clearing heat and promoting dieresis, invigorating qi, and blood, etc. Volatile oil is an important component of Angelicae sinensis Radix [2,3]. The content of neutral oil in the volatile oil of Angelicae sinensis Radix is the highest. The volatile oil of Angelicae sinensis Radix mainly contains Z-ligustilide and other components [4,5]. Volatile oil is also considered as the main effective component of Angelicae sinensis Radix, which has the effects of treating hypertension [6], analgesic and anti-inflammatory [7], and anti-tumor [8]. The preparation process of the Angelicae sinensis Radix volatile oil exists in the production of many Chinese medicines such as Danggui Tiaojing granules, Ruhe Sanjie tablets, Compound Herba Leonuri capsules, and Yangxueyin oral liquid, which are all included in the 2020 edition of the Chinese Pharmacopoeia (1st Part).
At present, the reported preparation methods of angelica volatile oil include steam distillation, supercritical CO2 extraction, organic solvent extraction, microwave assisted extraction, ultrasound assisted extraction, and so on [9]. Microwave assisted extraction has the advantages of rapidity [10], but a high microwave power may lead to a decrease in the yield of essential oil because of the loss of active compounds [9]. Ultrasound assisted extraction has the advantages of less solvent consumption and better efficiency [11], but there are difficulties in scale-up production. In industry, steam distillation is most widely used to obtain angelica volatile oil [4]. Steam distillation has the advantages of simple equipment, low cost, and a safe solvent [12,13,14]. In the Chinese Pharmacopoeia, there are more than 30 preparations that use the volatile oil of Angelicae sinensis Radix, where more than 80% used the steam distillation process. It is also the method used to determine the volatile oil content of drugs in the Chinese Pharmacopoeia.
However, there are some practical problems in the distillation of volatile oil from Angelicae sinensis Radix such as low yield, easy emulsification, and long-time consumption, which bring great challenges to the manufacturing processes and the quality control of Chinese medicines containing Angelicae sinensis Radix. Therefore, it is necessary to optimize the distillation process of volatile oil from Angelicae sinensis Radix.
In recent years, the quality-by-design concept has been widely applied to optimize the pharmaceutical processes of Chinese medicines [15]. Its implementation steps include: defining the critical process parameters (CPPs), establishing quantitative models of CPPs and pharmaceutical process evaluation indices, establishing the design space, determining the quality control strategies, and continuous improvement, etc. [16]. Before building mathematical models, response surface methodology (RSM) is usually used to collect data with small number of experiments [10]. The Box–Behnken design is often used in many studies [17] as it suggests how to select points from a three-level factorial arrangement, which allows for the efficient estimation of the first- and second-order coefficients of the mathematical model [18].
In this work, steam distillation was used to prepare the volatile oil from Angelicae sinensis Radix. A volatile oil extractor with enhanced cooling was used to collect more volatile oil. The CPPs of the distillation process were determined, and the distillation kinetic models were established. According to the results, the main mass transfer resistance of the volatile oil distillation was determined. Then, taking the volatile oil yield as the index, the Box–Behnken experimental design was used to study the effects of the distillation time, soaking time, and liquid–material ratio [10]. The design space of the distillation process of volatile oil from Angelicae sinensis Radix was obtained with a probability-based method, which was verified by experiments. Finally, the chemical constituents in the volatile oil of Angelicae sinensis Radix were determined by GC-MS, and the results were compared with those reported in the literature.

2. Materials and Methods

2.1. Materials, Reagents, and Instruments

Angelicae sinensis Radix was obtained from the Jiangxi Zhangshu Tianqitang Chinese Herbal Pieces Co., Ltd. (batch number: 2009007, origin: Gansu). All id the above medicinals were identified by Professor Ge Fei from the Identification Teaching and Research Section of the Jiangxi University of Chinese Medicine. The moisture content of Angelica sinensis was 10.0 ± 0.5%. Anhydrous sodium sulfate (batch number: 180408) and anhydrous ethanol (batch number: 200707) were purchased from Xilong Scientific Co., Ltd. (Shantou, China). Distilled water was self-made in the laboratory.

2.2. Distillation Method

Gansu Province is a genuine producing area of Angelicae sinensis Radix. Therefore, the Angelicae sinensis Radix from Gansu Province was used in this work. The steam distillation was used to prepare the volatile oil from Angelicae sinensis Radix. Figure 1 shows the essential oil extractor. With the homemade essential oil extractor, better cooling could be achieved. Compared to that of using a conventional extractor, more essential oil can be extracted.
A proper amount of Angelicae sinensis Radix was taken and then crushed with a crusher (FD200T, Shanghai Traditional Chinese Medicine Machinery Factory, Shanghai, China). A total of 100 g of Angelicae sinensis Radix powder was taken and added into a 2000 mL flask. Water and several glass beads were added into the flask. After that, the flask was shaken to wet the herbal powders and was connected with a condenser tube. Before the experiment, water was added into the volatile oil extractor. An electric heating jacket (KDM type, Shandong Zhencheng Hualu Electric Equipment Co., Ltd., Qingdao, China) was used to heat the flask. After boiling, the electric heating jacket was adjusted to keep boiling. The heating was stopped after a certain time of distillation. After more than 1 h, the height of the volatile oil was measured with a vernier caliper (0–150 mm, 3V type, Guilin Ganglu Digital Measurement and Control Co., Ltd., Guilin, China). The volume of volatile oil was calculated, then the volatile oil was collected. The obtained volatile oil was dehydrated and dried with anhydrous sodium sulfate overnight to obtain an oil-like product until the amount of volatile oil no longer increased. Thee Angelicae sinensis Radix volatile oil was sealed in brown reagent bottles and stored in a refrigerator at 4 °C. The volatile oil yield was calculated using Equation (1).
Volatile oil yield (%) = [volume of oil (mL)/weight of decoction pieces (g)] × 100%

2.3. Optimization of Distillation Process Parameters

After the literature review, the process parameters of distilling the volatile oil from herbal medicines by steam distillation were found and are shown in an Ishikawa diagram (Figure 2).
According to the pre-experiments, it was considered that the potential CPPs affecting the volatile oil yield by steam distillation were distillation time (A), soaking time (B), and the liquid–material ratio (C). In this study, the above three factors were investigated by single factor experiments. During the investigation, a total of 100 g of Angelicae sinensis Radix powder was taken, soaked for a certain time, and then heated to distill the volatile oil. After reaching the distillation time, the heating was stopped. After the volatile oil was cooled, the volume of the volatile oil was accurately read, and the volatile oil yield was calculated. When the distillation time was studied, the liquid–material ratio was 10:1 mL·g−1 and the soaking time was 2 h. When the soaking time was studied, the ratio of the liquid to solid was 10:1 mL·g−1 and the distillation time was 4 h. When the liquid–material ratio was studied, the distillation time was 4 h and the soaking time was 2 h.
Next, the Box–Behnken design was used to optimize the distillation process of the volatile oil from Angelicae sinensis Radix with the distillation time, soaking time, and liquid–material ratio as the factors. Compared with some other commonly used response surface designs, the number of experiments in the Box–Behnken design is relatively smaller when there are three 3-level factors. For example, there are 29, 17, and 13 runs for the full factorial design, central composite design, and the Box–Behnken design when there are three repetitions of the center point, respectively. The volatile oil yield was considered as the evaluation index. The factor levels are shown in Table 1 and Table 2.

2.4. Analysis of Chemical Constituents of Volatile Oil from Angelicae sinensis Radix

A total of 100 μL of volatile oil was measured accurately, diluted to 10 mL with ethanol, filtered with a 0.22 μm microporous filter membrane, and loaded into a sample bottle. The chemical constituents of the volatile oil samples were determined by gas chromatography-mass spectrometry (Agilent 7890A/5975C, Agilent, Santa Clara, CA, USA).
The gas chromatographic conditions were as follows [19]: HP-5MS capillary quartz column (30 m × 250 μm × 0.25 μm); injection volume: 1.0 μL; the inlet temperature: 230 °C; no diversion; heating procedure: 80 °C, keeping for 0 min; heating at 3 °C min−1 to 167 °C for 2.5 min, at 2 °C min−1 to 180 °C for 0 min, at 30 °C min−1 to 280 °C for 1 min; the transmission line temperature: 300 °C; the carrier gas was high purity helium, the flow rate was 1 mL min−1, and the injection volume was 1 μL. Mass spectrometry conditions: Ion source temperature: 230 °C, ionization source was EI, electron energy: 70 eV; the quadrupole temperature: 150 °C, and the scanning quality range was 30–550 amu.
The identification of volatile compounds was performed by computer matching their mass spectra with those stored in a digital library of mass spectral data (NIST 14). The identification results were tentatively identified by the EI-MS spectrum and further experiments are planned to confirm their identification by authentic standards.

2.5. Data Processing Method

2.5.1. Kinetic Models

Volatile oil yield was fitted by Equations (2)–(4), as shown below.
First-order kinetic model:
C = C e q 1 e k t
where C e q is the concentration of solution at equilibrium; k is the total distillation rate constant; and t   is time.
Peleg’s model:
C = t k 1 + k 2 t
where k1 is the mass transfer rate constant and k2 is the concentration of the solution at equilibrium.
The diffusion model of spherical particle:
C   = C e q 1 6 π 2 n = 1 1 n 2 exp D n 2 π 2 t R 2
where D is the solute diffusion coefficient in the solvent and R is the particle radius. In the calculation, only the first three terms (n = 3) were taken for the sum of the infinite order.
In order to evaluate the fitting effect of each model, R2 can be calculated according to Equation (5).
R 2 = 1 C a c t C f i t 2 C a c t C a c t ¯ 2
where Cact is the experimental value; C a c t ¯ is the average value of the measured value, and Cfit is the model fitting value. The software MATLAB R2019b (American Math Works Company) was used to analyze and process the data. Except for R2, other indices of RMSE, MSE, SSE, and MAE were also calculated according to Equations (6)–(9).
R M S E = C a c t C f i t 2 n
  M S E = C a c t C f i t 2 n
S S E = C a c t C f i t 2
A E = C a c t C f i t n
The software Microsoft Office Excel (American Microsoft Company) was used to analyze and process the data.

2.5.2. Data Processing of Volatile Oil Distillation Rate Obtained from Experimental Design

The software “Design-Expert 8.0.6” was used to analyze the experimental data collected from the Box–Behnken designed experiments. The volatile oil yield was taken as the evaluation index (Y), and a second-order polynomial equation fitting was carried out. Equation (10) was adopted as the mathematical model between the evaluation index and the three process parameters.
Y = b 0 + b 1 A + b 2 B + b 3 C + b 4 AB + b 5 AC + b 6 BC + b 7 A 2 + b 8 B 2 + b 9 C 2
where b1–9 is the partial regression coefficient and b0 is the intercept.

2.5.3. Optimization of Distillation Parameters of Volatile Oil by the Monte Carlo Method

The design space was calculated by the Monte Carlo method with a parameter optimization software compiled by MATLAB R2018b (Math Works Company of America) [20]. According to previous work, the design space calculated with this method is more reliable [20]. Three factors affecting the volatile oil yield were simulated randomly. The combination of parameters with the probability of no less than 0.80 of attaining the preset volatile oil yield value was taken as the design space. In the calculation, the step lengths of distillation time, soaking time, and liquid-–material ratio were set to 0.04, 0.02, and 0.08, respectively. The simulations were conducted 2000 times.

3. Results and Discussion

3.1. Critical Process Parameters of Volatile Oil Distillation

The distillation time was changed to 2, 4, 6, 8, and 10 h, respectively, and the volatile oil yield is shown in Figure 3. The volatile oil yield increased continuously within the first 8 h of distillation time. When the distillation time reached 8 h, the volatile oil yield reached 0.68%, and the yield tended to be stable when the distillation time was prolonged.
The soaking time was changed to 0, 1, 2, 3, and 4 h, the yield of volatile oil from Angelicae sinensis Radix is shown in Figure 4. When the soaking time was 3 h, a maximum oil yield of 0.67% was obtained. After soaking time reached 3 h, the oil yield decreased with the extension in the soaking time. The reason may be that the aqueous extract became viscous when soaking for too long, which was not conducive to the distillation of the volatile oil.
The liquid–material ratio was changed to 6:1, 8:1, 10:1, 12:1, and 14:1 (mL·g−1), respectively, and the volatile oil yield is shown in Figure 5. When the liquid–material ratio was not higher than 10:1, the volatile oil yield increased as the liquid–material ratio increased. The reason may lie in the more uniform heating when the liquid–material ratio increased, which was beneficial to the diffusion and dissolution of the volatile oil. In Figure 4, when the liquid–material ratio was greater than 10:1, the volatile oil yield showed an obvious downward trend, which may be due to the loss caused by the dissolution of the volatile oil in water.
The above results indicate that the distillation time, soaking time, and liquid–material ratio all had a great influence on the yield of the volatile oil from Angelicae sinensis Radix and were all indeed the CPPs.

3.2. Study on Kinetics of Volatile Oil Distillation

Three models were used to fit the kinetic data of the volatile oil distillation, and the fitting results are shown in Figure 2. The first-order kinetic model fitting results C e q and k were 0.6644% and 0.5861 h−1, respectively, and R2 was 0.8198. The Peleg’s model fitted k1 and k2 were 1.6156 h·%−1 and 1.3105 %−1, respectively, and the R2 was 0.9338. The diffusion model of spherical particle fitted C e q and D R 2 were 2.5317% and 0.0007 h−1, respectively, and the R2 was 0.9710. Among the three models, the diffusion model of spherical particle fitted R2 was the highest, which could explain more than 97% of the variance in experimental data.
Table 3 shows the errors between the predicted values and the actual values of different models. The prediction values of the diffusion model were closest to the actual values.
The diffusion model of spherical particle assumed that the major resistance of mass transfer in the distillation process of volatile oil was the diffusion of volatile oil from the inside to the surface of the medicinal materials. The mass transfer rate of volatile oil from the surface of the medicinal materials to the extraction solution and the mass transfer rate in the air inside the extractor were all much faster. These results indicate that a rapid distillation of volatile oil can probably be realized by lowering the average particle size of Angelicae sinensis Radix.

3.3. Optimization of Distillation Parameters of Volatile Oil

3.3.1. Data Processing and Model Fitting

By modeling the oil yield data, the multivariate binomial regression model obtained was as follows: Y = 0.710 + 0.026A + 0.010B − 0.026C + 2.50 × 10−3AB − 5.00 × 10−3AC − 7.50 × 10−3BC − 0.068A2 − 0.030B2 − 0.082C2. The model determination coefficient R2 was 0. 9917. The variance analysis of each item in the model is shown in Table 4. The F value of the model was 93.13, and p < 0.0001, which showed that the model was extremely significant. The model R2 exceeded 0.99, which showed that the model could well explain the data variation in the experiment. A (p < 0.0001) and C (p < 0.0001) were extremely significant, and B (p < 0.05) was significant. Among the quadratic terms, A2 (p < 0.0001), B2 (p = 0. 0003), and C2 (p < 0.0001) were all extremely significant terms, which showed that the influence of the three parameters on the oil yield was nonlinear. The results were consistent with the previous single factor experimental results. The p values of AB, BC, and AC were all greater than 0.10, which showed that the interaction between the three factors was small. Figure 6 shows the Pareto chart. It can be concluded that all the linear terms and quadratic terms were significant.

3.3.2. Response Surface Diagram and Contour Diagram

The response surface diagram and contour diagram of the distillation process of the volatile oil from Angelicae sinensis Radix are shown in Figure 7. The figure reflects the effects of the distillation time, soaking time, and liquid–material ratio on the oil yield. When the liquid–material ratio is fixed, the oil yield increases first and then decreases gradually with the increase in distillation time. When the distillation time is too long, some volatile components in the volatile oil may volatilize and be lost. When the soaking time is fixed, the volatile oil yield increases first and then decreases gently with the increase in the distillation time, and increases first and then decreases with the increase in the liquid–material ratio. When the distillation time is fixed, the oil yield increases first and then stabilizes with the increase in soaking time, and increases first and then decreases with the increase in the liquid–material ratio.

3.3.3. Design Space Calculation and Verification

The acceptable lower limit of volatile oil yield was set at 0.65%, and the lowest probability of target attainment was set at 0.80. The design space calculated by the probability-based method is shown in Figure 7. Figure 7a is a three-dimensional design space diagram, and Figure 8b is a two-dimensional space diagram after fixing the liquid–material ratio at 9.3:1. It can be seen from the figure that the design space was close to the shape of a football.
The design space was verified with a selected condition of 8.4 h of distillation time, 3.2 h of soaking time, and 9.3:1 (mL·g−1) of liquid–material ratio. Under this condition, the probability for target attainment was 1.0 and the predicted volatile oil yield was 0.715%. Three parallel validation tests were carried out, and the oil yield of Angelicae sinensis Radix volatile oil was 0.711%, 0.709%, and 0.712%, respectively. The average oil yield was 0.711% and the relative standard deviation was 0.21%. The density of the collected volatile oil was 1.0075 g/cm3 at room temperature. The mathematical model was developed according to the results of the Box–Behnken designed experiments. The measured values were close to the predicted values, indicating that the prediction of the model was accurate. The volatile oil yield was higher than the preset standard, which shows that the design space was reliable.

3.4. Qualitative Analysis of Chemical Constituents of Volatile Oil

The GC-MS analysis of the volatile oil from Angelicae sinensis Radix was carried out, and the total ion flow diagram is shown in Figure 9. The results of the mass spectrometry are shown in Table 5. Twenty chemical constituents were identified, accounting for 99.426% of the total volatile oil. Among them, Z-ligustilide had the highest relative content, and its relative content reached 85.385%. Its structural formula can be found in SciFinder [21]. Other components with relatively high contents were α-pinene, β-ocimene, 2-methyl−1,3-benzoxazole, butylidenephthalide, and E-ligustilide.

4. Discussion

Li Tao et al. [22] found that the content of volatile oil from wild Angelicae sinensis Radix (1.14%) was more than twice as high than that from the artificially cultivated Angelicae sinensis Radix (0.4%). Yan Hui et al. [23] found that longer sunshine was not conducive to the increase in volatile components that mainly consisted of ligustilide. Lin Haiming [24] found that the content of the volatile oil and alcohol-soluble extract decreased gradually with the increase in the drying temperature. Ji Peng et al. [25] determined and analyzed the chemical components of volatile oil in raw Angelicae sinensis Radix and different processed products of Angelicae sinensis Radix by GC-MS. It was found that different processing methods would affect the total amount of volatile oil in Angelicae sinensis Radix and the contents of ligustilide and butenyl phthalolactone in the volatile oil [25]. Li Runhong et al. [26] studied the difference in the volatile oil composition and content in Angelicae sinensis Radix from different producing areas. The results showed that the highest volatile oil content of wild Angelicae sinensis Radix in Linzhi (Tibet Province) was 2-hydrobutenyl phthalide (70.184%), followed by n-butenyl phthalide (9.288%) [26]. In the volatile oil of wild Angelicae sinensis Radix in Pingwu (Sichuan Province), the content of 2-hydrobutenyl phthalide (92.551%) was the highest, followed by butenyl phthalide (3.037%) [26]. These results were different compared with those reported in this study and most of the literature, in which ligustilide was found to be the most abundant constituent. Generally speaking, the composition and content of volatile oil in Angelicae sinensis Radix are affected by the growth environment, harvest time, drying method, processing method, extraction technology, etc. [14,27]. Therefore, the composition and content of volatile oil in Angelicae sinensis Radix may have significant differences.
The results of some published works [28] have shown that the yield of volatile oil obtained from Angelica sinensis by steam distillation was low (about 0.3–0.5%). In industry, the yield is even lower. Sometimes, only the volatile oil aqueous solution can be obtained. In this work, an improved steam distillation with enhanced cooling was used. Compared with that of using a traditional volatile oil extractor, the collected volatile oil could be increased by 30–58%. The collected volatile oil in Angelica sinensis was different to that in the literature [22,29,30], which may be attributed to germplasm resources, harvesting time, planting altitude, processing technology, and so on. However, the efficiency of steam distillation was not high. More volatile oil may be obtained when using other techniques such as supercritical fluid extraction [31].

5. Conclusions

In this work, an improved volatile oil extractor with better cooling was applied to collect the volatile oil from Angelicae sinensis Radix by steam distillation. The distillation process was optimized according to the concept of quality-by-design. The soaking time, distillation time, and liquid–material ratio were determined as CPPs. It was observed that the volatile oil yield increased with the increase in the distillation time. Furthermore, three models were used to fit the distillation kinetics data of the volatile oil. The fitting effect of the diffusion model of spherical particle was best, and the R2 exceeded 0.97. This indicates that a more rapid distillation process can probably be realized when smaller Angelicae sinensis Radix can be used. The volatile oil yield increased first and then decreased with the increase in the soaking time and liquid–material ratio. Then, the distillation time, soaking time, and liquid–material ratio were investigated with the Box–Behnken designed experiments. The R2 of the established second-order polynomial model exceeded 0.99. Then, the design space for distilling the volatile oil from Angelicae sinensis Radix was obtained with a probability-based method. The design space was verified at a condition shown as follows: Distillation time: 8.41 h, soaking time: 3.20 h, liquid–material ratio: 9.3:1 (mL·g−1). The probability of the target attainment was 1 at this condition. The actual volatile oil yield was 0.711%, which was close to the predicted value of 0.715% and higher than the preset standard, indicating that the model could predict accurately and the design space was reliable. The volatile oil of Angelicae sinensis Radix was also analyzed by GC-MS, and 21 chemical constituents were found, accounting for 99.426% of the total volatile oil. The content of Z-ligustilide in the volatile oil was the highest, and the relative content reached 85.38%.

Author Contributions

Conceptualization, X.G.; Investigation, N.W., J.L. and Z.W.; data curation, N.W., J.L., X.C. and Q.Z.; writing—original draft preparation, N.W., X.C., J.L. and X.G.; writing—review and editing, Q.Z. and X.G.; supervision, Q.Z. and X.G.; funding acquisition, X.G. and Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (82060720); the Jiangxi University of Traditional Chinese Medicine 1050 Youth Talent Project; the Open Fund Project of Key Laboratory of Modern Preparation of TCM; the Ministry of Education; the Jiangxi University of Chinese Medicine (zdsys-202110); and the Chinese Medicine Innovation Team and Talent Support Program of The National Administration of Traditional Chinese Medicine (ZYYCXTD-D-202002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Acknowledgments

The authors are grateful to Wanying Wang for her help during the manuscript preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, Z.; Wu, G.; Sun, M.; Du, L.; Ren, Y. Research progress of ligustilide in Angelica sinensis. J. Gansu Univ. Chin. Med. 2018, 35, 102–105. [Google Scholar]
  2. Zhao, X.; Wang, D.; Zhao, D.; Ji, H.; Pei, Y.; Bai, J. Isolation and identification of the chemical constituents from roots of Angelica sinenses (Oliv.) Diels. J. Shenyang Pharm. Univ. 2013, 30, 182–185. [Google Scholar]
  3. Wu, H.; Hua, Y.; Guo, Y.; Wei, Y. Extraction and Analysis of Volatile Oil in Different Parts of Radix Angelicae Sinensis from Min County in Gansu. Nat. Prod. Res. Dev. 2012, 24, 1225–1229. [Google Scholar]
  4. Sun, M.; Ma, Q.; Liu, F.; Ren, Y. New Progress in the Study of Volatile Oil from Angelica sinensis. Abstr. World Curr. Med. Inf. 2019, 19, 56–58. [Google Scholar]
  5. Cui, F.; Feng, L.; Hu, J. Factors affecting stability of z-ligustilide in the volatile oil of radix angelicae sinensis and ligusticum chuanxiong and its stability prediction. Drug Dev Ind Pharm. 2006, 32, 747–755. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, Y.; Zhang, X.; Zou, J.; Jia, Y.; Wang, C.; Shi, Y.; Guo, D. Network analysis of the pharmacological mechanism of Angelica sinensis volatile oil in the treatment of hypertension. Nat. Prod. Res. Dev. 2021, 33, 657–666. [Google Scholar]
  7. Liu, X.; Huang, X.; Zhang, X.; Wang, L.; Zhu, L.; Xu, J.; Chen, Q.; Yang, M.; Wang, F. Q-Marker prediction of volatile oil of Angelicae Sinensis Radix based on GC-MS analysis combined with network pharmacology. Chin. Tradit. Herb. Drugs 2021, 52, 2696–2706. [Google Scholar]
  8. Zhu, L.; Luo, J.; Song, R.; Wang, L.; Zhang, A.; Zang, K. Effect and mechanism of volatile oil of angelica on apoptosis and autophagy in human colorectal cancer SW480 cells. Chin. J. Clin. Pharmacol. 2021, 37, 3253–3256. [Google Scholar]
  9. Kusuma, H.S.; Mahfud, M. Chemical composition of essential oil of Indonesia sandalwood extracted by microwave-assisted hydrodistillation. AIP Conf. Proc. 2016, 1755, 50001. [Google Scholar]
  10. Kusuma, H.S.; Mahfud, M. Box-Behnken design for investigation of microwave-assisted extraction of patchouli oil. AIP Conf. Proc. 2015, 1699, 50014. [Google Scholar]
  11. Priyadarshi, S.; Kashyap, P.; Gadhave, R.K.; Jindal, N. Effect of ultrasound-assisted hydrodistillation on extraction kinetics, chemical composition, and antimicrobial activity of Citrus jambhiri peel essential oil. J. Food Process Eng. 2021, 44, e13904. [Google Scholar] [CrossRef]
  12. Lan, Z.; Wang, L.; Li, Q.; Wang, S.; Meng, J. Analysis of volatile oil components of different species of Curcumae rhizoma based on GC-MS and chemometrics. China J. Chin. Mater. Med. 2021, 46, 3614–3624. [Google Scholar]
  13. Zou, J.; Zhang, X.; Shi, Y.; Guo, D.; Cheng, J.; Cui, C.; Tai, J.; Liang, Y.; Wang, Y.; Wang, M. Kinetic study of extraction of volatile components from turmeric by steam distillation. China J. Tradit. Chin. Med. Pharm. 2020, 35, 1175–1180. [Google Scholar]
  14. Wu, Z.; Xie, L.; Li, Y.; Wang, Y.; Wang, X.; Wan, N.; Huang, X.; Zhang, X.; Yang, M. A novel application of the vacuum distillation technology in extracting Origanum vulgare L. essential oils. Ind. Crop. Prod. 2019, 139, 111516. [Google Scholar] [CrossRef]
  15. Gong, X.; Chen, T.; Qu, H. Research advances in secondary development of Chinese patentmedicines based on quality by design concept. China J. Chin. Mater. Med. 2017, 42, 1031–1036. [Google Scholar]
  16. Tai, Y.; Qu, H.; Gong, X. Design Space Calculation and Continuous Improvement Considering a Noise Parameter: A Case Study of Ethanol Precipitation Process Optimization for Carthami Flos Extract. Separations 2021, 8, 74. [Google Scholar] [CrossRef]
  17. El-Shamy, A.M.; El-Hadek, M.A.; Nassef, A.E.; El-Bindary, R.A. Box-Behnken design to enhance the corrosion resistance of high strength steel alloy in 3.5 wt.% NaCl solution. Moroc. J. Chem. 2020, 8, 788–800. [Google Scholar]
  18. Kusuma, H.S.; Sudrajat, R.G.M.; Febrilliant, D.; Susanto, D.F.; Gala, S.; Mahfud, M. Response Surface Methodology (RSM) Modeling of Microwave-Assisted Extraction of Natural Dye from Swietenia mahagony: A comparation Between Box-Behnken and Central Composite Design Method. AIP Conf. Proc. 2015, 1699, 50009. [Google Scholar]
  19. Zhao, M.; Yang, S.; Sun, Y.; Li, F.; Zhang, S.; Jiao, J. Gas Chromatography-mass Spectrometry Analysis of Volatile Oil of Angelica sinensis Root Dry Chemical Composition. Chem. World 2018, 59, 231–234. [Google Scholar]
  20. Shao, J.; Qu, H.; Gong, X. Comparison of two algorithms for development of design space-overlapping method and probability-based method. China J. Chin. Mater. Med. 2018, 43, 2074–2080. [Google Scholar]
  21. Zhang, Y.; Zhang, Y.; Han, Y.; Tian, Y.; Wu, P.; Xin, A.; Wei, X.; Shi, Y.; Zhang, Z.; Su, G.; et al. Pharmacokinetics, tissue distribution, and safety evaluation of a ligustilide derivative (LIGc). J. Pharm. Biomed. Anal. 2020, 182, 113140. [Google Scholar] [CrossRef]
  22. Li, T.; He, X. Analysis of essential oil from the roots of Angelica sinensis by GC-MS. West China J. Pharm. Sci. 2015, 30, 249–250. [Google Scholar]
  23. Yan, H.; Zhang, X.-B.; Zhu, S.-D.; Qian, D.-W.; Guo, L.-P.; Huang, L.-Q.; Duan, J.-A. Production regionalization study of Chinese angelica based on MaxEnt model. China J. Chin. Mater. Med. 2016, 41, 3139–3147. [Google Scholar]
  24. Liu, P.; Chen, J.; Zhou, B.; Xu, Y.; Qian, D.W.; Duan, J.A. Analysis of variation of coumarin and volatile compounds in Angelica Dahuricae radix in different drying methods and conditions. Zhongguo Zhong Yao Za Zhi 2014, 39, 2653–2659. [Google Scholar]
  25. Ji, P.; Hua, Y.; Xue, W.; Wu, H.; Guo, Y.; Wei, Y. Extraction and Composition Analysis of Essential Oil from Raw Radix Angelicae Sinensis and Its Different Processed Products. Nat. Prod. R D. 2012, 24, 1230–1234, 1238. [Google Scholar]
  26. Li, R. Analysis of Volatile Oil Ingredients in Wild Angelica sinensis from Different Origin. Chin. J. Ethnomed. Ethnopharm. 2019, 28, 41–44. [Google Scholar]
  27. Li, Z.; Cui, J.; Fu, Z.; Mu, J. Supercritical CO2 Fluid Extraction and Chemical Composition Analysis of Volatile Oil from Atractylodis macrocephalae rhizoma. Food Drug 2019, 21, 269–273. [Google Scholar]
  28. Wang, T.; Cheng, Z. Research on Angelica Naphtha Components form Different Regions. Pharm. Biotechnol. 2013, 20, 535–537. [Google Scholar]
  29. Wang, X.; Zhao, Z.; Zhang, Y.; Shi, X.; Cheng, F.; Guo, M.; Wen, X. Comparative analysis of volatile compounds in Angelica Sinensis by different acquisition methods with GC-MS technique. J. Tradit. Chin. Vet. Med. 2018, 37, 61–65. [Google Scholar]
  30. Li, D.; Ma, X.; Song, P.; Zhao, J.; Ding, Y. Analysis of Three Methods for Extracting Volatile Constituent from DangGui by Combination of Gas Chromatography and Mass Spectrometry. West. J. Tradit. Chin. Med. 2013, 26, 15–18. [Google Scholar]
  31. Yue, H.; Jia, H.; Wang, J. Study on the Extraction Technology of the Volatile Oil in Angelica sinensis (Oliv.) Diels Based on Supercritical CO2 Extraction. J. Anhui Agric. Sci. 2010, 1, 47–50. [Google Scholar]
Figure 1. Homemade essential oil extractor.
Figure 1. Homemade essential oil extractor.
Separations 09 00137 g001
Figure 2. The Ishikawa diagram of the steam distillation process for volatile oil distillation.
Figure 2. The Ishikawa diagram of the steam distillation process for volatile oil distillation.
Separations 09 00137 g002
Figure 3. The effect of distillation time on the volatile oil yield (n = 3).
Figure 3. The effect of distillation time on the volatile oil yield (n = 3).
Separations 09 00137 g003
Figure 4. The effect of soaking time on the volatile oil yield (n = 3).
Figure 4. The effect of soaking time on the volatile oil yield (n = 3).
Separations 09 00137 g004
Figure 5. The effect of solid–liquid ratio on the volatile oil yield (n = 3).
Figure 5. The effect of solid–liquid ratio on the volatile oil yield (n = 3).
Separations 09 00137 g005
Figure 6. The Pareto chart of the parameters.
Figure 6. The Pareto chart of the parameters.
Separations 09 00137 g006
Figure 7. The 3D response surface and 2D contour map for the distillation process of volatile oil.
Figure 7. The 3D response surface and 2D contour map for the distillation process of volatile oil.
Separations 09 00137 g007
Figure 8. The design space. (a) Panorama of the design space. (b) The design space when the liquid–material ratio = 9.3:1 (mL·g−1). The red star in the figure represents the condition of the verification experiments, and the color bar indicates the probability of the target attainment.
Figure 8. The design space. (a) Panorama of the design space. (b) The design space when the liquid–material ratio = 9.3:1 (mL·g−1). The red star in the figure represents the condition of the verification experiments, and the color bar indicates the probability of the target attainment.
Separations 09 00137 g008
Figure 9. The total ion chromatogram of the angelica volatile oil.
Figure 9. The total ion chromatogram of the angelica volatile oil.
Separations 09 00137 g009
Table 1. Factors and levels of the Box–Behnken design.
Table 1. Factors and levels of the Box–Behnken design.
Factor Level
Low (−1)Medium (0)High (1)
A: Distillation time (h)6810
B: Soaking time (h)234
C: Liquid–material ratio (mL·g−1)6:110:114:1
Table 2. The Box–Behnken designed experiments and results.
Table 2. The Box–Behnken designed experiments and results.
Serial NumberA (Distillation Time/h)B (Soaking Time/h)C [Liquid-Material Ratio (mL·g−1)]Y (Volatile Oil Yield/%)
1−1−100.58
21−100.63
3−1100.59
41100.65
5−10−10.55
610−10.61
7−1010.52
81010.56
90−1−10.61
1001−10.65
110−110.56
120110.57
130 (8 h)0 (3 h)0 (10:1)0.70
140 (8 h)0 (3 h)0 (10:1)0.71
150 (8 h)0 (3 h)0 (10:1)0.71
160 (8 h)0 (3 h)0 (10:1)0.72
170 (8 h)0 (3 h)0 (10:1)0.71
Table 3. The error results of the different functions.
Table 3. The error results of the different functions.
ModelRMSE
(10−2 mL/g)
SSE
(10−4 mL2/g2)
MSE
(10−4 mL2/g2)
MAE
(10−2 mL/g)
First-order kinetic model0.0320.00510.00100.031
Peleg’s model0.0190.00190.000370.019
Diffusion model of spherical particle0.0130.000820.000160.010
Table 4. The results of the analysis of variance.
Table 4. The results of the analysis of variance.
Variance SourceSum of SquareDegree of FreedomMean SquareF Valuep Value
Model0.06997.65 × 10−393.13<0.0001
A5.513 × 10−315.513 × 10−367.11<0.0001
B8 × 10−418 × 10−49.740.0168
C5.512 × 10−315.512 × 10−367.11<0.0001
AB2.5 × 10−512.5 × 10−50.30.5983
AC1 × 10−411 × 10−41.220.3064
BC2.25 × 10−412.25 × 10−42.740.1419
A220.01910.019233.55<0.0001
B223.789 × 10−313.789 × 10−346.130.0003
C220.02910.029348.88<0.0001
Residual error5.75 × 10−478.214 × 10−5
Misfit term3.75 × 10−431.25 × 10−42.50.1985
Pure error2 × 10−445 × 10−5
Total error0.06916
Table 5. The tentative identification results of the chemical composition of angelica volatile oil.
Table 5. The tentative identification results of the chemical composition of angelica volatile oil.
Serial NumberT/minCompoundChemical FormulaRelative Content/%m/z Value of the Fragment with Maximum AbundanceAbundance
13.85α-Pinene
2,6,6-trimethylbicyclo [3.1.1]hept-2-ene
C10H161.2193.1239,680
25.55β-Ocimene
3,7-dimethylocta−1,3,6-triene
C10H161.6793.1298,432
38.846-butylcyclohepta−1,4-dieneC11H180.5379.161,936
412.68Diamyl Ketone
undecan-6-one
C11H22O0.2171.122,096
514.334-ethenyl-2-methoxyphenolC9H10O20.30150.132,720
615.95Duraldehyde
2,4,5-trimethylbenzaldehyde
C10H12O0.46147.134,272
718.24β-Cedrene
(1S,2R,5S,7S)-2,6,6-trimethyl-8-methylidenetricyclo [5.3.1.01,5]undecane
C15H240.27161.27463
819.364-(1,2-dimethylcyclopent-2-en−1-yl)butan-2-oneC11H18O0.2593.112,845
919.88Sesquichamene
2,4a,8,8-tetramethyl−1,1a,4,5,6,7-hexahydrocyclopropa[j]naphthalene
C15H240.0293.14696
1021.45(1R,2R)−1-ethenyl−1-methyl−4-propan-2-ylidene-2-prop−1-en-2-ylcyclohexaneC15H240.75121.150,288
1121.90β-Bisabolene
(4S)−1-methyl−4-(6-methylhepta−1,5-dien-2-yl)cyclohexene
C15H240.1769.19638
1223.77Alloaromadendrene
1,1,7-trimethyl−4-methylidene-2,3,4a,5,6,7,7a,7b-octahydro−1aH-cyclopropa[e]azulene
C15H240.11121.14412
1324.57Spathulenol
(7S)−1,1,7-trimethyl−4-methylidene−1a,2,3,4a,5,6,7a,7b-octahydrocyclopropa[h]azulen-7-ol
C15H24O0.9691.122,384
1427.362-methyl−1,3-benzoxazoleC8H7NO1.06133.0240,832
1527.521-oxido−4-[2-(1-oxidopyridin−1-ium−4-yl)ethyl]pyridin−1-iumC12H12N2O20.71108.170,000
1627.75Cyclopentadiene
2,5,5-trimethylcyclopenta−1,3-diene
C8H120.09159.0373,632
1728.12Butylidenephthalide
(3E)-3-butylidene-2-benzofuran−1-one
C12H12O22.79159.056,168
1830.06Senkyunolide A
(3S)-3-butyl−4,5-dihydro-3H-2-benzofuran−1-one
C12H16O20.36107.127,088
1930.88Z-ligustilide
(3Z)-3-butylidene−4,5-dihydro-2-benzofuran−1-one
C12H14O285.4161.11,526,784
2032.34E-ligustilide
(3E)-3-butylidene−4,5-dihydro-2-benzofuran−1-one
C12H14O22.04161.186,504
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wan, N.; Lan, J.; Wu, Z.; Chen, X.; Zheng, Q.; Gong, X. Optimization of Steam Distillation Process and Chemical Constituents of Volatile Oil from Angelicaesinensis Radix. Separations 2022, 9, 137. https://doi.org/10.3390/separations9060137

AMA Style

Wan N, Lan J, Wu Z, Chen X, Zheng Q, Gong X. Optimization of Steam Distillation Process and Chemical Constituents of Volatile Oil from Angelicaesinensis Radix. Separations. 2022; 9(6):137. https://doi.org/10.3390/separations9060137

Chicago/Turabian Style

Wan, Na, Jing Lan, Zhenfeng Wu, Xinying Chen, Qin Zheng, and Xingchu Gong. 2022. "Optimization of Steam Distillation Process and Chemical Constituents of Volatile Oil from Angelicaesinensis Radix" Separations 9, no. 6: 137. https://doi.org/10.3390/separations9060137

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop