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).

Physiologically based kinetic modelling based prediction of in vivo rat and human acetylcholinesterase (AChE) inhibition upon exposure to diazinon.

The present study predicts in vivo human and rat red blood cell (RBC) acetylcholinesterase (AChE) inhibition upon diazinon (DZN) exposure using physiological based kinetic (PBK) modelling-facilitated reverse dosimetry. Due to the fact that both DZN and its oxon metabolite diazoxon (DZO) can inhibit AChE, a toxic equivalency factor (TEF) was included in the PBK model to combine the effect of DZN and DZO when predicting in vivo AChE inhibition. The PBK models were defined based on kinetic constants derived from in vitro incubations with liver fractions or plasma of rat and human, and were used to translate in vitro concentration-response curves for AChE inhibition obtained in the current study to predicted in vivo dose-response curves. The predicted dose-response curves for rat matched available in vivo data on AChE inhibition, and the benchmark dose lower confidence limits for 10% inhibition (BMDL10 values) were in line with the reported BMDL10 values. Humans were predicted to be 6-fold more sensitive than rats in terms of AChE inhibition, mainly because of inter-species differences in toxicokinetics. It is concluded that the TEF-coded DZN PBK model combined with quantitative in vitro to in vivo extrapolation (QIVIVE) provides an adequate approach to predict RBC AChE inhibition upon acute oral DZN exposure, and can provide an alternative testing strategy for derivation of a point of departure (POD) in risk assessment.

Defining in vivo dose-response curves for kidney DNA adduct formation of aristolochic acid I in rat, mouse and human by an in vitro and physiologically based kinetic modeling approach.

Aristolochic acid I (AAI) is a well-known genotoxic kidney carcinogen. Metabolic conversion of AAI into the DNA-reactive aristolactam-nitrenium ion is involved in the mode of action of tumor formation. This study aims to predict in vivo AAI-DNA adduct formation in the kidney of rat, mouse and human by translating the in vitro concentration-response curves for AAI-DNA adduct formation to the in vivo situation using physiologically based kinetic (PBK) modeling-based reverse dosimetry. DNA adduct formation in kidney proximal tubular LLC-PK1 cells exposed to AAI was quantified by liquid chromatography-electrospray ionization-tandem mass spectrometry. Subsequently, the in vitro concentration-response curves were converted to predicted in vivo dose-response curves in rat, mouse and human kidney using PBK models. Results obtained revealed a dose-dependent increase in AAI-DNA adduct formation in the rat, mouse and human kidney and the predicted DNA adduct levels were generally within an order of magnitude compared with values reported in the literature. It is concluded that the combined in vitro PBK modeling approach provides a novel way to define in vivo dose-response curves for kidney DNA adduct formation in rat, mouse and human and contributes to the reduction, refinement and replacement of animal testing.

Induction of peroxisome proliferator activated receptor γ (PPARγ) mediated gene expression and inhibition of induced nitric oxide production by Maerua subcordata (Gilg) DeWolf.

