A UK framework for the assessment and integration of different scientific evidence streams in chemical risk assessment.

BACKGROUND: Few methods are available for transparently combining different evidence streams for chemical risk assessment to reach an integrated conclusion on the probability of causation. Hence, the UK Committees on Toxicity (COT) and on Carcinogenicity (COC) have reviewed current practice and developed guidance on how to achieve this in a transparent manner, using graphical visualisation. METHODS/APPROACH: All lines of evidence, including toxicological, epidemiological, new approach methodologies, and mode of action should be considered, taking account of their strengths/weaknesses in their relative weighting towards a conclusion on the probability of causation. A qualitative estimate of the probability of causation is plotted for each line of evidence and a combined estimate provided. DISCUSSION/CONCLUSIONS: Guidance is provided on integration of multiple lines of evidence for causation, based on current best practice. Qualitative estimates of probability for each line of evidence are plotted graphically. This ensures a deliberative, consensus conclusion on likelihood of causation is reached. It also ensures clear communication of the influence of the different lines of evidence on the overall conclusion on causality. Issues on which advice from the respective Committees is sought varies considerably, hence the guidance is designed to be sufficiently flexible to meet this need.

PBTK modelling platforms and parameter estimation tools to enable animal-free risk assessment: recommendations from a joint EPAA--EURL ECVAM ADME workshop.

Information on toxicokinetics is critical for animal-free human risk assessment. Human external exposure must be translated into human tissue doses and compared with in vitro actual cell exposure associated to effects (in vitro-in vivo comparison). Data on absorption, distribution, metabolism and excretion in humans (ADME) could be generated using in vitro and QSAR tools. Physiologically-based toxicokinetic (PBTK) computer modelling could serve to integrate disparate in vitro and in silico findings. However, there are only few freely-available PBTK platforms currently available. And although some ADME parameters can be reasonably estimated in vitro or in silico, important gaps exist. Examples include unknown or limited applicability domains and lack of (high-throughput) tools to measure penetration of barriers, partitioning between blood and tissues and metabolic clearance. This paper is based on a joint EPAA--EURL ECVAM expert meeting. It provides a state-of-the-art overview of the availability of PBTK platforms as well as the in vitro and in silico methods to parameterise basic (Tier 1) PBTK models. Five high-priority issues are presented that provide the prerequisites for wider use of non-animal based PBTK modelling for animal-free chemical risk assessment.

A novel method to derive a human safety limit for PFOA by gene expression profiling and modelling.

Perfluorooctanoic acid (PFOA) is a persistent environmental contaminant that can accumulate in the human body due to its long half-life. This substance has been associated with liver, pancreatic, testicular and breast cancers, liver steatosis and endocrine disruption. PFOA is a member of a large group of substances also known as "forever chemicals" and the vast majority of substances of this group lack toxicological data that would enable their effective risk assessment in terms of human health hazards. This study aimed to derive a health-based guidance value for PFOA intake (ng/kg BW/day) from in vitro transcriptomics data. To this end, we developed an in silico workflow comprising five components: (i) sourcing in vitro hepatic transcriptomics concentration-response data; (ii) deriving molecular points of departure using BMDExpress3 and performing pathway analysis using gene set enrichment analysis (GSEA) to identify the most sensitive molecular pathways to PFOA exposure; (iii) estimating freely-dissolved PFOA concentrations in vitro using a mass balance model; (iv) estimating in vivo doses by reverse dosimetry using a PBK model for PFOA as part of a quantitative in vitro to in vivo extrapolation (QIVIVE) algorithm; and (v) calculating a tolerable daily intake (TDI) for PFOA. Fourteen percent of interrogated genes exhibited in vitro concentration-response relationships. GSEA pathway enrichment analysis revealed that "fatty acid metabolism" was the most sensitive pathway to PFOA exposure. In vitro free PFOA concentrations were calculated to be 2.9% of the nominal applied concentrations, and these free concentrations were input into the QIVIVE workflow. Exposure doses for a virtual population of 3,000 individuals were estimated, from which a TDI of 0.15 ng/kg BW/day for PFOA was calculated using the benchmark dose modelling software, PROAST. This TDI is comparable to previously published values of 1.16, 0.69, and 0.86 ng/kg BW/day by the European Food Safety Authority. In conclusion, this study demonstrates the combined utility of an "omics"-derived molecular point of departure and in silico QIVIVE workflow for setting health-based guidance values in anticipation of the acceptance of in vitro concentration-response molecular measurements in chemical risk assessment.

