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.

A systematic approach to evaluate plausible modes of actions for mouse lung tumors in mice exposed to 4-methylimidozole.

The National Toxicology Program (NTP) reported that chronic dietary exposure to 4-methylimidazole (4-MeI) increased the incidence of lung adenomas/carcinomas beyond the normally high spontaneous rate in B6C3F1 mice. To examine plausible modes of action (MoAs) for mouse lung tumors (MLTs) upon exposure to high levels of 4-MeI, and their relevance in assessing human risk, a systematic approach was used to identify and evaluate mechanistic data (in vitro and in vivo) in the primary and secondary literature, along with high-throughput screening assay data. Study quality, relevance, and activity of mechanistic data identified across the evidence-base were organized according to key characteristics of carcinogens (KCCs) to identify potential key events in known or novel MLT MoAs. Integration of these evidence streams provided confirmation that 4-MeI lacks genotoxic and cytotoxic activity with some evidence to support a lack of mitogenic activity. Further evaluation of contextual and chemical-specific characteristics of 4-MeI was consequently undertaken. Due to lack of genotoxicity, along with transcriptomic and histopathological lung changes up to 28 and 90 days of exposure, the collective evidence suggests MLTs observed following exposure to high levels of 4-MeI develop at a late stage in the mouse chronic bioassay, albeit the exact MoA remains unclear.

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 AC50s 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 toxicogenomic approach for the risk assessment of the food contaminant acetamide.

Acetamide (CAS 60-35-5) is detected in common foods. Chronic rodent bioassays led to its classification as a group 2B possible human carcinogen due to the induction of liver tumors in rats. We used a toxicogenomics approach in Wistar rats gavaged daily for 7 or 28 days at doses of 300 to 1500 mg/kg/day (mkd) to determine a point of departure (POD) and investigate its mode of action (MoA). Ki67 labeling was increased at doses ≥750 mkd up to 3.3-fold representing the most sensitive apical endpoint. Differential gene expression analysis by RNA-Seq identified 1110 and 1814 differentially expressed genes in male and female rats, respectively, following 28 days of treatment. Down-regulated genes were associated with lipid metabolism while up-regulated genes included cell signaling, immune response, and cell cycle functions. Benchmark dose (BMD) modeling of the Ki67 labeling index determined the BMD10 lower confidence limit (BMDL10) as 190 mkd. Transcriptional BMD modeling revealed excellent concordance between transcriptional POD and apical endpoints. Collectively, these results indicate that acetamide is most likely acting through a mitogenic MoA, though specific key initiating molecular events could not be elucidated. A POD value of 190 mkd determined for cell proliferation is suggested for risk assessment purposes.

An in vitro approach for prioritization and evaluation of chemical effects on glucocorticoid receptor mediated adipogenesis.

Rising obesity rates worldwide have socio-economic ramifications. While genetics, diet, and lack of exercise are major contributors to obesity, environmental factors may enhance susceptibility through disruption of hormone homeostasis and metabolic processes. The obesogen hypothesis contends that chemical exposure early in development may enhance adipocyte differentiation, thereby increasing the number of adipocytes and predisposing for obesity and metabolic disease. We previously developed a primary human adipose stem cell (hASC) assay to evaluate the effect of environmental chemicals on PPARG-dependent adipogenesis. Here, the assay was modified to determine the effects of chemicals on the glucocorticoid receptor (GR) pathway. In differentiation cocktail lacking the glucocorticoid agonist dexamethasone (DEX), hASCs do not differentiate into adipocytes. In the presence of GR agonists, adipocyte maturation was observed using phenotypic makers for lipid accumulation, adipokine secretion, and expression of key genes. To evaluate the role of environmental compounds on adipocyte differentiation, progenitor cells were treated with 19 prioritized compounds previously identified by ToxPi as having GR-dependent bioactivity, and multiplexed assays were used to confirm a GR-dependent mode of action. Five chemicals were found to be strong agonists. The assay was also modified to evaluate GR-antagonists, and 8/10 of the hypothesized antagonists inhibited adipogenesis. The in vitro bioactivity data was put into context with extrapolated human steady state concentrations (Css) and clinical exposure data (Cmax). These data support using a human adipose-derived stem cell differentiation assay to test the potential of chemicals to alter human GR-dependent adipogenesis.

