Rapid Determination of Pachymic Acid Content by Near-Infrared Spectroscopy

National R & D Center for Edible Fungus Processing Technology, Henan University, Kaifeng 475004, China Functional Food Engineering Technology Research Center of Henan Province, Kaifeng 475004, China Huaihe Hospital, Henan University, Kaifeng 475004, China Joint International Research Laboratory of Food & Medicine Resource Function, Kaifeng 475004, Henan, China Kaifeng Key Laboratory of Functional Components in Health Food, Kaifeng 475004, Henan, China


Introduction
Poria cocos, one important medicinal fungus, has been used in China for a long time [1,2]. Pachymic acid, a rich component in Poria, is the key active ingredient of Poria, with the effects of anti-inflammatory [3], antitumor [4], and hypoglycemic [5]. Usually, pachymic acid was used as the mark component to reflect the quality of Poria [6]. Nowadays, the common determination methods for pachymic acid include the colorimetric method and the HPLC method. For its high accuracy and stability, HPLC is usually taken as a reference method, but its analysis process is labor-intensive and time-consuming, which has hampered its wide applicability in the market. Moreover, the 2020 edition of the "Chinese Pharmacopoeia" has not included one quantitative determination method for this component. erefore, a method with high efficiency and throughput for the determination of pachymic acid content in Poria is necessary.
NIR spectroscopy is one rapid determination method widely used in the research field of traditional Chinese medicine recently for the advantages of fast speed [7,8], high efficiency, high sensitivity, simple sample processing, and convenient instrument operation [9][10][11]. NIR spectroscopy has been widely used in food [12,13], medicine [14][15][16], and other industries [17,18]. Gao Hongbin et al. [19] used vector normalization and first derivative to process the NIR spectrum data of Glycyrrhizae Radix et Rhizoma Yinpian and established the corresponding rapid NIR analysis model for the determination of glycyrrhizic acid and glycyrrhizin with the PLS method. e correlation coefficients of the models were 0.981 and 0.919, and the mean square errors of cross-validation were 0.184 μg/mL and 0.144 μg/mL, respectively. Ying Liang [20] used the PLS method combined with cross-validation methods with regression analysis and MSC + SD data processing methods to develop a NIR analysis model for rapid determination of naringin content. e model has good indexes and can rapidly evaluate the content of naringin in Exocarpium Citri Grandis.
According to references, it can be concluded that the NIR spectrum majorly reflects the doubling-and combination-frequency absorption of the movement of some characteristic groups (such as C-H, O-H, C�O, -CH 2 , and -CH 3 ), and these groups have a fixed vibration frequency which can excite to generate resonance when irradiated by infrared rays [21]. e quantity and connection of these groups of pachymic acid could form a specific NIR spectrum in the instrument, and the amount of energy in the nearinfrared light absorbed by the pachymic acid, which can be mined in the NIR spectrum by the Poria sample, could reflect the quantity of pachymic acid. In this study, NIR were used with PLS regression to establish a rapid quantitative analysis model to detect pachymic acid, which was expected to provide technical reference for quality control and market management of Poria.

Materials and Reagents.
A total of 108 batches of Poria samples were collected in small cube shapes from the market in 2019. e collected regions covered the main production provinces and the Dabie Mountains of China, including Hunan, Anhui, Yunnan, Hubei, Sichuan, Guizhou, Chongqing, Fujian, Guangdong, Henan, Jiangxi, Zhejiang, and Guangxi Provinces.
During this research, acetonitrile (HPLC grade) and methanol (HPLC grade) were purchased from Longxi Science Co., Ltd., acetic acid (HPLC grade) was purchased from Tianjin Fuyu Fine Chemical Co., Ltd., purified water were from Hangzhou Wahaha Baili Food Co., Ltd., and pachymic acid (Batch No. 18062101) was purchased from Chengdu Pufei Biotechnology Co., Ltd.

Preparation of the Reference Solution.
Pachymic acid was accurately weighed and dissolved in methanol to prepare a 0.49 mg/mL reference substance solution, which was filtered through a 0.22 μm microporous membrane and reserved for testing.

