Uncertainty factors for chemical risk assessment. human variability in the pharmacokinetics of CYP1A2 probe substrates

Current therapeutic targets and multifaceted physiological impacts of caffeine

Caffeine, which shares consubstantial structural similarity with purine adenosine, has been demonstrated as a nonselective adenosine receptor antagonist for eliciting most of the biological functions at physiologically relevant dosages. Accumulating evidence supports caffeine's beneficial effects against different disorders, such as total cardiovascular diseases and type 2 diabetes. Conversely, paradoxical effects are also linked to caffeine ingestion in humans including hypertension-hypotension and tachycardia-bradycardia. These observations suggest the association of caffeine action with its ingested concentration and/or concurrent interaction with preferential molecular targets to direct explicit events in the human body. Thus, a coherent analysis of the functional targets of caffeine, relevant to normal physiology, and disease pathophysiology, is required to understand the pharmacology of caffeine. This review provides a broad overview of the experimentally validated targets of caffeine, particularly those of therapeutic interest, and the impacts of caffeine on organ-specific physiology and pathophysiology. Overall, the available empirical and epidemiological evidence supports the dose-dependent functional activities of caffeine and advocates for further studies to get insights into the caffeine-induced changes under specific conditions, such as asthma, DNA repair, and cancer, in view of its therapeutic applications.

Variability of caffeine content in green and roasted Coffea arabica regarding the origin, post-harvest processing, and altitude, and overview of recommended daily allowance

Caffeine content is a crucial attribute of coffee. Its concentration and thus maximum cups of Coffea arabica from Africa, Asia, Central America, and South America from different altitudes of growing areas, altitude, and process using different post-harvest processing (dry, wet, and pulped natural). Our results suggest that geographical origin might affect the alkaloid concentration in C. arabica. The caffeine concentration pattern in green samples was as follows: Central America > South America > Asia > Africa. Altitude affected the concentrations, lowlands > midlands > highlands, however, not significantly. Given caffeine is thermostable, the medium roasting process did not affect the concentration of caffeine directly, but a small increase was observed. Scientific opinion on the safety of habitual caffeine consumption of up to 400 mg per day does not raise safety concerns for non-pregnant adults. A cup (7 g coffee in 120 mL of water) was used for recalculation. Results suggest that mostly highlands and midlands coffee from Africa reached levels of caffeine that might be consumed in more than 5.5 cups a day.

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.

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.

OpenCYP: An open source database exploring human variability in activities and frequencies of polymophisms for major cytochrome P-450 isoforms across world populations

The open source database "OpenCYP database" has been developed based on the results of extensive literature searches from the peer-reviewed literature. OpenCYP provides data on human variability on baseline of activities and polymophism frequencies for selected cytochrome P-450 isoforms (CYP1A2, CYP2A6, CYP2D6, CYP3A4/3A5 and CYP3A7) in healthy adult populations from world populations. CYP enzymatic activities were generally expressed as the metabolic ratio (MR) between an unchanged probe drug and its metabolite(s) in urine or plasma measured in healthy adults. Data on other age groups were very limited and fragmented, constituting an important data gap. Quantitative comparisons were often hampered by the different experimental conditions used. However, variability was quite limited for CYP1A2, using caffeine as a probe substrate, with a symmetrical distribution of metabolic activity values. For CYP3A4, human variability was dependent on the probe substrate itself with low variability when data considering the dextromethorphan/demethilathed metabolite MR were used and large variability when the urinary 6β-hydroxycortisol/cortisol ratio was used. The largest variability in CYP activity was shown for CYP2D6 activity, after oral dosing of dextromethorphan, for which genetic polymorphisms are well characterised and constitute a significant source of variability. It is foreseen that the OpenCYP database can contribute to promising tools to support the further development of QIVIVE and PBK models for human risk assessment of chemicals particularly when combined with information on isoform-specific content in cells using proteomic approaches.

An approach to evaluate metabolite-related phototoxicity with combined use of photochemical properties and skin deposition

Some chemicals have been reported to cause metabolite-related phototoxicity, and this study aimed to verify the applicability of photosafety assessment based on photochemical and pharmacokinetic properties to evaluate the metabolite-related phototoxicity risk. The phototoxic risk of imipramine (IMI) and its metabolite, desipramine (DMI), was evaluated by photochemical and pharmacokinetic analyses. IMI and DMI were found to have similar photoreactivities based on the generation of reactive oxygen species. The skin concentrations of IMI and DMI reached maximal levels at approximately 1 and 4 h, respectively, after oral administration of IMI (10 mg/kg), and DMI showed high skin deposition compared with IMI. According to the results, DMI was identified as a contributor to phototoxicity induced by orally-taken IMI. In in vivo phototoxicity testing, ultraviolet A irradiation from 3 to 6 h after oral administration of IMI (100 mg/kg) caused more potent phototoxic reactions compared with that from 0 to 3 h, and DMI yielded by metabolism of IMI would be associated with phototoxic reactions caused by orally-administered IMI. In addition to the data on IMI, a parent chemical, photochemical and pharmacokinetic profiling of its metabolite, DMI, led to reliable phototoxicity prediction of orally-administered IMI. Thus, characterization of the photosafety of metabolites would generate reliable information on the phototoxicity risk of parent chemicals, and the proposed strategy may facilitate comprehensive photosafety assessment of drug candidates in pharmaceutical development.

Human Variability in Carboxylesterases and carboxylesterase-related Uncertainty Factors for Chemical Risk Assessment

