Pervasive environmental chemicals impair oligodendrocyte development

Exposure to environmental chemicals can impair neurodevelopment, and oligodendrocytes may be particularly vulnerable, as their development extends from gestation into adulthood. However, few environmental chemicals have been assessed for potential risks to oligodendrocytes. Here, using a high-throughput developmental screen in cultured cells, we identified environmental chemicals in two classes that disrupt oligodendrocyte development through distinct mechanisms. Quaternary compounds, ubiquitous in disinfecting agents and personal care products, were potently and selectively cytotoxic to developing oligodendrocytes, whereas organophosphate flame retardants, commonly found in household items such as furniture and electronics, prematurely arrested oligodendrocyte maturation. Chemicals from each class impaired oligodendrocyte development postnatally in mice and in a human 3D organoid model of prenatal cortical development. Analysis of epidemiological data showed that adverse neurodevelopmental outcomes were associated with childhood exposure to the top organophosphate flame retardant identified by our screen. This work identifies toxicological vulnerabilities for oligodendrocyte development and highlights the need for deeper scrutiny of these compounds' impacts on human health.

Adverse effects in traditional and alternative toxicity tests

Chemical safety assessment begins with defining the lowest level of chemical that alters one or more measured endpoints. This critical effect level, along with factors to account for uncertainty, is used to derive limits for human exposure. In the absence of data regarding the specific mechanisms or biological pathways affected, non-specific endpoints such as body weight and non-target organ weight changes are used to set critical effect levels. Specific apical endpoints such as impaired reproductive function or altered neurodevelopment have also been used to set chemical safety limits; however, in test guidelines designed for specific apical effect(s), concurrently measured non-specific endpoints may be equally or more sensitive than specific endpoints. This means that rather than predicting a specific toxicological response, animal data are often used to develop protective critical effect levels, without assuming the same change would be observed in humans. This manuscript is intended to encourage a rethinking of how adverse chemical effects are interpreted: non-specific endpoints from in vivo toxicological studies data are often used to derive points of departure for use with safety assessment factors to create recommended exposure levels that are broadly protective but not necessarily target-specific.

Screening for drinking water contaminants of concern using an automated exposure-focused workflow

BACKGROUND

The number of chemicals present in the environment exceeds the capacity of government bodies to characterize risk. Therefore, data-informed and reproducible processes are needed for identifying chemicals for further assessment. The Minnesota Department of Health (MDH), under its Contaminants of Emerging Concern (CEC) initiative, uses a standardized process to screen potential drinking water contaminants based on toxicity and exposure potential.

OBJECTIVE

Recently, MDH partnered with the U.S. Environmental Protection Agency (EPA) Office of Research and Development (ORD) to accelerate the screening process via development of an automated workflow accessing relevant exposure data, including exposure new approach methodologies (NAMs) from ORD's ExpoCast project.

METHODS

The workflow incorporated information from 27 data sources related to persistence and fate, release potential, water occurrence, and exposure potential, making use of ORD tools for harmonization of chemical names and identifiers. The workflow also incorporated data and criteria specific to Minnesota and MDH's regulatory authority. The collected data were used to score chemicals using quantitative algorithms developed by MDH. The workflow was applied to 1867 case study chemicals, including 82 chemicals that were previously manually evaluated by MDH.

RESULTS

Evaluation of the automated and manual results for these 82 chemicals indicated reasonable agreement between the scores although agreement depended on data availability; automated scores were lower than manual scores for chemicals with fewer available data. Case study chemicals with high exposure scores included disinfection by-products, pharmaceuticals, consumer product chemicals, per- and polyfluoroalkyl substances, pesticides, and metals. Scores were integrated with in vitro bioactivity data to assess the feasibility of using NAMs for further risk prioritization.

SIGNIFICANCE

This workflow will allow MDH to accelerate exposure screening and expand the number of chemicals examined, freeing resources for in-depth assessments. The workflow will be useful in screening large libraries of chemicals for candidates for the CEC program.

Exploring the effects of experimental parameters and data modeling approaches on in vitro transcriptomic point-of-departure estimates

Multiple new approach methods (NAMs) are being developed to rapidly screen large numbers of chemicals to aid in hazard evaluation and risk assessments. High-throughput transcriptomics (HTTr) in human cell lines has been proposed as a first-tier screening approach for determining the types of bioactivity a chemical can cause (activation of specific targets vs. generalized cell stress) and for calculating transcriptional points of departure (tPODs) based on changes in gene expression. In the present study, we examine a range of computational methods to calculate tPODs from HTTr data, using six data sets in which MCF7 cells cultured in two different media formulations were treated with a panel of 44 chemicals for 3 different exposure durations (6, 12, 24 hr). The tPOD calculation methods use data at the level of individual genes and gene set signatures, and compare data processed using the ToxCast Pipeline 2 (tcplfit2), BMDExpress and PLIER (Pathway Level Information ExtractoR). Methods were evaluated by comparing to in vitro PODs from a validated set of high-throughput screening (HTS) assays for a set of estrogenic compounds. Key findings include: (1) for a given chemical and set of experimental conditions, tPODs calculated by different methods can vary by several orders of magnitude; (2) tPODs are at least as sensitive to computational methods as to experimental conditions; (3) in comparison to an external reference set of PODs, some methods give generally higher values, principally PLIER and BMDExpress; and (4) the tPODs from HTTr in this one cell type are mostly higher than the overall PODs from a broad battery of targeted in vitro ToxCast assays, reflecting the need to test chemicals in multiple cell types and readout technologies for in vitro hazard screening.

Reproducibility of organ-level effects in repeat dose animal studies

