Abstract

The natural product goldenseal is a clinical inhibitor of CYP3A activity, as evidenced by a 40%–60% increase in midazolam area under the plasma concentration versus time curve (AUC) after coadministration with goldenseal. The predominant goldenseal alkaloids berberine and (−)-β-hydrastine were previously identified as time-dependent CYP3A inhibitors using human liver microsomes. Whether these alkaloids contribute to the clinical interaction, as well as the primary anatomic site (hepatic vs. intestinal) and mode of CYP3A inhibition (reversible vs. time-dependent), remain uncharacterized. The objective of this study was to mechanistically assess the pharmacokinetic goldenseal-midazolam interaction using an integrated in vitro–in vivo–in silico approach. Using human intestinal microsomes, (−)-β-hydrastine was a more potent time-dependent inhibitor of midazolam 1’-hydroxylation than berberine (KI and kinact: 8.48 μM and 0.041 minutes−1, respectively, vs. >250 μM and ∼0.06 minutes−1, respectively). Both the AUC and Cmax of midazolam increased by 40%–60% after acute (single 3-g dose) and chronic (1 g thrice daily × 6 days) goldenseal administration to healthy adults. These increases, coupled with a modest or no increase (≤23%) in half-life, suggested that goldenseal primarily inhibited intestinal CYP3A. A physiologically based pharmacokinetic interaction model incorporating berberine and (−)-β-hydrastine successfully predicted the goldenseal-midazolam interaction to within 20% of that observed after both chronic and acute goldenseal administration. Simulations implicated (−)-β-hydrastine as the major alkaloid precipitating the interaction, primarily via time-dependent inhibition of intestinal CYP3A, after chronic and acute goldenseal exposure. Results highlight the potential interplay between time-dependent and reversible inhibition of intestinal CYP3A as the mechanism underlying natural product-drug interactions, even after acute exposure to the precipitant.

SIGNIFICANCE STATEMENT Natural products can alter the pharmacokinetics of an object drug, potentially resulting in increased off-target effects or decreased efficacy of the drug. The objective of this work was to evaluate fundamental mechanisms underlying the clinically observed goldenseal-midazolam interaction. Results support the use of an integrated approach involving established in vitro assays, clinical evaluation, and physiologically based pharmacokinetic modeling to elucidate the complex interplay between multiple phytoconstituents and various pharmacokinetic processes driving a drug interaction.

Introduction

