Use of physiologically based biokinetic (PBBK) modeling to study estragole bioactivation and detoxification in humans as compared with male rats

Toxins in Botanical Drugs and Plant-derived Food and Feed - from Science to Regulation: A Workshop Review

In September 2022, the 3rd International Workshop on pyrrolizidine alkaloids (PAs) and related phytotoxins was held on-line, entitled 'Toxins in botanical drugs and plant-derived food and feed - from science to regulation'. The workshop focused on new findings about the occurrence, exposure, toxicity, and risk assessment of PAs. In addition, new scientific results related to the risk assessment of alkenylbenzenes, a distinct class of herbal constituents, were presented. The presence of PAs and alkenylbenzenes in plant-derived food, feed, and herbal medicines has raised health concerns with respect to their acute and chronic toxicity but mainly related to the genotoxic and carcinogenic properties of several congeners. The compounds are natural constituents of a variety of plant families and species widely used in medicinal, food, and feed products. Their individual occurrence, levels, and toxic properties, together with the broad range of congeners present in nature, represent a striking challenge to modern toxicology. This review tries to provide an overview of the current knowledge on these compounds and indicates needs and perspectives for future research.

FEMA GRAS assessment of derivatives of basil, nutmeg, parsley, tarragon and related allylalkoxybenzene-containing natural flavor complexes

In 2015, the Expert Panel of the Flavor and Extract Manufacturers Association (FEMA) initiated a program for the re-evaluation of the safety of over 250 natural flavor complexes (NFCs) used as flavoring ingredients in food. In this publication, tenth in the series, NFCs containing a high percentage of at least one naturally occurring allylalkoxybenzene constituent with a suspected concern for genotoxicity and/or carcinogenicity are evaluated. In a related paper, ninth in the series, NFCs containing anethole and/or eugenol and relatively low percentages of these allylalkoxybenzenes are evaluated. The Panel applies the threshold of toxicological concern (TTC) concept and evaluates relevant toxicology data on the NFCs and their respective constituent congeneric groups. For NFCs containing allylalkoxybenzene constituent(s), the estimated intake of the constituent is compared to the TTC for compounds with structural alerts for genotoxicity and when exceeded, a margin of exposure (MOE) is calculated. BMDL₁₀ values are derived from benchmark dose analyses using Bayesian model averaging for safrole, estragole and methyl eugenol using EPA's BMDS software version 3.2. BMDL₁₀ values for myristicin, elemicin and parsley apiole were estimated by read-across using relative potency factors. Margins of safety for each constituent congeneric group and MOEs for each allylalkoxybenzene constituent for each NFC were determined that indicate no safety concern. The scope of the safety evaluation contained herein does not include added use in dietary supplements or any products other than food. Ten NFCs, derived from basil, estragon (tarragon), mace, nutmeg, parsley and Canadian snakeroot were determined or affirmed as generally recognized as safe (GRAS) under their conditions of intended use as flavor ingredients based on an evaluation of each NFC and the constituents and congeneric groups therein.

Identifying xenobiotic metabolites with in silico prediction tools and LCMS suspect screening analysis

Understanding the metabolic fate of a xenobiotic substance can help inform its potential health risks and allow for the identification of signature metabolites associated with exposure. The need to characterize metabolites of poorly studied or novel substances has shifted exposure studies towards non-targeted analysis (NTA), which often aims to profile many compounds within a sample using high-resolution liquid-chromatography mass-spectrometry (LCMS). Here we evaluate the suitability of suspect screening analysis (SSA) liquid-chromatography mass-spectrometry to inform xenobiotic chemical metabolism. Given a lack of knowledge of true metabolites for most chemicals, predictive tools were used to generate potential metabolites as suspect screening lists to guide the identification of selected xenobiotic substances and their associated metabolites. Thirty-three substances were selected to represent a diverse array of pharmaceutical, agrochemical, and industrial chemicals from Environmental Protection Agency's ToxCast chemical library. The compounds were incubated in a metabolically-active in vitro assay using primary hepatocytes and the resulting supernatant and lysate fractions were analyzed with high-resolution LCMS. Metabolites were simulated for each compound structure using software and then combined to serve as the suspect screening list. The exact masses of the predicted metabolites were then used to select LCMS features for fragmentation via tandem mass spectrometry (MS/MS). Of the starting chemicals, 12 were measured in at least one sample in either positive or negative ion mode and a subset of these were used to develop the analysis workflow. We implemented a screening level workflow for background subtraction and the incorporation of time-varying kinetics into the identification of likely metabolites. We used haloperidol as a case study to perform an in-depth analysis, which resulted in identifying five known metabolites and five molecular features that represent potential novel metabolites, two of which were assigned discrete structures based on in silico predictions. This workflow was applied to five additional test chemicals, and 15 molecular features were selected as either reported metabolites, predicted metabolites, or potential metabolites without a structural assignment. This study demonstrates that in some-but not all-cases, suspect screening analysis methods provide a means to rapidly identify and characterize metabolites of xenobiotic chemicals.

