Human variability in xenobiotic metabolism and pathway-related uncertainty factors for chemical risk assessment: a review

Development of an LC-HRMS non-targeted method for comprehensive profiling of the exposome of nicotine and tobacco product users - A showcase for cigarette smokers

The global prevalence of electronic cigarettes, heated tobacco products, and other smokeless alternatives has grown significantly in the last ten years. These products have been suggested as combustion-free alternatives for conventional tobacco products like cigarettes, aiming to reduce the negative health impacts associated with smoking. However, the impact of those products on the health and safety of the general population are still unclear, as the absolute exposure from those products has not been thoroughly studied, yet. In this project, a non-targeted LC-HRMS method was developed comprising four different analytical modes for the investigation of the exposure profile in urine of the product users. The method is characterized by its high sensitivity and reproducibility, as shown during method validation. As a proof of concept, we first applied this method to detect significant differences in biomarkers of exposure (BoEs) between smokers and non-smokers. We observed a total of 171 BoEs significantly elevated in smokers, including several well-known biomarkers of smoke exposure like nicotine and its metabolites, mercapturic acid derivatives, and phenolic compounds. Some of the detected biomarkers are present at low ng/mL concentrations in urine, proving the high sensitivity needed for a holistic exploration of the exposome. Moreover, we were able to identify BoEs that have not been reported previously for smoking, such as 2,6-dimethoxyphenol and 7-methyl-1-naphthol glucuronide.

Cytochrome P450 isoforms contribution, plasma protein binding, toxicokinetics of enniatin A in rats and in vivo clearance prediction in humans

Emerging mycotoxins, such as enniatin A (ENNA), are becoming a worldwide concern owing to their presence in different types of food and feed. However, comprehensive toxicokinetic data that links intake, exposure and toxicological effects of ENNA has not been elucidated yet. Therefore, the present study investigated the in vitro (rat and human) and in vivo (rat) toxicokinetic properties of ENNA. Towards this, an easily applicable and sensitive bioanalytical method was developed and validated for the estimation of ENNA in rat plasma. ENNA exhibited high plasma protein binding (99%), high hepatic clearance and mainly underwent metabolism via CYP3A4 (74%). The in-house predicted hepatic clearance (54 mL/min/kg) and observed in vivo rat clearance (55 mL/min/kg) were comparable. The predicted in vivo human hepatic clearance was 18 mL/min/kg. ENNA underwent slow absorption (T_max = 4 h) and rapid elimination following oral administration to rats. The absolute oral bioavailability was 47%. The toxicokinetic findings for ENNA from this study will help in designing and interpreting toxicological studies in rats. Besides, these findings could be used in physiologically based toxicokinetic (PBTK) model development for exposure predictions and risk assessment for ENNA in humans.

Distributions for time, inter‐ and intraspecies extrapolation for deriving occupational exposure limits

This work aimed at improving the empirical database of time (i.e., exposure duration), interspecies and intraspecies extrapolation when deriving occupational exposure limits (OELs). For each extrapolation step, a distribution was derived, which can be used to model the associated uncertainties. For time and interspecies extrapolation, distributions of ratios of dose descriptors were derived from studies of different length or species. National Toxicology Program (NTP) study data were manually assessed, and data from REACH (Registration, Evaluation and Authorisation of Chemicals) registration dossiers were evaluated semi‐automatically. Intraspecies extrapolation was investigated by compiling published studies on human toxicokinetic and toxicodynamic variability. A new database was established for toxicokinetic differences in interindividual susceptibility, including many inhalation studies. Using NTP data produced more reliable results than using REACH data. The geometric mean (GM) for time extrapolation subacute/chronic agreed with previous evaluations (GM = 4.11), whereas the GM for subchronic/chronic extrapolation was slightly higher (GM = 2.93) than the GMs found by others. No significant differences were observed between systemically and locally acting substances. Observed interspecies differences confirmed the suitability of allometric scaling, with the derived distribution describing remaining uncertainty. Distributions of intraspecies variability at the 1% and 5% incidence level had medians of 7.25 and 3.56, respectively. When compared with assessment factors (AFs) currently used in the EU, probabilities that these AFs are protective enough span a wide range from 10% to 95%, depending on the extrapolation step. These results help to select AFs in a transparent and informed way and, by allowing to compare protection levels achieved, to harmonise methods for deriving OELs.

This work aimed at improving the empirical database of time (i.e., exposure duration), interspecies and intraspecies extrapolation. Distributions were derived, which can be used to model the associated uncertainties. When compared with assessment factors (AFs) currently used in the EU, probabilities that these AFs are protective enough span a wide range from 10% to 95%, depending on the extrapolation step. These results help to select AFs in a transparent and informed way and to harmonise methods for deriving OELs.

Skin Barrier Function in Infants: Update and Outlook

A good understanding of infant skin should provide a rationale for optimum management of the health of this integument. In this review, we discuss the skin barrier function of infants, particularly with reference to the use of diapers and baby wipes. The skin barrier of newborns continues to develop with age. Two years after birth, the barrier properties of infant skin closely resemble those of adult skin. However, several risk factors may contribute to impaired skin barrier and altered skin permeability in infants. Problems may arise from the use of diapers and baby wipes. The skin covered by a diaper is effectively an occluded environment, and thus is vulnerable to over-hydration. To date there has been no published information regarding dermal absorption of ingredients contained in baby wipes. Similarly, dermal absorption of topical ingredients in infants with underlying skin conditions has not been widely explored. Clearly, there are serious ethical concerns related to conducting skin permeation studies on infant skin. However, the increasing availability of non-invasive methods for in vivo studies is encouraging and offers new directions for studying this important patient group.

Human variability in polymorphic CYP2D6 metabolism: Implications for the risk assessment of chemicals in food and emerging designer drugs