Preparation of the Test Solution.
Each sample was crushed into powder of 80 mesh fineness and placed in a 10 mL conical flask. e powder was dissolved with methanol to make a 0.125 g/mL solution and was disposed of by ultrasonic (1200 W) for 40 min. e supernatant was taken and reserved for testing after filtering through a 0.22 μm microporous filter membrane.

e HPLC Conditions.
e pachymic acid content of the reference solution and test solution was determined by HPLC with the chromatographic column being ermo BDS Hypersil C18 (250 × 4.6 mm, 5 μm). e LC time program is shown in Table 1. e flow rate was 1.4 mL/min, the column temperature was 30°C, the detection wavelength was 210 nm, and the sample size was 20 μL. Based on this condition, the corresponding chromatogram is shown in Figure 1. e peak with the retention time of 26 min was the pachymic acid.

Evaluation of HPLC Methods. Linear Relation Investigation.
In this part, the chromatograph and corresponding peak areas of 5 injections of reference solutions with different concentrations (0.0245, 0.049, 0.49, 0.98, and 1.47 mg/mL) were recorded. e linear regression was established with the content of reference solutions (x) and corresponding chromatographic peak areas (y) and shown as the equation y � 3 × 10 8 x +4162.2. e linear regression coefficient of this equation was 0.9999, indicating that the peak area of pachymic acid had a good linear correlation with the content in the range of 0.0245∼1.4700 mg/mL. Reproducibility. A total of 6 testing solutions prepared from the same sample were tested separately. e RSD of the test was 0.48%, indicating that this method had a good repeatability.
Precision. In this part, one reference solution was injected 6 times continuously for the precision testing of this method. e RSD was 0.31%, which indicated that the instrument and method could determine the content of pachymic acid accurately.
Stability. Samples prepared according to 2.3.2 were injected at 0, 4, 8, 16, 20, and 24 h, respectively, for testing the stability of the testing solution. e results showed that the RSD of this method was 0.28%, indicating that pachymic acid in the test solution was basically stable within 24 h.
Recovery. In this part, 6 samples with a known content of pachymic acid were added with a certain amount of standard pachymic acid, respectively. e average recovery rate of this

NIR Spectrum Acquisition and Analysis.
e Poria samples were crushed into powder, sieved through 80 mesh, and evenly spread on the bottom of the test dish (its diameter is 5 cm). e resolution was set to 8 cm −1 , scanning was 64 times per average spectrum within the spectral range of 4000-10000 cm −1 , the scanning type was interleaved mode, and the internal reference was automatically used as the background. Based on the determination parameters, the spectrum of each sample could be obtained (Figure 2).

Establishment of the NIR Model.
PLS, artificial neural networks (ANN), support vector machines (SVM), and other nonlinear correction methods are used to establish the model with their own advantages and disadvantages [22], and in this research, PLS was adopted to build the NIR analysis model with the software Spectrum Quant. In addition, 4 parameters including r 2 c , RMSEC, RMSEP, and r 2 v were used to evaluate the established model. r 2 c is used to evaluate the degree of linearity between the predicted and measured values of the samples. e closer r 2 c is to 1, the better the regression or prediction result is. RMSEC is mainly used to evaluate the feasibility and prediction ability of modeling methods. RMSEP is mainly used to evaluate the prediction ability of the model for external samples. e smaller the value is, the better the prediction ability of the model for external samples is. r 2 v is used to evaluate the linearity between the predicted value and the measured value of the model for the unknown sample during the external validation of the model [21]. Its value is approximately close to 1, indicating that the model predicts better results. Spectrum Quant was used for modeling in this study.  Figure 3 is a scatter plot about the pachymic acid content of 108 batches of Poria. A total of 91 Poria samples were randomly selected from 108 samples according to the sequence of each sample by using the random function in Excel for composing the calibration set, and the remaining samples were included in the validation set. At the same time, the maximum and minimum values of pachymic acid content from the calibration set should cover the validation set. e information characteristics of the samples from the two sets are shown in Table 2.

Selection of Spectral Interval of NIR Analysis Model.
e absorption effect of substances varies in different NIR spectra, and it is necessary to determine the proper spectral range for the PA determination [23]. In this study, 4 potential spectral intervals were adopted, and the evaluation results of the corresponding analysis model are shown in Table 3.
According to the results, the spectral range is 4000-10000 cm −1 , the full spectral range with the largest r 2 c , the smallest RMSEC and RMSEP, as well as the proper ratio of latter two indices were selected for modeling.