The botanical natural product goldenseal [Hydrastis canadensis L. (Ranunculaceae)], a perennial herb native to North America, is commonly used to self-treat the common cold and other respiratory infections, allergic rhinitis, and digestive complications such as diarrhea and constipation (Hobbs, 1990) (https://www.nccih.nih.gov/health/goldenseal). Gurley and colleagues previously reported goldenseal to precipitate clinically significant interactions with midazolam, a drug that undergoes extensive metabolism via cytochrome P450 (P450) 3A–mediated hydroxylation (Kronbach et al., 1989; Tian et al., 2019). The metabolic ratio, as determined by the plasma 1’-hydroxymidazolam/midazolam ratio measured at a single blood collection time point, decreased by 40% after a 4-week administration of goldenseal root extract (900 mg thrice daily) to 12 healthy adult participants (Gurley et al., 2005). The decreased CYP3A activity was recapitulated by a 1.6-fold increase in midazolam area under the plasma concentration versus time curve (AUC) after 2 weeks of goldenseal root extract administration (1.3 g thrice daily) to 16 healthy adults (Gurley et al., 2008). These observations clearly indicate that goldenseal is a clinical inhibitor of CYP3A. However, the primary precipitant, mode of inhibition, and site of interaction remain uncharacterized.

Major phytoconstituents of goldenseal include the isoquinoline alkaloids berberine and (−)-β-hydrastine. These alkaloids contain methylenedioxyphenyl moieties, which are known structural alerts for potential time-dependent inhibition (TDI) of P450s (Kalgutkar et al., 2007). Both berberine and (−)-β-hydrastine demonstrated TDI of CYP3A activity (midazolam 1’-hydroxylation) in human liver microsomes (HLMs) (Chatterjee and Franklin, 2003; McDonald et al., 2020). The kinact/KI ratios for berberine and (−)-β-hydrastine were 1.3 and 2.0 ml/min/μmol, respectively (McDonald et al., 2020). Berberine also activated midazolam 1’-hydroxylation in both HLMs and recombinant CYP3A5, but not recombinant CYP3A4, as shown by a reduction in Km(app). In addition to the liver, the intestine represents a potential site for first-pass drug interactions. Based on the presumed high concentrations of exogenous precipitants at this anatomic site, along with the fact that midazolam undergoes extensive first-pass metabolism (Paine et al., 1996; Thummel et al., 1996), modulation of intestinal CYP3A may contribute to the observed goldenseal-midazolam interaction.

Physiologically based pharmacokinetic (PBPK) modeling is a dynamic tool that can provide mechanistic insights into biologic processes governing an observed pharmacokinetic interaction by integrating organ physiology with compound-dependent physicochemical and pharmacokinetic parameters (Sager et al., 2015). Unlike for drug-drug interactions, the application of PBPK modeling to natural product-drug interactions remains underutilized due to the phytochemical complexity of natural products, the contents of phytoconstituents in marketed products varying between different manufacturers, and the lack of robust physicochemical and pharmacokinetic data for most phytoconstituents. A PBPK model for goldenseal was previously developed and verified to simulate the observed goldenseal-midazolam interaction (Adiwidjaja et al., 2022). Although the model provided a reliable estimate of the magnitude of the interaction, elucidating mechanisms underlying this natural product-drug interaction was not a premise of the study.

The objective of the current study was to elucidate the fundamental mechanism(s) underlying the goldenseal-midazolam interaction using a comprehensive integrated approach involving established in vitro assays, clinical evaluation, and PBPK modeling. The aims were to 1) describe the inhibition kinetics of berberine and (−)-β-hydrastine toward intestinal CYP3A activity using human intestinal microsomes (HIMs), 2) evaluate the pharmacokinetics of midazolam before and after oral administration of a well-characterized goldenseal product to healthy adult participants, and 3) develop a PBPK model to gain further mechanistic insights into the goldenseal-mediated drug interaction (i.e., contributions between reversible inhibition and TDI by berberine and/or (−)-β-hydrastine and the predominant anatomic site of interaction). Results illustrate how this integrated approach can be used to elucidate the complex interplay between multiple phytoconstituents and multiple pharmacokinetic processes driving an observed natural product-drug interaction.

Materials and MethodsChemicals and Reagents

Midazolam, 1’-hydroxymidazolam, 4-hydroxymidazolam, midazolam 1’-O-glucuronide, midazolam 4-O-glucuronide, midazolam N-glucuronide, alprazolam, β-nicotinamide adenine dinucleotide phosphate (NADPH), ketoconazole, atenolol, berberine chloride, (−)-β-hydrastine (USP reference standard), and penicillin/streptomycin (10,000 units/10 mg per ml) were purchased from Sigma-Aldrich (St. Louis, MO). HIMs (10 mg/ml; H0610.I, lot 1610314), pooled from 15 donors of mixed sex, were purchased from Xenotech (Kansas City, KS). Fasted state simulated intestinal fluid (FaSSIF) was acquired from Biorelevant (London, UK). Caco-2 cells were procured from American Type Culture Collection (ATCC) (Manassas, VA). All other chemicals and reagents were of analytical grade.

IC50 Shift Experiment for the Predominant Goldenseal Alkaloids Using HIMs

The primary mixture contained HIMs (0.05 mg/ml), potassium phosphate buffer (100 mM, pH 7.4), and either berberine or (−)-β-hydrastine (0.1–500 μM) incubated in the presence and absence of NADPH (1 mM). After 30 minutes at 37°C, 196 μl of the mixtures were transferred to a secondary reaction, which contained 2 μl midazolam (final concentration: 4 μM) and either 2 μl potassium phosphate buffer or 2 μl NADPH (final concentration: 1 mM). The final methanol concentration in all incubation mixtures was ≤1% (v/v). Reactions were terminated after 10 minutes by removing 100 μl of the mixture and adding to 200 μl of ice-cold methanol containing alprazolam (100 nM) as internal standard. After centrifugation (3800 g × 10 minutes), the supernatant was analyzed for 1'-hydroxymidazolam using ultrahigh-performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) (Supplemental Table 1). IC50 values were recovered from percent control activity versus log-transformed concentration data via nonlinear regression using GraphPad Prism 9 (GraphPad Software Inc., La Jolla, CA).

