Comprehensive chemical profiling of volatile constituents of Angong Niuhuang Pill in vitro and in vivo based on gas chromatography coupled with mass spectrometry

Background Angong Niuhuang Pill (ANP), a renowned precious traditional Chinese medicine prescription, is extensively utilized for the clinical treatment of stroke, meningitis and encephalorrhagia in China. As a classic resuscitation-inducing aromatic prescription, ANP has been investigated for its pharmacological effects in recent years, while the volatile composition in ANP still lacks comprehensive elucidation. Method To better explore the volatile constituents in ANP, a qualitative analysis method was developed based on gas chromatography coupled with mass spectrometry. Furthermore, a validated quantitative method was established to determine 21 main compounds in 8 batches of commercially available ANP samples by gas chromatography-tandem mass spectrometry. The quantitative data were successively subjected to Pearson correlation coefficient analysis. Additionally, the absorbed volatile constituents in rat plasma after single oral administration of ANP have also been characterized. Results A total of 93 volatile constituents including 29 sesquiterpenoids, 28 monoterpenoids, 13 fatty acids and their esters, 7 alkanes, 6 ketones, 3 phenols, 3 aldehydes, 2 benzoate esters, and 2 other types, were preliminarily characterized, which primarily originated from Borneolum, Moschus, Curcumae Radix, and Gardeniae Fructus. d-Borneol, isoborneol and muscone were the top three abundant ingredients (> 600 μg/g) in 8 batches of ANP samples. Subsequently, the average Pearson correlation coefficient of the contents of 21 analytes was 0.993, inferring the high batch-to-batch similarity among 8 batches. After oral administration of ANP, d-borneol, isoborneol, muscone and camphor were the main volatile constituents absorbed in the rat plasma. Conclusion This research may be helpful for the comprehensive quality control study of ANP, and provide for guarantee the clinical efficacy of ANP. Supplementary Information The online version contains supplementary material available at 10.1186/s13020-022-00659-8.


recorded in Treatise on Differentiation and Treatment of Epidemic Febrile Diseases in the Qing Dynasty and has a history of over 200 years of application in China.
According to the TCM theory and clinical practice, ANP has been proved to be beneficial to the treatment of various central nervous system diseases, such as stroke coma, meningitis and intracerebral hemorrhage [1]. A recent meta-analysis of 18 trials involving 1601 patients reported that adjuvant treatment with ANP could significantly improve the total response rate and reduce the neurologic deficit score in patients with acute cerebral infarction and acute intracerebral hemorrhage [2].
Generally, ANP contains volatile and non-volatile constituents, contributing to an integral part of the overall efficacy of ANP. The volatile substances of aromatic drugs are mainly deemed to process the important properties of resuscitation and awakening in traditional Chinese medicines. Since ANP is a typical resuscitation-inducing aromatic prescription, its therapeutic effect on cerebrovascular diseases may be closely related to the volatile components. For instance, muscone, d-borneol, and isoborneol, the representative components from Moschus, Borneolum, respectively, showed various biological activities in vitro and in vivo. Muscone might treat myocardial infarction and protect cardiovascular and cerebrovascular system [3,4], regulate neuroprotective system [5] and improve osteoarticular injuries [6]. Isoborneol and d-borneol had the characteristics of anti-inflammatory [7], antiatherosclerosis [8] and promoting penetration [9], and played significant roles in the treatment of cerebrovascular diseases [12]. However, there are few researches on systematical characterization the volatile chemical composition in ANP.
In this study, the volatile constituents in ANP were comprehensively characterized based on the sensitive gas chromatography coupled with mass spectrometry (GC-MS) in vitro and in vivo. The GC-MS method was complemented for the volatile compounds identification based on database matching, retention indices and standard references confirmation. Then, 21 compounds were quantified in 8 batches of ANP by gas chromatography-tandem mass spectrometry (GC-MS/ MS), and the batch-to-batch similarity was evaluated by Pearson correlation coefficient (PCC) analysis. Furthermore, the volatile constituent of ANP in vivo were characterized in the plasma of rats which was single oral administrated. To our knowledge, the determination of main volatile compounds in ANP was reported for the first time. It is hoped that the results can provide a valuable reference for quality control and clinical efficacy research of ANP.
Anhydrous ethanol and ethyl acetate (HPLC grade) were purchased from Yonghua Chemical Technology Co., Ltd (Suzhou, China). Distilled water was prepared using a Milli-Q Integral water purification system (Millipore, Bedford, MA, USA). Other reagents were analytical grade.