BACKGROUND: The health benefits of botanicals is linked to their phytochemicals that often exert pleiotropic effects via targeting multiple molecular signaling pathways such as the peroxisome proliferator-activated receptors (PPARs) and the nuclear factor kappaB (NFκB). The PPARs are transcription factors that control metabolic homeostasis and inflammation while the NF-κB is a master regulator of inflammatory genes such as the inducible nitric-oxide synthase that result in nitric oxide (NO) overproduction. METHODS: Extracts of Maerua subcordata (MS) and selected candidate constituents thereof, identified by liquid chromatography coupled to mass spectroscopy, were tested for their ability to induce PPARγ mediated gene expression in U2OS-PPARγ cells using luciferase reporter gene assay and also for their ability to inhibit lipopolysaccharide (LPS) induced NO production in RAW264.7 macrophages. While measuring the effect of test samples on PPARγ mediated gene expression, a counter assay that used U2OS-Cytotox cells was performed to monitor cytotoxicity or any non-specific changes in luciferase activity. RESULTS: The results revealed that the fruit, root, and seed extracts were non-cytotoxic up to a concentration of 30 g dry weight per litre (gDW/L) and induced PPARγ mediated gene expression but the leaf extract showed some cytotoxicity and exhibited minimal induction. Instead, all extracts showed concentration (1-15 gDW/L) dependent inhibition of LPS induced NO production. The root extract showed weaker inhibition. Among the candidate constituents, agmatine, stachydrine, trigonelline, indole-3-carboxyaldehyde, plus ethyl-, isobutyl-, isopropyl, and methyl-isothiocyanates showed similar inhibition, and most showed increased inhibition with increasing concentration (1-100 µM) although to a lesser potency than the positive control, aminoguanidine. CONCLUSION: The present study demonstrated for the first time the induction of PPARγ mediated gene expression by MS fruit, root, and seed extracts and the inhibition of LPS induced NO production by MS fruit, leaf, root, and seed extracts and some candidate constituents thereof.

Microfluidic chip for culturing intestinal epithelial cell layers: Characterization and comparison of drug transport between dynamic and static models.

Dynamic flow in vitro models are currently widely explored for their applicability in drug development research. The application of gut-on-chip models in toxicology is lagging behind. Here we report the application of a gut-on-chip model for biokinetic studies and compare the observed biokinetics of reference compounds with those obtained using a conventional static in vitro model. Intestinal epithelial Caco-2 cells were cultured on a porous membrane assembled between two glass flow chambers for the dynamic model, or on a porous membrane in a Transwell model. Confocal microscopy, lucifer yellow translocation, and alkaline phosphatase activity evaluation revealed that cells cultured in the gut-on-chip model formed tight, differentiated, polarized monolayers like in the static cultures. In the dynamic gut-on-chip model the transport of the high permeability compounds antipyrine, ketoprofen and digoxin was lower (i.e. 4.2-, 2.7- and 1.9-fold respectively) compared to the transport in the static Transwell model. The transport of the low permeability compound, amoxicillin, was similar in both the dynamic and static in vitro model. The obtained transport values of the compounds are in line with the compound Biopharmaceuticals Classification System. It is concluded that the gut-on-chip provides an adequate model for transport studies of chemicals.

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.

The role of metabolism in the developmental toxicity of polycyclic aromatic hydrocarbon-containing extracts of petroleum substances.

In vitro assays presently used for prenatal developmental toxicity (PDT) testing only assess the embryotoxic potential of parent substances and not that of potentially embryotoxic metabolites. Here we combined a biotransformation system, using hamster liver microsomes, with the ES-D3 cell differentiation assay of the embryonic stem cell test (EST) to compare the in vitro PDT potency of two 5-ring polycyclic aromatic hydrocarbons (PAHs), benzo[a]pyrene (BaP) and dibenz[a,h]anthracene (DBA), and dimethyl sulfoxide extracts from five PAH-containing petroleum substances (PS) and a gas-to-liquid base oil (GTLb), with and without bioactivation. In the absence of bioactivation, DBA, but not BaP, inhibited the differentiation of ES-D3 cells into beating cardiomyocytes in a concentration-dependent manner. Upon bioactivation, BaP induced in vitro PDT, while its major metabolite 3-hydroxybenzo[a]pyrene was shown to be active in the EST as well. This means BaP needs biotransformation to exert its embryotoxic effects. GTLb extracts tested negative in the EST, with and without bioactivation. The PS-induced PDT in the EST was not substantially changed following bioactivation, implying that metabolism may not play a crucial role for the PS extracts under study to exert the in vitro PDT effects. Altogether, these results indicate that although some PAH require bioactivation to induce PDT, some do not and this latter appears to hold for the (majority of) the PS constituents responsible for the in vitro PDT of these complex substances.

Hazard assessment of Maerua subcordata (Gilg) DeWolf. for selected endpoints using a battery of in vitro tests.

