Abstract

Background. Total flavonoids content (TFC) is one of the most important quality indexes of Ginkgo biloba leaf, and it is concerned with total antioxidant activity. Near-infrared spectroscopy (NIR) method has showed its advantages in fast, accurate, qualitative, and quantitative analysis of various components in many quality control researches. In this study, a calibration model was built by partial least squares regression (PLSR) coupling with NIR spectrum to quantitatively analyze the TFC and total antioxidant activity of Ginkgo biloba leaf. Results. During the model establishing, some spectrum pretreatment and outlier diagnosis methods were optimized to establish the final model. The coefficients of determination (R2) for TFC and total antioxidant activity prediction were 0.8863 and 0.8486, respectively; and the root mean square errors of prediction (RMSEP) were 2.203 mg/g and 0.2211 mM/g, respectively. Conclusion. These results showed that NIR method combined with chemometrics is suitable for quantitative analysis of main components and their activities and might be applied to quality control of relevant products.

1. Introduction

Ginkgo biloba L. (Ginkgoaceae) is an ancient tree growing in China for thousands of years. In recent decades, Ginkgo biloba L. occupies a prominent position among the best-selling natural products owing to its reliable and remarkable biological activities [1, 2]. Some studies and clinical trials have observed that Ginkgo biloba L. shows potent actions on cardiovascular system and cerebral vascular activity. Furthermore, due to its antioxidant properties, it has been used in Alzheimer’s patients [35].

The main active components in Ginkgo leaf are flavonoids and terpene trilactones. Because only several kinds of terpene trilactones were found in Ginkgo, the pharmacological effects of them are relatively clear and the corresponding quality evaluation is simply achieved [6, 7], while more than 70 kinds of flavonoids were identified, which associate with various kinds of pharmacological actives [811]. Therefore, more researches have been conducted which focused on flavonoids in Ginkgo as their broad-spectrum of antioxidant and free-radical scavenging activity. Thus, total flavonoids content is often considered as an important quality index of Ginkgo biloba leaf. Determination of total flavonoids contents in Ginkgo biloba leaf and further estimating their antioxidant are of great importance to their qualities [7, 12].

In recent years, flavonoids in Ginkgo biloba L. have received considerable attention in various literatures, especially due to their widely recognized free-radical scavenging activity. TFC is often considered an important quality index of Ginkgo biloba products and samples, while DPPH radical methods are a stable N-centered radical at room temperature that is widely employed to assess the radical scavenging properties of antioxidants, such as total flavonoids. The traditional methods for determining total flavonoids in botanical materials are based on chemical extraction and couple with various analytical techniques, such as HPLC [1315], GC [16], and ultraviolet spectrometry [17]. These methods are precise, and some methods are used as the reference methods for TFC detection. Yet, these methods all needed the sample pretreat or extraction process, which are often time-consuming and destructive. Therefore, a rapid, accurate, and even nondestructive analytical method is needed to identify the TFC and further determine their actives for the quality control of Ginkgo biloba L.

Near-infrared spectroscopy (NIR) has been used for various applications, such as quality estimation and quality control of various food, agriculture, and pharmaceutical products. There are also many researches based on NIR in herbal quality control and main contents analysis, and they show the advantages, including simple sample preparation and rapid and simultaneous analysis of several analytes in a large number of samples. The NIR spectra combining with appropriate mathematical models and pattern recognition techniques can be used to qualitatively and quantitatively determine quality of various products [12, 1820].

Some studies have been published by coupling near-infrared spectroscopy with chemometrics methods to qualitative and quantitative analysis of flavonoids concentrations in various botanical leaves and relevant products, including Ginkgo leaf. Shi et al. determined the TFC in fresh Ginkgo leaves with different colors by using the NIR spectroscopy; furthermore they also analyzed the basic structure of flavonoids and relationship of wavelength regions [12]. Liu et al. published their reviews about the roles of flavonoid and its broad-spectrum free-radical scavenging activities in Ginkgo biloba chemical analysis and quality control [8]. Geng et al. established a quantitative near-infrared diffuse reflectance spectroscopy method for the simultaneous determination of three flavonol aglycones in Ginkgo biloba extracts [21]. Yet, no research has been reported to quantitatively analyze the TFC and their total antioxidant activity in Ginkgo biloba samples, simultaneously.

