Next Article in Journal
Discrimination of Adzuki Bean (Vigna angularis) Geographical Origin by Targeted and Non-Targeted Metabolite Profiling with Gas Chromatography Time-of-Flight Mass Spectrometry
Previous Article in Journal
Color Mutations Alter the Biochemical Composition in the San Marzano Tomato Fruit
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 10
Issue 3
10.3390/metabo10030111
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compound Characterization and Metabolic Profile Elucidation after In Vitro Gastrointestinal and Hepatic Biotransformation of an Herniaria hirsuta Extract Using Unbiased Dynamic Metabolomic Data Analysis

by
Laura Peeters
1,*,
Anastasia Van der Auwera
1,
Charlie Beirnaert
2,
Sebastiaan Bijttebier
1,
Kris Laukens
2,
Luc Pieters
1,
Nina Hermans
1 and
Kenn Foubert
1
1
Natural Products & Food Research and Analysis (NatuRA), Department of Pharmaceutical Sciences, University of Antwerp, Universiteitsplein 1, 2610 Antwerp, Belgium
2
Adrem Data Lab, Department of Mathematics—Computer Sciences, University of Antwerp, Middelheimlaan 1, 2020 Antwerp, Belgium
*
Author to whom correspondence should be addressed.
Metabolites 2020, 10(3), 111; https://doi.org/10.3390/metabo10030111
Submission received: 5 February 2020 / Revised: 11 March 2020 / Accepted: 12 March 2020 / Published: 16 March 2020
Browse Figures
Versions Notes

Abstract

:
Herniaria hirsuta L. (Caryophyllaceae) is used for treatment of urinary stones and as a diuretic. Little is known about the active compounds and the mechanism of action. The phytochemical composition of H. hirsuta was comprehensively characterized using UHPLC-UV-HRMS (Ultrahigh-Performance Liquid Chromatography-Ultraviolet-High Resolution Mass Spectrometry) data. An in vitro gastrointestinal model was used to simulate biotransformation, which allowed the monitoring of the relative abundances of individual compounds over time. To analyze the longitudinal multiclass LC–MS data, XCMS, a platform that enables online metabolomics data processing and interpretation, and EDGE, a statistical method for time series data, were used to extract significant differential profiles from the raw data. An interactive Shiny app in R was used to rate the quality of the resulting features. These ratings were used to train a random forest model. The most abundant aglycone after gastrointestinal biotransformation was subjected to hepatic biotransformation using human S9 fractions. A diversity of compounds was detected, mainly saponins and flavonoids. Besides the known saponins, 15 new saponins were tentatively identified as glycosides of medicagenic acid, acetylated medicagenic acid and zanhic acid. It is suggested that metabolites of phytochemicals present in H. hirsuta, most likely saponins, are responsible for the pharmaceutical effects. It was observed that the relative abundance of saponin aglycones increased, indicating loss of sugar moieties during colonic biotransformation, with medicagenic acid as the most abundant aglycone. Hepatic biotransformation of this aglycone resulted in different metabolites formed by phase I and II reactions.

Graphical Abstract

1. Introduction

Urinary stone disease is considered as an important healthcare problem that affects 10–15% of the population in the developed world, but the incidence can be as high as 20–25% in the Middle East, with a peak at ages 20 to 40 years. In addition, the disease is characterized by its high recurrence rate, about 50% in 10 years and 75% in 20 years [1,2,3,4,5,6].
Over the years, many remedies and surgical treatments have been described, varying from dietary recommendations to interventional procedures. Endoscopic management, both ureteroscopic and percutaneous, offers an efficient and efficacious way to treat stones [6]. Nowadays, extracorporeal shock wave lithotripsy (ESWL) is widely used for stone fragmentation without invasive instrumentation of the body. However, the aforementioned therapies show some significant side effects such as residual stone fragments that can act as a nidus for new stone formation [5,7]. Existing oral treatments with evidence that supports their long-term efficacy to prevent kidney stones only include a handful of drugs. Thiazide diuretics, potassium citrate and allopurinol have been used as therapy for more than 30 years. However, their adverse effects prevent long term consumption [2,8].
Phytotherapy might be useful as an alternative treatment. Some herbal remedies have been used for centuries and have been shown to be effective against urolithiasis, although the mechanism of action is often not well established through systematic pharmacological and clinical studies. Many plants have been and are still used to treat urinary stones, such as Phyllanthus niruri, Zea mays, Agropyron repens and Herniaria hirsuta [5,9].
An aqueous extract of the aerial parts of H. hirsuta (hairy rupturewort) is an herbal medicine widely used against urolithiasis and which contains diuretic properties. The Arabic name, herras lehjer, is translated as stonebreaker and refers to its use as a traditional herbal remedy for urinary stones. The beneficial effects of the extract have been demonstrated in several studies [9,10,11,12,13,14]. Even though a lot of research has been done to prove the activity of H. hirsuta against urolithiasis, little is known about the active compounds and the exact mechanism of action. Previous phytochemical research on Herniaria species revealed the presence of saponins, flavonoids and coumarins [15,16,17,18]. Four herniariasaponins have been reported in H. hirsuta: herniariasaponins E, F and G, three monodesmosidic derivatives of medicagenic acid and herniariasaponin H, a bidesmosidic derivative of medicagenic acid [17,19]. Literature suggests that the antilithiatic potential of H. hirsuta is attributed to saponins or metabolites thereof [17,20,21].
It is well known that many natural products are so-called prodrugs, which become active after biotransformation [22]. Nevertheless, this aspect is usually overlooked when searching for new therapeutic agents using classical approaches. After oral administration, an herbal extract is inevitably brought into contact with gastrointestinal enzymes and intestinal microflora, which might lead to biotransformation of the compounds. Saponins are high-molecular-weight glycosides, containing a triterpenoid or steroidal sapogenin aglycone covalently linked to one or two sugar chains via an ether or ester glycosidic bond, referred to as monodesmosidic or bidesmosidic saponins, respectively. Before absorption in the gastrointestinal tract, saponins are presumably hydrolyzed with loss of sugar chains [23]. The aglycones or metabolites thereof may be absorbed and may further be biotransformed in the liver to the ultimate active components [24]. In this study, the biotransformation of a well-characterized H. hirsuta extract was investigated in an in vitro gastrointestinal model followed by hepatic biotransformation. The longitudinal multiclass data were subjected to different data analysis workflows to screen and tentatively identify metabolites formed after in vitro biotransformation.

2. Results and Discussion

2.1. Identification of Compounds

A dereplication workflow performed on the extract of H. hirsuta resulted in a molecular network containing different clusters, composed of different nodes connected by edges which define the degree of similarity between the MS/MS spectra. Compounds that are structurally related are linked to each other. Mass differences between the nodes were studied and provided information about hitherto unidentified compounds. Major clusters were analyzed to tentatively identify unknown compounds. As an example, the cluster containing herniariasaponin H is shown in Figure 1. Herniariasaponin H (Figure 2) is the most abundant saponin previously observed in H. hirsuta and served as starting point for identification of other compounds in this cluster.
Structures were assigned to peaks only when both the mass/charge (m/z) ratios and the molecular formulae of the precursor ion and product ions were in agreement. Additional information about the proposed structures was often provided by photodiode array (PDA) spectra and retention times. However, the information was not always sufficient for full elucidation of the structural composition at an acceptable confidence level. Additional information was obtained from published data.
All amMS and PDA data used for peak identification are listed in Table 1. The table also specifies the literature consulted to obtain additional information in order to confirm compound identity. In accordance with previous research, mainly flavonoids and saponins were identified [15,16,17,18,19,20,21]. However, many of the identified compounds have not been reported previously in H. hirsuta.
Flavonoids were mainly observed as glycoconjugates, i.e., a flavonoid aglycone attached to one or more sugar moieties. Several flavonoid compounds were identified using analytical standards. Others were tentatively identified with ddMS2. Collision-induced fragmentation resulted in product ions caused by elimination of sugar moieties, indicating glycosidic O-linkages [24]. Among others, rutin, narcissin and quercetin-3-O-(2″-O-α-L-rhamnopyranosyl)-β-D-glucuronopyranoside were identified, in agreement with previous studies [17,20]. However, isorhamnetin-3-robinobioside and isorhamnetin-3-[3-feruloylrhamnosyl-(1→6)-galactoside], previously reported to be present in H. fontanesii, were not detected in H. hirsuta, which is probably due to variances in phytochemical profile between different plant species [18].
A rich diversity of saponins was detected, differing in the presence of various aglycones, a different composition of sugar chains and varied linkages of sugar moieties. Fragmentation led to the distinct presentation of product ions caused by elimination of sugar moieties, indicating glycosidic O-linkages. The product ions in the MS2 spectra provide additional structural information about the sugar residues and the aglycones of the fragmented saponins. Losses of m/z 132, 146, 162 and 176 were observed, corresponding to pentose, deoxyhexose, hexose and uronic acid moieties, respectively. Three different aglycone moieties were observed with fragment ions at m/z 543.33234, 541.31787 and 501.32074, corresponding to acetylated medicagenic acid, zanhic acid and medicagenic acid, respectively [25,26].
Several other compounds, such as free phenolic and hydroxycinnamic acids, were detected. Amongst others, caffeic, p-coumaric and chlorogenic acids were identified using analytical standards. The coumarins herniarin and umbelliferone, previously reported to be present in Herniaria species, were not detected. This variation in compounds is probably due to differences in origin or regional climate [17,20]. More than one chromatographic peak was often present for the same precursor ion, indicating the presence of structural isomers with different linkage positions and/or different sugar moieties.
Figure 3 illustrates the identification of a hitherto unknown saponin. As shown in Figure 3A, the unknown saponin has a retention time of 15.33 min. Figure 3B shows a molecular ion at m/z 1703.71838 [M-H]−. Since saponins contain the elements C, H and O, 19 molecular formulas are possible within an error range of 5 ppm.
The sugar residues and the aglycone were identified from the MS2 spectra as shown in Figure 3C and Table 2. The product ion of m/z 1381.62830 [M-H-176-146]− results from the loss of uronic acid (176 u) and deoxyhexose (146 u).
The ion at m/z 879.30701 originates from a sugar chain containing two deoxyhexosyl (146 u), two pentosyl (132 u) and two hexosyl (162 u) units. The ion at m/z 823.41089 [aglycone-H+146+176]− represents the aglycone attached to a deoxyhexose unit (146 u) and a uronic acid moiety (176 u). The ion at m/z 501.32275 is characteristic for medicagenic acid as aglycone, resulting from the loss of sugar residues. The ion at m/z 439.31934 [aglycone-H-18-44]− results from an additional loss of water (18 u) and CO2 (44 u). This suggests a bidesmosidic saponin with R1 comprising uronic acid and deoxyhexose and R2 comprising two pentosyl, two deoxyhexosyl and two hexosyl moieties, with a hexosyl moiety in a terminal position. Combining all the information with the retention time and literature data, the saponin with m/z 1703.71734 supports a molecular formula of C76H120O42 and was tentatively identified as medicagenic acid attached to three deoxyhexosyl, two hexosyl, two pentosyl and a uronic acid moiety [15,16,17,18,19,20,21]. However, MSn does not provide enough information for absolute structural characterization. Nevertheless, the use of a hybrid orbitrap mass analyzer enabled the tentative identification of 15 saponins that have not been reported before in H. hirsuta (indicated with e in Table 2). Full structural characterization is not the main goal in this project, as the complete herbal extract is subjected to biotransformation studies.