In vitro and in silico study on consequences of combined exposure to the food-borne alkenylbenzenes estragole and safrole

Potential consequences of combined exposure to the selected food-borne alkenylbenzenes safrole and estragole or their proximate carcinogenic 1'-hydroxy metabolites were evaluated in vitro and in silico. HepG2 cells were exposed to 1'-hydroxyestragole and 1'-hydroxysafrole individually or in equipotent combination subsequently detecting cytotoxicity and DNA adduct formation. Results indicate that concentration addition adequately describes the cytotoxic effects and no statistically significant differences were shown in the level of formation of the major DNA adducts. Furthermore, physiologically based kinetic modeling revealed that at normal dietary intake the concentration of the parent compounds and their 1'-hydroxymetabolites remain substantially below the Km values for the respective bioactivation and detoxification reactions providing further support for the fact that the simultaneous presence of the two carcinogens or of their proximate carcinogenic 1'-hydroxy metabolites may not affect their DNA adduct formation. Overall, these results point at the absence of interactions upon combined exposure to selected food-borne alkenylbenzenes at realistic dietary levels of intake.

A Physiological-Based Model for Simulating the Bioavailability and Kinetics of Sulforaphane from Broccoli Products

There are no known physiological-based digestion models that depict glucoraphanin (GR) to sulforaphane (SR) conversion and subsequent absorption. The aim of this research was to make a physiological-based digestion model that includes SR formation, both by endogenous myrosinase and gut bacterial enzymes, and to simulate the SR bioavailability. An 18-compartment model (mouth, two stomach, seven small intestine, seven large intestine, and blood compartments) describing transit, reactions and absorption was made. The model, consisting of differential equations, was fit to data from a human intervention study using Mathwork’s Simulink and Matlab software. SR urine metabolite data from participants who consumed different broccoli products were used to estimate several model parameters and validate the model. The products had high, medium, low, and zero myrosinase content. The model’s predicted values fit the experimental values very well. Parity plots showed that the predicted values closely matched experimental values for the high (r² = 0.95), and low (r² = 0.93) products, but less so for the medium (r² = 0.85) and zero (r² = 0.78) myrosinase products. This is the first physiological-based model to depict the unique bioconversion processes of bioactive SR from broccoli. This model represents a preliminary step in creating a predictive model for the biological effect of SR, which can be used in the growing field of personalized nutrition.

Alkenylbenzenes in Foods: Aspects Impeding the Evaluation of Adverse Health Effects

Alkenylbenzenes are naturally occurring secondary plant metabolites, primarily present in different herbs and spices, such as basil or fennel seeds. Thus, alkenylbenzenes, such as safrole, methyleugenol, and estragole, can be found in different foods, whenever these herbs and spices (or extracts thereof) are used for food production. In particular, essential oils or other food products derived from the aforementioned herbs and spices, such as basil-containing pesto or plant food supplements, are often characterized by a high content of alkenylbenzenes. While safrole or methyleugenol are known to be genotoxic and carcinogenic, the toxicological relevance of other alkenylbenzenes (e.g., apiol) regarding human health remains widely unclear. In this review, we will briefly summarize and discuss the current knowledge and the uncertainties impeding a conclusive evaluation of adverse effects to human health possibly resulting from consumption of foods containing alkenylbenzenes, especially focusing on the genotoxic compounds, safrole, methyleugenol, and estragole.

Assessment of the predictive capacity of a physiologically based kinetic model using a read-across approach

•

Potential regulatory application of PBK modelling information to assist read-across.

•

Presents workflow to read across PBK model information from data-rich to data-poor chemicals.

•

Describes appropriate analogue selection based on a set of specific criteria.

•

Uses estragole and safrole as source chemicals for a target chemical - methyleugenol.

•

Example of PBK model validation where in vivo kinetic data are lacking.