Based on these reasons, we aimed to establish a calibration model to quantitatively analyze the TFC in Ginkgo biloba leaves and further quantitatively estimate their antioxidant properties. During the model establishing process, some NIR signal pretreat methods were adopted to optimize the calibration model. The feasibility of combining NIR spectroscopy with chemometrics methods to rapid and nondestructive determination of TFC and their antioxidant properties was investigated.

2. Material and Methods

2.1. Chemicals and Materials

1,1-Diphenyl-2-picrylhydrazyl (DPPH) and Trolox were purchased from Sigma-Aldrich Chemical Co. (St. Louis, MO, USA). NaNO2, NaOH, and Al (NO3)3 were analytical grade and acquired from Shanghai Macklin Biochemical Co., Ltd. (Shanghai, China). Rutin was obtained from the National Institution for Food and Drug Control (Beijing, China). Distilled water was filtered by using a Milli-Q water-purification system (Millipore, Bedford, MA, USA). 113 batches of samples were collected from Zhejiang Province. All samples were powdered by a grinder mill after dried and passed through 60-mesh place before analysis.

2.2. Total Flavonoids Content (TFC)

The concentrations of flavonoids were quantified based on a colorimetric assay method [12], with slight modifications. Briefly, rutin was used as a standard to establish calibration linear with function:

0.90-1.00 g samples were weighted, and 10 mL 60% ethanol aqueous was used to extract flavonoids from theses samples with supersonic (KQ-300DE, Kunshan Ultrasonic Equipment Co., China) for 30 min. These samples were further centrifuged at 3000 g. All the supernatant was transferred to 25 mL volumetric flask and then was fixed to 25 mL with 60% ethanol aqueous. 1.5 mL of each extracts and 4.5 mL of distilled water were pipetted into a 25 mL tube and then mixed with 1 mL 5% (w v−1) NaNO2 solutions. After incubation for 6 min, 1 mL of the 10% (w v−1) Al (NO3)3 solutions was added to the mixture. The mixture was kept for 6 min before adding 10 mL 4% (w v−1) NaOH solutions and fixed to 25 mL with 60% ethanol aqueous. Finally, the mixture was reacted for 15 min and the absorbance of the mixture solution was measured with a spectrophotometer (SP-1901, Shanghai Spectrum Instruments Co., China) at 510 nm against a blank containing 5 mL of extraction solvent. Samples were independently analyzed in triplicate times, and the mean of three tests were used and the total flavonoid content was expressed as mg rutin equivalent per g dry weight (DW).

2.3. Determination of the Total Antioxidant Activity

DPPH radical scavenging activity was determined as described by Okawa et al. [22] with a slight modification. Solutions of known Trolox concentration were used for calibration. 2 μL of samples or Trolox was mixed with 250 μL of methanolic DPPH. The homogenate was shaken vigorously and kept in darkness for 30 min. Absorption of the samples was measured on the spectrophotometer at 515 nm. Results were expressed as Trolox equivalent per g of dry weight (mM TE g dried extract−1).

2.4. NIR Spectra Acquisition and Preprocessing

The NIR spectra were measured in a diffuse reflectance mode by Antaris II FT-NIR spectrophotometer (Thermo EIectron Co., USA) equipped with an integrating sphere. The spectra (4000 to 8000 cm−1 were analyzed, and total 4150 points/spectrum) were collected in the log (1/R) mode which was converted by the reflectance value (R). Each sample (0.5 g) was placed in the sample cup, each sample was measured three times, and the mean of the three spectra was used for further statistical analysis. The temperature was kept around at 25°C.

2.5. Multivariate Calibration Methods Establishing

The whole establishing process of calibration model has been described as follow steps.