Carboxylesterases (CES) are an important class of enzymes involved in the hydrolysis of a range of chemicals and show large inter-individual variability in vitro. An extensive literature search was performed to identify in vivo probe substrates for CES1 and CES2 together with their protein content and enzymatic activity. Human pharmacokinetic (PK) data on Cmax, clearance, and AUC were extracted from 89 publications and Bayesian meta-analysis was performed using a hierarchical model to derive CES-related variability distributions and related uncertainty factors (UF). The CES-related variability indicated that 97.5% of healthy adults are covered by the kinetic default UF (3.16), except for clopidogrel and dabigatran etexilate. Clopidogrel is metabolised for a small amount by the polymorphic CYP2C19, which can have an impact on the overall pharmacokinetics, while the variability seen for dabigatran etexilate might be due to differences in the absorption, since this can be influenced by food intake. The overall CES-related variability was moderate to high in vivo (<CV 50%), which might be due to possible polymorphism in the enzyme but also to the small sample size available per chemical. The presented CES-related variability can be used in combination with in vitro data to derive pathway-specific distributions.

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.

Precision Dosing Management with Intelligent Computing in Digital Health

Pediatric dosing is not only critical for successful pediatric trials in drug development but also paramount to safety and effective treatment at bedside. Due to the complex pharmacokinetic of children compared to adults, several challenges are posed in managing dosing precisely during drug development and after drug approval to clinicians. In particular, given the real-world practice, understanding the impact of development on the dose-exposure-response relationship is essential in optimizing the dosing to children of different ages. In this paper we propose a novel intelligent computing framework to examine how the growth and maturation create size- and age-dependent variability in pharmacokinetics and pharmacodynamics, and summarize the use of modeling-based approaches for dose finding in pediatric drug development, allowing clinicians to anticipate probable treatment effects and to have a higher likelihood of achieving optimal dose regimens early, as well as reducing the drug development cycling time and cost.

Two Approaches for Evaluating the Effects of Galangin on the Activities and mRNA Expression of Seven CYP450

Galangin is a marker compound of honey and Alpinia officinarum Hance that exhibits great potential for anti-microbial, anti-diabetic, anti-obesity, anti-tumour and anti-inflammatory applications. Galangin is frequently consumed in combination with common clinical drugs. Here, we evaluated the effects of galangin on cytochrome P450 (CYP)-mediated metabolism, using two different approaches, to predict drug⁻drug interactions. Male Sprague Dawley rats were administered galangin daily for 8 weeks. A "cocktail-probes" approach was employed to evaluate the activities of different CYP450 enzymes. Blood samples of seven probe drugs were analysed using liquid chromatography-tandem mass spectrometry in positive and negative electrospray-ionisation modes. Pharmacokinetic parameters were calculated to identify statistical differences. CYP mRNA-expression levels were investigated in real-time quantitative polymerase chain reaction experiments. The galangin-treated group showed significantly decreased AUC_0⁻∞ and C_max values for CYP1A2, and CYP2B3. The galangin-treated group showed significantly increased AUC_0⁻∞ and C_max values for CYP2C13 and CYP3A1. No significant influences were observed in the pharmacokinetic profiles of CYP2C11, CYP2D4 and CYP2E1. The mRNA-expression results were consistent with the pharmacokinetic results. Thus, CYP450 enzyme activities may be altered by long-term galangin administration, suggesting galangin to be a promising candidate molecule for enhancing oral drug bioavailability and chemoprevention and reversing multidrug resistance.

The clinical toxicology of caffeine: A review and case study

Caffeine is a widely recognized psychostimulant compound with a long history of consumption by humans. While it has received a significant amount of attention there is still much to be learned with respect to its toxicology in humans, especially in cases of overdose. A review of the history of consumption and the clinical toxicology of caffeine including clinical features, pharmacokinetics, toxicokinetics, a thorough examination of mechanism of action and management/treatment strategies are undertaken. While higher (i.e., several grams) quantities of caffeine are known to cause toxicity and potentially lethality, cases of mainly younger individuals who have experienced severe side effects and death despite consuming doses not otherwise known to cause such harm is troubling and deserves further study. An attempted case reconstruction is performed in an effort to shed light on this issue with a focus on the pharmacokinetics and pharmacodynamics of caffeine.

Caffeine toxicity in forensic practice: possible effects and under-appreciated sources

Caffeine is considered a very safe stimulant and is widely consumed in a variety of forms, from pure caffeine to beverages and foods. Typically, death is only seen when gram quantities of caffeine are consumed, usually in suicide attempts. Even in this scenario, death is rare. However, there are special populations that need to be considered in forensic presentations, who may be at greater risk. These include poor metabolizers, people with liver disease, and people with cardiac conditions, who can die as a result of caffeine intake at levels well below what is ordinarily considered toxic. Also, caffeine intake may be hidden. For example, herbal medicines with substantial caffeine content may not disclose these concentrations on their product label. The role of caffeine in medicolegal deaths is yet to be defined, however, herbal medicines and herbal weight loss supplements may represent an underappreciated source of caffeine in this context.

Bayesian Population Physiologically-Based Pharmacokinetic (PBPK) Approach for a Physiologically Realistic Characterization of Interindividual Variability in Clinically Relevant Populations

Interindividual variability in anatomical and physiological properties results in significant differences in drug pharmacokinetics. The consideration of such pharmacokinetic variability supports optimal drug efficacy and safety for each single individual, e.g. by identification of individual-specific dosings. One clear objective in clinical drug development is therefore a thorough characterization of the physiological sources of interindividual variability. In this work, we present a Bayesian population physiologically-based pharmacokinetic (PBPK) approach for the mechanistically and physiologically realistic identification of interindividual variability. The consideration of a generic and highly detailed mechanistic PBPK model structure enables the integration of large amounts of prior physiological knowledge, which is then updated with new experimental data in a Bayesian framework. A covariate model integrates known relationships of physiological parameters to age, gender and body height. We further provide a framework for estimation of the a posteriori parameter dependency structure at the population level. The approach is demonstrated considering a cohort of healthy individuals and theophylline as an application example. The variability and co-variability of physiological parameters are specified within the population; respectively. Significant correlations are identified between population parameters and are applied for individual- and population-specific visual predictive checks of the pharmacokinetic behavior, which leads to improved results compared to present population approaches. In the future, the integration of a generic PBPK model into an hierarchical approach allows for extrapolations to other populations or drugs, while the Bayesian paradigm allows for an iterative application of the approach and thereby a continuous updating of physiological knowledge with new data. This will facilitate decision making e.g. from preclinical to clinical development or extrapolation of PK behavior from healthy to clinically significant populations.

