- Research
- Open access
- Published:
Metabolite identification of salvianolic acid A in rat using post collision-induced dissociation energy-resolved mass spectrometry
Chinese Medicine volume 19, Article number: 64 (2024)
Abstract
Background
As one of the most famous natural products, salvianolic acid A (SAA) is undergoing clinical trials for the treatments of angina pectoris and coronary heart disorders. However, the in vivo metabolites of SAA have only been tentatively identified, leading to a barrier for precise therapeutical drug monitoring.
Methods
Ultra-high performance liquid chromatography coupled with quadrupole time of flight tandem mass spectrometry (UPLC–Qtof-MS/MS) was firstly employed to acquire high-resolution MS1 and MS2 spectra for all metabolites. Through paying special attention onto the features of ester bond dissociation, metabolism sites were restricted at certain regions. To further determine the metabolism site, such as the monomethylated products (M23, M25, and M26), post collision-induced dissociation energy-resolved mass spectrometry (post-CID ER-MS) was proposed through programming progressive exciting energies to the second collision chamber of hybrid triple quadrupole-linear ion trap mass spectrometry (Qtrap-MS) device.
Results
After SAA oral administration, 29 metabolites (M1–M29), including five, thirteen, and sixteen ones in rat plasma, urine, and feces, respectively, were detected in rats. The metabolism route was initially determined by applying well-defined mass fragmentation pathways to those HR-m/z values of precursor and fragment ions. Metabolism site was limited to SAF- or DSS-unit based on the fragmentation patterns of ester functional group. Through matching the dissociation trajectories of concerned 1st-generation fragment ions with expected decomposition product anions using post-CID ER-MS strategy, M23 and M25 were unequivocally assigned as 3'-methyl-SAA and 3''-methyl-SAA, and M26 was identified as 2-methyl-SAA or 3-methyl-SAA. Hydrolysis, methylation, glucuronidation, sulfation, and oxidation were the primary metabolism channels being responsible for the metabolites' generation.
Conclusion
Together, the metabolism regions and sites of SAA metabolites were sequentially identified based on the ester bond dissociation features and post-CID ER-MS strategy. Importantly, the present study provided a promising way to elevate the structural identification confidence of natural products and metabolites.
Graphical abstract
Background
Salvianolic acid A (SAA), sourced from Miltiorrhizae Radix et. Rhizoma (Chinese name: Danshen), is reputed as one of the most famous natural products [1]. As the primary contributor for the therapeutic benefit of Danshen, SAA exhibits a broad pharmacological spectrum, including anti-tumor effects [2, 3], myocardial protection [4, 5], hepatoprotective activity [6, 7], and modulation of gut microbiota imbalance [8]. It is worth noting that a new drug candidate derived from SAA is currently undergoing phase II clinical trials for the treatments of angina pectoris and coronary heart disorders. Because of the great absorption property and nonetheless, low oral bioavailability, the metabolites of SAA are believed to be the active forms in vivo and responsible for the therapeutic outcomes via interactions with drug targets [9, 10]. To our knowledge, the metabolism of SAA in vivo remains insufficiently investigated in comparison to its well-defined pharmacological effects. Therefore, in-depth clarification of SAA metabolism is a prerequisite to disclose the underlying therapeutic mechanisms and facilitate the development of potentially active drug candidates from metabolites.
SAA (Fig. 1) contains a total of six hydroxy moieties and one carboxyl moiety, all of which serve as the metabolism labile sites [11,12,13]. Notably, various metabolism sites lead to the widespread presence of isomers, resulting in that the confirmative annotation is crucial and challenging. Using nuclear magnetic resonance (NMR), confirmative identification of isomers can be achieved [9], but the application is constrained by limited sensitivity and labor- and time-consuming shortcomings, as well as risks of metabolite lability. LC–HR-MS/MS enables the determination of metabolic type by applying mass fragmentation pathways to those HR-m/z values obtained from MS1 and MS2 spectra [14,15,16]. However, definitive identification of SAA metabolites is difficult to achieve, because the isomers usually share identical MS1 and indistinguishable MS2 spectra. For example, methylated SAA were ambiguously assigned without specifying the substitution site [17, 18]. Online energy-resolved (ER)-MS [19, 20] and full collision energy ramp-MS2 (FCER-MS2) spectrum [21,22,23,24,25] have been proposed to investigate the differential MS/MS behaviors among isomers. In addition to the m/z values of precursor and fragment ions, the FCER-MS2 spectrum contains orthogonal information, such as the collision energy at 50% survival yield (CE50) of the precursor ion, optimal collision energy (OCE), and maximum relative ion intensity at OCE (RIImax) of each fragment ion. Combined with quantum chemical calculations, OCE and RIImax of identical ion transitions for isomers are positively and negatively correlated with bond dissociation energy (BDE) and thermal enthalpy, respectively. Isomeric differentiation relies on aligning the ranking patterns of OCEs and RIImax of detected compounds with the corresponding BDEs and thermal enthalpy values of the candidate structures. Unfortunately, not all candidate structures can be detected simultaneously in practical samples, especially when extensive substitution sites exist, which means that FCER-MS2 strategy still has certain limitations in confirmatively identifying isomeric structures. Therefore, it is necessary to develop novel strategies for acquiring auxiliary structural information and advancing the identification of SAA metabolites in vivo.
The difficulty of structural identification can be reduced by hydrolyzing conjugated compound into decomposition products (i.e., biosynthetic pioneers) and then comparing them with appropriate authentic simple compounds [26, 27]. We aimed to implement a similar concept in MS with improved sensitivity, specificity, and throughput [28]. The hybrid triple quadrupole-linear ion trap mass spectrometry (Qtrap-MS) is equipped with two tandem-in-space collision cells, i.e., q2 collision cell and linear ion trap (LIT) chamber, and the latter one allows for MS3 spectra acquisition. Inspired by that FCER-MS2 spectrum contains trajectories from precursor ion to 1st-generation fragment ions, we proposed post collision-induced dissociation energy-resolved MS (post-CID ER-MS) to obtain trajectories from selected 1st-generation fragment ion to 2nd-generation fragment ions [29, 30]. Following the post-CID ER-MS strategy, intensities of residual 1st-generation fragment ion and all its 2nd-generation fragment ions were collected at different exciting energies (EEs). After normalization and curve fitting, full exciting energy ramp-MS3 (FEER-MS3) spectrum of selected 1st-generation fragment ion was obtained, which included exciting energy at 50% survival yield (EE50) of 1st-generation fragment ion, optimal exciting energy (OEE), and maximum relative ion intensity (RIIOEE) of each 2nd-generation fragment ion [31, 32]. Theoretically, FEER-MS3 spectrum of 1st-generation fragment ion for conjugated structure should be same as FEER-MS2 spectrum of the (de)protonated ion for decomposition product bearing the expectational geometry, and vice versa.
SAA is the esterification product of salvianolic acid F (SAF) and danshensu (DSS), and similarly, hydrolysis of ester bond can yield the two biosynthetic pioneers, SAF and DSS (Fig. 1). Fortunately, the ester bond of SAA is also susceptible to cleavage in collision-induced dissociation (CID), resulting in a series of related fragment ions, including a−, b−, c−, x−, y−, and z− ions according to the nomenclature rule [22]. Therefore, based on the ester bond dissociation features and post-CID ER-MS strategy, we proposed an analytical pipeline aiming to enhance the structural identification confidence and clarify the metabolic profile of SAA, following several progressive steps: (1) capturing HR-MS/MS information of SAA to analyze the features of ester bond dissociation; (2) obtaining HR-MS/MS information of SAA metabolites to determine the metabolism route and region (SAF- or DSS-unit) based on the mass fragmentation pathways and the features of ester bond dissociation, respectively; (3) acquiring FEER-MS3 spectra of targeted 1st-generation fragment ions for monomethylated SAA isomers using post-CID ER-MS strategy and comparing them with FEER-MS2 spectra of monomethylated SAF/DSS anions, to confirmatively identify the structures of monomethylated metabolites; and (4) proposing the possible metabolism channels for SAA. The obtained findings are expected to acquire a comprehensive metabolism map of SAA, and also demonstrate the significance of post-CID ER-MS strategy for elevating the structural annotation confidence.
Materials and methods
Chemicals and reagents
SAA, sodium DSS, and SAF were commercially supplied by Standard Technology Co., Ltd. (Shanghai, China). The purity of each compound was determined to be greater than 98% using LC–UV–IT-TOF–MS (Shimadzu, Kyoto, Japan). LC–MS grade formic acid, methanol, and acetonitrile (ACN) were purchased from Thermo-Fisher (Pittsburgh, PA). Deionized water was prepared in-house on a Milli-Q water purification instrument (Millipore, Bedford, MA). The other chemicals and reagents were of analytical grade and purchased from Beijing Chemical Co., Ltd. (Beijing, China).
Sample collection and preparation
The animal study protocols were approved by the Animal Ethics Committee of Beijing University of Chinese Medicine (Beijing, China). Twenty-four male Sprague–Dawley rats (220 ± 20 g) were provided by Vital River Laboratory Animal Technology Co., Ltd. (Beijing, China). All animals were acclimated in an animal breeding room for 3 days with the temperature of 22 ± 2 °C, a relative humidity of 60 ± 5% and a 12 h light/dark cycle. Standard chow and Milli-Q water were provided ad libitum. All animals were fasted overnight and given free access to tap water prior to treatment. All rats were randomly divided into four groups (n = 6) that orally administrated (100 mg/kg for each compound) with SAA, sodium DSS, SAF, and vehicle (tap water), accordingly. The animals were individually housed in metabolic cages.
A 200 μL aliquot of blood was sampled from the orbital venous plexus of each rat at 0.5, 1, 2, 6, 12, and 24 h after oral administration. Urine and feces were collected individually from metabolic cages over 0–24 h in an ice water bath. Drug-free biological samples were collected from vehicle group in parallel. Plasma and urine were centrifuged under 3 000 rpm and 12 000 rpm, respectively, for 10 min at 4 °C, and supernatants were retained. Feces underwent lyophilization and grinding, successively. Thereafter, all samples were stored at −80 °C until usage.
After being thawed in ice-water bath, each sample type was pooled within each group to generate a single homogenized sample. Three volumes of ACN were fortified into 300 μL homogenized plasma or urine samples for protein precipitation and desalting. Ten volumes of methanol were utilized to extract pulverized feces with an ultrasonic-assisted manner. The supernatants of plasma samples were evaporated to dryness at 30 °C under gentle nitrogen stream, and the residues were reconstituted with 50 μL of 25% aqueous ACN, vortexed for 2 min, and centrifugated at 12 000 rpm for 15 min at 4 °C, successively. Following thorough vortex and centrifugation, the supernatants of urine and feces were individually filtered with 0.22 μm nylon membranes.
LC–MS/MS measurements
All biological samples participated in LC–MS/MS measurements. All chromatographic separations were conducted on Shimadzu LC modular system (Kyoto, Japan) equipped with a Waters ACQUITY UPLC HSS T3 column (2.1 × 100 mm, 1.8 μm, Milford, CT). Mobile phase consisting of 0.1% aqueous formic acid (A) and ACN (B) was delivered in gradient at a total flow rate of 0.2 mL/min by operating the program as follows: 0–4 min, 2–10% B; 4–20 min, 10–30% B; 20–25 min, 30–65% B; 25–28 min, 65–100% B; and 28–32 min, 100% B. After each analytical run, 2% B was delivered for another five minutes for the re-equilibration. Column oven and injection volume were set as 40 °C and 2 μL, respectively.
The column outlet was directly connected to ESI interface of either Qtof-MS (TripleTOF 6600+, SCIEX, Foster City, CA) or Qtrap-MS (5500 Qtrap-MS, SCIEX) device. Actually, the whole measurement consisted of three progressive steps: UPLC–Qtof-MS/MS for acquiring HR-m/z values of precursor and fragment ions, online ER-MS for obtaining OCEs of the targeted 1st-generation fragment ions, and post-CID ER-MS for gaining the trajectories from selected 1st-generation fragment ions to 2nd-generation fragment ions.
Qtof-MS assays
The ion source parameters were set as follows: ion-spray needle voltage (IS), −4500 V; nebulizer (GS1), heater (GS2), and curtain gas flow rate, 55, 55, and 30 units, respectively; heater gas temperature (TEM), 500 °C. Mass range as 100–1000 Da was defined to record MS1 and MS2 spectra. MS2 spectra were acquired by data-dependent acquisition (DDA) manner with CE as −30 eV, collision energy spread (CES) as 15 eV. SCIEX PeakView1.2 software was utilized for data processing.
Qtrap-MS assays
The online ER-MS and post-CID ER-MS were configured using MRM and MS3 scan modes of Qtrap-MS. A methylation product namely M25 was taken as a case to clarify the detailed settings. Based on the fragmentation patterns of ester functional group in CID [33], the fragment ions at m/z 313.07, 295.06, 211.06, and 193.05 were assigned as c−, b−, y−, and z−, respectively, in MS2 spectrum of m/z 507.13 ([M–H]−) (Table 1). The c− and b− ions had identical m/z values with [MSAF–H]− and [MSAF–H–H2O]−, while y− and z− ions showed an increase of 14 Da (CH2 fortification) compared to [MDSS–H]− and [MDSS–H–H2O]−. Hence, it could be deduced that the methylation site was located at DSS-unit. Because our concern was the comparability between certain fragment ions and the deprotonated ions of those decomposition products, subsequent experiments focused only on the m/z 313.0676 (c−) and 211.0638 (y−) ions. The data recording and processing were performed using SCIEX Analyst software (Version 1.6.3).
The previously described online ER-MS [21] was strictly followed to gain OCEs for the concerned MS2 fragment ions (c− and y− ions). For the ion transition m/z 507.13 > 211.06, a set of pseudo-ion transitions (PITs) was derived, such as m/z 507.100 > 211.100, 507.100 > 211.101, 507.100 > 211.102, etc., corresponding to progressive CEs ranging from −5 to −81 eV. The responses under different CEs were imported into GraphPad Prism 8.0 software (San Diego, CA) to plot Gaussian-shaped breakdown graphs. The obtained OCE was used for subsequent MS3 data acquisition.
The follow-up post-CID ER-MS strategy in MS3 mode involved two CID, which occurred sequentially in q2 and LIT, corresponding to energies called CE and EE, respectively. OCE was set in q2 to accumulate enough fragment ions for MS3 acquisition. To distinguish fragment ions generated in q2 and LIT, we referred to them as 1st-generation fragment ions and 2nd-generation fragment ions, respectively [34]. Precursor ion, concerned 1st-generation fragment ion, and 2nd-generation fragment ions of each compound were involved to construct Q1 > Q3 > QLIT ion transitions [29,30,31]. Further fragmentation trajectories of m/z 211.06 for M25 were obtained at progressive absolute values of EE ranging from 0.05 to 0.30 V in LIT. The breakdown graphs of residual 1st-generation fragment ion (m/z 211.06) and 2nd-generation ions (m/z 211.06, 193.00, 178.00, 165.00, 150.00, 149.10, and 134.00) were fitted using sigmoid plot and multiple iterations with OriginPro9 software (OriginLab, Northampton, MA), respectively. Thereafter, breakdown graphs were assembled and normalized to obtain the FEER-MS3 spectrum of m/z 211.06 for M25, which implied molecular descriptors (MDs), including EE50, OEE, and RIIOEE.
The related decomposition products, monomethylated DSS (D14 and D15), which had the same m/z value as the y− of monomethylated SAA (M25), were also involved in post-CID ER-MS measurements. The ion transition m/z 211.06 > 211.06 > QLIT of D14 and D15 was employed to produce breakdown graphs. Under such status, the MS2 spectrum generated through CID in the LIT cell and the resultant breakdown graph set was termed as FEER-MS2 spectrum [31].
Results and discussion
Mass fragmentation behaviors of SAA
The fragmentation behaviors of SAA were investigated by employing UPLC–Qtof-MS/MS. The deprotonated ion ([M–H]−) of SAA was observed at m/z 493.11, corresponding to the elemental composition as C26H22O10. After converting HR-m/z to elemental compositions, the fragment ions in MS2 spectrum were assigned as m/z 313.07 [M–H–C9H8O4]−, 295.06 [M–H–C9H10O5]−, 269.08 [M–H–C10H8O6]−, 203.04 [M–H–C9H8O4–C6H6O2]−, 197.05 [M–H–C17H12O5]−, 185.02 [M–H–C9H10O5–C6H6O2]−, 179.03 [M–H–C17H14O6]−, 159.05 [M–H–C10H8O6–C6H6O2]−, 135.04 [M–H–C18H14O8]−, and 109.03 [M–H–C20H16O8]− (Additional file 1: Fig. S1A). Referring to the chemical structure of SAA, the dissociations around the ester bond were regarded as the primary contributions towards the entire mass fragmentation pattern. Similar to the nomenclature rule of peptides, a nomenclature rule was proposed for fragment ions related to the ester bond dissociation, including a−, b−, c−, x−, y−, and z− [22]. Based on the mathematical relationships between complementary ions (such as a− and x−, b− and y−, as well as c− and z−), along with the unique mass differences between a− and b− (CO, 25.98 Da), b− and c− (H2O, 18.01 Da) [33], the featured ions m/z 269.08, 295.06, 313.07, 197.05, and 179.03 were assigned as a−, b−, c−, y−, and z−, respectively. The fragmentation pathways of SAA are proposed in Additional file 1: Fig. S1B, which provide references for the follow-up identification of in vivo metabolites.
Metabolite identification for SAA by applying ester bond dissociation pattern
After carefully comparing each homogenized biological sample with the respective drug-free sample, a total of 29 metabolites (M1–M29) were detected. Thereof, five ones (M17, M18, M20, M22, and M24), thirteen ones (M2–M5, M7, M11, M12, M14, M16, M25, M26, M28, and M29), and sixteen ones (M1, M5–M10, M13, M15, M19, M21, M23, and M26–M29) were distributed in plasma, urine, and feces, respectively. The detailed chromatographic and mass spectrometric information is summarized in Table 1, and representative chromatograms are illustrated in Fig. 2.
Among the metabolites, M1–M4 were derived from hydrolysis of SAA and identified as DSS, SAF glucuronide, hydrogenated SAF, and SAF, respectively. M1 and M4 were identified by comparing with literature and reference substances [35, 36]. The M3 anion (tR, 19.73 min) existed at m/z 315.09 (C17H16O6), which was 2 Da greater than SAF anion. Additionally, typical fragment ions were observed at m/z 271.10 [M–H–CO2]−, 255.06 [M–H–C2H4O2]−, 135.05 [M–H–C9H8O4]−, and 109.02 [M–H–C11H10O4]−; therefore, M3 was finally identified as hydrogenated SAF. The deprotonated ion ([M–H]−) of M2 (tR, 18.15 min) was observed at m/z 489.11, corresponding to a molecular formula as C23H22O12, which was 176 Da (C6H8O6) greater than SAF anion. The fragment ions observed after a neutral loss of 176 Da were consistent with those of SAF, confirming the identification of M2 as a glucuronidation product of SAF.
All other metabolites (M5–M29) were direct derivatives of SAA, and the metabolism regions, either SAF-unit or DSS-unit, were judged based on MS2 spectral information and the ester bond dissociation pattern. The presence of substitution at SAF-unit was determined by calculating the mass difference between c−, b−, and a− ions of metabolites and the corresponding ions of SAA, and similarly, the substitution site at DSS-unit was determined by the mass difference between z− and y− ions of metabolites and their relevant ions of SAA.
Deprotonated ions of M5 (tR, 20.54 min), M7 (tR, 21.52 min), and M8 (tR, 21.94 min) were observed at m/z 669.15, suggesting the molecular formula as C32H30O16 for each. Compared to SAA (molecular formula: C26H22O10), glucuronidation should be responsible for GluA (C6H8O6, 176 Da) fortification. In addition, the signals at m/z 493.11 ([M–H–GluA]−) in MS2 spectra of M5, M7, and M8 further consolidated that these were glucuronidation products of SAA. Subsequent efforts were towards the identification of metabolism regions based on the pattern of ester bond dissociation. The b− (m/z 471.10) of M5 and M7 exhibited 176 Da greater than b− (m/z 295.06) of SAA, while c− (m/z 489.10) and a− (m/z 445.10) ions of M8 also showed a similar increase by 176 Da compared to c− (m/z 313.07) and a− (m/z 269.08) ions of SAA, indicating that glucuronidation occurred at the SAF-unit ([SAF-glucuronide]-[DSS]) for all three ones.
M10 (tR, 22.55 min), M12 (tR, 23.14 min), and M13 (tR, 23.32 min) were all assigned as SAA sulfate (molecular formula: C26H22O13S) attributing to the signals at m/z 573.07 ([M–H]−) in MS1 spectra and m/z 493.11 ([M–H–SO3]−) in MS2 spectra. The c− (m/z 393.04 [M–H–C9H8O4]−) of M10 and M13 showed 80 Da greater than c− of SAA, indicating that sulfation sites might located at the SAF-unit ([SAF-sulfate]-[DSS]). Nevertheless, because the MS2 spectrum of M12 contained few characteristic ions, the conjugation region was not confirmed in the current study.
Three isomers, M23 (tR, 25.61 min), M25 (tR, 25.78 min), and M26 (tR, 25.96 min), were all annotated as the methylation products of SAA, because their deprotonated ions at m/z 507.13 ([M–H]−, C27H24O10) were 14 Da (CH2) greater than SAA anion. Thereof, the featured ions c− (m/z 327.09 [M–H–C9H8O4]−) and b− (m/z 309.07 [M–H–C9H10O5]−) of M23 and M26 exhibited a mass increase of 14 Da compared to c− and b− ions of SAA, respectively. Meanwhile, the featured ions y− (m/z 197.05 [M–H–C18H14O5]−) and z− (m/z 179.03 [M–H–C18H16O6]−) of M23 and M26 were consistent with those observed in SAA. The above demonstrated that methylation region was located at SAF-unit ([methyl-SAF]-[DSS]). In addition, m/z 294.06 ([M–H–C9H10O5–CH3·]−·) of M26 was generated by methyl radical cleavage (CH3·) of b−, further confirming the deduction about methylation region. While y− (m/z 211.06) and z− (m/z 193.05) of M25 were found to be 14 Da greater than their counterparts of SAA, this proved that its methylation site was at the DSS-unit ([SAF]-[methyl-DSS]).
M11 (tR, 22.75 min) and M14 (tR, 23.68 min) shared identical precursor ions at m/z 683.16, and the molecular formula was consequently calculated as C33H32O16, implying that either could be the product of methylation plus glucuronidation of SAA. The fragment ions at m/z 507.10 in MS2 spectra were generated by the natural loss of 176 Da. M11 and M14 possessed almost identical fragment ions to M5 ([SAF-glucuronide]-[DSS]), including m/z 471.10 (b−), 313.07 (c−–176 Da), and 295.06 (b−–176 Da), suggesting that glucuronidation site of M11 and M14 was located at the SAF-unit. Subsequently, based on the complementarity of b− and y− ion, the mass of y− was presumed to be m/z 211.06, which was consistent with [M–H]− of monomethylated DSS, indicating that the methylation site was at the DSS-unit. Together, M11 and M14 were identified as [methyl-DSS]-[SAF-glucuronide].
In the case of a pair of isomers, M28 (tR, 26.73 min) and M29 (tR, 26.89 min), their empirical formulas were all suggested as C28H26O10 based on the deprotonated ions at m/z 521.14 ([M–H]−). Compared to SAA, di-methylation should be responsible for two CH2 (28 Da) fortification. The MS2 spectra of M28 and M29 exhibited an identical set of fragment ions, including c− (m/z 327.09), b− (m/z 309.08), y− (m/z 211.06), and z− (m/z 193.05). Thereof, c− and b− were identical with c− and b− of M23 and M26 ([methyl-SAF]-[DSS]), respectively; meanwhile, y− and z− corresponded to those of M25 ([SAF]-[methyl-DSS]). Hence, two methylation sites separately occurred at either unit of SAA, and finally, M28 and M29 were termed as [methyl-SAF]-[methyl-DSS].
The deprotonated ions ([M–H]–) of M16 (tR, 24.66 min), M17 (tR, 24.68 min), M18 (tR, 25.07 min), M20 (tR, 25.42 min), M22 (tR, 25.54 min), and M24 (tR, 25.70 min) occurred at m/z 697.18 (molecular formula: C34H34O16), indicating that they were resulted from fortifying two CH2 together with a glucuronyl group to SAA. In MS2 spectra of M16–M18, M20, and M24, the product ions at m/z 327.09, 309.08, 294.06, 211.07, and 193.05 corresponded to the featured ions of the di-methylation metabolites of SAA described above. Therefore, two methylation sites separately located at SAF-unit and DSS-unit, i.e., [methyl-SAF]-[methyl-DSS]-glucuronide. However, in the MS2 spectrum of M22, there was not enough characteristic fragment ions to support the determination of the methylation sites; therefore, M22 was tentatively identified as di-methyl-SAA glucuronide. Unfortunately, neutral loss of GluA was observed prior to ester bond dissociation in CID, resulting in the absence of featured fragment ions, and as a consequence, it was challenging to determine the glucuronidation site.
The deprotonated molecular ion ([M–H]–) of M21 (tR, 25.52 min) was observed at m/z 477.12 (C26H22O9), which was 16 Da less than SAA. Other featured ions in MS2 spectrum, such as c− (m/z 313.07) and b− (m/z 295.06) were consistent with the fragment ions of the prototype. Therefore, M21 was identified as the dehydroxylated derivative of SAA and dehydroxylation occurred on the DSS-unit ([SAF]-[DSS-dehydroxylated]).
M19 (tR, 25.32 min) possessed a [M–H]− signal as m/z 491.10 (molecular formula: C26H20O10), 2 Da less than SAA, indicating it might be dehydrogenation product of SAA. The metabolism site was determined to be at the SAF-unit by comparing the featured fragment ions of M19 and SAA, including c− (m/z 311.05 vs. 313.07), b− (m/z 293.05 vs. 295.06), a− (m/z 267.06 vs. 269.08), y− (m/z 197.05 vs. 197.05), and z− (m/z 179.04 vs. 179.04). Through referring the descriptions in the literature [17, 36, 37], M19 was definitely determined as salvianolic acid C (SAC).
M6 (tR, 20.91 min), M9 (tR, 22.42 min), and M15 (tR, 24.06 min) shared an identical [M–H]− ion at m/z 667.13 (C32H28O16), which was 176 Da greater than that of SAC, and their generation should be attribute to the glucuronidation for SAC. Meanwhile, in MS2 spectra of M6 and M15, the featured ions including b− (m/z 469.07) and c− (m/z 487.10), were 176 Da greater than the b− (m/z 293.05) and c− (m/z 311.05) of SAC, respectively, which reminded that the glucuronidation sites might locate at the SAF-unit ([tournefolic acid A-glucuronide]-[DSS]). Nevertheless, the MS2 spectrum of M9 displayed a characteristic ion of y− (m/z 355.07), indicating that the conjugation site on the DSS-unit ([tournefolic acid A]-[DSS-glucuronide]).
The [M–H]− ion of M27 was observed at m/z 475.10 with a retention time of 26.42 min, and the molecular formula thereafter came out as C26H20O9. M27 was 16 Da less than SAC, and as a result, dehydroxylation should be responsible for the metabolic process. Moreover, fragment ions c− (m/z 311.07) and b− (m/z 293.04) were detected in the MS2 spectrum, corresponding to c− and b− of SAC, implying that the metabolism region occurred at the DSS-unit. Consequently, M27 was tentatively identified as dehydroxylated SAC ([tournefolic acid A]-[DSS-dehydroxylated]).
Isomeric identification of methylation metabolites
Through analyzing HR-MS/MS spectral information, the metabolites of SAA in vivo were tentatively annotated. Thereafter, based on the ester bond dissociation features, one step forward was achieved through placing the metabolism site(s) onto SAF-unit and/or DSS-unit. To further advance structural identification confidence, three isomeric metabolites of methylation, such as M23, M25, and M26, were employed as representatives to elaborate structural identification details. Noteworthily, to aid structural identification for SAA metabolites, the metabolites of SAF (F1–F30) and DSS (D1–D15) were first characterized carefully through matching with the information archived in the literature [38,39,40]. Attributing to the relatively simple structures, confirmative identification was usually achieved, and their identities together with key MS/MS spectral information were summarized in Additional file 1: Tables S1, S2.
According to the annotation of SAA methylation region mentioned above, M23 and M26 occurred at the SAF-unit, while M25 took place at the DSS-unit, but the specific metabolism site could not be determined, e.g., either the C-3-OH or C-4-OH of DSS. Based on the fact that the c− and y− ions produced by ester bond dissociation should show identical structures to those of the hydrolysis product anions, we thereby employed the post-CID ER-MS strategy. The strategy was performed to determine the structures of the c− and y− ions by comparing the FEER-MS3 spectra of c− and y− ions for esters with FEER-MS2 spectra of [M–H]– for candidate hydrolysis products. Therefore, the methylation metabolites (F28 and F29, as well as D14 and D15) of both SAF and DSS were involved.
Through programming online ER-MS to the constructed two ion transitions (m/z 507.13 > 211.06 for M25 and 507.13 > 327.09 for M23 and M26), the OCEs were −38.2 eV and −29.1 eV, respectively. When constructing Q1 > Q3 > QLIT ion transitions, OCE was defined in q2, and progressive EEs were implemented in LIT chamber to acquire the FEER-MS3 spectra. Meanwhile, [M–H]– of each candidate hydrolysis product was permitted to pass q2 cell without dissociation through defining low CE as −10 eV, and thereby acquiring the FEER-MS2 spectra of [M–H]– for candidate hydrolysis products.
FEER-MS3 spectra of m/z 507.13 > 211.06 for M25 together with m/z 211.06 > 211.06 for D14 and D15 are illustrated in Fig. 3A. FEER-MS3 spectrum of y− (m/z 211) for M25 involved the breakdown graphs of m/z 211 (EE50, −0.05626 V), 193 (OEE, −0.06746 V; RIIOEE, 100%), 178 (−0.07292 V; 3.87%), 165 (−0.07239 V; 12.46%), 150 (−0.07918 V; 10.25%), 149 (−0.07172 V; 7.51%), and 134 (−0.07870 V; 43.54%). The EE50 and OEEs of M25 did not exhibited consistency to D14 and D15, while RIIOEE values demonstrated the similarities and the differences. The RIIOEE values for each 2nd-generation fragment ion of M25 were consistent with those of D14 anion (m/z 193, 100%; m/z 178, 3.57%; m/z 165, 13.22%; m/z 150, 10.08%; m/z 149, 7.48% and m/z 134, 38.30%), but did not exhibit such agreement with the D15 anion (m/z 193, 100%; m/z 178, 10.29%; m/z 165, 1.61%; m/z 150, 4.77%; m/z 149, 8.42% and m/z 134, 32.97%). The detailed information and the RIIOEE-trend curves of M25, D14, and D15 are shown in Additional file 1: Table S3 and Fig. S2. The FEER-MS3 spectrum of y− for M25 matched well with the FEER-MS2 spectrum for D14 (3-methyl-DSS) anion. In conclusion, M25 was configured as 3″-methyl-SAA.
In the same vein, M23, M26, F28, and F29 were also assayed using post-CID ER-MS strategy to produce FEER-MS3 and FEER-MS2 spectra (Fig. 3B). Great consistence occurred between FEER-MS3 spectrum of c− for M23 and FEER-MS2 spectrum for F28 anion, either of which was made up of the post-CID breakdown graphs of m/z 327 (EE50, −0.04542 V for M23 and −0.04541 V for F28), 268 (OEE and RIIOEE, −0.06916 V and 91.51% for M23 vs. −0.05896 V and 55.66% for F28), 159 (−0.05461 V and 100% for M23 vs. −0.04939 V and 100% for F28), 146 (−0.06836 V and 6.75% for M23 vs. −0.06971 V and 12.86% for F28), 123 (−0.04614 V and 26.11% for M23 vs. −0.04619 V and 26.12% for F28) and 108 (−0.06773 V and 18.58% for M23 vs. −0.06082 V and 25.38% for F28). Similarly, there exhibited a significant similarity between FEER-MS3 spectrum of y− for M26 and FEER-MS2 spectrum for F29 anion, either of which was composed of breakdown graphs of m/z 327 (−0.03645 V for M26 and − 0.03427 V for F29), 268 (−0.05664 V and 100% for M26 vs. −0.05594 V and 100% for F29), 159 (−0.06762 V and 8.71% for M26 vs. −0.07407 V and 6.00% for F29, 146 (−0.06419 V and 29.29% for M26 vs. −0.07071 V and 30.10% for F29), 123 (−0.05786 V and 0.08% for M26 vs. −0.07866 V and 0.60% for F29) and 108 (−0.04591 V and 1.28% for M26 vs. −0.04017 V and 1.08% for F29). The information above is listed in Additional file 1: Table S4 and trend curves of RIIOEE are shown in Additional file 1: Fig. S3. The OEE values of M23, M26, F28, and F29 exhibited no discernible distinction, but EE50 value for M23 closely resembled that of F28, while M26 consistently aligned with the F29. Finally, M23 was tentatively identified as 3'-methyl-SAA and M26 as either 2-methyl-SAA or 3-methyl-SAA.
Above all, FEER-MS3 and FEER-MS2 spectra matching is meaningful for structural annotation because RIIOEE feature is eligible to capturing diagnostic fragment ions and post CID ER-MS should be a credible strategy to identified fragment ions.
Noteworthily, the post-CID ER-MS strategy could be applicable for identifying any other metabolites, not limited to monomethylated metabolites. In theory, the acquisition of the FEER-MS3 spectrum of the concerned 1st-generation fragment ion can be achieved by employing the post-CID ER-MS strategy whenever the c–and z– ions generated from ester bond dissociation are detectable in the MS2 spectrum. However, glucuronidated and sulfonated substituents are susceptible to neutral loss under CID, and when they are fragmented preferentially over the ester bond, c− or y− ions are unavailable to acquire FEER-MS3 spectra. Subsequent attempts can be made to retain the glucuronosyl and sulfonyl by altering to other dissociation patterns, such as electron capture dissociation (ECD) and electron transfer dissociation (ETD). To some extent, it is extremely challenging to elucidate the exact structures of all metabolites for SAA. Therefore, the current study focused solely on discussing and analyzing methylation as a typical example, which effectively demonstrated the feasibility and practicality of this approach called FEER-MS3 and FEER-MS2 spectra matching.
Metabolism channels of SAA in rats
After the confidence-enhanced structural annotation for all metabolites (M1–M29), the metabolic pathways for SAA were proposed (Fig. 4). SAA underwent ester bond hydrolysis and further metabolism to produce M1–M4. SAA underwent glucuronidation, sulfation, methylation, or a combination of these to generate a series of metabolites (M5, M7, M8, M10–M14, M16–M18, M20, M22–M26, M28, and M29). SAA could directly undergo dihydroxylation to produce M21. SAA was able to give birth to SAC (M19) via oxidation, and subsequently to generate M27, M6, M9, and M15 through dehydroxylation and glucuronidation. The absence of SAA in plasma indicated that extensive metabolism occurred. All five metabolites detected in plasma were di-methylation plus glucuronidation products. The lower number of metabolites in rat plasma may be attributed to the lower bioconversion activity or the significantly low oral bioavailability and short plasma elimination half-life of SAA in rats [9]. SAC resulting from SAA dehydrogenation and further metabolites were only detected in feces, but not excreted via urine, presumably because feces and urine were closely related biliary and renal excretion, respectively.
Conclusion
Confirmative identification of SAA metabolites in vivo is highly challenging due to the variety of possible metabolism sites. In the present study, one step forward was achieved for structural identification of most metabolites by utilizing the features of ester bond dissociation and the mass differences between featured ions. Post-CID ER-MS strategy was proposed to specify the metabolism sites. And fortunately, the metabolism sites of monomethylated SAA (M23, M25, and M26) were confirmatively determined by matching FEER-MS3 spectra of 1st-generation fragment ions with FEER-MS2 spectra of monomethylated SAF and monomethylated DSS (biosynthetic pioneers). As a result, a total of 29 metabolites of SAA were preliminarily identified in rats. The primary metabolism channels of SAA encompassed hydrolysis, methylation, glucuronidation, sulfation, and oxidation. Overall, this study provides substantial evidence for comprehensively understanding the various metabolic forms of SAA in vivo after oral administration, and the developed strategy should be able to strengthen the identification confidence of complex compounds.
Availability of data and materials
All the data used to support the findings of this study are available from the corresponding author upon reasonable request.
Abbreviations
- CE50 :
-
The collision energy (CE) at 50% survival yield
- CID:
-
Collision-induced dissociation
- DSS:
-
Danshensu
- EE50 :
-
The exciting energy (EE) at 50% survival yield
- EEs:
-
Exciting energies
- FCER-MS2 :
-
Full collision energy ramp-MS2
- FEER-MS3 :
-
Full exciting energy ramp-MS3
- GluA:
-
Glucuronidation
- LIT:
-
Linear ion trap
- MDs:
-
Molecular descriptors
- OCE:
-
Optimal collision energy
- OEE:
-
Optimal exciting energy
- Online ER-MS:
-
Online energy-resolved mass spectrometry
- PITs:
-
Pseudo-ion transitions
- Post-CID ER-MS:
-
Post collision-induced dissociation energy-resolved mass spectrometry
- RIImax :
-
Maximum relative ion intensity at OCE
- RIIOE E :
-
Maximum relative ion intensity at OEE
- SAA:
-
Salvianolic acid A
- SAF:
-
Salvianolic acid F
References
Li L, Tan R, Chen W. Salvianolic acid A, a new depside from roots of Salvia miltiorrhiza. Planta Med. 1984;50(3):227–8.
Pei R, Si T, Lu Y, et al. Salvianolic acid A, a novel PI3K/Akt inhibitor, induces cell apoptosis and suppresses tumor growth in acute myeloid leukemia. Leuk Lymphoma. 2018;59(8):1959–67.
Chuang C, Ho Y, Lin C, et al. Salvianolic acid A suppresses MMP-2 expression and restrains cancer cell invasion through ERK signaling in human nasopharyngeal carcinoma. J Ethnopharmacol. 2020;252:112601.
Gong D, Sun S, Wang R, et al. Salvianolic acid A improve mitochondrial respiration and cardiac function via inhibiting apoptosis pathway through CRYAB in diabetic cardiomyopathy. Biomed Pharmacother. 2023;160:114382.
Zhou R, Gao J, Xiang C, et al. Salvianolic acid A attenuated myocardial infarction-induced apoptosis and inflammation by activating Trx. Naunyn Schmiedebergs Arch Pharmacol. 2020;393(6):991–1002.
Zhang Z, Qi D, Wang X, et al. Protective effect of salvianolic acid A on ischaemia-reperfusion acute kidney injury in rats through protecting against peritubular capillary endothelium damages. Phytother Res. 2018;32(1):103–14.
Li S, Qian Q, Ying N, et al. Activation of the AMPK-SIRT1 pathway contributes to protective effects of salvianolic acid A against lipotoxicity in hepatocytes and NAFLD in mice. Front Pharmacol. 2020;11:560905.
Wang K, Yang Q, Ma Q, et al. Protective effects of salvianolic acid A against dextran sodium sulfate-induced acute colitis in rats. Nutrients. 2018;10(6):791.
Xu H, Li Y, Che X, et al. Metabolism of salvianolic acid A and antioxidant activities of its methylated metabolites. Drug Metab Dispos. 2014;42(2):274–81.
Lv H, Jiang F, Guan D, et al. Metabolomics and its application in the development of discovering biomarkers for osteoporosis research. Int J Mol Sci. 2016;17(12):2018.
Xu J, Shen J, Cheng Y, et al. Simultaneous detection of seven phenolic acids in Danshen injection using HPLC with ultraviolet detector. J Zhejiang Univ Sci B. 2008;9(9):728–33.
Hou Y, Peng J, Chao R. Pharmacokinetic study of salvianolic acid A in rat after intravenous administration of Danshen injection. Biomed Chromatogr. 2007;21(6):598–601.
Pei L, Bao Y, Wang H, et al. A sensitive method for determination of salvianolic acid A in rat plasma using liquid chromatography/tandem mass spectrometry. Biomed Chromatogr. 2008;22(7):786–94.
Guo J, Yu H, Xing S, et al. Addressing big data challenges in mass spectrometry-based metabolomics. Chem Commun. 2022;58(72):9979–90.
Li Z, Ryenchindorj L, Liu B, et al. Chemical profiles and metabolite study of raw and processed Cistanche deserticola in rats by UPLC-Q-TOF-MSE. Chin Med. 2021;16(1):95.
Xu X, Li W, Li T, et al. Direct infusion-three-dimensional-mass spectrometry enables rapid chemome comparison among herbal medicines. Anal Chem. 2020;92(11):7646–56.
Zhou F, Teng L, Liu Y, et al. Elaboration of the comprehensive metabolic profile of salvianolic acid A in vivo and in vitro using UFLC-Q/TOF-MS. J Agric Food Chem. 2019;67(44):12199–207.
Shen Y, Wang X, Xu L, et al. Characterization of metabolites in rat plasma after intravenous administration of salvianolic acid A by liquid chromatography/time-of-flight mass spectrometry and liquid chromatography/ion trap mass spectrometry. Rapid Commun Mass Spectrom. 2009;23(12):1810–6.
Song Q, Li J, Huo H, et al. Retention time and optimal collision energy advance structural annotation relied on LC-MS/MS: an application in metabolite identification of an antidementia agent namely echinacoside. Anal Chem. 2019;91(23):15040–8.
Cao Y, Chai C, Chang A, et al. Optimal collision energy is an eligible molecular descriptor to boost structural annotation: an application for chlorogenic acid derivatives-focused chemical profiling. J Chromatogr A. 2020;1609:460515.
Guan P, Liu W, Cao Y, et al. Full collision energy ramp-MS2 spectrum in structural analysis relying on MS/MS. Anal Chem. 2021;93(46):15381–9.
Liu W, Li W, Zhang P, et al. Quality structural annotation for the metabolites of chlorogenic acid in rat. Food Chem. 2022;379:132134.
Tang H, Cao Y, Liu L, et al. High-level structural analysis of proanthocyanidins using full collision energy ramp-MS2 spectrum. J Pharm Biomed Anal. 2022;211:114634.
Li W, Liao Y, Hong Q, et al. Widely quasi-quantitative analysis of both metabolites and tryptic peptides in animal-originated medicinal materials: Bufonis venenum as a case. J Pharm Biomed Anal. 2023;223:115143.
Song Y, Song Q, Liu W, et al. High-confidence structural identification of metabolites relying on tandem mass spectrometry through isomeric identification: a tutorial. TrAC Trend Anal Chem. 2023;160:116982.
Lai C, Tan T, Zeng S, et al. An enzymatic protocol for absolute quantification of analogues: application to specific protopanoxadiol-type ginsenosides. Green Chem. 2015;17(4):2580–6.
Zhang L, Pan S, Lin G, et al. Intestinal absorbability of three Radix Puerariae isoflavones including daidzein, daidzin and puerarin. Chin Med. 2011;6:41.
Li T, Zhang K, Niu X, et al. MS/MS fingerprint comparison between adjacent generations enables substructure identification: flavonoid glycosides as cases. J Pharm Biomed Anal. 2023;234:115559.
Cao Y, Li W, Chen W, et al. Squared energy-resolved mass spectrometry advances quantitative bile acid-submetabolome characterization. Anal Chem. 2022;94(44):15395–404.
Li W, Chen W, Niu X, et al. Characterization of metabolic correlations of ursodeoxycholic acid with other bile acid species through in vitro sequential metabolism and isomer-focused identification. Molecules. 2023;28(12):4801.
Li W, Zhang K, Chen W, et al. Positional isomers of caffeoylquinic acid generate identical fragment ion after collision-induced dissociation. Rapid Commun Mass Spectrom. 2023;37(21):e9636.
Shi J, Zhang K, Niu X, et al. Performance comparison of energy-resolved mass spectrometry between quadrupole cell and linear ion trap chamber. Rapid Commun Mass Spectrom. 2023;37(7):e9475.
Li H, Zhang K, Niu X, et al. Full collision energy ramp-MS2 spectral features of natural esters: salvianolic acid A as a case. Rapid Commun Mass Spectrom. 2022;36(16): e9344.
Guironnet A, Wiest L, Vulliet E. Advantages of MS/MS/MS (MRM3) vs. classic MRM quantification for complex environmental matrices: analysis of beta-lactams in WWTP sludge. Anal Chim Acta. 2022;1205:339773.
Zhu J, Yi X, Zhang J, et al. Chemical profiling and antioxidant evaluation of Yangxinshi tablet by HPLC-ESI-Q-TOF-MS/MS combined with DPPH assay. J Chromatogr B Analyt Technol Biomed Life Sci. 2017;1060:262–71.
Liu A, Guo H, Ye M, et al. Detection, characterization and identification of phenolic acids in Danshen using high-performance liquid chromatography with diode array detection and electrospray ionization mass spectrometry. J Chromatogr A. 2007;1161(1–2):170–82.
Xu J, Zeng S, Chen X, et al. Isolation and identification of degradation products of salvianolic acid A by NMR and LC-MS. Fitoterapia. 2011;82(2):260–6.
Zhang K, Li H, Shi J, et al. Strategy strengthens structural identification through hyphenating full collision energy ramp-MS2 and full exciting energy ramp-MS3 spectra: an application for metabolites identification of rosmarinic acid. Anal Chim Acta. 2024;1296:342346.
Gu J, Feng L, Zhang M, et al. New metabolite profiles of danshensu in rats by ultraperformance liquid chromatography/quadrupole-time-of-flight mass spectrometry. J Chromatogr B. 2014;955–956:20–5.
Mei X, Wang Y, Wang Z, et al. Identification of metabolites of danshensu in vivo in rats. Zhongguo Zhong Yao Za Zhi. 2018;43(19):3933–9.
Funding
This study was financially supported by the National Natural Science Foundation of China (No. 81973444) and the Open Research Project Programme of the State Key Laboratory of Quality Research in Chinese Medicine (University of Macau) (No. SKLQRCM-OP21011).
Author information
Authors and Affiliations
Contributions
Han Li: methodology, investigation, writing—original draft. Ke Zhang: methodology, investigation, writing—original draft. Wei Chen: formal analysis, writing—original draft. Yuxuan Zhou: formal analysis, writing—original draft. Jun Li supervision, writing—review & editing. Yunfang Zhao (corresponding author): writing—review & editing. Yuelin Song (corresponding author): conceptualization, supervision, funding acquisition, writing—review & editing. All authors contributed to data interpretation and preparation of the manuscript for publication and they approved the final version. Han Li and Ke Zhang contributed equally to this article.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The entire animal study protocols were approved by the Animal Ethics Committee of Beijing University of Chinese Medicine (No. BUCM-2023060701–0007).
Consent for publication
Not applicable.
Competing interests
This work has no competing interests that might be relevant to the contents of this manuscript.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
Chromatographic and spectrometric information of in vivo metabolites (F1–F30) for salvianolic acid F (SAF). Table S2. Chromatographic and spectrometric information of in vivo metabolites (D1–D15) for sodium danshensu (sodium DSS). Table S3. FEER-MS3/MS2 spectral information of m/z 211 for M25, D14, and D15. Table S4. FEER-MS3/MS2 spectral information of m/z 327 for M23, M26, F28, and F29. Fig. S1. High-resolution MS/MS spectrum of SAA (A) and the proposed fragmentation pathways (B-1 and B-2). Fig. S2. RIIOEE-trend curves of M25, D14, and D15. Fig. S3. RIIOEE-trend curves of M23, M26, F28, and F29.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, H., Zhang, K., Chen, W. et al. Metabolite identification of salvianolic acid A in rat using post collision-induced dissociation energy-resolved mass spectrometry. Chin Med 19, 64 (2024). https://doi.org/10.1186/s13020-024-00931-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13020-024-00931-z