Use of Physiologically Based Pharmacokinetic Modeling to Predict Human Gut Microbial Conversion of Daidzein to S-Equol.

A physiologically based pharmacokinetic (PBPK) model was developed for daidzein and its metabolite S-equol. Anaerobic in vitro incubations of pooled fecal samples from S-equol producers and nonproducers allowed definition of the kinetic constants. PBPK model-based predictions for the maximum daidzein plasma concentration (Cmax) were comparable to literature data. The predictions also revealed that the Cmax of S-equol in producers was only up to 0.22% that of daidzein, indicating that despite its higher estrogenicity, S-equol is likely to contribute to the overall estrogenicity upon human daidzein exposure to a only limited extent. An interspecies comparison between humans and rats revealed that the catalytic efficiency for S-equol formation in rats was 210-fold higher than that of human S-equol producers. The described in vitro-in silico strategy provides a proof-of-principle on how to include microbial metabolism in humans in PBPK modeling as part of the development of new approach methodologies (NAMs).

Use of Physiologically Based Kinetic Modeling to Predict Rat Gut Microbial Metabolism of the Isoflavone Daidzein to S-Equol and Its Consequences for ERα Activation.

SCOPE: To predict gut microbial metabolism of xenobiotics and the resulting plasma concentrations of metabolites formed, an in vitro-in silico-based testing strategy is developed using the isoflavone daidzein and its gut microbial metabolite S-equol as model compounds. METHODS AND RESULTS: Anaerobic rat fecal incubations are optimized and performed to derive the apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for gut microbial conversion of daidzein to dihydrodaidzein, S-equol, and O-desmethylangolensin, which are input as parameters for a physiologically based kinetic (PBK) model. The inclusion of gut microbiota in the PBK model allows prediction of S-equol concentrations and slightly reduced predicted maximal daidzein concentrations from 2.19 to 2.16 µm. The resulting predicted concentrations of daidzein and S-equol are comparable to in vivo concentrations reported. CONCLUSION: The optimized in vitro approach to quantify kinetics for gut microbial conversions, and the newly developed PBK model for rats that includes gut microbial metabolism, provide a unique tool to predict the in vivo consequences of daidzein microbial metabolism for systemic exposure of the host to daidzein and its metabolite S-equol. The predictions reveal a dominant role for daidzein in ERα-mediated estrogenicity despite the higher estrogenic potency of its microbial metabolite S-equol.

Physiologically Based Kinetic Modeling-Facilitated Reverse Dosimetry to Predict In Vivo Red Blood Cell Acetylcholinesterase Inhibition Following Exposure to Chlorpyrifos in the Caucasian and Chinese Population.

Organophosphates have a long history of use as insecticides over the world. The aim of the present study was to investigate the interethnic differences in kinetics, biomarker formation, and in vivo red blood cell acetylcholinesterase inhibition of chlorpyrifos (CPF) in the Chinese and the Caucasian population. To this purpose, physiologically based kinetic models for CPF in both the Chinese and Caucasian population were developed, and used to study time- and dose-dependent interethnic variation in urinary biomarkers and to convert concentration-response curves for red blood cell acetylcholinesterase inhibition to in vivo dose-response curves in these 2 populations by reverse dosimetry. The results obtained revealed a marked interethnic difference in toxicokinetics of CPF, with lower urinary biomarker levels at similar dose levels and slower CPF bioactivation and faster chlorpyrifos-oxon detoxification in the Chinese compared with the Caucasian population, resulting in 5- to 6-fold higher CPF sensitivity of the Caucasian than the Chinese population. These differences might be related to variation in the frequency of single-nucleotide polymorphisms for the major biotransformation enzymes involved. To conclude, the interethnic variation in kinetics of CPF may affect both its biomarker-based exposure assessment and its toxicity and risk assessment and physiologically based kinetic modeling facilitates the characterization and quantification of these interethnic variations.

In vitro-in silico-based analysis of the dose-dependent in vivo oestrogenicity of the soy phytoestrogen genistein in humans.

BACKGROUND AND PURPOSE: The in vivo oestrogenicity of genistein and its glycoside genistin is still under debate. The present study aimed to develop a physiologically based kinetic (PBK) model that provides insight in dose-dependent plasma concentrations of genistein aglycone and its metabolites and enables prediction of in vivo oestrogenic effective dose levels of genistein and genistin in humans. EXPERIMENTAL APPROACH: A PBK model for genistein and genistin in humans was developed based on in vitro metabolic parameters. The model obtained was used to translate in vitro oestrogenic concentration-response curves of genistein to in vivo oestrogenic dose-response curves for intake of genistein and genistin. KEY RESULTS: The model predicted that genistein-7-O-glucuronide was the major circulating metabolite and that levels of the free aglycone were generally low [0.5-17% of total plasma genistein at oral doses from 0.01 to 50 mg (kg·bw)-1 ]. The predicted in vivo benchmark dose for 5% response values for oestrogenicity varied between 0.06 and 4.39 mg kg-1 genistein. For genistin, these values were 1.3-fold higher. These values are in line with reported human data and show that oestrogenic responses can be expected at an Asian dietary and a supplementary intake, while intake resulting from a Western diet may not be effective. CONCLUSIONS AND IMPLICATIONS: The present study shows how plasma concentrations of genistein and its metabolites and oestrogenic dose levels of genistein in humans can be predicted by combining in vitro oestrogenicity with PBK model-based reverse dosimetry, eliminating the need for human intervention studies.