This work estimates benchmarks for new approach method (NAM) performance in predicting organ-level effects in repeat dose studies of adult animals based on variability in replicate animal studies. Treatment-related effect values from the Toxicity Reference database (v2.1) for weight, gross, or histopathological changes in the adrenal gland, liver, kidney, spleen, stomach, and thyroid were used. Rates of chemical concordance among organ-level findings in replicate studies, defined by repeated chemical only, chemical and species, or chemical and study type, were calculated. Concordance was 39 - 88%, depending on organ, and was highest within species. Variance in treatment-related effect values, including lowest effect level (LEL) values and benchmark dose (BMD) values when available, was calculated by organ. Multilinear regression modeling, using study descriptors of organ-level effect values as covariates, was used to estimate total variance, mean square error (MSE), and root residual mean square error (RMSE). MSE values, interpreted as estimates of unexplained variance, suggest study descriptors accounted for 52-69% of total variance in organ-level LELs. RMSE ranged from 0.41 - 0.68 log10-mg/kg/day. Differences between organ-level effects from chronic (CHR) and subchronic (SUB) dosing regimens were also quantified. Odds ratios indicated CHR organ effects were unlikely if the SUB study was negative. Mean differences of CHR - SUB organ-level LELs ranged from -0.38 to -0.19 log10 mg/kg/day; the magnitudes of these mean differences were less than RMSE for replicate studies. Finally, in vitro to in vivo extrapolation (IVIVE) was employed to compare bioactive concentrations from in vitro NAMs for kidney and liver to LELs. The observed mean difference between LELs and mean IVIVE dose predictions approached 0.5 log10-mg/kg/day, but differences by chemical ranged widely. Overall, variability in repeat dose organ-level effects suggests expectations for quantitative accuracy of NAM prediction of LELs should be at least ± 1 log10-mg/kg/day, with qualitative accuracy not exceeding 70%.

Assessing utility of thyroid in vitro screening assays through comparisons to observed impacts in vivo

To better understand endocrine disruption, the U.S. Environmental Protection Agency's (USEPA) Endocrine Disruptor Screening Program (EDSP) utilizes a two-tiered approach to investigate the potential of a chemical to interact with the estrogen, androgen, or thyroid systems. As in vivo testing lacks the throughput to address data gaps on endocrine bioactivity for thousands of chemicals, in vitro high-throughput screening (HTS) methods are being developed to screen larger chemical libraries. The primary objective of this work was to investigate for how many of the 52 chemicals with weight-of-evidence (WoE) determinations from EDSP Tier 1 screening there are available in vitro HTS data supporting a thyroid impact. HTS data from the USEPA ToxCast program and the EDSP WoE were collected for this analysis. Considering the complexity of endocrine disruption and interpreting HTS data, concordance between in vitro activity and in vivo effects ranges from 58 to 78%. Based on this evaluation, we conclude that the current suite of HTS assays is beneficial for prioritizing chemicals for further inquiry; however, without a more detailed analysis, one cannot conclude whether HTS results are the primary mode-of-action. Furthermore, development of in vitro assays for additional thyroid-relevant molecular initiating events is required to effectively predict in vivo thyroid impacts.

Evaluation of a high-throughput H295R homogenous time resolved fluorescence assay for androgen and estrogen steroidogenesis screening

The H295R test guideline assay evaluates the effect of test substances on synthesis of 17β-estradiol (E2) and testosterone (T). The objective of this study was to leverage commercial immunoassay technology to develop a more efficient H295R assay to measure E2 and T levels in 384-well format. The resulting Homogenous Time Resolved Fluorescence assay platform (H295R-HTRF) was evaluated against a training set of 36 chemicals derived from the OECD inter-laboratory validation study, EPA guideline 890.1200 aromatase assay, and azole fungicides active in the HT-H295R assay. Quality control performance criteria were met for all conditions except E2 synthesis inhibition where low basal hormone synthesis was observed. Five proficiency chemicals were active for both the E2 and T endpoints, consistent with guideline classifications. Of the nine OECD core reference chemicals, 9/9 were concordant with outcomes for E2 and 7/9 for T. Likewise, 9/13 and 11/13 OECD supplemental chemicals were concordant with anticipated effects for E2 and T, respectively. Of the 10 azole fungicides screened, 7/10 for E2 and 8/10 for T exhibited concordant outcomes for inhibition. Generally, all active chemicals in the training set demonstrated equivalent or greater potency in the H295R-HTRF assay, supporting the sensitivity of the platform. The adaptation of HTRF technology to the H295R model provides an efficient way to evaluate E2 and T modulators in accordance with guideline specifications.

ToxRefDB v2.1: update to curated in vivo study data in the Toxicity Reference Database

The Toxicity Reference Database (ToxRefDB) contains in vivo study data from over 5,900 guideline or guideline-like studies for over 1,100 chemicals. The database includes information regarding study design, chemical treatment, dosing, treatment group parameters, treatment-related (significantly different from control) and critical (adverse) effects, guided by a controlled effect vocabulary, as well as endpoint testing status according to health effects guideline requirements. ToxRefDB v2.1 is an update to address a compilation error found in ToxRefDB v2.0 that resulted in some effects being inadvertently omitted from the database. Though effect data has been recovered, no new studies were added. The recovered data improves the utility of ToxRefDB as a resource for curated legacy in vivo information, which enhances scientific confidence in vitro high-throughput screening of chemicals and supports retrospective and predictive toxicology applications for which outcomes in traditional regulatory toxicology studies serve as reference information.

Application of cell painting for chemical hazard evaluation in support of screening-level chemical assessments

'Cell Painting' is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to ~100 μM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.

Evaluation of Per- and Polyfluoroalkyl Substances (PFAS) In Vitro Toxicity Testing for Developmental Neurotoxicity

Per- and polyfluoroalkyl substances (PFAS) are a diverse set of commercial chemicals widely detected in humans and the environment. However, only a limited number of PFAS are associated with epidemiological or experimental data for hazard identification. To provide developmental neurotoxicity (DNT) hazard information, the work herein employed DNT new approach methods (NAMs) to generate in vitro screening data for a set of 160 PFAS. The DNT NAMs battery was comprised of the microelectrode array neuronal network formation assay (NFA) and high-content imaging (HCI) assays to evaluate proliferation, apoptosis, and neurite outgrowth. The majority of PFAS (118/160) were inactive or equivocal in the DNT NAMs, leaving 42 active PFAS that decreased measures of neural network connectivity and neurite length. Analytical quality control indicated 43/118 inactive PFAS samples and 10/42 active PFAS samples were degraded; as such, careful interpretation is required as some negatives may have been due to loss of the parent PFAS, and some actives may have resulted from a mixture of parent and/or degradants of PFAS. PFAS containing a perfluorinated carbon (C) chain length ≥8, a high C:fluorine ratio, or a carboxylic acid moiety were more likely to be bioactive in the DNT NAMs. Of the PFAS positives in DNT NAMs, 85% were also active in other EPA ToxCast assays, whereas 79% of PFAS inactives in the DNT NAMs were active in other assays. These data demonstrate that a subset of PFAS perturb neurodevelopmental processes in vitro and suggest focusing future studies of DNT on PFAS with certain structural feature descriptors.