ETHNOPHARMACOLOGICAL RELEVANCE: Maerua subcordata (Gilg) DeWolf is a medicinal and wild food plant growing mainly in east Africa. Especially its root tuber is widely used in traditional medicine to treat several infectious and chronic diseases but also in some toxicity implications like use as abortifacient. AIM OF THE STUDY: the present study applied in silico and in vitro tests to identify possible hazards of M. subcordata (fruit, leaf, root, seed) methanol extracts focussing on developmental toxicity. MATERIALS AND METHODS: Ames test, estrogen receptor alpha (ERα) assay, aryl hydrocarbon receptor (AhR) assay, embryonic stem cell test (EST), and zebrafish embryotoxicity test (ZET) were employed. Besides, a Derek Nexus toxicity prediction was performed on candidate structures obtained from metabolomics profiling of the extracts using liquid chromatography coupled to multistage mass spectroscopy (LC/MSn) and a MAGMa software based structural annotation. RESULTS: Glucosinolates, which degrade to isothiocyanates, and biogenic amines were among the candidate molecules identified in the extracts by LC/MSn - MAGMa software structural annotation. Isothiocyanates and some other candidate molecules suggested a positive mutagenicity alert in Derek toxicity predictions. All the extracts showed negative mutagenicity in the Ames test. However, the Derek predictions also identified endocrine and developmental toxicity as possible endpoints of concern. This was further assessed using in vitro tests. Results obtained reveal that leaf extract shows AhR and ERα agonist activities, inhibited differentiation of ES-D3 stem cells into contracting cardiomyocytes in the EST (p < 0.001) as well as inhibited hatching (p < 0.01) and showed acute toxicity (p < 0.01) in the ZET. Also, the fruit extract showed toxicity (p < 0.05) towards zebrafish embryos and both fruit and seed extracts showed AhR agonist activities while root extract was devoid of activity in all in vitro assays. CONCLUSION: The leaf extract tests positive in in vitro tests that may point towards a developmental toxicity hazard. The current evaluations did not raise concerns of genotoxicity or developmental toxicity for the fruit, seed and root extracts. This is important given the use of especially these parts of M. subcordata, in traditional medicine and/or as (famine) food.

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.

Effects of Maerua subcordata (Gilg) DeWolf on electrophile-responsive element (EpRE)-mediated gene expression in vitro.

Plant extracts and phytochemicals may prevent chronic diseases via activation of adaptive cellular stress response pathways including induction of antioxidant and phase II detoxifying enzymes. The regulatory regions of these inducible genes encode the electrophile-response element (EpRE). This study tested the EpRE induction ability of Maerua subcordata (fruit, leaf, root, seed) methanol extracts and selected candidate constituents thereof, identified by liquid chromatography coupled with multistage mass spectroscopy, employing an EpRE luciferase reporter gene assay using hepa-1c1c7 mouse hepatoma cells. A parallel Cytotox CALUX assay using human osteosarcoma U2OS cells was used to monitor any non-specific changes in luciferase activity or cytotoxicity. Results showed that fruit, root, and seed extracts were non-cytotoxic up to a concentration of 30 gram dry weight per litre but the leaf extract exhibited some cytotoxicity and that the leaf (despite some cytotoxicity), fruit, and seed extracts showed strong induction of EpRE mediated gene expression while induction by the root extract was minimal. Selected candidates included glucosinolates, isothiocyanates, and some biogenic amines. Subsequent studies showed that methyl-, ethyl-, isopropyl-, isobutyl- isothiocyanates, and sec-butyl thiocyanate as well as glucobrassicin induced concentration (1-100 µM) dependent EpRE-mediated gene expression while the biogenic amines stachydrine and trigonelline acted as inhibitors of EpRE-mediated gene expression at 100 µM. The identification of glucolepidiin, glucobrassicin, glucocapparin, stachydrine, and trigonelline in all extracts was confirmed using standards and based on multiple reaction monitoring; yet, glucobrassicin level in the root extract was negligible. In conclusion, this study provided a first report on EpRE mediated gene expression effects of M. subcordata; and despite detection of different glucosinolates in all extracts, those containing glucobrassicin particularly displayed high EpRE induction. Because EpRE inducers are cytoprotective and potential chemopreventive agents while inhibitors are suggested adjuvants of chemotherapy, results of this study imply that process manipulation of this plant may result in herbal preparations that may be used as chemopreventive agents or adjuvants of chemotherapies.