Time-Dependent Inhibition of Midazolam 1’-Hydroxylation in HIMs by Goldenseal Alkaloids

HIMs (0.05 mg/ml) diluted in potassium phosphate buffer were incubated with berberine (25–500 μM), (−)-β-hydrastine (3–90 μM), or 6,7-dihydroxybergamottin (1, 0.5 μM); the final concentration of methanol was ≤1% (v/v). Primary incubation mixtures were equilibrated to 37°C, after which the reactions were initiated by adding NADPH (1 mM). After 0, 2, 5, 10, 15, and 20 minutes, 198 μl of the mixture were added to a tube containing 2 μl of midazolam (final concentration: 40 μM). The secondary incubations were terminated 5 minutes later with one volume of ice-cold methanol containing alprazolam (100 nM) as internal standard. After centrifugation (3800 g × 10 minutes), the supernatant was analyzed for 1'-hydroxymidazolam using UHPLC-MS/MS (Supplemental Table 1). KI and kinact estimates were recovered from kinact,app versus [I] data, which were analyzed via nonlinear regression using Phoenix WinNonlin (v7.0; Certara USA, Princeton, NJ) as described (Brantley et al., 2013; Tanna et al., 2021).

Allosterism of Midazolam Hydroxylation in HIMs by Berberine

HIMs (0.05 mg/ml) diluted in potassium phosphate buffer were incubated with berberine (5–100 μM) and midazolam (0.5–250 μM). Control mixtures contained ≤1% (v/v) methanol in place of berberine. Reactions proceeded for 5 minutes after initiating with NADPH (1 mM). Samples were processed as described above and analyzed for 1’-hydroxymidazolam and 4-hydroxymidazolam using UHPLC-MS/MS (Supplemental Table 1). Relative changes in Vmax and Km were visualized from initial velocity versus [S] data.

Dissolution Profiles of the Predominant Goldenseal Alkaloids

The goldenseal product administered to the clinical study participants was extensively characterized as described (Kellogg et al., 2020; Wallace et al., 2020). In brief, each gram of the product (Now Foods, Bloomingdale, IL) contained 29.8 mg of berberine and 25.4 mg of (−)-β-hydrastine. The dissolution profiles for berberine and (−)-β-hydrastine were determined in fasted state simulated intestinal fluid (FaSSIF) using Apparatus 2 (708-DS Dissolution Apparatus; Agilent Technologies, Santa Clara, CA). To mimic the physiologic conditions established in the previously published clinical study (Nguyen et al., 2021), two capsules of the goldenseal product (500 mg/capsule) were added to 250 ml of FaSSIF in glass vessels attached with a paddle rotating at 50 rpm. Samples (∼1 ml) of the mixture were collected at 0, 5, 10, 15, 20, 30, 60, and 120 minutes and filtered through Millex hydrophilic PTFE 0.20 μm membrane filters (Merck Millipore Ltd, Tullagreen, Ireland). Samples were diluted 1:500 in methanol containing the internal standard alprazolam (100 nM) prior to UHPLC-MS/MS analysis for berberine and (−)-β-hydrastine (Supplemental Table 1).

Permeability Assessment of Major Goldenseal Alkaloids

The apparent permeability (Papp) for berberine and (−)-β-hydrastine (10 μM) was determined using Caco-2 cell monolayers, which were developed with cells of low passage count (12–30) using established methods (Hubatsch et al., 2007). In brief, Caco-2 cells were seeded onto 12-mm Transwell polyester membrane inserts (Corning 3460; Corning, NY) at a density of 3 × 105 cells/cm2 and grown for 21–28 days in Dulbecco’s modified Eagle’s medium (Corning 10-009-CV) supplemented with 20% fetal bovine serum (ATCC 30-2020) and penicillin/streptomycin (final concentration: 100 units/0.1 mg/ml). Papp in the apical to the basolateral direction was assessed for each test article. Papp for berberine, (−)-β-hydrastine, and the low permeability control atenolol (10 μM) was determined after 30 and 60 minutes of incubation, whereas Papp for the high permeability control midazolam (10 μM) was determined after 15 and 30 minutes. Samples were diluted in methanol containing alprazolam (100 nM; internal standard) prior to UHPLC-MS/MS analysis (Supplemental Table 1). To ensure the integrity of cell monolayers, only those with a net transepithelial electrical resistance measurement ≥220 were used.

Clinical Pharmacokinetic Goldenseal-Midazolam Interaction Study

