Spectrum-Effect Relationships Between Chemical Fingerprints and Antibacterial Effects of Lonicerae Japonicae Flos and Lonicerae Flos Base on UPLC and Microcalorimetry

The traditional Chinese medicines Lonicerae Japonicae Flos (LJF, Jinyinhua in Chinese) and Lonicerae Flos (LF, Shanyinhua in Chinese) refer to the flower buds of five plants belonging to the Caprifoliaceae family. Until 2000, all of these were officially listed as a single item, LJF (Jinyinhua), in the Chinese Pharmacopoeia. However, there have recently been many academic controversies concerning the separation and combination of LJF and LF in administrative regulation. Till now there has been little work completed evaluating the relationships between biological activity and chemical properties among these drugs. Microcalorimetry and UPLC were used along with principal component analysis (PCA), hierarchical cluster analysis (HCA), and canonical correlation analysis (CCA) to investigate the relationships between the chemical ingredients and the antibacterial effects of LJF and LF. Using multivariate statistical analysis, LJF and LF could be initially separated according to their chemical fingerprints, and the antibacterial effects of the two herbal drugs were divided into two clusters. This result supports the disaggregation of LJF and LF by the Pharmacopoeia Committee. However, the sample of Lonicera fulvotomentosa Hsu et S. C. Cheng turned out to be an intermediate species, with similar antibacterial efficacy as LJF. The results of CCA indicated that chlorogenic acid and 3,4-Dicaffeoylquinic acid were the major components generating antibacterial effects. Furthermore, 3,4-Dicaffeoylquinic acid could be used as a new marker ingredient for quality control of LJF and LF.

The traditional Chinese medicines Lonicerae Japonicae Flos (LJF, Jinyinhua in Chinese) and Lonicerae Flos (LF, Shanyinhua in Chinese) refer to the flower buds of five plants belonging to the Caprifoliaceae family. Until 2000, all of these were officially listed as a single item, LJF (Jinyinhua), in the Chinese Pharmacopoeia. However, there have recently been many academic controversies concerning the separation and combination of LJF and LF in administrative regulation. Till now there has been little work completed evaluating the relationships between biological activity and chemical properties among these drugs. Microcalorimetry and UPLC were used along with principal component analysis (PCA), hierarchical cluster analysis (HCA), and canonical correlation analysis (CCA) to investigate the relationships between the chemical ingredients and the antibacterial effects of LJF and LF. Using multivariate statistical analysis, LJF and LF could be initially separated according to their chemical fingerprints, and the antibacterial effects of the two herbal drugs were divided into two clusters. This result supports the disaggregation of LJF and LF by the Pharmacopoeia Committee. However, the sample of Lonicera fulvotomentosa Hsu et S. C. Cheng turned out to be an intermediate species, with similar antibacterial efficacy as LJF. The results of CCA indicated that chlorogenic acid and 3,4-Dicaffeoylquinic acid were the major components generating antibacterial effects. Furthermore, 3,4-Dicaffeoylquinic acid could be used as a new marker ingredient for quality control of LJF and LF.