Theoretical study on the N-demethylation mechanism of theobromine catalyzed by P450 isoenzyme 1A2

Theobromine, a widely consumed pharmacological active substance, can cause undesirable muscle stiffness, nausea and anorexia in high doses ingestion. The main N-demethylation metabolic mechanism of theobromine catalyzed by P450 isoenzyme 1A2 (CYP1A2) has been explored in this work using the unrestricted hybrid density functional method UB3LYP in conjunction with the LACVP(Fe)/6-31G (H, C, N, O, S, Cl) basis set. Single-point calculations including empirical dispersion corrections were carried out at the higher 6-311++G** basis set. Two N-demethylation pathways were characterized, i.e., 3-N and 7-N demethylations, which involve the initial N-methyl hydroxylation to form carbinolamines and the subsequent carbinolamines decomposition to yield monomethylxanthines and formaldehydes. Our results have shown that the rate-limiting N-methyl hydroxylation occurs via a hydrogen atom transfer (HAT) mechanism, which proceeds in a spin-selective mechanism (SSM) in the gas phase. The carbinolamines generated are prone to decomposition via the contiguous heteroatom-assisted proton-transfer. Strikingly, 3-N demethylation is more favorable than 7-N demethylation due to its lower free energy barrier and 7-methylxanthine therefore is the optimum product reported for the demethylation of theobromine catalyzed by CYP1A2, which are in good agreement with the experimental observation. This work has first revealed the detail N-demethylation mechanisms of theobromine at the theoretical level. It can offer more significant information for the metabolism of purine alkaloid.

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.

Developmental pharmacokinetics in pediatric populations

Information on drug absorption and disposition in infants and children has increased considerably over the past 2 decades. However, the impact of specific age-related effects on pharmacokinetics, pharmacodynamics, and dose requirements remains poorly understood. Absorption can be affected by the differences in gastric pH and stomach emptying time that have been observed in the pediatric population. Low plasma protein concentrations and a higher body water composition can change drug distribution. Metabolic processes are often immature at birth, which can lead to a reduced clearance and a prolonged half-life for those drugs for which metabolism is a significant mechanism for elimination. Renal excretion is also reduced in neonates due to immature glomerular filtration, tubular secretion, and reabsorption. Limited data are available on the pharmacodynamic behavior of drugs in the pediatric population. Understanding these age effects provide a mechanistic way to identify initial doses for the pediatric population. The various factors that impact pharmacokinetics and pharmacodynamics mature towards adult values at different rates, thus requiring continual modification of drug dose regimens in neonates, infants, and children. In this paper, the age-related changes in drug absorption, distribution, metabolism, and elimination in infants and children are reviewed, and the age-related dosing regimens for this population are discussed.

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.

Assessing the impact of child/adult differences in hepatic first-pass effect on the human kinetic adjustment factor for ingested toxicants

The objective of this study was to evaluate the impact of interindividual differences in hepatic first-pass effect (FPE) on the magnitude of the human kinetic adjustment factor (HKAF) for ingested toxicants. This factor aims at replacing a default value of 3.2 used in non-cancer risk assessment. Coupled with Monte Carlo simulations, steady-state equations that account for FPE were used to obtain distributions of arterial blood concentrations (CAss) and rates of metabolism in adults, neonates, infants and toddlers continuously exposed to an oral dose of 1 μg/kg/d of theoretical CYP2E1 and CYP1A2 substrates. For such substrates exhibiting a range of blood:air partition coefficients (Pb: 1-10,000) and hepatic extraction ratios in an average adult (E(ad): 0.01-0.99), HKAFs were computed as the ratio of the 95th percentile of dose metrics for each subpopulation over the 50th percentile value in adults. The reduced hepatic enzyme content in neonates as compared to adults resulted in correspondingly diminished FPE. Consequently, HKAFs greater than 3.2 could be observed, based on CAss only, in the following cases: for some CYP2E1 substrates with E(ad) ≤ 0.3, in neonates (max.: 6.3); and for some CYP1A2 substrates with E(ad) ≤ 0.1 and 0.7, in, respectively, neonates and infants (max.: 28.3). Overall, this study pointed out the importance of accounting for child/adult differences in FPE when determining the HKAF for oral exposure.

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.

Extrapolation, uncertainty factors, and the precautionary principle

This essay examines the relationship between the precautionary principle and uncertainty factors used by toxicologists to estimate acceptable exposure levels for toxic chemicals from animal experiments. It shows that the adoption of uncertainty factors in the United States in the 1950s can be understood by reference to the precautionary principle, but not by cost-benefit analysis because of a lack of relevant quantitative data at that time. In addition, it argues that uncertainty factors continue to be relevant to efforts to implement the precautionary principle and that the precautionary principle should not be restricted to cases involving unquantifiable hazards.

Options for incorporating children's inhaled dose into human health risk assessment

Increasing attention has been placed on inhalation dosimetry in children because of children's greater air intake rate and unique windows of vulnerability for various toxicants and health outcomes. However, risk assessments have not incorporated this information because dosimetric adjustments have focused upon extrapolation across species rather than across age groups within the human population. The objectives of this study were to synthesize information regarding child/adult intake and dosimetry differences for particles and gases for potential application to risk assessment. Data and models gathered at a 2006 workshop and more recent studies were reviewed to better understand lung development and inhaled dose in children. The results show that child/adult differences exist both on a chemical intake basis and on a deposited or systemic dose basis. These differences can persist for several years and are not captured by standard intraspecies uncertainty factors or by USEPA's reference concentration (RfC) methodology. Options for incorporating children's inhalation exposures into human risk assessments include (1) 3-fold default air intake adjustment for the first 3 years of life with a reduced factor for older children; (2) superseding this default via simplified dosimetry models akin to USEPA's RfC methodology modified for children; (3) utilizing more sophisticated models with better anatomical and air flow descriptions; (4) running these models with input distributions to reflect interchild variability; (5) developing more advanced approaches involving imaging techniques and computational fluid dynamic (CFD) models. These options will enable children's inhaled dose to have a quantitative role in risk assessment that has been lacking and will establish a basis for ongoing research.