The major human cytochrome P450 CYP2D6 isoform enzyme plays important roles in the liver and in the brain with regards to xenobiotic metabolism. Xenobiotics as CYP2D6 substrates include a whole range of pharmaceuticals, pesticides and plant alkaloids to cite but a few. In addition, a number of endogenous compounds have been shown to be substrates of CYP2D6 including trace amines in the brain such as tyramine and 5-methoxytryptamine as well as anandamide and progesterone. Because of the polymorphic nature of CYP2D6, considerable inter-phenotypic and inter-ethnic differences in the pharmaco/toxicokinetics (PK/TK) and metabolism of CYP2D6 substrates exist with potential consequences on the pharmacology and toxicity of chemicals. Here, large extensive literature searches have been performed to collect PK data from published human studies for a wide range of pharmaceutical probe substrates and investigate human variability in CYP2D6 metabolism. The computed kinetic parameters resulted in the largest open source database, quantifying inter-phenotypic differences for the kinetics of CYP2D6 probe substrates in Caucasian and Asian populations, to date. The database is available in supplementary material (CYPD6 DB) and EFSA knowledge junction (DOI to added). Subsequently, meta-analyses using a hierarchical Bayesian model for markers of chronic oral exposure (oral clearance, area under the plasma concentration time curve) and acute oral exposure (maximum plasma concentration (Cmax) provided estimates of inter-phenotypic differences and CYP2D6-related uncertainty factors (UFs) for chemical risk assessment in Caucasian and Asian populations classified as ultra-rapid (UM), extensive (EMs), intermediate (IMs) and poor metabolisers (PMs). The model allowed the integration of inter-individual (i.e. inter-phenotypic and inter-ethnic), inter-compound and inter-study variability together with uncertainty in each PK parameter. Key findings include 1. Higher frequencies of PMs in Caucasian populations compared to Asian populations (>8% vs 1-2%) for which EM and IM were the most frequent phenotype. 2. Large inter-phenotypic differences in PK parameters for Caucasian EMs (coefficients of variation (CV) > 50%) compared with Caucasian PMs and Asian EMs and IMs (i.e CV < 40%). 3. Inter-phenotypic PK differences between EMs and PMs in Caucasian populations increase with the quantitative contribution of CYP2D6 for the metabolism (fm) for a range of substrates (fm_CYP2D6 range: 20-95% of dose) (range: 1-54) to a much larger extent than those for Asian populations (range: 1-4). 4. Exponential meta-regressions between Fm_CYP2D6 in EMs and inter-phenotypic differences were also shown to differ between Caucasian and Asian populations as well as CYP2D6-related UFs. Finally, implications of these results for the risk assessment of food chemicals and emerging designer drugs of public health concern, as CYP2D6 substrates, are highlighted and include the integration of in vitro metabolism data and CYP2D6-variability distributions for the development of quantitative in vitro in vivo extrapolation models.

Body fluids from the rat exposed to chlorpyrifos induce cytotoxicity against the corresponding tissue-derived cells in vitro

BACKGROUND

This study aims to establish an in vitro monitoring approach to evaluate the pesticide exposures. We studied the in vitro cytotoxicity of three different body fluids of rats to the respective corresponding tissue-derived cells.

METHODS

Wistar rats were orally administrated daily with three different doses of chlorpyrifos (1.30, 3.26, and 8.15 mg/kg body weight/day, which is equal to the doses of 1/125, 1/50, and 1/20 LD50, respectively) for consecutive 90 days. Blood samples as well as 24-hour urine and fecal samples were collected and processed. Then, urine, serum, and feces samples were used to treat the correspondent cell lines, i.e., T24 bladder cancer cells, Jurkat lymphocytes, and HT-29 colon cancer cells respectively, which derived from the correspondent tissues that could interact with the respective corresponding body fluids in organism. Cell viability was determined by using MTT or trypan blue staining.

RESULTS

The results showed that urine, serum, and feces extract of the rats exposed to chlorpyrifos displayed concentration- and time-dependent cytotoxicity to the cell lines. Furthermore, we found that the cytotoxicity of body fluids from the exposed animals was mainly due to the presence of 3, 4, 5-trichloropyrindinol, the major toxic metabolite of chlorpyrifos.

CONCLUSIONS

These findings indicated that urine, serum, and feces extraction, especially urine, combining with the corresponding tissue-derived cell lines as the in vitro cell models could be used to evaluate the animal exposure to pesticides even at the low dose with no apparent toxicological signs in the animals. Thus, this in vitro approach could be served as complementary methodology to the existing toolbox of biological monitoring of long-term and low-dose exposure to environmental pesticide residues in practice.

Finger sweat analysis enables short interval metabolic biomonitoring in humans

Metabolic biomonitoring in humans is typically based on the sampling of blood, plasma or urine. Although established in the clinical routine, these sampling procedures are often associated with a variety of compliance issues, which are impeding time-course studies. Here, we show that the metabolic profiling of the minute amounts of sweat sampled from fingertips addresses this challenge. Sweat sampling from fingertips is non-invasive, robust and can be accomplished repeatedly by untrained personnel. The sweat matrix represents a rich source for metabolic phenotyping. We confirm the feasibility of short interval sampling of sweat from the fingertips in time-course studies involving the consumption of coffee or the ingestion of a caffeine capsule after a fasting interval, in which we successfully monitor all known caffeine metabolites as well as endogenous metabolic responses. Fluctuations in the rate of sweat production are accounted for by mathematical modelling to reveal individual rates of caffeine uptake, metabolism and clearance. To conclude, metabotyping using sweat from fingertips combined with mathematical network modelling shows promise for broad applications in precision medicine by enabling the assessment of dynamic metabolic patterns, which may overcome the limitations of purely compositional biomarkers.

The codification of hazard and its impact on the hazard versus risk controversy

The long running controversy about the relative merits of hazard-based versus risk-based approaches has been investigated. There are three levels of hazard codification: level 1 divides chemicals into dichotomous bands of hazardous and non-hazardous; level 2 divides chemicals into bands of hazard based on severity and/or potency; and level 3 places each chemical on a continuum of hazard based on severity and/or potency. Any system which imposes compartments onto a continuum will give rise to issues at the boundaries, especially with only two compartments. Level 1 schemes are only justifiable if there is no variation in severity, or potency or if there is no threshold. This is the assumption implicit in GHS/EU classification for carcinogenicity, reproductive toxicity and mutagenicity. However, this assumption has been challenged. Codification level 2 hazard assessments offer a range of choices and reduce the built-in conflict inherent in the level 1 process. Level 3 assessments allow a full range of choices between the extremes and reduce the built-in conflict even more. The underlying reason for the controversy between hazard and risk is the use of level 1 hazard codification schemes in situations where there are ranges of severity and potency which require the use of level 2 or level 3 hazard codification. There is not a major difference between level 2 and level 3 codification, and they can both be used to select appropriate risk management options. Existing level 1 codification schemes should be reviewed and developed into level 2 schemes where appropriate.

Human variability in isoform-specific UDP-glucuronosyltransferases: markers of acute and chronic exposure, polymorphisms and uncertainty factors

UDP-glucuronosyltransferases (UGTs) are involved in phase II conjugation reactions of xenobiotics and differences in their isoform activities result in interindividual kinetic differences of UGT probe substrates. Here, extensive literature searches were performed to identify probe substrates (14) for various UGT isoforms (UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B7 and UGT2B15) and frequencies of human polymorphisms. Chemical-specific pharmacokinetic data were collected in a database to quantify interindividual differences in markers of acute (Cmax) and chronic (area under the curve, clearance) exposure. Using this database, UGT-related uncertainty factors were derived and compared to the default factor (i.e. 3.16) allowing for interindividual differences in kinetics. Overall, results show that pharmacokinetic data are predominantly available for Caucasian populations and scarce for other populations of different geographical ancestry. Furthermore, the relationships between UGT polymorphisms and pharmacokinetic parameters are rarely addressed in the included studies. The data show that UGT-related uncertainty factors were mostly below the default toxicokinetic uncertainty factor of 3.16, with the exception of five probe substrates (1-OH-midazolam, ezetimibe, raltegravir, SN38 and trifluoperazine), with three of these substrates being metabolised by the polymorphic isoform 1A1. Data gaps and future work to integrate UGT-related variability distributions with in vitro data to develop quantitative in vitro-in vivo extrapolations in chemical risk assessment are discussed.

Association Between the rs1229984 Polymorphism in the Alcohol Dehydrogenase 1B Gene and Risk for Restless Legs Syndrome

Background/Objectives

Several studies have raised the possibility of an association between alcohol consumption and the risk of developing restless legs syndrome (RLS). Moreover, an important percentage of patients under alcohol detoxification therapy develop RLS symptoms that fulfil the criteria for idiopathic RLS during alcohol withdrawal. We have aimed to establish the possible association between two common single nucleotide polymorphisms (SNPs) in the alcohol-dehydrogenase 1B (ADH1B) gene and the risk for RLS.

Methods

We studied, using specific TaqMan assays, the genotype and allelic variant frequencies of ADH1B rs1229984 and ADH1B rs6413413 SNPs in 205 RLS patients and 505 gender-matched healthy controls.

Results

The sum of the frequencies of rs1229984CT and rs1229984TT genotypes, as well as the frequency of the rs1229984T allelic variant, was significantly higher in RLS patients than in controls, both in the whole group and in females. The frequencies of genotypes and allelic variants of the rs6413413 SNP were similar between the two groups. RLS patients with the rs1229984CT genotype were younger, and those with the rs122984TT genotype older, at onset of RLS symptoms than those with the rs1229984CC genotype. None of the studied SNPs were related either with positivity of family history for RLS or with RLS severity.

Conclusions

These results suggest an association between rs1229984 SNP and the risk for RLS.

Association between the missense alcohol dehydrogenase rs1229984T variant with the risk for Parkinson's disease in women

BACKGROUND/OBJECTIVE

Several meta-analyses including retrospective case-control studies have shown that the risk of developing Parkinson's disease (PD) correlates inversely with alcohol consumption and (PD), although the results of prospective longitudinal studies are far from being conclusive. The reasons for this inverse association are not well-known. Because alcohol dehydrogenase is one of the most important alcohol-detoxification enzymes, we tried to replicate a putative association of the risk of developing PD with two missense gene variations affecting the alcohol dehydrogenase 1B (ADH1B) gene (one of them related with aversive effects to alcohol).

METHODS

In a cohort composed of 629 PD patients and 865 age- and gender-matched healthy individuals, we analyzed genotypes and allele frequencies for two common missense ADH1B single nucleotide polymorphisms (SNPs), namely rs1229984 (His48Arg) and rs6413413 (Thr60Ser) using specifically designed TaqMan assays.

RESULTS

The frequency of individuals carrying rs1229984T alleles in homozygosity or in heterozygosity was higher in PD than in controls in the whole study cohort (P < 0.001 and P = 0.005, respectively), and in women (P < 0.001 and P < 0.001, respectively). The genotypes for rs6413413 were similar in PD patients and control subjects. Age at onset of PD patients was not statistically related to rs1229984 or rs6413413 genotypes.

CONCLUSIONS

The missense variant rs1229984T is statistically associated with the risk of developing PD mainly in women, which could explain differences in alcohol consumption in this gender.

Reverse dosimetry modeling of toluene exposure concentrations based on biomonitoring levels from the Canadian health measures survey

Biomonitoring might provide useful estimates of population exposure to environmental chemicals. However, data uncertainties stemming from interindividual variability are common in large population biomonitoring surveys. Physiologically based pharmacokinetic (PBPK) models might be used to account for age- and gender-related variability in internal dose. The objective of this study was to reconstruct air concentrations consistent with blood toluene measures reported in the third Canadian Health Measures Survey using reverse dosimetry PBPK modeling techniques. Population distributions of model's physiological parameters were described based upon age, weight, and size for four subpopulations (12-19, 20-39, 40-59, and 60-79 years old). Monte Carlo simulations applied to PBPK modeling allowed converting the distributions of venous blood measures of toluene obtained from CHMS into related air levels. Based upon blood levels observed at the 50^th, 90^th and 95^th percentiles, corresponding air toluene concentrations were estimated for teenagers aged 12-19 years as being, respectively, 0.009, 0.04 and 0.06 ppm. Similarly, values were computed for adults aged 20-39 years (0.007, 0.036, and 0.06 ppm), 40-59 years (0.007, 0.036 and 0.06 ppm) and 60-79 years (0.006, 0.022 and 0.04 ppm). These estimations are well below Health Canada's maximum recommended chronic air guidelines for toluene. In conclusion, PBPK modeling and reverse dosimetry may be combined to help interpret biomonitoring data for chemical exposure in large population surveys and estimate the associated toxicological health risk.

Pharmacokinetic data reduce uncertainty in the acceptable daily intake for benzoic acid and its salts

The current acceptable daily intake (ADI) for benzoic acid and its salts as food additives is 0-5 mg/kg body weight. This accounts for a total uncertainty factor (UF) of 100, which includes a default factor of 10 for interspecies differences. Based on pharmacokinetic data in rodents and humans, we derived a chemical-specific adjustment factor (CSAF) of 2 for the pharmacokinetic component of the interspecies UF. Additional analyses indicate that this CSAF is conservative and interspecies differences between rats and humans are likely closer to unity. Human clinical studies indicate that the pharmacokinetics of benzoic acid and its salts are similar in children and adults, and that there is a lack of adverse events in humans at doses comparable to the no observed adverse effect level (NOAEL) in rodents; this suggests that the pharmacokinetic UF for intraspecies variability, as well as the pharmacodynamic components of the UFs, may also be reduced, although we did not calculate to what degree. In conclusion, the total UF can be reduced to 50 (2 for interspecies differences in pharmacokinetics, 2.5 for interspecies differences in pharmacodynamics, and 10 for intraspecies variability), which would increase the ADI to 0-10 mg/kg body weight.

Guidance on the risk assessment of substances present in food intended for infants below 16 weeks of age

Following a request from the European Commission to EFSA, the EFSA Scientific Committee (SC) prepared a guidance for the risk assessment of substances present in food intended for infants below 16 weeks of age. In its approach to develop this guidance, the EFSA SC took into account, among others, (i) an exposure assessment based on infant formula as the only source of nutrition; (ii) knowledge of organ development in human infants, including the development of the gut, metabolic and excretory capacities, the brain and brain barriers, the immune system, the endocrine and reproductive systems; (iii) the overall toxicological profile of the substance identified through the standard toxicological tests, including critical effects; (iv) the relevance for the human infant of the neonatal experimental animal models used. The EFSA SC notes that during the period from birth up to 16 weeks, infants are expected to be exclusively fed on breast milk and/or infant formula. The EFSA SC views this period as the time where health-based guidance values for the general population do not apply without further considerations. High infant formula consumption per body weight is derived from 95th percentile consumption. The first weeks of life is the time of the highest relative consumption on a body weight basis. Therefore, when performing an exposure assessment, the EFSA SC proposes to use the high consumption value of 260 mL/kg bw per day. A decision tree approach is proposed that enables a risk assessment of substances present in food intended for infants below 16 weeks of age. The additional information needed when testing substances present in food for infants below 16 weeks of age and the approach to be taken for the risk assessment are on a case-by-case basis, depending on whether the substance is added intentionally to food and is systemically available.

Derivation of an occupational exposure limit for an inhalation analgesic methoxyflurane (Penthrox(®))

Methoxyflurane (MOF) a haloether, is an inhalation analgesic agent for emergency relief of pain by self administration in conscious patients with trauma and associated pain. It is administered under supervision of personnel trained in its use. As a consequence of supervised use, intermittent occupational exposure can occur. An occupational exposure limit has not been established for methoxyflurane. Human clinical and toxicity data have been reviewed and used to derive an occupational exposure limit (referred to as a maximum exposure level, MEL) according to modern principles. The data set for methoxyflurane is complex given its historical use as anaesthetic. Distinguishing clinical investigations of adverse health effects following high and prolonged exposure during anaesthesia to assess relatively low and intermittent exposure during occupational exposure requires an evidence based approach to the toxicity assessment and determination of a critical effect and point of departure. The principal target organs are the kidney and the central nervous system and there have been rare reports of hepatotoxicity, too. Methoxyflurane is not genotoxic based on in vitro bacterial mutation and in vivo micronucleus tests and it is not classifiable (IARC) as a carcinogenic hazard to humans. The critical effect chosen for development of a MEL is kidney toxicity. The point of departure (POD) was derived from the concentration response relationship for kidney toxicity using the benchmark dose method. A MEL of 15 ppm (expressed as an 8 h time weighted average (TWA)) was derived. The derived MEL is at least 50 times higher than the mean observed TWA (0.23 ppm) for ambulance workers and medical staff involved in supervising use of Penthrox. In typical treatment environments (ambulances and treatment rooms) that meet ventilation requirements the derived MEL is at least 10 times higher than the modelled TWA (1.5 ppm or less) and the estimated short term peak concentrations are within the MEL. The odour threshold for MOF of 0.13-0.19 ppm indicates that the odour is detectable well below the MEL. Given the above considerations the proposed MEL is health protective.

Proteomic changes in Corbicula fluminea exposed to wastewater from a psychiatric hospital

The increase use of pharmaceutical compounds in veterinary practice and human population results in the ubiquitous presence of these compounds in aquatic ecosystems. Because pharmaceuticals are highly bioactive, there is concern about their toxicological effects in aquatic organisms. Therefore, the aim of this study was to assess the effects of an effluent from a psychiatric hospital (containing a complex mixture of 25 pharmaceutical compounds from eleven therapeutic classes) on the freshwater clam Corbicula fluminea using a proteomic approach. The exposure of C. fluminea to this complex effluent containing anxiolytics, analgesics, lipid regulators, beta blockers, antidepressants, antiepileptics, antihistamines, antihypertensives, antiplatelets and antiarrhythmics induced protein changes after 1 day of exposure in clam gills and digestive gland more evident in the digestive gland. These changes included increase in the abundance of proteins associated with structural (actin and tubulin), cellular functions (calreticulin, proliferating cell nuclear antigen (PCNA), T complex protein 1 (TCP1)) and metabolism (aldehyde dehydrogenase (ALDH), alcohol dehydrogenase, 6 phosphogluconate dehydrogenase). Results from this study indicate that calreticulin, PCNA, ALDH and alcohol dehydrogenase in the digestive gland and T complex protein 1 (TCP1)) and 6 phosphogluconate dehydrogenase in the gills represent useful biomarkers for the ecotoxicological characterization of psychiatric hospital effluents in this species.