2.2. Gastrointestinal Biotransformation

Herbal extracts comprise a mixture of compounds, covering a wide range of bioactive constituents, aside from active compounds or prodrugs. As the main compounds are classified as flavonoids and saponins, extensive biotransformation after oral intake can be expected. Therefore, gastrointestinal biotransformation of the lyophilized extract was simulated in vitro to monitor the levels of the identified compounds using amMS and UV data. However, gastrointestinal enzymes and fecal microflora present in the samples cause a lot of matrix interference, increasing the complexity of the data. The multiclass samples are measured as a function of time, adding a longitudinal aspect to the complex data and impeding the interpretation of the data. The concentration of precursor compounds and metabolites can increase, decrease or show any combination of these patterns during biotransformation.
Figure 4 shows the chromatogram of the H. hirsuta extract before and after biotransformation. Before biotransformation (t0), peaks are mostly attributed to the tentatively identified compounds. After the colonic phase of gastrointestinal biotransformation, the chromatogram contains peaks of compounds, metabolites and matrix interferences, originating from enzymes and fecal microflora, which are absent at t0.
Processing the chromatograms of every time point for every compound allows monitoring of the intensity over time. However, this manual approach is time-consuming and will only provide information about the abundance of previously identified compounds. Manual screening for metabolites is comparable to looking for a needle in a haystack because of the complexity of the chromatograms and matrix interferences.
Generally, for most tentatively identified compounds, a decrease in relative abundance is observed mainly at the colonic phase, caused by bacterial degradation, suggesting formation of metabolites. However, some compounds, mainly saponins, show an increase in relative abundance over time. It is suggested that this increase is induced by conversion of one compound into another, caused by loss of one or more sugar moieties. For example, compound 46 (medicagenic acid + three deoxyhexosyl + two pentosyl + one glucuronyl moiety) shows an increase in relative abundance and reaches a maximum after 24 h colon phase, which is not observed in the negative control samples (NCHEX) (Figure 5). In addition, the signal was absent in the blank (MB), as expected. At the same time, a decrease of compound 45 (medicagenic acid + one hexosyl + four deoxyhexosyl + two pentosyl + one glucuronyl moiety) is observed. The increase of compound 46 can be related to the decrease of compound 45, caused by the loss of a hexosyl and a deoxyhexosyl moiety. The decrease is observed in both samples and negative control during the stomach phase and the small intestinal phase. At these time points, there is no difference in samples and negative control as the negative controls only lack fecal bacteria during the colonic phase. Since identification based on MS-data is not absolute, only tentative conclusions can be drawn from these observations. However, the data support the theory that gastrointestinal biotransformation leads to a loss of sugar moieties in saponins.
An automated data analysis workflow, previously described by Peeters et al. and Beirnaert et al., was used to screen for metabolites [23,27]. The workflow resulted in 11,072 features in negative ion mode, of which 4846 were annotated as isotopes. In positive mode, 5242 features were detected after exclusion of isotopes.
Firstly, the metabolic pathway of herniariasaponin H (m/z 1541.6675, rt: 17.3 min) was studied in detail to verify if the workflow was suitable for herbal extracts. The time profile of herniariasaponin H shows a decrease in intensity over time, especially during the colonic phase, which is not observed in the negative control samples (NCHEX) (Figure 6). In addition, the signal is absent in the method blank (MB). The feature is rated as notable based on the results of tinderesting, the interactive Shiny app developed in R, with a score of 1 and a ranking of 242 out of 6226 features. During biotransformation, stepwise elimination of sugar moieties is observed resulting in formation of the aglycone (Figure 7). All intermediate metabolites show a tinderesting score of at least 0.98, meaning that they are all rated as interesting. Finally, the aglycone medicagenic acid (m/z 501.3213; rt: 21.68) is formed, showing an increase in signal over time in the samples, with a tinderesting score of 0.998 (Figure 6).
Features with a tinderesting score higher than 0.8 were selected for in-depth exploration with the following inclusion criteria: a reasonable intensity with good-quality MS/MS spectra and a retention time between 2 and 25 min.
The main challenge in the analysis of a large amount of analytical data remains in the identification and annotation of the features of interest. However, the data analysis workflow ranks features according to one single score taking the longitudinal aspect and the three sample classes into account, which facilitates selection of features of interest. On the other hand, determination of the molecular formula of compounds with a molecular mass above 1000 Dalton becomes more challenging. Therefore, this step in the data analysis process is still performed manually.
As expected, the features included in the in-depth exploration could be attributed to saponins, flavonoids and their metabolites. Saponins (compounds 35–41, 43–45, 47–48, 51, 58, 59 and 61) show the same biotransformation pathway as the studied example, herniariasaponin H. In brief, the intensity of the signal decreases over time due to stepwise elimination of sugar moieties, leading to formation of the aglycone. The observed biotransformation of flavonoids is similar. In nature, flavonoids mainly occur as glycosides. Gastrointestinal biotransformation causes deglycosylation, leading to formation of the aglycone. For example, quercetin glycosides are abundantly and ubiquitously present flavonoids in nature. The aglycone is formed after biotransformation of several precursors, such as isoquercetin and rutin. Isoquercetin (m/z 463.0871; rt: 11.67 min) has a tinderesting score of 0.86 and shows a decrease in intensity over time (Figure 8). The decrease is already observed in the stomach phase and continues in the small intestinal phase and colonic phases until the signal is completely absent. In the negative control samples, the signal stays stable during the colonic phase, because there are no bacteria present. The signal is absent in the method blank. For quercetin (m/z 301.0345; rt: 16.54 min), an increase in intensity is observed over time, which is not observed in the negative control or method blank (Figure 8). The feature is rated as interesting with a tinderesting score of 1. This suggests deglycosylation of isoquercetin leading to the formation of quercetin.
A feature with m/z 181.0497 and rt: 4.9 min was included in the in-depth exploration as it showed a tinderesting score of 0.992. The feature shows an increase in intensity over time, mainly during the colonic phase (Figure 8). The metabolite was tentatively identified as 3,4-dihydroxyphenylpropionic acid (A) or 3-methoxy-4-hydroxyphenylacetic acid (B). Both compounds were previously described as metabolites of quercetin [28,29]. Based on the obtained spectra, absolute structural characterization is not possible. Nonetheless, the product ion could be identified as a metabolite of quercetin. This shows that despite the increase in intensity of quercetin, the compound is further biotransformed (Figure 9).

2.3. Liver Biotransformation

Figure 10 shows an overview of the proposed metabolic pathway of medicagenic acid after S9 biotransformation. For metabolites 1–3 (M1–3), three different chromatographically separated peaks were observed (rt 17.39; 17.60 and 18.45 min), suggesting the presence of three different structural isomers. This was confirmed by the prediction provided by Meteor. The software proposed hydroxylation at C-11, C-24, or C-29. Double hydroxylation (M4), oxidation of a hydroxy group (M5) and a combination of hydroxylation and oxidation (M9–10) were also predicted by Meteor. The nontarget screening approach confirmed the formation of the predicted metabolites, and one additional metabolite was detected. Metabolite 8 (M8) shows a reduction of the carboxylic acid moiety. However, limited by the poor MS/MS fragmentation of the aglycones and lack of NMR data, identification remains tentative. Further elucidation of the molecular structures was not possible. The only phase II products observed were glucuronidated metabolites of medicagenic acid (M6–7). Sulfation reactions were not observed, although they were predicted by Meteor and Biotransformer. To confirm that the in vitro model was capable to form sulfated conjugates, biotransformation of testosterone was examined as positive control. Testosterone was analyzed after hepatic biotransformation and proved the capability of the in vitro model to form both sulfated and glucuronidated conjugates.
The phase II hydrophilic metabolites may be excreted via the kidneys and are often considered detoxifying as they usually possess lower activity compared to their parent aglycones. However, a number of glucuronide conjugates are known to be active. The 6-O-glucuronide of morphine, a widely used opioid analgesic, is the most well-known example of a glucuronidated metabolite possessing pharmacological activity greater than the parent aglycone [30]. Based on the observation of in vitro formation of glucuronidated hepatic metabolites, it is a possible hypothesis that these metabolites may be responsible for the urolithiatic activity [31]. A possible mode of action for prevention and dissolution of urinary stones could be the emulsifying properties of the glucuronidated metabolites, facilitating dissolution of crystals.