Systems toxicology approaches for understanding the joint effects of environmental chemical mixtures

Environmental mixtures of chemicals constitute a prevalent issue in ecotoxicology and the development of new methods to reduce the uncertainties associated with their ecological risk assessment is a critical research need. Historically, a number of models have been explored to predict the potential combined effects of chemicals on species. These models, especially concentration addition and the independent action, have been applied to a number of mixtures. While often providing a good prediction of joint effect, there are cases where these models can have limitations: notably in cases where there are interactions for which they fail to adequately predict joint effects. To support the better mechanistic understanding of interactions in mixture toxicology a framework to support experimental studies to investigate the basis of observed interactions is proposed. The conceptual framework is derived from the extension of a three stage scheme which has previously been applied to understand chemical bioavailability. The framework considers that interactions in mixtures result from processes related to 1) the speciation, binding and transport of chemicals in the exposure medium (external exposure); 2) the adsorption, distribution, metabolism and excretion of chemicals within the organisms (toxicokinetics); 3) associations governing the binding and toxicity of the chemical(s) at the target site (toxicodynamics). The current state of the art in (eco)toxicology in relation to investigation of the mechanisms of interactions between chemicals is discussed with particular emphasis towards the multi-disciplinary tools and techniques within environmental chemistry; toxicology; biochemistry and systems biology that can be used to address such effects.

Metabolism, variability and risk assessment