Firstly, Kennard and Stone algorithm (K-S) [23] was adopted to split samples into calibration dataset (80%) and prediction dataset (20%), respectively. Then, the calibration dataset was used to develop calibration model for TFC and its antioxidant activity by using the partial least squares regression (PLSR). In calibration model, the number of PLS factors were optimized by 10-fold cross-validation method. It was performed as follows: 90% of the calibration dataset samples were used to form the calibration model, and the remaining 10% samples were used to validate this model, and the procedure was repeated by 10 times; the root mean square error of cross-validation (RMSECV) was then calculated as follows:

where n is the number of samples in the calibration set, is the measured result for sample , and is the predicted value of sample . The performances of the optimal model were evaluated according to root mean square error of calibration (RMSEC). RMSEC is calculated as follows:

where is the number of samples in the calibration set, is the measured result for sample, and is the predicted value of sample .

Then, the optimal model which was validated by prediction samples in the prediction dataset RMSEP is calculated as follows:

where is the number of samples in the prediction dataset, is the measured result for sample , and is the predicted value of sample .

Correction coefficient between the predicted value of PLSR model and the measurement value is calculated as follows for both the calibration and prediction set:where n is the number of samples in the calibration or the prediction set and is the mean of measurement value for the calibration or the prediction set.

2.6. Spectral Signal Preprocessing

In the model establishing process, some signal pretreat need be optimized for achieving the best calibration model. In this study, some data preprocessing methods were used to process these NIR signals, such as multiplicative scattering correction (MSC), standard normal transformation (SNV), moving window smoothing, and Savitzky-Golay first derivative or second derivative (S/G 1st/2nd der). The detailed descriptions of these process methods can be found in previous researches [2325]. All of the algorithms were implemented in MATLAB 8.0.1, and all the programs were written by own group (Mathworks).

3. Results and Discussion

3.1. TFC and Antioxidant Activity

The reduction capability of DPPH radical was determined by the decrease in absorbance by reduced DPPH to by plant antioxidants [26]. Therefore, getting systematic knowledge of flavonoids in Ginkgo biloba L. and their activities are highly important for the research and development of this plant. The results of TFC of these samples and their antioxidant activities were listed in Table 1.

3.2. Near-Infrared Spectra

Original NIR spectra of all samples are similar and broad; they consist of many overlapping narrow bands of different vibrational modes, as showed in Figure 1(a), the raw NIR spectra of 113 samples. It can be seen that the intensive spectral peaks are mainly in the region of 4000-8000 cm−1. The multiplicative scattering correction processed spectra of all samples (Figure 1(b)) showed the most intensive band in the spectrum belonging to the vibration of the second overtone of the carbonyl group (5352 cm-1); these are caused by the stretch or deformation vibration of C-H, O-H, and N-H groups, the first two of which are abundant in the flavonoids. Also, it might be caused by the combination of stretching and deformation of the O-H group in water for the spectral peaks intense absorption bands at 6900 cm−1 and 5180 cm−1.

3.3. NIR Calibration Model Establish

In this section, we aimed to establish a reliable and accurate calibration model for TFC and their antioxidant activities quantitative estimation. Thus, some process algorithms were taken into account. The whole calibration model established process in this study contains these steps.

Firstly, the PCA method was adopted to analyze these data for exposing cluster trends in the samples information. Secondly, some anomalous spectra were detected by using the Mahalanobis distance and hat matrix method to detect the outliers. Thirdly, after removing these anomalous spectra, the remaining samples were divided into a calibration set and a prediction set by using the Kennard-Stone (K-S) algorithm. The calibration set was used to optimize the model pretreat processes and establish the calibration model; the prediction set was used as external set to validate model. In the calibration model establish process, some spectra of pretreat methods and variable selective approach were optimized.

3.4. Principle Component Analysis

Principle component analysis (PCA) can be used as a preprocessing step to show samples distribution and visually descript the samples similarities and dissimilarities. In this section, we planned to use PCA to check the cluster of leaves’ NIR spectra. Plotting PCA scores in two or three dimensions provides an effortless way to observe the data distribution. In the PC1-PC2 plot, the first two PCs contain about 95% (PC1: 90.32%, PC2: 5.17%) information of the raw data, and samples are unevenly distributed without obvious cluster which can be found (Figure 2). All samples were located at the whole positions of PCA scores plot. However, we still could find some samples, such as 99, 100, and 113 which are located far from other samples. These samples might be the outlier samples, which can affect the model calibration.