INTRODUCTION
Lonicerae Japonicae Flos (LJF, Lonicera japonica Thunb., Jinyinhua in Chinese) is one of the most commonly used traditional Chinese medicines. Use of LJF was first recorded in Shen-Nong's Herbals, one of the world's earliest pharmacopeias, and it has been widely used throughout China for disease prevention and treatment. Specifically, LJF has a variety of bioactive effects, including antibacterial, anti-inflammatory (Tae et al., 2003), and antiviral properties as well as liver protection activity (Jiang et al., 2014). LJF is officially listed in the 2015 version of the Chinese pharmacopeia, where it is listed simultaneously with Lonicerae Flos (LF, Shanyinhua in Chinese). They have the same descriptions of flavor, meridian tropism, functions, dosage, and indications in the pharmacopeia (Commission, 2015). LF involves the flowers or buds of four plants: Lonicera macranthoides Hand.-Mazz., Lonicera hypoglauca Miq., Lonicera confusa DC., and Lonicera fulvotnetosa Hsu et S. C. Cheng. Few studies of LF exist, and most of them focus only on L. macranthoides Hand.-Mazz. (Li et al., 2015). Before the 2000 version of the Chinese pharmacopoeia, both LJF and LF were listed under the same category of LJF (Jinyinhua in Chinese). However, for medical safety, especially in regards to drug injection, the Pharmacopoeia Commission separated these into two herbal drugs. Based on the divergence of chemical ingredients, only Lonicera japonica Thunb. is listed in LJF (Jinyinhua in Chinese) in the 2015 version of the Chinese pharmacopoeia, while the other plant sources are listed in LF (Shanyinhua in Chinese). However, due to the similarity of these two drugs in clinical applications, there are many disputes about the administrative regulation of LJF and LF. It is therefore urgent to do further research to determine the pharmaceutical activity of the two herbal drugs.
Sensitive biological methods are needed to explore the diversity and bioactivity of these two medicines. Traditional bioassays are time-consuming and also consume large amounts of materials, and more importantly, they are not able to easily and clearly identify the different bioactivities of closely related species. Microcalorimetry is an ideal method to overcome these challenges because it is a real-time and quantitative means of measuring thermal dynamic changes across many fields, including studies of microorganisms (Kabanova et al., 2012) and pharmacological analysis. Essentially, useful bioactivity information can be obtained by measuring the thermokinetics of the drug-microorganism complex. P. aeruginosa was chosen as a model microorganism in this study because it is widespread in nature and a common bacteria found in wound infections (Wu et al., 2015). Furthermore, P. aeruginosa exhibited large thermo-kinetic changes, which improved testing and analysis.
Ultra Performance Liquid Chromatography (UPLC) is a highly efficient and accurate technique and is widely used to quantify the components of herbal drugs and generate chemical fingerprints. Chromatographic fingerprint methods can be used to characterize the holistic chemical profiles of herbal medicine. Chemical fingerprint is a useful tool for evaluation of herbal drugs' quality, differentiation of origin, identification of authenticity, and so on (Deng and Yang, 2013). Especially, it plays an important role in controlling the quality of herbal medicine without a reasonable strategy for ensuring the safety and efficacy of this herbal medicine .
In the present study, microcalorimetry, UPLC, PCA, hierarchical cluster analysis (HCA), and canonical correlation analysis (CCA) were used in combination in order to distinguish the differences in bioactivity and chemical properties between LJF and LF, with the ultimate goal of characterizing their spectrum-effect relationships.

Chemicals, Reagents, and Materials
Sixteen batches of Lonicera flowers buds were collected from various sources in China, as shown in Supplementary Table 1. Each species was authenticated by Professor Chunsheng Liu (School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing, China).

Sample Preparation
Fifty grams of fine powder from each batch of flower buds was weighed and then added to 500 ml of 70% ethanol solution. It was then extracted in a 100 Hz ultrasonic bath for 60 min, the solution was filtered, and the residue was extracted using the methods described above. After two rounds of extraction, the filtrate was combined and concentrated to dryness by evaporation. The residue was further dried in a vacuum freeze drier and then ground into powder.

Instruments
UPLC fingerprints were measured using a Waters Acquity UPLC ™ system (Waters, Milford, MA, USA) equipped with a binary solvent delivery pump, auto-sampler manager, column compartment, and a Photo-Diode Array (PDA) detector, connected to Waters Empower 2 software. The real-time metabolic HFP-time curves of P. aeruginosa were detected using a TAM 3114/3236 Bio-activity monitor (TA, Sweden) at 37 • C connected to Picolog software. The baseline stability of the microcalorimeter was 0.2 mW for 24 h.

Bacterial Strain and Recovery
P. aeruginosa (ATCC 27853) was supplied by the Clinical Examination Center of 302 Hospital of People's Liberation Army (PLA), Beijing, China and was stored in glycerinum solution at −80 • C. Prior to experiments the stock solution was inoculated into 25 mL Luria-Bertani (LB) culture medium. It was then placed in a 37 • C incubator shaker for 4 h with a rotational speed of 110 rpm and finally stored in a refrigerator at 4 • C.

Microcalorimetric Measurements
All experimental manipulations were performed under sterile conditions. An ampoule method was used for microcalorimetry experiments and the ampoules were cleaned and autoclaved before use. The suspension of P. aeruginosa was added into equal volumes of culture medium containing different sample solutions of extract powder dissolved in medium. Each 20 mL ampoule contained 10 mL LB culture medium that included sample and bacteria. These were hermetically incubated at 37 • C. The concentration of Lonicera flowers bud extract was 10 mg/mL. When the temperature of the sealed ampoules reached 37 • C (this took about 30 min) a stable baseline was recorded, and the bacterial thermogenic curves were recorded until the output returned to the initial baseline.