The probable future of toxicology - probabilistic risk assessment.

Both because of the shortcomings of existing risk assessment methodologies, as well as newly available tools to predict hazard and risk with machine learning approaches, there has been an emerging emphasis on probabilistic risk assessment. Increasingly sophisticated AI models can be applied to a plethora of exposure and hazard data to obtain not only predictions for particular endpoints but also to estimate the uncertainty of the risk assessment outcome. This provides the basis for a shift from deterministic to more probabilistic approaches but comes at the cost of an increased complexity of the process as it requires more resources and human expertise. There are still challenges to overcome before a probabilistic paradigm is fully embraced by regulators. Based on an earlier white paper (Maertens et al., 2022), a workshop discussed the prospects, challenges and path forward for implementing such AI-based probabilistic hazard assessment. Moving forward, we will see the transition from categorized into probabilistic and dose-dependent hazard outcomes, the application of internal thresholds of toxicological concern for data-poor substances, the acknowledgement of user-friendly open-source software, a rise in the expertise of toxicologists required to understand and interpret artificial intelligence models, and the honest communication of uncertainty in risk assessment to the public.

Probabilistic risk assessment, initially from engineering, is applied in toxicology to understand chemical-related hazards and their consequences. In toxicology, uncertainties abound ­ unclear molecular events, varied proposed outcomes, and population-level assessments for issues like neurodevelopmental disorders. Establishing links between chemical exposures and diseases, especially rare events like birth defects, often demands extensive studies. Existing methods struggle with subtle effects or those affecting specific groups. Future risk assessments must address developmental disease origins, presenting challenges beyond current capabilities. The intricate nature of many toxicological processes, lack of consensus on mechanisms and outcomes, and the need for nuanced population-level assessments highlight the complexities in understanding and quantifying risks associated with chemical exposures in the field of toxicology.

Refinement and calibration of a human PBPK model for the plasticiser, Di-(2-propylheptyl) phthalate (DPHP) using in silico, in vitro and human biomonitoring data.

An existing physiologically based pharmacokinetic model for Di-(2-propylheptyl) phthalate (DPHP) was refined to improve the simulations of the venous blood concentrations of the primary monoester metabolite, mono-(2-propylheptyl) phthalate (MPHP). This was considered a significant deficiency that should be addressed because the primary metabolite of other high molecular weight phthalates has been associated with toxicity. The various processes that influence the concentration of DPHP and MPHP in blood were re-evaluated and modified. A few simplifications of the existing model were made, including the removal of enterohepatic recirculation (EHR) of MPHP. However, the primary development was describing the partial binding of MPHP to plasma proteins following uptake of DPHP and metabolism in the gut affording better simulation of the trends observed in the biological monitoring data. Secondly, the relationship between blood concentrations and the urinary excretion of secondary metabolites was explored further because the availability of two data streams provides a better understanding of the kinetics than reliance on just one. Most human studies are conducted with few volunteers and generally with the absence of blood metabolite measurements which would likely imply an incomplete understanding of the kinetics. This has important implications for the "read across" approach proposed as part of the development of New Approach Methods for the replacement of animals in chemical safety assessments. This is where the endpoint of a target chemical is predicted by using data for the same endpoint from another more "data rich" source chemical. Validation of a model parameterized entirely with in vitro and in silico derived parameters and calibrated against several data streams would constitute a data rich source chemical and afford more confidence for future evaluations of other similar chemicals using the read-across approach.