3.5. Anomalous Spectra Detection

In this section, two outlier measure methods were applied to accurately explore the sample information. Techniques based on the Mahalanobis distance (MD) [27] and hat matrix [28] were applied in different fields of chemometrics such as multivariate calibration, pattern recognition, and process control. In the original variable space, the MD considered the correlation in the data, since it is calculated using the inverse of the variance-covariance matrix of the data set. The “Mahalanobis distance” between all the pairs of samples was calculated. As can be seen from Figure 3, 4 samples (2, 5, 99, and 100) were defined as the outliers as the four samples were significant far away other samples.

Furthermore, the hat matrix method was used to estimate the similarities of these samples. We can find that the 4 samples (2, 5, 99, and 100) were also chosen as the outliers (Figure 4). Based on these analysis, the four samples were deleted, and the remaining 108 samples were used for establishing the calibration model.

3.6. Signal Pretreatment and Prediction Model Establish

Kennard-Stone (K-S) algorithm was used to split the dataset into calibration dataset and prediction dataset with split ratio 80%. Thus, 88 samples were used to optimize the calibration model, and remaining 20 samples were used to estimate the established model. In the application of PLS algorithm, it is generally known that the spectral preprocessing methods and the number of PLS factors are critical parameters. Here, their effects on the results are discussed. The optimum number of factors is determined by the lowest root mean square error cross-validation (RMSECV).

Pretreating spectra are a procedure to optimize data and avoid disturbance due to a changing baseline. Common used pretreatments method is averaging, smoothing and normalizing with first and second derivative spectra. The first derivative can eliminate shift errors and the second derivative eliminate tilt errors. Other methods such as multiplicative scatter correction (MSC), Savizaky-Golay method (SG), and standard normal variate (SNV) are also widely used in the NIR spectra. The number of PLS factors included in the model is chosen according to the lowest RMSECV. For RMSECV, a 10-fold cross-validation was performed. Figures 5(a) and 5(b) showed RMSECV plotted versus relevant PLS factors for determining TFC and their antioxidant activity with different spectral preprocessing methods, respectively. Standard normal variate spectral preprocessing method is obviously superior to other methods with lowest RMSECV values. Therefore, standard normal variate (SNV) spectral preprocessing method and corresponding optimized factor were selected to establish the calibration model.

3.7. Calibration Model Validation

The robustness of the method obtained by NIR technology was validated with the 20 prediction samples. The performance of the final PLS model was evaluated in terms of root mean square error of cross-validation (RMSECV), the root mean square error of prediction (RMSEP), and the square of correlation coefficient (R2). As can be seen from Figure 6, the coefficients of determination (R2) for TFC and total antioxidant activity prediction were 0.8863 and 0.8486, respectively, and the root mean square errors of prediction (RMSEP) were 2.203 mg g−1 and 0.2211 mM g−1, respectively.

4. Conclusion

In this study, a method was proposed to quantitatively analyze the TFC and their antioxidant activity of Ginkgo biloba L. by combining NIR spectroscopy coupled with chemometrics methods. The results verified that NIR spectroscopy was a suitable tool for quantification of TFC and their antioxidant activity, simultaneously. Comparing with other analysis methods, the NIR method has its advantages such as being simply pretreated, fast analysis speed, and being nondestructive; these make this approach has the potential of high sample throughput analysis and low costs and widely applied to products’ quality control.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Ling-jia Zhao and Wei Liu have equal contribution to this article.

Acknowledgments

This study was funded by the National Natural Science Foundation of China (81503041), Natural Science Foundation of Zhejiang Province (LY14H270008; LY15H280007), Natural Science Foundation of Hunan Province (2017JJ4045), and Project Foundation of Changsha Science and Technology Bureau (kq1701073).