The clinical protocol and consent form were approved by the Washington State University (WSU) Institutional Review Board (IRB 16620) and registered in the ClinicalTrials.gov database (NCT03772262). All clinical activities were conducted at the WSU Clinical Research Unit on the Health Sciences Campus in accordance with 45 CFR 46 and adherence to Good Clinical Practice guidelines. The study was initially conducted in an open-label, two-arm (baseline and chronic goldenseal exposure), fixed sequence crossover manner as described (Nguyen et al., 2021). In brief, during the baseline arm, healthy adult participants (eight males, eight females) were orally administered 2.5 mg midazolam syrup (West-Ward Pharmaceuticals, Eatontown, NJ) as a sensitive CYP3A substrate in a probe drug cocktail. After at least a 1-week washout, participants were administered 1 g goldenseal thrice daily for 6 days (chronic exposure). On the sixth day, they were administered the probe drug cocktail containing midazolam ∼30 minutes after the first dose of goldenseal. For both arms, plasma and urine were collected up to 96 and 24 hours, respectively, after midazolam administration. A subset of these participants (four males, four females) completed an exploratory study to determine whether a single large dose of goldenseal (acute exposure) would precipitate a pharmacokinetic interaction with midazolam. During this arm, participants were administered a single 3-g dose of goldenseal ∼30 minutes prior to midazolam syrup (2.5 mg). Plasma and urine were collected up to 12 hours after midazolam administration in the acute exposure arm. During the three study arms, vital signs (blood pressure, pulse, oxygen saturation) were recorded periodically throughout the inpatient day. Plasma and urine were analyzed for midazolam, midazolam metabolites, berberine, and (−)-β-hydrastine via UHPLC-MS/MS (Supplemental Table 1).

Pharmacokinetic Analysis

The pharmacokinetics of midazolam, midazolam metabolites, berberine, and (−)-β-hydrastine were determined by noncompartmental analysis using Phoenix WinNonlin as described (Tian et al., 2019). The following were recovered: plasma AUC from time zero to the last measured concentration (AUC0–tlast), AUC from time zero to infinity (AUCinf), terminal half-life (t1/2), maximum plasma concentration (Cmax), time to reach Cmax (tmax), oral clearance (Cl/F), and renal clearance (CLR). Midazolam AUC0–tlast was determined up to 12 hours after administration due to concentrations being below the limit of quantification.

Statistical and Power Analysis

The clinical study was initially designed and powered to assess the effects of chronic goldenseal exposure on midazolam pharmacokinetics as previously reported (Nguyen et al., 2021); the protocol was later amended to include the exploratory acute goldenseal exposure arm. The primary endpoints included the chronic goldenseal exposure-to-baseline ratios of midazolam AUC0–tlast and AUCinf with a predefined no effect range of 0.80–1.25. Secondary endpoints included Cmax, t1/2, Cl/F, and CLR of midazolam after chronic goldenseal exposure; AUC0–tlast, AUCinf, Cmax, t1/2, and CLR of midazolam metabolites after chronic goldenseal exposure; and all of these pharmacokinetic measures for midazolam and metabolites after acute goldenseal exposure. Statistical analyses were performed on all secondary endpoints using Wilcoxon matched-pairs signed-rank test.

PBPK Modeling of the Goldenseal-Midazolam Interaction

The PBPK model was developed using the modeling and simulation platform Simcyp (v22; Certara USA, Princeton, NJ); the model structure has been extensively described elsewhere (Jamei et al., 2009, 2013; Pathak et al., 2017). The default midazolam model within the software was applied to evaluate CYP3A-mediated drug interactions (Supplemental Methods). The inhibitor models for berberine and (−)-β-hydrastine were developed using physicochemical properties, in vitro–derived kinetic parameters, and external pharmacokinetic data (Supplemental Table 2) and verified with the current clinical pharmacokinetic study. Based on the dissolution and permeability results, the first-order absorption model was selected for the berberine PBPK model, whereas the advanced dissolution, absorption, and metabolism (ADAM) model was adapted for the (−)-β-hydrastine PBPK model. Due to limitations of the modeling and simulation platform, only one inhibition value (obtained from either HIMs or HLMs) can be used to characterize the interaction potential for a given inhibitor; thus, the value correlating with the highest inhibition potency was used for model development. Simulations were conducted in a virtual population with demographics matching those of the participants enrolled in the current clinical study (Supplemental Table 3). Model predictions were considered successful if the predicted AUC0–tlast and Cmax were within 2-fold of observed values.