Development, testing, parameterisation, and calibration of a human PBK model for the plasticiser, di (2-ethylhexyl) adipate (DEHA) using in silico, in vitro and human biomonitoring data.

Introduction: A physiologically based biokinetic model for di (2-ethylhexyl) adipate (DEHA) based on a refined model for di-(2-propylheptyl) phthalate (DPHP) was developed to interpret the metabolism and biokinetics of DEHA following a single oral dosage of 50 mg to two male and two female volunteers. Methods: The model was parameterized using in vitro and in silico methods such as, measured intrinsic hepatic clearance scaled from in vitro to in vivo and algorithmically predicted parameters such as plasma unbound fraction and tissue:blood partition coefficients (PCs). Calibration of the DEHA model was achieved using concentrations of specific downstream metabolites of DEHA excreted in urine. The total fractions of ingested DEHA eliminated as specific metabolites were estimated and were sufficient for interpreting the human biomonitoring data. Results: The specific metabolites of DEHA, mono-2-ethyl-5-hydroxyhexyl adipate (5OH-MEHA), mono-2-ethyl-5-oxohexyl adipate (5oxo-MEHA), mono-5-carboxy-2-ethylpentyl adipate (5cx-MEPA) only accounted for ∼0.45% of the ingested DEHA. Importantly, the measurements of adipic acid, a non-specific metabolite of DEHA, proved to be important in model calibration. Discussion: The very prominent trends in the urinary excretion of the metabolites, 5cx-MEPA and 5OH-MEHA allowed the important absorption mechanisms of DEHA to be modelled. The model should be useful for the study of exposure to DEHA of the general human population.

Development, testing, parameterisation, and calibration of a human PBPK model for the plasticiser, di-(2-ethylhexyl) terephthalate (DEHTP) using in silico, in vitro and human biomonitoring data.

A physiologically based pharmacokinetic model for di-(2-ethylhexyl) terephthalate (DEHTP) based on a refined model for di-(2-propylheptyl) phthalate (DPHP) was developed to interpret the metabolism and biokinetics of DEHTP following a single oral dose of 50 mg to three male volunteers. In vitro and in silico methods were used to generate parameters for the model. For example, measured intrinsic hepatic clearance scaled from in vitro to in vivo and plasma unbound fraction and tissue:blood partition coefficients (PCs) were predicted algorithmically. Whereas the development and calibration of the DPHP model was based upon two data streams, blood concentrations of parent chemical and first metabolite and the urinary excretion of metabolites, the model for DEHTP was calibrated against a single data stream, the urinary excretion of metabolites. Despite the model form and structure being identical significant quantitative differences in lymphatic uptake between the models were observed. In contrast to DPHP the fraction of ingested DEHTP entering lymphatic circulation was much greater and of a similar magnitude to that entering the liver with evidence for the dual uptake mechanisms discernible in the urinary excretion data. Further, the absolute amounts absorbed by the study participants, were much higher for DEHTP relative to DPHP. The in silico algorithm for predicting protein binding performed poorly with an error of more than two orders of magnitude. The extent of plasma protein binding has important implications for the persistence of parent chemical in venous blood-inferences on the behaviour of this class of highly lipophilic chemicals, based on calculations of chemical properties, should be made with extreme caution. Attempting read across for this class of highly lipophilic chemicals should be undertaken with caution since basic adjustments to PCs and metabolism parameters would be insufficient, even when the structure of the model itself is appropriate. Therefore, validation of a model parameterized entirely with in vitro and in silico derived parameters would need to be calibrated against several human biomonitoring data streams to constitute a data rich source chemical to afford confidence for future evaluations of other similar chemicals using the read-across approach.

Alternative (non-animal) methods for cosmetics testing: current status and future prospects-2010.

