Search

Calibration of a kinetic model for the transfer of PCDD/Fs and dl-PCBs from feed to the hen's body and eggs was thus far restricted to the total TEQ concentration, i.e. the summed concentrations of PCDD/Fs and dl-PCBs expressed in terms of equivalents of 2,3,7,8-TCDD. However, this approach may lead to over- or underestimation of the transfer if the mixture contains congeners with kinetic characteristics which differ considerably from those used in such a model. This paper extends a previous transfer model of PCDD/Fs and dl-PCBs from feed to egg yolk fat and abdominal fat of high production laying hens, based on the total TEQ approach, to the level of individual congeners. Both modelling approaches are compared and the new approach is presented as a webtool application. This congener-specific approach enabled the calibration of 25 of the 29 relevant PCDD/F and dl-PCB congeners with respect to their individual transfer characteristics to body fat and egg yolk fat and their clearance from the body. Limitations of the available experimental data prevented the calibration of 1,2,3,4,6,7,8-HpCDD, OCDD, OCDF and PCB 123. The fraction transferred to egg yolk fat after long-term daily intake of contaminated feed was found to be at least 0.78 for 2,3,7,8-TCDD, 0.75 for PeCDD, 0.42-0.61 for HxCDDs, 0.70 for 2,3,7,8-TCDF, 0.71 for PeCDF, 0.54-0.60 for HxCDFs, 0.18-0.24 for HpCDFs and 0.89-1.00 for dl-PCBs. Various experimental and feed incident mixtures were used to compare the total TEQ- model with the congener-specific approach. An overestimation of the transfer by the total TEQ method was shown in particular for mixtures with a substantial contribution of hexa-, hepta- and octa-PCDD/Fs to the total TEQ level.

Contains fulltext : 291310.pdf (Publisher’s version ) (Open Access) Poly- and perfluoroalkyl substances (PFASs) are omnipresent in the environment and have been shown to accumulate in humans. Most PFASs are not biotransformed in animals and humans, so that elimination is largely dependent on non-metabolic clearance via bile and urine. Accumulation of certain PFASs in humans may relate to their reabsorption from the pre-urine by transporter proteins in the proximal tubules of the kidney, such as URAT1 and OAT4. The present study assessed the in vitro transport of 7 PFASs (PFHpA, PFOA, PFNA, PFDA, PFBS, PFHxS and PFOS) applying URAT1- or OAT4-transfected human embryonic kidney (HEK) cells. Virtually no transport of PFASs could be measured in URAT1-transfected HEK cells. All PFASs, except PFBS, showed clear uptake in OAT4-transfected HEK cells. In addition, these in vitro results were further supported by in silico docking and molecular dynamic simulation studies assessing transporter-ligand interactions. Information on OAT4-mediated transport may provide insight into the accumulation potential of PFASs in humans, but other kinetic aspects may play a role and should also be taken into account. Quantitative information on all relevant kinetic processes should be integrated in physiologically based kinetic (PBK) models, to predict congener-specific accumulation of PFASs in humans in a more accurate manner.

With an increasing need to incorporate new approach methodologies (NAMs) in chemical risk assessment and the concomitant need to phase out animal testing, the interpretation of in vitro assay readouts for quantitative hazard characterisation becomes more important. Physiologically based kinetic (PBK) models, which simulate the fate of chemicals in tissues of the body, play an essential role in extrapolating in vitro effect concentrations to in vivo bioequivalent exposures. As PBK-based testing approaches evolve, it will become essential to standardise PBK modelling approaches towards a consensus approach that can be used in quantitative in vitro-to-in vivo extrapolation (QIVIVE) studies for regulatory chemical risk assessment based on in vitro assays. Based on results of an ECETOC expert workshop, steps are recommended that can improve regulatory adoption: (1) define context and implementation, taking into consideration model complexity for building fit-for-purpose PBK models, (2) harmonise physiological input parameters and their distribution and define criteria for quality chemical-specific parameters, especially in the absence of in vivo data, (3) apply Good Modelling Practices (GMP) to achieve transparency and design a stepwise approach for PBK model development for risk assessors, (4) evaluate model predictions using alternatives to in vivo PK data including read-across approaches, (5) use case studies to facilitate discussions between modellers and regulators of chemical risk assessment. Proof-of-concepts of generic PBK modelling approaches are published in the scientific literature at an increasing rate. Working on the previously proposed steps is, therefore, needed to gain confidence in PBK modelling approaches for regulatory use.