With current progress in science, there is growing interest in developing and applying Physiologically Based Kinetic (PBK) models in chemical risk assessment, as knowledge of internal exposure to chemicals is critical to understanding potential effects in vivo. In particular, a new generation of PBK models is being developed in which the model parameters are derived from in silico and in vitro methods. To increase the acceptance and use of these “Next Generation PBK models”, there is a need to demonstrate their validity. However, this is challenging in the case of data-poor chemicals that are lacking in kinetic data and for which predictive capacity cannot, therefore, be assessed. The aim of this work is to lay down the fundamental steps in using a read across framework to inform modellers and risk assessors on how to develop, or evaluate, PBK models for chemicals without in vivo kinetic data. The application of a PBK model that takes into account the absorption, distribution, metabolism and excretion characteristics of the chemical reduces the uncertainties in the biokinetics and biotransformation of the chemical of interest. A strategic flow-charting application, proposed herein, allows users to identify the minimum information to perform a read-across from a data-rich chemical to its data-poor analogue(s). The workflow analysis is illustrated by means of a real case study using the alkenylbenzene class of chemicals, showing the reliability and potential of this approach. It was demonstrated that a consistent quantitative relationship between model simulations could be achieved using models for estragole and safrole (source chemicals) when applied to methyleugenol (target chemical). When the PBK model code for the source chemicals was adapted to utilise input values relevant to the target chemical, simulation was consistent between the models. The resulting PBK model for methyleugenol was further evaluated by comparing the results to an existing, published model for methyleugenol, providing further evidence that the approach was successful. This can be considered as a “read-across” approach, enabling a valid PBK model to be derived to aid the assessment of a data poor chemical.

Estragole DNA adduct accumulation in human liver HepaRG cells upon repeated in vitro exposure

