Search

The goal of the present study was to develop an online web-based toolbox that contains generic physiologically based kinetic (PBK) models for rats and humans, including underlying calculation tools to predict plasma protein binding and tissue:plasma distribution, to be used for quantitative in-vitro-to-in-vivo extrapolations (QIVIVE). The PBK models within the toolbox allow first estimations of internal plasma and tissue concentrations of chemicals to be made, based on the logP and pKa of the chemicals and values for intestinal uptake and intrinsic hepatic clearance. As a case study, the toolbox was used to predict oral equivalent doses of in vitro ToxCast bioactivity data for the food additives methylparaben, propyl gallate, octyl gallate, and dodecyl gallate. These oral equivalent doses were subsequently compared with human exposure estimates, as a low tier assessment allowing prioritization for further assessment. The results revealed that daily intake levels of especially propyl gallate can lead to internal plasma concentrations that are close to in vitro biological effect concentrations, particularly with respect to the inhibition of human thyroid peroxidase (TPO). Estrogenic effects were not considered likely to be induced by the food additives, as daily exposure levels of the different compounds remained 2 orders of magnitude below the oral equivalent doses for in vitro estrogen receptor activation. Overall, the results of the study show how the toolbox, which is freely accessible through www.qivivetools.wur.nl, can be used to obtain initial internal dose estimates of chemicals and to prioritize chemicals for further assessment, based on the comparison of oral equivalent doses of in vitro biological activity data with human exposure levels.

In vitro models are widely used to study the biotransformation of xenobiotics and to provide input parameters to physiologically based kinetic models required to predict the kinetic behavior in vivo. For farm animals this is not common practice yet. The use of slaughterhouse-derived tissue material may provide opportunities to study biotransformation reactions in farm animals. The goal of the present study was to explore the potential of slaughterhouse-derived bovine liver S9 (S9) and precision cut liver slices (PCLSs) to capture observed biotransformation reactions of lidocaine in cows. The in vitro data obtained with both S9 and PCLSs confirm in vivo findings that 2,6-dimethylaniline (DMA) is an important metabolite of lidocaine in cows, being for both PCLSs and S9 the end-product. In case of S9, also conversion of lidocaine to lidocaine-N-oxide and monoethylglycinexylidine (MEXG) was observed. MEGX is considered as intermediate for DMA formation, given that this metabolite was metabolized to DMA by both PLCSs and S9. In contrast to in vivo, no in vitro conversion of DMA to 4-OH-DMA was observed. Further work is needed to explain this lack of conversion and to further evaluate the use of slaughterhouse-derived tissue materials to predict the biotransformation of xenobiotics in farm animals.

Background and purpose In vivo estrogenicity 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 estrogenic 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 estrogenic concentration-response curves of genistein to in vivo estrogenic 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-50 mg (kg bw)-1). The predicted in vivo BMD05 (benchmark dose for 5% response) values for estrogenicity varied between 0.06 and 4.39 mg (kg bw)-1 genistein. For genistin, these values were 1.3-fold higher. These values are in line with reported human data and show that estrogenic responses can be expected at an Asian dietary and a supplementary intake, while intake resulting from a Western diet may not be effective. Conclusion and implications The present study shows how plasma concentrations of genistein and its metabolites and estrogenic dose levels of genistein in humans can be predicted by combining in vitro estrogenicity with PBK model based reverse dosimetry, eliminating the need for human intervention studies.

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., This study describes a physiologically based kinetic (PBK) modeling‐based reverse dosimetry approach to convert in vitro concentration‐response curves for aristolochic acid I‐DNA adduct formation quantified in LLC‐PK1 cells to in vivo dose‐response curves in the kidney of rat, mouse and human. The predicted DNA adduct levels were within an order of magnitude from values reported in the literature. In conclusion, the in vitro PBK approach predicts DNA adduct formation with a similar level of uncertainty, as observed in the experimental animal studies.

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. Electronic supplementary material The online version of this article (doi:10.1007/s00204-016-1881-x) contains supplementary material, which is available to authorized users.

The goal of the present study was to obtain an in vivo relevant prioritization method for the endocrine potencies of different polycarbonate monomers, by combining in vitro bioassay data with physiologically based kinetic (PBK) modelling. PBK models were developed for a selection of monomers, including bisphenol A (BPA), two bisphenol F (BPF) isomers and four different bisphenol A diglycidyl ethers (BADGEs), using in vitro input data. With these models, the plasma concentrations of the compounds were simulated, providing means to estimate the dose levels at which the in vitro endocrine effect concentrations are reached. The results revealed that, whereas the in vitro relative potencies of different BADGEs (predominantly anti-androgenic effects) can be up to fourfold higher than BPA, the estimated in vivo potencies based on the oral equivalent doses are one to two orders of magnitude lower than BPA because of fast detoxification of the BADGEs. In contrast, the relative potencies of 2,2-BPF and 4,4-BPF increase when accounting for the in vivo availability. 4,4-BPF is estimated to be fivefold more potent than BPA in humans in vivo in inducing estrogenic effects and both 2,2-BPF and 4,4-BPF are estimated to be, respectively, 7 and 11-fold more potent in inducing anti-androgenic effects. These relative potencies were considered to be first-tier estimates, particularly given that the potential influence of intestinal metabolism on the in vivo availability was not accounted for. Overall, it can be concluded that both 2,2-BPF and 4,4-BPF are priority compounds.