Preparation of standard and sample solutions
The 21 reference standards were dissolved with anhydrous ethanol respectively to prepare the corresponding stock solutions at the concentration of 1 mg/mL. An appropriate amount of individual standard stock solutions was mixed into a 10 mL volumetric flask to prepare a mixed standard stock solution, which was diluted with anhydrous ethanol to obtain a series of working solutions at proper concentrations for calibration curves. All working solutions were stored at − 20 °C until analysis.
After ground into fine powder, 1.5 g ANP was accurately weighed and extracted using steam distillation for 4 h with 100 mL distilled water. The volatile oils were collected and residual water was removed with anhydrous sodium sulfate, finally were dissolved in 2 mL ethyl acetate. The solution was directly analyzed after adding n-tridecane as an internal standard (IS). Since some compounds of high contents such as isoborneol, d-borneol, β-caryophyllene, ar-turmerone and muscone were overloaded, the solution was also analyzed after diluted by 40 times with ethyl acetate and added with IS. The n-tridecane in all injection samples was at the final concentration of 14.20 μg/mL.

Collection and precipitation procedure of plasma samples
Male Wistar rats (Animal certificate number: SCXK [HU]-2018-0016) were purchased from Shanghai Lab. Animal Research Center (Shanghai, China). All rats were housed at 24 ± 2 °C on a 12 h light/dark cycle, and fed a standard diet and water for 1 week before the experiment. Then all rats (200 ± 20 g) were randomly divided into two groups: ANP group and control group. The rats were fasted for 12 h prior to experiments but water was provided ad libitum. ANP suspension dissolved in physiological saline was orally administrated to nine rats (ANP group) at a dosage of 8.1 g/kg, and the rest two rats (control group) were orally administered with the same dose of saline respectively. All procedures were carried out in accordance with Guide for the Care and Use of Laboratory Animals (National Institutes of Health).
The blood samples were collected in heparinized 1.5-mL polythene tubes at 0, 0.25, 0.5, 0.75, 1, 1.5, 2 and 4 h after oral administration. At each time point, the blood samples of 3 rats (ANP group) or 2 rats (control group) were mixed into one sample, and were centrifuged immediately at 4500 rpm for 10 min at 4 °C to collect the plasma, which were stored at − 80 °C before analysis. For GC-MS analysis, a 100 μL aliquot of plasma sample was added with 100 μL of ethyl acetate, followed by vortex-mixing for 1 min and centrifugation for 10 min at 13,000 rpm. An aliquot of supernatant was transferred into an injection vial for analysis.

Chromatographic and mass spectrometric conditions
The qualitative analysis was performed on an Agilent 7890B gas chromatography system coupled to an Agilent 5977A quadrupole mass spectrometer (GC-MS, Agilent Technologies, USA). The GC separation was achieved on an Agilent DB-5 MS capillary column (60 m × 0.25 mm i.d.) coated with a 0.25 μm film of 5% phenyl polymethyl siloxane. High-purity helium gas was used as carrier gas, with a flow rate of 1.0 mL/min. The injection and interface temperatures were set to 280 °C and 250 °C, respectively. The column temperature was programmed as follows: the oven was initially maintained at 60 °C for 3 min, sequentially increased to 100 °C at 10 °C/min, increased to 135 °C at 2 °C/min, then increased to 165 °C at 5 °C/min, increased to 168 °C at 1 °C/min, increased to 200 °C at 5 °C/min, held for 2 min, and then increased to 295 °C at 10 °C/min, held constant for 5 min. The split ratio was set to 10:1. The injection volume was 1 μL. The electron energy was set to 70 eV. The source temperature was 230 °C and the quadrupole temperature was 150 °C. The mass data was acquired using full scan mode with a mass range of m/z 50-600 after a solvent delay of 9.1 min. Data acquisition was obtained by Agilent MassHunter GC-MS Acquisition Software Version B.07.03.2129.
For quantitative analysis, the Agilent 7890B GC system equipped with an Agilent 7000D triple quadrupole mass spectrometer (GC-MS/MS, Agilent Technologies, USA) was performed. The Agilent DB-5 MS UI capillary column (30 m × 0.25 mm i.d., 0.25 μm) were used for separation, the initial oven temperature was 60 °C, and raised to 76 °C with 8 °C/min, then increased to 82 °C at the rate of 2 °C/min, and ramped with 15 °C/min to 130 °C, further rose at 20 °C/min to 230 °C and held for 2 min. The injection and interface temperatures were both kept at 250 °C. The carrier gas (helium, > 99.999%) flow rate was set at 1.0 mL/min. The electron energy and source temperature were set to 70 eV and 230 °C, respectively. The flow rate of quenching gas (helium, > 99.999%) and collision gas (nitrogen, > 99.999%) was 2.25 mL/min and 1.5 mL/min, separately. Data acquisition was achieved on Agilent MassHunter Workstation GC/MS Data Acquisition Software Version 10.0.368.

