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

Development of an allele-specific PCR (AS-PCR) method for identifying high-methyl eugenol-containing purple Tulsi (Ocimum tenuiflorum L.) in market samples

BACKGROUND

Ocimum tenuiflorum L. is a highly traded medicinal with several therapeutic values. Green Tulsi and purple Tulsi are two subtypes in O. tenuiflorum and both have the same medicinal properties. Recent reports have revealed that purple Tulsi contains higher quantities of methyl eugenol (ME), which is moderately toxic and potentially carcinogenic. Therefore, we developed an allele-specific PCR (AS-PCR) method to distinguish the green and purple Tulsi.

METHODS AND RESULT

Using the green Tulsi as a reference, 12 single nucleotide polymorphisms (SNPs) and 10 insertions/deletions (InDels) were identified in the chloroplast genome of the purple Tulsi. The C > T SNP at the 1,26,029 position in the ycf1 gene was selected for the development of the AS-PCR method. The primers were designed to amplify 521 bp and 291 bp fragments specific to green and purple Tulsi, respectively. This AS-PCR method was validated in 10 accessions from each subtype and subsequently verified using Sanger sequencing. Subsequently, 30 Tulsi powder samples collected from the market were subjected to molecular identification by AS-PCR. The results showed that 80% of the samples were purple Tulsi, and only 3.5% were green Tulsi. About 10% of the samples were a mixture of both green and purple Tulsi. Two samples (6.5%) did not contain O. tenuiflorum and were identified as O. gratissimum.

CONCLUSION

The market samples of Tulsi were predominantly derived from purple Tulsi. The AS-PCR method will be helpful for quality control and market surveillance of Tulsi herbal powders.

Urinary Excretion of Mercapturic Acids of the Rodent Carcinogen Methyleugenol after a Single Meal of Basil Pesto: A Controlled Exposure Study in Humans

Methyleugenol (ME), found in numerous plants and spices, is a rodent carcinogen and is classified as "possibly carcinogenic to humans". The hypothesis of a carcinogenic risk for humans is supported by the observation of ME-derived DNA adducts in almost all human liver and lung samples examined. Therefore, a risk assessment of ME is needed. Unfortunately, biomarkers of exposure for epidemiological studies are not yet available. We hereby present the first detection of N-acetyl-l-cysteine conjugates (mercapturic acids) of ME in human urine samples after consumption of a popular ME-containing meal, pasta with basil pesto. We synthesized mercapturic acid conjugates of ME, identified the major product as N-acetyl-S-[3'-(3,4-dimethoxyphenyl)allyl]-l-cysteine (E-3'-MEMA), and developed methods for its extraction and LC-MS/MS quantification in human urine. For conducting an exposure study in humans, a basil cultivar with a suitable ME content was grown for the preparation of basil pesto. A defined meal containing 100 g of basil pesto, corresponding to 1.7 mg ME, was served to 12 participants, who collected the complete urine at defined time intervals for 48 h. Using d6-E-3'-MEMA as an internal standard for LC-MS/MS quantification, we were able to detect E-3'-MEMA in urine samples of all participants collected after the ME-containing meal. Excretion was maximal between 2 and 6 h after the meal and was completed within about 12 h (concentrations below the limit of detection). Excreted amounts were only between 1 and 85 ppm of the ME intake, indicating that the ultimate genotoxicant, 1'-sulfooxy-ME, is formed to a subordinate extent or is not efficiently detoxified by glutathione conjugation and subsequent conversion to mercapturic acids. Both explanations may apply cumulatively, with the ubiquitous detection of ME DNA adducts in human lung and liver specimens arguing against an extremely low formation of 1'-sulfooxy-ME. Taken together, we hereby present the first noninvasive human biomarker reflecting an internal exposure toward reactive ME species.

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.

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.

Human Family 1-4 cytochrome P450 enzymes involved in the metabolic activation of xenobiotic and physiological chemicals: an update