ResultsIC50 Shift for Berberine and (−)-β-Hydrastine in HIMs

An IC50 shift for berberine was not observed due to apparent activation of midazolam 1’-hydroxylation in the absence of NADPH in the primary incubation (Fig. 1A). In contrast, an ∼15-fold leftward shift in IC50 was observed for (−)-β-hydrastine after preincubation with NADPH (Fig. 1B).

Fig. 1.

IC50 for berberine (A) and (−)-β-hydrastine (B) preincubated with human intestinal microsomes in the presence and absence of NADPH. The index reaction, midazolam 1’-hydroxylation (4 μM midazolam), was used to assess CYP3A activity. Symbols and error bars denote means and standard deviations, respectively, of triplicate incubations. Lines denote nonlinear least-squares regression of the percent control activity versus log-transformed concentration data.

Time-Dependent Inhibition of Midazolam 1’-Hydroxylation by Berberine and (−)-β-Hydrastine in HIMs

Both berberine and (−)-β-hydrastine exhibited concentration-dependent and time-dependent inhibition of midazolam 1’-hydroxylation (Fig. 2). KI, kinact, and kinact/KI for (−)-β-hydrastine were 8.48 ± 3.09 μM, 0.041 ± 0.0034 minutes−1, and 4.83 ml/min/μmol, respectively. A robust estimate of KI and kinact for berberine could not be obtained because testing concentrations higher than 500 μM exceeded the solubility limit. As such, the maximum kinact,app (0.061 ± 0.001 minutes−1) was used as a conservative estimate of kinact, of which the corresponding KI was 250 μM. Using these values, kinact/KI for berberine was estimated to be 0.244 ml/min/μmol.

Fig. 2.

Time-dependent inhibitory effects of berberine (A) and (−)-β-hydrastine (B) on midazolam 1’-hydroxylation in human intestinal microsomes. Symbols and lines denote means of duplicate incubations and linear regression of initial monoexponential decline in activity, respectively. Corresponding graphs on the right (C and D) depict the rate of CYP3A inactivation as a function of inhibitor concentration. Symbols denote observed inactivation rates at each inhibitor concentration, and the curve (D) denotes nonlinear least-squares regression of the observed data. A robust estimate of KI and kinact for berberine (C) could not be recovered because testing concentrations higher than 500 μM exceeded the solubility limit.

Allosterism of Midazolam Hydroxylation by Berberine in HIMs

Using a reversible inhibition experimental design, berberine activated 1’-hydroxymidazolam formation at midazolam concentrations ≤5 μM but inhibited the reaction at concentrations ≥50 μM (Fig. 3A). In contrast, berberine inhibited midazolam 4-hydroxylation throughout the tested concentration range (5–100 μM) (Fig. 3B). The 1’-hydroxymidazolam/4-hydroxymidazolam ratio decreased with increasing midazolam concentrations (Fig. 3C). Berberine showed a concentration-dependent increase in this ratio, particularly at lower midazolam concentrations (≤5 μM).

Fig. 3.

Effects of berberine (5–100 μM) on the rate of midazolam 1’-hydroxylation (A) and 4-hydroxylation (B) in human intestinal microsomes at midazolam concentrations ranging from 5 to 250 μM. The corresponding ratios of 1’-hydroxymidazolam to 4-hydroxymidazolam (C) are shown for midazolam concentrations up to 250 μM.

Dissolution and Permeability of Berberine and (−)-β-Hydrastine

More than 80% of total berberine and (−)-β-hydrastine content (29.8 and 25.4 mg, respectively) in two goldenseal capsules dissolved in FaSSIF within 20 minutes (Fig. 4, A and B). The Papp for berberine and (−)-β-hydrastine was 0.18 × 10−6 and 9.72 × 10−6 cm/s, respectively (Fig. 4C).

Fig. 4.

Dissolution profiles for berberine (A) and (−)-β-hydrastine (B) observed in FaSSIF over the course of 2 hours and apical-to-basolateral translocation (C) of berberine and (−)-β-hydrastine in Caco-2 cell monolayers grown on Transwell inserts. Midazolam and atenolol served as high and low permeability controls, respectively. Symbols and error bars denote means and standard deviations, respectively, of three replicates (A and B); areas circled by dashes represent the methylenedioxyphenyl moiety, a structural alert for time-dependent inhibition. Horizontal lines (C) denote means of two to three replicates from three separate experiments.