UPLC Conditions
Samples were separated on an ACQUITY UPLC HSS T3 C18 column (2.1 × 100 mm, 1.8 µm) at 28 • C using acetonitrile (A) and a 0.3% solution of formic acid in water (B). The gradient elution program utilized began with 10% A at 0-8.5 min, then 10-25% B at 8.5-9.5 min, 25% B at 9.5-16 min, 25-100% B at 16-25 min, and then maintained for 5 min. The flow-rate was kept constant at 0.2 mL/min. The effluent was monitored at 330 nm and the sample injection volume was 1.0 µL.

Hierarchical Clustering Analysis
HCA is a statistic method used to determine the natural structure of objects in multidimensional profiles. This analysis can divide individuals or objects into specific clusters in a way that maximizes homogeneity within each cluster while also maximizing heterogeneity between clusters. In this study HCA was used to assess correlations in the UPLC fingerprints.

Principal Components Analysis
Principal components analysis (PCA) is a statistical technique used to simplify a complex issue by transforming multiple variables into a small number of integrated variables while losing very little information. These integrated variables are known as principal components, and are uncorrelated linear combinations of the original variables, making these principal components superior compared to the original variables for many applications. In our study, PCA was used to process the thermo-kinetics parameters obtained from the power-time curve of the P. aeruginosa growth with different solutions of sample extracts.

Canonical Correlation Analysis
CCA is a multivariate analysis used to study the correlation between two sets of variables. It is used for the dimension reduction of PCA and to extract the main principal components, and then describes the whole linear relationship of two sets of variables by the relevance of two principal components. In our study, CCA was used to analyze the relevance between the peak area values from the UPLC fingerprints and thermo-kinetic data from microcalorimetry.

RESULT UPLC Fingerprints
Under optimized conditions chromatograms were generated for all batches of bud extracts (Figure 1) as well as for a mixture of reference substances (Figure 2). There were similar chemical profiles across different bud sources, and seven common peaks were found by comparing ultraviolet spectra and UPLC retention times for the 16 chromatograms. The methods used to identify common peaks have been referenced in similar reports (Zheng et al., 2014).
Seven common characteristic peaks with large areas and good resolution from consecutive peaks were collected systematically from the 16 chromatograms. Because the S15 chromatogram lacked peak X7, this peak area was defined as 0. By comparing the chromatograms of the samples to that of the mixture of reference substances, peaks X1, X2, X4, X5, and X6 were identified as 5-caffeoylquinic acid, Chlorogenic acid, Luteoloside, 3,4-Dicaffeoylquinic acid, and 3,5-Dicaffeoylquinic acid, respectively. Peaks X3 and X7 remain unknown compounds. Some differences were noted among data collected for the common peaks (Table 1), such that the coefficients of variance (C.V.%) for all common peaks were larger than 31.71%. This is due to the diversity in the content of chemical compounds contained in C.V.% = σ /µ*100; σ was the standard deviation, µ was the average value of peak area.
different batches. As a next step, HCA was used to further elucidate the relative contents of the samples.

HCA
Results of HCA showed that the fingerprints separated into three clusters, as shown in Figure 3: cluster one included S1-S8, cluster 2 contained only S9, and cluster three contained S10-S16. The samples separated into cluster one were batches of LJF from different places of origin, samples in cluster three were LF, and S9 in cluster two was the LF from Guizhou province. S9 is the species of LF closest to LJF in chemical profile, and is from the plant Lonicera fulvotomentosa Hsu et S. C. Cheng, one of the various species in the LF samples. The results of HCA showed that the extracts of LJF and LF from different origins had similar chemical fingerprints, yet HCA could initially separate the LJF and LF samples at the chemical level. Based on these results, a combination of HCA and UPLC with a PDA detector could, at least roughly, separate the bud extracts from two herbal drugs.