This is an overview of the metabolic activation of drugs, natural products, physiological compounds, and general chemicals by the catalytic activity of cytochrome P450 enzymes belonging to Families 1-4. The data were collected from > 5152 references. The total number of data entries of reactions catalyzed by P450s Families 1-4 was 7696 of which 1121 (~ 15%) were defined as bioactivation reactions of different degrees. The data were divided into groups of General Chemicals, Drugs, Natural Products, and Physiological Compounds, presented in tabular form. The metabolism and bioactivation of selected examples of each group are discussed. In most of the cases, the metabolites are directly toxic chemicals reacting with cell macromolecules, but in some cases the metabolites formed are not direct toxicants but participate as substrates in succeeding metabolic reactions (e.g., conjugation reactions), the products of which are final toxicants. We identified a high level of activation for three groups of compounds (General Chemicals, Drugs, and Natural Products) yielding activated metabolites and the generally low participation of Physiological Compounds in bioactivation reactions. In the group of General Chemicals, P450 enzymes 1A1, 1A2, and 1B1 dominate in the formation of activated metabolites. Drugs are mostly activated by the enzyme P450 3A4, and Natural Products by P450s 1A2, 2E1, and 3A4. Physiological Compounds showed no clearly dominant enzyme, but the highest numbers of activations are attributed to P450 1A, 1B1, and 3A enzymes. The results thus show, perhaps not surprisingly, that Physiological Compounds are infrequent substrates in bioactivation reactions catalyzed by P450 enzyme Families 1-4, with the exception of estrogens and arachidonic acid. The results thus provide information on the enzymes that activate specific groups of chemicals to toxic metabolites.

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

SCOPE

To predict gut microbial metabolism of xenobiotics and the resulting plasma concentrations of metabolites formed, an in vitro-in silico-based testing strategy is developed using the isoflavone daidzein and its gut microbial metabolite S-equol as model compounds.

METHODS AND RESULTS

Anaerobic rat fecal incubations are optimized and performed to derive the apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for gut microbial conversion of daidzein to dihydrodaidzein, S-equol, and O-desmethylangolensin, which are input as parameters for a physiologically based kinetic (PBK) model. The inclusion of gut microbiota in the PBK model allows prediction of S-equol concentrations and slightly reduced predicted maximal daidzein concentrations from 2.19 to 2.16 µm. The resulting predicted concentrations of daidzein and S-equol are comparable to in vivo concentrations reported.

CONCLUSION

The optimized in vitro approach to quantify kinetics for gut microbial conversions, and the newly developed PBK model for rats that includes gut microbial metabolism, provide a unique tool to predict the in vivo consequences of daidzein microbial metabolism for systemic exposure of the host to daidzein and its metabolite S-equol. The predictions reveal a dominant role for daidzein in ERα-mediated estrogenicity despite the higher estrogenic potency of its microbial metabolite S-equol.

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

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

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

Methyleugenol DNA adducts in human liver are associated with SULT1A1 copy number variations and expression levels

Methyleugenol is a rodent hepatocarcinogen occurring in many herbs and spices as well as essential oils used for flavoring. Following metabolic activation by cytochromes P450 (CYPs) and sulfotransferases (SULTs), methyleugenol can form DNA adducts. Previously, we showed that DNA adduct formation by methyleugenol in mouse liver is dependent on SULT1A1 expression and that methyleugenol DNA adducts are abundant in human liver specimens. In humans, SULT1A1 activity is affected by genetic polymorphisms, including single-nucleotide polymorphisms (SNPs) and copy number variations (CNVs). Here we investigated the relationship between individual methyleugenol DNA adduct levels and SULT1A1 in human liver samples. Using isotope-dilution ultraperformance liquid chromatography coupled with tandem mass spectrometry, we quantified methyleugenol DNA adducts in 121 human surgical liver samples. Frequent CNVs, including deletions (f = 3.3%) and duplications (f = 36.4%) of SULT1A1, were identified using qPCR and TaqMan assays in the donors' genomic DNA. SULT1A1 mRNA and protein levels were quantified using microarray data and Western blot analysis, respectively. Methyleugenol DNA adducts were detected in all 121 liver samples studied. Their levels varied 122-fold between individuals and were significantly correlated to both mRNA and protein levels of SULT1A1 (r _s = 0.43, and r _s = 0.44, respectively). Univariate and multivariate statistical analysis identified significant associations of SULT1A1 CNVs with mRNA (p = 1.7 × 10^-06) and protein (p = 4.4 × 10^- 10) levels as well as methyleugenol DNA adduct levels (p = 0.003). These data establish the importance of SULT1A1 genotype for hepatic methyleugenol DNA adducts in humans, and they confirm a strong impact of SULT1A1 CNVs on SULT1A1 hepatic phenotype.