The 7th amendment to the EU Cosmetics Directive prohibits to put animal-tested cosmetics on the market in Europe after 2013. In that context, the European Commission invited stakeholder bodies (industry, non-governmental organisations, EU Member States, and the Commission's Scientific Committee on Consumer Safety) to identify scientific experts in five toxicological areas, i.e. toxicokinetics, repeated dose toxicity, carcinogenicity, skin sensitisation, and reproductive toxicity for which the Directive foresees that the 2013 deadline could be further extended in case alternative and validated methods would not be available in time. The selected experts were asked to analyse the status and prospects of alternative methods and to provide a scientifically sound estimate of the time necessary to achieve full replacement of animal testing. In summary, the experts confirmed that it will take at least another 7-9 years for the replacement of the current in vivo animal tests used for the safety assessment of cosmetic ingredients for skin sensitisation. However, the experts were also of the opinion that alternative methods may be able to give hazard information, i.e. to differentiate between sensitisers and non-sensitisers, ahead of 2017. This would, however, not provide the complete picture of what is a safe exposure because the relative potency of a sensitiser would not be known. For toxicokinetics, the timeframe was 5-7 years to develop the models still lacking to predict lung absorption and renal/biliary excretion, and even longer to integrate the methods to fully replace the animal toxicokinetic models. For the systemic toxicological endpoints of repeated dose toxicity, carcinogenicity and reproductive toxicity, the time horizon for full replacement could not be estimated.

Derivation of a Human In Vivo Benchmark Dose for Bisphenol A from ToxCast In Vitro Concentration Response Data Using a Computational Workflow for Probabilistic Quantitative In Vitro to In Vivo Extrapolation.

A computational workflow which integrates physiologically based kinetic (PBK) modelling; global sensitivity analysis (GSA), Approximate Bayesian Computation (ABC), Markov Chain Monte Carlo (MCMC) simulation and the Virtual Cell Based Assay (VCBA) for the estimation of the active, free in vitro concentration of chemical in the reaction medium was developed to facilitate quantitative in vitro to in vivo extrapolation (QIVIVE). The workflow was designed to estimate parameter and model uncertainty within a computationally efficient framework. The workflow was tested using a human PBK model for bisphenol A (BPA) and high throughput screening (HTS) in vitro concentration-response data, for estrogen and pregnane X receptor activation determined in human liver and kidney cell lines, from the ToxCast/Tox21 database. In vivo benchmark dose 10% lower confidence limits (BMDL10) for oral uptake of BPA (ng/kg BW/day) were calculated from the in vivo dose-responses and compared to the human equivalent dose (HED) BMDL10 for relative kidney weight change in the mouse derived by European Food Safety Authority (EFSA). Three from four in vivo BMDL10 values calculated in this study were similar to the EFSA values whereas the fourth was much smaller. The derivation of an uncertainty factor (UF) to accommodate the uncertainties associated with measurements using human cell lines in vitro, extrapolated to in vivo, could be useful for the derivation of Health Based Guidance Values (HBGV).

Development, Testing, Parameterisation and Calibration of a Human PBPK Model for the Plasticiser, Di-(2-propylheptyl) Phthalate (DPHP) Using in Silico, in vitro and Human Biomonitoring Data.

A physiologically based pharmacokinetic model for Di-(2-propylheptyl) phthalate (DPHP) was developed to interpret the biokinetics in humans after single oral doses. The model was parameterized with in vitro and in silico derived parameters and uncertainty and sensitivity analysis was used during the model development process to assess structure, biological plausibility and behaviour prior to simulation and analysis of human biological monitoring data. To provide possible explanations for some of the counter-intuitive behaviour of the biological monitoring data the model included a simple lymphatic uptake process for DPHP and enterohepatic recirculation (EHR) for DPHP and the mono ester metabolite mono-(2-propylheptyl) phthalate (MPHP). The model was used to simultaneously simulate the concentration-time profiles of blood DPHP, MPHP and the urinary excretion of two metabolites, mono-(2-propyl-6-hydroxyheptyl) phthalate (OH-MPHP) and mono-(2-propyl-6-carboxyhexyl) phthalate (cx-MPHP). The availability of blood and urine measurements permitted a more robust qualitative and quantitative investigation of the importance of EHR and lymphatic uptake. Satisfactory prediction of blood DPHP and urinary metabolites was obtained whereas blood MPHP was less satisfactory. However, the delayed peak of DPHP concentration relative to MPHP in blood and second order metabolites in urine could be explained as a result of three processes: 1) DPHP entering the systemic circulation from the lymph, 2) rapid and very high protein binding and 3) the efficiency of the liver in removing DPHP absorbed via the hepatic route. The use of sensitivity analysis is considered important in the evaluation of uncertainty around in vitro and in silico derived parameters. By quantifying their impact on model output sufficient confidence in the use of a model should be afforded. This approach could expand the use of PBPK models since parameterization with in silico techniques allows for rapid model development. This in turn could assist in reducing the use of animals in toxicological evaluations by enhancing the utility of "read across" techniques.