Development of a Combined In Vitro Physiologically Based Kinetic (PBK) and Monte Carlo Modelling Approach to Predict Interindividual Human Variation in Phenol-Induced Developmental Toxicity.

With our recently developed in vitro physiologically based kinetic (PBK) modelling approach, we could extrapolate in vitro toxicity data to human toxicity values applying PBK-based reverse dosimetry. Ideally information on kinetic differences among human individuals within a population should be considered. In the present study, we demonstrated a modelling approach that integrated in vitro toxicity data, PBK modelling and Monte Carlo simulations to obtain insight in interindividual human kinetic variation and derive chemical specific adjustment factors (CSAFs) for phenol-induced developmental toxicity. The present study revealed that UGT1A6 is the primary enzyme responsible for the glucuronidation of phenol in humans followed by UGT1A9. Monte Carlo simulations were performed taking into account interindividual variation in glucuronidation by these specific UGTs and in the oral absorption coefficient. Linking Monte Carlo simulations with PBK modelling, population variability in the maximum plasma concentration of phenol for the human population could be predicted. This approach provided a CSAF for interindividual variation of 2.0 which covers the 99th percentile of the population, which is lower than the default safety factor of 3.16 for interindividual human kinetic differences. Dividing the dose-response curve data obtained with in vitro PBK-based reverse dosimetry, with the CSAF provided a dose-response curve that reflects the consequences of the interindividual variability in phenol kinetics for the developmental toxicity of phenol. The strength of the presented approach is that it provides insight in the effect of interindividual variation in kinetics for phenol-induced developmental toxicity, based on only in vitro and in silico testing.

Study on inter-ethnic human differences in bioactivation and detoxification of estragole using physiologically based kinetic modeling.

Considering the rapid developments in food safety in the past decade in China, it is of importance to obtain insight into what extent safety and risk assessments of chemicals performed for the Caucasian population apply to the Chinese population. The aim of the present study was to determine physiologically based kinetic (PBK) modeling-based predictions for differences between Chinese and Caucasians in terms of metabolic bioactivation and detoxification of the food-borne genotoxic carcinogen estragole. The PBK models were defined based on kinetic constants for hepatic metabolism derived from in vitro incubations using liver fractions of the two ethnic groups, and used to evaluate the inter-ethnic differences in metabolic activation and detoxification of estragole. The models predicted that at realistic dietary intake levels, only 0.02% of the dose was converted to the ultimate carcinogenic metabolite 1'-sulfooxyestragole in Chinese subjects, whereas this amounted to 0.09% of the dose in Caucasian subjects. Detoxification of 1'-hydroxyestragole, mainly via conversion to 1'-oxoestragole, was similar within the two ethnic groups. The 4.5-fold variation in formation of the ultimate carcinogenic metabolite of estragole accompanied by similar rates of detoxification may indicate a lower risk of estragole for the Chinese population at similar levels of exposure. The study provides a proof of principle for how PBK modeling can identify differences in ethnic sensitivity and provide a more refined risk assessment for a specific ethnic group for a compound of concern.

Integrating in vitro data and physiologically based kinetic (PBK) modelling to assess the in vivo potential developmental toxicity of a series of phenols.