Use of Physiologically Based Kinetic Modeling-Based Reverse Dosimetry to Predict in Vivo Toxicity from in Vitro Data

The development of reliable nonanimal based testing strategies, such as in vitro bioassays, is the holy grail in current human safety testing of chemicals. However, the use of in vitro toxicity data in risk assessment is not straightforward. One of the main issues is that concentration-response curves from in vitro models need to be converted to in vivo dose-response curves. These dose-response curves are needed in toxicological risk assessment to obtain a point of departure to determine safe exposure levels for humans. Recent scientific developments enable this translation of in vitro concentration-response curves to in vivo dose-response curves using physiologically based kinetic (PBK) modeling-based reverse dosimetry. The present review provides an overview of the examples available in the literature on the prediction of in vivo toxicity using PBK modeling-based reverse dosimetry of in vitro toxicity data, showing that proofs-of-principle are available for toxicity end points ranging from developmental toxicity, nephrotoxicity, hepatotoxicity, and neurotoxicity to DNA adduct formation. This review also discusses the promises and pitfalls, and the future perspectives of the approach. Since proofs-of-principle available so far have been provided for the prediction of toxicity in experimental animals, future research should focus on the use of in vitro toxicity data obtained in human models to predict the human situation using human PBK models. This would facilitate human- instead of experimental animal-based approaches in risk assessment.

The influence of the SULT1A status - wild-type, knockout or humanized - on the DNA adduct formation by methyleugenol in extrahepatic tissues of mice

Methyleugenol, present in herbs and spices, has demonstrated carcinogenic activity in the liver and, to a lesser extent, in extrahepatic tissues of rats and mice. It forms DNA adducts after hydroxylation and sulphation. As previously reported, hepatic DNA adduct formation by methyleugenol in mice is strongly affected by their sulphotransferase (SULT) 1A status. Now, we analysed the adduct formation in extrahepatic tissues. The time course of the adduct levels was determined in transgenic (tg) mice, expressing human SULT1A1/2, after oral administration of methyleugenol (50 mg per kg body mass). Nearly maximal adduct levels were observed 6 h after treatment. They followed the order: liver > caecum > kidney > colon > stomach > small intestine > lung > spleen. We then selected liver, caecum, kidney and stomach for the main study, in which four mouse lines [wild-type (wt), Sult1a1-knockout (ko), tg, and humanized (ko-tg)] were treated with methyleugenol at varying dose levels. In the liver, caecum and kidney, adduct formation was nearly completely dependent on the expression of SULT1A enzymes. In the liver, human SULT1A1/2 led to higher adduct levels than mouse Sult1a1, and the effects of both enzymes were approximately additive. In the caecum, human SULT1A1/2 and mouse Sult1a1 were nearly equally effective, again with additive effects in tg mice. In the kidney, only human SULT1A1/2 played a role: no adducts were detected in wt and ko mice even at the highest dose tested and the adduct levels were similar in tg and ko-tg mice. In the stomach, adduct formation was unaffected by the SULT1A status.

IN CONCLUSION

(i) the SULT1A enzymes only affected adduct formation in those tissues in which they are highly expressed (mouse Sult1a1 in the liver and caecum, but not in the kidney and stomach; human SULT1A1/2 in the liver, caecum and kidney, not in the stomach of tg mice and humans), indicating a dominating role of local bioactivation; (ii) the additivity of the effects of both enzymes in the liver and caecum implies that the enzyme level was limiting in the adduct formation; (iii) SULT1A forms dominated the activation of methyleugenol in several tissues, but non-Sult1a1 forms or SULT-independent mechanisms were involved in its adduct formation in the stomach.

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.