Environmental mixtures and breast cancer: identifying co-exposure patterns between understudied vs breast cancer-associated chemicals using chemical inventory informatics

BACKGROUND

Although evidence linking environmental chemicals to breast cancer is growing, mixtures-based exposure evaluations are lacking.

OBJECTIVE

This study aimed to identify environmental chemicals in use inventories that co-occur and share properties with chemicals that have association with breast cancer, highlighting exposure combinations that may alter disease risk.

METHODS

The occurrence of chemicals within chemical use categories was characterized using the Chemical and Products Database. Co-exposure patterns were evaluated for chemicals that have an association with breast cancer (BC), no known association (NBC), and understudied chemicals (UC) identified through query of the Silent Spring Institute's Mammary Carcinogens Review Database and the U.S. Environmental Protection Agency's Toxicity Reference Database. UCs were ranked based on structure and physicochemical similarities and co-occurrence patterns with BCs within environmentally relevant exposure sources.

RESULTS

A total of 6793 chemicals had data available for exposure source occurrence analyses. 50 top-ranking UCs spanning five clusters of co-occurring chemicals were prioritized, based on shared properties with co-occuring BCs, including chemicals used in food production and consumer/personal care products, as well as potential endocrine system modulators.

SIGNIFICANCE

Results highlight important co-exposure conditions that are likely prevalent within our everyday environments that warrant further evaluation for possible breast cancer risk.

IMPACT STATEMENT

Most environmental studies on breast cancer have focused on evaluating relationships between individual, well-known chemicals and breast cancer risk. This study set out to expand this research field by identifying understudied chemicals and mixtures that may occur in everyday environments due to their patterns of commercial use. Analyses focused on those that co-occur alongside chemicals associated with breast cancer, based upon in silico chemical database querying and analysis. Particularly in instances when understudied chemicals share physicochemical properties and structural features with carcinogens, these chemical mixtures represent conditions that should be studied in future clinical, epidemiological, and toxicological studies.

Environmental mixtures and breast cancer: identifying co-exposure patterns between understudied vs breast cancer-associated chemicals using chemical inventory informatics

BACKGROUND

Although evidence linking environmental chemicals to breast cancer is growing, mixtures-based exposure evaluations are lacking.

OBJECTIVE

This study aimed to identify environmental chemicals in use inventories that co-occur and share properties with chemicals that have association with breast cancer, highlighting exposure combinations that may alter disease risk.

METHODS

The occurrence of chemicals within chemical use categories was characterized using the Chemical and Products Database. Co-exposure patterns were evaluated for chemicals that have an association with breast cancer (BC), no known association (NBC), and understudied chemicals (UC) identified through query of the Silent Spring Institute's Mammary Carcinogens Review Database and the U.S. Environmental Protection Agency's Toxicity Reference Database. UCs were ranked based on structure and physicochemical similarities and co-occurrence patterns with BCs within environmentally relevant exposure sources.

RESULTS

A total of 6793 chemicals had data available for exposure source occurrence analyses. 50 top-ranking UCs spanning five clusters of co-occurring chemicals were prioritized, based on shared properties with co-occuring BCs, including chemicals used in food production and consumer/personal care products, as well as potential endocrine system modulators.

SIGNIFICANCE

Results highlight important co-exposure conditions that are likely prevalent within our everyday environments that warrant further evaluation for possible breast cancer risk.

IMPACT STATEMENT

Most environmental studies on breast cancer have focused on evaluating relationships between individual, well-known chemicals and breast cancer risk. This study set out to expand this research field by identifying understudied chemicals and mixtures that may occur in everyday environments due to their patterns of commercial use. Analyses focused on those that co-occur alongside chemicals associated with breast cancer, based upon in silico chemical database querying and analysis. Particularly in instances when understudied chemicals share physicochemical properties and structural features with carcinogens, these chemical mixtures represent conditions that should be studied in future clinical, epidemiological, and toxicological studies.

Evaluation of 147 perfluoroalkyl substances for immunotoxic and other (patho)physiological activities through phenotypic screening of human primary cells

A structurally diverse set of 147 per- and polyfluoroalkyl substances (PFAS) was screened in a panel of 12 human primary cell systems by measuring 148 biomarkers relevant to (patho)physiological pathways to inform hypotheses about potential mechanistic effects of data-poor PFAS in human model systems. This analysis focused on immunosuppressive activity, which was previously reported as an in vivo effect of perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), by comparing PFAS responses to four pharmacological immunosuppressants. The PFOS response profile had little correlation with reference immunosuppressants, suggesting in vivo activity does not occur by similar mechanisms. The PFOA response profile did share features with the profile of dexamethasone, although some distinct features were lacking. Other PFAS, including 2,2,3,3-tetrafluoropropyl acrylate, demonstrated more similarity to the reference immunosuppressants but with additional activities not found in the reference immunosuppressive drugs. Correlation of PFAS profiles with a database of environmental chemical responses and pharmacological probes identified potential mechanisms of bioactivity for some PFAS, including responses similar to ubiquitin ligase inhibitors, deubiquitylating enzyme (DUB) inhibitors, and thioredoxin reductase inhibitors. Approximately 21% of the 147 PFAS with confirmed sample quality were bioactive at nominal testing concentrations in the 1-60 micromolar range in these human primary cell systems. These data provide new hypotheses for mechanisms of action for a subset of PFAS and may further aid in development of a PFAS categorization strategy useful in safety assessment.

Advances in computational methods along the exposure to toxicological response paradigm