Strain-related differences in mouse lung gene expression over a two-year period of inhalation exposure to styrene: Relevance to human risk assessment.

Both CD-1 and C57BL/6 wildtype (C57BL/6-WT) mice show equivalent short-term lung toxicity from exposures to styrene, while long-term tumor responses are greater in CD-1 mice. We analyzed lung gene expression from styrene exposures lasting from 1-day to 2-years in male mice from these two strains, including a Cyp2f2(-/-) knockout (C57BL/6-KO) and a Cyp2F1/2A13/2B6 transgenic mouse (C57BL/6-TG). With short term exposures (1-day to 1-week), CD-1 and C57BL/6-WT mice had thousands of differentially expressed genes (DEGs), consistent with changes in pathways for cell proliferation, cellular lipid metabolism, DNA-replication and inflammation. C57BL/6-WT mice responded within a single day; CD-1 mice required several days of exposure. The numbers of exposure related DEGs were greatly reduced at longer times (4-weeks to 2-years) with enrichment only for biological oxidations in C57BL/6-WT and metabolism of lipids and lipoproteins in CD-1. Gene expression results indicate a non-genotoxic, mouse specific mode of action for short-term styrene responses related to activation of nuclear receptor signaling and cell proliferation. Greater tumor susceptibility in CD-1 mice correlated with the presence of the Pas1 loci, differential Cytochrome P450 gene expression, down-regulation of Nr4a, and greater inflammatory pathway activation. Very few exposure-related responses occurred at any time in C57BL/6-KO or -TG mice indicating that neither the short term nor long term responses of styrene in mice are relevant endpoints for assessing human risks.

Assessing molecular initiating events (MIEs), key events (KEs) and modulating factors (MFs) for styrene responses in mouse lungs using whole genome gene expression profiling following 1-day and multi-week exposures.

Styrene increased lung tumors in mice at chronic inhalation exposures of 20ppm and greater. MIEs, KEs and MFs were examined using gene expression in three strains of male mice (the parental C57BL/6 strain, a CYP2F2(-/-) knock out and a CYP2F2(-/-) transgenic containing human CYP2F1, 2A13 and 2B6). Exposures were for 1-day and 1, 4 and 26weeks. After 1-day exposures at 1, 5, 10, 20, 40 and 120ppm significant increases in differentially expressed genes (DEGs) occurred only in parental strain lungs where there was already an increase in DEGs at 5ppm and then many thousands of DEGs by 120ppm. Enrichment for 1-day and 1-week exposures included cell cycle, mitotic M-M/G1 phases, DNA-synthesis and metabolism of lipids and lipoproteins pathways. The numbers of DEGs decreased steadily over time with no DEGs meeting both statistical significance and fold-change criteria at 26weeks. At 4 and 26weeks, some key transcription factors (TFs) - Nr1d1, Nr1d2, Dbp, Tef, Hlf, Per3, Per2 and Bhlhe40 - were upregulated (|FC|>1.5), while others - Npas, Arntl, Nfil3, Nr4a1, Nr4a2, and Nr4a3 - were down-regulated. At all times, consistent changes in gene expression only occurred in the parental strain. Our results support a MIE for styrene of direct mitogenicity from mouse-specific CYP2F2-mediated metabolites activating Nr4a signaling. Longer-term MFs include down-regulation of Nr4a genes and shifts in both circadian clock TFs and other TFs, linking circadian clock to cellular metabolism. We found no gene expression changes indicative of cytotoxicity or activation of p53-mediated DNA-damage pathways.

Combining transcriptomics and PBPK modeling indicates a primary role of hypoxia and altered circadian signaling in dichloromethane carcinogenicity in mouse lung and liver.

