Abstract
Plant-based therapeutics, including herbal products, continue to represent a growing facet of the contemporary health care market. Mechanistic descriptions of the pharmacokinetics and pharmacodynamics of constituents composing these products remain nascent, particularly for metabolites produced following herbal product ingestion. Generation and characterization of authentic metabolite standards are essential to improve the quantitative mechanistic understanding of herbal product disposition in both in vitro and in vivo systems. Using the model herbal product, milk thistle, the objective of this work was to biosynthesize multimilligram quantities of glucuronides of select constituents (flavonolignans) to fill multiple knowledge gaps in the understanding of herbal product disposition and action. A partnership between clinical pharmacology and natural products chemistry expertise was leveraged to optimize reaction conditions for efficient glucuronide formation and evaluate alternate enzyme and reagent sources to improve cost effectiveness. Optimized reaction conditions used at least one-fourth the amount of microsomal protein (from bovine liver) and cofactor (UDP glucuronic acid) compared with typical conditions using human-derived subcellular fractions, providing substantial cost savings. Glucuronidation was flavonolignan-dependent. Silybin A, silybin B, isosilybin A, and isosilybin B generated five, four, four, and three monoglucuronides, respectively. Large-scale synthesis (40 mg of starting material) generated three glucuronides of silybin A: silybin A-7-O-β-d-glucuronide (15.7 mg), silybin A-5-O-β-d-glucuronide (1.6 mg), and silybin A-4´´-O-β-d-glucuronide (11.1 mg). This optimized, cost-efficient method lays the foundation for a systematic approach to synthesize and characterize herbal product constituent glucuronides, enabling an improved understanding of mechanisms underlying herbal product disposition and action.
Introduction
Plant-based therapeutics, including herbal products, have existed for thousands of years, predating recorded history. Today, herbal products represent an ever-flourishing facet of the contemporary health care market, which has grown steadily in the United States since passage of the Dietary Supplement Health and Education Act in 1994 (Lindstrom et al., 2014). Despite millennia of usage, mechanistic descriptions of the pharmacokinetics and pharmacodynamics of constituents composing these products remain nascent.
The top-selling herbal product, milk thistle [Silybum marianum (L.) Gaertn. (Asteraceae)], is an extensively studied hepatoprotectant, as well as a perpetrator of herb-drug interactions (Wellington and Jarvis, 2001; Gurley et al., 2006, 2012; Abenavoli et al., 2010; Brantley et al., 2010, 2013, 2014a,b; Loguercio and Festi, 2011; Gurley, 2012; Polyak et al., 2013a,b; Gufford et al., 2014a, 2015). Despite considerable investigation, quantitative data describing mechanisms underlying the disposition of milk thistle constituents remain elusive and are virtually nonexistent when considering herbal constituent metabolites generated in vivo. Failing to account for the disposition and possible bioactivity of these metabolites may be a source of frequently cited in vitro–in vivo disconnects when translating preclinical “hits” to clinical application (Won et al., 2010, 2012). Generation and characterization of authentic metabolite standards of these constituents are essential to improve the quantitative mechanistic understanding of their disposition in both in vitro and in vivo systems.
Synthesis of metabolite standards of herbal products poses challenges beyond those with drug metabolites (Walker et al., 2011, 2014; Di and Obach, 2015). For example, because herbal products are plant-derived, product variation can occur at multiple stages of production, ranging from plant growth and sourcing to harvest, drying, and manufacturing (Won et al., 2010, 2012; Brantley et al., 2014a; Gufford et al., 2014b). These challenges in consistency necessitate careful consideration of key individual constituents prior to synthesis of metabolite standards. In addition, the supply of individual constituents to use as starting materials for the generation of metabolite standards can range from limited to nonexistent, a challenge not frequently encountered when synthesizing a metabolite standard of a conventional small-molecule pharmaceutical. Addressing these additional hurdles requires collaboration between the divergent yet complementary skill sets of clinical pharmacologists and natural products chemists. This partnership can unravel the complicated composition and pharmacology of herbal products to improve the understanding of the disposition and action of these products.
Chemoenzymatic generation of major β-glucuronides of the flavonolignans silybin A and silybin B, the primary constituents of the semipurified milk thistle extract silibinin, has been reported (Kren et al., 2000; Han et al., 2004; Jančová et al., 2011; Křen et al., 2013; Charrier et al., 2014). However, approaches to enhance the efficiency of these procedures remain unexplored. Moreover, the crude milk thistle extract silymarin contains at least seven flavonolignan diastereoisomers, including isosilybin A and isosilybin B, which also undergo extensive glucuronidation. Challenges in isolating sufficient quantities of individual parent flavonolignans from silymarin have likely hindered generation of glucuronide metabolites for these structurally related constituents. Methods to isolate gram quantities of pure flavonolignans, including silybin A, silybin B, isosilybin A, and isosilybin B, have been developed (Graf et al., 2007; Monti et al., 2010). Chemical synthesis of glucuronide conjugates could be pursued but would require synthesis of an extensive library of orthogonally protected precursors of each flavonolignan (Stachulski and Meng, 2013; Zhang et al., 2013). Chemoenzymatic processes using isolated constituents as starting materials, followed by purification and structural characterization using natural products chemistry protocols, can be used to efficiently generate flavonolignan glucuronides on a multimilligram scale.
The cost-effective availability of authentic glucuronide metabolites of herbal product constituents will facilitate evaluation of pharmacologic activity, quantitative assessment of metabolic clearance pathways, direct quantification in biologic matrices to determine systemic exposure, and development of additional assays to describe herbal product constituent disposition. Using milk thistle as the model herbal product, the objective of this work was to generate glucuronides of select flavonolignans to fill multiple knowledge gaps in the understanding of herbal product disposition and action. The aims were to 1) enhance efficiency by optimizing reaction conditions for maximal glucuronide formation and 2) evaluate alternate enzyme and reagent sources to improve cost effectiveness. The results form the foundation for a systematic approach to synthesize and characterize herbal product metabolites (Fig. 1).
Schematic of proposed approach to generate and characterize authentic herbal product glucuronides. PD, pharmacodynamics; PK, pharmacokinetics.
Materials and Methods
Materials and Chemicals.
Human liver microsomes (HLMs; pooled from 50 donors, mixed gender) and bovine liver microsomes (BLMs; pooled from two donors, mixed gender, >12 months old) and bovine liver S9 fraction (pooled from two donors, mixed gender, >12 months old) were purchased from Xenotech, LLC (Lenexa, KS). Bovine serum albumin (BSA), MgCl2, silibinin, and UDP-glucuronic acid (UDPGA) were purchased from Sigma-Aldrich (St. Louis, MO). Alamethicin was purchased from Cayman Chemical Company (Ann Arbor, MI). The individual milk thistle flavonolignans were purified as described in detail previously (Graf et al., 2007) and were >97% pure as determined by ultra-performance liquid chromatography (UPLC) (Napolitano et al., 2013). Methanol [liquid chromatography/mass spectrometry (MS) grade], ethanol, Tris-HCl, Tris base, and formic acid were purchased from Fisher Scientific (Waltham, MA). Strata-X 33-µm polymeric reversed-phase solid phase extraction columns (60 mg of medium each) were purchased from Phenomenex (Torrance, CA). Acetone-d6 was purchased from Cambridge Isotope Laboratories Inc. (Tewksbury, MA).
Initial Optimization of Conditions.
Glucuronidation reaction conditions evaluating a range of microsomal and UDPGA concentrations were examined and compared with typical conditions (Walsky et al., 2012; Walker et al., 2014) to optimize yield and cost. A series of 12 different conditions (Table 1) were examined at six different time points over a 24-hour period for silybin A, silybin B, isosilybin A, and isosilybin B. Conditions examined included enzyme source (BLMs versus S9 fraction), enzyme concentration (1, 0.5, and 0.25 mg/ml for BLMs), and UDPGA concentration relative to flavonolignan concentration (10-fold and 2-fold molar excess). Optimization of alamethicin concentration and source were evaluated as potential cost-saving measures using a series of reactions with varying alamethicin concentration and an in-house alamethicin source (alamethicin F50, which was isolated from a filamentous fungal culture, MSX70741, as detailed previously) (Ayers et al., 2012).
Comparison between typical and screened reaction conditions to optimize cost, efficiency, and yield of milk thistle glucuronides
A master mixture of MgCl2, BSA, and alamethicin was prepared in Tris-HCl buffer (100 mM, pH 7.4) at 3 times the final reaction concentration. This mixture was split into two aliquots, and UDPGA was added at either 10- or 2-fold molar excess. Each of these mixtures was split into four aliquots, and the four flavonolignans were added separately from dimethylsulfoxide (DMSO) stock solutions (19.28 mg/ml). Each flavonolignan-containing mixture was divided further into six microcentrifuge tubes. Three different volumes of BLMs (2 mg/ml) or bovine S9 (4 mg/ml) were added, depending on the desired concentration, with Tris-HCl buffer added to produce a final reaction volume of 0.65 ml. The 48 mixtures were incubated at 37°C, after which 100-μl aliquots were removed and quenched with 300 µl of ice-cold MeOH at each time point (0, 2, 4, 8, 12, and 24 hours), generating 288 samples, which were stored at 4°C until analysis.
Samples were prepared for UPLC-MS analysis by centrifugation (15,000g × 2 minutes), and 60 µl of the supernatant was removed and mixed with 90 µl of H2O in a 96-well plate (final concentration of 30% MeOH and 40 µM total flavonolignans). The resultant mixture (3-µl injection volume) was separated using an Acquity UPLC HSS C18 (1.8 µm, 50 × 2.1 mm) column with an Acquity UPLC system (Waters, Milford, MA) coupled to a Q Exactive Plus high-resolution electrospray ionization mass spectrometer (Thermo Fisher Scientific, San Jose, CA) operated in positive mode.
Large-Scale Reaction of Silybin A and Purification of Glucuronides.
After determining the optimal reaction conditions for silybin A (0.8 mM UDPGA, 0.25 mg/ml BLMs, 37°C × 12 hours), a large-scale reaction was initiated. Tris-HCl (204 ml, 100 mM, pH 7.4), MgCl2∙6H2O (212 mg, 5 mM), UDPGA (108 mg, 0.8 mM), BSA (1.05 g, 0.5% w/v), silybin A (40 mg in 1.7 ml of DMSO, 0.4 mM), alamethicin (8 mg in 0.4 ml of DMSO, 38 μg/ml, 19.5 μM), and BLMs (2.5 ml of 20 mg/ml, 0.25 mg/ml) were combined to yield a total reaction volume of 209 ml (1.0% DMSO final). The reaction was incubated at 37°C in a rotary incubator at 100 rpm to prevent settling of the BLMs during the incubation and was terminated after 12 hours by adding 3 volumes (630 ml) of ice-cold MeOH and stored at 4°C. The resulting mixture was filtered [qualitative grade Whatman filter paper (GE Healthcare Life Sciences, Little Chalfont, Buckinghamshire, UK) in a Büchner funnel], and the supernatant was evaporated to a volume of ∼20 ml. Wash solution (10% MeOH) was added to yield a volume of ∼50 ml (5% MeOH). The reaction concentrate was purified further using six 60-mg Strata-X cartridges Phenomenex, loading 1 ml at a time over eight successive purification cycles. Before the first cycle, each cartridge was conditioned with 4 ml of 100% MeOH and equilibrated with 3 ml of 10% MeOH. For each cycle, the cartridges were pre-equilibrated with 3 ml of 10% MeOH, loaded with 1.0 ml of the reaction concentrate, and washed with 3 ml of 10% MeOH; the combined load and wash fractions were designated pool 1. The sample-loaded cartridges were eluted with 3 ml of 50% MeOH, generating pool 2, then eluted with 4 ml of 100% MeOH, generating pool 3. Each of the pools was dried under vacuum, yielding 4707 mg of pool 1, 51 mg of pool 2, and 24 mg of pool 3.
Pool 2 (a white solid), containing the desired reaction products, was dissolved in 500 µl of 4:1 DMSO:dioxane (producing a fluorescent yellow-green solution) and purified in two successive separations (250 µl each) by preparative reverse-phase high-performance liquid chromatography (HPLC) using a Gemini-NX C18 column (5 µm, 250 × 21 mm) with a gradient from 20:80 to 40:60 CH3CN:H2O (0.1% formic acid) over 30 minutes at a flow rate of 21.2 ml/min with peak elution monitored at 288 nm (Phenomenex, Torrance, CA). Separations were accomplished using a Varian ProStar HPLC system (Varian, Palo Alto, CA) equipped with ProStar 210 pumps and a ProStar 335 photodiode array detector and Galaxie Chromatography Workstation (version 1.9.3.2; Varian).
Characterization and Structural Confirmation.
All NMR experiments were conducted in acetone-d6 using either an Agilent 700 NMR spectrometer equipped with a cryoprobe (700 and 175 MHz for 1H and 13C, respectively; Agilent Technologies, Santa Clara, CA) or a JEOL ECS-400 NMR spectrometer (399.78 MHz for 1H and 100.53 MHz for 13C; JEOL, Tokyo, Japan) equipped with an auto tune 5 mm field gradient tunable Royal probe (NM-03810RO5/UPG). High-resolution electrospray ionization mass spectrometry data were collected using an electrospray ionization source coupled to a Q Exactive Plus system (Thermo Fisher Scientific, San Jose, CA) in positive and negative ionization modes via a liquid chromatography/autosampler system composed of an Acquity UPLC system (Waters).
Comparison of BLM- versus HLM-Generated Glucuronides.
Flavonolignan glucuronidation in BLMs was compared with HLMs using the previously determined optimal reaction conditions for silybin A in BLMs (Table 2). A mixture of MgCl2, BSA, alamethicin, and UDPGA was prepared in Tris-HCl buffer (100 mM, pH 7.4) at 2.5 times the final reaction concentration and split into four aliquots. The four flavonolignans were added separately from DMSO stock solutions (19.28 mg/ml) and divided to yield duplicate microcentrifuge tubes. BLMs or HLMs (0.424 mg/ml) in Tris-HCl buffer were added to produce a final reaction volume of 0.65 ml. Final concentrations of all reaction ingredients were the same as those described earlier for the large-scale reaction of silybin A. The eight mixtures were incubated at 37°C, after which 100-μl aliquots were removed and quenched with 300 µl of ice-cold MeOH at each time point (0, 2, 4, 8, 12, and 24 hours), generating 48 samples, which were stored at 4°C until analysis. Samples were prepared for and analyzed by UPLC-MS as described earlier.
Summary of optimal reaction conditions and projected cost savings to generate milk thistle flavonolignan glucuronides
Results
Synthesis of Authentic Glucuronides of Herbal Products Is Achievable.
Milk thistle flavonolignans demonstrated compound-dependent glucuronidation (Figs. 2 and 3); typically, three to five glucuronidated products were observed by mass spectrometry for each reaction. Silybin A generated five monoglucuronides, silybin B generated four monoglucuronides, isosilybin A generated four monoglucuronides, and isosilybin B generated three monoglucuronides (Figs. 2 and 3). Diglucuronide flavonolignans were detected only in trace amounts. Optimization and scaling up of incubation procedures were focused on monoglucuronide formation due to the paucity of diglucuronide conjugates. Optimal incubation times based upon maximal product formation (0.8 mM UDPGA) were 8, 8, 4, and 12 hours for silybin A, silybin B, isosilybin A, and isosilybin B, respectively (Figs. 2 and 3). Twofold molar excess of UDPGA and 0.25 mg/ml protein was optimal for all flavonolignans. Only silybin B was conjugated more rapidly by S9 than BLMs, representing a potential cost-savings measure for generating silybin B glucuronides.
Optimization of incubation time, UDPGA concentration, and BLM concentration. Silybin A (A), silybin B (B), isosilybin A (C), or isosilybin B (D) (0.4 mM) was incubated with 1 (black), 0.5 (dashed), or 0.25 (crimson) mg/ml BLMs with 4 (upper) or 0.8 (lower) mM UDPGA. Disappearance of flavonolignan starting material (left column) and appearance of glucuronides (right column) were monitored over time. Glucuronide data are presented as the sum of all monoglucuronide products formed at each time point.
Optimization of incubation time, UDPGA concentration, and S9 concentration. Silybin A (A), silybin B (B), isosilybin A (C), or isosilybin B (D) (0.4 mM) was incubated with 2 (gray), 1 (dashed), or 0.5 (black) mg/ml S9 with 4 (upper) or 0.8 (lower) mM UDPGA. Disappearance of flavonolignan starting material (left column) and appearance of glucuronides (right column) were monitored over time. Flavonolignan disappearance and glucuronide formation data using BLMs (0.25 mg/ml) are included as a comparison with microsomal systems (crimson). Glucuronide data are presented as the sum of all monoglucuronide products formed at each time point.
BLMs Represent a Cost-Effective and Efficient System for Generating Glucuronides.
Optimized reaction conditions required 4-fold less enzyme and 5-fold less UDPGA cofactor compared with standard conditions (Table 1). Use of BLMs and optimized reaction conditions yielded substantial (at least 5-fold) cost savings compared with standard reaction conditions using human-derived subcellular fractions (Table 2). Further cost savings can be realized using an alternate alamethicin source [i.e., alamethicin F50 from fungal culture MSX70741 (Ayers et al., 2012)] for synthetic reactions (∼$150 for a 25-mg scale reaction). Relatively low protein concentrations of optimized conditions would be expected to reduce compound loss to nonspecific binding and simplify purification procedures.
Large-Scale Generation of Glucuronides Facilitates Structural Characterization.
Large-scale synthesis generated three glucuronides of silybin A, which were purified by preparative HPLC (Fig. 4A). Compound 1 eluted at 13.5 minutes (15.7 mg), compound 2 at 15.1 minutes (1.6 mg), and compound 3 at 18.5 minutes (11.1 mg); the weights indicate the amounts isolated based on a reaction that started with 40 mg of silybin A. The purity of the glucuronides was measured by UPLC-MS, and the UV data were indicative of the conjugated core of flavonolignans (Fig. 4B). The structures of the three silybin A glucuronides (compounds 1–3) were assigned primarily using 1H, 13C, correlation spectroscopy, heteronuclear single-quantum correlation, and heteronuclear multiple-bond correlation (HMBC) NMR data (Fig. 5; Table 3); the data for the major product (compound 1) are discussed later, and the other glucuronides (compounds 2 and 3) were assigned in an analogous manner (Supplemental Figs. S1–S6).
Analysis of the large-scale reaction for silybin A. (A) Purification of silybin A glucuronides. Preparative HPLC UV trace (288 nm) of the reaction mixture applied to a Gemini-NX C18 column eluted using a 20–40% gradient of CH3CN/H2O (0.1% formic acid). Isolated silybin A glucuronides (1–3) are indicated in the chromatogram. (B) Stacked plots showing the selected ion chromatogram (left), UV absorbance at 288 nm (middle), and UV spectrum from 190 to 500 nm (right) of isolated compounds silybin A-7-O-β-d-glucuronide (1), silybin A-5-O-β-d-glucuronide (2), and silybin A-4′′-O-β-d-glucuronide (3).
Selected HMBC spectrum of compound 1 (400/100 MHz) in acetone-d6. The HMBC correlation from H-1′′′ to C-7 supports the assigned site of attachment of the glucuronic acid to silybin A.
NMR spectroscopic data for silybin A glucuronides (1–3) in acetone-d6
Values of δ are in ppm.
The high-resolution electrospray ionization mass spectrometry (Fig. 6), 13C NMR (Table 3), and heteronuclear single-quantum correlation data for the major product of the silybin A glucuronidation reaction (compound 1 in Figs. 4–7) indicated a molecular formula of C31H30O16, suggesting an index of hydrogen deficiency of 17. The 1H NMR spectrum (Supplemental Fig. S1) of compound 1 displayed signals for a methoxy group, eight aromatic protons, nine oxymethines, and a pair of nonequivalent protons for the oxymethylene group, whereas the 13C NMR spectrum revealed 31 carbons [20 sp2 (eight protonated) and 11 sp3 (nine tertiary, one secondary, and one primary)]. Two 1,2,4-trisubstituted aromatic rings were identified based on 1H-1H coupling constants (2JH,H = 8.3 Hz and 4JH,H = 1.8 Hz), and the corresponding 13C NMR shifts (δC 144.7, 145.2, 148.5, and 148.0 for C-3′, C-4′, C-3′′, and C-4′′, respectively) required that they each be 1,2-dioxygenated. An additional tetrasubstituted aromatic ring, which was 1,3,5-trioxygenated, was also observed. Finally, two carbonyl carbons (δC 199.0 and 170.0 for C-4 and C-7′′′, respectively) accounted for the remaining sp2 carbons.
Stacked plots showing the high-resolution electrospray ionization mass spectrometry for positive (left) and negative (right) ionization modes of isolated compounds silybin A-7-O-β-d-glucuronide (1), silybin A-5-O-β-d-glucuronide (2), and silybin A-4′′-O-β-d-glucuronide (3). Values in parentheses represent the difference between the measured versus calculated mass for the indicated molecular formula; values within 5 ppm are considered valid.
Structures and diagnostic HMBC correlations of silybin A glucuronides (1–3).
The structure of the silybin A portion of compound 1 was established by one-dimensional and two-dimensional NMR data and through comparisons to the literature (Lee and Liu, 2003). The structural data from the 1H NMR spectrum of compound 1 (Supplemental Fig. S1; Table 3) were consistent with the characteristics of a 5,7-dioxygenated flavanonol by the signals at δH 6.20 (1H, s, H-6), 6.20 (1H, s, H-8), 5.18 (1H, d, J = 11.7 Hz, H-2), and 4.75 (1H, d, J = 11.7 Hz, H-3). The aromatic proton signals at δH 7.17 (1H, d, J = 1.8 Hz, H-2′), 6.97 (1H, d, J = 8.3 Hz, H-5′), and 7.10 (1H, dd, J = 8.3, 1.8 Hz, H-6′) could be attributed to the B-ring unit of silybin A. The aromatic proton signals at δH 7.15 (1H, d, J = 1.8 Hz, H-2′′), 6.88 (1H, d, J = 8.3 Hz, H-5′′), and 6.98 (1H, dd, J = 8.3, 1.8 Hz, H-6′′) could be attributed to the E-ring unit of silybin A. The proton signals at δH 5.01 (1H, d, J = 8.3 Hz, H-7′′), δH 4.17 (1H, ddd, J = 8.3, 4.4, 2.4 Hz, H-8′′), δH 3.76 (1H, dd, J = 12.4, 2.4 Hz, H-9′′a), and δH 3.52 (1H, dd, J = 12.4, 4.4 Hz, H-9′′b) along with the aromatic proton signals for the E-ring unit could be attributed to the coniferyl alcohol moiety of silybin A.
Additional proton (5 oxymethines) and carbon signals in the spectra of compound 1 indicated condensation of silybin A with a glucuronide moiety. The anomeric H-1′′′ of the glucuronide portion of compound 1 was readily identifiable due to its chemical shift, coupling pattern, and coupling constant (δH 5.27, d, 7.8 Hz). The two main structural units of compound 1 were linked on the basis of an HMBC correlation of H-1′′′ to C-7 (Figs. 5 and 6; Supplemental Fig. S2). The presence of a sharp peak at approximately δH 11.5, due to intramolecular H-bonding between the phenolic proton at C-5 to the C-4 carbonyl, further verified the linkage at the C-7 position.
The structures of compounds 2 and 3 were elucidated in a similar manner (Fig. 6; Supplemental Figs. S3–S6; Table 3). With respect to the characterization of flavonolignans by NMR, in general (Napolitano et al., 2013), this study varied in one major aspect. Based on earlier research that derivatized the flavonolignans (Sy-Cordero et al., 2013; Althagafy et al., 2013), the phenol moieties were hypothesized to be glucuronidated readily; indeed, other researchers have reported this observation (Kren et al., 2000; Jančová et al., 2011; Křen et al., 2013; Charrier et al., 2014). Determining the position of glucuronidation required an HMBC correlation from the anomeric proton back to the carbon attached to the phenol (Fig. 6). However, such a correlation would be somewhat inconclusive if the phenolic protons at the nonglucuronidated positions were not observed, which was a challenge at the 5 position due to chelation with the carbonyl and proton exchange in protic solvents, e.g., MeOH-d4. Thus, all NMR spectra were acquired in acetone-d6 such that the phenolic protons could be observed readily. For compounds 2 and 3, the 13C NMR data (Supplemental Figs. S3 and S5; Table 3) displayed 31 signals, which were six more than silybin A due to the resonances from the glucuronic acid motif at C-1′′′ through C-7′′′. The anomeric C-1′′′ appeared at approximately δC 102 and the anomeric H-1′′′ at approximately δH 5.0 as a doublet with a coupling constant at approximately 7 Hz, characteristic of the β-anomer (Pearson et al., 2005; Walker et al., 2007). Data from the HMBC experiments (Supplemental Figs. S4 and S6) demonstrated the connectivity of the glucuronide moiety, where correlations were observed between H-1′′′ to C-5 and H-1′′′ to C-4′′ for compounds 2 and 3, respectively.
BLMs Generate Human-Relevant Glucuronides with Enhanced Efficiency.
Compared with HLMs, BLMs converted all four flavonolignans to their respective glucuronides more rapidly and completely (Fig. 8). Silybin A and silybin B glucuronides were produced similarly by both BLMs and HLMs. The predominant glucuronides of isosilybin A and isosilybin B generated using BLMs (elution time 3.1 minutes) were opposite of the predominant glucuronides generated using HLMs (elution time 2.8 minutes) (Fig. 8). Overall, BLMs generated glucuronide profiles similar to those generated by HLMs (Fig. 8).
Comparison of BLM- versus HLM-generated glucuronides. Silybin A (A), silybin B (B), isosilybin A (C), or isosilybin B (D) (0.4 mM) was incubated with 0.25 mg/ml BLMs (crimson) or HLMs (black) with 0.8 mM UDPGA. (Left panel) Disappearance of flavonolignan starting material (dashed lines) and appearance of glucuronides (solid lines) were monitored over time. Glucuronide data are presented as the sum of all monoglucuronide products formed at each time point. (Right panel) UPLC-MS chromatograms of flavonolignan glucuronides (659 m/z) generated following 12-hour incubation with BLMs (upper) or HLMs (lower).
Discussion
Partnerships between clinical pharmacologists and natural products chemists with divergent yet complementary skill sets are imperative to advance the understanding of herbal product constituent disposition and pharmacology. Complete characterization of herbal product metabolites will enable identification of sites of metabolic lability that could be targeted to optimize pharmacokinetics of therapeutically promising constituents. These studies would enable structure-activity relationships that lay the groundwork for application of rational drug design concepts to compounds of natural origin.
Metabolite characterization can identify the potential for bioactivation and toxicity of a molecule that may occur in vivo. Identifying and characterizing metabolites with appreciable human exposure is critical to a complete understanding of the pharmacology of a particular herbal product. Minor metabolites may also be important to overall activity, as they can be unique in both disposition and action. Individual constituents within complex herbal product mixtures can be differentially metabolized. Milk thistle flavonolignans, although structurally similar (i.e., constitutional isomers), exemplify this consideration, as evidenced by the in vitro depletion curves demonstrating compound-dependent metabolism in BLMs (Figs. 2 and 3) and differential metabolism in BLMs compared with HLMs (Fig. 8). Although evaluation of isolated constituents can simplify investigation, the possibility exists that one constituent may modulate metabolism of another both in vitro and in vivo. This possibility underscores the importance of evaluating both individual constituents and complex mixtures to ascertain the metabolic fates in each scenario; such experiments would be exceedingly difficult without the availability of authentic metabolite standards.
Typical chemoenzymatic glucuronidation reactions include a 10-fold molar excess of UDPGA and ≥1 mg/ml of microsomal protein. Saturating UDPGA concentrations (5–10 mM) are routinely used when recovering in vitro kinetic parameters using HLMs (Km ∼ 0.2–1.3 mM) (Court et al., 2001) to eliminate UDPGA availability as a potential rate-limiting step. Application of these relatively high concentrations to synthetic reactions assumes that excess cofactor and enzyme enhance yield. The current work demonstrated that using less of these reagents for chemoenzymatic synthesis is an effective cost-saving measure, and that excess enzyme and reagents can reduce yields via what appear to be reverse catalases. Saccharolactone, a β-glucuronidase inhibitor, is often added to microsomal incubations to prevent enzymatic degradation of glucuronides. However, hepatic microsomal preparations are not expected to be substantially impacted by adding saccharolactone (Oleson and Court, 2008), and degradation of glucuronides under optimized conditions was minimal (Fig. 2). Consequently, saccharolactone was not included in the optimized reaction. Detergents could be explored as an additional cost-savings measure to replace alamethicin, but the impact of detergent selection and concentration on glucuronide yield would need to be evaluated (Soars et al., 2003). Albumin was added to all microsomal incubations based upon reports of substrate-dependent decreases in Km and increases in Vmax (Rowland et al., 2008; Manevski et al., 2011, 2013; Gill et al., 2012; Walsky et al., 2012), which are believed to reflect binding to inhibitory fatty acids (Rowland et al., 2007). Selection of additives that can modulate microsomal UDP-glucuronosyl transferase activity is an important consideration that can be addressed in microscale optimization reactions prior to scaling up synthetic procedures.
Generation of flavonolignan glucuronides on the desired scale initially appeared to be overtly cost-prohibitive based on literature protocols. However, relatively inexpensive microscale optimization reduced the total costs by at least an order of magnitude, making these previously cost-prohibitive experiments (thousands of U.S. dollars per reaction) feasible (Table 2). Cost-effective generation of large (multimilligram) quantities are needed to support detailed investigation of the kinetics and pharmacology of herbal product metabolites across multiple laboratories. However, cost savings could be realized with microscale (<1 mg) incubations followed by quantitative NMR to generate small, yet useful, quantities of metabolites using approaches similar to those described for the synthesis of drug metabolites (Walker et al., 2011, 2014). Flavonolignan cost was not factored into total cost estimates, as these compounds are isolated in house, precluding meaningful comparisons to alternate herbal product constituents that can vary widely in cost and availability. Synthetic processes using isolated constituent starting materials enhance optimization of specific reaction conditions and simplify metabolite purification, leading to enhanced yield. Enhanced efficiency is crucial to reduce waste of limited individual constituents that may take weeks to months to isolate and purify in quantities sufficient for metabolite synthesis. Deconvolution of complex parent and metabolite mixtures encountered in both in vitro and clinical samples can be overcome partially using isolated constituents to simplify the diversity of generated metabolites in a particular system.
Bioanalysis of herbal product glucuronides historically has relied on indirect methods using enzymatic (e.g., β-glucuronidase) incubations to cleave the conjugates, followed by analysis of the parent compound. Indirect quantification methods increase sample preparation time and frequently are plagued by poor accuracy and reproducibility resulting from incomplete enzymatic hydrolysis and instability of metabolites under enzymatic conditions. Moreover, simultaneous cleavage of multiple glucuronides present in a given system precludes quantitative description of individual glucuronides. A lack of commercially available glucuronide standards necessitates chemoenzymatic synthesis and purification of glucuronide metabolites for analytical development and application. The reaction scales described in the current work are capable of supporting bioanalytical method development, metabolism inquiries, and pharmacology evaluations.
In summary, techniques to efficiently and cost-effectively generate and purify authentic glucuronide metabolites of a model herbal product have been developed (Fig. 1). This approach was bolstered by a partnership between clinical pharmacology and natural products chemistry expertise. These techniques are essential to improve the understanding of mechanisms underlying herbal product disposition and action. Authentic metabolites can be used to provide quantitative description of metabolite formation rates, exposure, and activity. Quantitative description of herbal product disposition and action could lead to improved translation of natural product discoveries from bench to bedside.
Acknowledgments
The authors thank Dr. Cedric Pearce of Mycosynthetix, Inc. for access to fungal culture MSX70741 and Dr. Daneel Ferreira of the University of Mississippi for helpful discussions. M.F.P. dedicates this article to Dr. David P. Paine.
Authorship Contributions
Participated in research design: Gufford, Graf, Oberlies, Paine.
Conducted experiments: Gufford, Graf, Paguigan.
Contributed new reagents or analytic tools: Graf, Oberlies.
Performed data analysis: Gufford, Graf, Paguigan, Oberlies, Paine.
Wrote or contributed to the writing of the manuscript: Gufford, Graf, Paguigan, Oberlies, Paine.
Footnotes
- Received June 26, 2015.
- Accepted August 26, 2015.
This work was supported by the National Institutes of Health National Institute of General Medical Sciences [Grant R01-GM077482-S1]. B.T.G. was supported by fellowships awarded by the American Foundation for Pharmaceutical Education and the James and Diann Robbers Student Research Fund. B.T.G. is currently supported by the National Institute of General Medical Sciences [Grant T32-GM008425]. Alamethicin F50 was isolated as part of Program Project Grant P01-CA125066 from the National Institutes of Health National Cancer Institute. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of General Medical Sciences or the National Institutes of Health.
↵
This article has supplemental material available at dmd.aspetjournals.org.
Abbreviations
- BLM
- bovine liver microsome
- BSA
- bovine serum albumin
- DMSO
- dimethylsulfoxide
- HLM
- human liver microsome
- HMBC
- heteronuclear multiple-bond correlation
- HPLC
- high-performance liquid chromatography
- MS
- mass spectrometry
- UDPGA
- UDP glucuronic acid
- UPLC
- ultra-performance liquid chromatography
- Copyright © 2015 by The American Society for Pharmacology and Experimental Therapeutics