Toxicity outcomes derived in vitro do not always reflect in vivo toxicity values, which was previously observed for a series of phenols tested in the embryonic stem cell test (EST). Translation of in vitro data to the in vivo situation is therefore an important, but still limiting step for the use of in vitro toxicity outcomes in the safety assessment of chemicals. The aim of the present study was to translate in vitro embryotoxicity data for a series of phenols to in vivo developmental toxic potency values for the rat by physiologically based kinetic (PBK) modelling-based reverse dosimetry. To this purpose, PBK models were developed for each of the phenols. The models were parameterised with in vitro-derived values defining metabolism and transport of the compounds across the intestinal and placental barrier and with in silico predictions and data from the literature. Using PBK-based reverse dosimetry, in vitro concentration-response curves from the EST were translated into in vivo dose-response curves from which points of departure (PoDs) were derived. The predicted PoDs differed less than 3.6-fold from PoDs derived from in vivo toxicity data for the phenols available in the literature. Moreover, the in vitro PBK-based reverse dosimetry approach could overcome the large disparity that was observed previously between the in vitro and the in vivo relative potency of the series of phenols. In conclusion, this study shows another proof-of-principle that the in vitro PBK approach is a promising strategy for non-animal-based safety assessment of chemicals.

Evaluation of the interindividual human variation in bioactivation of methyleugenol using physiologically based kinetic modeling and Monte Carlo simulations.

The present study aims at predicting the level of formation of the ultimate carcinogenic metabolite of methyleugenol, 1'-sulfooxymethyleugenol, in the human population by taking variability in key bioactivation and detoxification reactions into account using Monte Carlo simulations. Depending on the metabolic route, variation was simulated based on kinetic constants obtained from incubations with a range of individual human liver fractions or by combining kinetic constants obtained for specific isoenzymes with literature reported human variation in the activity of these enzymes. The results of the study indicate that formation of 1'-sulfooxymethyleugenol is predominantly affected by variation in i) P450 1A2-catalyzed bioactivation of methyleugenol to 1'-hydroxymethyleugenol, ii) P450 2B6-catalyzed epoxidation of methyleugenol, iii) the apparent kinetic constants for oxidation of 1'-hydroxymethyleugenol, and iv) the apparent kinetic constants for sulfation of 1'-hydroxymethyleugenol. Based on the Monte Carlo simulations a so-called chemical-specific adjustment factor (CSAF) for intraspecies variation could be derived by dividing different percentiles by the 50th percentile of the predicted population distribution for 1'-sulfooxymethyleugenol formation. The obtained CSAF value at the 90th percentile was 3.2, indicating that the default uncertainty factor of 3.16 for human variability in kinetics may adequately cover the variation within 90% of the population. Covering 99% of the population requires a larger uncertainty factor of 6.4. In conclusion, the results showed that adequate predictions on interindividual human variation can be made with Monte Carlo-based PBK modeling. For methyleugenol this variation was observed to be in line with the default variation generally assumed in risk assessment.

Conversion of major soy isoflavone glucosides and aglycones in in vitro intestinal models.

SCOPE: This study compares conversion of three major soy isoflavone glucosides and their aglycones in a series of in vitro intestinal models. METHODS AND RESULTS: In an in vitro human digestion model isoflavone glucosides were not deconjugated, whereas studies in a Caco-2 transwell model confirmed that deconjugation is essential to facilitate transport across the intestinal barrier. Deconjugation was shown upon incubation of the isoflavone glucosides with rat as well as human intestinal S9. In incubations with rat intestinal S9 lactase phlorizin hydrolase, glucocerebrosidase, and cytosolic broad-specific ß-glucosidase all contribute significantly to deconjugation, whereas in incubations with human intestinal S9 deconjugation appeared to occur mainly through the activity of broad-specific ß-glucosidase. Species differences in glucuronidation and sulfation were limited and generally within an order of magnitude with 7-O-glucuronides being the major metabolites for all three isoflavone aglycones and the glucuronidation during first pass metabolism being more efficient in rats than in humans. Comparison of the catalytic efficiencies reveals that deconjugation is less efficient than conjugation confirming that aglycones are unlikely to enter the systemic circulation. CONCLUSION: Altogether, the data point at possible differences in the characteristics for intestinal conversion of the major soy isoflavones between rat and human, especially with respect to their deconjugation.