Dichloromethane (DCM) is a lung and liver carcinogen in mice at inhalation exposures≥2000ppm. The modes of action (MOA) of these responses have been attributed to formation of genotoxic, reactive metabolite(s). Here, we examined gene expression in lung and liver from female B6C3F1 mice exposed to 0, 100, 500, 2000, 3000 and 4000ppm DCM for 90days. We also simulated dose measures - rates of DCM oxidation to carbon monoxide (CO) in lung and liver and expected blood carboxyhemoglobin (HbCO) time courses with a PBPK model inclusive of both conjugation and oxidation pathways. Expression of large numbers of genes was altered at 100ppm with maximal changes in the numbers occurring by 500 or 2000ppm. Most changes in genes common to the two tissues were related to cellular metabolism and circadian clock. At the lower concentrations, the changes in metabolism-related genes were discordant - up in liver and down in lung. These processes included organelle biogenesis, TCA cycle, and respiratory electron transport. Changes in circadian cycle genes - primarily transcription factors - showed strong concentration-related response at higher concentrations (Arntl, Npas2, and Clock were down-regulated; Cry2, Wee1, Bhlhe40, Per3, Nr1d1, Nr1d2 and Dbp) were up-regulated with similar directionality in both tissues. Overall, persistently elevated HbCO from DCM oxidation appears to cause extended periods of hypoxia, leading to altered circadian coupling to cellular metabolism. The dose response for altered circadian processes correlates with the cancer outcome. We found no evidence of changes in genes indicative of responses to cytotoxic, DNA-reactive metabolites.

Multiple receptors shape the estrogen response pathway and are critical considerations for the future of in vitro-based risk assessment efforts.

Current in life toxicity testing paradigms are being challenged as the future of risk assessment moves towards more comprehensive mode of action/adverse outcome pathway based approaches. In particular, endocrine disruption screening is now a global activity and key initiatives in the United States focus on the use of high throughput in vitro assays to prioritize compounds for further testing of estrogen, androgen or thyroid disruption. Of these pathways, much of the emphasis to date has been on high-throughput methods for estrogenic activity primarily using ligand binding and trans-activation assays. However, as the knowledge regarding estrogen receptor signaling pathways continues to evolve, it is clear that the assumption of a simple one-receptor pathway underlying current in vitro screening assays is out of date. To develop more accurate models for estrogen-initiated pathways useful for quantitative safety assessments, we must design assays that account for the key signaling processes driving cellular dose response based on up-to-date understanding of the biological network. In this review, we summarize the state of the science for the estrogen receptor signaling network, particularly with regard to proliferative effects, and highlight gaps in current high throughput approaches. From the sum of this literature, we propose a model for the estrogen-signaling pathway that should serve as a starting point for development of new in vitro methods fit for the purpose of predicting dose response for estrogenic chemicals in the human.

Information-dependent enrichment analysis reveals time-dependent transcriptional regulation of the estrogen pathway of toxicity.

The twenty-first century vision for toxicology involves a transition away from high-dose animal studies to in vitro and computational models (NRC in Toxicity testing in the 21st century: a vision and a strategy, The National Academies Press, Washington, DC, 2007). This transition requires mapping pathways of toxicity by understanding how in vitro systems respond to chemical perturbation. Uncovering transcription factors/signaling networks responsible for gene expression patterns is essential for defining pathways of toxicity, and ultimately, for determining the chemical modes of action through which a toxicant acts. Traditionally, transcription factor identification is achieved via chromatin immunoprecipitation studies and summarized by calculating which transcription factors are statistically associated with up- and downregulated genes. These lists are commonly determined via statistical or fold-change cutoffs, a procedure that is sensitive to statistical power and may not be as useful for determining transcription factor associations. To move away from an arbitrary statistical or fold-change-based cutoff, we developed, in the context of the Mapping the Human Toxome project, an enrichment paradigm called information-dependent enrichment analysis (IDEA) to guide identification of the transcription factor network. We used a test case of activation in MCF-7 cells by 17ß estradiol (E2). Using this new approach, we established a time course for transcriptional and functional responses to E2. ERα and ERß were associated with short-term transcriptional changes in response to E2. Sustained exposure led to recruitment of additional transcription factors and alteration of cell cycle machinery. TFAP2C and SOX2 were the transcription factors most highly correlated with dose. E2F7, E2F1, and Foxm1, which are involved in cell proliferation, were enriched only at 24 h. IDEA should be useful for identifying candidate pathways of toxicity. IDEA outperforms gene set enrichment analysis (GSEA) and provides similar results to weighted gene correlation network analysis, a platform that helps to identify genes not annotated to pathways.