The goal of the present study was to assess the predictive performance of a minimal generic rat physiologically based kinetic (PBK) model based on in vitro and in silico input data to predict peak plasma concentrations (Cmax) upon single oral dosing. To this purpose, a dataset was generated of 3960 Cmax predictions for 44 compounds, applying different combinations of in vitro and in silico approaches for chemical parameterization, and comparison of the predictions to reported in vivo data. Best performance was obtained when (1) the hepatic clearance was parameterized based on in vitro measured intrinsic clearance values, (2) the method of Rodgers and Rowland for calculating partition coefficients, and (3) in silico calculated fraction unbound plasma and Papp values (the latter especially for very lipophilic compounds). Based on these input data, the median Cmax of 32 compounds could be predicted within 10-fold of the observed Cmax, with 22 out of these 32 compounds being predicted within 5-fold, and 8 compounds within 2-fold. Overestimations of more than 10-fold were observed for 12 compounds, whereas no underestimations of more than 10-fold occurred. Median Cmax predictions were frequently found to be within 10-fold of the observed Cmax when the scaled unbound hepatic intrinsic clearance (Clint,u) was either higher than 20 l/h or lower than 1 l/h. Similar findings were obtained with a test set of 5 in-house BASF compounds. Overall, this study provides relevant insights in the predictive performance of a minimal PBK model based on in vitro and in silico input data.

Poly- and perfluoroalkyl substances (PFASs) are omnipresent in the environment and have been shown to accumulate in humans. Most PFASs are not biotransformed in animals and humans, so that elimination is largely dependent on non-metabolic clearance via bile and urine. Accumulation of certain PFASs in humans may relate to their reabsorption from the pre-urine by transporter proteins in the proximal tubules of the kidney, such as URAT1 and OAT4. The present study assessed the in vitro transport of 7 PFASs (PFHpA, PFOA, PFNA, PFDA, PFBS, PFHxS and PFOS) applying URAT1- or OAT4-transfected human embryonic kidney (HEK) cells. Virtually no transport of PFASs could be measured in URAT1-transfected HEK cells. All PFASs, except PFBS, showed clear uptake in OAT4-transfected HEK cells. In addition, these in vitro results were further supported by in silico docking and molecular dynamic simulation studies assessing transporter-ligand interactions. Information on OAT4-mediated transport may provide insight into the accumulation potential of PFASs in humans, but other kinetic aspects may play a role and should also be taken into account. Quantitative information on all relevant kinetic processes should be integrated in physiologically based kinetic (PBK) models, to predict congener-specific accumulation of PFASs in humans in a more accurate manner.

The intention of this study was to determine the utility of high-throughput screening (HTS) data, as exemplified by ToxCast and Tox21, for application in toxicological read-across in food-relevant chemicals. Key questions were addressed on the extent to which the HTS data could provide information enabling (1) the elucidation of underlying bioactivities associated with apical toxicological outcomes, (2) the closing of existing toxicological data gaps, and (3) the definition of the boundaries of chemical space across which bioactivity could reliably be extrapolated. Results revealed that many biological targets apparently activated within the chemical groupings lack, at this time, validated toxicity pathway associations. Therefore, as means of providing proof-of-principle, a comparatively well-characterized end point - estrogenicity - was selected for evaluation. This was facilitated through the preparation of two exploratory case studies, focusing upon groupings of paraben-gallates and pyranone-type compounds (notably flavonoids). Within both, the HTS data were seen to reflect estrogenic potencies in a manner which broadly corresponded to established structure-activity group relationships, with parabens and flavonoids displaying greater estrogen receptor affinity than benzoate esters and alternative pyranone-containing molecules, respectively. As such, utility in the identification of out-of-domain compounds was demonstrated, indicating potential for application in addressing point (3) as detailed above.

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.