Approaches to the Development of Human Health Toxicity Values for Active Pharmaceutical Ingredients in the Environment

Management of active pharmaceutical ingredients (API) in the environment is challenging because these substances represent a large and diverse group of compounds. Advanced wastewater treatment technologies that can remove API tend to be costly. Because of the potential resources required to address API in the environment, there is a need to establish environmental benchmarks that can serve as targets for treatment and release. To date, there are several different approaches that have been taken to derive human health toxicity values for API. These methods include traditional risk assessment approaches that calculate "safe" doses using experimental data and uncertainty (safety) factors; point of departure (POD), which starts from a therapeutic human dose and applies uncertainty factors; and threshold of toxicological concern (TTC), a generic approach that establishes threshold values across broad classes of chemicals based on chemical structure. To evaluate the use of these approaches, each of these methods was applied to three API commonly encountered in the environment: acetaminophen, caffeine, and chlorpromazine. The results indicate that the various methods of estimating toxicity values produce highly varying doses. Associated doses are well below typical intakes, or toxicity thresholds cannot be derived due to a lack of information. No uniform approach can be applied to establishing thresholds for multiple substances. Rather, an individualized approach will need to be applied to each target API.

Uncertainty and variability in human exposure limits - a chemical-specific approach for ciprofloxacin and methotrexate