Pathway Based Toxicology and Fit-for-Purpose Assays.

The field of toxicity testing for non-pharmaceutical chemicals is in flux with multiple initiatives in North America and the EU to move away from animal testing to mode-of-action based in vitro assays. In this arena, there are still obstacles to overcome, such as developing appropriate cellular assays, creating pathway-based dose-response models and refining in vitro-in vivo extrapolation (IVIVE) tools. Overall, it is necessary to provide assurances that these new approaches are adequately protective of human and ecological health. Another major challenge for individual scientists and regulatory agencies is developing a cultural willingness to shed old biases developed around animal tests and become more comfortable with mode-of-action based assays in human cells. At present, most initiatives focus on developing in vitro alternatives and assessing how well these alternative methods reproduce past results related to predicting organism level toxicity in intact animals. The path forward requires looking beyond benchmarking against high dose animal studies. We need to develop targeted cellular assays, new cell biology-based extrapolation models for assessing regions of safety for chemical exposures in human populations, and mode-of-action-based approaches which are constructed on an understanding of human biology. Furthermore, it is essential that assay developers have the flexibility to 'validate' against the most appropriate mode-of-action data rather than against apical endpoints in high dose animal studies. This chapter demonstrates the principles of fit-for-purpose assay development using pathway-targeted case studies. The projects include p53-mdm2-mediated DNA-repair, estrogen receptor-mediated cell proliferation and PPARα receptor-mediated liver responses.

Identification of a tissue-selective heat shock response regulatory network.

The heat shock response (HSR) is essential to survive acute proteotoxic stress and has been studied extensively in unicellular organisms and tissue culture cells, but to a lesser extent in intact metazoan animals. To identify the regulatory pathways that control the HSR in Caenorhabditis elegans, we performed a genome-wide RNAi screen and identified 59 genes corresponding to 7 positive activators required for the HSR and 52 negative regulators whose knockdown leads to constitutive activation of the HSR. These modifiers function in specific steps of gene expression, protein synthesis, protein folding, trafficking, and protein clearance, and comprise the metazoan heat shock regulatory network (HSN). Whereas the positive regulators function in all tissues of C. elegans, nearly all of the negative regulators exhibited tissue-selective effects. Knockdown of the subunits of the proteasome strongly induces HS reporter expression only in the intestine and spermatheca but not in muscle cells, while knockdown of subunits of the TRiC/CCT chaperonin induces HS reporter expression only in muscle cells. Yet, both the proteasome and TRiC/CCT chaperonin are ubiquitously expressed and are required for clearance and folding in all tissues. We propose that the HSN identifies a key subset of the proteostasis machinery that regulates the HSR according to the unique functional requirements of each tissue.

Using gene expression profiling to evaluate cellular responses in mouse lungs exposed to V2O5 and a group of other mouse lung tumorigens and non-tumorigens.

Many compounds test positive for lung tumors in two-year NTP carcinogenicity bioassays in B6C3F1 mice. V2O5 was identified as a lung carcinogen in this assay, leading to its IARC (International Agency for Research on Cancer) classification as group 2b or a "possible" human carcinogen. To assess potential tumorigenic mode of action of V2O5, we compared gene expression and gene ontology enrichment in lung tissue of female B6C3F1 mice exposed for 13 weeks to a V2O5 particulate aerosol at a tumorigenic level (2.0 mg/m(3)). Relative to 12 other compounds also tested for carcinogenicity in 2-year bioassays in mice, there were 1026 differentially expressed genes with V2O5, of which 483 were unique to V2O5. Ontology analysis of the 1026 V2O5 differentially expressed genes showed enrichment for hyaluronan and sphingolipid metabolism, adenylate cyclase functions, c-AMP signaling and PKA activation/signaling. Enrichment of lipids/lipoprotein metabolism and inflammatory pathways were consistent with previously reported clinical findings. Enrichment of c-AMP and PKA signaling pathways may arise due to inhibition of phosphatases, a known biological action of vanadate. We saw no enrichment for DNA-damage, oxidative stress, cell cycle, or apoptosis pathway signaling in mouse lungs exposed to V2O5 which is in contrast with past studies evaluating in vivo gene expression in target tissues of other carcinogens (arsenic, formaldehyde, naphthalene and chloroprene).

Integrating gene expression and splicing dynamics across dose-response oxidative modulators.

Toxicological risk assessment increasingly utilizes transcriptomics to derive point of departure (POD) and modes of action (MOA) for chemicals. One essential biological process that allows a single gene to generate several different RNA isoforms is called alternative splicing. To comprehensively assess the role of splicing dysregulation in toxicological evaluation and elucidate its potential as a complementary endpoint, we performed RNA-seq on A549 cells treated with five oxidative stress modulators across a wide dose range. Differential gene expression (DGE) showed limited pathway enrichment except at high concentrations. However, alternative splicing analysis revealed variable intron retention events affecting diverse pathways for all chemicals in the absence of significant expression changes. For instance, diazinon elicited negligible gene expression changes but progressive increase in the number of intron retention events, suggesting splicing alterations precede expression responses. Benchmark dose modeling of intron retention data highlighted relevant pathways overlooked by expression analysis. Systematic integration of splicing datasets should be a useful addition to the toxicogenomic toolkit. Combining both modalities paint a more complete picture of transcriptomic dose-responses. Overall, evaluating intron retention dynamics afforded by toxicogenomics may provide biomarkers that can enhance chemical risk assessment and regulatory decision making. This work highlights splicing-aware toxicogenomics as a possible additional tool for examining cellular responses.

Macro-level modeling of the response of C. elegans reproduction to chronic heat stress.

A major goal of systems biology is to understand how organism-level behavior arises from a myriad of molecular interactions. Often this involves complex sets of rules describing interactions among a large number of components. As an alternative, we have developed a simple, macro-level model to describe how chronic temperature stress affects reproduction in C. elegans. Our approach uses fundamental engineering principles, together with a limited set of experimentally derived facts, and provides quantitatively accurate predictions of performance under a range of physiologically relevant conditions. We generated detailed time-resolved experimental data to evaluate the ability of our model to describe the dynamics of C. elegans reproduction. We find considerable heterogeneity in responses of individual animals to heat stress, which can be understood as modulation of a few processes and may represent a strategy for coping with the ever-changing environment. Our experimental results and model provide quantitative insight into the breakdown of a robust biological system under stress and suggest, surprisingly, that the behavior of complex biological systems may be determined by a small number of key components.

Physically grounded approach for estimating gene expression from microarray data.

High-throughput technologies, including gene-expression microarrays, hold great promise for the systems-level study of biological processes. Yet, challenges remain in comparing microarray data from different sources and extracting information about low-abundance transcripts. We demonstrate that these difficulties arise from limitations in the modeling of the data. We propose a physically motivated approach for estimating gene-expression levels from microarray data, an approach neglected in the microarray literature. We separately model the noises specific to sample amplification, hybridization, and fluorescence detection, combining these into a parsimonious description of the variability sources in a microarray experiment. We find that our model produces estimates of gene expression that are reproducible and unbiased. While the details of our model are specific to gene-expression microarrays, we argue that the physically grounded modeling approach we pursue is broadly applicable to other molecular biology technologies.

Beyond pharmaceuticals: Fit-for-purpose new approach methodologies for environmental cardiotoxicity testing.

Environmental factors play a substantial role in determining cardiovascular health, but data informing the risks presented by environmental toxicants is insufficient. In vitro new approach methodologies (NAMs) offer a promising approach with which to address the limitations of traditional in vivo and in vitro assays for assessing cardiotoxicity. Driven largely by the needs of pharmaceutical toxicity testing, considerable progress in developing NAMs for cardiotoxicity analysis has already been made. As the scientific and regulatory interest in NAMs for environmental chemicals continues to grow, a thorough understanding of the unique features of environmental cardiotoxicants and their associated cardiotoxicities is needed. Here, we review the key characteristics of as well as important regulatory and biological considerations for fit-for-purpose NAMs for environmental cardiotoxicity. By emphasizing the challenges and opportunities presented by NAMs for environmental cardiotoxicity we hope to accelerate their development, acceptance, and application.

Biological system considerations for application of toxicogenomics in next-generation risk assessment and predictive toxicology.

There is increasing interest in using modern 'omics technologies, such as whole transcriptome sequencing, to inform decisions about human health safety and chemical toxicity hazard. High throughput methodologies using in vitro assays offer a path forward in reducing or eliminating animal testing. However, many aspects of these technologies need assessment before they will gain the trust of regulators and the public as viable alternative test methods for human health and safety. We used a high throughput whole transcriptome sequence assay (TempO-Seq) to assess the use of three widely used cancer cell lines (HepG2, MCF7, and Ishikawa cells) as in vitro systems for determination of cellular modes of action for two well studied compounds with canonical liver responses: ketoconazole and phenobarbital. We evaluated transcriptomic data to infer points of departure for use in risk analyses of compounds. Both compounds displayed shortcomings in evidence for canonical liver-related responses in any cell line, despite a strong dose response in all three. This raises questions about the competence of simple, mono-cultured cancer cell lines as appropriate surrogates for some adverse effects or toxic endpoints. Points of departure derived from benchmark doses were highly consistent across all three cell lines however, indicating the use of transcriptomic BMD analyses for such purposes would be a reliable and consistent approach.

NAM-based prediction of point-of-contact toxicity in the lung: A case example with 1,3-dichloropropene.

Time, cost, ethical, and regulatory considerations surrounding in vivo testing methods render them insufficient to meet existing and future chemical safety testing demands. There is a need for the development of in vitro and in silico alternatives to replace traditional in vivo methods for inhalation toxicity assessment. Exposures of differentiated airway epithelial cultures to gases or aerosols at the air-liquid interface (ALI) can assess tissue responses and in vitro to in vivo extrapolation can align in vitro exposure levels with in-life exposures expected to give similar tissue exposures. Because the airway epithelium varies along its length, with various regions composed of different cell types, we have introduced a known toxic vapor to five human-derived, differentiated, in vitro airway epithelial cell culture models-MucilAir of nasal, tracheal, or bronchial origin, SmallAir, and EpiAlveolar-representing five regions of the airway epithelium-nasal, tracheal, bronchial, bronchiolar, and alveolar. We have monitored toxicity in these cultures 24 h after acute exposure using an assay for transepithelial conductance (for epithelial barrier integrity) and the lactate dehydrogenase (LDH) release assay (for cytotoxicity). Our vapor of choice in these experiments was 1,3-dichloropropene (1,3-DCP). Finally, we have developed an airway dosimetry model for 1,3-DCP vapor to predict in vivo external exposure scenarios that would produce toxic local tissue concentrations as determined by in vitro experiments. Measured in vitro points of departure (PoDs) for all tested cell culture models were similar. Calculated rat equivalent inhaled concentrations varied by model according to position of the modeled tissue within the airway, with nasal respiratory tissue being the most proximal and most sensitive tissue, and alveolar epithelium being the most distal and least sensitive tissue. These predictions are qualitatively in accordance with empirically determined in vivo PoDs. The predicted PoD concentrations were close to, but slightly higher than, PoDs determined by in vivo subchronic studies.

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 2 = 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.