Human health risk assessment is a function of chemical toxicity, bioavailability to reach target biological tissues, and potential environmental exposure. These factors are complicated by many physiological, biochemical, physical and lifestyle factors. Furthermore, chemical health risk assessment is challenging in view of the large, and continually increasing, number of chemicals found in the environment. These challenges highlight the need to prioritize resources for the efficient and timely assessment of those environmental chemicals that pose greatest health risks. Computational methods, either predictive or investigative, are designed to assist in this prioritization in view of the lack of cost prohibitive in vivo experimental data. Computational methods provide specific and focused toxicity information using in vitro high throughput screening (HTS) assays. Information from the HTS assays can be converted to in vivo estimates of chemical levels in blood or target tissue, which in turn are converted to in vivo dose estimates that can be compared to exposure levels of the screened chemicals. This manuscript provides a review for the landscape of computational methods developed and used at the U.S. Environmental Protection Agency (EPA) highlighting their potentials and challenges.

An expert-driven literature review of "negative" chemicals for developmental neurotoxicity (DNT) in vitro assay evaluation

To date, approximately 200 chemicals have been tested in US Environmental Protection Agency (EPA) or Organization for Economic Co-operation and Development (OECD) developmental neurotoxicity (DNT) guideline studies, leaving thousands of chemicals without traditional animal information on DNT hazard potential. To address this data gap, a battery of in vitro DNT new approach methodologies (NAMs) has been proposed. Evaluation of the performance of this battery will increase the confidence in its use to determine DNT chemical hazards. One approach to evaluate DNT NAM performance is to use a set of chemicals to evaluate sensitivity and specificity. Since a list of chemicals with potential evidence of in vivo DNT has been established, this study aims to develop a curated list of "negative" chemicals for inclusion in a "DNT NAM evaluation set". A workflow, including a literature search followed by an expert-driven literature review, was used to systematically screen 39 chemicals for lack of DNT effect. Expert panel members evaluated the scientific robustness of relevant studies to inform chemical categorizations. Following review, the panel discussed each chemical and made categorical determinations of "Favorable", "Not Favorable", or "Indeterminate" reflecting acceptance, lack of suitability, or uncertainty given specific limitations and considerations, respectively. The panel determined that 10, 22, and 7 chemicals met the criteria for "Favorable", "Not Favorable", and "Indeterminate", for use as negatives in a DNT NAM evaluation set. Ultimately, this approach not only supports DNT NAM performance evaluation but also highlights challenges in identifying large numbers of negative DNT chemicals.

Integration of toxicodynamic and toxicokinetic new approach methods into a weight-of-evidence analysis for pesticide developmental neurotoxicity assessment: A case-study with DL- and L-glufosinate

DL-glufosinate ammonium (DL-GLF) is a registered herbicide for which a guideline Developmental Neurotoxicity (DNT) study has been conducted. Offspring effects included altered brain morphometrics, decreased body weight, and increased motor activity. Guideline DNT studies are not available for its enriched isomers L-GLF acid and L-GLF ammonium; conducting one would be time consuming, resource-intensive, and possibly redundant given the existing DL-GLF DNT. To support deciding whether to request a guideline DNT study for the L-GLF isomers, DL-GLF and the L-GLF isomers were screened using in vitro assays for network formation and neurite outgrowth. DL-GLF and L-GLF isomers were without effects in both assays. DL-GLF and L-GLF (1-100 μM) isomers increased mean firing rate of mature networks to 120-140% of baseline. In vitro toxicokinetic assessments were used to derive administered equivalent doses (AEDs) for the in vitro testing concentrations. The AED for L-GLF was ∼3X higher than the NOAEL from the DL-GLF DNT indicating that the available guideline study would be protective of potential DNT due to L-GLF exposure. Based in part on the results of these in vitro studies, EPA is not requiring L-GLF isomer guideline DNT studies, thereby providing a case study for a useful application of DNT screening assays.

tcplfit2: an R-language general purpose concentration-response modeling package

SUMMARY

Many applications of chemical screening are performed in concentration or dose-response mode, and it is necessary to extract appropriate parameters, including whether the chemical/assay pair is active and if so, what are concentrations where activity is seen. Typically, multiple mathematical models or curve shapes are tested against the data to assess the best fit. There are several commercial programs used for this purpose as well as open-source libraries. A widely used system for managing high-throughput screening (HTS) concentration-response data is tcpl (ToxCast Pipeline). The current implementation of tcpl has the concentration-response modeling code tightly integrated with the data management and databasing aspects of HTS data processing. Tcplfit2 is a stand-alone version of the curve-fitting and hitcalling core of tcpl that has been extended to include a large number of standard curve classes and to use benchmark dose modeling. This package will be useful for HTS concentration-response data such as high-throughput whole genome transcriptomics.

AVAILABILITY AND IMPLEMENTATION

tcplfit2 is written in R and is available from CRAN.

OUP accepted manuscript

In vivo developmental neurotoxicity (DNT) testing is resource intensive and lacks information on cellular processes affected by chemicals. To address this, DNT new approach methodologies (NAMs) are being evaluated, including: the microelectrode array neuronal network formation assay; and high-content imaging to evaluate proliferation, apoptosis, neurite outgrowth, and synaptogenesis. This work addresses 3 hypotheses: (1) a broad screening battery provides a sensitive marker of DNT bioactivity; (2) selective bioactivity (occurring at noncytotoxic concentrations) may indicate functional processes disrupted; and, (3) a subset of endpoints may optimally classify chemicals with in vivo evidence for DNT. The dataset was comprised of 92 chemicals screened in all 57 assay endpoints sourced from publicly available data, including a set of DNT NAM evaluation chemicals with putative positives (53) and negatives (13). The DNT NAM battery provides a sensitive marker of DNT bioactivity, particularly in cytotoxicity and network connectivity parameters. Hierarchical clustering suggested potency (including cytotoxicity) was important for classifying positive chemicals with high sensitivity (93%) but failed to distinguish patterns of disrupted functional processes. In contrast, clustering of selective values revealed informative patterns of differential activity but demonstrated lower sensitivity (74%). The false negatives were associated with several limitations, such as the maximal concentration tested or gaps in the biology captured by the current battery. This work demonstrates that this multi-dimensional assay suite provides a sensitive biomarker for DNT bioactivity, with selective activity providing possible insight into specific functional processes affected by chemical exposure and a basis for further research.

Utility of In Vitro Bioactivity as a Lower Bound Estimate of In Vivo Adverse Effect Levels and in Risk-Based Prioritization

