Next generation risk assessment for occupational chemical safety - A real world example with sodium-2-hydroxyethane sulfonate.

Next Generation Risk Assessment (NGRA) is an exposure-led approach to safety assessment that uses New Approach Methodologies (NAMs). Application of NGRA has been largely restricted to assessments of consumer use of cosmetics and is not currently implemented in occupational safety assessments, e.g. under EU REACH. By contrast, a large proportion of regulatory worker safety assessments are underpinned by toxicological studies using experimental animals. Consequently, occupational safety assessment represents an area that would benefit from increasing application of NGRA to safety decision making. Here, a workflow for conducting NGRA under an occupational safety context was developed, which is illustrated with a case study chemical; sodium 2-hydroxyethane sulphonate (sodium isethionate or SI). Exposures were estimated using a standard occupational exposure model following a comprehensive life cycle assessment of SI and considering factory-specific data. Outputs of this model were then used to estimate internal exposures using a Physiologically Based Kinetic (PBK) model, which was constructed with SI specific Absorption, Distribution, Metabolism and Excretion (ADME) data. PBK modelling indicated a worst-case plasma maximum concentration (Cmax) of 0.8 µM across the SI life cycle. SI bioactivity was assessed in a battery of NAMs relevant to systemic, reproductive, and developmental toxicity; a cell stress panel, high throughput transcriptomics in three cell lines (HepG2, HepaRG and MCF-7 cells), pharmacological profiling and specific assays relating to developmental toxicity (Reprotracker and devTOX quickPredict). Points of Departure (PoDs) for SI ranged from 104 to 5044 µM. Cmax values obtained from PBK modelling of occupational exposures to SI were compared with PoDs from the bioactivity assays to derive Bioactivity Exposure Ratios (BERs) which demonstrated the safety for workers exposed to SI under current levels of factory specific risk management. In summary, the tiered and iterative workflow developed here represents an opportunity for integrating non animal approaches for a large subset of substances for which systemic worker safety assessment is required. Such an approach could be followed to ensure that animal testing is only conducted as a "last resort" e.g. under EU REACH.

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

Mode of action based risk assessment of the botanical food-borne alkenylbenzene apiol from parsley using physiologically based kinetic (PBK) modelling and read-across from safrole.

The present study developed physiologically-based kinetic (PBK) models for the alkenylbenzene apiol in order to facilitate risk assessment based on read-across from the related alkenylbenzene safrole. Model predictions indicate that in rat liver the formation of the 1'-sulfoxy metabolite is about 3 times lower for apiol than for safrole. These data support that the lower confidence limit of the benchmark dose resulting in a 10% extra cancer incidence (BMDL10) that would be obtained in a rodent carcinogenicity study with apiol may be 3-fold higher for apiol than for safrole. These results enable a preliminary risk assessment for apiol, for which tumor data are not available, using a BMDL10 value of 3 times the BMDL10 for safrole. Based on an estimated BMDL10 for apiol of 5.7-15.3 mg/kg body wt per day and an estimated daily intake of 4 × 10(-5) mg/kg body wt per day, the margin of exposure (MOE) would amount to 140,000-385,000. This indicates a low priority for risk management. The present study shows how PBK modelling can contribute to the development of alternatives for animal testing, facilitating read-across from compounds for which in vivo toxicity studies on tumor formation are available to compounds for which these data are unavailable.

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.

A physiologically based in silico model for trans-2-hexenal detoxification and DNA adduct formation in rat.

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.