Derivation of a Human In Vivo Benchmark Dose for Perfluorooctanoic Acid From ToxCast In Vitro Concentration-Response Data Using a Computational Workflow for Probabilistic Quantitative In Vitro to In Vivo Extrapolation.

A computational workflow which integrates physiologically based kinetic (PBK) modeling, global sensitivity analysis (GSA), approximate Bayesian computation (ABC), and Markov Chain Monte Carlo (MCMC) simulation was developed to facilitate quantitative in vitro to in vivo extrapolation (QIVIVE). The workflow accounts for parameter and model uncertainty within a computationally efficient framework. The workflow was tested using a human PBK model for perfluorooctanoic acid (PFOA) and high throughput screening (HTS) in vitro concentration-response data, determined in a human liver cell line, from the ToxCast/Tox21 database. In vivo benchmark doses (BMDs) for PFOA intake (ng/kg BW/day) and drinking water exposure concentrations (µg/L) were calculated from the in vivo dose responses and compared to intake values derived by the European Food Safety Authority (EFSA). The intake benchmark dose lower confidence limit (BMDL5) of 0.82 was similar to 0.86 ng/kg BW/day for altered serum cholesterol levels derived by EFSA, whereas the intake BMDL5 of 6.88 was six-fold higher than the value of 1.14 ng/kg BW/day for altered antibody titer also derived by the EFSA. Application of a chemical-specific adjustment factor (CSAF) of 1.4, allowing for inter-individual variability in kinetics, based on biological half-life, gave an intake BMDL5 of 0.59 for serum cholesterol and 4.91 (ng/kg BW/day), for decreased antibody titer, which were 0.69 and 4.31 the EFSA-derived values, respectively. The corresponding BMDL5 for drinking water concentrations, for estrogen receptor binding activation associated with breast cancer, pregnane X receptor binding associated with altered serum cholesterol levels, thyroid hormone receptor α binding leading to thyroid disease, and decreased antibody titer (pro-inflammation from cytokines) were 0.883, 0.139, 0.086, and 0.295 ng/ml, respectively, with application of no uncertainty factors. These concentrations are 5.7-, 36-, 58.5-, and 16.9-fold lower than the median measured drinking water level for the general US population which is approximately, 5 ng/ml.

Recombinant expression of aryl hydrocarbon receptor for quantitative ligand-binding analysis.