Use of high-throughput, in vitro bioactivity data in setting a point-of-departure (POD) has the potential to accelerate the pace of human health safety evaluation by informing screening-level assessments. The primary objective of this work was to compare PODs based on high-throughput predictions of bioactivity, exposure predictions, and traditional hazard information for 448 chemicals. PODs derived from new approach methodologies (NAMs) were obtained for this comparison using the 50th (PODNAM, 50) and the 95th (PODNAM, 95) percentile credible interval estimates for the steady-state plasma concentration used in in vitro to in vivo extrapolation of administered equivalent doses. Of the 448 substances, 89% had a PODNAM, 95 that was less than the traditional POD (PODtraditional) value. For the 48 substances for which PODtraditional < PODNAM, 95, the PODNAM and PODtraditional were typically within a factor of 10 of each other, and there was an enrichment of chemical structural features associated with organophosphate and carbamate insecticides. When PODtraditional < PODNAM, 95, it did not appear to result from an enrichment of PODtraditional based on a particular study type (eg, developmental, reproductive, and chronic studies). Bioactivity:exposure ratios, useful for identification of substances with potential priority, demonstrated that high-throughput exposure predictions were greater than the PODNAM, 95 for 11 substances. When compared with threshold of toxicological concern (TTC) values, the PODNAM, 95 was greater than the corresponding TTC value 90% of the time. This work demonstrates the feasibility, and continuing challenges, of using in vitro bioactivity as a protective estimate of POD in screening-level assessments via a case study.

Extrapolating In Vitro Screening Assay Data for Thyroperoxidase Inhibition to Predict Serum Thyroid Hormones in the Rat

Thyroperoxidase (TPO) is an enzyme essential for thyroid hormone (TH) synthesis and a target site for a number of xenobiotics that disrupt TH homeostasis. An in vitro high-throughput screening assay for TPO inhibition, the Amplex UltraRed-TPO (AUR-TPO), has been used to screen the ToxCast chemical libraries for this action. Output from this assay would be most useful if it could be readily translated into an in vivo response, namely a reduction of TH in serum. To this end, the relationship between TPO inhibition in vitro and serum TH decreases was examined in rats exposed to 2 classic TPO inhibitors, propylthiouracil (PTU) and methimazole (MMI). Serum and gland PTU, MMI, and TH levels were quantified using tandem liquid chromatography mass spectrometry. Thyroperoxidase activity was determined in thyroid gland microsomes treated with PTU or MMI in vitro and ex vivo from thyroid gland microsomes prepared from exposed animals. A quantitative model was constructed by contrasting in vitro and ex vivo AUR-TPO results and the in vivo time-course and dose-response analysis. In vitro:ex vivo correlations of AUR-TPO outputs indicated that less than 30% inhibition of TPO in vitro was sufficient to reduce serum T4 by 20%, a degree of regulatory significance. Although further testing of model estimates using other TPO inhibitors is essential for verification of these initial findings, the results of this study provide a means to translate in vitro screening assay results into predictions of in vivo serum T4 changes to inform risk assessment.

Structure-based QSAR Models to Predict Repeat Dose Toxicity Points of Departure

Human health risk assessment for environmental chemical exposure is limited by a vast majority of chemicals with little or no experimental in vivo toxicity data. Data gap filling techniques, such as quantitative structure activity relationship (QSAR) models based on chemical structure information, can predict hazard in the absence of experimental data. Risk assessment requires identification of a quantitative point-of-departure (POD) value, the point on the dose-response curve that marks the beginning of a low-dose extrapolation. This study presents two sets of QSAR models to predict POD values (POD_QSAR) for repeat dose toxicity. For training and validation, a publicly available in vivo toxicity dataset for 3592 chemicals was compiled using the U.S. Environmental Protection Agency's Toxicity Value database (ToxValDB). The first set of QSAR models predict point-estimates of POD values (POD_QSAR) using structural and physicochemical descriptors for repeat dose study types and species combinations. A random forest QSAR model using study type and species as descriptors showed the best performance, with an external test set root mean square error (RMSE) of 0.71 log₁₀-mg/kg/day and coefficient of determination (R²) of 0.53. The second set of QSAR models predict the 95% confidence intervals for POD_QSAR using a constructed POD distribution with a mean equal to the median POD value and a standard deviation of 0.5 log₁₀-mg/kg/day, based on previously published typical study-to-study variability that may lead to uncertainty in model predictions. Bootstrap resampling of the pre-generated POD distribution was used to derive point-estimates and 95% confidence intervals for each POD prediction. Enrichment analysis to evaluate the accuracy of POD_QSAR showed that 80% of the 5% most potent chemicals were found in the top 20% of the most potent chemical predictions, suggesting that the repeat dose POD QSAR models presented here may help inform screening level human health risk assessments in the absence of other data.

Selecting a minimal set of androgen receptor assays for screening chemicals

Screening certain environmental chemicals for their ability to interact with endocrine targets, including the androgen receptor (AR), is an important global concern. We previously developed a model using a battery of eleven in vitro AR assays to predict in vivo AR activity. Here we describe a revised mathematical modeling approach that also incorporates data from newly available assays and demonstrate that subsets of assays can provide close to the same level of predictivity. These subset models are evaluated against the full model using 1820 chemicals, as well as in vitro and in vivo reference chemicals from the literature. Agonist batteries of as few as six assays and antagonist batteries of as few as five assays can yield balanced accuracies of 95% or better relative to the full model. Balanced accuracy for predicting reference chemicals is 100%. An approach is outlined for researchers to develop their own subset batteries to accurately detect AR activity using assays that map to the pathway of key molecular and cellular events involved in chemical-mediated AR activation and transcriptional activity. This work indicates in vitro bioactivity and in silico predictions that map to the AR pathway could be used in an integrated approach to testing and assessment for identifying chemicals that interact directly with the mammalian AR.

Variability in in vivo studies: Defining the upper limit of performance for predictions of systemic effect levels