Background Upon ingestion, nanoparticles can interact with the intestinal epithelial barrier potentially resulting in systemic uptake of nanoparticles. Nanoparticle properties have been described to influence the protein corona formation and subsequent cellular adhesion, uptake and transport. Here, we aimed to study the effects of nanoparticle size and surface chemistry on the protein corona formation and subsequent cellular adhesion, uptake and transport. Caco-2 intestinal cells, were exposed to negatively charged polystyrene nanoparticles (PSNPs) (50 and 200 nm), functionalized with sulfone or carboxyl groups, at nine nominal concentrations (15–250 μg/ml) for 10 up to 120 min. The protein coronas were analysed by LC–MS/MS. Results Subtle differences in the protein composition of the two PSNPs with different surface chemistry were noted. High-content imaging analysis demonstrated that sulfone PSNPs were associated with the cells to a significantly higher extent than the other PSNPs. The apparent cellular adhesion and uptake of 200 nm PSNPs was not significantly increased compared to 50 nm PSNPs with the same surface charge and chemistry. Surface chemistry outweighs the impact of size on the observed PSNP cellular associations. Also transport of the sulfone PSNPs through the monolayer of cells was significantly higher than that of carboxyl PSNPs. Conclusions The results suggest that the composition of the protein corona and the PSNP surface chemistry influences cellular adhesion, uptake and monolayer transport, which might be predictive of the intestinal transport potency of NPs.

Following legislation, European Member States should have multi-annual control programs for contaminants, such as for mycotoxins, in feed and food. These programs need to be risk based implying the checks are regular and proportional to the estimated risk for animal and human health. This study aimed to prioritize feed products in the Netherlands for deoxynivalenol and aflatoxin B1 monitoring. Historical mycotoxin monitoring results from the period 2007–2016 were combined with data from other sources. Based on occurrence, groundnuts had high priority for aflatoxin B1 monitoring; some feed materials (maize and maize products and several oil seed products) and complete/complementary feed excluding dairy cattle and young animals had medium priority; and all other animal feeds and feed materials had low priority. For deoxynivalenol, maize by-products had a high priority, complete and complementary feed for pigs had a medium priority and all other feed and feed materials a low priority. Also including health consequence estimations showed that feed materials that ranked highest for aflatoxin B1 included sunflower seed and palmkernel expeller/extracts and maize. For deoxynivalenol, maize products were ranked highest, followed by various small grain cereals (products); all other feed materials were of lower concern. Results of this study have proven to be useful in setting up the annual risk based control program for mycotoxins in animal feed and feed materials.

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.

Absorption, distribution, metabolism and excretion (ADME) of food-borne toxic compounds may be influenced by other compounds or constituents present in the food. The present review presents an overview of evidence currently available on food matrix-derived combination effects influencing the ADME characteristics of food-borne toxic compounds and the possible implications for risk assessment. The results obtained indicate that interactions may occur at all levels of ADME and that the interactions may decrease but also increase the bioavailability and/or toxicity of the compounds of interest. The overview also illustrates that food matrix-derived combination effects should be considered on a case-by-case basis, taking into account especially the mode of action underlying the interactions and the dose dependency of the effects. Especially food matrix-derived combination effects that proceed by a reversible mode of action, such as for example binding to biotransformation enzymes or transport proteins, may be detected at concentrations used in in vitro assays and at dose levels used in animal bioassays but may be absent at dose levels representing realistic human intake. It is concluded that although food matrix-derived combination effects may exist, their detection in in vitro assays or in animal bioassays at high dose levels may not improve risk assessment practice because interactions observed may not be maintained at low realistic levels of intake. Insight in the mode of action underlying the interactions combined with physiologically based kinetic (PBK) modelling may prove a way to obtain better insight in whether interactions detected at high dose levels will still be relevant at more realistic lower intake levels, and thus to what extent these effects should be taken into account in the risk assessment for human exposure.

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. © The Author 2017. Published by Oxford University Press on behalf of the Society of Toxicology.

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.5 h at different dose levels ranging from 0.1 to 50 mg/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.

ScopeThe present work investigates whether the previous observation that the basil flavonoid nevadensin is able to inhibit sulfotransferase (SULT)-mediated estragole DNA adduct formation in primary rat hepatocytes could be validated in vivo. Methods and resultsEstragole and nevadensin were co-administered orally to Sprague-Dawley rats, at a ratio reflecting their presence in basil. Moreover, previously developed physiologically based biokinetic (PBBK) models to study this inhibition in rat and in human liver were refined by including a submodel describing nevadensin kinetics. Nevadensin resulted in a significant 36% reduction in the levels of estragole DNA adducts formed in the liver of rats. The refined PBBK model predicts the formation of estragole DNA adducts in the liver of rat with less than twofold difference compared to in vivo data and suggests more potent inhibition in the liver of human compared to rat due to less efficient metabolism of nevadensin in human liver and intestine. ConclusionGiven the role of the SULT-mediated DNA adduct formation in the hepatocarcinogenicity of estragole, the results of the present study suggest that the likelihood of bioactivation and subsequent adverse effects in rodent bioassays may be lower when estragole is dosed with nevadensin compared to dosing of pure estragole.