3. Materials and Methods

3.1. Chemicals

Ultrapure water with a resistivity of 18.2 MΩ·cm at 25 °C was generated with a Millipore™-purification system. Ultrahigh-Performance Liquid Chromatography (UHPLC)-grade MeOH, acetonitrile and formic acid were purchased from Biosolve (Dieuze, France), dichloromethane by Merck (Darmstadt, Germany). The following analytical standards were provided by Sigma–Aldrich (St. Louis, MO, USA): apigenin, benzoic acid, caffeic acid, +/− catechin, chlorogenic acid, cinnamic acid, coumarin, emodin, epicatechin, ferulic acid, isorhamnetin, naringenin, p-coumaric acid, protocatechuic acid, quercetin, quercitrin, rutin, salicylic acid, sinapic acid, β-sitosterol, stigmasterol, syringic acid, tannic acid, taxifolin, theophylline and vanillic acid. Luteolin and procyanidin B2 were provided by Santa Cruz Biotechnology (Santa Cruz, CA, USA). Gallic acid and p-hydroxybenzoic acid were provided by Carl Roth (Karlsruhe, Germany). Human liver S9 fraction and NADPH RS (Regenerating System) were purchased from XenoTech. Medicagenic acid was provided by Phytolab (Vestenbergsgreuth, Germany). All other chemicals and biochemicals were purchased from Sigma-Aldrich (St. Louis, MO, USA).

3.2. Preparation of Standard Solutions

Standard stock solutions of the analytical standards were prepared at a concentration of 1 mg mL−1 in UHPLC-grade MeOH for each analyte separately and stored in the dark at −80 °C. Dilutions of these solutions were prepared in 60:40 (v:v) MeOH:ammonium formate buffer (40 mM, aqueous). Standard stock and working solutions were stored at −80 °C in the dark.

3.3. Sample Preparation

Aerial parts from H. hirsuta were collected in d’Oujda, Morocco. A voucher specimen of the plant is kept at the Muséum National d’Histoire Naturelle—Institut Scientifique-Université Mohamed V Agdal, Morocco—(Reference number: 5902). The air-dried plant material of H. hirsuta was ground prior to extraction with a PF 10 basic Microfine grinder drive (IKA-Werke GmbH & Co. KG, Staufen, Germany) using a sieve mesh size 0.5 mm.
To cover the full range of constituents, a comprehensive extraction protocol was applied, previously described by Bijttebier et al. with minor adaptations [32,33]. Combination of extraction with H2O:EtOAc and CHCl3:MeOH:H2O afforded both polar and apolar compounds. Approximately 18 g of sample was subsequently wetted with 54 mL of water and 225 mL of EtOAc. The mixture was stirred for 2 h followed by 1 h of ultrasound-assisted extraction. The solvent was removed, and the residue was washed three times with 20 mL EtOAc. The extract was evaporated to dryness. CHCl3:MeOH:H2O extraction was performed by adding 18 g of sample to 180 mL CHCl3:MeOH:H2O (4:4:2). The mixture was stirred for 5 min, followed by 5 min of ultrasound-assisted extraction and stirred for 15 more min. The solvent was removed, and the residue was washed 3 times with 20 mL CHCl3:MeOH:H2O (4:4:2). The dried H2O:EtOAc extract was redissolved with the extract from the CHCl3:MeOH:H2O extraction method, leading to an apolar and a polar phase. Both phases were separated and dried.

3.4. Gastrointestinal Biotransformation

Many natural products are pro-drugs that are biotransformed and activated after oral administration. A previously developed and validated in vitro gastrointestinal biotransformation model was used to mimic human biotransformation processes [34]. This model allows the study of biotransformation processes in the stomach, small intestine and colon using a culture of pooled human feces, avoiding extensive in vivo studies.
The digestive juices and fecal suspension were formulated corresponding to the human conditions, previously described by Breynaert et al. and Peeters et al. [23,34]. Small adaptations to the method were made [35]. Briefly, the pepsin solution was prepared by dissolving 16% (w/v) of pepsin powder in 0.1 M HCl (622,000 FIP-U 100 mL−1). To obtain a pancreatin-bile mixture, 0.4% (w/v) pancreatin and 2.5% (w/v) bile were dissolved in 0.1 M NaHCO3 (32,000 FIP-U lipase, 143,600 FIP-U amylase, 16,400 FIP-U protease and 58.4 mmol bile L−1). A suspension of 10% (w/v) feces was prepared by homogenizing 3 human fecal samples with a sterile phosphate-buffer solution (0.1 M, pH 7.0). The phosphate-buffer solution consisted of NaH2PO4.2H2O (1.03% w/v), sodium thioglycolate broth (3.45% v/v) and glycerol (17%). Prior to use, the fecal suspension was cultivated. To 10% (v/v) pooled fecal suspension in phosphate buffer 90% (v/v) phosphate buffer was added in an anaerobic glove-box (5% H2, 5% CO2 and 90% N2) (Jacomex Glove Box T3, TCPS, Belgium). The bacteria were cultivated in an anaerobic environment for 1 h at 37 °C by continuously stirring.
Three groups were included in the experiment: (1) samples containing the Herniaria extract, which were treated with digestive enzymes and fecal microflora (HEX); (2) negative control samples, also containing the extract with addition of digestive enzymes but not of fecal bacteria (NCHEX); and (3) method blanks, comprising an equal volume of solvent instead of extract and undergoing treatment with digestive enzymes and fecal bacteria (MB). For the samples, an amount of 300 mg lyophilized polar phase of the comprehensive extract of H. hirsuta was accurately weighed in triplicate.
Gastrointestinal biotransformation was simulated in vitro. An aliquot of each sample was taken at several time points: before biotransformation (t0), after the gastric stage (S), after the small intestinal phase (SI) and at several time points during the colonic phase (after 2, 4, 6, 10, 14, 18, 22, 24, 32, 40, 48 and 72 h). Every sample was diluted with MeOH (1:2) and centrifuged for 5 min at 3500 rpm. The supernatant was collected, and samples were diluted 10 times with MeOH:H2O (60:40) before analysis.

3.5. Liver Biotransformation

Liver biotransformation mimicking phase I and II reactions was performed using medicagenic acid, the most abundant aglycone after gastrointestinal biotransformation. Hepatic biotransformation was simulated in vitro by using pooled S9 fractions, previously described by Van den Eede et al. [36]. Small adaptations to the method were made. Briefly, a mixture of 65 mM TRIS buffer (pH adjusted to 7.4 at 37 °C), human liver S9 fraction (1 mg mL−1 final protein concentration) and medicagenic acid (100 µM final concentration) was prepared in a total volume of 0.5 mL and preincubated in a shaking waterbath at 37 °C. The reaction was initiated by addition of an NADPH regenerating system (1 mM NADP, final concentration). Cofactors were added to expose the samples to phase II glucuronidation and sulfation: UDPGA (2 mM final concentration), GSH (2 mM final concentration) and PAPS (0.1 mM final concentration). The reaction was inhibited after 1 h by adding 0.5 mL of acetonitrile and storing the tubes on ice. Thereafter, the tubes were centrifuged for 5 min at 10,000 rpm (4 °C). The supernatant was collected and analyzed by LC–MS-MS. Blank samples were prepared as described above, replacing substrate by solvent. Negative control samples were prepared by immediately adding acetonitrile to quench the biotransformation reactions. A positive control was included by incubating testosterone (100 µM final concentration) [36].

3.6. Instrumental Analysis