HFP-Time Curves of P. aeruginosa
An HFP-time curve was used to record the power change of the ampoules containing medium inoculated with P. aeruginosa at 37 • C. Because the ampoules were airtight, two peaks were recorded in the HFP-time curve, one showing the aerobic metabolism of the bacteria and the other due to the anaerobic metabolism. The InP-t curve reflected the metabolic character of P. aeruginosa. As is well known, bacterial growth can be divided into four phases, i.e., the initial exponential phase, the lag phase, the second exponential phase, and the decline phase. These phases are displayed in the HFP-time plot shown in the  By analyzing the HFP-time curve of the growth of P. aeruginosa, eight parameters (k 1 , k 2 , P m1 , P m2 , T m1 , T m2 , Q 1 , and Q 2 ) were extracted from the four phases. These could be used to quantify the effects of various bud extracts on P. aeruginosa, as shown in Supplementary Table 2. The values of k 1 and k 2 are the maximum growth rate constants of the two peaks, P m1 , P m2 , T m1 , and T m2 are the maximum heat output and the appearance time, respectively, of the two peaks, and Q 1 and Q 2 are the FIGURE 4 | HFP-time curves of P. aeruginosa growth without any extracts (blue), and the corresponding InP-time curve (red). The curves were divided into four phases (I-IV).The two curves reflect characteristic growth of P. aeruginosa at 37 • C. total quantity of heat in the aerobic and anaerobic stages of P. aeruginosa metabolism. Among the eight parameters, P m1 , P m2 , T m1 , T m2 , Q 1 , and Q 2 could be obtained directly from the HFP-time curve. However, k 1 and k 2 were calculated using the following equations (Xie et al., 1988): where P 0 is the Heatflow-power at time t = 0 and P t is the power at time t.
The HFP-time curves of P. aeruginosa after the addition of different bud extracts (Figure 5) were different from the control curve. This indicates that the bud extracts were able to influence bacterial growth to different extents. For example, the point of maximum heat of the P. aeruginosa growth was postponed for all samples, indicating that the medium that included extracts might inhibit the bacterial growth. This means that the bacterial culture would need more time to produce the defined number of bacteria needed to generate enough heat to achieve the intense signal. Furthermore, some curves measured using medium with extracts were severely deformed due to strong bioactivity on P. PCA was introduced as a further tool to analyze the thermokinetic parameters and further elucidate the extracts influence on bacteria.
PCA Through use of a PCA procedure it was possible to determine the most important variables, in other words, those variables that represent a majority of the total information. In this process, the values of thermo-kinetic parameters (k 1 , k 2 , T m1 , T m2 , P m1 , P m2 , Q 1 , and Q 2 ) were analyzed using SPSS statistical analysis software. The results showed that the first three principal components accounted for about 80% of the information contained in the original eight parameters.
The relationships of the three principal components from eight indexes referencing the eigenvalue of the correlation matrix were expressed as the following equations: These result demonstrate that parameters k 2 , P m2 , and T m2 might be the principle parameters, as they account for 79.37% of all components in the thermokinetic parameters (Supplementary Table 3).
Analysis of the thermokinetic parameters in the PCA score plots generated by the Simca-P software (Figure 6) showed that the 16 batches of bud extracts were clearly divided into two groups, LF and LJF. These results demonstrated that there was obviously different effects on P. aeruginosa due to LJF vs. LF. From Figure 6, different batches of LJF showed little diversity, illustrating that the bioactivity of LJF was relatively consistent. In contrast, there was a large divergence among different batches of LF. S2 and S13 had similar bioactivity effects on P. aeruginosa, and S9 appeared to have effects similar to those of LJF. Overall, these results were similar to those from the HCA of UPLC data. In order to elucidate the relationships between thermokinetic data and UPLC fingerprints, CCA was introduced as a next step.

Canonical Correlation Analysis
CCA was used to determine spectrum-effect relationships between area values of seven common peaks in the UPLC data and three principal components (k 2 , P m2 , and T m2 ) extracted from the previously calculated thermodynamic parameters. The outcome of this analysis is shown in Table 2. Five peaks: X1, X2, X3, X5, and X6 were highly correlated (|R|> 0.6) with the major thermodynamic parameters. The correlation coefficients also showed that P m2 was negatively correlated with X1, X2, X3, X4, X6, and X7. This indicates that the thermokinetic parameter P m2 could be inhibited by these compounds. The X6 peak, previously identified as Isochlorogenic acid A was the most correlated (R = −0.7489) to biological effect. Other correlated peaks X1, X2, and X5 had been previously identified as neochlorogenic acid (5-caffeoylquinic acid), Chlorogenic acid, and 3,4-Dicaffeoylquinic acid. However, the X3 peak, which also showed good correlation, remains an unknown compound. Furthermore, parameter T m2 had a marked positive correlation (R > 0.6) with peaks X2 and X5, demonstrating that chlorogenic acid and 3,4-Dicaffeoylquinic acid, respectively, could delay the appearance of the second peak. Inhibition of P m2 and delaying T m2 are both antibacterial FIGURE 6 | PCA plot of microcalorimetry parameters generated using Simca-P software. Samples of LJF appear as black points, while samples of LF appear as red points.
effects . Therefore, the above results indicated that chlorogenic acid and 3,4-Dicaffeoylquinic acid are the major antibacterial compounds in LJF and LF. Next the content levels of these compounds in LJF and LF were analyzed, as shown in Figure 7. The content of 3,4-Dicaffeoylquinic acid was significantly different between the two herbal drug samples, and could therefore be a marker for quality control in LJF and LF.