Validation of quantitative method Calibration curves, LODs and LOQs
After added with an equal amount of internal standard (IS) n-tridecane (final concentration: 14.20 μg/mL), a series of working solutions at multiple concentrations were analyzed. The calibration curves were constructed by plotting the relationships between peak area and concentration of the analytes and IS.
The mixed standard stock solution was further diluted with anhydrous ethanol and analyzed by GC-MS/MS. The concentration of the analyte with signal-to-noise ratio (S/N) of 10 was defined as the limit of quantification (LOQ), and the concentration of the analyte with S/N of 3 was assigned as the limit of detection (LOD).

Precision, repeatability, stability and recovery
The precision was evaluated by the determination of intra-and inter-day variances. The mixed standard solutions at three different concentration levels (low, medium, high) were analyzed for consecutive 3 days and continuous six times per day, and the peak area ratio of each compound and IS was recorded to calculate the relative standard deviation (RSD) value, respectively. Six parallel ANP samples from the same batch were prepared and analyzed for the repeatability test. To confirm the stability, a single sample solution was stored in sample chamber and analyzed at 0, 2, 4, 8, 12 and 24 h, respectively. Recovery test was used to verify the accuracy of the established method. A known amount of mixed reference solution at middle concentration was added in the same sample for six parallel extraction and analysis.

Data analysis
The qualitative analysis was realized on Agilent Mass-Hunter Workstation Qualitative Analysis Software 10.0. The retention indices of all chromatographic peaks were calculated based on the data of C7-C40 n-alkanes acquired with the same GC-MS method. The acquired components were tentatively identified by comparison with mass spectra and retention indices in the National Institute of Standards and Technology (NIST) 2017 library, and some were unambiguously determined by direct comparison with reference standards. Agilent MassHunter Workstation Quantitative Analysis Software Version 10.0 was used for quantitative analysis. Related graphical analysis was conducted on Graphpad Prism 8.02 (San Diego, USA).
The pairwise PCC between batches was based on the vectors of concentrations from each batch as Formula (1). Here, B1 and B2 represented the vectors of ingredients' concentration of the first and second batches respectively. Finally, the PCCs were displayed through the correlation heat map, the batch-to-batch similarity was (1) evaluated according to the average value of PCCs among all pairs from batches.