New approach methodologies (NAMs) for chemical hazard assessment are often evaluated via comparison to animal studies; however, variability in animal study data limits NAM accuracy. The US EPA Toxicity Reference Database (ToxRefDB) enables consideration of variability in effect levels, including the lowest effect level (LEL) for a treatment-related effect and the lowest observable adverse effect level (LOAEL) defined by expert review, from subacute, subchronic, chronic, multi-generation reproductive, and developmental toxicity studies. The objectives of this work were to quantify the variance within systemic LEL and LOAEL values, defined as potency values for effects in adult or parental animals only, and to estimate the upper limit of NAM prediction accuracy. Multiple linear regression (MLR) and augmented cell means (ACM) models were used to quantify the total variance, and the fraction of variance in systemic LEL and LOAEL values explained by available study descriptors (e.g., administration route, study type). The MLR approach considered each study descriptor as an independent contributor to variance, whereas the ACM approach combined categorical descriptors into cells to define replicates. Using these approaches, total variance in systemic LEL and LOAEL values (in log₁₀-mg/kg/day units) ranged from 0.74 to 0.92. Unexplained variance in LEL and LOAEL values, approximated by the residual mean square error (MSE), ranged from 0.20-0.39. Considering subchronic, chronic, or developmental study designs separately resulted in similar values. Based on the relationship between MSE and R-squared for goodness-of-fit, the maximal R-squared may approach 55 to 73% for a NAM-based predictive model of systemic toxicity using these data as reference. The root mean square error (RMSE) ranged from 0.47 to 0.63 log₁₀-mg/kg/day, depending on dataset and regression approach, suggesting that a two-sided minimum prediction interval for systemic effect levels may have a width of 58 to 284-fold. These findings suggest quantitative considerations for building scientific confidence in NAM-based systemic toxicity predictions.

Respirometric Screening and Characterization of Mitochondrial Toxicants Within the ToxCast Phase I and II Chemical Libraries

Mitochondrial toxicity drives several adverse health outcomes. Current high-throughput screening assays for chemically induced mitochondrial toxicity typically measure changes to mitochondrial structure and may not detect known mitochondrial toxicants. We adapted a respirometric screening assay (RSA) measuring mitochondrial function to screen ToxCast chemicals in HepG2 cells using a tiered testing strategy. Of 1042 chemicals initially screened at a singlemaximal concentration, 243 actives were identified and rescreened at 7 concentrations. Concentration-response data for 3 respiration phases confirmed activity and indicated a mechanism for 193 mitochondrial toxicants: 149 electron transport chain inhibitors (ETCi), 15 uncouplers and 29 adenosine triphosphate synthase inhibitors. Subsequently, an electron flow assay was used to identify the target complex for 84 of the 149 ETCi. Sixty reference chemicals were used to compare the RSA to existing ToxCast and Tox21 mitochondrial toxicity assays. The RSA was most predictive (accuracy = 90%) of mitochondrial toxicity. The Tox21 mitochondrial membrane potential assay was also highly predictive (accuracy = 87%) of bioactivity but underestimated the potency of well-known ETCi and provided no mechanistic information. The tiered RSA approach accurately identifies and characterizes mitochondrial toxicants acting through diverse mechanisms and at a throughput sufficient to screen large chemical inventories. The electron flow assay provides additional confirmation and detailed mechanistic understanding for ETCi, the most common type of mitochondrial toxicants among ToxCast chemicals. The mitochondrial toxicity screening approach described herein may inform hazard assessment and the in vitro bioactive concentrations used to derive relevant doses for screening level chemical assessment using new approach methodologies.

Evaluation of Chemical Effects on Network Formation in Cortical Neurons Grown on Microelectrode Arrays

Thousands of chemicals to which humans are potentially exposed have not been evaluated for potential developmental neurotoxicity (DNT), driving efforts to develop a battery of in vitro screening approaches for DNT hazard. Here, 136 unique chemicals were evaluated for potential DNT hazard using a network formation assay (NFA) in cortical cells grown on microelectrode arrays. The effects of chemical exposure from 2 h postplating through 12 days in vitro (DIV) on network formation were evaluated at DIV 5, 7, 9, and 12, with cell viability assessed at DIV 12. Only 82 chemicals altered at least 1 network development parameter. Assay results were reproducible; 10 chemicals tested as biological replicates yielded qualitative results that were 100% concordant, with consistent potency values. Toxicological tipping points were determined for 58 chemicals and were similar to or lower than the lowest 50% effect concentrations (EC50) for all parameters. When EC50 and tipping point values from the NFA were compared to the range of potencies observed in ToxCast assays, the NFA EC50 values were less than the lower quartile for ToxCast assay potencies for a subset of chemicals, many of which are acutely neurotoxic in vivo. For 13 chemicals with available in vivo DNT data, estimated administered equivalent doses based on NFA results were similar to or lower than administered doses in vivo. Collectively, these results indicate that the NFA is sensitive to chemicals acting on nervous system function and will be a valuable contribution to an in vitro DNT screening battery.

Python BMDS: A Python interface library and web application for the canonical EPA dose-response modeling software

Several primary sources of publicly available, quantitative dose-response data from traditional toxicology study designs relevant to predictive toxicology applications are now available, including the redeveloped U.S. Environmental Protection Agency's Toxicity Reference Database (ToxRefDB v2.0), the Health Assessment Workspace Collaborative (HAWC), and the National Toxicology Program's Chemical Program's Chemical Effects in Biological Systems (CEBS). These resources provide effect level information but modeling these data to a curve may be more informative for predictive toxicology applications. Benchmark Dose Software (BMDS) has been recognized broadly and used for regulatory applications at multiple agencies. However, the current BMDS software was not amenable to modeling large datasets. Herein we describe development and use of a Python package that implements a wrapper around BMDS, a software that requires manual input in the dose-response modeling process (i.e., best-fitting model-selection, reporting, and dose-dropping). In the Python BMDS, users can select the BMDS version, customize model recommendation logic, and export summaries of the resultant BMDS output. Further, using the Python interface, a web-based application programming interface (API) has been developed for easy integration into other software systems, pipelines, or databases. Software utility was demonstrated via modeling nearly 28,000 datasets in ToxRefDB v2.0, re-creation of an existing, published large-scale analysis, and demonstration of usage in software such as CEBS and HAWC. Python BMDS enables rapid-batch processing of dose-response datasets using a modeling software with broad acceptance in the toxicology community, thereby providing an important tool for leveraging the publicly available quantitative toxicology data in a reproducible manner.