A number of a,ß-unsaturated aldehydes are present in food both as natural constituents and as flavouring agents. Their reaction with DNA due to their electrophilic a,ß-unsaturated aldehyde moiety may result in genotoxicity as observed in some in vitro models, thereby raising a safety concern. A question that remains is whether in vivo detoxification would be efficient enough to prevent DNA adduct formation and genotoxicity. In this study, a human physiologically based kinetic/dynamic (PBK/D) model of trans-2-hexenal (2-hexenal), a selected model a,ß-unsaturated aldehyde, was developed to examine dose-dependent detoxification and DNA adduct formation in humans upon dietary exposure. The kinetic model parameters for detoxification were quantified using relevant pooled human tissue fractions as well as tissue fractions from 11 different individual subjects. In addition, a Monte Carlo simulation was performed so that the impact of interindividual variation in 2-hexenal detoxification on the DNA adduct formation in the population as a whole could be examined. The PBK/D model revealed that DNA adduct formation due to 2-hexenal exposure was 0.039 adducts/108 nucleotides (nt) at the estimated average 2-hexenal dietary intake (0.04 mg 2-hexenal/kg bw) and 0.18 adducts/108 nt at the 95th percentile of the dietary intake (0.178 mg 2-hexenal/kg bw) in the most sensitive people. These levels are three orders of magnitude lower than natural background DNA adduct levels that have been reported in disease-free humans (6.8–110 adducts/108 nt), suggesting that the genotoxicity risk for the human population at realistic dietary daily intakes of 2-hexenal may be negligible.

cope 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