Characterization of volatile constituents in ANP
For the qualitative analysis, the temperature-programmed conditions of GC-MS were optimized in terms of resolution, symmetricity and capacity of chromatographic peaks. Then, the mixed reference standard solution and ANP sample solution were analyzed by the optimized GC-MS method, and the total ion chromatogram of the essential oil from ANP is shown in Fig. 1. The volatile constituents of ANP were identified by comparing their mass spectra and retention indices with those of NIST library, and accurately identified by matching with authentic standards. As illustrated in Fig. 1 and Table 1, a total of 93 components, including 29 sesquiterpenoids, 28 monoterpeneoids, 13 fatty acids and their esters, 7 alkanes, 6 ketones, 3 phenols, 3 aldehydes, 2 benzoate esters and 2 other type components, were preliminarily characterized. Among them, 21 compounds were determined by comparison with reference standards. The essential oil was dominated by d-borneol, isoborneol and muscone, making up 90.69% of the identified oil composition.   Among the 93 preliminarily identified volatile compounds, 11 components might be derived from Moschus, including 1 ketone, 1 benzoate ester and 9 fatty acids and their esters. The most abundant composition in Moschus was muscone (Peak 75) [13], which was also one of the principal constituents in ANP. Furthermore, 19 compounds comprising 16 monoterpeneoids, 2 fatty acid esters and 1 alkane, were tentatively assigned to Borneolum, of which the highest were d-borneol (Peak 29) and isoborneol (Peak 28) [7]. Additionally, 46 components were originated from Curcumae Radix, consisting of 29 sesquiterpenoids, 13 monoterpeneoids, 1 phenol, 1 ketone, 1 fatty acid ester, and 1 alcohol. Chinese Pharmacopoeia recorded its four botanical origins, including Curcuma wenyujin Y. H. Chen et C. Ling, Curcuma longa L., Curcuma phaeocaulis Val. and Curcuma kwangsiensis S. G. Lee et C. F. Liang. Sesquiterpenoid was one of the most abundant types in Curcumae Radix. Baesd on the characteristic constituents including ar-Tumerone (Peak 62), β-caryophyllene (Peak 42) and α-curcumene (Peak 48), the Curcumae Radix in those ANP could be traced to Curcuma longa L. [14][15][16]. Furthermore, 32 volatile components in ANP could be sourced from Gardeniae Fructus, containing 10 fatty acids and their esters, 6 alkanes, 5 monoterpeneoids, 4 ketones, 3 aldehydes, 2 phenols, and 1 furan and 1 benzoate ester. The constituents containing "-cyclohexane-trimethyl-" group represented by 4-methyleneisophorone (Peak 34) and β-Isophorone (Peak 14) might be the characteristic aromatic compounds in Gardeniae Fructus [17], which were detected in ANP. The detailed sources of volatile constituents in ANP are shown in Table 1 and Additional file 2: Fig. S1.

Quantitative analysis of representative volatile constituents in ANP
For the quantitative analysis, a temperature-programmed process of GC-MS/MS for 21 volatile constituents was developed within 15 min. A dynamic multiple reaction monitoring (dMRM) pattern was selected for quantification. As an alternative for multiple reaction monitoring, dMRM can monitor analytes automatically around the expected retention time, without defining a specific time period for the selected transitions, thereby decreasing concurrent transitions and improving the sensitivity [18,19]. The segmented dMRM chromatogram of 21 target compounds is shown in Fig. 2. The dMRM parameters for all analytes and IS are presented in Additional file 1: Table S2.

Method validation
The results of methodological study are shown in Additional file 1: Tables S3, S4. For the calibration curve established, all analytes showed good linearity within their linear ranges and their correlation coefficient R 2 values varied from 0.9918 to 0.9997. The LODs and LOQs ranged from 0.60 to 150.40 ng/mL and from 6.00 to 451.20 ng/mL, respectively, indicating this method was sensitive enough for the quantitative analysis of the major volatile constituents of ANP. The RSD values of the intra-and inter-day were in the ranges of 0.04-5.84% and 0.48-6.83%, respectively. The RSD values for repeatability of all compounds ranged from 1.13 to 9.00%, which showed that the repeatability of tested samples was good.
The RSD values for stability of all analytes were with the range of 1.39-6.70%, indicating that the samples in sample chamber were relatively stable within 24 h. The overall recoveries of 21 compounds were in a range of 90.80-109.83% with RSD values less than 8.68%, which

Contents of representative volatile constituents in ANP
In total, eight batches of ANP were analyzed by the developed and verified GC-MS/MS method, and the quantitative results are presented in Fig. 3 and Additional file 1: Table S5. It was found that d-borneol (6727.63-7734.92 μg/g), isoborneol (4614.45-5620.61 μg/g) and muscone (637.90-694.74 μg/g) were the top three abundant ingredients in all batches of ANP samples, making up 95% of the total contents at least. Notably, the characteristic ar-tumerone (186.33-324.22 μg/g), β-caryophyllene (106.81-244.95 μg/g), and α-curcumene (20.16-80.29 μg/g) were of high content level. The content of each component varied greatly, while the content distribution of each component in the eight batches was relatively consistent.