Recombinant expression of the aryl hydrocarbon receptor (AhR) yields small amounts of ligand-binding-competent AhR. Therefore, Spodoptera frugiperda (Sf9) cells and baculovirus have been evaluated for high-level and functional expression of AhR. Rat and human AhR were expressed as soluble protein in significant amounts. Expression of ligand-binding-competent AhR was sensitive to the protein concentration of Sf9 extract, and coexpression of the chaperone p23 failed to affect the yield of functional ligand-binding AhR. The expression system yielded high levels of functional protein, with the ligand-binding capacity (Bmax) typically 20-fold higher than that obtained with rat liver cytosol. Quantitative estimates of the ligand-binding affinity of human and rat AhR were obtained; the Kd for recombinant rat AhR was indistinguishable from that of native rat AhR, thereby validating the expression system as a faithful model for native AhR. The human AhR bound TCDD with significantly lower affinity than the rat AhR. These findings demonstrate high-level expression of ligand-binding-competent AhR, and sufficient AhR for quantitative analysis of ligand binding.

Development, Testing, Parameterization, and Calibration of a Human Physiologically Based Pharmacokinetic Model for the Plasticizer, Hexamoll® Diisononyl-Cyclohexane-1, 2-Dicarboxylate Using In Silico, In Vitro, and Human Biomonitoring Data.

A physiologically based pharmacokinetic model for Hexamoll® diisononyl-cyclohexane-1, 2-dicarboxylate was developed to interpret the biokinetics in humans after single oral doses. The model was parameterized with in vitro and in silico derived parameters and uncertainty and sensitivity analysis was used during the model development process to assess structure, biological plausibility and behavior prior to simulation and analysis of human biological monitoring data. The model provided good simulations of the urinary excretion (Curine) of two metabolites; cyclohexane-1,2-dicarboxylic acid mono hydroxyisononyl ester (OH-MINCH) and cyclohexane-1, 2-dicarboxylic acid mono carboxyisononyl ester (cx-MINCH) from the biotransformation of mono-isononyl-cyclohexane-1, 2-dicarboxylate (MINCH), the monoester metabolite of di-isononyl-cyclohexane-1,2-dicarboxylate. However, good simulations could be obtained, with and without, a lymphatic compartment. Selection of an appropriate model structure was informed by sensitivity analysis which could identify and quantify the contribution to variability in Curine by parameters, such as, the fraction of oral dose that directly entered the lymphatic compartment and therefore by-passed the liver and the fraction of MINCH bio-transformed to cx-MINCH and OH-MINCH. By constraining these parameters within biologically plausible limits the presence of a lymphatic compartment was deemed an important component of model structure. Furthermore, the use of sensitivity analysis is important in the evaluation of uncertainty around in silico derived parameters. By quantifying their impact on model output sufficient confidence in the use of a model should be afforded. This type of approach could expand the use of physiologically based pharmacokinetic models since parameterization with in silico techniques allows for rapid model development. This in turn could assist in reducing the use of animals in toxicological evaluations by enhancing the utility of "read across" techniques.

Use cases, best practice and reporting standards for metabolomics in regulatory toxicology.

Metabolomics is a widely used technology in academic research, yet its application to regulatory science has been limited. The most commonly cited barrier to its translation is lack of performance and reporting standards. The MEtabolomics standaRds Initiative in Toxicology (MERIT) project brings together international experts from multiple sectors to address this need. Here, we identify the most relevant applications for metabolomics in regulatory toxicology and develop best practice guidelines, performance and reporting standards for acquiring and analysing untargeted metabolomics and targeted metabolite data. We recommend that these guidelines are evaluated and implemented for several regulatory use cases.

A Computational Workflow for Probabilistic Quantitative in Vitro to in Vivo Extrapolation.

