Novel mathematical model for the assessment of similarity of chromatographic fingerprints of volatile oil from Houttuynia cordata
Jin Zhou1, Qimeng Fan1, Yutian Zhang1, Roxanne Castillo2, Meifeng Xiao1, Hui Liu1, Zhifei Zhu1, Youzhi Liu1, Yantao Yang1, Yiqun Zhou1, Xue Pan1, Fuyuan He1
1 College of Pharmacy, Hunan University of Chinese Medicine; Hunan Key Laboratory of Druggability and Preparation Modification for Traditional Chinese Medicine; Supramolecular Mechanism and Mathematic-Physics Characterization for Chinese Materia Medicine, Changsha, Hunan, China
2 University of California Los Angeles, Los Angeles, USA
|Date of Submission||22-May-2020|
|Date of Decision||06-Jul-2020|
|Date of Acceptance||15-Dec-2020|
|Date of Web Publication||15-Apr-2021|
College of Pharmacy, Hunan University of Chinese Medicine, Changsha, Hunan 410208
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background: The analysis of similarities among fingerprints of Chinese herbal medicines is an important quality control tool to determine the authenticity of the herbal medicines. Objectives: In this study, we aimed to develop a novel mathematical model to analyze the similarity of the chromatographic fingerprints of Houttuynia cordata (HC). Materials and Methods: Total quantum statistical moment similarity (TQSMS) expressions were deduced to evaluate the similarities between two chromatographic fingerprints. The volatile oil samples of HC were analyzed with gas chromatography-mass spectrometry, and the fingerprints were constructed by the area under the peak of the chromatograms. Results: There were nine peaks in common, and a total of 733 chemical constituents observed among 15 batches of samples. The number of peaks in the chromatographic fingerprints of the 15 batches of HC was 49–137, with a relative standard deviation (RSD) of 30.13%. The sum of area under the peak was 1.159 × 107–3.437 × 108 μv × s, with an RSD 174.56%; MCRTT was 9.410–18.602 min, with an RSD of 20.79%; and VCRTT was 37.549–81.504, with an RSD of 23.27%. The volatile oil composition and content of HC showed strong fluctuation. Therefore, its quality control from the variety and content of the components is impractical. Since TQSMS method can characterize the sample similarity, we can quantitate the correct probability of positive and negative conclusions regardless of the population origin of the samples. Conclusion: Our results show that TQSMS can be an additional method that can be used to assess the similarity of two chromatographic fingerprints.
Keywords: Chromatographic fingerprints, gas chromatography-mass spectrometry, herbal medicine, Houttuynia cordata, total quantum statistical moment similarity, volatile oil
|How to cite this article:|
Zhou J, Fan Q, Zhang Y, Castillo R, Xiao M, Liu H, Zhu Z, Liu Y, Yang Y, Zhou Y, Pan X, He F. Novel mathematical model for the assessment of similarity of chromatographic fingerprints of volatile oil from Houttuynia cordata. Phcog Mag 2021;17:154-62
|How to cite this URL:|
Zhou J, Fan Q, Zhang Y, Castillo R, Xiao M, Liu H, Zhu Z, Liu Y, Yang Y, Zhou Y, Pan X, He F. Novel mathematical model for the assessment of similarity of chromatographic fingerprints of volatile oil from Houttuynia cordata. Phcog Mag [serial online] 2021 [cited 2021 Aug 4];17:154-62. Available from: http://www.phcog.com/text.asp?2021/17/73/154/313489
- The volatile oil composition and content of Houttuynia cordata showed strong fluctuation
- The total quantum statistical moment similarity can characterize the sample similarity, and we can quantitate the correct probability of positive and negative conclusions regardless of the origin of the samples.
Abbreviations used: TQSMS: Total quantum statistical moment similarity; TQSM: Total quantum statistical moment; GC-MS: Gas chromatography-mass spectrometry; AUCT: Area under the curve of total quantum; MCRTT: Mean chromatographic retention time of total quantum; VCRTT: Variance of mean chromatographic retention time of total quantum; D: Deviation; pV: Variable probability; 1 − β: Confidence of probability ST: Similarity of total quantum statistical moment; α: Confidence coefficient; HC: Houttuynia cordata; RSD: Relative standard deviation.
| Introduction|| |
Chinese herbal medicines have been successfully used to treat various diseases since thousands of years, yet there is a lack of adequate evidence to substantiate the quality of herbal medicines. Recently, chromatographic fingerprints have been proposed to identify the authenticity of herbal medicines, which provide quality control measures for complex systems with multiple components.,, Currently, many mathematical models have been applied to assess the similarities or differences between various herbal medicines, such as angle cosine, correlation coefficient, fuzzy cusp T distribution, and Euclidean distance. The methods used in the fingerprint similarity analysis usually are divided into the fingerprint characteristic peak response value into discrete data information, adopted to a corresponding characteristic peak of a multidimensional vector method, and calculated to find the similarity, often causing an unstable result.
Although great progress has been made in the field of fingerprint analysis, some problems still need attention., First, under the same conditions of analysis, researchers have obtained variations in peak retention times and heights of a given sample. However, it is noteworthy that the peak areas of a given concentration were relatively stable. Second, the differences among the chromatographic fingerprints are influenced by factors such as analytical technics, geographical factors, and harvest time.,,,, In other words, for one sample of herbal medicine, any change in conditions, analytical or environmental, will lead to differences in their chromatographic fingerprints. Composition and characteristics of herbal medicines vary with their collection time. Unfortunately, existing methods usually do not compare samples with an objective statistic to evaluate the similarities or differences between two fingerprints. Thus, an objective method is needed to evaluate the similarities among fingerprints for quality control and routine authenticity of herbal medicines.
Houttuynia cordata (HC, Yuxingcao in Chinese), the dried aerial part of HC Thunb. (Saururaceae), is one of the best-known natural herbs used since thousands of years in the history of traditional Chinese medicine. Modern pharmacological studies have shown that it shows pharmacological activities, such as antimutagenic, anti-inflammatory, antiviral, antibacterial, antiallergic, antidiabetic, antioxidant, and antiobesity. The primary active components of HC include volatile oils, flavonoids, alkaloids, organic acids, and trace elements. The volatile oils are considered as the major functional components in HC. As is known, the composition of HC changed greatly in the variety and contents with planting environment and harvesting time. However, these existing methods of similarity analysis, including the angle cosine method and correlation coefficient method, cannot accurately reflect the reality of multidimensional vector deviation degree and are more sensitive to the change in peaks with a larger response and are less sensitive to small peak.
To address this issue, we established a novel qualitative and quantitative mathematical model for chromatographic fingerprint analysis, namely, the total quantum statistical moment similarity (TQSMS). TQSMS can be used to characterize the characteristic information of chromatographic fingerprints. Specifically, the area under the curve of total quantum (AUCT) can be used in quantitative analysis, whereas the mean retention time of total quantum (MCRTT) and variance of mean retention time of total quantum (VCRTT) can be used in qualitative analysis. According to the properties of the normal distribution probability density function and the total quantum statistical moment, the TQSMS can then be established and elucidated. To test the model, we analyzed gas chromatographic fingerprints of the volatile oil of HC.
| Materials and Methods|| |
Total quantum statistical moment similarity established for chromatographic fingerprint
Statistical moment methods and theory are powerful techniques for characterizing the chromatographic peaks of fingerprints. In detail, the area of a peak is defined as the zeroth moment, the retention time is defined as the first moment, the variance in retention time is defined as the second moment, whereas the higher moments are associated with the peak shape. However, for a more in-depth analysis, further elaboration of zeroth, first, and second moment of a chromatographic fingerprint is essential.
Zeroth moment of total quanta (AUCT)
Each chromatographic peak in a fingerprint can be considered as a Gaussian curve. Therefore, a complete chromatographic fingerprint can be regarded as the superposition of m characteristic peak response curves. The area under curve (AUC) of the chromatographic fingerprint spectrum (AUCT) was defined as the zeroth moment of total quanta, i.e., all chromatographic peak areas integrated versus time from zero to infinity, and can be calculated by Equation (1).
First moment of total quanta (MCRTT/tT)
The first moment of total quanta is defined as the mean chromatographic retention time (MCRTT). It can be calculated by Equation (2).
Second moment of total quanta (VCRTT/)
The second moment of total quanta was defined as the variance of the chromatographic retention time (VCRTT), a degree of residence as Equation (3).
Total quantum statistical moment similarity
The two TQSM parameters, first moment (MCRTT), i.e., mean, and second moment (VCRTT), i.e., variance, illustrate the mean chromatographic retention time center and discrete degree for the chromatographic fingerprints [Figure 1]. The two parameters can be converted into a normal distribution probability density function,, as shown in the following Equation (4).
|Figure 1: Total quantum statistical moment similarity modeling process for chromatographic fingerprint. (a and b) Chromatographic fingerprints of two samples. (c-e) Three scenarios of total quantum statistical moment similarity model: (c) two chromatographic fingerprints were overlapped completely and the similarity was 1; (d) only one cross-point exists; (e) two cross-point exists|
Click here to view
Where t means the chromatographic retention time,tT is the MCRTT, and is the VCRTT. Assume that the first moments for two chromatographic fingerprints are ta and tb and second moments are and . The intersection point of two normal curves of distribution represents t1and t2, yielding Equation (5).
Then, Equation (5) was regulated, and Equation (6) was obtained.
Finally, two solutions following Equation (7) are given.,
The TQSMS for two chromatographic fingerprints can then be defined as the overlapping area for two probability density functions depicted in [Figure 1], as in Equation (8).
There were three scenarios as follows:
- When as shown in [Figure 1]b in [Figure 1]c, two chromatographic fingerprints were overlapped completely and the similarity ST was 1
- When as shown in [Figure 1]a in [Figure 1]d, only one cross-point t1 exists as in Equation (9)
Then, TQSMS can be calculated by Equation (10).
- When as shown in [Figure 1]c in [Figure 1]e, two cross-points both t1 and t2 exist in two probability density functions where similarity ST is calculated by Equation (8).
Standard total quantum statistical moment similarity, deviation, variable probability, positive or negative judgment, and critical values for standard normal distribution
Similar to our previous study on TQSMS of pharmacokinetics, the parameters of the statistical test for TQSMS of chromatographic fingerprints can be obtained. Under test size uα, the parameters, standard TQSMS (TQSMSu), deviation (D), variable probability (pV), and the confidence of probability (1 − β), can be calculated and shown in [Supplementary Table 1] for negative judgment and [Supplementary Table 2] for positive judgment, and finally, their critical values can also be ascertained. With this, it is convenient to make judgment about similarity or differences in the chromatographic fingerprints.
When D = 1.96, TQSMSu would be 0.05, and the pV is 95%, whereas the confidence of probability (1 − β) varied with the confidence coefficient α. If a negative judgment was made that the two chromatographic fingerprints were different, then [Supplementary Table 1] should be adopted. When α value is 0.05, 1 − β is 0.5, i.e., there is 50% confidence of probability to make a judgment on 95% of chromatographic fingerprints being different. When 1 - β value is more than 0.75 and TQSMSu is less than 0.008 (α = 0.05, significance level) or 0.001 (α = 0.01), it can be considered as the critical value of negative conclusion that two chromatographic fingerprints are from different populations. If a positive judgment was made that the two chromatographic fingerprints were similar, then [Supplementary Table 2] should be adopted. When 1 − α is valued 0.95, β = 0.95, TQSMSu = 0.803, i.e., there is 95% confidence of probability to make a judgment on 95% of two fingerprints are similar, and then it can be considered as a critical value to a positive conclusion that the two samples were from a same population; as the routine requirement that 1 − α value was less than 0.95, β value is more than 0.900, and TQSMSu is more than 0.803, then there is 90% confidence of probability to make a judgment on 95% of two samples are similar. In other words, the critical value of TQSMSu to determine whether the two fingerprints are similar is 0.8030, which is an important and objective statistical parameter to judge the similarity between two samples. Within 25%–75% confidence of probability, a precise judgment would be made that the two chromatographic fingerprints are from the same population, while the risk of these conclusions can also be estimated. The availability of statistical test parameters is an important feature of the TQSMS model superior to other similarity methods.
Chemicals and materials
Standards 2-undecanone (Lot. 110834-200502, purity ≥ 99.8%), α-pinene (Lot. 897-2000001, purity ≥ 98.0%), and internal standard n-pentadecane (Lot. 11677-200401, purity ≥ 100.0%) were obtained from the China National Institute for the Control of Pharmaceutical and Biological Products (Beijing, China). Ethyl acetate (Lot. 013092701, purity ≥ 99.5%) was purchased from Chengdu Chron Chemicals (Chengdu, China), n-hexane (Lot. K46764991 524, purity ≥ 99.9%) from Merck Ltd. (Darmstadt, German), and ethanol (Lot. 0907360, purity ≥ 99.7%) from Anhui Ante Biochemistry Co. Ltd. (Suzhou, China). All samples of fresh HC were collected from Changsha Hunan from April 2017 to August 2017, and Prof. Ji-Lian Shi (Professor from Hunan University of Chinese Medicine, Changsha, China) authenticated the plants. n-Pentadecane (0.1358 g) was dissolved in 10 mL of n-hexane to and diluted to obtain a final concentration of 0.679 mg/mL. Briefly, 2-undecanone (0.0726 g) and α-pinene (0.1393 g) were dissolved in n-pentadecane internal standard solution (0.679 mg/mL) to a final volume of 10 mL and final concentrations of 7.29 and 13.93 mg/mL. All these solutions were kept at 4°C until for gas chromatography-mass spectrometry (GC-MS) analysis.
Apparatus and chromatographic conditions
The procedure for the extraction of HC volatile oil and GC-MS analysis is presented in our previous study. Volatile oil of HC was extracted by an extractor apparatus purchased from Sichuan Shubo (Group) Co., Ltd. (Chongzhou, China). GC-MS-QP2010 (SHIMADZU, Japan) was used for GC-MS analysis, and a quartz capillary column SE-30 (the stationary phase: AT SE-30, 0.25 mm × 30 m × 0.25 μm, Dalian Physiochemical Institute, China) was used for chromatographic separation. The injection volume of 2.0 μL was used for analysis. The injection port temperature was held at 250°C. The oven temperature was kept at 60°C for 3 min and then increased up to 140°C at 2°C/min, held at 140°C for 5 min and then increased up to 200°C at 10°C/min, and finally held at 200°C for 5 min. The total flow rate was kept at 37.1 mL/min, the column pressure was maintained at 65.2 kPa, and the temperature of transfer line and source was maintained at 230°C, with split mode (ratio 30:1). The temperature of ESI ion source was held at 230°C, the interface at 280°C, and the quadruple temperature at 150°C; the electron energy was maintained at 70 eV. SCAN mode was used for detection. The solvent cutoff time was 3.5 min. Constituents were identified by contrasting their mass spectra of chromatographic peak from NIST08.
A total of 15 batches of fresh HC were weighed and chopped into 2–3 cm pieces. The HC volatile oil was extracted five times by steam distillation with water. The distillate was collected (100%) (2.0 mL) and sequentially labeled as HC-01–HC-15 and stored at −20°C until GC-MS analysis. To perform GC-MS analysis, 0.5 μL of the sample was diluted to 1.0 mL volume with ethyl acetate.
From previously reported systems, integral conditions for the chromatographic fingerprint of HC volatile oil were obtained as follows: peak height of 200 μv, peak area of 4000 μv × s, and drift value of 15 μv.
Standard solution was determined for precision by five continuous times. The RSD of the peak retention time was <2.56% and the RSD of the peak area was <2.83%, indicating that the precision of the method meets standard requirements.
Stability of the sample solutions was tested by comparing sample solutions that were kept at room temperature with standard solutions in 1, 2, 4, 6, 12, and 24 h. We found that the sample solutions were stable within 24 h (the RSD of peak retention time <2.45%, the RSD of peak area <4.62%).
Five independent samples were prepared and analyzed to determine the repeatability. The RSD of peak retention time was <0.44% and the RSD of peak area was <4.10%, indicating that the repeatability of the method met standard requirements.
The chromatographic fingerprints obtained were analyzed with TQSMS method established in this study according to Equations (1)–(10). Supplementary [Table 1] and [Table 2] are applied for the statistical tests. Furthermore, the included angle cosine similarities and correlation coefficient were also used to compare with the TQSMS method.
|Table 1: Retention time, peak area, and relative content of common chemical compounds of volatile oils from 15 batches Houttuynia cordata|
Click here to view
|Table 2: Relative content (%) of common chemical compounds of volatile oils from Houttuynia cordata of different samples|
Click here to view
| Results and Discussion|| |
Identification of volatile compounds in Houttuynia cordata
The volatile oils from the 15 batches of HC samples were analyzed by GC-MS. [Figure 2] shows the representative chromatograms. The peak area of each component was chosen as the analytical signal for the relative content, and these identified components from the 15 batches of HC are listed in Supplementary [Table 3], [Table 4], [Table 5], [Table 6], [Table 7], [Table 8], [Table 9], [Table 10], [Table 11], [Table 12], [Table 13], [Table 14], [Table 15], [Table 16], [Table 17]. The composition and relative content of volatile oil obtained by steam distillation of 15 batches of HC were all different. There were nine peaks in common, and a total of 733 chemical constituents were observed among 15 batches of samples. [Figure 3] shows the mass scan spectra and chemical structures of the nine common compounds [Figure 2]. [Table 1] shows the retention time, peak area, and relative content of common chemical compounds of essential oils from 15 batches of HC. The RSD of retention time of nine common chemical compounds ranged from 0.47% to 1.90%, indicating that the precision of the method met the standard requirements. The RSD of peak area and relative content of nine common chemical compounds ranged from 33.79% to 182.04% and from 37.45% to 101.04%, respectively. As outlined in [Table 2] and [Figure 4], the common components are represented from 20.74% to 74.13% of the extracted volatile oils. The identified volatile constituents consisted of aromatic, aliphatic, and terpenoid compounds. The most abundant component of the HC volatile oils was 2-undecanone (4.24–29.88%). Furthermore, the volatile oils also contained β-myrcene (5.41%–27.00%), β-pinene (0.83–17.90%), α-pinene (0.37%–10.46%); 2-tridecanone (0.76%–8.80%); bicyclo[3.1.0]hexane, 4-methylene-1-(1-methylethyl)- (0.02%–8.09%); n-decanoic acid (0.53%–6.62%), camphene (0.13%–1.86%), and n-hexadecanoic acid (0.26%–1.65%). However, the peak area percentages of other volatile constituents ranged from 79.26% to 25.87%. These results show that the composition of volatile oil and content of HC have a strong fluctuation, hinting that its quality control cannot be considered solely from the variety and content of the components. In situations like these, TQSMS method was a better tool to assess the similarities of these chromatographic fingerprints.
|Figure 2: Gas chromatography fingerprint of Houttuynia cordata volatile oil samples. S1–S15, respectively, represent the gas chromatography fingerprint of batch 1 to batch 15 Houttuynia cordata volatile oil samples. (a) .alpha.-pinene; (b) camphene; (c) bicyclo[3.1.0]hexane, 4-methylene-1-(1-methylethyl)-; (d) .beta.-pinene; (e) .beta.-myrcene; (f) 2-undecanone; (g) n-decanoic acid; (h) 2-tridecanone; (i) n-hexadecanoic acid|
Click here to view
|Table 3: Total quantum statistical moment parameters of gas chromatography fingerprints of 15 batches of Houttuynia cordata volatile oil|
Click here to view
|Table 4: Total quantum statistical moment similarity of gas chromatography fingerprints of 15 batches of Houttuynia cordata volatile oil|
Click here to view
|Table 5: Correlation coefficient of gas chromatography fingerprints of 15 batches of Houttuynia cordata volatile oil|
Click here to view
|Table 6: Angle cosine of gas chromatography fingerprints of 15 batches of Houttuynia cordata volatile oil|
Click here to view
|Figure 3: Ion fragmentation patterns and chemical structures of (a) .alpha.-pinene; (b) camphene; (c) bicyclo[3.1.0]hexane, 4-methylene-1-(1-methylethyl)-; (d) .beta.-pinene; (e) .beta.-myrcene; (f) 2-undecanone; (g) n-decanoic acid; (h) 2-tridecanone; (i) n-hexadecanoic acid|
Click here to view
|Figure 4: Score plot of relative content of common chemical compositions of volatile oils from 15 batches of Houttuynia cordata|
Click here to view
Total quantum statistical moment similarity analysis of gas chromatography-mass spectrometry fingerprint of volatile oil from Houttuynia cordata
Chromatographic fingerprints of volatile oil from HC were evaluated by the TQSMS method established in this study. The TQSMS parameters of the chromatographic fingerprints of the 15 batches of HC volatile oil are presented in [Table 3] and [Figure 5]. The peak numbers of chromatographic fingerprints of 15 batches of HC volatile oil were 49–137, with an RSD value of 30.13%; the sums of peak area were 1.159 × 107–3.437 × 108 μv × s, with an RSD value of 174.56%; MCRTT was 9.410–18.602 min, with an RSD value of 20.79%; VCRTT was 37.549–81.504, with an RSD value of 23.27%. The TQSMS of the chromatographic fingerprints of the 15 batches of HC volatile oil was also obtained [Table 4]. TQSMS close to 1 suggests a high similarity between the two chromatographic fingerprints. As outlined in [Table 4], the TQSMS was ranging from 0.4973 (HC-7 and HC-12) to 0.9905 (HC-4 and HC-8). According to the statistical test for TQSMS and our previous work, α error (Type-I error) = 0.05, β error (Type-II error) = 0.95, and ST = 0.8030, TQSMS of chromatographic fingerprints for HC volatile oil was shown as significant deviation among the batches of HC volatile oil. The reason for the occurrence of these differences is that the biosynthesis of secondary metabolites during the growth of the plant is closely related to soil, temperature, water quality, ecological environment, and other factors. Based on the values in [Table 4], samples from HC-01 to HC-09 were to have similarity (range from 0.8752 to 0.9905), whereas samples from HC-10 is similar to HC-15 (TQSMS = 0.9893), HC-11 is similar to HC-02 (TQSMS = 0.9705), HC-12 is similar to HC-14 (TQSMS = 0.8282), HC-13 is similar to HC-14 (TQSMS = 0.9264), HC-14 is similar to HC-13 (TQSMS = 0.9264), and HC-15 is similar to HC-10 (TQSMS = 0.9893). Otherwise, based on the mean TQSMS of each HC sample to other 14 samples, the order from largest to smallest is HC-08 (TQSMS = 0.9036), HC-11 (TQSMS = 0.9025), HC-04 (TQSMS = 0.9020), HC-03 (TQSMS = 0.8982), HC-02 (TQSMS = 0.8972), HC-05 (TQSMS = 0.8913), HC-06 (TQSMS = 0.8887), HC-15 (TQSMS = 0.8878), HC-01 (TQSMS = 0.8832), HC-10 (TQSMS = 0.8827), HC-09 (TQSMS = 0.8591), HC-07 (TQSMS = 0.8383), HC-13 (TQSMS = 0.8249), HC-14 (TQSMS = 0.7703), and HC-12 (TQSMS = 0.6510). In other words, samples HC-08, HC-11, and HC-04 have high similarities to other 14 samples. However, samples HC-12 and HC-14 showed significant differences from that of other 14 samples. These results indicate that the volatile oils in HC varied slightly. The TQSMS of samples HC-12 and 14 is different from other samples, probably due to some environmental factors., Therefore, the specific reasons for these differences and whether these differences are related to their efficacy need further research.
|Figure 5: Heatmap of total quantum statistical moment similarity of gas chromatography fingerprints of 15 batches of Houttuynia cordata volatile oil|
Click here to view
The included angle cosine and correlation coefficient of these fingerprints are, respectively, listed in [Table 5] and [Table 6], and the corresponding heatmaps are shown in [Figure 6] and [Figure 7]. As outlined in [Table 5], the correlation coefficient ranged between 0.3172 (HC-05 to HC-13) and 0.9987 (HC-02 to HC-05). Meanwhile, as outlined in [Table 6], the correlation coefficient ranged between 0.1183 (HC-11 to HC-12) and 0.9988 (HC-10 to HC-15). According to the principle of each similarity calculation method, the included angle cosine method and correlation coefficient method cannot accurately reflect the reality of multidimensional vector deviation degree. The two methods are more sensitive to the change of the peaks with a larger response and less sensitive to small peak. Therefore, the aforementioned two methods for fingerprints of herbal medicine with different batches, less characteristic peak, and more fingerprint peak are less effective. The composition and content of herbal medicine are complicated, which will be affected by a series of factors such as variety, origin, and processing conditions. Therefore, the similarity analysis of fingerprint is more important to the analysis of the whole spectrum for quality control of herbal medicine. The TQSMS method significantly reduces the requirements on the test method so that its similarity mainly reflects the similarity degree of the components, which can be used to analyze the fingerprint with no obvious characteristic peak. To sum up, TQSMS is an effective method to analyze the chromatographic fingerprints with some major outstanding characteristics, such as reducing the requirement of test method and having a coupling ability to couple with multiple varies to form a multidimensional functional curve and calculate the multi-dimensional TQSMS parameters. TQSMS model has been elucidated and established according to the TQSM parameters and normal distribution probability density function properties. Thus, we can quantitative analysis the correct probability to make positive and negative conclusions regardless of the origin of the samples with any confident coefficient α.
|Figure 6: Heatmap of correlation coefficient of gas chromatography fingerprints of 15 batches of Houttuynia cordata volatile oil|
Click here to view
|Figure 7: Heatmap of included angle cosine of gas chromatography fingerprints of 15 batches of Houttuynia cordata volatile oil|
Click here to view
| Conclusion|| |
Using statistical moment theories, the TQSMS was applied to analyze chromatographic fingerprints in this study. The volatile oil composition and content of HC have strong fluctuations. As a result, its quality control cannot be considered solely from the variety and content of its components. The TQSMS of chromatographic fingerprints of the 15 batches of HC volatile oil ranged from 0.4973 to 0.9905. Except for samples HC-12 to HC-14, the 13 other samples were found to be highly similar to each other, whereas samples HC-12 and 14 were found to be significantly different from others samples. Due to the simplicity of the TQSMS similarity method, it can be easily adopted and applied. Thus, TQSMS can be an additional method applied in the assessment of the similarity of two chromatographic fingerprints of herbal medicine or other complex systems with multiple components.
The authors are thankful to the help from colleagues of Hunan Key Laboratory of Druggability and Preparation Modification for Traditional Chinese Medicine, Yunfeng Lu research group of University of California Los Angeles, and China Scholarship.
Financial support and sponsorship
National Natural Science Foundation of China (Grant No. 81903759, 81874507, 81703824, 81573691, 81803729), the Natural Science Foundation of Hunan Province (Grant No. 2017JJ3236), the Youth Foundation of Hunan Province Department of Education (Grant No. 17B200), The First-Class Discipline of Pharmaceutical Science of Hunan (Grant No. 2018XY09), Changsha Science and Technology Bureau project (Grant No. kq2004059), the Open Fund of Hunan Key Laboratory of Druggability and Preparation Modification for Traditional Chinese Medicine supported the study.
Conflicts of interest
There are no conflicts of interest.
| References|| |
Tistaert C, Dejaegher B, Vander Heyden Y. Chromatographic separation techniques and data handling methods for herbal fingerprints: A review. Anal Chim Acta 2011;690:148-61.
Fan XH, Cheng YY, Ye ZL, Lin RC, Qian ZZ. Multiple chromatographic fingerprinting and its application to the quality control of herbal medicines. Anal Chim Acta 2006;555:217-24.
Xu CJ, Liang YZ, Chau FT, Heyden YV. Pretreatments of chromatographic fingerprints for quality control of herbal medicines. J Chromatogr A 2006;1134:253-9.
Gan F, Ye R. New approach on similarity analysis of chromatographic fingerprint of herbal medicine. J Chromatogr A 2006;1104:100-5.
Goodarzi M, Russell PJ, Vander Heyden Y. Similarity analyses of chromatographic herbal fingerprints: A review. Anal Chim Acta 2013;804:16-28.
Alaerts G, Van Erps J, Pieters S, Dumarey M, van Nederkassel AM, Goodarzi M, et al
. Similarity analyses of chromatographic fingerprints as tools for identification and quality control of green tea. J Chromatogr B Analyt Technol Biomed Life Sci 2012;910:61-70.
Yang B, Wang Y, Shan L, Zou J, Wu Y, Yang F, et al
. A novel and practical chromatographic “Fingerprint-ROC-SVM” strategy applied to quality analysis of traditional Chinese medicine injections: Using KuDieZi injection as a case study. Molecules 2017;22:1237.
Xie P, Chen S, Liang YZ, Wang X, Tian R, Upton R. Chromatographic fingerprint analysis—A rational approach for quality assessment of traditional Chinese herbal medicine. J Chromatogr A 2006;1112:171-80.
Xu C, Yang B, Zhu W, Li X, Tian J, Zhang L. Characterisation of polyphenol constituents of Linderae aggregate leaves using HPLC fingerprint analysis and their antioxidant activities. Food Chem 2015;186:83-9.
Xue S, Zhang R, Su J, Li X, Qiu M, Wang Y. 27 Fingerprint analysis of Clerodendranthus Spicatus
from different regions by high-performance liquid chromatography. Asian J Tradit Med 2016;11:1-10.
Zhong RX, Feng GF, Liu YJ, Zhong L, Wu CJ, Wan J. Origin discrimination and quality evaluation of Gastrodiae rhizoma
) by high-performance liquid chromatographic fingerprint. Trop J Pharm Res 2016;15:1535-41.
Zhang N, Song Y, Song Q, Shi S, Zhang Q, Zhao Y, et al
. Qualitative and quantitative assessments of Aconiti Lateralis
Radix Praeparata using high-performance liquid chromatography coupled with diode array detection and hybrid ion trap-time-of-flight mass spectrometry. J Chromatogr Sci 2016;54:888-901.
Committee SP. Pharmacopoeia of the People's Republic of China. Beijing: People's Medical Publishing House; 2015.
Chen YY, Liu JF, Chen CM, Chao PY, Chang TJ. A study of the antioxidative and antimutagenic effects of Houttuynia cordata
Thunb. using an oxidized frying oil-fed model. J Nutr Sci Vitaminol (Tokyo) 2003;49:327-33.
Park E, Kum S, Wang C, Park SY, Kim BS, Schuller-Levis G. Anti-inflammatory activity of herbal medicines: Inhibition of nitric oxide production and tumor necrosis factor-α secretion in an activated macrophage-like cell line. Am J Chinese Med 2005;33:415-24.
Chiang LC, Chang JS, Chen CC, Ng LT, Lin CC. Anti-herpes simplex virus activity of Bidens pilosa
and Houttuynia cordata.
Am J Chin Med 2003;31:355-62.
Kim GS, Kim DH, Lim JJ, Lee JJ, Han DY, Lee WM, et al
. Biological and antibacterial activities of the natural herb Houttuynia cordata
water extract against the intracellular bacterial pathogen salmonella within the RAW 264.7 macrophage. Biol Pharm Bull 2008;31:2012-7.
Lai KC, Chiu YJ, Tang YJ, Lin KL, Chiang JH, Jiang YL, et al. Houttuynia cordata
Thunb extract inhibits cell growth and induces apoptosis in human primary colorectal cancer cells. Anticancer Res 2010;30:3549-56.
Lau KM, Lee KM, Koon CM, Cheung CS, Lau CP, Ho HM, et al
. Immunomodulatory and anti-SARS activities of Houttuynia cordata
. J Ethnopharmacol 2008;118:79-85.
Kumar M, Prasad SK, Krishnamurthy S, Hemalatha S. Antihyperglycemic activity of Houttuynia cordata
Thunb. in streptozotocin-induced diabetic rats. Adv Pharmacol Sci 2014;2014:809438.
Miyata M, Koyama T, Yazawa K. Water extract of Houttuynia cordata
Thunb. leaves exerts anti-obesity effects by inhibiting fatty acid and glycerol absorption. J Nutr Sci Vitaminol (Tokyo) 2010;56:150-6.
Fu JG, Dai L, Lin Z, Lu HM. Houttuynia cordata
Thunb: A review of phytochemistry and pharmacology and quality control. Chinese Medicine 2013;4:101-23.
Xiao MF, Duan XP, Deng KW, Yang YT, Liu WL, He FY. Dynamic chromatopharmacokinetics model for components in Chinese medicine and validation in Buyang Huanwu Decoction. Zhongguo Zhong Yao Za Zhi 2019;44:574-81.
Cai CP, Wu NS. Statistical moment analysis and deconvolution of overlapping chromatographic peaks. Chromatographia 1991;31:595-9.
Martin B. An Introduction to Medical Statistics. 3rd
ed. Oxford: Oxford University Press; 2000. p. 45.
Aragao GM, Corradini MG, Normand MD, Peleg M. Evaluation of the Weibull and log normal distribution functions as survival models of Escherichia coli
under isothermal and non isothermal conditions. Int J Food Microbiol 2007;119:243-57.
Fan Q, Liu W, Yang Y, Zhou J, Tang Y, Xiao M, et al
. A new similarity method for assessment of pharmacokinetic interaction between flucloxacillin and midazolam. Pharmazie 2019;74:397-405.
Zhang YT, Xiao MF, Liao Q, Liu WL, Deng KW, Zhou YQ, et al
. Application of TQSM polypharmacokinetics and its similarity approach to ascertain Q-marker by analyses of transitivity in vivo
of five candidates in Buyanghuanwu injection. Phytomedicine 2018;45:18-25.
Liu WL, Zhang XL, Fan SQ, Zhu JP, Liang HH, Zhang YT, et al
. A novel concept of Q-markers: Molecular connectivity index. Phytomedicine 2018;45:36-40.
Zhou JF, Fang XQ, Zhang TM, Zhao Z, Zhu R, Xiang FQ, et al
. Quantitative similarity assessment of non-linear chemical fingerprint of traditional Chinese medicine by similarity system theory. J Cent South Univ T 2011;18:343-52.
Ashraf MA, Riaz M, Arif MS. Environmental Stress and Secondary Metabolites in Plants: An Overview. In: Plant Metabolites and Regulation under Environmental Stress. Elsevier: Amsterdam; 2018. p. 153-67.
Atyane LH, Lagram K, Ben EC, Lachheb M, Salaka L, Serghini MA, et al
. Study of the influence of geographical origin and environment conditions on the three secondary metabolites of Moroccan saffron by UV-visible spectrometry. Acta Horticulturae 2017;1184:267-72.
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7]
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5], [Table 6]