Combining in vitro embryotoxicity data with physiologically based kinetic (PBK) modelling to define in vivo dose-response curves for developmental toxicity of phenol in rat and human.

In vitro assays are often used for the hazard characterisation of compounds, but their application for quantitative risk assessment purposes is limited. This is because in vitro assays cannot provide a complete in vivo dose-response curve from which a point of departure (PoD) for risk assessment can be derived, like the no observed adverse effect level (NOAEL) or the 95 % lower confidence limit of the benchmark dose (BMDL). To overcome this constraint, the present study combined in vitro data with a physiologically based kinetic (PBK) model applying reverse dosimetry. To this end, embryotoxicity of phenol was evaluated in vitro using the embryonic stem cell test (EST), revealing a concentration-dependent inhibition of differentiation into beating cardiomyocytes. In addition, a PBK model was developed on the basis of in vitro and in silico data and data available from the literature only. After evaluating the PBK model performance, effective concentrations (ECx) obtained with the EST served as an input for in vivo plasma concentrations in the PBK model. Applying PBK-based reverse dosimetry provided in vivo external effective dose levels (EDx) from which an in vivo dose-response curve and a PoD for risk assessment were derived. The predicted PoD lies within the variation of the NOAELs obtained from in vivo developmental toxicity data from the literature. In conclusion, the present study showed that it was possible to accurately predict a PoD for the risk assessment of phenol using in vitro toxicity data combined with reverse PBK modelling.

Physiologically based kinetic models for the alkenylbenzene elemicin in rat and human and possible implications for risk assessment.

The present study describes physiologically based kinetic (PBK) models for the alkenylbenzene elemicin (3,4,5-trimethoxyallylbenzene) in rat and human, based on the PBK models previously developed for the structurally related alkenylbenzenes estragole, methyleugenol, and safrole. Using the newly developed models, the level of metabolic activation of elemicin in rat and human was predicted to obtain insight in species differences in the bioactivation of elemicin and read across to the other methoxy allylbenzenes, estragole and methyleugenol. Results reveal that the differences between rat and human in the formation of the proximate carcinogenic metabolite 1'-hydroxyelemicin and the ultimate carcinogenic metabolite 1'-sulfoxyelemicin are limited (<3.8-fold). In addition, a comparison was made between the relative importance of bioactivation for elemicin and that of estragole and methyleugenol. Model predictions indicate that compound differences in the formation of the 1'-sulfoxymetabolites are limited (<11-fold) in rat and human liver. The insights thus obtained were used to perform a risk assessment for elemicin using the margin of exposure (MOE) approach and read across to the other methoxy allylbenzene derivatives for which in vivo animal tumor data are available. This reveals that elemicin poses a lower priority for risk management as compared to its structurally related analogues estragole and methyleugenol. Altogether, the results obtained indicate that PBK modeling provides an important insight in the occurrence of species differences in the metabolic activation of elemicin. Moreover, they provide an example of how PBK modeling can facilitate a read across in risk assessment from compounds for which in vivo toxicity studies are available to a compound for which only limited toxicity data have been described, thus contributing to the development of alternatives for animal testing.

In vivo validation of DNA adduct formation by estragole in rats predicted by physiologically based biodynamic modelling.