Human exposure limits (HELs) for chemicals with a toxicological threshold are traditionally derived using default assessment factors that account for variations in exposure duration, species sensitivity and individual sensitivity. The present paper elaborates a probabilistic approach for human hazard characterization and the derivation of HELs. It extends the framework for evaluating and expressing uncertainty in hazard characterization recently proposed by WHO-IPCS, i.e. by the incorporation of chemical-specific data on human variability in toxicokinetics. The incorporation of human variability in toxicodynamics was based on the variation between adverse outcome pathways (AOPs). Furthermore, sources of interindividual variability and uncertainty are propagated separately throughout the derivation process. The outcome is a two-dimensional human dose distribution that quantifies the population fraction exceeding a pre-selected critical effect level with an estimate of the associated uncertainty. This enables policy makers to set separate standards for the fraction of the population to be protected and the confidence level of the assessment. The main sources of uncertainty in the human dose distribution can be identified in order to plan new research for reducing uncertainty. Additionally, the approach enables quantification of the relative risk for specific subpopulations. The approach is demonstrated for two pharmaceuticals, i.e. the antibiotic ciprofloxacin and the antineoplastic methotrexate. For both substances, the probabilistic HEL is mainly influenced by uncertainty originating from: (1) the point of departure (PoD), (2) extrapolation from sub-acute to chronic toxicity and (3) interspecies extrapolation. However, when assessing the tails of the two-dimensional human dose distributions, i.e. the section relevant for the derivation of human exposure limits, interindividual variability in toxicodynamics also becomes important.