Effects of Goldenseal on the Pharmacokinetics of Midazolam and Its Metabolites in Healthy Adult Participants

Midazolam pharmacokinetics observed in the baseline and chronic goldenseal exposure arms have been previously reported (Nguyen et al., 2021). Acute administration of goldenseal resulted in similar changes to midazolam pharmacokinetics as those after chronic administration, albeit the acute goldenseal exposure arm involved only eight of the 16 participants from the baseline/chronic arm (Fig. 5). The geometric mean ratio of midazolam AUC0–12h (AUC from 0 to 12 hours), AUCinf, and Cmax increased after chronic (30%–40%) and acute (40%–60%) administration (Table 1). Chronic, but not acute, administration modestly increased midazolam half-life. Goldenseal had minimal to no effects on midazolam tmax and renal clearance after both chronic and acute administration.

Fig. 5.

Plasma concentration-time profiles for midazolam (A), midazolam-N-glucuronide (B), 1’-hydroxymidazolam (C), midazolam 1’-O-glucuronide (D), 4-hydroxymidazolam (E), and midazolam 4-O-glucuronide (F) after oral administration of a single dose of midazolam (2.5 mg) to healthy adult participants in the absence of goldenseal (gray) and after chronic (blue) and acute (orange) exposure to goldenseal. Symbols denote geometric means of the indicated number of participants except for midazolam 4-O-glucuronide, which represents arithmetic means of two participants. Error bars are excluded for visual clarity; the reader is referred to Table 1 for the 90% confidence intervals of various pharmacokinetic parameters. The profiles for midazolam at baseline and after chronic goldenseal administration (A) are reprinted with permission from the American Society for Clinical Pharmacology and Therapeutics from Nguyen et al. (2021).

TABLE 1

Pharmacokinetics of midazolam and its metabolites after oral administration of midazolam (2.5 mg) alone and after chronic (1 g thrice daily for 6 days) or acute (a single dose of 3 g) exposure to goldenseal

Data represent geometric means [90% confidence intervals].

The effects of goldenseal on the pharmacokinetics of midazolam metabolites varied between chronic and acute administration (Table 1). Acute goldenseal increased the AUC of 1’-hydroxymidazolam, leading to no change in the metabolite-to-parent AUC ratio, whereas chronic goldenseal had no effect on the AUC of 1’-hydroxymidazolam, leading to a decrease in the metabolite-to-parent AUC ratio. Both chronic and acute goldenseal had no effect on the renal clearance of 1’-hydroxymidazolam. Acute goldenseal decreased the Cmax and metabolite-to-parent AUC ratio of midazolam 1’-O-glucuronide, whereas chronic goldenseal only decreased the midazolam 1’-O-glucuronide metabolite-to-parent AUC ratio. Chronic and acute goldenseal had variable effects on the minor metabolites 4-hydroxymidazolam, midazolam 4-O-glucuronide, and midazolam N-glucuronide. The change in midazolam pharmacokinetics after chronic goldenseal exposure for the eight participants who completed all three arms of the study (Supplemental Table 4) was comparable to the change after acute exposure.

CYP3A5 Genotype

The demographics of the participants in the clinical study are summarized (Supplemental Table 3). Of the 16 participants who completed the study, 14 were CYP3A5*3*3 homozygotes, one a CYP3A5*1*1 homozygote, and one a CYP3A5*1*3 heterozygote. Because none of the participants carried the nonfunctional variants CYP3A5*6 and CYP3A5*7, a subanalysis of the pharmacokinetic data was not necessary.

PBPK Modeling and Simulation of the Goldenseal-Midazolam Interaction

The developed PBPK models predicted the AUC0–12.5h and Cmax of berberine and (−)-β-hydrastine to within 2-fold of those observed after a single dose and multiple doses of goldenseal (Table 2). In addition, the simulated plasma concentration-time profiles for berberine and (−)-β-hydrastine accurately described the concentrations observed after chronic and acute goldenseal exposure in healthy adult participants (Fig. 6). However, the extrapolated area (AUCtlast–inf) for the observed berberine AUCinf exceeded 20% of AUC0–tlast, resulting in an underprediction of AUCinf (Table 2) due to lack of a robust estimate of the terminal half-life in model simulations.