For characterization of the H. hirsuta extract, an aliquot of 5 µL was injected with a CTC PAL™ autosampler (CTC Analytics) on a Waters Acquity UHPLC BEH SHIELD RP18 column (3.0 mm × 150 mm, 1.7 µm; Waters, Milford, MA, USA) and eluted with an Accela™ quaternary solvent manager. The temperature of the column was kept at 40 °C using a “Hot Pocket” column oven (Thermo Fisher Scientific, Waltham, MA, USA). The mobile phase solvents consisted of water + 0.1% formic acid (A) and acetonitrile + 0.1% formic acid (B). Compounds were eluted using the following gradient (min/B%): 0.00/0, 9.91/26, 18.51/65, 18.76/100, 20.76/100, 20.88/0, 23.00/0. For detection, an amMS (Q Exactive™; Thermo Fisher Scientific) was used with heated electrospray ionization (HESI). Spray voltage was set at −2.5 kV, sheath gas and auxiliary gas at 47 and 15 (adimensional), respectively, and capillary temperature at 350 °C. The PDA detector was set to scan from 190 to 800 nm. During the first analysis, full scan data were acquired in positive and negative ion mode over an m/z range of 120–1800 and resolving power was set at 70,000 at full width at half maximum (FWHM). Consecutively, a second analysis was performed to selectively fragment the generated ions by data-dependent fragmentation (ddMS2) in positive and negative ionization mode. Selective ion fragmentation was obtained by using a hybrid quadrupole-orbitrap MS analyzer (Q Exactive™, Thermo Fisher Scientific). The raw data were processed using XCalibur™ 3.0 software (Thermo Fisher Scientific).
For the qualitative UHPLC-UV-QTOF analyses of the biotransformation samples, an aliquot of 5 µL was injected on a Waters Acquity UHPLC BEH SHIELD RP18 column (3.0 mm × 150 mm, 1.7 µm; Waters). The temperature of the column was kept at 40 °C. The mobile phase solvents consisted of water + 0.1% formic acid (A) and acetonitrile + 0.1% formic acid (B), and the gradient was set as follows (min/B%): 0/2, 1/2, 14/26, 24/65, 26/100, 29/100, 31/2, 41/2. The flow rate was set at 0.4 mL min−1. For detection, accurate mass measurements were done using a Xevo G2-XS QTof spectrometer (Waters, Milford, MA, USA) coupled with an ACQUITY LC system equipped with MassLynx 4.1 software. During the first analysis, full scan data were recorded in ESI (+) and ESI (−) mode from m/z 50 to 2000, and the analyzer was used in sensitivity mode (approximate resolution: 22,000 FWHM). The spray voltage was set at either +1.5 kV and −1.0 kV; cone gas flow and desolvation gas flow at 50.0 L h−1 and 1000.0 L h−1 respectively; source temperature and desolvation temperature at 120 °C and 500 °C respectively. Data were also recorded using MSE in positive and negative ionization modes, providing separate MS and MSE data. A ramp collision energy from 20 V to 30 V was applied to obtain additional structural information. Leucine encephalin was used as lock mass. To monitor analytical drift and assess precision, quality control (QC) samples were injected after every time point. QC samples were prepared using a dilution of the standard solution mix (39 ng mL−1).

3.7. Data Analysis