Uncertainties in human health risk assessment of environmental contaminants: A review and perspective

Addressing uncertainties in human health risk assessment is a critical issue when evaluating the effects of contaminants on public health. A range of uncertainties exist through the source-to-outcome continuum, including exposure assessment, hazard and risk characterisation. While various strategies have been applied to characterising uncertainty, classical approaches largely rely on how to maximise the available resources. Expert judgement, defaults and tools for characterising quantitative uncertainty attempt to fill the gap between data and regulation requirements. The experiences of researching 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) illustrated uncertainty sources and how to maximise available information to determine uncertainties, and thereby provide an 'adequate' protection to contaminant exposure. As regulatory requirements and recurring issues increase, the assessment of complex scenarios involving a large number of chemicals requires more sophisticated tools. Recent advances in exposure and toxicology science provide a large data set for environmental contaminants and public health. In particular, biomonitoring information, in vitro data streams and computational toxicology are the crucial factors in the NexGen risk assessment, as well as uncertainties minimisation. Although in this review we cannot yet predict how the exposure science and modern toxicology will develop in the long-term, current techniques from emerging science can be integrated to improve decision-making.

Conservation of the coding regions of the glycine N-acyltransferase gene further suggests that glycine conjugation is an essential detoxification pathway

Thorough investigation of the glycine conjugation pathway has been neglected. No defect of the glycine conjugation pathway has been reported and this could reflect the essential role of glycine conjugation in hepatic metabolism. Therefore, we hypothesised that genetic variation in the open reading frame (ORF) of the GLYAT gene should be low and that deleterious alleles would be found at low frequencies. This hypothesis was investigated by analysing the genetic variation of the human GLYAT ORF using data available in public databases. We also sequenced the GLYAT ORF of a small cohort of South African Afrikaner Caucasian individuals. In total, data from 1537 individuals was analysed. The two most prominent GLYAT haplotypes in all populations analysed, were S156 (70%) and T17S156 (20%). The S156C199 and S156H131 haplotypes, which have a negative effect on the enzyme activity of a recombinant human GLYAT, were detected at very low frequencies. In the Afrikaner Caucasian cohort a novel Q61L SNP occurring at a high frequency (12%) was detected. The results of this study indicated that the GLYAT ORF is highly conserved and supported the hypothesis that the glycine conjugation pathway is an essential detoxification pathway. These findings emphasise the importance of future investigations to determine the in vivo capacity of the glycine conjugation pathway for the detoxification of benzoate and other xenobiotics.

Mechanisms of tolvaptan-induced toxicity in HepG2 cells

Tolvaptan, a vasopressin receptor 2 antagonist used to treat hyponatremia, has recently been reported to be associated with an increased risk of liver injury. In this study, we explored the underlying mechanisms of hepatotoxicity of tolvaptan using human HepG2 cells. Tolvaptan inhibited cell growth and caused cell death in a concentration- and time-dependent manner. Tolvaptan treatment led to delayed cell cycle progression, accompanied by decreased levels of several cyclins and cyclin-dependent kinases. Tolvaptan was found to cause DNA damage, as assessed by alkaline comet assays; this was confirmed by increased levels of 8-oxoguanine and phosphorylation of histone H2AX. Exposure of HepG2 cells to tolvaptan enhanced cytochrome C release and triggered apoptosis by modulating Bcl-2 family members. The activation of p38 contributed to tolvaptan-mediated apoptosis via down-regulation of Bcl-2. Proteasome inhibition altered tolvaptan-induced cell cycle deregulation and enhanced tolvaptan-induced apoptosis and cytotoxicity. Moreover, tolvaptan treatment induced autophagy. Inhibition of autophagy by knocking-down an autophagy-related gene increased tolvaptan-induced apoptosis and cytotoxicity. Taken together, our findings suggest that the cytotoxicity of tolvaptan results from delayed cell cycle progression, the induction of DNA damage, and the execution of apoptosis. In addition, a number of signaling pathways were perturbed by tolvaptan and played an important role in its cytotoxicity.