Physiologically based kinetic modeling of hesperidin metabolism and its use to predict in vivo effective doses in humans.

SCOPE: To develop a physiologically based kinetic (PBK) model that describes the absorption, distribution, metabolism, and excretion of hesperidin in humans, enabling the translation of in vitro concentration-response curves to in vivo dose-response curves. METHODS AND RESULTS: The PBK model for hesperidin in humans was developed based on in vitro metabolic parameters. Hesperidin was predicted to mainly occur in the systemic circulation as different monoglucuronides. The plasma concentrations of hesperidin aglycone (hesperetin) was predicted to be <0.02 mg/L at an oral dose of 50 mg/kg bw. The developed PBK model allowed conversion of in vitro concentration-response curves for different effects to in vivo dose-response curves. The BMD05 (benchmark dose for 5% response) values for protein kinase A inhibition ranged between 135 and 529 mg/kg bw hesperidin, and for inhibition of endothelial cell migration and prostaglandin E2 and nitric oxide production ranged between 2.19 and 44 mg/kg bw hesperidin. These values are in line with reported human data showing in vivo effects by hesperidin and show that these effects may occur at Western dietary and supplementary intake of hesperidin. CONCLUSIONS: The developed PBK model adequately predicts absorption, distribution, metabolism, and excretion of hesperidin in humans and allows to evaluate the human in vivo situation without the need for human intervention studies.

Use of physiologically based kinetic (PBK) modeling to study interindividual human variation and species differences in plasma concentrations of quercetin and its metabolites.

Biological activities of flavonoids in vivo ultimately depend on the systemic bioavailability of the aglycones and their metabolites. We aimed to develop physiologically based kinetic (PBK) models to predict plasma concentrations of the flavonoid quercetin and its metabolites in individual human subjects and to define species differences compared with male rat. The human models were developed based on in vitro metabolic parameters derived from incubations with pooled and 20 individual human tissue fractions and by fitting kinetic parameters to available in vivo data. The outcomes obtained were compared to a previously developed model for quercetin and its metabolites formation in male rat. Quercetin-3'-O-glucuronide was predicted to be the major circulating metabolite in 19 out of 20 individuals, while in male rat di- and tri-conjugates of quercetin containing a glucuronic acid, sulfate and/or methyl moieties are the major metabolites. Significant species differences occur in major circulating metabolites of quercetin suggesting that rat is not an adequate model to study effects of quercetin in man. The defined PBK models can be used to guide the experimental design of in vitro experiments with flavonoids, especially to better take into account the relevance of metabolism and the contribution of metabolites to the biological activity in humans.

A physiologically based kinetic (PBK) model describing plasma concentrations of quercetin and its metabolites in rats.

Biological activities of flavonoids in vivo are ultimately dependent on the systemic bioavailability of the aglycones as well as their metabolites. In the present study, a physiologically based kinetic (PBK) model was developed to predict plasma concentrations of the flavonoid quercetin and its metabolites and to tentatively identify the regiospecificity of the major circulating metabolites. The model was developed based on in vitro metabolic parameters and by fitting kinetic parameters to literature available in vivo data. Both exposure to quercetin aglycone and to quercetin-4'-O-glucoside, for which in vivo data were available, were simulated. The predicted plasma concentrations of different metabolites adequately matched literature reported plasma concentrations of these metabolites in rats exposed to 4'-O-glucoside. The bioavailability of aglycone was predicted to be very low ranging from 0.004%-0.1% at different oral doses of quercetin or quercetin-4'-O-glucoside. Glucuronidation was a crucial pathway that limited the bioavailability of the aglycone, with 95-99% of the dose being converted to monoglucuronides within 1.5-2.5h at different dose levels ranging from 0.1 to 50mg/kg bw quercetin or quercetin-4'-O-glucoside. The fast metabolic conversion to monoglucuronides allowed these metabolites to further conjugate to di- and tri-conjugates. The regiospecificity of major circulating metabolites was observed to be dose-dependent. As we still lack in vivo kinetic data for many flavonoids, the developed model has a great potential to be used as a platform to build PBK models for other flavonoids as well as to predict the kinetics of flavonoids in humans.

Malabaricone C-containing mace extract inhibits safrole bioactivation and DNA adduct formation both in vitro and in vivo.

Safrole, present in mace and its essential oils, causes liver tumors in rodents at high dose levels due to formation of a DNA reactive 1'-sulfooxysafrole. The present study identifies malabaricone C as a mace constituent able to inhibit safrole DNA adduct formation at the level of sulfotransferase mediated bioactivation. This inhibition was incorporated into physiologically based biokinetic rat and human models. Dosing safrole at 50mg/kg body weight and malabaricone C-containing mace extract at a ratio reflecting the relative presence in mace, and assuming 100% or 1% uptake of malabaricone C-containing mace extract, the model predicted inhibition of 1'-sulfooxysafrole formation for rats and humans by 90% and 100% or 61% and 91%, respectively. To validate the model, mace extract and safrole were co-administered orally to Sprague-Dawley rats. LC-ECI-MS/MS based quantification of DNA adduct levels revealed a significant (p<0.01) 55% reduction of safrole DNA adduct formation by malabaricone C-containing mace extract in the liver of rats exposed to safrole. The data obtained were used to perform a refined risk assessment of safrole. Overall, the results suggest a lower tumor incidence when safrole would be tested within a relevant food matrix containing sulfotransferase inhibitors compared to dosing pure safrole.