Progress in data interoperability to support computational toxicology and chemical safety evaluation

New approach methodologies (NAMs) in chemical safety evaluation are being explored to address the current public health implications of human environmental exposures to chemicals with limited or no data for assessment. For over a decade since a push toward "Toxicity Testing in the 21st Century," the field has focused on massive data generation efforts to inform computational approaches for preliminary hazard identification, adverse outcome pathways that link molecular initiating events and key events to apical outcomes, and high-throughput approaches to risk-based ratios of bioactivity and exposure to inform relative priority and safety assessment. Projects like the interagency Tox21 program and the US EPA ToxCast program have generated dose-response information on thousands of chemicals, identified and aggregated information from legacy systems, and created tools for access and analysis. The resulting information has been used to develop computational models as viable options for regulatory applications. This progress has introduced challenges in data management that are new, but not unique, to toxicology. Some of the key questions require critical thinking and solutions to promote semantic interoperability, including: (1) identification of bioactivity information from NAMs that might be related to a biological process; (2) identification of legacy hazard information that might be related to a key event or apical outcomes of interest; and, (3) integration of these NAM and traditional data for computational modeling and prediction of complex apical outcomes such as carcinogenesis. This work reviews a number of toxicology-related efforts specifically related to bioactivity and toxicological data interoperability based on the goals established by Findable, Accessible, Interoperable, and Reusable (FAIR) Data Principles. These efforts are essential to enable better integration of NAM and traditional toxicology information to support data-driven toxicology applications.

ToxRefDB version 2.0: Improved utility for predictive and retrospective toxicology analyses

The Toxicity Reference Database (ToxRefDB) structures information from over 5000 in vivo toxicity studies, conducted largely to guidelines or specifications from the US Environmental Protection Agency and the National Toxicology Program, into a public resource for training and validation of predictive models. Herein, ToxRefDB version 2.0 (ToxRefDBv2) development is described. Endpoints were annotated (e.g. required, not required) according to guidelines for subacute, subchronic, chronic, developmental, and multigenerational reproductive designs, distinguishing negative responses from untested. Quantitative data were extracted, and dose-response modeling for nearly 28,000 datasets from nearly 400 endpoints using Benchmark Dose (BMD) Modeling Software were generated and stored. Implementation of controlled vocabulary improved data quality; standardization to guideline requirements and cross-referencing with United Medical Language System (UMLS) connects ToxRefDBv2 observations to vocabularies linked to UMLS, including PubMed medical subject headings. ToxRefDBv2 allows for increased connections to other resources and has greatly enhanced quantitative and qualitative utility for predictive toxicology.

Evaluating Chemicals for Thyroid Disruption: Opportunities and Challenges with in Vitro Testing and Adverse Outcome Pathway Approaches

BACKGROUND

Extensive clinical and experimental research documents the potential for chemical disruption of thyroid hormone (TH) signaling through multiple molecular targets. Perturbation of TH signaling can lead to abnormal brain development, cognitive impairments, and other adverse outcomes in humans and wildlife. To increase chemical safety screening efficiency and reduce vertebrate animal testing, in vitro assays that identify chemical interactions with molecular targets of the thyroid system have been developed and implemented.

OBJECTIVES

We present an adverse outcome pathway (AOP) network to link data derived from in vitro assays that measure chemical interactions with thyroid molecular targets to downstream events and adverse outcomes traditionally derived from in vivo testing. We examine the role of new in vitro technologies, in the context of the AOP network, in facilitating consideration of several important regulatory and biological challenges in characterizing chemicals that exert effects through a thyroid mechanism.

DISCUSSION

There is a substantial body of knowledge describing chemical effects on molecular and physiological regulation of TH signaling and associated adverse outcomes. Until recently, few alternative nonanimal assays were available to interrogate chemical effects on TH signaling. With the development of these new tools, screening large libraries of chemicals for interactions with molecular targets of the thyroid is now possible. Measuring early chemical interactions with targets in the thyroid pathway provides a means of linking adverse outcomes, which may be influenced by many biological processes, to a thyroid mechanism. However, the use of in vitro assays beyond chemical screening is complicated by continuing limits in our knowledge of TH signaling in important life stages and tissues, such as during fetal brain development. Nonetheless, the thyroid AOP network provides an ideal tool for defining causal linkages of a chemical exerting thyroid-dependent effects and identifying research needs to quantify these effects in support of regulatory decision making. https://doi.org/10.1289/EHP5297.

Novel application of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene sets associated with disease: use case in breast carcinogenesis

Advances in technology within biomedical sciences have led to an inundation of data across many fields, raising new challenges in how best to integrate and analyze these resources. For example, rapid chemical screening programs like the US Environmental Protection Agency's ToxCast and the collaborative effort, Tox21, have produced massive amounts of information on putative chemical mechanisms where assay targets are identified as genes; however, systematically linking these hypothesized mechanisms with in vivo toxicity endpoints like disease outcomes remains problematic. Herein we present a novel use of normalized pointwise mutual information (NPMI) to mine biomedical literature for gene associations with biological concepts as represented by Medical Subject Headings (MeSH terms) in PubMed. Resources that tag genes to articles were integrated, then cross-species orthologs were identified using UniRef50 clusters. MeSH term frequency was normalized to reflect the MeSH tree structure, and then the resulting GeneID-MeSH associations were ranked using NPMI. The resulting network, called Entity MeSH Co-occurrence Network (EMCON), is a scalable resource for the identification and ranking of genes for a given topic of interest. The utility of EMCON was evaluated with the use case of breast carcinogenesis. Topics relevant to breast carcinogenesis were used to query EMCON and retrieve genes important to each topic. A breast cancer gene set was compiled through expert literature review (ELR) to assess performance of the search results. We found that the results from EMCON ranked the breast cancer genes from ELR higher than randomly selected genes with a recall of 0.98. Precision of the top five genes for selected topics was calculated as 0.87. This work demonstrates that EMCON can be used to link in vitro results to possible biological outcomes, thus aiding in generation of testable hypotheses for furthering understanding of biological function and the contribution of chemical exposures to disease.

High-Throughput H295R Steroidogenesis Assay: Utility as an Alternative and a Statistical Approach to Characterize Effects on Steroidogenesis