Diversity Outbred Mice Identify Population-Based Exposure Thresholds and Genetic Factors that Influence Benzene-Induced Genotoxicity

BACKGROUND

Inhalation of benzene at levels below the current exposure limit values leads to hematotoxicity in occupationally exposed workers.

OBJECTIVE

We sought to evaluate Diversity Outbred (DO) mice as a tool for exposure threshold assessment and to identify genetic factors that influence benzene-induced genotoxicity.

METHODS

We exposed male DO mice to benzene (0, 1, 10, or 100 ppm; 75 mice/exposure group) via inhalation for 28 days (6 hr/day for 5 days/week). The study was repeated using two independent cohorts of 300 animals each. We measured micronuclei frequency in reticulocytes from peripheral blood and bone marrow and applied benchmark concentration modeling to estimate exposure thresholds. We genotyped the mice and performed linkage analysis.

RESULTS

We observed a dose-dependent increase in benzene-induced chromosomal damage and estimated a benchmark concentration limit of 0.205 ppm benzene using DO mice. This estimate is an order of magnitude below the value estimated using B6C3F1 mice. We identified a locus on Chr 10 (31.87 Mb) that contained a pair of overexpressed sulfotransferases that were inversely correlated with genotoxicity.

CONCLUSIONS

The genetically diverse DO mice provided a reproducible response to benzene exposure. The DO mice display interindividual variation in toxicity response and, as such, may more accurately reflect the range of response that is observed in human populations. Studies using DO mice can localize genetic associations with high precision. The identification of sulfotransferases as candidate genes suggests that DO mice may provide additional insight into benzene-induced genotoxicity.

Assessment of health risks resulting from early-life exposures: Are current chemical toxicity testing protocols and risk assessment methods adequate?

Abstract Over the last couple of decades, the awareness of the potential health impacts associated with early-life exposures has increased. Global regulatory approaches to chemical risk assessment are intended to be protective for the diverse human population including all life stages. However, questions persist as to whether the current testing approaches and risk assessment methodologies are adequately protective for infants and children. Here, we review physiological and developmental differences that may result in differential sensitivity associated with early-life exposures. It is clear that sensitivity to chemical exposures during early-life can be similar, higher, or lower than that of adults, and can change quickly within a short developmental timeframe. Moreover, age-related exposure differences provide an important consideration for overall susceptibility. Differential sensitivity associated with a life stage can reflect the toxicokinetic handling of a xenobiotic exposure, the toxicodynamic response, or both. Each of these is illustrated with chemical-specific examples. The adequacy of current testing protocols, proposed new tools, and risk assessment methods for systemic noncancer endpoints are reviewed in light of the potential for differential risk to infants and young children.

Toward systems epidemiology of coffee and health

PURPOSE OF REVIEW

Coffee is one of the most widely consumed beverages in the world and has been associated with many health conditions. This review examines the limitations of the classic epidemiological approach to studies of coffee and health, and describes the progress in systems epidemiology of coffee and its correlated constituent, caffeine. Implications and applications of this growing body of knowledge are also discussed.

RECENT FINDINGS

Population-based metabolomic studies of coffee replicate coffee-metabolite correlations observed in clinical settings but have also identified novel metabolites of coffee response, such as specific sphingomyelin derivatives and acylcarnitines. Genome-wide analyses of self-reported coffee and caffeine intake and serum levels of caffeine support an overwhelming role for caffeine in modulating the coffee consumption behavior. Interindividual variation in the physiological exposure or response to any of the many chemicals present in coffee may alter the persistence and magnitude of their effects. It is thus imperative that future studies of coffee and health account for this variation.

SUMMARY

Systems epidemiological approaches promise to inform causality, parse the constituents of coffee responsible for health effects, and identify the subgroups most likely to benefit from increasing or decreasing coffee consumption.

Ecogenetics of mercury: from genetic polymorphisms and epigenetics to risk assessment and decision-making

The risk assessment of mercury (Hg), in both humans and wildlife, is made challenging by great variability in exposure and health effects. Although disease risk arises following complex interactions between genetic ("nature") and environmental ("nurture") factors, most Hg studies thus far have focused solely on environmental factors. In recent years, ecogenetic-based studies have emerged and have started to document genetic and epigenetic factors that may indeed influence the toxicokinetics or toxicodynamics of Hg. The present study reviews these studies and discusses their utility in terms of Hg risk assessment, management, and policy and offers perspectives on fruitful areas for future research. In brief, epidemiological studies on populations exposed to inorganic Hg (e.g., dentists and miners) or methylmercury (e.g., fish consumers) are showing that polymorphisms in a number of environmentally responsive genes can explain variations in Hg biomarker values and health outcomes. Studies on mammals (wildlife, humans, rodents) are showing Hg exposures to be related to epigenetic marks such as DNA methylation. Such findings are beginning to increase understanding of the mechanisms of action of Hg, and in doing so they may help identify candidate biomarkers and pinpoint susceptible groups or life stages. Furthermore, they may help refine uncertainty factors and thus lead to more accurate risk assessments and improved decision-making.

Characterization of the human kinetic adjustment factor for the health risk assessment of environmental contaminants

A default uncertainty factor of 3.16 (√10) is applied to account for interindividual variability in toxicokinetics when performing non-cancer risk assessments. Using relevant human data for specific chemicals, as WHO/IPCS suggests, it is possible to evaluate, and replace when appropriate, this default factor by quantifying chemical-specific adjustment factors for interindividual variability in toxicokinetics (also referred to as the human kinetic adjustment factor, HKAF). The HKAF has been determined based on the distributions of pharmacokinetic parameters (e.g., half-life, area under the curve, maximum blood concentration) in relevant populations. This article focuses on the current state of knowledge of the use of physiologically based algorithms and models in characterizing the HKAF for environmental contaminants. The recent modeling efforts on the computation of HKAF as a function of the characteristics of the population, chemical and its mode of action (dose metrics), as well as exposure scenario of relevance to the assessment are reviewed here. The results of these studies, taken together, suggest the HKAF varies as a function of the sensitive subpopulation and dose metrics of interest, exposure conditions considered (route, duration, and intensity), metabolic pathways involved and theoretical model underlying its computation. The HKAF seldom exceeded the default value of 3.16, except in very young children (i.e., <≈ 3 months) and when the parent compound is the toxic moiety. Overall, from a public health perspective, the current state of knowledge generally suggest that the default uncertainty factor is sufficient to account for human variability in non-cancer risk assessments of environmental contaminants.