Estragole is a naturally occurring food-borne genotoxic compound found in a variety of food sources, including spices and herbs. This results in human exposure to estragole via the regular diet. The objective of this study was to quantify the dose-dependent estragole-DNA adduct formation in rat liver and the urinary excretion of 1'-hydroxyestragole glucuronide in order to validate our recently developed physiologically based biodynamic (PBBD) model. Groups of male outbred Sprague Dawley rats (n = 10, per group) were administered estragole once by oral gavage at dose levels of 0 (vehicle control), 5, 30, 75, 150, and 300mg estragole/kg bw and sacrificed after 48h. Liver, kidney and lungs were analysed for DNA adducts by LC-MS/MS. Results obtained revealed a dose-dependent increase in DNA adduct formation in the liver. In lungs and kidneys DNA adducts were detected at lower levels than in the liver confirming the occurrence of DNA adducts preferably in the target organ, the liver. The results obtained showed that the PBBD model predictions for both urinary excretion of 1'-hydroxyestragole glucuronide and the guanosine adduct formation in the liver were comparable within less than an order of magnitude to the values actually observed in vivo. The PBBD model was refined using liver zonation to investigate whether its predictive potential could be further improved. The results obtained provide the first data set available on estragole-DNA adduct formation in rats and confirm their occurrence in metabolically active tissues, i.e. liver, lung and kidney, while the significantly higher levels found in liver are in accordance with the liver as the target organ for carcinogenicity. This opens the way towards future modelling of dose-dependent estragole liver DNA adduct formation in human.

Relative embryotoxic potency of p-substituted phenols in the embryonic stem cell test (EST) and comparison to their toxic potency in vivo and in the whole embryo culture (WEC) assay.

The applicability of the embryonic stem cell test (EST) as an alternative for in vivo embryotoxicity testing was evaluated for a series of five p-substituted phenols. To this purpose, the potency ranking for this class of compounds derived from the inhibition of cardiomyocyte differentiation in the EST was compared to in vivo embryotoxic potency data obtained from literature and to the potency ranking defined in the in vitro whole embryo culture (WEC) assay. From the results obtained it appears that the EST was able to identify the embryotoxic potential for p-substituted phenols, providing an identical potency ranking compared to the WEC assay. However, the EST was not able to predict an accurate ranking for the phenols compared to their potency observed in vivo. Only phenol, the least potent compound within this series, was correctly ranked. Furthermore, p-mercaptophenol was correctly identified as a relative potent congener of the phenols tested, but its ranking was distorted by p-heptyloxyphenol, of which the toxicity was overestimated in the EST. It is concluded that when attempting to explain the observed disparity in potency rankings between in vitro and in vivo embryotoxicity, the in vitro models should be combined with a kinetic model describing in vivo absorption, distribution, metabolism and excretion processes of the compounds.

Physiologically based kinetic modeling of bioactivation and detoxification of the alkenylbenzene methyleugenol in human as compared with rat.

This study defines a physiologically based kinetic (PBK) model for methyleugenol (ME) in human based on in vitro and in silico derived parameters. With the model obtained, bioactivation and detoxification of methyleugenol (ME) at different doses levels could be investigated. The outcomes of the current model were compared with those of a previously developed PBK model for methyleugenol (ME) in male rat. The results obtained reveal that formation of 1'-hydroxymethyleugenol glucuronide (1'HMEG), a major metabolic pathway in male rat liver, appears to represent a minor metabolic pathway in human liver whereas in human liver a significantly higher formation of 1'-oxomethyleugenol (1'OME) compared with male rat liver is observed. Furthermore, formation of 1'-sulfooxymethyleugenol (1'HMES), which readily undergoes desulfonation to a reactive carbonium ion (CA) that can form DNA or protein adducts (DA), is predicted to be the same in the liver of both human and male rat at oral doses of 0.0034 and 300 mg/kg bw. Altogether despite a significant difference in especially the metabolic pathways of the proximate carcinogenic metabolite 1'-hydroxymethyleugenol (1'HME) between human and male rat, the influence of species differences on the ultimate overall bioactivation of methyleugenol (ME) to 1'-sulfooxymethyleugenol (1'HMES) appears to be negligible. Moreover, the PBK model predicted the formation of 1'-sulfooxymethyleugenol (1'HMES) in the liver of human and rat to be linear from doses as high as the benchmark dose (BMD10) down to as low as the virtual safe dose (VSD). This study shows that kinetic data do not provide a reason to argue against linear extrapolation from the rat tumor data to the human situation.