Tox21 and ToxCast are high-throughput in vitro screening programs coordinated by the U.S. National Toxicology Program and the U.S. Environmental Protection Agency, respectively, with the goal of forecasting biological effects in vivo based on bioactivity profiling. The present study investigated whether mechanistic insights in the biological targets of food-relevant chemicals can be obtained from ToxCast results when the chemicals are grouped according to structural similarity. Starting from the 556 direct additives that have been identified in the ToxCast database by Karmaus et al. [Karmaus, A. L., Trautman, T. D., Krishan, M., Filer, D. L., and Fix, L. A. (2017). Curation of food-relevant chemicals in ToxCast. Food Chem. Toxicol. 103, 174–182.], the results showed that, despite the limited number of assays in which the chemical groups have been tested, sufficient results are available within so-called “DNA binding” and “nuclear receptor” target families to profile the biological activities of the defined chemical groups for these targets. The most obvious activity identified was the estrogen receptor-mediated actions of the chemical group containing parabens and structurally related gallates, as well the chemical group containing genistein and daidzein (the latter 2 being particularly active toward estrogen receptor β as a potential health benefit). These group effects, as well as the biological activities of other chemical groups, were evaluated in a series of case studies. Overall, the results of the present study suggest that high-throughput screening data could add to the evidence considered for regulatory risk assessment of food chemicals and to the evaluation of desirable effects of nutrients and phytonutrients. The data will be particularly useful for providing mechanistic information and to fill data gaps with read-across.

For almost fifteen years, the availability and regulatory acceptance of new approach methodologies (NAMs) to assess the absorption, distribution, metabolism and excretion (ADME/biokinetics) in chemical risk evaluations are a bottleneck. To enhance the field, a team of 24 experts from science, industry, and regulatory bodies, including new generation toxicologists, met at the Lorentz Centre in Leiden, The Netherlands. A range of possibilities for the use of NAMs for biokinetics in risk evaluations were formulated (for example to define species differences and human variation or to perform quantitative in vitro-in vivo extrapolations). To increase the regulatory use and acceptance of NAMs for biokinetics for these ADME considerations within risk evaluations, the development of test guidelines (protocols) and of overarching guidance documents is considered a critical step. To this end, a need for an expert group on biokinetics within the Organisation of Economic Cooperation and Development (OECD) to supervise this process was formulated. The workshop discussions revealed that method development is still required, particularly to adequately capture transporter mediated processes as well as to obtain cell models that reflect the physiology and kinetic characteristics of relevant organs. Developments in the fields of stem cells, organoids and organ-on-a-chip models provide promising tools to meet these research needs in the future.

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.

Non-animal methods for toxicokinetics, such as in vitro hepatic metabolic clearance studies, play an important role in chemical risk evaluations. To gain regulatory acceptance of such clearance data, the development of a test guideline for performing in vitro hepatic clearance studies is crucial. The aim of the present study was to obtain insight in the experimental conditions of clearance studies that influence obtained intrinsic clearance (CLint) values. To that end, in vitro hepatic CLint data obtained with rat or human hepatocytes and methodological aspects of the experiments, were collected from 42 different suitable studies published between 1995 and 2018. The CLint values for the majority of chemicals differed by more than one order of magnitude. We estimated the systematic effect of different experimental setups on the CLint values using a random forest regression analysis, revealing that 'hepatocyte concentration', 'species' (rat or human hepatocytes) and 'culture medium' have the largest impact. Calculating unbound CLint (CLint,u) values slightly reduced the variation for most chemicals. Given that in vivo clearance is in general underpredicted based on in vitro clearance data, a harmonized protocol is preferably based on a protocol that provides relatively high in vitro CLint values.

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.

Rigorous peer-review is the corner-stone of high-quality academic publishing [...]

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.

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.

Quantitative predictions of in vivo chemical levels based on in vitro data will become a cornerstone of next generation nonanimal risk evaluations. Both regulatory and scientific experience with quantitative toxicokinetics must increase now for this transition to happen.

The goal of the current work was to implement and verify previously published rat and human PBPK modelling codes for TCDD into Berkeley Madonna. The US‐EPA has used these PBPK models in the reassessment of TCDD. A procurement contract has been set up to explore the possibilities to adequately run the models and reproduce previously published results. The implementation of the available codes in Berkeley Madonna was carried out at RIKILT‐WUR under the framework agreement with RIVM. The results obtained with a newly coded PBPK model for TCDD in the rat in Berkeley Madonna revealed good reproducibility of previously reported graphs. A selection of key figures could be adequately simulated with the recoded models. Altogether it can be concluded that the PBPK codes for TCDD in rat and human can be adequately run in Berkeley Madonna.

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.