Batch-to-batch similarity evaluation of ANP samples
In order to evaluate the batch-to-batch similarity of multiple ANP samples, the quantitative data of 8 batches were subjected to PCC analysis. For the sake of reducing the concentration distribution difference of different components, the concentrations were processed with logarithm first. Furthermore, the concentration heat map of eight batches was displayed in Additional file 2: Fig.  S2, indicating that the concentration distribution of these batches are generally similar. Subsequently, the batch-tobatch similarity was quantified by constructing the correlation coefficient matrix, and the heat map of it was shown in Fig. 4. It could be found that the average PCC was 0.993, and even the lowest PCC between S2 and S5 was also more than 0.986, further revealing the high similarity among the eight batches.

Characterization of volatile constituents of ANP in vivo
To clarify the absorbed volatile component in vivo, the rat plasma after oral administration of ANP were analyzed by the developed GC-MS method. Since ANP was a classic first-aid Chinese patent medicine, single administration was performed. As shown in Fig. 5A, B, four absorbed components from ANP were characterized in 0.25 h, including three prototypes (P1-P3) and one metabolite (M1). By comparison with the retention time and mass spectral data of the reference standards, M1, P1, P2, and P3 were assigned to camphor, isoborneol (28), d-borneol (29) and muscone (75), respectively (Table 2), and their relative plasma concentration profiling in 0.25 to 4 h were depicted in Fig. 5C. The relative plasma concentration profiling of muscone could not be displayed because the content of muscone was extremely low at 0.25-4 h. The volatile constituents of most two   high content in ANP were quickly absorbed after oral administration, and then decreased rapidly. While the content of camphor, the 2-position oxide of isoborneol and borneol [20,21], remained relatively stable within 4 h, suggesting that isoborneol and borneol may be gradually metabolized into camphor. Curiously, no metabolites of muscone were detected, which were consistent with previously reports [22,23]. Additionally, few metabolites were characterized partly because that single administration not accumulated the drug in vivo.
Most notably, isoborneol, d-borneol and muscone played significant roles in the treatment of cerebrovascular diseases, corresponding to the main efficacy of ANP. Dong et al. [12] found that isoborneol and d-borneol may improve the function of neurovascular units through anti-apoptotic and anti-inflammatory properties so as to exert their protective effects against cerebral ischemia injury. Moreover, borneol and muscone acted as absorption enhancers in the blood-brain barrier by increasing paracellular and transcellular transport [11,24,25]. Muscone have been reported that could mediate neuroprotective effects against cerebral ischemia by preventing oxidative stress and Ca 2+ influx [10]. When combined with other drugs, they can effectively promote the penetration of drugs through the blood-brain barrier, which is of great significance for the development of treating cerebrovascular diseases.

Conclusion
In this study, GC-MS and GC-MS/MS were used to conduct a comprehensive qualitative and quantitative analysis of the representative volatile constituents in commercial ANP. Combined with retention indices and reference standards, the chemical characteristics of volatile constituents in ANP were described comprehensively and systematically for the first time. A total of 93 volatile components, assigned to four aromatic Chinese medicine, were preliminarily identified in ANP. Among them, 21 volatile components were selected and accurately quantified within 15 min based on the established GC-MS/MS method, in which d-borneol, isoborneol and muscone were the top three abundant ingredients in the ANP. According to the contents of 21 determined constituents, 8 batches of ANP were generally chemical similar based on PCC analysis. In addition, d-borneol, isoborneol and muscone and one metabolite (camphor) were found in rat plasma after single oral administration, suggesting that these prototypes and their metabolites could be potential bioactive substances of ANP. Collectively, this study provides a comprehensive information of volatile component of ANP in vitro and in vivo, which may helpful for the quality evaluation and pharmacological research of ANP. Meanwhile, the non-volatile constituents were also important for its holistic efficacy, which need further investigated using the appropriate analytical method.
Additional file 1: Table S1. The information on eight batches of commercially available ANP samples. Table S2. Dynamic multiple reaction monitoring parameters of all analytes and internal standard. Table S3. Calibration curves, LODs and LOQs of 21 volatile analytes in ANP samples. Table S4. Intra-day precision, inter-day precision, repeatability, stability and recovery of 21 volatile analytes in ANP samples. Table S5. Contents of 21 volatile analytes in ANP samples (μg/g, mean ± SD, n = 4).
Additional file 2: Figure S1. The sources of volatile constituents in ANP samples. Figure S2. The concentration heat map of 21 volatile analytes in ANP samples from eight batches.