Dereplication was used as an untargeted workflow for rapid identification of the major compounds present in the extract, regardless of their chemical class. Global Natural Products Social Molecular Networking (GNPS, a web-based mass spectrometry open-access knowledge base) was used to facilitate identification [37]. Raw files were uploaded into the METABOLOMICS-SNETS-V2 workflow from GNPS (Global Natural Products Social Molecular Networking) in order to cluster the data based on their MS/MS fragmentation pattern. The precursor ion mass tolerance and the fragment ion mass tolerance were set at 2.0 Da and 0.5 Da respectively, and the cosine score was set at 0.7. The generated molecular network was visualized using Cytoscape 3.7.1 software.
To encounter the dynamic nature of the biotransformation process, a data analysis workflow for in vitro gastrointestinal biotransformation was implemented, previously developed and validated using hederacoside C as a model compound [23,27]. Briefly, data were converted to the open-source mzData format to allow further processing [38]. XCMS was used to convert the raw data into features via peak-picking, using following parameters: ppm = 10, peakwidth = c(5, 25, snthresh = 10, noise = 1000, mzdiff = 0.01, prefilter = c(3, 5000), integrate = 1. XCMS was followed by grouping, and EDGE was used for the extraction of significant differential profiles [39,40]. Tinderesting, an interactive Shiny app developed in R, was used to rate the quality of the resulting features. These ratings were used to train a random forest model for predicting experts response. The machine learning model provided a single score for each feature, referred to as tinderesting score, which allowed ranking of all features based on the difference over time between the three groups. The maximal score of 1 corresponds to the model labeling this feature as interesting, and the minimal score of 0 defines an uninteresting feature [23,27].
Several peaks in the output correspond to the same type of sample molecule, reflecting its isotope pattern. Isotopes were removed, using following parameters: IsoMass = 1.003355, IsoMassMaxDiff = 0.05, RTdiff = 10, corrThreshold = 0.8, ppmIsoDiffThreshold = 5.
The metabolic pathway of herniariasaponin H was studied since it is the most abundant saponin present in the extract and it is structurally related to hederacoside C, which was used for optimization of the workflow [23]. Both compounds contain a sapogenin linked to sugar moieties via the hydroxy groups at C-3 and C-28 of the aglycone. Herniariasaponin H comprises medicagenic acid as aglycone, which only differs from hederagenin by an additional hydroxy group at C-2 and a C-24 hydroxycarbonyl group instead of a hydroxymethyl group. However, the challenge in revealing the metabolic pathway of herniariasaponin H was that these experiments were not conducted with a single saponin, but with an herbal extract, comprising a complex mixture of compounds of high structural variety.
A suspect screening method, based on in silico metabolite prediction, was combined with a nontarget screening workflow to enhance the identification of products formed by in vitro liver biotransformation [41].
A list of potential biotransformation products was generated using Meteor Nexus 2.1 (Lhasa Limited, Leeds, UK). For phase I biotransformation, all redox and non-redox biotransformations were selected. For phase II biotransformation, O-glucuronidation, O-sulfonylation, acetylation and conjugation with amino acids were selected. The maximum depth was set at 3, and the maximum number of biotransformation products at 1000.
Biotransformer was used as an additional software tool to predict phase II metabolites. The tool uses both a knowledge-based approach and a machine-learning-based approach to predict biotransformation [42]. The SMILES string of medicagenic acid was uploaded, and “Phase I Transformation” and “Phase II Transformation” were selected separately. The generated csv file contained InChIKey, synonyms, major isotope mass, molecular formula, type of biotransformation reaction and precursor ID.
A nonarget screening workflow was conducted with MZmine 2.29 software (pluskal, mzmine 2) after converting the raw data files to the open-source mzXML format. Firstly, m/z features were detected using the centroid algorithm followed by a chromatogram building step. Resulting chromatograms were deconvoluted using the noise amplitude algorithm. Further data reduction occurred by deisotoping, keeping the lowest m/z value as the representative isotope. Meanwhile, chromatographic peaks were also filtered according to peak width: only peaks with a width between 0.05 and 1.00 min were retained. Next, peaks were aligned across samples using the random sample consensus (RANSAC) alignment algorithm correcting for nonlinear shifts in retention times. Finally, using the same rt and m/z range gap filler algorithm, any missing peaks were reiteratively extracted. The obtained m/z features were imported in R and processed as described by Vervliet et al. [43]. In brief, fold changes between samples and negative controls were calculated. A volcano plot was constructed to plot the p-value from a student t-test as a function of the calculated fold change for every feature. Features with a p-value lower than 0.05 and a log 10-fold change higher than 10 were selected for in depth investigation [41].

4. Conclusions

Since there is no prophylactic therapy for urinary stone disease, there is an urgent need for new drugs. An integrated strategy, based on natural prodrugs and their metabolites, is developed to characterize new lead compounds derived from H. hirsuta. The combination of a comprehensive extraction method and biotransformation via in vitro gastrointestinal and hepatic models, followed by metabolomics profiling, is an innovative concept to identify and further develop new lead compounds for drugs.
Comprehensive extraction led to the tentative identification of 15 saponins that have not been reported before in H. hirsuta. The automated data analysis workflow used for unbiased screening for metabolites revealed the gastrointestinal biotransformation pathway of several compounds. A decrease in relative abundance over time was observed for the majority of all identified compounds, especially during the colon phase. Additionally, the relative abundance of saponin aglycones increased, illustrating the biotransformation of possible saponin prodrugs to their respective aglycones, which was also observed for flavonoids. Further hepatic biotransformation of the most abundant saponin aglycone revealed further biotransformation of medicagenic acid to different phase I and phase II metabolites, which can be used as a possible target for in vitro activity studies against urinary stones.

Author Contributions

Conceptualization, L.P. (Laura Peeters), A.V.d.A., C.B., S.B., K.L., L.P. (Luc Pieters), N.H. and K.F.; methodology, L.P. (Laura Peeters), A.V.d.A.; software, C.B., L.P. (Laura Peeters); validation, L.P. (Laura Peeters), A.V.d.A. and C.B.; formal analysis, L.P. (Laura Peeters), A.V.d.A.; investigation, L.P. (Laura Peeters).; resources, K.L., L.P. (Luc Pieters) and N.H.; data curation, L.P. (Laura Peeters).; writing—original draft preparation, L.P. (Laura Peeters); writing—review and editing, L.P. (Laura Peeters), L.P. (Luc Pieters) and K.F.; visualization, L.P. (Laura Peeters); supervision, K.L., L.P. (Luc Pieters) and N.H.; project administration, L.P. (Luc Pieters) and N.H.; funding acquisition, L.P. (Luc Pieters) and N.H. All authors have read and agreed to the published version of the manuscript.

Funding

Research was funded by Special Fund for Research of the University of Antwerp (Concerted Action), project ID no. 30732 and by the Research Foundation (FWO)—Flanders Project no. G032415N.

Acknowledgments

Y. Cherrah (University of Rabbat, Morocco) is kindly acknowledged for providing the original plant material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eisner, B.H.; Goldfarb, D.S. Pharmacologic Treatment of Kidney Stone Disease. Urol. Clin. N. Am. 2013, 40, 21–30. [Google Scholar] [CrossRef] [PubMed]
  2. Pachaly, M.A.; Baena, C.P. Therapy of Nephrolithiasis: Where Is the Evidence from Clinical Trials? J. Bras. Nefrol. 2016, 38, 99–106. [Google Scholar] [CrossRef] [PubMed]
  3. Rule, A.D.; Lieske, J.C. The Roks Nomogram for Predicting a Second Symptomatic Stone Episode. J. Am. Soc. Nephrol. 2014, 25, 2878–2886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Trinchieri, A. Epidemiology of Urolithiasis: An Update. Clin. Cases Miner. Bone Metab. 2008, 5, 101–106. [Google Scholar]
  5. Gurocak, S.; Kupeli, B. Consumption of Historical and Current Phytotherapeutic Agents for Urolithiasis: A Critical Review. J. Urol. 2006, 176, 450–455. [Google Scholar] [CrossRef]
  6. Moe, O.W. Kidney Stones: Pathophysiology and Medical Management. Lancet 2006, 367, 333–344. [Google Scholar] [CrossRef]
  7. Kirkali, Z.; Rasooly, R. Urinary Stone Disease: Progress, Status, and Needs. Urology 2015, 86, 651–653. [Google Scholar] [CrossRef] [Green Version]
  8. Erickson, S.B.; Vrtiska, T.J. Effect of Cystone(R) on Urinary Composition and Stone Formation over a One Year Period. Phytomedicine 2011, 18, 863–867. [Google Scholar] [CrossRef] [Green Version]
  9. Atmani, F.; Slimani, Y. Effect of Aqueous Extract from Herniaria hirsuta L. On Experimentally Nephrolithiasic Rats. J. Ethnopharmacol. 2004, 95, 87–93. [Google Scholar] [CrossRef]
  10. Atmani, F.; Slimani, Y. Prophylaxis of Calcium Oxalate Stones by Herniaria hirsuta on Experimentally Induced Nephrolithiasis in Rats. BJU Int. 2003, 92, 137–140. [Google Scholar] [CrossRef] [Green Version]
  11. Aggarwal, A. Urolithiasis: Phytotherapy as an Adjunct Therapy. Indian J. Exp. Biol. 2014, 53, 103–111. [Google Scholar]
  12. Fouada, A.; Yamina, S. In Vitro and in Vivo Antilithiasic Effect of Saponin Rich Fraction Isolated from Herniaria hirsuta. J. Bras. Nefrol. 2006, 28, 199–203. [Google Scholar]
  13. Atmani, F. Medical Management of Urolithiasis, What Opportunity for Phytotherapy. Front. Biosci. 2003, 8, 507–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Atmani, F.; Farell, G. Extract from Herniaria hirsuta Coats Calcium Oxalate Monohydrate Crystals and Blocks Their Adhesion to Renal Epithelial Cells. J. Urol. 2004, 172, 1510–1514. [Google Scholar] [CrossRef]
  15. Charrouf, Z.; Nait-Mbark, A. Herniaria Saponin B, a Novel Triterpenoid Saponin from Herniaria Fontanesii. Adv. Exp. Med. Biol. 1996, 405, 241–245. [Google Scholar]
  16. Mbark, A.N.; Charrouf, Z. Herniaria Saponin a, a Novel Saponin Fromherniaria Fontanesii. Nat. Prod. Lett. 2006, 6, 233–240. [Google Scholar] [CrossRef]
  17. van Dooren, I.; Foubert, K. Saponins and Flavonoids from an Infusion of Herniaria hirsuta. Planta Med. 2016, 82, 1576–1583. [Google Scholar] [CrossRef] [Green Version]
  18. MBark, A.N.; Guillaume, D. Triterpenoid Saponins from Herniaria Fontanesii. Phytochemistry 1996, 43, 1075–1077. [Google Scholar] [CrossRef]
  19. Mbark, A.N.; Charouf, Z. Monodesmosidic Saponins from Herniaria hirsuta. Pharmazie 2000, 55, 690–692. [Google Scholar]
  20. Wagner, H.; Bladt, S. Drogenanalyse: Dünnschichtchromatographische Analyse Von Arzneidrogen; Springer: Berlin/Heidelberg, Germany, 1983; pp. 225–230. [Google Scholar]
  21. Schröder, H.; Manfred, S.-Z. A triterpene saponin from Herniaria glabra. Phytochemistry 1993, 34, 1609–1613. [Google Scholar] [CrossRef]
  22. Rautio, J.; Kumpulainen, H. Prodrugs: Design and Clinical Applications. Nat. Rev. Drug Discov. 2008, 7, 255–270. [Google Scholar] [CrossRef] [PubMed]
  23. Peeters, L.; Beirnaert, C. Revelation of the Metabolic Pathway of Hederacoside C Using an Innovative Data Analysis Strategy for Dynamic Multiclass Biotransformation Experiments. J. Chromatogr. A 2019, 1595, 240–247. [Google Scholar] [CrossRef] [PubMed]
  24. Osbourn, A.E. Preformed Antimicrobial Compounds and Plant Defense against Fungal Attack. Plant Cell 1996, 8, 1821–1831. [Google Scholar] [CrossRef] [PubMed]
  25. Cuyckens, F.; Claeys, M. Mass Spectrometry in the Structural Analysis of Flavonoids. J. Mass Spectrom. 2004, 39, 1–15. [Google Scholar] [CrossRef]
  26. Mikolajczyk-Bator, K.; Blaszczyk, A. Identification of Saponins from Sugar Beet (Beta Vulgaris) by Low and High-Resolution Hplc-Ms/Ms. J. Chromatogr. B 2016, 1029, 36–47. [Google Scholar] [CrossRef]
  27. Beirnaert, C.; Peeters, L. Using Expert Driven Machine Learning to Enhance Dynamic Metabolomics Data Analysis. Metabolites 2019, 9, 54. [Google Scholar] [CrossRef] [Green Version]
  28. Aura, A.-M.; O’leary, K. Quercetin Derivatives Are Deconjugated and Converted to Hydroxyphenylacetic Acids but Not Methylated by Human Fecal Flora in Vitro. J. Agric. Food Chem. 2002, 50, 1725–1730. [Google Scholar] [CrossRef]
  29. Justesen, U.; Arrigoni, E. Electrospray Ionisation Mass Spectrometric Study of Degradation Products of Quercetin, Quercetin-3-Glucoside and Quercetin-3-Rhamnoglucoside, Produced by in Vitro Fermentation with Human Faecal Flora. Rapid Commun. Mass Spectrom. 2001, 15, 477–483. [Google Scholar] [CrossRef]
  30. Osborne, R.; Thompson, P. The Analgesic Activity of Morphine-6-Glucuronide. Br. J. Clin. Pharmacol. 1992, 34, 130–138. [Google Scholar] [CrossRef]
  31. Ritter, J.K. Roles of Glucuronidation and Udp-Glucuronosyltransferases in Xenobiotic Bioactivation Reactions. Chem. Biol. Interact. 2000, 129, 171–193. [Google Scholar] [CrossRef]
  32. ijttebier, S.; Van der Auwera, A. Bridging the Gap between Comprehensive Extraction Protocols in Plant Metabolomics Studies and Method Validation. Anal. Chim. Acta 2016, 935, 136–150. [Google Scholar] [CrossRef] [PubMed]
  33. Theodoridis, G.; Gika, H. Lc-Ms Based Global Metabolite Profiling of Grapes: Solvent Extraction Protocol Optimisation. Metabolomics 2012, 8, 175–185. [Google Scholar] [CrossRef]
  34. Breynaert, A.; Bosscher, D. Development and Validation of an in Vitro Experimental Gastrointestinal Dialysis Model with Colon Phase to Study the Availability and Colonic Metabolisation of Polyphenolic Compounds. Planta Med. 2015, 81, 1075–1083. [Google Scholar] [CrossRef] [PubMed]
  35. Mortelé, O.; Iturrospe, E. Optimization of an in Vitro Gut Microbiome Biotransformation Platform with Chlorogenic Acid as Model Compound: From Fecal Sample to Biotransformation Product Identification. J. Pharm. Biomed. Anal. 2019, 112768. [Google Scholar] [CrossRef]
  36. Van den Eede, N.; Maho, W. First Insights in the Metabolism of Phosphate Flame Retardants and Plasticizers Using Human Liver Fractions. Toxicol. Lett. 2013, 223, 9–15. [Google Scholar] [CrossRef]
  37. Wang, M.; Carver, J.J. Sharing and Community Curation of Mass Spectrometry Data with Global Natural Products Social Molecular Networking. Nat. Biotechnol. 2016, 34, 828–837. [Google Scholar] [CrossRef] [Green Version]
  38. Stanstrup, J.; Gerlich, M. Metabolite Profiling and Beyond: Approaches for the Rapid Processing and Annotation of Human Blood Serum Mass Spectrometry Data. Anal. Bioanal. Chem. 2013, 405, 5037–5048. [Google Scholar] [CrossRef]
  39. Storey, J.D.; Xiao, W. Significance Analysis of Time Course Microarray Experiments. Proc. Natl. Acad. Sci. USA 2005, 102, 12837–12842. [Google Scholar] [CrossRef] [Green Version]
  40. Leek, J.T.; Monsen, E. Edge: Extraction and Analysis of Differential Gene Expression. Bioinformatics 2006, 22, 507–508. [Google Scholar] [CrossRef] [Green Version]
  41. Vervliet, P.; Mortele, O. Suspect and Non-Target Screening Workflows to Investigate the in Vitro and in Vivo Metabolism of the Synthetic Cannabinoid 5cl-Thj-018. Drug Test. Anal. 2019, 11, 479–491. [Google Scholar] [CrossRef]
  42. Djoumbou-Feunang, Y.; Fiamoncini, J. Biotransformer: A Comprehensive Computational Tool for Small Molecule Metabolism Prediction and Metabolite Identification. J. Cheminform. 2019, 11, 1–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Vervliet, P. Invitro-in Vitro Hlm Assay Feature Prioritization (Version V1.0.1). Zenodo 2018. [Google Scholar] [CrossRef]
Metabolites 10 00111 g001 550
Figure 1. Molecular network in negative ion mode with the cluster containing herniariasaponin H shown in detail.
Figure 1. Molecular network in negative ion mode with the cluster containing herniariasaponin H shown in detail.
Metabolites 10 00111 g001
Metabolites 10 00111 g002 550
Figure 2. Chemical structure of herniariasaponin H.
Figure 2. Chemical structure of herniariasaponin H.
Metabolites 10 00111 g002
Metabolites 10 00111 g003 550
Figure 3. LC–MS data for an unknown saponin. (A) Chromatogram obtained for the extract of Herniaria hirsuta. (B) MS scan of peak at 15.33 min in the heated electrospray ionization (HESI) negative mode. (C) MS2 fragmentation spectrum of the ion at m/z 1703.71734.
Figure 3. LC–MS data for an unknown saponin. (A) Chromatogram obtained for the extract of Herniaria hirsuta. (B) MS scan of peak at 15.33 min in the heated electrospray ionization (HESI) negative mode. (C) MS2 fragmentation spectrum of the ion at m/z 1703.71734.
Metabolites 10 00111 g003
Metabolites 10 00111 g004 550
Figure 4. Total ion chromatogram in negative ion mode of H. hirsuta before and after biotransformation.
Figure 4. Total ion chromatogram in negative ion mode of H. hirsuta before and after biotransformation.
Metabolites 10 00111 g004
Metabolites 10 00111 g005 550
Figure 5. Time profile of the biotransformation of compounds 46 (left) and 45 (right) in samples (HEX), negative controls (NCHEX) and blank (MB).
Figure 5. Time profile of the biotransformation of compounds 46 (left) and 45 (right) in samples (HEX), negative controls (NCHEX) and blank (MB).
Metabolites 10 00111 g005
Metabolites 10 00111 g006 550
Figure 6. Time profile of gastrointestinal biotransformation of herniariasaponin H (m/z 1541.6675 [M-H]−) and medicagenic acid (m/z 501.3213 [M-H]−).
Figure 6. Time profile of gastrointestinal biotransformation of herniariasaponin H (m/z 1541.6675 [M-H]−) and medicagenic acid (m/z 501.3213 [M-H]−).
Metabolites 10 00111 g006
Metabolites 10 00111 g007 550
Figure 7. Metabolic pathway of herniariasaponin H by in vitro gastrointestinal biotransformation.
Figure 7. Metabolic pathway of herniariasaponin H by in vitro gastrointestinal biotransformation.
Metabolites 10 00111 g007
Metabolites 10 00111 g008 550
Figure 8. Time profiles of biotransformation of isoquercetin (m/z 463.0871 [M-H]−), quercetin (m/z 301.0345 [M-H]−) and a metabolite (m/z 181.0497 [M-H]−) tentatively identified as 3,4-dihydroxyphenylpropionic acid or 3-methoxy-4-hydroxyphenylacetic acid.
Figure 8. Time profiles of biotransformation of isoquercetin (m/z 463.0871 [M-H]−), quercetin (m/z 301.0345 [M-H]−) and a metabolite (m/z 181.0497 [M-H]−) tentatively identified as 3,4-dihydroxyphenylpropionic acid or 3-methoxy-4-hydroxyphenylacetic acid.
Metabolites 10 00111 g008
Metabolites 10 00111 g009 550
Figure 9. Biotransformation pathway of isoquercetin showing quercetin as intermediate metabolite which is further transformed into 3,4-dihydroxyphenylpropionic acid (A) or 3-methoxy-4-hydroxyphenylacetic acid (B).
Figure 9. Biotransformation pathway of isoquercetin showing quercetin as intermediate metabolite which is further transformed into 3,4-dihydroxyphenylpropionic acid (A) or 3-methoxy-4-hydroxyphenylacetic acid (B).
Metabolites 10 00111 g009
Metabolites 10 00111 g010 550
Figure 10. B Metabolic pathway of medicagenic acid after in vitro hepatic biotransformation.
Figure 10. B Metabolic pathway of medicagenic acid after in vitro hepatic biotransformation.
Metabolites 10 00111 g010
Table
Table 1. Chromatographic and Spectrometric Data of the Tentatively Identified Compounds in H. hirsuta Detected with a Generic LC-PDA-amMS Method for Moderately Polar Phytochemicals.
Table 1. Chromatographic and Spectrometric Data of the Tentatively Identified Compounds in H. hirsuta Detected with a Generic LC-PDA-amMS Method for Moderately Polar Phytochemicals.
Com-Pound NumberCompound aMolecular FormulaRT (min) HESI Full MS bHESI ddMS2Neutral LossesMax Abs. (nm)Literature
1quinic acid cC7H12O61.54neg191.05611-
pos--
2gallic acid cC7H6O54.06neg169.01412-
pos--
3protocatechuic acid [3,4-dihydroxybenzoic acid] cC7H6O46.46neg153.01933109.02950109.02950: C6H6O2259; 294
pos--
4caffeyl hexoseC15H18O97.68neg341.08817179.03439; 135.04439179.03439: C9H8O4; 135.04439: C9H8O4 - CO2
pos365.08380 [M+Na]+; 360.12839 [M+NH4]+181.04962; 163.03906; 145.02861; 135.04443181.04962: C9H8O4; 163.03906: C9H8O4 - H2O; 145.02861: C9H6O4 - 2H2O; 135.04443: C8H6O2
5asperulosidic acidC18H24O128.01neg431.11975; 499.10757 [M-H+NaFA]−; 863.24756 [2M-H]−329.09064; 161.04463329.09064: C14H18O9; 161.04463: C6H10O5
pos433.13428271.08188; 187.06091; 145.04990; 127.03919271.08188: C18H24O12 - C6H10O5; 187.06091: C8H10O5; 145.04990: C6H8O4; 127.03919: C6H6O3
6dihydrocaffeic acid cC9H10O48.07neg181.05041-
pos--
71-O-feruloylquinic acid dC17H20O98.76neg367.10349; 435.09108 [M-H+NaFA]−-
pos369.11795-
8chlorogenic acid [3.4-dihydroxycinnamoylquinic acid; 5-Caffeoylquinic acid] cC16H18O98.94neg353.08685; 707.18355 [2M-H]−191.05545191.05611: C7H12O6326
pos355.10160; 377.08360 [M+Na]+; 731.17920 [2M+Na]+163.03897163.03897: C9H6O3
9isorhamnetin-3-O-rutinoside-7-O-glucoside dC34H42O219.27neg785.21631; 853.20349 [M-H+NaFA]−623.16479; 476.09396; 315.05124; 300.02750; 271.02466 623.16479: C34H42O21 - C6H10O5; 315.05124: C16H12O7; 300.02750: C15H9O7; 271.02466: C14H7O6
pos787.22931641.17114; 479.11838641.17114: C34H42O21 - C6H10O4; 479.11838: C34H42O21 - C6H10O4 - C6H10O5
10chlorogenic acid cC16H18O99.36neg353.08781191.05611191.05611: C7H12O6
pos355.10236163.03897
11Catechin cC15H14O611.21neg289.07176- 279Mbark et al.
pos--
12epicatechin cC15H14O611.25neg289.07176- 280
pos--
13procyanidin B2 cC30H26O129.68neg577.13515-
pos579.14970-
14caffeic acid cC9H8O49.7neg179.03498135.04515135.04515: C9H8O4 - CO2
pos181.04954163.03915; 145.02861; 135.04443163.03897: C9H8O4 - H2O
151-O-feruloylquinic acid dC17H20O9 10.51neg367.10358; 435.09106 [M-H+NaFA]−193.05003; 173.04489; 134.03622193.05003: C10H10O4; 173.04485: C7H10O5
pos391.10001 [M+Na]+145.02855; 149.05952; 117.03381; 177.05475; 194.05750145.02855: C9H4O2; 149.05952: C9H8O2; 117.03381: C8H4O; 177.05475: C10H10O4 - H2O
16quercetin 3-rhamnosyl-(1->2)-rhamnosyl-(1->6)-glucoside or quercetin 3-rhamninosideC33H40O2011.05neg755.20471; 823.19204 [M-H+NaFA]−; 1511.41418 [2M-H]−300.02771; 271.02499; 255.02994300.02771: C15H9O7; 271.02499: C14H8O6; 255.02994: C14H8O5
pos757.21765; 611.16085465.1033; 303.05008465.1033: C33H40O20 - 2C6H10O4; 303.05008: C33H40O20 - 2C6H10O4 - C6H10O5
17matairesinosideC26H32O1111.3neg565.19336 [M-H+FA]−339.12378; 329.13812; 327.12405; 324.09967; 312.09885; 309.07953; 297.07669339.12378: C20H20O5; 329.13812: C19H22O5; 327.12405: C19H20O5; 324.09967: C19H17O5; 312.09885: C18H17O5; 309.07953: C18H14O5; 297.07669: C17H14O5355
pos543.18420 [M+Na]+; 538.22786 [M+NH4]+503.19199; 341.13672; 323.12738; 311.12756; 175.07526; 137.05983503.19199: C26H30O10; 341.13672: C20H20O5; 323.12738: C20H18O4; 311.12756: C19H18O4; 175.07526: C11H10O2; 137.05983: C8H8O2
18isorhamnetin-O-glycoside (one hexosyl and two deoxyhexosyl units)C34H42O2011.74neg769.22040 [M-H]; 1539.44527 [2M-H]−; 632.32239357.06104; 314.04349; 271.02512; 299.01981; 285.04065632.32239: C34H42O20 - C6H10O4; 357.06104: C18H14O8; 314.04349: C16H11O7; 299.01981: C15H8O7; 285.04065: C15H10O6; 271.02512: C14H8O6355
pos771.23315; 625.17635; 479.11849 625.17743; 479.11874; 317.06540625.17743: C34H42O20 - C6H10O4; 479.11874: C34H42O20 - 2C6H10O4; 317.06540: C34H42O20 - 2C6H10O4 - C6H10O5
19quercetin 3-O-β-(2″-O-α-L-rhamnopyranosyl)-glucopyranosideC27H28O1711.96neg623.12537; 1247.25657 [2M-H]−301.03543301.03543: C27H28O17 - C6H10O4 - C6H8O6 255; 355van Dooren et al. (2016)
pos625.13943; 479.08180479.08240; 303.04999479.08180: C27H28O17 - C6H10O4; 303.04999: C27H28O17 - C6H10O4 - C6H8O6
20rutin cC27H30O1612.08neg609.14667; 677.13349 [M-H+NaFA]−300.02762; 271.02481; 255.02990; 243.02990300.02755: C15H9O7; 271.02481: C14H8O6; 255.02990: C14H8O5; 243.02990: C13H8O5256; 355van Dooren et al. (2016)
pos611.15973465.10267; 303.04999465.10275: C21H22O12; 303.04993: C15H12O7
21cynaroside cC21H20O1112.2neg447.09364- Wagner et al. (2013)
pos--
22hesperidin cC28H34O1512.24neg609.18153-
pos611.19705-
23ferulic acid cC10H10O412.34neg193.05063178.02643178.02643: C9H7O4
pos195.06519177.0547; 145.02855; 135.04424; 117.03378177.0547: C10H10O4 - H2O; 149.05984: C9H12O2; 145.02855: C9H4O2; 135.04424: C8H6O2; 117.03378: C8H4O
24herniariasaponin 1C54H56O3212.59neg1215.2698-
pos1239.26343 [M+Na]+-
25isoquercetin cC21H20O1212.61neg463.0882300.02725; 271.02539; 255.02905; 243.02859; 151.00230300.02725: C15H9O7; 271.02539: C14H8O6; 255.02905: C14H8O5; 243.02859: C13H8O5; 151.00230: C7H4O4
pos465.10275303.05008; 287.05533; 257.04483303.05008: C15H10O7; 287.05533: C15H10O6; 257.04483: C14H8O5
26p-coumaric acid cC9H8O312.62neg163.04007119.05024119.05024: C8H8O
pos165.05462147.04414147.04414: C9H6O2
272-(3.4-Dihydroxyphenyl)-5.7-dihydroxy-4-oxo-4H-chromen-3-yl-6-O-(2-O-acetyl-α-L-arabinopyranosyl)-β-D-glucopyranosideC28H30O1712.67neg637.14125; 1275.28809 [2M-H]−463.09000; 300.02728463.09000: C21H20O12; 300.02728: C15H9O7
pos639.15543; 493.09749465.10071; 303.04972493.09749: C28H30O17 – C6H10O4; 465.10071: C21H20O12; 303.04972: C15H10O7
28miquelianin [quercetin-3-O-glucuronide]cC21H18O1312.89neg477.06746301.03538; 271.02481 255.02990; 178.99860; 151.00368301.03538: C15H10O7; 271.02481: C14H8O6; 255.02990: C14H8O5; 178.99860: C8H4O5; 151.00368: C7H4O4
pos479.08202303.04993; 257.04445303.04993: C15H10O7; 257.04445: C14H8O5
29narcissin [isorhametin-3-O-rutinoside]C28H32O1612.89neg623.16211; 477.06749; 1247.33026 [2M-H]−315.05072477.06749: C28H32O16 - C6H10O4; 315.05072: C28H32O16 - C6H10O4 - C6H10O5 254; 355van Dooren et al. (2016) Wagner et al. (2013)
pos625.14606479.11804; 317.06586479.11804: C28H32O16 - C6H10O4; 317.06586: C28H32O16 - C6H10O4 - C6H10O5
30avicularin [quercetin-3-O-arabinoside] cC20H18O1113.16neg433.0782301.03586 256; 355
pos435.09219-
31apigetrin cC21H24O913.27neg431.09837268.03827268.03827: C15H9O5
pos433.11292-
32astragalin [kaempferol 3-O-glucoside] cC21H20O1113.34neg447.09381284.03314; 255.03024; 227.03455284.03314: C15H9O6; 255.03024: C14H8O5; 227.03455: C13H8O4265; 348
pos449.10814-
33quercitrin cC21H20O1113.35neg447.09328- 265; 346
pos449.10784-
34salicylic acid cC7H6O314.07neg137.0244293.0345993.03459: C6H6O236; 302
pos--
35herniariasaponin FC59H94O2814.17neg1249.58482; 647.29342501.32303; 439.32358647.29342: C30H46O6 + C6H10O4; 501.32303: C30H46O6; 439.32358: C30H46O6 - H2O - CO2 Mbark et al. (2000)
pos1273.58240 [M+Na]+; 1268.62600 [M+NH4]+; 1122.56885; 645.27747; 457.33133771.25336457.33133: C29H44O4; 771.25336: 2C6H10O5 + 2C6H10O4 + C5H8O4 + Na
36acetylated medicagenic acid + two hexosyl + two deoxyhexosyl + one pentosyl moieties eC61H96O2914.53neg1291.59644 1249.59241; 1231.57312; 821.43506; 543.33441; 501.32294; 483.31107; 439.32251821.43506: C61H96O29 - 2C6H10O5 - C6H10O4; 543.33441: C32H48O7; 501.32294: C30H46O6; 483.31107: C30H46O6 - H2O; 439.32251: C30H46O6 - H2O - CO2
pos1310.63818 [M+NH4]+ -
37acetylated medicagenic acid + three hexosyl + one deoxyhexosyl moieties eC56H88O2614.61neg1175.549 705.38928; 543.33337; 501.32281; 439.32205705.38928: C32H48O7 + C6H10O5; 543.33337: C32H48O7; 501.32281: C30H46O6; 439.32205: C30H46O6 - H2O - CO2 Schröder et al. (1993)
pos1194.59088 [M+NH4]+-
38medicagenic acid + two hexosyl + two deoxyhexosyl + two pentosyl + 1 glucuronic acid moiety eC70H110O3814.89neg1557.6595; 1411.603171235.56738; 983.57922; 796.23248; 439.320431411.60317: C70H110O38 - C6H10O4; 1235.56738: C70H110O38 - C6H10O4 - C6H8O6; 439.32043: C29H44O3 Mbark et al. (2000)
pos1581.65795 [M+Na]+1259.56189; 741.24219; 595.18250; 445.13116
39zanhic acid + one hexose + four deoxyhexoses + two pentoses + one glucuronic acid moiety eC76H120O4214.92neg1703.71708; 851.35632 [M-2H]2−1381.62830; 823.41089; 517.31879; 455.31766; 839.408451381.62830: C76H120O42 - C6H10O4 - C6H8O6; 517.31879: C30H46O7; 455.31766: C30H46O7 - H2O - CO2; 839.40845: C30H47O7 + C6H10O4 + C6H8O6
pos1727.71399 [M+Na]+; 1722.75891 [M+NH4]+ -
40herniariasaponin EC55H86O2515.03neg1145.5387; 1213.52771 [M-H+NaFA]−; 572.26715; 595.26941 [M-2H-FA]2−543.33105; 483.31180; 439.32227; 601.19806; 439.32208543.33105: C32H48O7; 483.31180: C30H46O6 - H2O; 439.32227: C30H46O6 - H2O - CO2; 601.19806: C5H8O4 + C6H10O4 + 2C6H10O5; 439.32208: C5H8O4 + C6H10O4 + C6H10O5 Mbark et al. (2000)
pos1169.53381 [M+Na]+; 1164.57869 [M+NH4]+545.34430545.34430: C32H48O7
41acetylated medicagenic acid + three hexosyl moieties eC50H78O2215.06neg1029.49191; 1075.49754 [M-H+FA]−543.33105; 483.31180; 439.32227; 791.07910; 585.34424543.33105: C32H48O7; 483.31180: C30H44O5; 439.32227; 791.07910: C29H28O26
pos1048.53173 [M+NH4]+; 1053.48612 [M+Na]+ -
pos1398.65376 [M+NH4]+-
42medicagenic acid + three deoxyhexosyl + one pentosyl + one glucuronic acid moiety eC59H92O2815.12neg1247.569925.48315; 501.31976
pos1271.56687 [M+Na]+1127.49792; 949.47644; 701.34851; 525.32062; 447.14737
43zanhic acid + three deoxyhexosyl + one pentosyl + one glucuronic acid moiety eC59H92O2915.24neg1263.56571163.15430; 1081.94128
pos1282.60583 [M+NH4]+ -
44medicagenic acid + two hexosyl + three deoxyhexosyl + two pentosyl + one glucuronic acid moiety eC76H120O4215.33neg1703.71741381.62830; 1219.22852; 879.30701; 843.52264; 823.41089; 501.32275; 439.31934
pos1727.71399 [M+Na]+; 1722.75891 [M+NH4]*+-
45medicagenic acid + one hexosyl + four deoxyhexosyl + two pentoses + one glucuronic acid moiety eC76H120O41 15.52neg1687.72291365.63477; 998.50537; 823.42078
pos1711.71936 [M+Na]+ -
46medicagenic acid + three deoxyhexosyl + two pentosyl + one glucuronic acid moiety d,eC64H100O3215.12; 15.62neg1379.61271057.52454; 925.48077; 823.41821; 501.32251; 439.320501057.52454: C64H100O32 - C6H8O6 - C6H10O4; 925.48077: C64H100O32 - C6H8O6 - C6H10O4 – C5H8O4; 823.41821: C35H68O21; 501.32251: C30H46O6; 439.32050: C30H46O6 - H2O - CO2
pos1403.60791 [M+Na]+1259.54272; 1081.51782; 847.40808; 579.18964
47medicagenic acid + one hexosyl + two deoxyhexosyl + two pentosyl + one glucuronic acid moiety eC64H100O3315.67neg1395.6069439.32248; 1219.58301439.32248: C30H46O6 - H2O - CO2; 1219.58301: C64H100O33 - C6H8O6 Charrouf et al. (1995)
pos1419.60059 [M+Na]+ -
48herniariasaponin BC65H102O3315.83neg1409.62441087.53479; 823.41272; 501.32300; 439.323431087.53479: C65H102O33 - C6H10O4 - C6H8O6; 823.41272; 501.32300: C30H46O6; 439.32343: C30H46O6 - H2O - CO2 Charrouf et al. (1995)
pos1433.61982 [M+Na]+; 1428.66283 [M+NH4]+; 1411.63638 -
49acetylated medicagenic acid + one deoxyhexosyl + two pentosyl moieties eC47H74O19 16neg941.47546501.32312; 439.32037; 821.43097; 543.33167821.43097: C47H74O19 - C5H8O4; 543.33167: C32H48O7; 501.32312: C30H46O6; 439.32037: C30H46O6 - H2O - CO2
pos960.51593 [M+NH4]+; 965.47021 [M+Na]+463.14227: C17H27O13Na
50quercetin cC15H10O716.4neg301.03538271.02539; 178.99774; 151.00273; 121.02832271.02539: C14H8O6; 178.99774: C8H4O5; 151.00273: C7H4O4; 121.02832: C7H6O2255; 372
pos303.04993285.03964; 257.04440; 229.04947; 201.05476; 165.01811; 153.01846; 137.02357285.03964: C15H8O6; 257.04440: C14H8O5; 229.04947: C13H8O4; 201.05476: C12H8O3; 165.01811: C8H4O4; 153.01846: C7H4O4; 137.02357: C7H4O3
51acetylated medicagenic acid + one hexosyl + two deoxyhexosyl + one pentosyl + one glucuronic acid moietyeC61H94O3016.76neg652.28625 [M-2H]2−; 1305.57764 439.32321; 585.20691439.32321: C29H43O3; 585.20691: C5H8O4 + 2C6H10O4 + C6H10O5
pos1329.57288 [M+Na]+; 1324.61719 [M+NH4]+-
52naringenin cC15H12O516.8neg271.06120253.05069253.05069: C15H10O4
pos273.07575-
53apigenin cC15H10O517.73neg269.04555225.05510; 151.00255; 117.03339225.05510: C14H10O3; 151.00255: C7H4O4
pos271.06010153.01845153.01845: C7H4O4
54isorhamnetin cC16H12O717.87neg315.05103300.02762; 271.02399; 243.02707; 216.04128; 164.01071; 151.00258300.02762: C15H9O7
pos317.06555302.04205; 285.03882; 153.01846153.01846: C7H4O4
55luteolin cC15H10O618.09neg285.04046199.03967; 151.00212; 133.02821 199.03967: C12H8O3; 151.00212: C7H4O4; 133.02821: C8H6O2
pos287.05501153.01787; 135.04414153.01787: C7H4O4; 135.04414: C8H6O2
56kaempferol cC15H10O618.09neg285.04046- 265; 364
pos287.05501258.05115; 231.06519; 165.01807; 153.01816; 121.02866258.05115: C14H9O5; 231.06519: C13H10O4; 165.01807: C8H4O4; 153.01816: C7H4O4; 121.02866: C7H4O2
57miquelianin [quercetin 3-O-glucuronide] cC21H18O1312.89–13.49neg477.0675301.03552; 271.02505; 255.03020;178.99788; 151.00258 257; 357
pos479.0824303.04993; 257.04370
58herniariasaponin GC53H84O24 14.61neg1103.52829; 574.26304 [M-2H+FA]−821.43842; 501.32306; 439.32156821.43842: C36H70O20; 501.32306: C30H46O6; 439.32156: C30H46O6 -H2O - CO2 van Dooren et al. (2016)
pos1127.52344 [M+Na]+; 1122.56841 [M+NH4]+; 811.44629; 503.33636; 457.33087973.50128; 811.45264; 649.39624; 503.33704; 457.33188811.44629: C53H84O24 - C5H8O4 - C6H10O5; 503.33636: C30H46O6; 973.50128: C53H84O24 - C5H8O4; 649.39624: C53H84O24 - C5H8O4 - 2C6H10O5
59medicagenic acid + two hexosyl + two deoxyhexosyl + two pentosyl + one glucuronic acid moiety eC70H110O3814.89–15.38neg1557.6595; 1411.603171235.56738; 983.57922; 796.23248; 439.320431411.60317: C70H110O38 - C6H10O4
pos1581.65795 [M+Na]+1259.56189; 741.24219; 595.18250; 445.13116
60medicagenic acid + one hexosyl + two deoxyhexosyl + two pentosyl + one glucuronic acid moiety eC64H100O3314.92–15.67neg1395.60691219.583
pos - -
61herniariasaponin H dC70H110O3715.12; 15.66; 15.72neg1541.66571446.38293; 1219.58069; 823.41711; 501.325871219.58069: C70H110O37 - C5H8O4 - C6H8O6; 823.41711: C70H110O37 - C6H10O5 - 2C5H8O4 - 2C6H10O4; 501.32587: C30H46O6 van Dooren et al. (2016)
pos1565.66118 [M+Na]+; 1560.70778 [M+NH4]+1243.57349; 905.31189; 741.24255
62chlorogenic acid dC16H18O97.31; 8.88neg353.08685191.05544; 179.03427; 173.04494191.05544: C7H12O6; 179.03427: C9H8O4; 173.04494: C7H12O6 - H2O
pos355.1016; 377.08360 [M+Na]+163.03915; 145.02861; 135.04443163.03915: C9H8O4 - H2O; 145.02861: C9H8O4 - 2H2O; 135.04443: C8H8O2
631-O-feruloylquinic acid dC17H20O99.20–10.62neg367.10349; 435.09108 [M-H+NaFA]−193.05003; 173.04485; 134.03622 Jaiswal et al. (2011)
pos369.11748; 391.09949 [M+Na]+145.02855; 149.05952; 117.03381; 177.05475; 194.05750
a If not specified otherwise: tentative identification based on accurate mass b Deprotonated and protonated molecules in negative and positive mode, respectively, unless stated otherwise c Identification with an analytical standard d Tentatively identified compound or an isomer e Tentative identification of saponins that have never been reported before in H. hirsuta.
Table
Table 2. Most probable sum formula and ppm deviation for the fragment ions of m/z 1703.71734.
Table 2. Most probable sum formula and ppm deviation for the fragment ions of m/z 1703.71734.
m/z [M-H]−Sum Formula∆ ppm
1703.71838C76H120O420.141
1381.62830C64H102O320.113
879.30701C34H56O269.445
823.41089C42H64O16−1.541
501.32275C30H45O61.172
439.31934C29H44O3−5.528

Share and Cite

MDPI and ACS Style

Peeters, L.; Van der Auwera, A.; Beirnaert, C.; Bijttebier, S.; Laukens, K.; Pieters, L.; Hermans, N.; Foubert, K. Compound Characterization and Metabolic Profile Elucidation after In Vitro Gastrointestinal and Hepatic Biotransformation of an Herniaria hirsuta Extract Using Unbiased Dynamic Metabolomic Data Analysis. Metabolites 2020, 10, 111. https://doi.org/10.3390/metabo10030111

AMA Style

Peeters L, Van der Auwera A, Beirnaert C, Bijttebier S, Laukens K, Pieters L, Hermans N, Foubert K. Compound Characterization and Metabolic Profile Elucidation after In Vitro Gastrointestinal and Hepatic Biotransformation of an Herniaria hirsuta Extract Using Unbiased Dynamic Metabolomic Data Analysis. Metabolites. 2020; 10(3):111. https://doi.org/10.3390/metabo10030111

Chicago/Turabian Style

Peeters, Laura, Anastasia Van der Auwera, Charlie Beirnaert, Sebastiaan Bijttebier, Kris Laukens, Luc Pieters, Nina Hermans, and Kenn Foubert. 2020. "Compound Characterization and Metabolic Profile Elucidation after In Vitro Gastrointestinal and Hepatic Biotransformation of an Herniaria hirsuta Extract Using Unbiased Dynamic Metabolomic Data Analysis" Metabolites 10, no. 3: 111. https://doi.org/10.3390/metabo10030111

APA Style

Peeters, L., Van der Auwera, A., Beirnaert, C., Bijttebier, S., Laukens, K., Pieters, L., Hermans, N., & Foubert, K. (2020). Compound Characterization and Metabolic Profile Elucidation after In Vitro Gastrointestinal and Hepatic Biotransformation of an Herniaria hirsuta Extract Using Unbiased Dynamic Metabolomic Data Analysis. Metabolites, 10(3), 111. https://doi.org/10.3390/metabo10030111

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop