Comparison of Rat Hepatocyte 2D-Monocultures and Hepatocytes Non-Parenchymal Cell Co-Cultures for Assessing Chemical Toxicity

Liver responses are the most common endpoints used as the basis for setting exposure standards. Liver hepatocytes play a vital role in biotransformation of xenobiotics, but non-parenchymal cells (NPCs) in the liver are also involved in certain liver responses. Development of in vitro systems that more faithfully capture liver responses to reduce reliance on animals is a major focus of New Approach Methodology (NAMs). Since rodent regulatory studies are frequently the sole source safety assessment data, mode-of-action data, and used for risk assessments, in vitro rodent models that reflect in vivo responses need to be developed to reduce reliance on animal models. In the work presented in this paper, we developed a 2-D hepatocyte monoculture and 2-D liver cell co-culture system using rat liver cells. These models were assessed for conditions for short-term stability of the cultures and phenotypic and transcriptomic responses of 2 prototypic hepatotoxicants compounds - acetaminophen and phenobarbital. The optimized multi-cellular 2-D culture required use of freshly prepared hepatocytes and NPCs from a single rat, a 3:1 ratio of hepatocytes to NPCs and growth medium using 50% Complete Williams E medium (WEM) and 50% Endothelial Cell Medium (ECM). The transcriptomic responses of the 2 model systems to PB were compared to previous studies from TG-Gates on the gene expression changes in intact rats and the co-culture model responses were more representative of the in vivo responses. Transcriptomic read-outs promise to move beyond conventional phenotypic evaluations with these in vitro NAMs and provide insights about modes of action.

Estimating provisional margins of exposure for data-poor chemicals using high-throughput computational methods

Current computational technologies hold promise for prioritizing the testing of the thousands of chemicals in commerce. Here, a case study is presented demonstrating comparative risk-prioritization approaches based on the ratio of surrogate hazard and exposure data, called margins of exposure (MoEs). Exposures were estimated using a U.S. EPA’s ExpoCast predictive model (SEEM3) results and estimates of bioactivity were predicted using: 1) Oral equivalent doses (OEDs) derived from U.S. EPA’s ToxCast high-throughput screening program, together with in vitro to in vivo extrapolation and 2) thresholds of toxicological concern (TTCs) determined using a structure-based decision-tree using the Toxtree open source software. To ground-truth these computational approaches, we compared the MoEs based on predicted noncancer TTC and OED values to those derived using the traditional method of deriving points of departure from no-observed adverse effect levels (NOAELs) from in vivo oral exposures in rodents. TTC-based MoEs were lower than NOAEL-based MoEs for 520 out of 522 (99.6%) compounds in this smaller overlapping dataset, but were relatively well correlated with the same (r² = 0.59). TTC-based MoEs were also lower than OED-based MoEs for 590 (83.2%) of the 709 evaluated chemicals, indicating that TTCs may serve as a conservative surrogate in the absence of chemical-specific experimental data. The TTC-based MoE prioritization process was then applied to over 45,000 curated environmental chemical structures as a proof-of-concept for high-throughput prioritization using TTC-based MoEs. This study demonstrates the utility of exploiting existing computational methods at the pre-assessment phase of a tiered risk-based approach to quickly, and conservatively, prioritize thousands of untested chemicals for further study.

Considerations for Improving Metabolism Predictions for In Vitro to In Vivo Extrapolation

High-throughput (HT) in vitro to in vivo extrapolation (IVIVE) is an integral component in new approach method (NAM)-based risk assessment paradigms, for rapidly translating in vitro toxicity assay results into the context of in vivo exposure. When coupled with rapid exposure predictions, HT-IVIVE supports the use of HT in vitro assays for risk-based chemical prioritization. However, the reliability of prioritization based on HT bioactivity data and HT-IVIVE can be limited as the domain of applicability of current HT-IVIVE is generally restricted to intrinsic clearance measured primarily in pharmaceutical compounds. Further, current approaches only consider parent chemical toxicity. These limitations occur because current state-of-the-art HT prediction tools for clearance and metabolite kinetics do not provide reliable data to support HT-IVIVE. This paper discusses current challenges in implementation of IVIVE for prioritization and risk assessment and recommends a path forward for addressing the most pressing needs and expanding the utility of IVIVE.

Time-dependent genomic response in primary human uroepithelial cells exposed to arsenite for up to 60 days

Evidence from both in vivo and in vitro studies suggests that gene expression changes from long-term exposure to arsenite evolve markedly over time, including reversals in the direction of expression change in key regulatory genes. In this study, human uroepithelial cells from the ureter segments of 4 kidney-donors were continuously treated in culture with arsenite at concentrations of 0.1 or 1 μM for 60 days. Gene expression at 10, 20, 30, 40, and 60 days was determined using Affymetrix human genome microarrays and signal pathway analysis was performed using GeneGo Metacore. Arsenic treated cells continued to proliferate for the full 60-day period, whereas untreated cells ceased proliferating after approximately 30 days. A peak in the number of gene changes in the treated cells compared to untreated controls was observed between 30 and 40 days of exposure, with substantially fewer changes at 10 and 60 days, suggesting remodeling of the cells over time. Consistent with this possibility, the direction of expression change for a number of key genes was reversed between 20 and 30 days, including CFOS and MDM2. While the progression of gene changes was different for each subject, a common pattern was observed in arsenic treated cells over time, with early upregulation of oxidative stress responses (HMOX1, NQ01, TXN, TXNRD1) and down-regulation of immune/inflammatory responses (IKKα). At around 30 days, there was a transition to increased inflammatory and proliferative signaling (AKT, CFOS), evidence of epithelial-to-mesenchymal transition (EMT), and alterations in DNA damage responses (MDM2, ATM). A common element in the changing response of cells to arsenite over time appears to involve up-regulation of MDM2 by inflammatory signaling (through AP-1 and NF-κB), leading to inhibition of P53 function.

Development of a physiologically based pharmacokinetic model of diisononyl phthalate (DiNP) in pregnant rat and human

A physiologically based pharmacokinetic (PBPK) model for di-isononyl phthalate (DiNP) was developed by adapting the existing models for di(2-ethylhexyl) phthalate (DEHP) and di-butylphthalate (DBP). Both pregnant rat and human time-course plasma and urine data were used to address the hydrolysis of DiNP in intestinal tract, plasma, and liver as well as hepatic oxidative metabolism and conjugation of the monoester and primary oxidative metabolites. Data in both rats and humans were available to inform the uptake and disposition of mono-isononyl phthalate (MiNP) as well as the three primary oxidative metabolites including hydroxy (7-OH)-, oxo (7-OXO)-, and carboxy (7-COX)-monoisononyl phthalate in plasma and urine. The DiNP model was reliable over a wide range of exposure levels in the pregnant rat as well as the two low exposure levels in humans including capturing the nonlinear behavior in the pregnant rat after repeated 750 mg/kg/day dosing. The presented DiNP PBPK model in pregnant rat and human, based upon an extensive kinetic dataset in both species, may provide a basis for assessing human equivalent exposures based upon either rodent or in vitro points of departure.

The role of fit-for-purpose assays within tiered testing approaches: A case study evaluating prioritized estrogen-active compounds in an in vitro human uterotrophic assay

Chemical risk assessment relies on toxicity tests that require significant numbers of animals, time and costs. For the >30,000 chemicals in commerce, the current scale of animal testing is insufficient to address chemical safety concerns as regulatory and product stewardship considerations evolve to require more comprehensive understanding of potential biological effects, conditions of use, and associated exposures. We demonstrate the use of a multi-level new approach methodology (NAMs) strategy for hazard- and risk-based prioritization to reduce animal testing. A Level 1/2 chemical prioritization based on estrogen receptor (ER) activity and metabolic activation using ToxCast data was used to select 112 chemicals for testing in a Level 3 human uterine cell estrogen response assay (IKA assay). The Level 3 data were coupled with quantitative in vitro to in vivo extrapolation (Q-IVIVE) to support bioactivity determination (as a surrogate for hazard) in a tissue-specific context. Assay AC₅₀s and Q-IVIVE were used to estimate human equivalent doses (HEDs), and HEDs were compared to rodent uterotrophic assay in vivo-derived points of departure (PODs). For substances active both in vitro and in vivo, IKA assay-derived HEDs were lower or equivalent to in vivo PODs for 19/23 compounds (83%). Activity exposure relationships were calculated, and the IKA assay was as or more protective of human health than the rodent uterotrophic assay for all IKA-positive compounds. This study demonstrates the utility of biologically relevant fit-for-purpose assays and supports the use of a multi-level strategy for chemical risk assessment.

A Qualitative Modeling Approach for Whole Genome Prediction Using High-Throughput Toxicogenomics Data and Pathway-Based Validation

Efficient high-throughput transcriptomics (HTT) tools promise inexpensive, rapid assessment of possible biological consequences of human and environmental exposures to tens of thousands of chemicals in commerce. HTT systems have used relatively small sets of gene expression measurements coupled with mathematical prediction methods to estimate genome-wide gene expression and are often trained and validated using pharmaceutical compounds. It is unclear whether these training sets are suitable for general toxicity testing applications and the more diverse chemical space represented by commercial chemicals and environmental contaminants. In this work, we built predictive computational models that inferred whole genome transcriptional profiles from a smaller sample of surrogate genes. The model was trained and validated using a large scale toxicogenomics database with gene expression data from exposure to heterogeneous chemicals from a wide range of classes (the Open TG-GATEs data base). The method of predictor selection was designed to allow high fidelity gene prediction from any pre-existing gene expression data set, regardless of animal species or data measurement platform. Predictive qualitative models were developed with this TG-GATES data that contained gene expression data of human primary hepatocytes with over 941 samples covering 158 compounds. A sequential forward search-based greedy algorithm, combining different fitting approaches and machine learning techniques, was used to find an optimal set of surrogate genes that predicted differential expression changes of the remaining genome. We then used pathway enrichment of up-regulated and down-regulated genes to assess the ability of a limited gene set to determine relevant patterns of tissue response. In addition, we compared prediction performance using the surrogate genes found from our greedy algorithm (referred to as the SV2000) with the landmark genes provided by existing technologies such as L1000 (Genometry) and S1500 (Tox21), finding better predictive performance for the SV2000. The ability of these predictive algorithms to predict pathway level responses is a positive step toward incorporating mode of action (MOA) analysis into the high throughput prioritization and testing of the large number of chemicals in need of safety evaluation.

Assessing bioactivity-exposure profiles of fruit and vegetable extracts in the BioMAP profiling system

The ToxCast program has generated in vitro screening data on over a thousand chemicals to assess potential disruption of important biological processes and assist in hazard identification and chemical testing prioritization. Few results have been reported for complex mixtures. To extend these ToxCast efforts to mixtures, we tested extracts from 30 organically grown fruits and vegetables in concentration-response in the BioMAP® assays. BioMAP systems use human primary cells primed with endogenous pathway activators to identify phenotypic perturbations related to proliferation, inflammation, immunomodulation, and tissue remodeling. Clustering of bioactivity profiles revealed separation of these produce extracts and ToxCast chemicals. Produce extracts elicited 87 assay endpoint responses per item compared to 20 per item for ToxCast chemicals. On a molar basis, the produce extracts were 10 to 50-fold less potent and when constrained to the maximum testing concentration of the ToxCast chemicals, the produce extracts did not show activity in as many assay endpoints. Using intake adjusted measures of dose, the bioactivity potential was higher for produce extracts than for agrichemicals, as expected based on the comparatively small amounts of agrichemical residues present on conventionally grown produce. The evaluation of BioMAP readouts and the dose responses for produce extracts showed qualitative and quantitative differences from results with single chemicals, highlighting challenges in the interpretation of bioactivity data and dose-response from complex mixtures.

Evaluating opportunities for advancing the use of alternative methods in risk assessment through the development of fit-for-purpose in vitro assays

An evolving regulatory, scientific, and legislative landscape is driving a fundamental change in how chemical safety decisions are made. As we move to implement changes, regulatory agencies and industry are beginning to adopt tiered approaches, which leverage high-throughput screening technologies for prioritization and read across, followed by interrogation of "hit chemicals" with more rigorous dose-response assessment either in fit-for-purpose human cell-based assays or with traditional in vivo tests. However, to date, suitable in vitro alternatives do not exist for the vast majority of the organ toxicities that form the basis of current regulatory decisions. To successfully support safety decisions, biologically relevant, quantitative, cell-based assays that evaluate dose-response and identify regions of safety for chemical exposure are required. This review evaluates the current state of the science in the development of such assays, identifies key gaps in the current tests, and recommends areas where research efforts may be focused to help move the risk assessment community towards more wide-spread use of in vitro methods. Our analysis suggests that a key shortcoming in the current efforts is the ability to test volatile compounds and to predict pulmonary toxicity. We present a mechanistically-based path forward for the development of a fit-for-purpose lung toxicity assay.

Contribution of ATM and ATR kinase pathways to p53-mediated response in etoposide and methyl methanesulfonate induced DNA damage

p53 is a key integrator of cellular response to DNA damage, supporting post-translational repair and driving transcription-mediated responses including cell cycle arrest, apoptosis, and repair. DNA damage sensing kinases recognize different types of DNA damage and initiate specific responses through various post-translational modifications of p53. This study evaluated chemical specificity of the p53 pathway response by manipulating p53 or its upstream kinases and assessing the effect on DNA damage and cellular responses to prototype chemicals: etoposide (ETP, topoisomerase II inhibitor) and methyl methane sulfonate (MMS, alkylating agent). p53-deficient cells demonstrated reduced accumulation of the p53 target proteins MDM2, p21, and Wip1; reduced apoptotic response; and increased DNA damage (p-H2AX and micronuclei) with both chemicals. However, p53 was not essential for cell cycle arrest in HT1080 or HCT116 cells. The two chemicals induced different patterns of kinase activation, particularly in terms of Chk 1, Chk 2, p38, and ERK 1/2. However, inhibition of the ATM pathway showed a greater effect on p53 activtation, apoptosis, and accumulation of DNA damage than ATR-Chk 1 or the MAP kinases regardless of the chemical used. These results indicate that ATM is the predominant upstream kinase responsible for activation of the p53-mediated DNA damage response for both MMS and ETP, though the downstream kinase response is markedly different. Environ. Mol. Mutagen. 58:72-83, 2017. © 2017 Wiley Periodicals, Inc.

Editor's Highlight: Development of an In vitro Assay Measuring Uterine-Specific Estrogenic Responses for Use in Chemical Safety Assessment

A toxicity pathway approach was taken to develop an in vitro assay using human uterine epithelial adenocarcinoma (Ishikawa) cells as a replacement for measuring an in vivo uterotrophic response to estrogens. The Ishikawa cell was determined to be fit for the purpose of recapitulating in vivo uterine response by verifying fidelity of the biological pathway components and the dose-response predictions to women of child-bearing age. Expression of the suite of estrogen receptors that control uterine proliferation (ERα66, ERα46, ERα36, ERβ, G-protein coupled estrogen receptor (GPER)) were confirmed across passages and treatment conditions. Phenotypic responses to ethinyl estradiol (EE) from transcriptional activation of ER-mediated genes, to ALP enzyme induction and cellular proliferation occurred at concentrations consistent with estrogenic activity in adult women (low picomolar). To confirm utility of this model to predict concentration-response for uterine proliferation with xenobiotics, we tested the concentration-response for compounds with known uterine estrogenic activity in humans and compared the results to assays from the ToxCast and Tox21 suite of estrogen assays. The Ishikawa proliferation assay was consistent with in vivo responses and was a more sensitive measure of uterine response. Because this assay was constructed by first mapping the key molecular events for cellular response, and then ensuring that the assay incorporated these events, the resulting cellular assay should be a reliable tool for identifying estrogenic compounds and may provide improved quantitation of chemical concentration response for in vitro-based safety assessments.

Editor's Highlight: Screening ToxCast Prioritized Chemicals for PPARG Function in a Human Adipose-Derived Stem Cell Model of Adipogenesis

The developmental origins of obesity hypothesis posits a multifaceted contribution of factors to the fetal origins of obesity and metabolic disease. Adipocyte hyperplasia in gestation and early childhood may result in predisposition for obesity later in life. Rodent in vitro and in vivo studies indicate that some chemicals may directly affect adipose progenitor cell differentiation, but the human relevance of these findings is unclear. The nuclear receptor peroxisome proliferator-activated receptor gamma (PPARG) is the master regulator of adipogenesis. Human adipose-derived stem cells (hASC) isolated from adipose tissue express endogenous isoforms of PPARG and represent a biologically relevant cell-type for evaluating activity of PPARG ligands. Here, a multi-endpoint approach based on a phenotypic adipogenesis assay was applied to screen a set of 60 chemical compounds identified in ToxCast Phase I as PPARG active (49) or inactive (11). Chemicals showing activity in the adipogenesis screen were further evaluated in a series of 4 orthogonal assays representing 7 different key events in PPARG-dependent adipogenesis, including gene transcription, protein expression, and adipokine secretion. An siRNA screen was also used to evaluate PPARG-dependence of the adipogenesis phenotype. A universal concentration-response design enabled inter-assay comparability and implementation of a weight-of-evidence approach for bioactivity classification. Collectively, a total of 14/49 (29%) prioritized chemicals were identified with moderate-to-strong activity for human adipogenesis. These results provide the first integrated screening approach of prioritized ToxCast chemicals in a human stem cell model of adipogenesis and provide insight into the capacity of PPARG-activating chemicals to modulate early life programming of adipose tissue.

Using Transcriptomics to Evaluate Thresholds in Genotoxicity Dose–Response

Several genotoxic chemicals have been reported to produce threshold-shaped dose–response curves for mutation and genotoxicity assays, both in vivo and in vitro, challenging the current default practice for risk assessment of genotoxic chemicals, which assumes a linear dose–response below the lowest tested dose. Statistical methods cannot determine whether a biological threshold exists with sufficient confidence to overturn this assumption of linearity. Indeed, to truly define the shape of the dose–response curves, we must look to the underlying biology and develop targeted experiments to identify and measure the key processes governing the response of the cell to DNA damage. This chapter describes a series of studies aimed at defining the key transcriptional responses. Two approaches were taken to evaluate transcriptional responses preventing micronucleus induction: (1) comparison of gene signatures for several prototype compounds at a single chemical dose that led to a similar activation of the p53-DNA damage pathway (i.e. 1.5-fold increase in total p53); and (2) evaluation of a subset of chemicals with in-depth dose–response studies. The goal of these efforts was to determine the transcriptional pathways responsible for maintaining homeostasis at low levels of DNA damage, i.e., the biological underpinning of threshold-shaped dose–response curves for mutagenicity.

Adaptive Posttranslational Control in Cellular Stress Response Pathways and Its Relationship to Toxicity Testing and Safety Assessment

Although transcriptional induction of stress genes constitutes a major cellular defense program against a variety of stressors, posttranslational control directly regulating the activities of preexisting stress proteins provides a faster-acting alternative response. We propose that posttranslational control is a general adaptive mechanism operating in many stress pathways. Here with the aid of computational models, we first show that posttranslational control fulfills two roles: (1) handling small, transient stresses quickly and (2) stabilizing the negative feedback transcriptional network. We then review the posttranslational control pathways for major stress responses-oxidative stress, metal stress, hyperosmotic stress, DNA damage, heat shock, and hypoxia. Posttranslational regulation of stress protein activities occurs by reversible covalent modifications, allosteric or non-allosteric enzymatic regulations, and physically induced protein structural changes. Acting in feedback or feedforward networks, posttranslational control may establish a threshold level of cellular stress. Sub-threshold stresses are handled adequately by posttranslational control without invoking gene transcription. With supra-threshold stress levels, cellular homeostasis cannot be maintained and transcriptional induction of stress genes and other gene programs, eg, those regulating cell metabolism, proliferation, and apoptosis, takes place. The loss of homeostasis with consequent changes in cellular function may lead to adverse cellular outcomes. Overall, posttranslational and transcriptional control pathways constitute a stratified cellular defense system, handling stresses coherently across time and intensity. As cell-based assays become a focus for chemical testing anchored on toxicity pathways, examination of proteomic and metabolomic changes as a result of posttranslational control occurring in the absence of transcriptomic alterations deserves more attention.

Approaches for characterizing threshold dose-response relationships for DNA-damage pathways involved in carcinogenicity in vivo and micronuclei formation in vitro

Assessing the shape of dose-response curves for DNA-damage in cellular systems and for the consequences of DNA damage in intact animals remains a controversial topic. This overview looks at aspects of the pharmacokinetics (PK) and pharmacodynamics (PD) of cellular DNA-damage/repair and their role in defining the shape of dose-response curves using an in vivo example with formaldehyde and in vitro examples for micronuclei (MN) formation with several test compounds. Formaldehyde is both strongly mutagenic and an endogenous metabolite in cells. With increasing inhaled concentrations, there were transitions in gene changes, from activation of selective stress pathway genes at low concentrations, to activation of pathways for cell-cycle control, p53-DNA damage, and stem cell niche pathways at higher exposures. These gene expression changes were more consistent with dose-dependent transitions in the PD responses to formaldehyde in epithelial cells in the intact rat rather than the low-dose linear extrapolation methods currently used for carcinogens. However, more complete PD explanations of non-linear dose response for creation of fixed damage in cells require detailed examination of cellular responses in vitro using measures of DNA damage and repair that are not easily accessible in the intact animal. In the second section of the article, we illustrate an approach from our laboratory that develops fit-for-purpose, in vitro assays and evaluates the PD of DNA damage and repair through studies using prototypical DNA-damaging agents. Examination of a broad range of responses in these cells showed that transcriptional upregulation of cell cycle control and DNA repair pathways only occurred at doses higher than those causing overt damage fixed damage-measured as MN formation. Lower levels of damage appear to be handled by post-translational repair process using pre-existing proteins. In depth evaluation of the PD properties of one such post-translational process (formation of DNA repair centers; DRCs) has indicated that the formation of DRCs and their ability to complete repair before replication are consistent with threshold behaviours for mutagenesis and, by extension, with chemical carcinogenesis.

Contributions of DNA repair and damage response pathways to the non-linear genotoxic responses of alkylating agents

From a risk assessment perspective, DNA-reactive agents are conventionally assumed to have genotoxic risks at all exposure levels, thus applying a linear extrapolation for low-dose responses. New approaches discussed here, including more diverse and sensitive methods for assessing DNA damage and DNA repair, strongly support the existence of measurable regions where genotoxic responses with increasing doses are insignificant relative to control. Model monofunctional alkylating agents have in vitro and in vivo datasets amenable to determination of points of departure (PoDs) for genotoxic effects. A session at the 2013 Society of Toxicology meeting provided an opportunity to survey the progress in understanding the biological basis of empirically-observed PoDs for DNA alkylating agents. Together with the literature published since, this review discusses cellular pathways activated by endogenous and exogenous alkylation DNA damage. Cells have evolved conserved processes that monitor and counteract a spontaneous steady-state level of DNA damage. The ubiquitous network of DNA repair pathways serves as the first line of defense for clearing of the DNA damage and preventing mutation. Other biological pathways discussed here that are activated by genotoxic stress include post-translational activation of cell cycle networks and transcriptional networks for apoptosis/cell death. The interactions of various DNA repair and DNA damage response pathways provide biological bases for the observed PoD behaviors seen with genotoxic compounds. Thus, after formation of DNA adducts, the activation of cellular pathways can lead to the avoidance of a mutagenic outcome. The understanding of the cellular mechanisms acting within the low-dose region will serve to better characterize risks from exposures to DNA-reactive agents at environmentally-relevant concentrations.

Profiling dose-dependent activation of p53-mediated signaling pathways by chemicals with distinct mechanisms of DNA damage

As part of a larger effort to provide proof-of-concept in vitro-only risk assessments, we have developed a suite of high-throughput assays for key readouts in the p53 DNA damage response toxicity pathway: double-strand break DNA damage (p-H2AX), permanent chromosomal damage (micronuclei), p53 activation, p53 transcriptional activity, and cell fate (cell cycle arrest, apoptosis, micronuclei). Dose-response studies were performed with these protein and cell fate assays, together with whole genome transcriptomics, for three prototype chemicals: etoposide, quercetin, and methyl methanesulfonate. Data were collected in a human cell line expressing wild-type p53 (HT1080) and results were confirmed in a second p53 competent cell line (HCT 116). At chemical concentrations causing similar increases in p53 protein expression, p53-mediated protein expression and cellular processes showed substantial chemical-specific differences. These chemical-specific differences in the p53 transcriptional response appear to be determined by augmentation of the p53 response by co-regulators. More importantly, dose-response data for each of the chemicals indicate that the p53 transcriptional response does not prevent micronuclei induction at low concentrations. In fact, the no observed effect levels and benchmark doses for micronuclei induction were less than or equal to those for p53-mediated gene transcription regardless of the test chemical, indicating that p53's post-translational responses may be more important than transcriptional activation in the response to low dose DNA damage. This effort demonstrates the process of defining key assays required for a pathway-based, in vitro-only risk assessment, using the p53-mediated DNA damage response pathway as a prototype.

Use of mode of action data to inform a dose-response assessment for bladder cancer following exposure to inorganic arsenic

In the recent National Research Council report on conducting a dose-response assessment for inorganic arsenic, the committee remarked that mode of action data should be used, to the extent possible, to extrapolate below the observed range for epidemiological studies to inform the shape of the dose-response curve. Recent in vitro mode of action studies focused on understanding the development of bladder cancer following exposure to inorganic arsenic provide data to inform the dose-response curve. These in vitro data, combined with results of bladder cancer epidemiology studies, inform the dose-response curve in the low-dose region, and include values for both pharmacokinetic and pharmacodynamic variability. Integration of these data provides evidence of a range of concentrations of arsenic for which no effect on the bladder would be expected. Specifically, integration of these results suggest that arsenic exposures in the range of 7-43 ppb in drinking water are exceedingly unlikely to elicit changes leading to key events in the development of cancer or noncancer effects in bladder tissue. These findings are consistent with the lack of evidence for bladder cancer following chronic ingestion of arsenic water concentrations <100 ppb in epidemiological studies.

Physiologically based pharmacokinetic/toxicokinetic modeling

Physiologically based pharmacokinetic (PBPK) models differ from conventional compartmental pharmacokinetic models in that they are based to a large extent on the actual physiology of the organism. The application of pharmacokinetics to toxicology or risk assessment requires that the toxic effects in a particular tissue are related in some way to the concentration time course of an active form of the substance in that tissue. The motivation for applying pharmacokinetics is the expectation that the observed effects of a chemical will be more simply and directly related to a measure of target tissue exposure than to a measure of administered dose. The goal of this work is to provide the reader with an understanding of PBPK modeling and its utility as well as the procedures used in the development and implementation of a model to chemical safety assessment using the styrene PBPK model as an example.

Perchlorate and radioiodide kinetics across life stages in the human: using PBPK models to predict dosimetry and thyroid inhibition and sensitive subpopulations based on developmental stage

Perchlorate (ClO4(-)) is a drinking-water contaminant, known to disrupt thyroid hormone homeostasis in rats. This effect has only been seen in humans at high doses, yet the potential for long term effects from developmental endocrine disruption emphasizes the need for improved understanding of perchlorate's effect during the perinatal period. Physiologically based pharmacokinetic/dynamic (PBPK/PD) models for ClO4(-) and its effect on thyroid iodide uptake were constructed for human gestation and lactation data. Chemical specific parameters were estimated from life-stage and species-specific relationships established in previously published models for various life-stages in the rat and nonpregnant adult human. With the appropriate physiological descriptions, these kinetic models successfully simulate radioiodide data culled from the literature for gestation and lactation, as well as ClO4(-) data from populations exposed to contaminated drinking water. These models provide a framework for extrapolating from chemical exposure in laboratory animals to human response, and support a more quantitative understanding of life-stage-specific susceptibility to ClO4(-). The pregnant and lactating woman, fetus, and nursing infant were predicted to have higher blood ClO4(-) concentrations and greater thyroid iodide uptake inhibition at a given drinking-water concentration than either the nonpregnant adult or the older child. The fetus is predicted to receive the greatest dose (per kilogram body weight) due to several factors, including placental sodium-iodide symporter (NIS) activity and reduced maternal urinary clearance of ClO4(-). The predicted extent of iodide inhibition in the most sensitive population (fetus) is not significant (approximately 1%) at the U.S. Environmental Protection Agency reference dose (0.0007 mg/kg-d).

PBPK model for radioactive iodide and perchlorate kinetics and perchlorate-induced inhibition of iodide uptake in humans

Detection of perchlorate (ClO4-) in several drinking water sources across the U.S. has lead to public concern over health effects from chronic low-level exposures. Perchlorate inhibits thyroid iodide (I-) uptake at the sodium (Na+)-iodide (I-) symporter (NIS), thereby disrupting the initial stage of thyroid hormone synthesis. A physiologically based pharmacokinetic (PBPK) model was developed to describe the kinetics and distribution of both radioactive I- and cold ClO4- in healthy adult humans and simulates the subsequent inhibition of thyroid uptake of radioactive I- by ClO4-. The model successfully predicts the measured levels of serum and urinary ClO4- from drinking water exposures, ranging from 0.007 to 12 mg ClO4-/kg/day, as well as the subsequent inhibition of thyroid 131I- uptake. Thyroid iodine, as well as total, free, and protein-bound radioactive I- in serum from various tracer studies, are also successfully simulated. This model's parameters, in conjunction with corresponding model parameters established for the male, gestational, and lactating rat, can be used to estimate parameters in a pregnant or lactating human, that have not been or cannot be easily measured to extrapolate dose metrics and correlate observed effects in perchlorate toxicity studies to other human life stages. For example, by applying the adult male rat:adult human ratios of model parameters to those parameters established for the gestational and lactating rat, we can derive a reasonable estimate of corresponding parameters for a gestating or lactating human female. Although thyroid hormones and their regulatory feedback are not incorporated in the model structure, the model's successful prediction of free and bound radioactive I- and perchlorate's interaction with free radioactive I- provide a basis for extending the structure to address the complex hypothalamic-pituitary-thyroid feedback system. In this paper, bound radioactive I- refers to I- incorporated into thyroid hormones or iodinated proteins, which may or may not be bound to plasma proteins.

Evidence for competitive inhibition of iodide uptake by perchlorate and translocation of perchlorate into the thyroid

Various published data sets that investigate the potential effect of exogenous perchlorate (ClO4-) on the uptake of iodide in the thyroid and subsequent changes in thyroid hormone levels are available. In order to best use the data towards the prediction of human health effects resulting from ClO4- exposure, the available literature data must be integrated into a self-consistent, coherent, and parsimonious quantitative model based on the most likely mode of action of perchlorate effect on thyroid function. We submit that the simplest mode of action for ClO4- in the thyroid that remains consistent with all available data involves competitive inhibition of iodide transport into the thyroid follicle, transport of perchlorate into the thyroid follicle against a concentration gradient, further transport into the thyroid lumen (where it may again interfere with iodide transport), and, finally, passive diffusion back into the blood. We believe this description of perchlorate's kinetic behavior should serve as the foundation for predictive physiologically based pharmacokinetic (PBPK) models and as a working hypothesis for further experimental exploration.

Predicting neonatal perchlorate dose and inhibition of iodide uptake in the rat during lactation using physiologically-based pharmacokinetic modeling

Perchlorate (ClO4-), a contaminant in drinking water, competitively inhibits active uptake of iodide (I-) into various tissues, including mammary tissue. During postnatal development, inhibition of I- uptake in the mammary gland and neonatal thyroid and the active concentration ClO4- in milk indicate a potentially increased susceptibility of neonates to endocrine disruption. A physiologically based pharmacokinetic (PBPK) model was developed to reproduce measured ClO4- distribution in the lactating and neonatal rat and predict resulting effects on I- kinetics from competitive inhibition at the sodium iodide symporter (NIS). Kinetic I- and ClO4- behavior in tissues with NIS (thyroid, stomach, mammary gland, and skin) was simulated with multiple subcompartments, Michaelis-Menten (M-M) kinetics and competitive inhibition. Physiological and kinetic parameters were obtained from literature and experiment. Systemic clearance and M-M parameters were estimated by fitting simulations to tissue and serum data. The model successfully describes maternal and neonatal thyroid, stomach, skin, and plasma, as well as maternal mammary gland and milk data after ClO4- exposure (from 0.01 to 10 mg/kg-day ClO4-) and acute radioiodide (2.1 to 33,000 ng/kg I-) dosing. The model also predicts I- uptake inhibition in the maternal thyroid, mammary gland, and milk. Model simulations predict a significant transfer of ClO4- through milk after maternal exposure; approximately 50% to 6% of the daily maternal dose at doses ranging from 0.01 to 10.0 mg ClO4-/kg-day, respectively. Comparison of predicted dosimetrics across life-stages in the rat indicates that neonatal thyroid I- uptake inhibition is similar to the adult and approximately tenfold less than the fetus.

PBPK predictions of perchlorate distribution and its effect on thyroid uptake of radioiodide in the male rat

Due to perchlorate's (ClO4-) ability to competitively inhibit thyroid iodide (I-) uptake through the sodium-iodide symporter (NIS), potential human health risks exist from chronic exposure via drinking water. Such risks may include hypothyroidism, goiter, and mental retardation (if exposure occurs during critical periods in neurodevelopment). To aid in predicting perchlorate's effect on normal I- kinetics, we developed a physiologically-based pharmacokinetic (PBPK) model for the adult male rat. The model structure describes simultaneous kinetics for both anions together with their interaction at the NIS, in particular, the inhibition of I- uptake by ClO4-. Subcompartments and Michaelis-Menten (M-M) kinetics were used to describe active uptake of both anions in the thyroid, stomach, and skin. Separate compartments for kidney, liver, plasma, and fat were described by passive diffusion. The model successfully predicts both 36ClO4- and 125I- kinetics after iv doses of 3.3 mg/kg and 33 mg/kg, respectively, as well as inhibition of thyroid 125I- uptake by ClO4- after iv doses of ClO4- (0.01 to 3.0 mg/kg). The model also predicts serum and thyroid ClO4- concentrations from 14-day drinking water exposures (0.01 to 30.0 mg ClO4-/kg/day) and compensation of perchlorate-induced inhibition of radioiodide uptake due to upregulation of the thyroid. The model can be used to extrapolate dose metrics and correlate observed effects in perchlorate toxicity studies to other species and life stages, such as rat gestation (Clewell et al., 2003). Because the model successfully predicts perchlorate's interaction with iodide, it provides a sound basis for future incorporation of the complex hypothalamic-pituitary-thyroid feedback system.

Predicting fetal perchlorate dose and inhibition of iodide kinetics during gestation: a physiologically-based pharmacokinetic analysis of perchlorate and iodide kinetics in the rat

Perchlorate (ClO4-) disrupts endocrine homeostasis by competitively inhibiting the transport of iodide (I-) into the thyroid. The potential for health effects from human exposure to ClO4- in drinking water is not known, but experimental animal studies are suggestive of developmental effects from ClO4- induced iodide deficiency during gestation. Normal hormone-dependent development relies, in part, on synthesis of hormones in the fetal thyroid from maternally supplied iodide. Although ClO4- crosses the placenta, the extent of inhibition in the fetal thyroid is unknown. A physiologically-based pharmacokinetic (PBPK) model was developed to simulate ClO4- exposure and the resulting effect on iodide kinetics in rat gestation. Similar to concurrent model development for the adult male rat, this model includes compartments for thyroid, stomach, skin, kidney, liver, and plasma in both mother and fetus, with additional compartments for the maternal mammary gland, fat, and placenta. Tissues with active uptake are described with multiple compartments and Michaelis-Menten (M-M) kinetics. Physiological and kinetic parameters were obtained from literature and experiment. Systemic clearance, placental-fetal transport, and M-M uptake parameters were estimated by fitting model simulations to experimental data. The PBPK model is able to reproduce maternal and fetal iodide data over five orders of magnitude (0.36 to 33,000 ng/kg 131I-), ClO4- distribution over three orders of magnitude (0.01 to 10 mg/kg-day ClO4-) and inhibition of maternal thyroid and total fetal I- uptake. The model suggests a significant fetal ClO4- dose in late gestation (up to 82% of maternal dose). A comparison of model-predicted internal dosimetrics in the adult male, pregnant, and fetal rat indicates that the fetal thyroid is more sensitive to inhibition than that of the adult.

Pharmacokinetics of toxic chemicals in breast milk: use of PBPK models to predict infant exposure

Factors controlling the transfer of potentially toxic chemicals in the breast milk of nursing mothers include both chemical characteristics, such as lipophilicity, and physiologic changes during lactation. Physiologically based pharmacokinetic (PBPK) models can aid in the prediction of infant exposure via breast milk. Benefits of these quantitative models include the ability to account for changing maternal physiology and transfer kinetics, as well as the chemical-specific characteristics, in order to produce more accurate estimates of neonatal risk. A recently developed PBPK model for perchlorate and iodide kinetics in the lactating and neonatal rat demonstrates the utility of PBPK modeling in predicting maternal and neonatal distribution of these two compounds. This model incorporates time-dependent changes in physiologic characteristics and includes interactions between iodide and perchlorate that alter the distribution and kinetics of iodide.

The use of physiologically based models to integrate diverse data sets and reduce uncertainty in the prediction of perchlorate and iodide kinetics across life stages and species

The effects of perchlorate on the incorporation of iodide into thyroid hormones have been studied for more than 40 years in many species and under varying exposure conditions. Nevertheless, the database for this drinking water contaminant is still incomplete, particularly with regard to human developmental risk. A method for integrating the available data and forming meaningful conclusions for risk assessment is needed. To this end, an initial suite of physiologically based pharmacokinetic (PBPK) models has been developed, which incorporates physiological data for the relevant species and life stages and kinetic data for perchlorate and iodide, as well as the interaction between the two anions. The validated models successfully describe perchlorate-induced inhibition of thyroid iodide uptake and perchlorate and iodide kinetics in the male, pregnant, lactating, fetal, and neonatal rats and the adult humans. The relationships of model-predicted internal dose metrics and kinetic parameters allow a direct comparison of internal dose metrics across life stages in rats and humans. By incorporating all the available data, these models provide a framework for species and life stage extrapolation where the lack of specific data sets would otherwise limit predictive capability. This paper demonstrates two approaches for calculating life stage-specific equivalent doses in a risk assessment for perchlorate: the direct combination of validated model predictions, and the development of preliminary PBPK models for the human-sensitive populations based on the relationship of the parameters in the validated rat and human models. Either approach can be used to perform the needed dosimetry. However, the second approach provides the advantage of a preliminary human life stage-specific PBPK model that can be used for identification of key data gaps and estimation of uncertainty.