TABLE 2

Predicted and observed pharmacokinetics for berberine and (−)-β-hydrastine

Data represent geometric mean [90% confidence intervals].

Fig. 6.

Observed (colored circles and triangles) and simulated (gray lines) plasma concentration-time profiles for berberine (A and C) and (−)-β-hydrastine (B and D) after oral administration of multiple doses (1 g thrice daily × 6 days) of goldenseal to 16 healthy adult participants (A and B) or a single dose of goldenseal (3 g) to eight healthy adult participants (C and D). Plasma collection for the chronic goldenseal exposure arm (A and B) began after the first dose of goldenseal was administered on day 6; goldenseal administrations are indicated with orange arrows. Each unique symbol-color pair represents a single participant, whereas green-shaded regions denote fifth and 95th percentiles for model simulations. Each gram of the goldenseal product contained 29.8 mg berberine and 25.4 mg (−)-β-hydrastine.

Simulations using the midazolam PBPK model in conjunction with the berberine and (−)-β-hydrastine models resulted in predicted pharmacokinetic interactions that were within 1.2-fold of those observed after chronic (Fig. 7A) and acute (Fig. 7B) goldenseal administration (Table 3). Simulations of active CYP3A enzyme after chronic (Fig. 7C) and acute (Fig. 7D) goldenseal administration showed that goldenseal primarily inhibited CYP3A in the small intestine, with maximum inhibition occurring at ∼2 hours after administration of the precipitant. In addition, modeling and simulations indicated that the drug interaction risk (as determined by midazolam AUC ratio >1.25) is sustained for ∼24 hours after goldenseal exposure (Fig. 7, E and F). Simulations with intravenously administered midazolam resulted in no drug interactions (Supplemental Table 5); similarly, no interactions were predicted when either the (−)-β-hydrastine model or time-dependent inhibition component was removed from the simulations.

Fig. 7.

Observed (blue dots) and simulated (gray lines) concentration-time profiles for an orally administered dose of midazolam (2.5 mg) after chronic (A) and acute (B) goldenseal exposure. Each observed data point represents the geometric mean of 16 and 8 participants for the chronic and acute goldenseal exposure groups, respectively. Green-shaded region denotes fifth and 95th percentiles for model simulations, whereas blue error bars denote the 90% confidence intervals for the observed data. Simulated time course of active CYP3A in the small intestine (SI) and liver before and after chronic (C) and acute (D) goldenseal (GS) exposure. Simulated drug interaction risk of the goldenseal-midazolam interaction with respect to time of midazolam administration after chronic (E) and acute (F) goldenseal exposure. Chronic exposure indicates multiple doses (1 g thrice daily × 6 days) of goldenseal, whereas acute exposure indicates a single dose (3 g) of goldenseal. Red dashes (E and F) indicate AUC ratio of 1.25, which represents the threshold for clinically relevant pharmacokinetic drug interactions.

TABLE 3

Predicted and observed changes in oral midazolam (2.5 mg) pharmacokinetics after chronic (1 g thrice daily for 6 days) or acute (a single 3-g dose) goldenseal administration

Data represent geometric means [90% confidence intervals]. Ratios denote the AUC0–12h or Cmax of midazolam in the presence to absence of goldenseal.

Discussion

Natural product-drug interactions remain a longstanding public health concern as more patients seek natural alternatives to supplement their pharmacotherapeutic regimens (Johnson et al., 2018; Safari et al., 2022). Goldenseal is one natural product that has been shown to inhibit CYP3A activity in humans. However, the major precipitant(s), mode of inhibition, and site of the interaction remain largely unknown. Three aspects of the pharmacokinetic goldenseal-midazolam interaction were evaluated: berberine versus (−)-β-hydrastine, reversible inhibition versus TDI, and gut versus liver. Physicochemical, CYP3A inhibition kinetic, and clinical pharmacokinetic data were used to develop a robust PBPK interaction model. Results from this integrated approach indicated that (−)-β-hydrastine, TDI of CYP3A, and the intestine are the predominant precipitant, mechanism, and site of the interaction.