Selection of Spectral Pretreatment Methods for NIR Analysis Models.
During the data analysis of NIR spectra, some random factors, such as sample particle size and light interference, could cause the baseline drift, affecting the reliability and stability of model building [24][25][26]. Single and combination data preprocessing methods were used for the modeling and the results are listed in Table 4. MSC mainly eliminates the scattering effect caused by uneven sample size and particle size. SNV is mainly used to eliminate the influence of particle size, surface scattering light, and optical path change on the near-infrared spectrum. Derivatives can effectively eliminate baseline and other background interference [21]. It can be seen from Table 4 that SD, SNV + SD, and MSC + SD are the better preprocessing methods for the spectral information of the NIR analysis model. e r 2 c of the three pretreatment methods is    greater than 0.9, and RMSEC and RMSEP meet the basic requirements of model building. Among them, the values of the three indexes (r 2 c , RMSEC, and RMSEP) of SNV + SD and MSC + SD pretreatment methods have a slight difference, so these three pretreatment methods are taken into consideration.

Determination of PLS Principal Factor Number in the NIR Analysis Model.
When PLS was used to establish the NIR analysis model, the number of PLS principal factors was also a key factor affecting the quality of the model. Either too many or too less principal factors would cause underfitting or overfitting problems, respectively, increasing the prediction error [27]. In this part, the corresponding results of the models with specific data preprocessing method and principle factors are listed in Table 5.
It can be seen from Table 5 that all alternative data preprocessing methods with more principal factors could cause higher r 2 c and smaller RMSEC and RMSEP, and the models with over 9 principal factors would have a suitable r 2 c higher than 0.90. So 9 principal factors were the basic requirements for the modeling.

Validation and Evaluation of the Model
On the basis of the former parts of this research, 6 aforementioned models (No. 1, 2, 4, 5, 7, and 8 in Table 5) with r 2 c over 0.90 were further verified by using validation set samples in this part, and it was found that the maximum r 2 v of model 2 is 0.917. From this, the optimal NIR analysis model (No. 2 model in Table 5) for predicting the content of pachymic acid of Poria was obtained, and its parameters are shown in Table 6. e distribution of the predicted and experimental values was scattered intuitively in a line graph (Figure 4) [28]. At the same time, SPSS statistical software (SPSS Inc., Chicago, USA) was used to analyze the difference between the two values by comparing the ratio of them with 1, and there was no significant difference found between them (P � 0.512).
is result indicated that the NIR analysis model had good prediction accuracy for pachymic acid of Poria.

Conclusion
ere are already many national standards by using NIR for qualitative and quantitative detection now [29,30]. NIR has also been widely recognized in the field of traditional Chinese medicine. It is widely used in the production process of traditional Chinese medicines, including qualitative identification of the types and authenticity of traditional Chinese medicines [31], quantitative determination of effective components, and monitoring of the drug production process [32,33]. In this study, different NIR analysis models for pachymic acid were evaluated with different parameters by r 2 c , RMSEC, RMSEP, and r 2 v values. e results showed that the NIR analysis model with pretreatment method of standard normal variate (SNV) + second derivative (SD), the spectral interval of 4000-10000 cm −1 , and 9 PLS major factors could predict the pachymic acid content in Poria exactly. Besides, there was no significant difference   Journal of Food Quality 5 between the predicted value and the measured value of the prediction set (P � 0.512). According to Williams et al.'s results, this NIR model is more precise than the standard for most applications and almost reaches the standard for quality assurance (r 2 ≥ 0.92) [34]. Hence, many improvements still need to be made to this model in further study, in order to promote its applicability and practicability. e pachymic acid content of Poria was successively determined by NIR combined with stoichiometry, expanding the applicability of NIR. Compared with other analytical methods, the NIR method has the advantages of simple pretreatment, fast analytical speed, high multisample processing capacity, and low cost, which is beneficial to the development promotion of the Poria industry [35].

Data Availability
e data used to support the findings of this study are included within the article.

Ethical Approval
Not applicable.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.