Dispelling urban myths about default uncertainty factors in chemical risk assessment--sufficient protection against mixture effects?

Assessing the detrimental health effects of chemicals requires the extrapolation of experimental data in animals to human populations. This is achieved by applying a default uncertainty factor of 100 to doses not found to be associated with observable effects in laboratory animals. It is commonly assumed that the toxicokinetic and toxicodynamic sub-components of this default uncertainty factor represent worst-case scenarios and that the multiplication of those components yields conservative estimates of safe levels for humans. It is sometimes claimed that this conservatism also offers adequate protection from mixture effects. By analysing the evolution of uncertainty factors from a historical perspective, we expose that the default factor and its sub-components are intended to represent adequate rather than worst-case scenarios. The intention of using assessment factors for mixture effects was abandoned thirty years ago. It is also often ignored that the conservatism (or otherwise) of uncertainty factors can only be considered in relation to a defined level of protection. A protection equivalent to an effect magnitude of 0.001-0.0001% over background incidence is generally considered acceptable. However, it is impossible to say whether this level of protection is in fact realised with the tolerable doses that are derived by employing uncertainty factors. Accordingly, it is difficult to assess whether uncertainty factors overestimate or underestimate the sensitivity differences in human populations. It is also often not appreciated that the outcome of probabilistic approaches to the multiplication of sub-factors is dependent on the choice of probability distributions. Therefore, the idea that default uncertainty factors are overly conservative worst-case scenarios which can account both for the lack of statistical power in animal experiments and protect against potential mixture effects is ill-founded. We contend that precautionary regulation should provide an incentive to generate better data and recommend adopting a pragmatic, but scientifically better founded approach to mixture risk assessment.

Optimization of trans-Resveratrol bioavailability for human therapy

We have developed an innovative soluble galenic form to overcome the low absorption of trans-Resveratrol (t-Res) as a dry powder. We present here data on pharmacokinetics, bioavailability, and toxicity of t-Res in human volunteers treated with this soluble form, plus additional data on biological effects in rodents. Fifteen healthy volunteers of both sexes received 40 mg of t-Res in two forms, the soluble formulation (caplets) and the original powder (capsules), in a crossover design. Blood samples were collected at 15 min, 30 min, and every hour for 5 h. Plasma concentrations of t-Res and its metabolites were analyzed by liquid chromatography and mass spectrometry. The single dose (40 mg) of the soluble t-Res was well absorbed and elicited biologically efficient blood levels (0.1-6 μM) for several hours, despite metabolization into glucuronide and sulfate conjugates coupled to renal elimination. In contrast, t-Res administered as a dry powder barely elicited efficient blood levels for a short duration. The new formulation led to 8.8-fold higher t-Res levels in plasma versus the powder. t-Res metabolism was not modified and neither intolerance nor toxicity were observed during the study and the following week. The soluble formulation elicited a robust anti-inflammatory effect in various tissues of mice fed a high-fat diet, while dry powder t-Res was almost inactive. Our data suggest that significant improvements in t-Res bioavailability and efficiency can be obtained by this soluble galenic form, also allowing lower doses. The use of t-Res in human therapy is thus greatly facilitated and the toxicity risk is reduced.

Preparation and culture of precision-cut organ slices from human and animal

1.Human and animal precision-cut organ slices are being widely used to obtain drug metabolism and toxicity profiles in vitro. These data are then used to predict what might be seen in human patients. The accuracy of this prediction and extrapolation of the findings based on human or animal in vitro systems to the findings that occur in vivo is dependent on both the quality of the tissue itself and the quality of the in vitro system. 2.The quality of human organs used in research is dependent on procurement methods, warm ischaemia time, preservation solutions, cold ischaemia time, and donor-specific factors. It is important to confirm that the organs being used are highly viable and fully functional before using them in scientific studies. 3.The optimal preparation and incubation of organ slices is also essential in maintaining slice viability and function. It is important to prepare the slices in a cold preservation solution, to prepare the slices at a correct thickness, and to incubate the slices in a system where the slice rotates in out of the oxygen atmosphere and medium. 4.Meeting the criteria outlined here will lead to successful organ slice cultures for investigating drug-induced mechanisms and organ-specific toxicity.

Comparisons between in vitro whole cell imaging and in vivo zebrafish-based approaches for identifying potential human hepatotoxicants earlier in pharmaceutical development

Drug-induced liver injury (DILI) is a major cause of attrition during both the early and later stages of the drug development and marketing process. Reducing or eliminating drug-induced severe liver injury, especially those that lead to liver transplants or death, would be tremendously beneficial for patients. Therefore, developing new pharmaceuticals that have the highest margins and attributes of hepatic safety would be a great accomplishment. Given the current low productivity of pharmaceutical companies and the high costs of bringing new medicines to market, any early screening assay(s) to identify and eliminate pharmaceuticals with the potential to cause severe liver injury in humans would be of economic value as well. The present review discusses the background, proof-of-concept, and validation studies associated with high-content screening (HCS) by two major pharmaceutical companies (Pfizer Inc and Jansen Pharmaceutical Companies of Johnson & Johnson) for detecting compounds with the potential to cause human DILI. These HCS assays use fluorescent-based markers of cell injury in either human hepatocytes or HepG2 cells. In collaboration with Evotec, an independent contract lab, these two companies also independently evaluated larval zebrafish as an early-stage in vivo screen for hepatotoxicity in independently conducted, blinded assessments. Details about this model species, the need for bioanalysis, and, specifically, the outcome of the phenotypic-based zebrafish screens are presented. Comparing outcomes in zebrafish against both HCS assays suggests an enhanced detection for hepatotoxicants of most DILI concern when used in combination with each other, based on the U.S. Food and Drug Administration DILI classification list.

Modeling the Human Kinetic Adjustment Factor for Inhaled Volatile Organic Chemicals: Whole Population Approach versus Distinct Subpopulation Approach