For non-genotoxic carcinogens, "thresholded toxicants", Acceptable/Tolerable Daily Intakes (ADI/TDI) represent a level of exposure "without appreciable health risk" when consumed everyday or weekly for a lifetime and are derived by applying an uncertainty factor of a 100-fold to a no-observed-adverse-effect-levels (NOAEL) or to a benchmark dose. This UF allows for interspecies differences and human variability and has been subdivided to take into account toxicokinetics and toxicodynamics with even values of 10(0.5) (3.16) for the human aspect. Ultimately, such refinements allow for chemical-specific adjustment factors and physiologically based models to replace such uncertainty factors. Intermediate to chemical-specific adjustment factors are pathway-related uncertainty factors which have been derived for phase I, phase II metabolism and renal excretion. Pathway-related uncertainty factors are presented here as derived from the result of meta-analyses of toxicokinetic variability data in humans using therapeutic drugs metabolised by a single pathway in subgroups of the population. Pathway-related lognormal variability was derived for each metabolic route. The resulting pathway-related uncertainty factors showed that the current uncertainty factor for toxicokinetics (3.16) would not cover human variability for genetic polymorphism and age differences (neonates, children, the elderly). Latin hypercube (Monte Carlo) models have also been developed using quantitative metabolism data and pathway-related lognormal variability to predict toxicokinetics variability and uncertainty factors for compounds handled by several metabolic routes. For each compound, model results gave accurate predictions compared to published data and observed differences arose from data limitations, inconsistencies between published studies and assumptions during model design and sampling. Finally, under the 6(th) framework EU project NOMIRACLE (http://viso.jrc.it/nomiracle/), novel methods to improve the risk assessment of chemical mixtures were explored (1) harmonization of the use of uncertainty factors for human and ecological risk assessment using mechanistic descriptors (2) use of toxicokinetics interaction data to derive UFs for chemical mixtures. The use of toxicokinetics data in risk assessment are discussed together with future approaches including sound statistical approaches to optimise predictability of models and recombinant technology/toxicokinetics assays to identify metabolic routes for chemicals and screen mixtures of environmental health importance.

The French approach to deriving toxicity reference values: an example using reprotoxic effects

Following the French health and environment action plan, the French Agency for Environmental and Occupational Health and Safety set up a workgroup to standardise a method of deriving toxicity reference values (TRVs). Over the last few decades, there has been increasing concern about the effect of exposure to chemicals on reproductive function, leading the group to take an interest in reprotoxic effects. This article presents the recommendations of the workgroup regarding specific reprotoxic effects. Abnormal development of foetuses and infants, together with impairment of reproduction were considered to be critical effects. Where critical windows of exposure were concerned, quantitative analysis suggested the need for several types of toxicity reference value, as a function of exposure duration: reprotoxic effects may result from acute or chronic exposure at any time of life, whilst developmental effects may occur after exposure during the pregnancy or during the lactation period. The choice of a critical study is based on epidemiological or toxicological quality criteria. The working group recommends the use of the benchmark dose approach in estimating the critical dose. Finally, the working group considered the application of uncertainty factors typically used to take into account the variability between animal and human, between different individuals, and the availability of the data.

Physiology-based simulations of a pathological condition: prediction of pharmacokinetics in patients with liver cirrhosis

BACKGROUND

Liver cirrhosis is a progressive disease characterized by loss of functional hepatocytes with concomitant connective tissue and nodule formation in the liver. The morphological and physiological changes associated with the disease substantially affect drug pharmacokinetics. Whole-body physiologically based pharmacokinetic (WB-PBPK) modelling is a predictive technique that quantitatively relates the pharmacokinetic parameters of a drug to such (patho-)physiological conditions.

OBJECTIVE

To extend an existing WB-PBPK model, based on the physiological changes associated with liver cirrhosis, which allows for prediction of drug pharmacokinetics in patients with liver cirrhosis.

METHODS

The literature was searched for quantitative measures of the physiological changes associated with the presence of Child-Pugh class A through C liver cirrhosis. The parameters that were included were the organ blood flows, cardiac index, plasma binding protein concentrations, haematocrit, functional liver volume, hepatic enzymatic activity and glomerular filtration rate. Predictions of pharmacokinetic profiles and parameters were compared with literature data for the model compounds alfentanil, lidocaine (lignocaine), theophylline and levetiracetam.

RESULTS

The predicted versus observed plasma concentration-time profiles for alfentanil and lidocaine were similar, such that the pharmacokinetic changes associated with Child-Pugh class A, B and C liver cirrhosis were adequately described. The theophylline elimination half-life was greatly increased in Child-Pugh class B and C patients compared with controls, as predicted by the model. Levetiracetam urinary excretion was consistently reduced with disease progression and very closely resembled observed values.

CONCLUSION

Consideration of the physiological differences between healthy individuals and patients with liver cirrhosis was important for the simulation of drug pharmacokinetics in this compromised group. The WB-PBPK model was altered to incorporate these physiological differences with the result of adequate simulation of drug pharmacokinetics. The information provided in this study will allow other researchers to further validate this liver cirrhosis model within a WB-PBPK model.

Approaches for applications of physiologically based pharmacokinetic models in risk assessment

Physiologically based pharmacokinetic (PBPK) models are particularly useful for simulating exposures to environmental toxicants for which, unlike pharmaceuticals, there is often little or no human data available to estimate the internal dose of a putative toxic moiety in a target tissue or an appropriate surrogate. This article reviews the current state of knowledge and approaches for application of PBPK models in the process of deriving reference dose, reference concentration, and cancer risk estimates. Examples drawn from previous U.S. Environmental Protection Agency (EPA) risk assessments and human health risk assessments in peer-reviewed literature illustrate the ways and means of using PBPK models to quantify the pharmacokinetic component of the interspecies and intraspecies uncertainty factors as well as to conduct route to route, high dose to low dose and duration extrapolations. The choice of the appropriate dose metric is key to the use of the PBPK models for the various applications in risk assessment. Issues related to whether uncertainty factors are most appropriately applied before or after derivation of human equivalent dose (or concentration) continue to be explored. Scientific progress in the understanding of life stage and genetic differences in dosimetry and their impacts on variability in susceptibility, as well as ongoing development of analytical methods to characterize uncertainty in PBPK models, will make their use in risk assessment increasingly likely. As such, it is anticipated that when PBPK models are used to express adverse tissue responses in terms of the internal target tissue dose of the toxic moiety rather than the external concentration, the scientific basis of, and confidence in, risk assessments will be enhanced.

Potential new approaches for children's inhalation risk assessment

The U.S. Environmental Protection Agency (EPA) practice of risk assessment is moving toward more thoroughly considering children's unique susceptibilities and exposure potential. Childhood is assessed as a sequence of life stages that reflects the fact that as humans develop, windows of susceptibility may appear that lead to enhanced sensitivity to exposure of environmental agents, while changes in behavior and physiology may increase exposure and dose. The U.S. EPA developed guidance in the past few years that addresses some aspects of increased susceptibility and exposure and dose. However, when it comes to considering inhalation exposure, dose, and risk, current U.S. EPA practice does not explicitly address children. The purpose here is to begin studying the adequacy of practice for children's health and to explore possible next steps in developing new methods to more accurately assess life-stage-specific differences. The existing guidelines and policies used to address potentially unique susceptibilities of children for inhaled environmental chemicals were considered, as well as what may be learned from examples of approaches that have been applied by state agencies (such as the California Environmental Protection Agency) or in the literature, to incorporate potentially unique susceptibilities and exposures to children. Finally, there is a discussion of possible approaches for considering inhalation exposure and susceptibility in U.S. EPA risk assessments.

Approaches to acrylamide physiologically based toxicokinetic modeling for exploring child-adult dosimetry differences

Dietary exposure to acrylamide is common as a result of its formation during the cooking of carbohydrate foods. This leads to widespread human exposure in adults and children alike. Acrylamide is neurotoxic and is metabolized by cytochrome P-450 (CYP) 2E1 to a mutagenic epoxide, glycidamide. This article describes a modeling framework for assessing acrylamide and glycidamide dosimetry in rats and human adults and children. The challenges in building a physiologically based toxicokinetic (PBTK) model that is compatible with existing rat and human data are described, with an emphasis on calibration against the hemoglobin adduct database. This exploratory PBTK model was adapted to children by incorporating life-stage-specific parameters consistent with children's changing physiology and metabolic capacity for processes involved in acrylamide disposition in terms of CYP2E1, glutathione conjugation, and epoxide hydrolase. Monte Carlo analysis was used to simulate the distribution of internal doses to gain an initial understanding of the range of child/adult differences possible. This analysis suggests modest dosimetry differences between children and adults, with area-under-the-curve (AUC) doses for the 99th percentile child up to fivefold greater than the median adult for both acrylamide and glycidamide. Early life immaturities tended to exert a greater effect on acrylamide than glycidamide dosimetry because immaturities in CYP2E1 and glutathione counteract one another for glycidamide AUC, but both lead to greater acrylamide dose. The analysis points toward glutathione conjugation parameters as being particularly influential and uncertain in early life, making this a key area for future research.

Assessment factors--applications in health risk assessment of chemicals

We review the scientific basis for default assessment factors used in risk assessment of nongenotoxic chemicals including the use of chemical- and pathways specific assessment factors, and extrapolation approaches relevant to species differences, age and gender. One main conclusion is that the conventionally used default factor of 100 does not cover all inter-species and inter-individual differences. We suggest that a species-specific default factor based on allometric scaling should be used for inter-species extrapolation (basal metabolic rate). Regarding toxicodynamic and remaining toxicokinetic differences we suggest that a percentile from a probabilistic distribution is chosen to derive the assessment factor. Based on the scarce information concerning the human-to-human variability it is more difficult to suggest a specific assessment factor. However, extra emphasis should be put on sensitive populations such as neonates and genetically sensitive subgroups, and also fetuses and children which may be particularly vulnerable during development and maturation. Factors that also need to be allowed for are possible gender differences in sensitivity, deficiencies in the databases, nature of the effect, duration of exposure, and route-to-route extrapolation. Since assessment factors are used to compensate for lack of knowledge we feel that it is prudent to adopt a "conservative" approach, erring on the side of protectiveness.

Development of a Physiology-Based Whole-Body Population Model for Assessing the Influence of Individual Variability on the Pharmacokinetics of Drugs

In clinical development stages, an a priori assessment of the sensitivity of the pharmacokinetic behavior with respect to physiological and anthropometric properties of human (sub-) populations is desirable. A physiology-based pharmacokinetic (PBPK) population model was developed that makes use of known distributions of physiological and anthropometric properties obtained from the literature for realistic populations. As input parameters, the simulation model requires race, gender, age, and two parameters out of body weight, height and body mass index. From this data, the parameters relevant for PBPK modeling such as organ volumes and blood flows are determined for each virtual individual. The resulting parameters were compared to those derived using a previously published model (P(3)M). Mean organ weights and blood flows were highly correlated between the two models, despite the different methods used to generate these parameters. The inter-individual variability differed greatly especially for organs with a log-normal weight distribution (such as fat and spleen). Two exemplary population pharmacokinetic simulations using ciprofloxacin and paclitaxel as model drugs showed good correlation to observed variability. A sensitivity analysis demonstrated that the physiological differences in the virtual individuals and intrinsic clearance variability were equally influential to the pharmacokinetic variability but were not additive. In conclusion, the new population model is well suited to assess the influence of individual physiological variability on the pharmacokinetics of drugs. It is expected that this new tool can be beneficially applied in the planning of clinical studies.

Human variability in hepatic and renal elimination: implications for risk assessment

Hepatic metabolism and renal excretion constitute the main routes of xenobiotic elimination in humans. Improving human risk assessment for threshold contaminants requires the incorporation of quantitative data related to their elimination (toxicokinetics) and potential toxic effects (toxicodynamics). This type of data provides a scientific basis to replace the standard uncertainty factor (UF = 10) allowing for the consideration of human variability in toxicokinetics and toxicodynamics. This review focuses on recent research efforts aiming to incorporate human variability in hepatic and renal elimination (toxicokinetics) into the risk assessment process. A therapeutic drug database was developed to quantify pathway-related variability in human phase I and phase II hepatic metabolism as well as renal excretion in subgroups of the population (healthy adults, neonates and the elderly), using data on compounds cleared primarily through each route (> 60% dose). For each subgroup of the population and elimination route, pathway-related UFs were then derived to cover 95-99% of each subgroup. Overall, the default toxicokinetic UFs would not cover neonates, the elderly for most elimination routes and any subgroup of the population for compounds metabolized via polymorphic isozymes (such as CYP2C19 and CYP2D6). These pathway-related UFs allow the incorporation of in vivo metabolism and toxicokinetic data in the risk assessment process and provide a flexible intermediate option between the default UF and chemical-specific adjustment factors (CSAFs) derived from physiologically based pharmacokinetic models. Implications of human variability in hepatic metabolism and renal excretion for chemical risk assessment are discussed.

Human and environmental risk assessment of pharmaceuticals: differences, similarities, lessons from toxicology

The presence of human and veterinary pharmaceuticals in the environment has caused increasing concern due their effects on ecological receptors. Improving the risk assessment of these compounds necessitates a quantitative understanding of their metabolism and elimination in the target organism (toxicokinetics), particularly via the ubiquitous cytochrome P-450 (CYP) system and their mechanisms of toxicity (toxicodynamics). This review focuses on a number of pharmaceuticals and veterinary medicines of environmental concern, and the differences and similarities between ecological and human risk assessment. CYP metabolism is discussed with particular reference to its ubiquity in species of ecological relevance. The important issue of pharmaceutical mixtures is discussed to assess how emerging technologies such as ecotoxicogenomics may assist in moving towards a more mechanism-based environmental risk assessment of pharmaceuticals.

Effect of Caffeine Intake 12 or 24 Hours Prior to Melatonin Intake and CYP1A2*1F Polymorphism on CYP1A2 Phenotyping by Melatonin

Earlier evidence suggests that melatonin is almost exclusively metabolised by CYP1A2 and could serve as a probe drug for CYP1A2 phenotyping. However, caffeine inhibits the metabolism of melatonin by CYP1A2 and dietary caffeine could be a potential confounder for the measurement of CYP1A2 activity with melatonin. We undertook a 3-phase cross-over study in 12 healthy volunteers to examine whether caffeine (200 mg single dose), taken 12 hr or 24 hr prior to melatonin intake, would affect the results of CYP1A2 phenotyping results as assessed by a spot sample melatonin concentration 1.5 hr after intake of 6 mg of melatonin orally. In addition we examined the influence of the CYP1A2*1F polymorphism on the phenotyping results by combining the present material with another 12 persons from a previous study. Caffeine, co-administered 12 or 24 hr prior to melatonin intake, did not have any significant effect on the 1.5 hr melatonin concentration (P=0.086 for ANOVA), but in two volunteers about 4 times increase in melatonin concentration was observed after caffeine intake 12 hr (but not 24 hr) before phenotyping with melatonin. Also, individuals homozygous for the CYP1A2*1A allele had clearly higher 1.5 hr melatonin concentration compared with the *1F/*1F or the *1F/*1A genotypes. Abstinence from caffeine for 24 hr prior to melatonin intake should be enough to overcome the possible confounding effect of caffeine on the CYP1A2 phenotyping with melatonin. Also, melatonin may be a sensitive probe to detect phenotypic differences with regard to CYP1A2*1F polymorphism. Melatonin might be, thus, advantageous for CYP1A2 phenotyping compared to the standard probe caffeine.

Human response to dioxin: aryl hydrocarbon receptor (AhR) molecular structure, function, and dose-response data for enzyme induction indicate an impaired human AhR

The aryl hydrocarbon receptor (AhR) mediates nearly all studied adverse effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and many related compounds. Binding of TCDD or related ligands to AhR is the key initiating event in downstream biochemical responses. The binding affinity of AhR for TCDD is specific to species and strain, and studies of human AhR demonstrate binding affinities approximately an order of magnitude or more lower than those observed in the most sensitive laboratory strains and species. Molecular genetic studies confirmed that human AhR shares key mutations with the DBA mouse strain that result in an "impaired" AhR (with respect to TCDD binding and responsiveness). Despite a number of polymorphisms in human AhR, the key "DBA-type" mutations appear to be a constant feature of the human AhR, and no polymorphisms have been identified that compensate for the impaired binding function conferred by these mutations. Consistent with the impaired binding status of the human AhR, human cells have consistently required approximately 10-fold higher concentrations of TCDD in vitro than rodent cells to respond with enzyme induction. Recent studies of in vivo enzyme induction-related endpoints in human populations with moderately and highly increased TCDD body burdens detected no relationship between these endpoints and TCDD body burdens at body-burden levels up to 250 ng TEQ/kg body weight, or approximately 25 times above the upper range of current general population background body burdens, while marked elevations in enzyme activity were observed in persons with body burdens above 750 ng TEQ/kg. In contrast, the more sensitive laboratory rodent strains and species exposed to TCDD exhibit significant enzyme induction at body burdens below 50 ng/kg. These interspecies data on the most sensitive and best understood response to binding of TCDD and related compounds to the AhR are consistent with the binding affinity and molecular structure data and support the hypothesis that the human AhR is less functional than the AhR of the more sensitive laboratory animals at a molecular level. Quantitative risk assessments involving interspecies extrapolation from sensitive laboratory species and strains should take these fundamental differences into account when margins of exposure and safety factors are considered.

A Mechanistic Approach for the Scaling of Clearance in Children

BACKGROUND AND OBJECTIVE

Clearance is an important pharmacokinetic concept for scaling dosage, understanding the risks of drug-drug interactions and environmental risk assessment in children. Accurate clearance scaling to children requires prior knowledge of adult clearance mechanisms and the age-dependence of physiological and enzymatic development. The objective of this research was to develop and evaluate ontogeny models that would provide an assessment of the age-dependence of clearance.

METHODS

Using in vitro data and/or in vivo clearance values for children for eight compounds that are eliminated primarily by one process, models for the ontogeny of renal clearance, cytochrome P450 (CYP) 3A4, CYP2E1, CYP1A2, uridine diphosphate glucuronosyltransferase (UGT) 2B7, UGT1A6, sulfonation and biliary clearance were developed. Resulting ontogeny models were evaluated using six compounds that demonstrated elimination via multiple pathways. The proportion of total clearance attributed to each clearance pathway in adults was delineated. Each pathway was individually scaled to the desired age, inclusive of protein-binding prediction, and summed to generate a total plasma clearance for the child under investigation. The paediatric age range included in the study was premature neonates to sub-adults.

RESULTS

There was excellent correlation between observed and predicted clearances for the model development (R2 = 0.979) and test sets (Q2 = 0.927). Clearance in premature neonates could also be well predicted (development R2 = 0.951; test Q2 = 0.899).

CONCLUSION

Paediatric clinical trial development could greatly benefit from clearance scaling, particularly in guiding dosing regimens. Furthermore, since the proportion of clearance via different elimination pathways is age-dependent, information could be gained on the developmental extent of drug-drug interactions.

Review and analysis of inhalation dosimetry methods for application to children's risk assessment

Young children have a greater ventilation rate per body weight or pulmonary surface area as compared to adults. The implications of this difference for inhalation dosimetry and children's risk assessment were evaluated in runs of the U.S. Environmental Protection Agency (U.S. EPA) 1994 reference concentration (RfC) methodology and the ICRP 1994 inhalation dosimetry model. Dosimetry estimates were made for 3-mo-old children and adults for particles and Category 1 and 2 reactive gases in the following respiratory-tract regions: extrathoracic (ET), tracheobronchial (BB), bronchioles (bb), and pulmonary (PU). Systemic dosimetry estimates were made for nonreactive (Category 3) gases. Results suggest similar ET dosimetry for children and adults for all types of inhaled materials. BB dosimetry was also similar across age groups except that the dosimetry of ultrafine particles in this region was twofold greater in 3-mo-old children than in adults. In contrast, the bb region generally showed higher dosimetry of particles and gases in adults than in children. Particle dose in the PU region was two- to fourfold higher in 3-mo-old children, with the greatest child/adult difference occurring for submicron size particles. Particulate dosimetry estimates with the default RfC methodology were below those found with the ICRP model for both adults and children for submicrometer sized particles. There were no cases in which reactive gas dosimetry was substantially greater in the respiratory regions of 3-mo-old children. Estimates of systemic dose of Category 3 gases were greater in 3-mo-old children than in adults, especially for liver dose of metabolite for rapidly metabolized gases. These analyses support the approach of assuming twofold greater inhalation dose in children than adults, although there are cases in which this differential can be greater and others where it can be less.

The refinement of uncertainty/safety factors in risk assessment by the incorporation of data on toxicokinetic variability in humans

The derivation of safe levels of exposure in humans for compounds that are assumed to cause threshold toxicity has relied on the application of a 100-fold uncertainty factor to a measure for the threshold, such as the no observed adverse effect level (NOAEL) or the benchmark dose (BMD). This 100-fold safety factor consists of the product of two 10-fold factors allowing for human variability and interspecies differences. The International Programme on Chemical Safety has suggested the subdivision of these 10-fold factors to allow for variability in toxicokinetics and toxicodynamics. This subdivision allows the replacement of the default uncertainty factors with a chemical-specific adjustment factor (CSAF) when suitable data are available. This short review describes potential options to refine safety factors used in risk assessment, with particular emphasis on pathway-related uncertainty factors associated with variability in kinetics. These pathway-related factors were derived from a database that quantified interspecies differences and human variability in phase I metabolism, phase II metabolism, and renal excretion. This approach allows metabolism and pharmacokinetic data in healthy adults and subgroups of the population to be incorporated in the risk-assessment process and constitutes an intermediate approach between simple default factors and chemical-specific adjustment factors.

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

This review provides an account of recent developments arising from a database that defined human variability in phase I metabolism (CYP1A2, CYP2A6, CYP2C9, CYP2C19, CYP2D6, CYP2E1, CYP3A4, hydrolysis, alcohol dehydrogenase), phase II metabolism (N-acetyltransferases, glucuronidation, glycine conjugation, sulphation) and renal excretion. This database was used to derive pathway-related uncertainty factors for chemical risk assessment that allow for human variability in toxicokinetics. Probe substrates for each pathway of elimination were selected on the basis that oral absorption was >95% and that the metabolic route was the primary route of elimination of the compound (60-100% of a dose). Intravenous data were used for compounds for which absorption was variable. Human variability in kinetics was quantified for each compound from published pharmacokinetic studies (after oral and intravenous dosing) in healthy adults and other subgroups of the population using parameters relating to chronic exposure (metabolic and total clearances, area under the plasma concentration-time curve (AUC)) and acute exposure (Cmax) (data not presented here). The pathway-related uncertainty factors were calculated to cover 95%, 97.5% and 99% of the population of healthy adults and of each subgroup. Pathway-related uncertainty factors allow metabolism data to be incorporated into the derivation of health-based guidance values. They constitute an intermediate approach between the general kinetic default factors (3.16) and a chemical-specific adjustment factor. Applications of pathway-related uncertainty factors for chemical risk assessment and future refinements of the approach are discussed. A knowledge-based framework to predict human variability in kinetics for xenobiotics showing a threshold dose below which toxic effects are not observed, is proposed to move away from default assumptions.

Risk-benefit analysis of micronutrients

Traditionally, different approaches have been used to determine the recommended dietary allowances for micronutrients, above which there is a low risk of deficiency, and safe upper levels, below which there is a negligible risk of toxicity. The advice given to risk managers has been in the form of point estimates, such as the recommended dietary allowance (RDA) and the tolerable upper level (UL). In future, the gap between the two intake-response curves may become narrower, as more sensitive indicators of deficiency and toxicity are used, and as health benefits above the recommended daily allowance are taken into account. This paper reviews the traditional approaches and proposes a novel approach to compare beneficial and adverse effects across intake levels. This model can provide advice for risk managers in a form that will allow the risk of deficiency or the risk of not experiencing the benefit to be weighed against the risk of toxicity. The model extends the approach used to estimate recommended dietary allowances to make it applicable to both beneficial and adverse effects and to extend the intake-incidence data to provide a range of estimates that can be considered by the risk manager. The data-requirements of the model are the incidence of a response at one or more levels of intake, and a suitable coefficient of variation to represent the person-to-person variations within the human population. A coefficient of variation of 10% or 15% has been used for established recommended dietary allowances and a value of 15% is proposed as default for considerations of benefit. A coefficient of variation of 45% is proposed as default for considerations of toxicity, based on analyses of human variability in the fate and effects of therapeutic drugs. Using this approach risk managers, working closely with risk assessors, will be able to define ranges of intake based on a balance between the risks of deficiency (or lack of benefit) and toxicity.

Impact of inter-individual differences in drug metabolism and pharmacokinetics on safety evaluation

Safety evaluation aims to assess the dose-response relationship to determine a dose/level of exposure for food contaminants below which no deleterious effect is measurable that is 'without appreciable health risk' when consumed daily over a lifetime. These safe levels, such as the acceptable daily intake (ADI) have been derived from animal studies using surrogates for the threshold such as the no-observed-adverse-effect-level (NOAEL). The extrapolation from the NOAEL to the human safe intake uses a 100-fold uncertainty factor, defined as the product of two 10-fold factors allowing for human variability and interspecies differences. The 10-fold factor for human variability has been further subdivided into two factors of 10(0.5) (3.16) to cover toxicokinetics and toxicodynamics and this subdivsion allows for the replacement of an uncertainty factor with a chemical-specific adjustment factor (CSAF) when compound-specific data are available. Recently, an analysis of human variability in pharmacokinetics for phase I metabolism (CYP1A2, CYP2A6, CYP2C9, CYP2C19, CYP2D6, CYP2E1, CYP3A4, hydrolysis, alcohol dehydrogenase), phase II metabolism (N-acetyltransferase, glucuronidation, glycine conjugation, sulphation) and renal excretion was used to derive pathway-related uncertainty factors in subgroups of the human population (healthy adults, effects of ethnicity and age). Overall, the pathway-related uncertainty factors (99th centile) were above the toxicokinetic uncertainty factor for healthy adults exposed to xenobiotics handled by polymorphic metabolic pathways (and assuming the parent compound was the proximate toxicant) such as CYP2D6 poor metabolizers (26), CYP2C19 poor metabolizers (52) and NAT-2 slow acetylators (5.2). Neonates were the most susceptible subgroup of the population for pathways with available data [CYP1A2 and glucuronidation (12), CYP3A4 (14), glycine conjugation (28)]. Data for polymorphic pathways were not available in neonates but uncertainty factors of up to 45 and 9 would allow for the variability observed in children for CYP2D6 and CYP2C19 metabolism, respectively. This review presents an overview on the history of uncertainty factors, the main conclusions drawn from the analysis of inter-individual differences in metabolism and pharmacokinetics, the development of pathway-related uncertainty factors and their use in chemical risk assessment.