DISCUSSION
Nowadays the strategy for quality control of herbal drugs depends mostly on measuring the content of constituents, which is correlated with the safety and efficacy of herbal medicine (Capasso et al., 2000). However, the relationships and efficacies of specific constituents in herbal drugs remain ambiguous. In this study, a bacterial model was used to evaluate the antibacterial effects of LJF and LF. The results indicated that the effects of LJF on P. aeruginosa were similar regardless of the geographical origins of LJF. This, suggested suggesting that the quality of LJF was more stable as compared to LF. Microcalorimetry could be used to provide further information regarding bioactivity . Furthermore, biological thermokinetic parameters can reflect antibacterial effects of whole compounds in the herbal extractions. Although chemical fingerprints could show the differences between most constituents in herbal medicine, it could not directly show the bioactivity of herbal drugs.  Based on these observations, it is a promising strategy to combine the use of chemical fingerprints and microcalorimetry in order to ensure the quality of herbal medicine. In our study which utilized statistical analysis, two major components were found to be primarily responsible for the antibacterial effects of herbal drugs. It was particularly interesting that the content of one of these components, 3,5-Dicaffeoylquinic acid, was significantly different between LJF and LF. This could be used in the future as a new marker for quality control of LJF and LF. Species and geographical origin influence the distinct content within plants (Canas et al., 2000). Only flower buds from Lonicera japonica Thunb. are defined as LJF. The LF samples were collected from three plant species for our study. The results of PCA showed that the samples of LJF were relatively concentrated although they were collected from different places. However, the LF samples were much more spread out. Most of the LF samples were from a single species, Lonicera macranthoides Hand.-Mazz., yet these samples were also quite discrete. This result demonstrates that when considering factors that can lead to different biological effects, geographical origin should obviously be included. After analyzing the geographical environment of the LJF and LF samples, it was determined that the sources of LJF were concentrated on the North China Plain, while most LF samples were widely distributed across mountainous Southwest China. Southwest China is quite different from the North China Plain in geographical characteristics, climate, and seasonality. As an adaptive consequence of the complicated ecology, the plants in the Southwest area conserve high levels of diversity, an effect which could contribute to the discrete bioactivities of the LF samples.
Ensuring the safety and efficacy of drugs involves multiple considerations and the quality of the drug must be fundamentally guaranteed. An important part of drug quality control is to ensure consistent medical and biological effects are delivered by the same drug dosage (Busse, 2000). In our research, there is a divergence in bioactivity between LJF and LF using the bacterial model. Separation of LF from LJF in the Chinese Pharmacopoeia could to a certain degree ensure the accuracy of clinical applications. However, microcalorimetry and chemical fingerprint results for L. fulvotomentosa Hsu et S. C. Cheng, a species of LF, closely approached those for LJF. This indicates that the buds of L. fulvotomentosa Hsu et S. C. Cheng may be used as LJF. Based on these complex results, there is an urgent need for further and more fundamental research on LJF and LF. There are many areas influenced by the modification of the Chinese Pharmacopoeia, such as history, economy, politics, civilization, and ethnic habits. Any change to the Chinese Pharmacopoeia could have far-reaching impacts. Therefore, the disaggregation of LJF and LF was not only a scientific question, but a much more comprehensive issue. The future will require that scientific researchers and administrators work together to study and develop relevant policies for promoting the development of the traditional Chinese medicine industry.

AUTHOR CONTRIBUTION
ZS did the writing of paper and UPLC experiments as well as microcalorimetric experiments; ZL did the UPLC experiments; CL and YZ did the plant identification; MW and HS did the microcalorimetric experiments; XM did the statistical analysis; JW, YZ, and XX supervised the project. All the authors read and approved the final manuscript.