A computational workflow was developed to facilitate the process of quantitative in vitro to in vivo extrapolation (QIVIVE), specifically the translation of in vitro concentration-response to in vivo dose-response relationships and subsequent derivation of a benchmark dose value (BMD). The workflow integrates physiologically based pharmacokinetic (PBPK) modeling; global sensitivity analysis (GSA), Approximate Bayesian Computation (ABC) and Markov Chain Monte Carlo (MCMC) simulation. For a given set of in vitro concentration and response data the algorithm returns the posterior distribution of the corresponding in vivo, population-based dose-response values, for a given route of exposure. The novel aspect of the workflow is a rigorous statistical framework for accommodating uncertainty in both the parameters of the PBPK model (both parameter uncertainty and population variability) and in the structure of the PBPK model itself recognizing that the model is an approximation to reality. Both these sources of uncertainty propagate through the workflow and are quantified within the posterior distribution of in vivo dose for a fixed representative in vitro concentration. To demonstrate this process and for comparative purposes a similar exercise to previously published work describing the kinetics of ethylene glycol monoethyl ether (EGME) and its embryotoxic metabolite methoxyacetic acid (MAA) in rats was undertaken. The computational algorithm can be used to extrapolate from in vitro data to any organism, including human. Ultimately, this process will be incorporated into a user-friendly, freely available modeling platform, currently under development, that will simplify the process of QIVIVE.

Toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin in the developing male Wistar(Han) rat. II: Chronic dosing causes developmental delay.

We have investigated whether fetal exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) causes defects in the male reproductive system of the rat using chronically exposed rats to ensure continuous exposure of the fetus. Five- to six-week-old rats were exposed to control diet, or diet containing TCDD, to attain an average dose of 2.4, 8, and 46 ng TCDD/kg/day for 12 weeks, whereupon the rats were mated and allowed to litter; rats were switched to control diet after parturition. Male offsprings were allowed to develop until kills on PND70 (25 per group) or PND120 (all remaining animals). Offspring from the high-dose group showed an increase in total litter loss, and the number of animals alive on postnatal day (PND)4 in the high-dose group was approximately 26% less than control. The high and medium dose offsprings showed decreased weights at various ages. Balano-preputial separation (BPS) was significantly delayed in all three dose groups compared to control. There were no significant effects of maternal treatment when the offsprings were subjected to a functional observational battery or learning tests, with the exception that the high-dose group showed a deficit in motor activity. Twenty rats per group were mated to females, and there were no significant effects of maternal treatment on the fertility of these rats or on the F1 or F2 sex ratio. Sperm parameters at PND70 and 120 showed no significant effect of maternal treatment, with the exception that there was an increase in the proportion of abnormal sperm in the high-dose group at PND70; this is associated with the developmental delay in puberty in this dose group. There were no remarkable findings of maternal treatment on organ weights, with the exception that testis weights were reduced by approximately 10% at PND70 (but not PND120), and although the experiment was sufficiently powered to detect small changes, ventral prostate weight was not reduced. There were no significant effects of maternal treatment upon histopathological comparison of high-dose and control group organs. These data confirm that developmental exposure to TCDD shows no potent effect on adult sperm parameters or accessory sexual organs, but show that delay in BPS occurs after exposure to low doses of TCDD, and this is dependent upon whether TCDD is administered acutely or chronically.

Toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin in the developing male Wistar(Han) rat. I: No decrease in epididymal sperm count after a single acute dose.

It has been reported that fetal exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) causes defects in the male reproductive system of the rat. We set out to replicate and extend these effects using a robust experimental design. Groups of 75 (control vehicle) or 55 (50, 200, or 1000 ng of TCDD/kg bodyweight) female Wistar(Han) rats were exposed to TCDD on gestational day (GD)15, then allowed to litter. The high-dose group dams showed no sustained weight loss compared to control, but four animals had total litter loss. Pups in the high-dose group showed reduced body weight up till day 21, and pups in the medium dose group showed reduced body weight in the first week postpartum. Balano-preputial separation was significantly delayed in the high-dose group male offspring. There were no significant effects of treatment when the offspring were subjected to a functional observational battery or mated with females to assess reproductive capability. Twenty-five males per group were killed on postnatal day (PND) 70, and approximately 60 animals per group (approximately 30 for the high-dose group) on PND120 to assess seminology and other end points. At PND120, the two highest dose groups showed a statistically significant elevation of sperm counts, compared to control; however, this effect was small (approximately 30%), within the normal range of sperm counts for this strain of rat, was not reflected in testicular spermatid counts nor PND70 data, and is therefore postulated to have no biological significance. Although there was an increase in the proportion of abnormal sperm at PND70, seminology parameters were otherwise unremarkable. Testis weights in the high-dose group were slightly decreased at PND70 and 120, and at PND120, brain weights were decreased in the high-dose group, liver to body weight ratios were increased for all three dose groups, with an increase in inflammatory cell foci in the epididymis in the high-dose group. These data show that TCDD is a potent developmental toxin after exposure of the developing fetus but that acute developmental exposure to TCDD on GD15 caused no decrease in sperm counts.