The abundant goldenseal alkaloids berberine and (−)-β-hydrastine were selected as potential precipitants of the pharmacokinetic goldenseal-midazolam interaction because both contain a methylenedioxyphenyl functional group (Fig. 4), a structural moiety associated with TDI (Kalgutkar et al., 2007). The CYP3A inhibition kinetics for both alkaloids were previously established using HLMs (McDonald et al., 2020). Because the parameters obtained from HLMs can vary from those obtained from HIMs (Lalovic et al., 2004), additional in vitro experiments were conducted. Regarding reversible inhibition, effects of a given alkaloid were consistent between HIMs and HLMs. That is, the IC50 for berberine was unattainable from both enzyme sources (due to solubility limitations when >500 μM), and the IC50 for (−)-β-hydrastine was comparable (46 vs. 58 μM in HIMs and HLMs, respectively). However, regarding TDI, the inactivation efficiency (kinact/KI) for berberine in HIMs was estimated to be 5× lower than that in HLMs (∼0.24 vs. 1.3 ml/min/μmol), whereas the kinact/KI for (−)-β-hydrastine in HIMs was moderately higher than that in HLMs (4.8 vs. 2.0 ml/min/μmol) (McDonald et al., 2020). These discrepancies may be attributed to differences in CYP3A4/5 structure, modulation, and/or expression between HIMs and HLMs (Lin et al., 2002), as well as differences in the protein milieu. Regardless, comparing the recovered kinact/KI values with clinically relevant time-dependent inhibitors (diltiazem and erythromycin: 2.7 and 3.6 ml/min/μmol, respectively) (Obach et al., 2007), these in vitro observations suggest that (−)-β-hydrastine is the primary alkaloid contributing to the goldenseal-midazolam interaction via TDI.

In addition to the clinical studies conducted by Gurley and colleagues, our group recently completed a pharmacokinetic goldenseal-drug interaction study that included midazolam (Nguyen et al., 2021). Results demonstrated that multiple dose administration of goldenseal for a shorter duration than previous studies led to a comparable increase in midazolam AUC (∼40%). One prior study (Gurley et al., 2008) showed a larger increase in midazolam half-life in the presence of goldenseal relative to baseline (3.1 vs. 2.0 hours) compared with our study (4.8 vs. 3.9 hours). This discrepancy is likely due to differences in blood collection intervals (6 vs. 12 hours), with the longer sampling duration providing a more robust estimate of half-life. The current work further showed that a single dose of goldenseal increased midazolam AUC and Cmax but had no effect on half-life (Table 1). Collectively, results suggest that goldenseal inhibited CYP3A primarily in the gut, consistent with a previous report involving weak CYP3A inhibitors (exposure-to-baseline AUC ratio <2) (Yamada et al., 2020).

To gain additional mechanistic insights into the goldenseal-midazolam interaction, a PBPK interaction model was developed using a middle-out approach, incorporating both in vitro and clinical pharmacokinetic data. PBPK models were developed for berberine and (−)-β-hydrastine independently, which were then used in combination with midazolam to simulate systemic exposure to the two major goldenseal alkaloids and to predict the goldenseal-midazolam interaction. Using a learn-and-confirm method, the first-order absorption model was selected for berberine because absorption is permeability-limited. In contrast, the ADAM model was selected for (−)-β-hydrastine because absorption may be dissolution rate-limited based on its relatively high Papp. The resulting PBPK models successfully simulated exposure (AUC0–12.5h and Cmax) to berberine and (−)-β-hydrastine to within 2-fold of that observed after chronic and acute goldenseal administration (Table 2). However, AUCinf for berberine after chronic goldenseal exposure was underpredicted due to lack of robust estimates of terminal half-life (Fig. 6A). Because the inhibition kinetics of each alkaloid varied between HIMs and HLMs, the KI and kinact values correlating with the highest recovered kinact/KI and the lowest reversible Ki values were applied to simulate worst-case interaction scenarios (i.e., highest midazolam AUC ratio). The predicted interactions (AUC0–12h and Cmax ratios) were within 1.2-fold of those observed after chronic and acute goldenseal administration (Table 3), indicating successful model predictions.

Modeling and simulations revealed that active intestinal and hepatic CYP3A decreased by 60%–80% and 3%–10%, respectively, in the presence of goldenseal (Fig. 7, C and D). Simulation of the interaction using intravenous midazolam resulted in no change in AUC (Supplemental Table 5), supporting the inference that goldenseal inhibits CYP3A primarily in the gut. This observation aligns with results from a previous clinical study in which goldenseal had no effect on the pharmacokinetics of indinavir, a drug primarily metabol