Accumulation of N²-(trans-isoestragol-3'-yl)-2'-deoxyguanosine (E-3'-N²-dG) DNA adducts derived from the alkenylbenzene estragole upon repeated dose exposure was investigated since the repair of this adduct was previously shown to be inefficient. To this end human HepaRG cells were exposed to repeating cycles of 2 h exposure to 50 μM estragole followed by 22 h repair to mimic daily exposure. The E-3'-N²-dG DNA adduct levels were quantified by LC-MS/MS after each cycle. The results show accumulation of E-3'-N²-dG DNA adducts at a rate of 17.53 adducts/10⁸ nts/cycle. This rate at the dose level calculated by physiologically based kinetic (PBK) modeling to result in 50 μM was converted to a rate expected at average human daily intake of estragole. The predicted time estimated to reach adduct levels reported at the BMD10 of the related alkenylbenzene methyleugenol of 10-100 adducts /10⁸ nts upon average human daily intake of estragole amounted to 8-80 (in rat) or 6-57 years (in human). It is concluded that the persistent nature of the E-3'-N²-dG DNA adducts may contribute to accumulation of substantial levels of DNA adducts upon prolonged dietary exposure.

Molecular Dynamics and In Vitro Quantification of Safrole DNA Adducts Reveal DNA Adduct Persistence Due to Limited DNA Distortion Resulting in Inefficient Repair

The formation and repair of N²-(trans-isosafrol-3'-yl)-2'-deoxyguanosine (S-3'-N²-dG) DNA adduct derived from the spice and herbal alkenylbenzene constituent safrole were investigated. DNA adduct formation and repair were studied in vitro and using molecular dynamics (MD) simulations. DNA adduct formation was quantified using liquid chromatography-mass spectrometry (LCMS) in wild type and NER (nucleotide excision repair) deficient CHO cells and also in HepaRG cells and primary rat hepatocytes after different periods of repair following exposure to safrole or 1'-hydroxysafrole (1'-OH safrole). The slower repair of the DNA adducts found in NER deficient cells compared to that in CHO wild type cells indicates a role for NER in repair of S-3'-N²-dG DNA adducts. However, DNA repair in liver cell models appeared to be limited, with over 90% of the adducts remaining even after 24 or 48 h recovery. In our further studies, MD simulations indicated that S-3'-N²-dG adduct formation causes only subtle changes in the DNA structure, potentially explaining inefficient activation of NER. Inefficiency of NER mediated repair of S-3'-N²-dG adducts points at persistence and potential bioaccumulation of safrole DNA adducts upon daily dietary exposure.

An in vitromodel to quantify interspecies differences in kinetics for intestinal microbial bioactivation and detoxification of zearalenone

Zearalenone (ZEN) is a mycotoxin known for its estrogenic activities. The metabolism of ZEN plays a role in the interspecies differences in sensitivity to ZEN, and is known to occur in the liver and via the intestinal microbiota, although the relative contribution of these two pathways remains to be characterized. In the present study a fecal in vitro model was optimized and used to quantify the interspecies differences in kinetics of the intestinal microbial metabolism of ZEN in rat, pig and human. V_max, K_m, and catalytic efficiencies (k_cat) were determined, and results obtained reveal that the k_cat values for formation of α-ZEL and β-ZEL amounted to 0.73 and 0.12 mL/h/kg bw for human microbiota, 2.6 and 1.3 mL/h/kg bw for rat microbiota and 9.4 and 6.3 mL/h/kg bw for pig microbiota showing that overall ZEN metabolism increased in the order human < rat < pig microbiota. Expressed per kg bw the k_cat for ZEN metabolism by the liver surpassed that of the intestinal microbiota in all three species. In conclusion, it is estimated that the activity of the intestinal colon microbiome may be up to 36 % of the activity of the liver, and that it can additionally contribute to the species differences in bioactivation and detoxification and thus the toxicity of ZEN in pigs and rats but not in humans. The results highlight the importance of the development of human specific models for the assessment of the metabolism of ZEN.

Mode of action-based risk assessment of genotoxic carcinogens

The risk assessment of chemical carcinogens is one major task in toxicology. Even though exposure has been mitigated effectively during the last decades, low levels of carcinogenic substances in food and at the workplace are still present and often not completely avoidable. The distinction between genotoxic and non-genotoxic carcinogens has traditionally been regarded as particularly relevant for risk assessment, with the assumption of the existence of no-effect concentrations (threshold levels) in case of the latter group. In contrast, genotoxic carcinogens, their metabolic precursors and DNA reactive metabolites are considered to represent risk factors at all concentrations since even one or a few DNA lesions may in principle result in mutations and, thus, increase tumour risk. Within the current document, an updated risk evaluation for genotoxic carcinogens is proposed, based on mechanistic knowledge regarding the substance (group) under investigation, and taking into account recent improvements in analytical techniques used to quantify DNA lesions and mutations as well as "omics" approaches. Furthermore, wherever possible and appropriate, special attention is given to the integration of background levels of the same or comparable DNA lesions. Within part A, fundamental considerations highlight the terms hazard and risk with respect to DNA reactivity of genotoxic agents, as compared to non-genotoxic agents. Also, current methodologies used in genetic toxicology as well as in dosimetry of exposure are described. Special focus is given on the elucidation of modes of action (MOA) and on the relation between DNA damage and cancer risk. Part B addresses specific examples of genotoxic carcinogens, including those humans are exposed to exogenously and endogenously, such as formaldehyde, acetaldehyde and the corresponding alcohols as well as some alkylating agents, ethylene oxide, and acrylamide, but also examples resulting from exogenous sources like aflatoxin B1, allylalkoxybenzenes, 2-amino-3,8-dimethylimidazo[4,5-f] quinoxaline (MeIQx), benzo[a]pyrene and pyrrolizidine alkaloids. Additionally, special attention is given to some carcinogenic metal compounds, which are considered indirect genotoxins, by accelerating mutagenicity via interactions with the cellular response to DNA damage even at low exposure conditions. Part C finally encompasses conclusions and perspectives, suggesting a refined strategy for the assessment of the carcinogenic risk associated with an exposure to genotoxic compounds and addressing research needs.

The safety evaluation of food flavoring substances: the role of genotoxicity studies

The Flavor and Extract Manufacturers Association (FEMA) Expert Panel relies on the weight of evidence from all available data in the safety evaluation of flavoring substances. This process includes data from genotoxicity studies designed to assess the potential of a chemical agent to react with DNA or otherwise cause changes to DNA, either in vitro or in vivo. The Panel has reviewed a large number of in vitro and in vivo genotoxicity studies during the course of its ongoing safety evaluations of flavorings. The adherence of genotoxicity studies to standardized protocols and guidelines, the biological relevance of the results from those studies, and the human relevance of these studies are all important considerations in assessing whether the results raise specific concerns for genotoxic potential. The Panel evaluates genotoxicity studies not only for evidence of genotoxicity hazard, but also for the probability of risk to the consumer in the context of exposure from their use as flavoring substances. The majority of flavoring substances have given no indication of genotoxic potential in studies evaluated by the FEMA Expert Panel. Examples illustrating the assessment of genotoxicity data for flavoring substances and the consideration of the factors noted above are provided. The weight of evidence approach adopted by the FEMA Expert Panel leads to a rational assessment of risk associated with consumer intake of flavoring substances under the conditions of use.

Physiologically based kinetic modelling-facilitated reverse dosimetry to predict in vivo red blood cell acetylcholinesterase inhibition following exposure to chlorpyrifos in the Caucasian and Chinese population

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

The safety evaluation of food flavouring substances: the role of metabolic studies

The safety assessment of a flavour substance examines several factors, including metabolic and physiological disposition data. The present article provides an overview of the metabolism and disposition of flavour substances by identifying general applicable principles of metabolism to illustrate how information on metabolic fate is taken into account in their safety evaluation. The metabolism of the majority of flavour substances involves a series both of enzymatic and non-enzymatic biotransformation that often results in products that are more hydrophilic and more readily excretable than their precursors. Flavours can undergo metabolic reactions, such as oxidation, reduction, or hydrolysis that alter a functional group relative to the parent compound. The altered functional group may serve as a reaction site for a subsequent metabolic transformation. Metabolic intermediates undergo conjugation with an endogenous agent such as glucuronic acid, sulphate, glutathione, amino acids, or acetate. Such conjugates are typically readily excreted through the kidneys and liver. This paper summarizes the types of metabolic reactions that have been documented for flavour substances that are added to the human food chain, the methodologies available for metabolic studies, and the factors that affect the metabolic fate of a flavour substance.

Natural occurrence of genotoxic and carcinogenic alkenylbenzenes in Indonesian jamu and evaluation of consumer risks

The consumer risks of jamu, Indonesian traditional herbal medicines, was assessed focussing on the presence of alkenylbenzene containing botanical ingredients. Twenty-three out of 25 samples contained alkenylbenzenes at levels ranging from 3.8 to 440 μg/kg, with methyleugenol being the most frequently encountered alkenylbenzene. The estimated daily intake (EDI) resulting from jamu consumption was estimated to amount to 0.2-171 μg/kg bw/day for individual alkenylbenzenes, to 0.9-203 μg/kg bw/day when adding up all alkenylbenzenes detected, and to 0.9-551 μg/kg bw/day when expressed in methyleugenol equivalents using interim relative potency (REP) factors. The margin of exposure (MOE) values obtained were generally <10,000 indicating a priority for risk management when assuming daily consumption during a lifetime. Using Haber's rule it was estimated that two weeks consumption of these jamu only once would not raise a concern (MOE >10,000). However, when considering use for two weeks every year during a lifetime, 5 samples still raise a concern. It is concluded that the consumption of alkenylbenzene containing jamu can be of concern especially when consumed on a daily basis for longer periods of time on a regular basis.

From in vitro to in vivo: Integration of the virtual cell based assay with physiologically based kinetic modelling

Physiologically based kinetic (PBK) models and the virtual cell based assay can be linked to form so called physiologically based dynamic (PBD) models. This study illustrates the development and application of a PBK model for prediction of estragole-induced DNA adduct formation and hepatotoxicity in humans. To address the hepatotoxicity, HepaRG cells were used as a surrogate for liver cells, with cell viability being used as the in vitro toxicological endpoint. Information on DNA adduct formation was taken from the literature. Since estragole induced cell damage is not directly caused by the parent compound, but by a reactive metabolite, information on the metabolic pathway was incorporated into the model. In addition, a user-friendly tool was developed by implementing the PBK/D model into a KNIME workflow. This workflow can be used to perform in vitro to in vivo extrapolation and forward as backward dosimetry in support of chemical risk assessment.

Investigating the state of physiologically based kinetic modelling practices and challenges associated with gaining regulatory acceptance of model applications

Physiologically based kinetic (PBK) models are used widely throughout a number of working sectors, including academia and industry, to provide insight into the dosimetry related to observed adverse health effects in humans and other species. Use of these models has increased over the last several decades, especially in conjunction with emerging alternative methods to animal testing, such as in vitro studies and data-driven in silico quantitative-structure-activity-relationship (QSAR) predictions. Experimental information derived from these new approach methods can be used as input for model parameters and allows for increased confidence in models for chemicals that did not have in vivo data for model calibration. Despite significant advancements in good modelling practice (GMP) for model development and evaluation, there remains some reluctance among regulatory agencies to use such models during the risk assessment process. Here, the results of a survey disseminated to the modelling community are presented in order to inform the frequency of use and applications of PBK models in science and regulatory submission. Additionally, the survey was designed to identify a network of investigators involved in PBK modelling and knowledgeable of GMP so that they might be contacted in the future for peer review of PBK models, especially in regards to vetting the models to such a degree as to gain a greater acceptance for regulatory purposes.

•

Physiologically Based kinetic (PBK) models are used widely in academia, industry, and government.

•

Good modelling practice (GMP) for model development and evaluation continues to expand.

•

Further guidance for establishing GMP is called for.

•

There remains some reluctance among regulatory agencies to use PBK models.

•

The next generation of PBK models could be developed using only data from in vitro and in silico methods.

Risk assessment of combined exposure to alkenylbenzenes through consumption of plant food supplements containing parsley and dill

A risk assessment was performed of parsley- and dill-based plant food supplements (PFS) containing apiol and related alkenylbenzenes. First, the levels of the alkenylbenzenes in the PFS and the resulting estimated daily intake (EDI) resulting from use of the PFS were quantified. Since most PFS appeared to contain more than one alkenylbenzene, a combined risk assessment was performed based on equal potency or using a so-called toxic equivalency (TEQ) approach based on toxic equivalency factors (TEFs) for the different alkenylbenzenes. The EDIs resulting from daily PFS consumption amount to 0.74-125 µg kg^-1 bw for the individual alkenylbenzenes, 0.74-160 µg kg^-1 bw for the sum of the alkenylbenzenes, and 0.47-64 µg kg^-1 bw for the sum of alkenylbenzenes when expressed in safrole equivalents. The margins of exposure (MOEs) obtained were generally below 10,000, indicating a priority for risk management if the PFS were to be consumed on a daily basis. Considering short-term use of the PFS, MOEs would increase above 10,000, indicating low priority for risk management. It is concluded that alkenylbenzene intake through consumption of parsley- and dill-based PFS is only of concern when these PFS are used for long periods of time.

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

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

Determination and risk assessment of naturally occurring genotoxic and carcinogenic alkenylbenzenes in basil-containing sauce of pesto

A risk assessment of basil-based pesto sauces containing methyleugenol and related alkenylbenzenes was performed based on their levels detected in a series of pesto sauces available on the Dutch market. The estimated daily intake (EDI) values of alkenylbenzenes as a result of consumption of the different pesto sauces amounted to 1.2–44.3 μg/kg bw for individual alkenylbenzenes, 14.3–43.5 μg/kg bw when adding up the alkenylbenzene levels assuming equal potency, and 17.3–62.9 μg/kg bw when expressed in methyleugenol equivalents using alkenylbenzenes defined toxic equivalency factors (TEF). The margin of exposure approach (MOE), used to evaluate the potential risks, resulted in MOE values that were generally lower than 10000 indicating a priority for risk management when assuming daily consumption. The levels of methyleugenol detected in the pesto sauces would allow consumption of 1.1–29.8, 7.5–208, 15.1–416.5, and 32.4–892.5 g of pesto sauce on a daily basis, once a week, once every two weeks, and once a month, respectively, to achieve MOE values above the 10000 limit indicating low priority for risk management. It is concluded that consumption of pesto sauces would only be of concern if consumed on a daily basis over longer periods of time.

Physiologically based kinetic modeling of the bioactivation of myristicin

The present study describes physiologically based kinetic (PBK) models for the alkenylbenzene myristicin that were developed by extension of the PBK models for the structurally related alkenylbenzene safrole in rat and human. The newly developed myristicin models revealed that the formation of the proximate carcinogenic metabolite 1′-hydroxymyristicin in liver is at most 1.8 fold higher in rat than in human and limited for the ultimate carcinogenic metabolite 1′-sulfoxymyristicin to (2.8–4.0)-fold higher in human. In addition, a comparison was made between the relative importance of bioactivation for myristicin and safrole. Model predictions indicate that for these related compounds, the formation of the 1′-sulfoxy metabolites in rat and human liver is comparable with a difference of <2.2-fold over a wide dose range. The results from this PBK analysis support that risk assessment of myristicin may be based on the BMDL₁₀ derived for safrole of 1.9–5.1 mg/kg bw per day. Using an estimated daily intake of myristicin of 0.0019 mg/kg bw per day resulting from the use of herbs and spices, this results in MOE values for myristicin that amount to 1000–2700, indicating a priority for risk management. The results obtained illustrate that PBK modeling provides insight into possible species differences in the metabolic activation of myristicin. Moreover, they provide an example of how PBK modeling can facilitate a read-across in risk assessment from a compound for which in vivo toxicity studies are available to a related compound for which tumor data are not reported, thus contributing to alternatives in animal testing.

Evaluation of Interindividual Human Variation in Bioactivation and DNA Adduct Formation of Estragole in Liver Predicted by Physiologically Based Kinetic/Dynamic and Monte Carlo Modeling

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.

Dose-dependent DNA adduct formation by cinnamaldehyde and other food-borne α,β-unsaturated aldehydes predicted by physiologically based in silico modelling

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.

Physiologically based pharmacokinetic modeling of impaired carboxylesterase-1 activity: effects on oseltamivir disposition

BACKGROUND AND OBJECTIVE

Human carboxylesterase-1 (CES1) is an enzyme that is primarily expressed in the liver, where it plays an important role in the metabolism of many commonly used medications. Ethanol (alcohol)-mediated inhibition of CES1 and loss-of-function polymorphisms in the CES1 gene can markedly reduce this enzyme's function. Such alterations in CES1 activity may have important effects on the disposition of substrate drugs. The aim of this study is to develop a physiologically based pharmacokinetic (PBPK) model to predict changes in CES1 substrate drug exposure in humans with CES1 activity impaired by ethanol or loss-of-function CES1 genetic polymorphisms.

METHODS

The antiviral drug oseltamivir, an ethyl ester prodrug that is rapidly converted in vivo to the active metabolite oseltamivir carboxylate (OSC) by CES1 was used as a probe drug for CES1 activity. Oseltamivir PBPK models integrating in vitro and in vivo data were developed and refined. Then the changes in oseltamivir and OSC exposure in humans with CES1 impaired by ethanol or polymorphisms were simulated using a PBPK model incorporating in vitro inhibition and enzyme kinetic data. Model assumptions were verified by comparison of simulations with observed and published data. A sensitivity analysis was performed to gain a mechanistic understanding of the exposure changes of oseltamivir and OSC.

RESULTS

The simulated changes in oseltamivir and OSC exposures in humans with CES1 impaired by ethanol or polymorphism were similar to the observed data. The observed exposures to oseltamivir were increased by 46 and 37 % for the area under the plasma concentration-time curve from time zero to 6 h (AUC6) and from time zero to 24 h (AUC24), respectively, with co-administration of ethanol 0.6 g/kg. In contrast, only a slight change was observed in OSC exposure. The simulated data show the same trend as evidenced by greater change in exposures to oseltamivir (27 and 26 % for AUC(6) and AUC(24), [corrected] respectively) than OSC (≤6 %).

CONCLUSIONS

The PBPK model of impaired CES1 activity correctly predicts observed human data. This model can be extended to predict the effects of drug interactions and other factors affecting the pharmacokinetics of other CES1 substrate drugs.

Matrix-derived combination effects influencing absorption, distribution, metabolism and excretion (ADME) of food-borne toxic compounds: implications for 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.

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

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

An integrated QSAR-PBK/D modelling approach for predicting detoxification and DNA adduct formation of 18 acyclic food-borne α,β-unsaturated aldehydes

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.

Investigation of the cytotoxic, genotoxic, and apoptosis-inducing effects of estragole isolated from fennel (Foeniculum vulgare)

The present study was undertaken to evaluate, in the HepG2 human hepatoma cell line, the in vitro cytotoxic, genotoxic, and apoptotic activities of estragole (1), contained in the essential oil of Foeniculum vulgare (fennel) and suspected to induce hepatic tumors in susceptible strains of mice. Toward this end, an MTT cytotoxicity assay, a trypan blue dye exclusion test, a double-staining (acridine orange and DAPI) fluorescence viability assay, a single-cell microgel-electrophoresis (comet) assay, a mitochondrial membrane potential (Δψm) assay, and a DNA fragmentation analysis were conducted. In terms of potential genotoxic effects, the comet assay indicated that estragole (1) was not able to induce DNA damage nor apoptosis under the experimental conditions used.

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

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

Fresh aromatic herbs containing methylchavicol did not exhibit the pro-oxidative effects of pure methylchavicol on a human hepatoma cell line, HepG2

Methylchavicol (CH(3)-CV), an important aromatic constituent of different plants like tarragon and basils, has been shown to be carcinogenic by a mechanism yet unclear, although it has been reported that carcinogenicity of CH(3)-CV in rodent might be linked to its metabolic conversion into a genotoxic electrophilic metabolite generated through a two steps bioactivation pathway catalyzed by cytochrome P450 enzymes and sulfotransferases. The induction of carcinogenesis by certain agents has been associated with the generation of oxidative stress. The aim of the present study was to determine whether pure methylchavicol applied on a human hepatoma cell line, HepG2, could promote oxidative stress and might alter the expression of procarcinogenic biomarkers such as the drug-metabolizing enzyme (CYP2E1), the inducible form of nitric oxide synthase (iNOS) and might induce the expression of Cu/Zn-superoxide dismutase (Cu/Zn-SOD) and Mn-SOD that control the redox equilibrium of the cells. CH(3)-CV was shown to cause a significant induction of oxidative stress, as revealed by luminol-dependent chemiluminescence (LDCL) and to alter dramatically the expression of CYP2E1, iNOS and Mn-SOD, indicating that the toxic effect of CH(3)-CV could be mediated through a nitric oxide dependent mechanism. Under similar experimental conditions, the extracts from tarragon, chervil and basil did not induce such biological changes. These results provide evidence that the generation of an oxidative stress may be a significant event occurring during CH(3)-CV-induced toxicity. It also suggests that natural extracts containing different amounts of CH(3)-CV (tarragon, chervil and basil) did not elicit such toxicity and might contain compounds able to counteract this detrimental property.

A physiologically based in silico kinetic model predicting plasma cholesterol concentrations in humans

Increased plasma cholesterol concentration is associated with increased risk of cardiovascular disease. This study describes the development, validation, and analysis of a physiologically based kinetic (PBK) model for the prediction of plasma cholesterol concentrations in humans. This model was directly adapted from a PBK model for mice by incorporation of the reaction catalyzed by cholesterol ester transfer protein and contained 21 biochemical reactions and eight different cholesterol pools. The model was calibrated using published data for humans and validated by comparing model predictions on plasma cholesterol levels of subjects with 10 different genetic mutations (including familial hypercholesterolemia and Smith-Lemli-Opitz syndrome) with experimental data. Average model predictions on total cholesterol were accurate within 36% of the experimental data, which was within the experimental margin. Sensitivity analysis of the model indicated that the HDL cholesterol (HDL-C) concentration was mainly dependent on hepatic transport of cholesterol to HDL, cholesterol ester transfer from HDL to non-HDL, and hepatic uptake of cholesterol from non-HDL-C. Thus, the presented PBK model is a valid tool to predict the effect of genetic mutations on cholesterol concentrations, opening the way for future studies on the effect of different drugs on cholesterol levels in various subpopulations in silico.

Evaluation of research activities and research needs to increase the impact and applicability of alternative testing strategies in risk assessment practice

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.

Quantitative comparison between in vivo DNA adduct formation from exposure to selected DNA-reactive carcinogens, natural background levels of DNA adduct formation and tumour incidence in rodent bioassays

This study aimed at quantitatively comparing the occurrence/formation of DNA adducts with the carcinogenicity induced by a selection of DNA-reactive genotoxic carcinogens. Contrary to previous efforts, we used a very uniform set of data, limited to in vivo rat liver studies in order to investigate whether a correlation can be obtained, using a benchmark dose (BMD) approach. Dose-response data on both carcinogenicity and in vivo DNA adduct formation were available for six compounds, i.e. 2-acetylaminofluorene, aflatoxin B1, methyleugenol, safrole, 2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline and tamoxifen. BMD(10) values for liver carcinogenicity were calculated using the US Environmental Protection Agency BMD software. DNA adduct levels at this dose were extrapolated assuming linearity of the DNA adduct dose response. In addition, the levels of DNA adducts at the BMD(10) were compared to available data on endogenous background DNA damage in the target organ. Although for an individual carcinogen the tumour response increases when adduct levels increase, our results demonstrate that when comparing different carcinogens, no quantitative correlation exists between the level of DNA adduct formation and carcinogenicity. These data confirm that the quantity of DNA adducts formed by a DNA-reactive compound is not a carcinogenicity predictor but that other factors such as type of adduct and mutagenic potential may be equally relevant. Moreover, comparison to background DNA damage supports the notion that the mere occurrence of DNA adducts above or below the level of endogenous DNA damage is neither correlated to development of cancer. These data strongly emphasise the need to apply the mode of action framework to understand the contribution of other biological effect markers playing a role in carcinogenicity.