Relationships between tissue levels of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), mRNAs, and toxicity in the developing male Wistar(Han) rat.

We compared the effects of a single acute dose, or chronic fetal exposure, to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) on the male reproductive system of the Wistar(Han) rat. Tissue samples were taken from dams on gestation day (GD)16 and GD21, and from offspring on postnatal days (PND)70 and 120. Steady-state concentration of TCDD was demonstrated in the chronic study: body burdens were comparable in both studies. Fetal TCDD concentrations were comparable after acute and chronic exposure, and demonstrate more potent toxicity after chronic versus acute dosing. In maternal liver, cytochrome P450 (CYP)1A1 and CYP1A2 RNA were induced. In fetus, there was induction of both CYP1A1 and CYP1A2 RNA at medium and high doses, but inadequate evidence for induction at low dose in either study. The low level induction of CYP1A1 RNA at low dose in fetus argues against AhR activation in fetus as a mechanism of toxicity of TCDD in causing delay in balanopreputial separation (BPS), and the greater induction of CYP1A1 RNA in PND70 offspring liver from chronically-dosed dams suggests that lactational transfer of TCDD is crucial to this toxicity. These data characterize the maternal and fetal disposition of TCDD, induction of CYP1A1 RNA as a measure of AhR activation, and suggest that lactational transfer of TCDD determines the difference in delay in BPS between the two studies.

Characterizing uncertainty and variability in physiologically based pharmacokinetic models: state of the science and needs for research and implementation.

Physiologically based pharmacokinetic (PBPK) models are used in mode-of-action based risk and safety assessments to estimate internal dosimetry in animals and humans. When used in risk assessment, these models can provide a basis for extrapolating between species, doses, and exposure routes or for justifying nondefault values for uncertainty factors. Characterization of uncertainty and variability is increasingly recognized as important for risk assessment; this represents a continuing challenge for both PBPK modelers and users. Current practices show significant progress in specifying deterministic biological models and nondeterministic (often statistical) models, estimating parameters using diverse data sets from multiple sources, using them to make predictions, and characterizing uncertainty and variability of model parameters and predictions. The International Workshop on Uncertainty and Variability in PBPK Models, held 31 Oct-2 Nov 2006, identified the state-of-the-science, needed changes in practice and implementation, and research priorities. For the short term, these include (1) multidisciplinary teams to integrate deterministic and nondeterministic/statistical models; (2) broader use of sensitivity analyses, including for structural and global (rather than local) parameter changes; and (3) enhanced transparency and reproducibility through improved documentation of model structure(s), parameter values, sensitivity and other analyses, and supporting, discrepant, or excluded data. Longer-term needs include (1) theoretical and practical methodological improvements for nondeterministic/statistical modeling; (2) better methods for evaluating alternative model structures; (3) peer-reviewed databases of parameters and covariates, and their distributions; (4) expanded coverage of PBPK models across chemicals with different properties; and (5) training and reference materials, such as cases studies, bibliographies/glossaries, model repositories, and enhanced software. The multidisciplinary dialogue initiated by this Workshop will foster the collaboration, research, data collection, and training necessary to make characterizing uncertainty and variability a standard practice in PBPK modeling and risk assessment.