Estragole is a known hepatocarcinogen in rodents at high doses following metabolic conversion to the DNA-reactive metabolite 1'-sulfooxyestragole. The aim of the present study was to model possible levels of DNA adduct formation in (individual) humans upon exposure to estragole. This was done by extending a previously defined PBK model for estragole in humans to include (i) new data on interindividual variation in the kinetics for the major PBK model parameters influencing the formation of 1'-sulfooxyestragole, (ii) an equation describing the relationship between 1'-sulfooxyestragole and DNA adduct formation, (iii) Monte Carlo modeling to simulate interindividual human variation in DNA adduct formation in the population, and (iv) a comparison of the predictions made to human data on DNA adduct formation for the related alkenylbenzene methyleugenol. Adequate model predictions could be made, with the predicted DNA adduct levels at the estimated daily intake of estragole of 0.01 mg/kg bw ranging between 1.6 and 8.8 adducts in 10(8) nucleotides (nts) (50th and 99th percentiles, respectively). This is somewhat lower than values reported in the literature for the related alkenylbenzene methyleugenol in surgical human liver samples. The predicted levels seem to be below DNA adduct levels that are linked with tumor formation by alkenylbenzenes in rodents, which were estimated to amount to 188-500 adducts per 10(8) nts at the BMD10 values of estragole and methyleugenol. Although this does not seem to point to a significant health concern for human dietary exposure, drawing firm conclusions may have to await further validation of the model's predictions.

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 (

The alkenylbenzene estragole is a constituent of several herbs and spices. It induces hepatomas in rodents at high doses following bioactivation by cytochrome P450s and sulfotransferases (SULTs) giving rise to the ultimate carcinogenic metabolite 1'-sulfooxyestragole which forms DNA adducts. Methanolic extracts from different alkenylbenzene-containing herbs and spices were able to inhibit SULT activity. Flavonoids including quercetin, kaempferol, myricetin, apigenin, and nevadensin were the major constituents responsible for this inhibition with Ki values in the nano to micromolar range. In human HepG2 cells exposed to the proximate carcinogen 1'-hydroxyestragole, the various flavonoids were able to inhibit estragole DNA adduct formation and shift metabolism in favor of glucuronidation which is a detoxification pathway for 1'-hydroxyestragole. In a next step, the kinetics for SULT inhibition were incorporated in physiologically based biokinetic (PBBK) models for estragole in rat and human to predict the effect of co-exposure to estragole and (mixtures of) the different flavonoids on the bioactivation in vivo. The PBBK-model-based predictions indicate that the reduction of estragole bioactivation in rat and human by co-administration of the flavonoids is dependent on whether the intracellular liver concentrations of the flavonoids can reach their Ki values. It is expected that this is most easily achieved for nevadensin which has a Ki value in the nanomolar range and is, due to its methyl ation, more metabolically stable than the other flavonoids.

trans-2-Hexenal (2-hexenal) is an α,β-unsaturated aldehyde that occurs naturally in a wide range of fruits, vegetables, and spices. 2-Hexenal as well as other α,β-unsaturated aldehydes that are natural food constituents or flavoring agents may raise a concern for genotoxicity due to the ability of the α,β-unsaturated aldehyde moiety to react with DNA. Controversy remains, however, on whether α,β-unsaturated aldehydes result in significant DNA adduct formation in vivo at realistic dietary exposure. In this study, a rat physiologically based in silico model was developed for 2-hexenal as a model compound to examine the time- and dose-dependent detoxification and DNA adduct formation of this selected α,β-unsaturated aldehyde. The model was developed based on in vitro and literature-derived parameters, and its adequacy was evaluated by comparing predicted DNA adduct formation in the liver of rats exposed to 2-hexenal with reported in vivo data. The model revealed that at an exposure level of 0.04 mg/kg body weight, a value reflecting estimated daily human dietary intake, 2-hexenal is rapidly detoxified predominantly by conjugation with glutathione (GSH) by glutathione S-transferases. At higher dose levels, depletion of GSH results in a shift to 2-hexenal oxidation and reduction as the major pathways for detoxification. The level of DNA adduct formation at current levels of human dietary intake was predicted to be more than 3 orders of magnitude lower than endogenous DNA adduct levels. These results support that rapid detoxification of 2-hexenal reduces the risk arising from 2-hexenal exposure and that at current dietary exposure levels, DNA adduct formation is negligible.

Estragole is a natural constituent of several herbs and spices including sweet basil. In rodent bioassays, estragole induces hepatomas, an effect ascribed to estragole bioactivation to 1'-sulfooxyestragole resulting in DNA adduct formation. The present paper identifies nevadensin as a basil constituent able to inhibit DNA adduct formation in rat hepatocytes exposed to the proximate carcinogen 1'-hydroxyestragole and nevadensin. This inhibition occurs at the level of sulfotransferase (SULT)-mediated bioactivation of 1'-hydroxyestragole. The Ki for SULT inhibition by nevadensin was 4 nM in male rat and human liver fractions. Furthermore, nevadensin up to 20 microM did not inhibit 1'-hydroxyestragole detoxification by glucuronidation and oxidation. The inhibition of SULT by nevadensin was incorporated into the recently developed physiologically based biokinetic (PBBK) rat and human models for estragole bioactivation and detoxification. The results predict that co-administration of estragole at a level inducing hepatic tumors in vivo (50mg/kg bw) with nevadensin at a molar ratio of 0.06, representing the ratio of their occurrence in basil, results in almost 100% inhibition of the ultimate carcinogen 1'-sulfooxyestragole when assuming 100% uptake of nevadensin. Assuming 1% uptake, inhibition would still amount to more than 83%. Altogether these data point at a nevadensin-mediated inhibition of the formation of the ultimate carcinogenic metabolite of estragole, without reducing the capacity to detoxify 1'-hydroxyestragole via glucuronidation or oxidation. These data also point at a potential reduction of the cancer risk when estragole exposure occurs within a food matrix containing SULT inhibitors compared to what is observed upon exposure to pure estragole.

The present study investigates interindividual variation in liver levels of the proximate carcinogenic metabolite of estragole, 1'-hydroxyestragole, due to variation in two key metabolic reactions involved in the formation and detoxification of this metabolite, namely 1'-hydroxylation of estragole and oxidation of 1'-hydroxyestragole. Formation of 1'-hydroxyestragole is predominantly catalyzed by P450 1A2, 2A6, and 2E1, and results of the present study support that oxidation of 1'-hydroxyestragole is catalyzed by 17beta-hydroxysteroid dehydrogenase type 2 (17beta-HSD2). In a first approach, the study defines physiologically based biokinetic (PBBK) models for 14 individual human subjects, revealing a 1.8-fold interindividual variation in the area under the liver concentration-time curve (AUC) for 1'-hydroxyestragole within this group of human subjects. Variation in oxidation of 1'-hydroxyestragole by 17beta-HSD2 was shown to result in larger effects than those caused by variation in P450 enzyme activity. In a second approach, a Monte Carlo simulation was performed to evaluate the extent of variation in liver levels of 1'-hydroxyestragole that could occur in the population as a whole. This analysis could be used to derive a chemical-specific adjustment factor (CSAF), which is defined as the 99th percentile divided by the 50th percentile of the predicted distribution of the AUC of 1'-hydroxyestragole in the liver. The CSAF was estimated to range between 1.6 and 4.0, depending on the level of variation that was taken into account for oxidation of 1'-hydroxyestragole. Comparison of the CSAF to the default uncertainty factor of 3.16 for human variability in biokinetics reveals that the default uncertainty factor adequately protects 99% of the population.

In chemical safety assessment, information on adverse effects after chronic exposure to low levels of hazardous compounds is essential for estimating human risks. Results from in vitro studies are often not directly applicable to the in vivo situation, and in vivo animal studies often have to be performed at unrealistic high levels of exposure. Physiologically based biokinetic (PBBK) modeling can be used as a platform for integrating in vitro metabolic data to predict dose- and species-dependent in vivo effects on biokinetics, and can provide a method to obtain a better mechanistic basis for extrapolations of data obtained in experimental animal studies to the human situation. Recently, we have developed PBBK models for the bioactivation of the alkenylbenzene estragole to its DNA binding ultimate carcinogenic metabolite 1'-sulfooxyestragole in both rat and human, as well as rat and human PBBK models for the bioactivation of coumarin to its hepatotoxic o-hydroxyphenylacetaldehyde metabolite. This article presents an overview of the results obtained so far with these in silico methods for PBBK modeling, focusing on the possible implications for risk assessment, and some additional considerations and future perspectives.

The present study defines a physiologically based biokinetic (PBBK) model for the alkenylbenzene estragole in rat based on in vitro metabolic parameters determined using relevant tissue fractions, in silico derived partition coefficients, and physiological parameters derived from the literature. The model consists of eight compartments including liver, lung and kidney as metabolizing compartments, and additional compartments for fat, arterial blood, venous blood, rapidly perfused tissue and slowly perfused tissue. Evaluation of the model was performed by comparing the PBBK predicted dose-dependent formation of the estragole metabolites 4-allylphenol and 1′-hydroxyestragole glucuronide to literature reported levels of these metabolites, which were demonstrated to be in the same order of magnitude. With the model obtained the relative extent of bioactivation and detoxification of estragole at different oral doses was examined. At low doses formation of 4-allylphenol, leading to detoxification, is observed to be the major metabolic pathway, occurring mainly in the lung and kidney due to formation of this metabolite with high affinity in these organs. Saturation of this metabolic pathway in the lung and kidney leads to a relative increase in formation of the proximate carcinogenic metabolite 1′-hydroxyestragole, occurring mainly in the liver. This relative increase in formation of 1′-hydroxyestragole leads to a relative increase in formation of 1′-hydroxyestragole glucuronide and 1′-sulfooxyestragole the latter being the ultimate carcinogenic metabolite of estragole. These results indicate that the relative importance of different metabolic pathways of estragole may vary in a dose-dependent way, leading to a relative increase in bioactiviation of estragole at higher doses. © 2008 Elsevier Inc. All rights reserved.

The effects of a basil extract on the sulfation and concomitant DNA adduct formation of the proximate carcinogen 1′-hydroxyestragole were studied using rat and human liver S9 homogenates and the human hepatoma cell line HepG2. Basil was chosen since it contains the procarcinogen estragole that can be metabolized to 1′-hydroxyestragole by cytochrome P450 enzymes. Basil extract addition to incubations of rat and human liver S9 homogenates with 1′-hydroxyestragole, the sulfotransferase cofactor PAPS, and 2′-deoxyguanosine resulted in a dose-dependent inhibition of N2-(trans-isoestragol-3′-yl)-2′-deoxyguanosine formation. Because the inhibition resembled the inhibition by the sulfotransferase inhibitor pentachlorophenol and since the inhibition was not observed in incubations with the direct electrophile 1′-acetoxyestragole it is concluded that the inhibition occurs at the level of the sulfotransferase mediated bioactivation step. Additional experiments in HepG2 cells revealed the same protective effect of basil extract in intact cells, demonstrating that the inhibitors are able to enter the cells. The results of this study suggest that bioactivation and subsequent adverse effects of 1′-hydroxyestragole might be lower in a matrix of other basil ingredients than what would be expected on the basis of experiments using 1′-hydroxyestragole as a single compound. © 2008 Elsevier Ltd. All rights reserved. Chemicals / CAS: 1' hydroxyestragole, 51410-44-7; cytochrome P450, 9035-51-2; deoxyguanosine, 961-07-9; DNA, 9007-49-2; estragole, 140-67-0; sulfotransferase, 9023-09-0; 1'-hydroxyestragole, 51410-44-7; Anisoles; Deoxyguanosine, 961-07-9; DNA Adducts; Enzyme Inhibitors; Plant Extracts; Sulfotransferases, EC 2.8.2.-; Tetrazolium Salts; Thiazoles; thiazolyl blue, 298-93-1

Aristolochic acids are naturally occurring nephrotoxins. This study aims to investigate whether physiologically based kinetic (PBK) model-based reverse dosimetry could convert in vitro concentration-response curves of aristolochic acid I (AAI) to in vivo dose response-curves for nephrotoxicity in rat, mouse and human. To achieve this extrapolation, PBK models were developed for AAI in these different species. Subsequently, concentration-response curves obtained from in vitro cytotoxicity models were translated to in vivo dose-response curves using PBK model-based reverse dosimetry. From the predicted in vivo dose-response curves, points of departure (PODs) for risk assessment could be derived. The PBK models elucidated species differences in the kinetics of AAI with the overall catalytic efficiency for metabolic conversion of AAI to aristolochic acid Ia (AAIa) being 2-fold higher for rat and 64-fold higher for mouse than human. Results show that the predicted PODs generally fall within the range of PODs derived from the available in vivo studies. This study provides proof of principle for a new method to predict a POD for in vivo nephrotoxicity by integrating in vitro toxicity testing with in silico PBK model-based reverse dosimetry.

To get more insight into possible species differences in the bioactivation of estragole, the kinetics for sulfonation of the proximate carcinogen 1'-hydroxyestragole were compared for male rat, male mouse, and mixed gender human liver S9 homogenates. In order to quantify sulfonation, 2'-deoxyguanosine was added to the incubation mixture in which sulfonation of 1'-hydroxyestragole was catalyzed to trap the reactive 1'-sulfooxyestragole. A method was developed with which the formation of the most abundant adduct with 2'-deoxyguanosine could be quantified using isotope dilution LC-ESI-MS/MS. Comparing the kinetics for sulfonation by liver S9 homogenates of male rat, male mouse, and humans revealed that sulfonation was about 30 times more efficient by male rat liver S9 than by human liver S9, whereas the catalytic efficiency by male mouse and human liver S9 was about the same. This indicates, as far as the bioactivation by sulfotransferase is concerned, that when extrapolating the cancer risk from laboratory animals to humans, using data from male rats may overestimate the cancer risk in humans, whereas using data from male mice may provide a better estimate of the cancer risk in humans.

Genotoxicity of α,β-unsaturated aldehydes shown in vitro raises a concern for the use of the aldehydes as food flavourings, while at low dose exposures the formation of DNA adducts may be prevented by detoxification. Unlike many α,β-unsaturated aldehydes for which in vivo data are absent, cinnamaldehyde was shown to be not genotoxic or carcinogenic in vivo. The present study aimed at comparing dose-dependent DNA adduct formation by cinnamaldehyde and 18 acyclic food-borne α,β-unsaturated aldehydes using physiologically based kinetic/dynamic (PBK/D) modelling. In rats, cinnamaldehyde was predicted to induce higher DNA adducts levels than 6 out of the 18 α,β-unsaturated aldehydes, indicating that these 6 aldehydes may also test negative in vivo. At the highest cinnamaldehyde dose that tested negative in vivo, cinnamaldehyde was predicted to form at least three orders of magnitude higher levels of DNA adducts than the 18 aldehydes at their respective estimated daily intake. These results suggest that for all the 18 α,β-unsaturated aldehydes DNA adduct formation at doses relevant for human dietary exposure may not raise a concern. The present study illustrates a possible use of physiologically based in silico modelling to facilitate a science-based comparison and read-across on the possible risks posed by DNA reactive agents.

In the present study, the effects of a large series of flavonoids on multidrug resistance proteins (MRPs) were studied in MRP1 and MRP2 transfected MDCKII cells. The results were used to define the structural requirements of flavonoids necessary for potent inhibition of MRP1- and MRP2-mediated calcein transport in a cellular model. Several of the methoxylated flavonoids are among the best MRP1 inhibitors (IC50 values, ranging between 2.7 and 14.3 μM) followed by robinetin, myricetin and quercetin (IC50 values ranging between 13.6 and 21.8 μM). Regarding inhibition of MRP2 activity especially robinetin and myricetin appeared to be good inhibitors (IC 50 values of 15.0 and 22.2 μM, respectively). Kinetic characterization revealed that the two transporters differ marginally in the apparent Km for the substrate calcein. For one flavonoid, robinetin, the kinetics of inhibition were studied in more detail and revealed competitive inhibition with respect to calcein, with apparent inhibition constants of 5.0 μM for MRP1 and 8.5 μM for MRP2. For inhibition of MRP1, a quantitative structure activity relationship (QSAR) was obtained that indicates three structural characteristics to be of major importance for MRP1 inhibition by flavonoids: the total number of methoxylated moieties, the total number of hydroxyl groups and the dihedral angle between the B- and C-ring. Regarding MRP2 mediated calcein efflux inhibition, only the presence of a flavonol B-ring pyrogallol group seems to be an important structural characteristic. Overall, this study provides insight in the structural characteristics involved in MRP inhibition and explores the differences between inhibitors of these two transporters, MRP1 and MRP2. Ultimately, MRP2 displays higher selectivity for flavonoid type inhibition than MRP1. © 2004 Elsevier Inc. All rights reserved.

Acyclic α,β-unsaturated aldehydes present in food raise a concern because the α,β-unsaturated aldehyde moiety is considered a structural alert for genotoxicity. However, controversy remains on whether in vivo at realistic dietary exposure DNA adduct formation is significant. The aim of the present study was to develop physiologically based kinetic/dynamic (PBK/D) models to examine dose-dependent detoxification and DNA adduct formation of a group of 18 food-borne acyclic α,β-unsaturated aldehydes without 2- or 3-alkylation, and with no more than one conjugated double bond. Parameters for the PBK/D models were obtained using quantitative structure-activity relationships (QSARs) defined with a training set of six selected aldehydes. Using the QSARs, PBK/D models for the other 12 aldehydes were defined. Results revealed that DNA adduct formation in the liver increases with decreasing bulkiness of the molecule especially due to less efficient detoxification. 2-Propenal (acrolein) was identified to induce the highest DNA adduct levels. At realistic dietary intake, the predicted DNA adduct levels for all aldehydes were two orders of magnitude lower than endogenous background levels observed in disease free human liver, suggesting that for all 18 aldehydes DNA adduct formation is negligible at the relevant levels of dietary intake. The present study provides a proof of principle for the use of QSAR-based PBK/D modelling to facilitate group evaluations and read-across in risk assessment.

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.

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 50 mg/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

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.

Methyleugenol (ME) occurs naturally in a variety of spices, herbs, including basil, and their essential oils. ME induces hepatomas in rodent bioassays following its conversion to a DNA reactive metabolite. In the present study, the basil constituent nevadensin was shown to be able to inhibit SULT-mediated DNA adduct formation in HepG2 cells exposed to the proximate carcinogen 1'-hydroxymethyleugenol in the presence of nevadensin. To investigate possible in vivo implications of SULT inhibition by nevadensin on ME bioactivation, the rat physiologically based kinetic (PBK) model developed in our previous work to describe the dose-dependent bioactivation and detoxification of ME in male rat was combined with the recently developed PBK model describing the dose-dependent kinetics of nevadensin in male rat. The resulting binary ME–nevadensin PBK model was used to predict the possible nevadensin mediated reduction in ME DNA adduct formation and resulting carcinogenicity at the doses of ME used by the NTP carcinogenicity study. Using these data an updated risk assessment using the Margin of Exposure (MOE) approach was performed. The results obtained point at a potential reduction of the cancer risk when rodents are orally exposed to ME within a relevant food matrix containing SULT inhibitors compared to exposure to pure ME.

A physiologically based biokinetic (PBBK) model for the alkenylbenzene safrole in humans was developed based on in vitro- and in silico-derived kinetic parameters. With the model obtained, the time- and dose-dependent formation of the proximate and ultimate carcinogenic metabolites, 1'-hydroxysafrole and 1'-sulfooxysafrole in human liver were estimated and compared with previously predicted levels of these metabolites in rat liver. In addition, Monte Carlo simulations were performed to predict interindividual variation in the formation of these metabolites in the overall population. For the evaluation of the model performance, a comparison was made between the predicted total amount of urinary metabolites of safrole and the reported total levels of metabolites in the urine of humans exposed to safrole, which adequately matched. The model results revealed no dose-dependent shifts in safrole metabolism and no relative increase in bioactivation at dose levels up to 100 mg/kg body weight/day. Species differences were mainly observed in the detoxification pathways of 1'-hydroxysafrole, with the formation of 1'-oxosafrole being a main detoxification pathway of 1'-hydroxysafrole in humans but a minor pathway in rats, and glucuronidation of 1'-hydroxysafrole being less important in humans than in rats. The formation of 1'-sulfooxysafrole was predicted to vary 4- to 17-fold in the population (fold difference between the 95th and median, and 95th and 5th percentile, respectively), with the median being three to five times higher in human than in rat liver. Comparison of the PBBK results for safrole with those previously obtained for the related alkenylbenzenes estragole and methyleugenol revealed that differences in 1'-sulfooxy metabolite formation are limited, being only twofold to fivefold.

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.

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.

A physiologically based biokinetic (PBBK) model for alkenylbenzene safrole in rats was developed using in vitro metabolic parameters determined using relevant tissue fractions. The performance of the model was evaluated by comparison of the predicted levels of 1,2-dihydroxy-4-allylbenzene and 1'-hydroxysafrole glucuronide to levels of these metabolites reported in the literature to be excreted in the urine of rats exposed to safrole and by comparison of the predicted amount of total urinary safrole metabolites to the reported levels of safrole metabolites in the urine of safrole exposed rats. These comparisons revealed that the predictions adequately match observed experimental values. Next, the model was used to predict the relative extent of bioactivation and detoxification of safrole at different oral doses. At low as well as high doses, P450 mediated oxidation of safrole mainly occurs in the liver in which 1,2-dihydroxy-4-allylbenzene was predicted to be the major P450 metabolite of safrole. A dose dependent shift in P450 mediated oxidation leading to a relative increase in bioactivation at high doses was not observed. Comparison of the results obtained for safrole with the results previously obtained with PBBK models for the related alkenylbenzenes estragole and methyleugenol revealed that the overall differences in bioactivation of the three alkenylbenzenes to their ultimate carcinogenic 1'-sulfooxy metabolites are limited. This is in line with the generally less than 4-fold difference in their level of DNA binding in in vitro and in vivo studies and their almost similar BMDL(10) values (lower confidence limit of the benchmark dose that gives 10% increase in tumor incidence over background level) obtained in in vivo carcinogenicity studies. It is concluded that in spite of differences in the rates of specific metabolic conversions, overall the levels of bioactivation of the three alkenylbenzenes are comparable which is in line with their comparable carcinogenic potential.

The present study defines a physiologically based biokinetic (PBBK) model for the alkenylbenzene methyleugenol in rat based on in vitro metabolic parameters determined using relevant tissue fractions, in silico derived partition coefficients, and physiological parameters derived from the literature. The model was based on the model previously developed for the related alkenylbenzene estragole and consists of eight compartments including liver, lung, and kidney as metabolizing compartments, and separate compartments for fat, arterial blood, venous blood, richly perfused and slowly perfused tissues. Evaluation of the model was performed by comparing the PBBK predicted concentration of methyleugenol in the venous compartment to methyleugenol plasma levels reported in the literature, by comparing the PBBK predicted dose-dependent percentage of formation of 2-hydroxy-4,5-dimethoxyallylbenzene, 3-hydroxy-4-methoxyallylbenzene, and 1'-hydroxymethyleugenol glucuronide to the corresponding percentage of metabolites excreted in urine reported in the literature, which were demonstrated to be in the same order of magnitude. With the model obtained the relative extent of bioactivation and detoxification of methyleugenol at different oral doses was examined. At low doses, formation of 3-(3,4-dimethoxyphenyl)-2-propen-1-ol and methyleugenol-2',3'-oxide leading to detoxification appear to be the major metabolic pathways, occurring in the liver. At high doses, the model reveals a relative increase in the formation of the proximate carcinogenic metabolite 1'-hydroxymethyleugenol, occurring in the liver. This relative increase in formation of 1'-hydroxymethyleugenol leads to a relative increase in formation of 1'-hydroxymethyleugenol glucuronide, 1'-oxomethyleugenol, and 1'-sulfooxymethyleugenol the latter being the ultimate carcinogenic metabolite of methyleugenol. These results indicate that the relative importance of different metabolic pathways of methyleugenol may vary in a dose-dependent way, leading to a relative increase in bioactiviation of methyleugenol at higher doses.

Estragole has been shown to be hepatocarcinogenic in rodent species at high-dose levels. Translation of these results into the likelihood of formation of DNA adducts, mutation, and ultimately cancer upon more realistic low-dose exposures remains a challenge. Recently we have developed physiologically based biokinetic (PBBK) models for rat and human predicting bioactivation of estragole. These PBBK models, however, predict only kinetic characteristics. The present study describes the extension of the PBBK model to a so-called physiologically based biodynamic (PBBD) model predicting in vivo DNA adduct formation of estragole in rat liver. This PBBD model was developed using in vitro data on DNA adduct formation in rat primary hepatocytes exposed to 1'-hydroxyestragole. The model was extended by linking the area under the curve for 1'-hydroxyestragole formation predicted by the PBBK model to the area under the curve for 1'-hydroxyestragole in the in vitro experiments. The outcome of the PBBD model revealed a linear increase in DNA adduct formation with increasing estragole doses up to 100 mg/kg bw. Although DNA adduct formation of genotoxic carcinogens is generally seen as a biomarker of exposure rather than a biomarker of response, the PBBD model now developed is one step closer to the ultimate toxic effect of estragole than the PBBK model described previously. Comparison of the PBBD model outcome to available data showed that the model adequately predicts the dose-dependent level of DNA adduct formation. The PBBD model predicts DNA adduct formation at low levels of exposure up to a dose level showing to cause cancer in rodent bioassays, providing a proof of principle for modeling a toxicodynamic in vivo endpoint on the basis of solely in vitro experimental data.

The extent of bioactivation of the herbal constituent estragole to its ultimate carcinogenic metabolite 1′-sulfooxyestragole depends on the relative levels of bioactivation and detoxification pathways. The present study investigated the kinetics of the metabolic reactions of both estragole and its proximate carcinogenic metabolite 1′-hydroxyestragole in humans in incubations with relevant tissue fractions. Based on the kinetic data obtained a physiologically based biokinetic (PBBK) model for estragole in human was defined to predict the relative extent of bioactivation and detoxification at different dose levels of estragole. The outcomes of the model were subsequently compared with those previously predicted by a PBBK model for estragole in male rat to evaluate the occurrence of species differences in metabolic activation. The results obtained reveal that formation of 1′-oxoestragole, which represents a minor metabolic route for 1′-hydroxyestragole in rat, is the main detoxification pathway of 1′-hydroxyestragole in humans. Due to a high level of this 1′-hydroxyestragole oxidation pathway in human liver, the predicted species differences in formation of 1′-sulfooxyestragole remain relatively low, with the predicted formation of 1′-sulfooxyestragole being twofold higher in human compared with male rat, even though the formation of its precursor 1′-hydroxyestragole was predicted to be fourfold higher in human. Overall, it is concluded that in spite of significant differences in the relative extent of different metabolic pathways between human and male rat there is a minor influence of species differences on the ultimate overall bioactivation of estragole to 1′-sulfooxyestragole.

Coumarin is bioactivated via 3,4-epoxidation resulting in formation of the hepatotoxic o-hydroxyphenylacetaldehyde (oHPA) and detoxified by cytochrome P450 2A6 (CYP2A6) hydroxylation leading to 7-hydroxycoumarin. The present study defines physiologically based biokinetic (PBBK) models to predict liver levels of the toxic oHPA metabolite of coumarin in rats and in human subjects with normal or deficient CYP2A6 catalyzed coumarin 7-hydroxylation. The results reveal that the predicted maximum tissue concentration (C(max)) of oHPA in the liver of wild type human subjects and of subjects deficient in CYP2A6 catalyzed 7-hydroxylation are, respectively, three and one order of magnitude lower than the values predicted for rat liver. Another difference between CYP2A6 deficient and wild type human subjects is a 500-fold difference in the area under the curve 0-24h (AUC(0-24h)) for the time-dependent oHPA liver concentration pointing at a relative higher percentage of the original dose converted in time through this pathway when CYP2A6 is deficient. For wild type human subjects and the subjects with completely deficient coumarin 7-hydroxylation the AUC(0-24h) values for oHPA in the liver are, respectively, three and one order of magnitude lower than that for rat liver. Even when 7-hydroxylation is deficient, the chances on formation of the hepatotoxic oHPA metabolite will be significantly lower in the liver of humans than those expected in the liver of rats when exposed to a similar dose on a body-weight basis. This conclusion should be taken into account when extrapolating data from experimental studies in sensitive animals, i.e., rats, to the general human population.