The objective of this study was to evaluate the impact of whole- and sub-population-related variabilities on the determination of the human kinetic adjustment factor (HKAF) used in risk assessment of inhaled volatile organic chemicals (VOCs). Monte Carlo simulations were applied to a steady-state algorithm to generate population distributions for blood concentrations (CAss) and rates of metabolism (RAMs) for inhalation exposures to benzene (BZ) and 1,4-dioxane (1,4-D). The simulated population consisted of various proportions of adults, elderly, children, neonates and pregnant women as per the Canadian demography. Subgroup-specific input parameters were obtained from the literature and P3M software. Under the “whole population” approach, the HKAF was computed as the ratio of the entire population's upper percentile value (99th, 95th) of dose metrics to the median value in either the entire population or the adult population. Under the “distinct subpopulation” approach, the upper percentile values in each subpopulation were considered, and the greatest resulting HKAF was retained. CAss-based HKAFs that considered the Canadian demography varied between 1.2 (BZ) and 2.8 (1,4-D). The “distinct subpopulation” CAss-based HKAF varied between 1.6 (BZ) and 8.5 (1,4-D). RAM-based HKAFs always remained below 1.6. Overall, this study evaluated for the first time the impact of underlying assumptions with respect to the interindividual variability considered (whole population or each subpopulation taken separately) when determining the HKAF.

Population effects and variability

Chemical risk assessment for human health requires a multidisciplinary approach through four steps: hazard identification and characterization, exposure assessment, and risk characterization. Hazard identification and characterization aim to identify the metabolism and elimination of the chemical (toxicokinetics) and the toxicological dose-response (toxicodynamics) and to derive a health-based guidance value for safe levels of exposure. Exposure assessment estimates human exposure as the product of the amount of the chemical in the matrix consumed and the consumption itself. Finally, risk characterization evaluates the risk of the exposure to human health by comparing the latter to with the health-based guidance value. Recently, many research efforts in computational toxicology have been put together to characterize population variability and uncertainty in each of the steps of risk assessment to move towards more quantitative and transparent risk assessment. This chapter focuses specifically on modeling population variability and effects for each step of risk assessment in order to provide an overview of the statistical and computational tools available to toxicologists and risk assessors. Three examples are given to illustrate the applicability of those tools: derivation of pathway-related uncertainty factors based on population variability, exposure to dioxins, dose-response modeling of cadmium.

Polybrominated diphenyl ethers, hydroxylated polybrominated diphenyl ethers, and measures of thyroid function in second trimester pregnant women in California

Prenatal exposure to polybrominated diphenyl ethers (PBDEs) may disrupt thyroid function and contribute to adverse neurodevelopmental outcomes. We conducted a pilot study to explore the relationship between serum concentrations of lower-brominated PBDEs (BDE-17 to -154), higher-brominated PBDEs (BDE-183 to -209), and hydroxylated PBDE metabolites (OH-PBDEs) with measures of thyroid function in pregnant women. Concentrations of PBDEs, OH-PBDEs, thyroid-stimulating hormone (TSH), total thyroxine (T(4)), and free T(4) were measured in serum samples collected between 2008 and 2009 from 25 second trimester pregnant women in California. Median concentrations of lower-brominated PBDEs and OH-PBDEs were the highest reported to date in pregnant women. Median concentrations of BDE-47 and the sum of lower-brominated PBDEs (ΣPBDE(5)) were 43.1 ng/g lipid and 85.8 ng/g lipid, respectively, and the sum of OH-PBDEs (ΣOH-PBDE(4)) was 0.084 ng/mL. We observed a positive association between the weighted sum of chemicals known to bind to transthyretin (ΣTTR binders) and TSH levels. We also found positive associations between TSH and ΣPBDE(5), ΣOH-PBDE(4), BDE-47, BDE-85, 5-OH-BDE47, and 4'-OH-BDE49, and an inverse association with BDE-207. Relationships with free and total T(4) were weak and inconsistent. Our results indicate that PBDE exposures are elevated in pregnant women in California and suggest a relationship with thyroid function. Further investigation is warranted to characterize the risks of PBDE exposures during pregnancy.

The impact of CYP2E1 genetic variability on risk assessment of VOC mixtures

Humans are simultaneously exposed to multiple chemicals in the environment. Many of the chemicals use the same enzymes in their metabolic pathways. Competitive inhibition may occur as one of the possible interactions between the xenobiotics in human body. For example, many volatile organic compounds (VOCs) are metabolized using P450 enzymes, specifically CYP2E1. Inheritable gene alterations may result in changes of function of the enzymes in different human subpopulations. Variations in quantity and/or quality of particular isoenzymes may cause differences in the metabolism of VOCs. These variations may cause higher sensitivity in certain populations. Using examples of three different mixtures, this review paper outlines the variances in CYP2E1 isoenzymes, effect of exposure to such mixtures on sensitive populations, and approaches to mixtures risk assessment.

Evaluation of the impact of the exposure route on the human kinetic adjustment factor

The objective of this study was to assess the impact of the exposure route on the human kinetic adjustment factor (HKAF), for which a default value of 3.16 is used in non-cancer risk assessment. A multi-route PBPK model was modified from the literature and used for computing the internal dose metrics in adults, neonates, children, elderly and pregnant women following three route-specific scenarios to chloroform, bromoform, tri- or per-chloroethylene (TCE or PERC). These include 24-h inhalation exposure, body-weight adjusted oral exposure and 30 min dermal exposure to contaminated drinking water. Distributions for body weight (BW), height (BH) and hepatic cytochrome P450 2E1 (CYP2E1) content were obtained from the literature, whereas model parameters (flows, volumes) were calculated from BW and BH. Monte Carlo simulations were performed and the HKAF was calculated as the ratio of the 95th percentile value of internal dose metrics in subpopulation to the 50th percentile value in adults. On the basis of the area under the parent compound's arterial blood concentration vs time curve (AUC(pc)), highest HKAFs were obtained in neonates for every scenario considered, and were the highest for bromoform (range: 3.6-7.4). Exceedance of the default value based on AUC(PC) was also observed for an oral exposure to chloroform in neonates (4.9). In all other cases, HKAFs remained below the default value. Overall, this study has pointed out the dependency of the HKAF on the exposure route, dose metrics and subpopulation considered, as well as characteristics of the chemicals investigated.