The objective of the present work was to review the availability and predictive value of non-animal toxicokinetic approaches and to evaluate their current use in European risk evaluations of food contaminants, additives and food contact materials, as well as pesticides and medicines. Results revealed little use of quantitative animal or human kinetic data in risk evaluations of food chemicals, compared with pesticides and medicines. Risk evaluations of medicines provided sufficient in vivo kinetic data from different species to evaluate the predictive value of animal kinetic data for humans. These data showed a relatively poor correlation between the in vivo bioavailability in rats and dogs versus that in humans. In contrast, in vitro (human) kinetic data have been demonstrated to provide adequate predictions of the fate of compounds in humans, using appropriate in vitro-in vivo scalers and by integration of in vitro kinetic data with in silico kinetic modelling. Even though in vitro kinetic data were found to be occasionally included within risk evaluations of food chemicals, particularly results from Caco-2 absorption experiments and in vitro data on gut-microbial conversions, only minor use of in vitro methods for metabolism and quantitative in vitro-in vivo extrapolation methods was identified. Yet, such quantitative predictions are essential in the development of alternatives to animal testing as well as to increase human relevance of toxicological risk evaluations. Future research should aim at further improving and validating quantitative alternative methods for kinetics, thereby increasing regulatory acceptance of non-animal kinetic data.

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

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.

The in vivo oestrogenicity of genistein and its glycoside genistin is still under debate. The present study aimed to develop a physiologically based kinetic (PBK) model that provides insight in dose-dependent plasma concentrations of genistein aglycone and its metabolites and enables prediction of in vivo oestrogenic effective dose levels of genistein and genistin in humans.A PBK model for genistein and genistin in humans was developed based on in vitro metabolic parameters. The model obtained was used to translate in vitro oestrogenic concentration-response curves of genistein to in vivo oestrogenic dose-response curves for intake of genistein and genistin.The model predicted that genistein-7-O-glucuronide was the major circulating metabolite and that levels of the free aglycone were generally low [0.5-17% of total plasma genistein at oral doses from 0.01 to 50 mg (kg·bw)The present study shows how plasma concentrations of genistein and its metabolites and oestrogenic dose levels of genistein in humans can be predicted by combining in vitro oestrogenicity with PBK model-based reverse dosimetry, eliminating the need for human intervention studies.

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.

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.

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.

The present paper aims at identifying strategies to increase the impact and applicability of alternative testing strategies in risk assessment. To this end, a quantitative and qualitative literature evaluation was performed on (a) current research efforts in the development of in vitro methods aiming for alternatives to animal testing, (b) the possibilities and limitations of in vitro methods for regulatory purposes and (c) the potential of physiologically-based kinetic (PBK) modeling to improve the impact and applicability of in vitro methods in risk assessment practice. Overall, the evaluation showed that the focus of state-of-the-art research activities does not seem to be optimally directed at developing in vitro alternatives for those endpoints that are most animal-demanding, such as reproductive and developmental toxicity, and carcinogenicity. A key limitation in the application of in vitro alternatives to such systemic endpoints is that in vitro methods do not provide so-called points of departure, necessary for regulators to set safe exposure limits. PBK-modeling could contribute to overcoming this limitation by providing a method that allows extrapolation of in vitro concentration-response curves to in vivo dose-response curves. However, more proofs of principle are required. (C) 2011 Elsevier Inc. All rights reserved.

Studies in the field of molecular nutrition and food research often aim at identifying effects of bioactive ingredients on living organisms. When data from human studies are difficult to obtain, effects are often studied in relevant animal or cellular in vitro models. This poses the need for adequate extrapolation from the in vitro to the in vivo situation, from high-dose levels to realistic low-dose levels and from experimental animals to humans. Furthermore, effects of genetic polymorphisms or lifestyle factors may have to be taken into account. Physiologically based kinetic (PBK) modeling provides a means to support these kinds of extrapolations. The present paper illustrates the basic concepts of PBK modeling. PBK modeling includes six steps: (i) definition of the conceptual model, (ii) translation into a mathematical model, (iii) defining parameter values, (iv) solving the equations, (v) evaluation of model performance and (vi) making predictions. The paper provides an overview of these basic steps and presents examples to illustrate how PBK modeling can be applied. This reveals that PBK modeling provides an important tool in the field of the 3Rs aiming at Replacement, Reduction and Refinement of animal studies and may also be a useful tool for risk assessment.

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

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.

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.