The U.S. Environmental Protection Agency Endocrine Disruptor Screening Program and the Organization for Economic Co-operation and Development (OECD) have used the human adrenocarcinoma (H295R) cell-based assay to predict chemical perturbation of androgen and estrogen production. Recently, a high-throughput H295R (HT-H295R) assay was developed as part of the ToxCast program that includes measurement of 11 hormones, including progestagens, corticosteroids, androgens, and estrogens. To date, 2012 chemicals have been screened at 1 concentration; of these, 656 chemicals have been screened in concentration-response. The objectives of this work were to: (1) develop an integrated analysis of chemical-mediated effects on steroidogenesis in the HT-H295R assay and (2) evaluate whether the HT-H295R assay predicts estrogen and androgen production specifically via comparison with the OECD-validated H295R assay. To support application of HT-H295R assay data to weight-of-evidence and prioritization tasks, a single numeric value based on Mahalanobis distances was computed for 654 chemicals to indicate the magnitude of effects on the synthesis of 11 hormones. The maximum mean Mahalanobis distance (maxmMd) values were high for strong modulators (prochloraz, mifepristone) and lower for moderate modulators (atrazine, molinate). Twenty-five of 28 reference chemicals used for OECD validation were screened in the HT-H295R assay, and produced qualitatively similar results, with accuracies of 0.90/0.75 and 0.81/0.91 for increased/decreased testosterone and estradiol production, respectively. The HT-H295R assay provides robust information regarding estrogen and androgen production, as well as additional hormones. The maxmMd from this integrated analysis may provide a data-driven approach to prioritizing lists of chemicals for putative effects on steroidogenesis.

ecancermedicalscience.. [PMID: 22276069 PMCID: PMC3239171 DOI: 10.3332/ecancer.2011.237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Details on the 28-year treatment history of a patient with an endocrine-responsive breast cancer are provided. She was originally diagnosed as having a T1N0M0 cancer after a modified radical mastectomy at age 41. Fifteen years later, in 1998, she presented with hemoptysis and pleuritic chest pain: a 10 cm right atrial tumor and estrogen receptor (ER) positive endobronchial and adjacent lung parenchyma adenocarcinoma were documented. Epithelial markers normalized as she manifested a partial response (PR) lasting 3 years with tamoxifen treatment. From 2001 to 2007 she benefitted from exemestane treatment. Upon progression in the previous lung area and left adrenal, exemestane withdrawal led to transient decrease in markers. Six months later (in July 2008), with growth in her adrenal tumor, laparoscopic adrenalectomy was performed: in addition to ER positivity, the tumor showed Her2 overexpression and amplification. She has subsequently had some control of disease with fulvestrant, letrozole + trastuzumab, and subsequently letrozole + lapatinib. In addition to the chronicity of disease, this history illustrates the expanding range of treatments available for endocrine-responsive breast cancer commensurate to our greater understanding of tumor biology.

A predictive data-driven framework for endocrine prioritization: a triazole fungicide case study

The US Environmental Protection Agency Endocrine Disruptor Screening Program (EDSP) is a tiered screening approach to determine the potential for a chemical to interact with estrogen, androgen, or thyroid hormone systems and/or perturb steroidogenesis. Use of high-throughput screening (HTS) to predict hazard and exposure is shifting the EDSP approach to (1) prioritization of chemicals for further screening; and (2) targeted use of EDSP Tier 1 assays to inform specific data needs. In this work, toxicology data for three triazole fungicides (triadimefon, propiconazole, and myclobutanil) were evaluated, including HTS results, EDSP Tier 1 screening (and other scientifically relevant information), and EPA guideline mammalian toxicology study data. The endocrine-related bioactivity predictions from HTS and information that satisfied the EDSP Tier 1 requirements were qualitatively concordant. Current limitations in the available HTS battery for thyroid and steroidogenesis pathways were mitigated by inclusion of guideline toxicology studies in this analysis. Similar margins (3-5 orders of magnitude) were observed between HTS-predicted human bioactivity and exposure values and between in vivo mammalian bioactivity and EPA chronic human exposure estimates for these products' registered uses. Combined HTS hazard and human exposure predictions suggest low priority for higher-tiered endocrine testing of these triazoles. Comparison with the mammalian toxicology database indicated that this HTS-based prioritization would have been protective for any potential in vivo effects that form the basis of current risk assessment for these chemicals. This example demonstrates an effective, human health protective roadmap for EDSP evaluation of pesticide active ingredients via prioritization using HTS and guideline toxicology information.

Tiered High-Throughput Screening Approach to Identify Thyroperoxidase Inhibitors Within the ToxCast Phase I and II Chemical Libraries

High-throughput screening for potential thyroid-disrupting chemicals requires a system of assays to capture multiple molecular-initiating events (MIEs) that converge on perturbed thyroid hormone (TH) homeostasis. Screening for MIEs specific to TH-disrupting pathways is limited in the U.S. Environmental Protection Agency ToxCast screening assay portfolio. To fill 1 critical screening gap, the Amplex UltraRed-thyroperoxidase (AUR-TPO) assay was developed to identify chemicals that inhibit TPO, as decreased TPO activity reduces TH synthesis. The ToxCast phase I and II chemical libraries, comprised of 1074 unique chemicals, were initially screened using a single, high concentration to identify potential TPO inhibitors. Chemicals positive in the single-concentration screen were retested in concentration-response. Due to high false-positive rates typically observed with loss-of-signal assays such as AUR-TPO, we also employed 2 additional assays in parallel to identify possible sources of nonspecific assay signal loss, enabling stratification of roughly 300 putative TPO inhibitors based upon selective AUR-TPO activity. A cell-free luciferase inhibition assay was used to identify nonspecific enzyme inhibition among the putative TPO inhibitors, and a cytotoxicity assay using a human cell line was used to estimate the cellular tolerance limit. Additionally, the TPO inhibition activities of 150 chemicals were compared between the AUR-TPO and an orthogonal peroxidase oxidation assay using guaiacol as a substrate to confirm the activity profiles of putative TPO inhibitors. This effort represents the most extensive TPO inhibition screening campaign to date and illustrates a tiered screening approach that focuses resources, maximizes assay throughput, and reduces animal use.