Profiling environmental chemicals for activity in the antioxidant response element signaling pathway using a high throughput screening approach

Effects of Chemicals in Reporter Gene Bioassays with Different Metabolic Activities Compared to Baseline Toxicity

High-throughput cell-based bioassays are used for chemical screening and risk assessment. Chemical transformation processes caused by abiotic degradation or metabolization can reduce the chemical concentration or, in some cases, lead to the formation of more toxic transformation products. Unaccounted loss processes may falsify the bioassay results. Capturing the formation and effects of transformation products is important for relating the in vitro effects to in vivo. Reporter gene cell lines are believed to have low metabolic activity, but inducibility of cytochrome P450 (CYP) enzymes has been reported. Baseline toxicity is the minimal toxicity a chemical can have and is caused by the incorporation of the chemical into cell membranes. In the present study, we improved an existing baseline toxicity model based on a newly defined critical membrane burden derived from freely dissolved effect concentrations, which are directly related to the membrane concentration. Experimental effect concentrations of 94 chemicals in three bioassays (AREc32, ARE-bla and GR-bla) were compared with baseline toxicity by calculating the toxic ratio (TR). CYP activities of all cell lines were determined by using fluorescence-based assays. Only ARE-bla showed a low basal CYP activity and inducibility and AREc32 showed a low inducibility. Overall cytotoxicity was similar in all three assays despite the different metabolic activities indicating that chemical metabolism is not relevant for the cytotoxicity of the tested chemicals in these assays. Up to 28 chemicals showed specific cytotoxicity with TR > 10 in the bioassays, but baseline toxicity could explain the effects of the majority of the remaining chemicals. Seven chemicals showed TR < 0.1 indicating inaccurate physicochemical properties or experimental artifacts like chemical precipitation, volatilization, degradation, or other loss processes during the in vitro bioassay. The new baseline model can be used not only to identify specific cytotoxicity mechanisms but also to identify potential problems in the experimental performance or evaluation of the bioassay and thus improve the quality of the bioassay data.

Prediction of drug-induced liver injury and cardiotoxicity using chemical structure and in vitro assay data

Drug-induced liver injury (DILI) and cardiotoxicity (DICT) are major adverse effects triggered by many clinically important drugs. To provide an alternative to in vivo toxicity testing, the U.S. Tox21 consortium has screened a collection of ∼10K compounds, including drugs in clinical use, against >70 cell-based assays in a quantitative high-throughput screening (qHTS) format. In this study, we compiled reference compound lists for DILI and DICT and compared the potential of Tox21 assay data with chemical structure information in building prediction models for human in vivo hepatotoxicity and cardiotoxicity. Models were built with four different machine learning algorithms (e.g., Random Forest, Naïve Bayes, eXtreme Gradient Boosting, and Support Vector Machine) and model performance was evaluated by calculating the area under the receiver operating characteristic curve (AUC-ROC). Chemical structure-based models showed reasonable predictive power for DILI (best AUC-ROC = 0.75 ± 0.03) and DICT (best AUC-ROC = 0.83 ± 0.03), while Tox21 assay data alone only showed better than random performance. DILI and DICT prediction models built using a combination of assay data and chemical structure information did not have a positive impact on model performance. The suboptimal predictive performance of the assay data is likely due to insufficient coverage of an adequately predictive number of toxicity mechanisms. The Tox21 consortium is currently expanding coverage of biological response space with additional assays that probe toxicologically important targets and under-represented pathways that may improve the prediction of in vivo toxicity such as DILI and DICT.

Application of ARE-reporter systems in drug discovery and safety assessment

The human body is continuously exposed to xenobiotics and internal or external oxidants. The health risk assessment of exogenous chemicals remains a complex and challenging issue. Alternative toxicological test methods have become an essential strategy for health risk assessment. As a core regulator of constitutive and inducible expression of antioxidant response element (ARE)-dependent genes, nuclear factor erythroid 2-related factor 2 (Nrf2) plays a critical role in maintaining cellular redox homeostasis. Consistent with the properties of Nrf2-mediated antioxidant response, Nrf2-ARE activity is a direct indicator of oxidative stress and thus has been used to identify and characterize oxidative stressors and redox modulators. To screen and distinguish chemicals or environmental insults that affect the cellular antioxidant activity and/or induce oxidative stress, various in vitro cell models expressing distinct ARE reporters with high-throughput and high-content properties have been developed. These ARE-reporter systems are currently widely applied in drug discovery and safety assessment. In the present review, we provide an overview of the basic structures and applications of various ARE-reporter systems employed for discovering Nrf2-ARE modulators and characterizing oxidative stressors.

Mechanism-driven modeling of chemical hepatotoxicity using structural alerts and an in vitro screening assay

Traditional experimental approaches to evaluate hepatotoxicity are expensive and time-consuming. As an advanced framework of risk assessment, adverse outcome pathways (AOPs) describe the sequence of molecular and cellular events underlying chemical toxicities. We aimed to develop an AOP that can be used to predict hepatotoxicity by leveraging computational modeling and in vitro assays. We curated 869 compounds with known hepatotoxicity classifications as a modeling set and extracted assay data from PubChem. The antioxidant response element (ARE) assay, which quantifies transcriptional responses to oxidative stress, showed a high correlation to hepatotoxicity (PPV=0.82). Next, we developed quantitative structure-activity relationship (QSAR) models to predict ARE activation for compounds lacking testing results. Potential toxicity alerts were identified and used to construct a mechanistic hepatotoxicity model. For experimental validation, 16 compounds in the modeling set and 12 new compounds were selected and tested using an in-house ARE-luciferase assay in HepG2-C8 cells. The mechanistic model showed good hepatotoxicity predictivity (accuracy = 0.82) for these compounds. Potential false positive hepatotoxicity predictions by only using ARE results can be corrected by incorporating structural alerts and vice versa. This mechanistic model illustrates a potential toxicity pathway for hepatotoxicity, and this strategy can be expanded to develop predictive models for other complex toxicities.

Identification of active or inactive agonists of tumor suppressor protein based on Tox21 library

Exposure of cells to xenobiotic human-made products can lead to genotoxicity and cause DNA damage. It is an urgent need to quickly identify the chemicals that cause DNA damage, and their toxicity should be predicted. In this study, recursive partitioning (RP), binary logistic regression, and one machine learning approach, namely, random forest (RF) classifier, were used to predict the active and inactive compounds of a total 5036 data based on the assay conducted by a β-lactamase reporter gene under control of the p53 response element (p53RE) from Tox21 library. Results show that the binary logistic regression model with a threshold of 0.5 has a high accuracy rate (83%) to distinguish active and inactive compounds. The RF classifier method has satisfactory results, with an accuracy rate (84.38%) approximately higher than that of binary logistic regression. The models established can identify compounds that induce DNA damage and activate p53, and provide a scientific basis for the risk assessment of organic chemicals in the environment.

Cell-Based Assays to Identify Modulators of Nrf2/ARE Pathway

The nuclear factor erythroid 2-related factor (Nrf2) and antioxidant response element (ARE) signaling pathway play an important role in the amelioration of cellular oxidative stress. Thus, assays that detect this pathway can be useful for identifying chemicals that induce or inhibit oxidative stress signaling. This chapter is to describe two cell-based Nrf2/ARE assays in a quantitative high-throughput screening (HTS) format to test a large collection of chemicals for oxidative stress induction ability. The assay descriptions involve cell handling, assay preparation, instrument usage, and assay procedure.

In Silico Approaches In Carcinogenicity Hazard Assessment: Current Status and Future Needs

Historically, identifying carcinogens has relied primarily on tumor studies in rodents, which require enormous resources in both money and time. In silico models have been developed for predicting rodent carcinogens but have not yet found general regulatory acceptance, in part due to the lack of a generally accepted protocol for performing such an assessment as well as limitations in predictive performance and scope. There remains a need for additional, improved in silico carcinogenicity models, especially ones that are more human-relevant, for use in research and regulatory decision-making. As part of an international effort to develop in silico toxicological protocols, a consortium of toxicologists, computational scientists, and regulatory scientists across several industries and governmental agencies evaluated the extent to which in silico models exist for each of the recently defined 10 key characteristics (KCs) of carcinogens. This position paper summarizes the current status of in silico tools for the assessment of each KC and identifies the data gaps that need to be addressed before a comprehensive in silico carcinogenicity protocol can be developed for regulatory use.

Application of In Vitro Metabolism Activation in High-Throughput Screening

In vitro methods which incorporate metabolic capability into the assays allow us to assess the activity of metabolites from their parent compounds. These methods can be applied into high-throughput screening (HTS) platforms, thereby increasing the speed to identify compounds that become active via the metabolism process. HTS was originally used in the pharmaceutical industry and now is also used in academic settings to evaluate biological activity and/or toxicity of chemicals. Although most chemicals are metabolized in our body, many HTS assays lack the capability to determine compound activity via metabolism. To overcome this problem, several in vitro metabolic methods have been applied to an HTS format. In this review, we describe in vitro metabolism methods and their application in HTS assays, as well as discuss the future perspectives of HTS with metabolic activity. Each in vitro metabolism method has advantages and disadvantages. For instance, the S9 mix has a full set of liver metabolic enzymes, but it displays high cytotoxicity in cell-based assays. In vitro metabolism requires liver fractions or the use of other metabolically capable systems, including primary hepatocytes or recombinant enzymes. Several newly developed in vitro metabolic methods, including HepaRG cells, three-dimensional (3D) cell models, and organ-on-a-chip technology, will also be discussed. These newly developed in vitro metabolism approaches offer significant progress in dissecting biological processes, developing drugs, and making toxicology studies quicker and more efficient.

Predicting oxidative stress induced by organic chemicals by using quantitative Structure-Activity relationship methods

Cellular exposure to xenobiotic human-made products will lead to oxidative stress that gives rise to DNA damage, as well as chemical or mechanical damage. Distinguishing the chemicals that will induce oxidative stress and predicting their toxicity is necessary. In the present study, 4270 compounds in the ARE-bla assay were investigated to predict active and inactive compounds by using simple algorithms, namely, recursive partitioning (RP) and binomial logistic regression, and to develop the quantitative structure-activity relationship (QSAR) models of chemicals that activate the ARE pathway to induce oxidative stress and exert toxic effects on cells. A decision tree based on scaffold-based fragments obtained through RP analysis showed the best identification accuracy. However, the overall identification accuracy of this model for active compounds was unsatisfactory due to limited fragments. Furthermore, a binomial logistic regression model was developed from 638 active compounds and 3632 inactive chemicals. The model with a cutoff of 0.15 could predict chemicals that were active or inactive with the prediction accuracy of 69.1%. Its area under the receiver operating characteristic (ROC) curve metric (AUROC) was 0.762, which indicated the acceptable predictive ability of this model. The parameters nBM (number of multiple bonds) and H% (percentage of H atom) played dominant roles in the prediction of the activity (inactive or active) of chemicals. A global QSAR model was developed to predict the toxicity of active chemicals. However, the model displayed an unsatisfactory result with R² = 0.316 and R²_ext = 0.090. Active chemicals were then classified on the basis of structure. A total of 79 compounds with carbon chains could be predicted with acceptable performance by using a QSAR model with six descriptors (R² = 0.722, R²_ext = 0.798, Q²_Loo = 0.654, Q²_Boot = 0.755, Q²_ext = 0.721). The simple models established here contribute to efforts on identification compounds inducing oxidative stress and provide the scientific basis for risk assessment to organisms in the environment.

Nrf2 in liver toxicology

Liver plays essential roles in the metabolism of many endogenous chemicals and exogenous toxicants. Mechanistic studies in liver have been at the forefront of efforts to probe the roles of bioactivation and detoxication of environmental toxins and toxicants in hepatotoxicity. Moreover, idiosyncratic hepatoxicity remains a key barrier in the clinical development of drugs. The now vast Nrf2 field emerged in part from biochemical and molecular studies on chemical inducers of hepatic detoxication enzymes and subsequent characterization of the modulation of drug/toxicant induced hepatotoxicities in mice through disruption of either Nrf2 or Keap1 genes. In general, loss of Nrf2 increases the sensitivity to such toxic chemicals, highlighting a central role of this transcription factor and its downstream target genes as a modifier to chemical stress. In this review, we summarize the impact of Nrf2 on the toxicology of multiple hepatotoxicants, and discuss efforts to utilize the Nrf2 response in predictive toxicology.

Baseline Toxicity and Volatility Cutoff in Reporter Gene Assays Used for High-Throughput Screening

Most studies using high-throughput in vitro cell-based bioassays tested chemicals up to a certain fixed concentration. It would be more appropriate to test up to concentrations predicted to elicit baseline toxicity because this is the minimal toxicity of every chemical. Baseline toxicity is also called narcosis and refers to nonspecific intercalation of chemicals in biological membranes, leading to loss of membrane structure and impaired functioning of membrane-related processes such as mitochondrial respiration. In cells, baseline toxicity manifests as cytotoxicity, which was quantified by a robust live-cell imaging method. Inhibitory concentrations for baseline toxicity varied by orders of magnitude between chemicals and were described by a simple quantitative structure activity relationship (QSAR) with the liposome-water partition constant as a sole descriptor. The QSAR equations were remarkably similar for eight reporter gene cell lines of different cellular origin, six of which were used in Tox21. Mass-balance models indicated constant critical membrane concentrations for all cells and all chemicals with a mean of 69 mmol·kg_lip^-1(95% CI: 49-89), which is in the same range as for bacteria and aquatic organisms and consistent with the theory of critical membrane burden of narcosis. The challenge of developing baseline QSARs for cell lines is that many confirmed baseline toxicants are rather volatile. We deduced from cytotoxicity experiments with semi-volatile chemicals that only chemicals with medium-air partition constants >10,000 L/L can be tested in standard robotic setups without appreciable loss of effect. Chemicals just below that cutoff showed crossover effects in neighboring wells, whereas the effects of chemicals with lower medium-air partition constants were plainly lost. Applying the "volatility cut-off" to >8000 chemicals tested in Tox21 indicated that approximately 20% of Tox21 chemicals could have partially been lost during the experiments. We recommend applying the baseline QSARs together with volatility cut-offs for experimental planning of reporter gene assays, that is, to dose only chemicals with medium-air partition constants >10,000 at concentrations up to the baseline toxicity level.

Antioxidant response elements: Discovery, classes, regulation and potential applications

Exposure to antioxidants and xenobiotics triggers the expression of a myriad of genes encoding antioxidant proteins, detoxifying enzymes, and xenobiotic transporters to offer protection against oxidative stress. This articulated universal mechanism is regulated through the cis-acting elements in an array of Nrf2 target genes called antioxidant response elements (AREs), which play a critical role in redox homeostasis. Though the Keap1/Nrf2/ARE system involves many players, AREs hold the key in transcriptional regulation of cytoprotective genes. ARE-mediated reporter constructs have been widely used, including xenobiotics profiling and Nrf2 activator screening. The complexity of AREs is brought by the presence of other regulatory elements within the AREs. The diversity in the ARE sequences not only bring regulatory selectivity of diverse transcription factors, but also confer functional complexity in the Keap1/Nrf2/ARE pathway. The different transcription factors either homodimerize or heterodimerize to bind the AREs. Depending on the nature of partners, they may activate or suppress the transcription. Attention is required for deeper mechanistic understanding of ARE-mediated gene regulation. The computational methods of identification and analysis of AREs are still in their infancy. Investigations are required to know whether epigenetics mechanism plays a role in the regulation of genes mediated through AREs. The polymorphisms in the AREs leading to oxidative stress related diseases are warranted. A thorough understanding of AREs will pave the way for the development of therapeutic agents against cancer, neurodegenerative, cardiovascular, metabolic and other diseases with oxidative stress.

Channel Interactions and Robust Inference for Ratiometric β-lactamase Assay Data: a Tox21 Library Analysis

Ratiometric β-lactamase (BLA) reporters are widely used to study transcriptional responses in a high-throughput screening (HTS) format. Typically, a ratio readout (background/target fluorescence) is used for toxicity assessment and structure-activity modeling efforts from BLA HTS data. This ratio readout may be confounded by channel-specific artifacts. To maximize the utility of BLA HTS data, we analyzed the relationship between individual channels and ratio readouts after fitting 10,000 chemical titration series screened in seven BLA stress-response assays from the Tox21 initiative. Similar to previous observations, we found that activity classifications based on BLA ratio readout alone are confounded by interference patterns for up to 85% (50 % on average) of active chemicals. Most Tox21 analyses adjust for this issue by evaluating target and ratio readout direction. In addition, we found that the potency and efficacy estimates derived from the ratio readouts may not represent the target channel effects and thus complicates chemical activity comparison. From these analyses we recommend a simpler approach using a direct evaluation of the target and background channels as well as the respective noise levels when using BLA data for toxicity assessment. This approach eliminates the channel interference issues and allows for straightforward chemical assessment and comparisons.

Omics-Based Platform for Studying Chemical Toxicity Using Stem Cells

The new strategy for chemical toxicity testing and modeling is to use in vitro human cell-based assays in conjunction with quantitative high-throughput screening (qHTS) technology, to identify molecular mechanisms and predict in vivo responses. Stem cells are more physiologically relevant than immortalized cell lines because of their unique proliferation and differentiation potentials. We established a robust two stem cells-two lineages assay system, encompassing human mesenchymal stem cells (hMSCs) along osteogenesis and human induced pluripotent stem cells (hiPSCs) along hepatogenesis. We performed qHTS phenotypic screening of LOPAC1280 and identified 38 preliminary hits for hMSCs. This was followed by validation of a selected number of hits and determination of their IC50 values and mechanistic studies of idarubicin and cantharidin treatments using proteomics and bioinformatics. In general, hiPSCs were more sensitive than hMSCs to chemicals, and differentiated progenies were less sensitive than their progenitors. We showed that chemical toxicity depends on both stem cell types and their differentiation stages. Proteomics identified and quantified over 3000 proteins for both stem cells. Bioinformatics identified apoptosis and G2/M as the top pathways conferring idarubicin toxicity. Our Omics-based assays of stem cells provide mechanistic insights into chemical toxicity and may help prioritize chemicals for in-depth toxicological evaluations.

Linear-No-Threshold Default Assumptions are Unwarranted for Cytotoxic Endpoints Independently Triggered by Ultrasensitive Molecular Switches

Crump's response in this issue to my critique of linear-no-threshold (LNT) default assumptions for noncancer and nongenotoxic cancer risks (Risk Analysis 2016; 36(3):589-604) is rebutted herein. Crump maintains that distinguishing between a low-dose linear dose response and a threshold dose response on the basis of dose-response data is impossible even for endpoints involving increased cytotoxicity. My rebuttal relies on descriptions and specific illustrations of two well-characterized ultrasensitive molecular switches that govern two key cytoprotective responses to cellular stress-heat shock response and antioxidant response element activation, respectively-each of which serve to suppress stress-induced apoptotic cell death unless overwhelmed. Because detailed dose-response data for each endpoint is shown to be J- or inverted-J-shaped with high confidence, and because independent pathways can explain background rates of apoptosis, LNT assumptions for this cytotoxic endpoint are unwarranted, at least in some cases and perhaps generally.

Exploring the oxidative stress response mechanism triggered by environmental water samples

Environmental waters can contain a wide range of micropollutants. Bioanalytical test batteries using assays indicative of different stages of cellular toxicity pathways, such as adaptive stress responses, have been applied to a range of water samples. Oxidative stress response assays have proven to be sensitive tools, but the mechanism by which water samples are inducing the oxidative stress response remains unclear because both electrophiles and reactive oxygen species (ROS) may activate the Nrf2-antioxidant response element (ARE) pathway. The current study aimed to explore the underlying mechanisms of the oxidative stress response triggered by exposure to surface water extracts previously shown to be active in the ARE GeneBLAzer oxidative stress response assay. ROS formation and changes in glutathione (GSH) concentration were assessed in human liver cells exposed to water extracts from a large river in addition to individual chemicals that were detected in these water extracts and reported to be active in the ARE GeneBLAzer assay in a previous study. Many of the surface water samples induced ROS formation and decreased the GSH to glutathione disulfide (GSSG) ratio, suggesting that the formation of ROS is an important mechanism. However, some of the most responsive samples in the ARE GeneBLAzer assay, as well as the individual chemicals, did not have a significant effect on either ROS formation or the GSH/GSSG ratio, suggesting a different mechanism. Antioxidants can also induce the Nrf2-ARE pathway and the ARE GeneBLAzer assay may also detect compounds that activate ARE by Nrf2-independent mechanisms, thus further research is required to characterise active chemicals in oxidative stress response assays. However, these tests are still useful for quantifying the integrated cellular response to multiple molecular initiating events and can be used complementary to assays indicative of specific effects, such as receptor-mediated assays.

Monohalogenated acetamide-induced cellular stress and genotoxicity are related to electrophilic softness and thiol/thiolate reactivity

Haloacetamides (HAMs) are cytotoxic, genotoxic, and mutagenic byproducts of drinking water disinfection. They are soft electrophilic compounds that form covalent bonds with the free thiol/thiolate in cysteine residues through an S_N2 reaction mechanism. Toxicity of the monohalogenated HAMs (iodoacetamide, IAM; bromoacetamide, BAM; or chloroacetamide, CAM) varied depending on the halogen substituent. The aim of this research was to investigate how the halogen atom affects the reactivity and toxicological properties of HAMs, measured as induction of oxidative/electrophilic stress response and genotoxicity. Additionally, we wanted to determine how well in silico estimates of electrophilic softness matched thiol/thiolate reactivity and in vitro toxicological endpoints. Each of the HAMs significantly induced nuclear Rad51 accumulation and ARE signaling activity compared to a negative control. The rank order of effect was IAM>BAM>CAM for Rad51, and BAM≈IAM>CAM for ARE. In general, electrophilic softness and in chemico thiol/thiolate reactivity provided a qualitative indicator of toxicity, as the softer electrophiles IAM and BAM were more thiol/thiolate reactive and were more toxic than CAM.

Individual expression features of GPX2, NQO1 and SQSTM1 transcript variants induced by hydrogen peroxide treatment in HeLa cells

Pathway activity assessment-based approaches are becoming highly influential in various fields of biology and medicine. However, these approaches mostly rely on analysis of mRNA expression, and total mRNA from a given locus is measured in the majority of cases. Notably, a significant portion of protein-coding genes produces more than one transcript. This biological fact is responsible for significant noise when changes in total mRNA transcription of a single gene are analyzed. The NFE2L2/AP-1 pathway is an attractive target for biomedical applications. To date, there is a lack of data regarding the agreement in expression of even classical target genes of this pathway. In the present paper we analyzed whether transcript variants of GPX2, NQO1 and SQSTM1 were characterized by individual features of expression when HeLa cells were exposed to pro-oxidative stimulation with hydrogen peroxide. We found that all the transcripts (10 in total) appeared to be significantly individually regulated under the conditions tested. We conclude that individual transcripts, rather than total mRNA, are best markers of pathway activation. We also discuss here some biological roles of individual transcript regulation.

Low-Dose Dose-Response for In Vitro Nrf2-ARE Activation in Human HepG2 Cells

Kelch-like ECH-associated protein 1 (Keap1), nuclear factor erythroid 2-like factor 2-related factor 2 (Nrf2), and the antioxidant response element (ARE) are interacting components of a master regulatory signaling pathway that coordinates redox homeostasis, cytoprotective responses, and shifts in stem cell state. This study reexamined detailed dose–response (DR) data reported for in vitro Nrf2-ARE activation in human hepatoblastoma HepG2 cell lines containing either a ARE-bla or ARE-luc reporter at 12 different concentrations of each of 15 chemicals. The normalized study data were combined among chemicals exhibiting a positive response, yielding n = 531 (179) DR data for 9 (7) chemicals using the ARE-bla (ARE-luc) assay. Three-parameter linear/kth-power regression fits obtained to each combined set of ARE-bla- or ARE-luc-assay response data provided good fits (R² = .99 or .91, respectively, P_fit > .99) that each incorporate a highly significant negative initial linear slope (P = 4 × 10⁻⁵ or .00025) and an overall J-shaped DR pattern. Results from this reanalysis of high-resolution ARE response data support the hypothesis that nonlinear ARE-mediated adaptive cellular responses to oxidative stress are governed by an ultrasensitive molecular switch.

Correlation between luminescence intensity and cytotoxicity in cell-based cytotoxicity assay using luciferase

The luciferase reporter assay has become one of the conventional methods for cytotoxicity evaluation. Typically, the decrease of luminescence expressed by a constitutive promoter is used as an index of cytotoxicity. However, to our knowledge, there have been no reports of the correlation between cytotoxicity and luminescence intensity. In this study, to accurately verify the correlation between them, beetle luciferase was stably expressed in human hepatoma HepG2 cells harboring the multi-integrase mouse artificial chromosome vector. We showed that the cytotoxicity assay using luciferase does not depend on the stability of luciferase protein and the kind of constitutive promoter. Next, HepG2 cells in which green-emitting beetle luciferase was expressed under the control of CAG promoter were exposed to 58 compounds. The luminescence intensity and cytotoxicity curves of cells exposed to 48 compounds showed similar tendencies, whereas those of cells exposed to 10 compounds did not do so, although the curves gradually approached each other with increasing exposure time. Finally, we demonstrated that luciferase expressed under the control of a constitutive promoter can be utilized both as an internal control reporter for normalizing a test reporter and for monitoring cytotoxicity when two kinds of luciferases are simultaneously used in the cytotoxicity assay.

Low-dose dose-response for reduced cell viability after exposure of human keratinocyte (HEK001) cells to arsenite

The in vitro arsenite (AsIII) cytotoxicity dose-response (DR) of human keratinocytes (HEK001) was examined at greater statistical resolution than ever previously reported using the MTT assay to determine cell viability. Fifty-four 96-well plates were treated with AsIII concentrations of 0.25, 0.5, 1, 2, 3, 4, 5, 7, 10, 15, 20, 25, or 30 μM. Because of unexpected variation in viability response patterns, a two-stage DR analysis was used in which data on plate-specific viability (%), estimated as 100% times the ratio of measured viability in exposed to unexposed cells, were fit initially to a generalized lognormal response function positing that HEK001 cells studied consisted of: a proportion P of relatively highly sensitive (HS) cells, a proportion Po of relatively resistant cells, and a remaining (1-P-Po) fraction of typical-sensitivity (TS) cells exhibiting the intermediate level of AsIII sensitivity characteristic of most cells in each assay. The estimated fractions P and Po were used to adjust data from all 54 plates (and from the 28 plates yielding the best fits) to reflect the condition that P = Po = 0 to provide detailed DR analysis specifically for TS cells. Four DR models fit to the combined adjusted data were each very predictive (R2 > 0.97) overall but were inconsistent with at least one of the data set examined (p < 10-5). Adjusted mean responses at ≤3 μM were approximately equal (p > 0.30) and exceeded 100% significance (p ≤ 10-6). A low-dose hormetic model provided the best fit to the combined adjusted data for TS cells (R2 = 0.995). Marked variability in estimates of P (the proportion of apparent HS cells) was unexpected, not readily explained, and warrants further study using additional cell lines and assay methods, and in vivo.

Effects of multi-component mixtures of polyaromatic hydrocarbons and heavy metal/loid(s) on Nrf2-antioxidant response element (ARE) pathway in ARE reporter-HepG2 cells

Exposure to polyaromatic hydrocarbons (PAHs) and heavy metal/loid(s) has been demonstrated to induce an oxidative stress response in mammalian cells. The combined effect of PAHs and heavy metal/loid(s) on the oxidative stress response has not been reported extensively. The Nrf2 antioxidant response pathway plays an important role in cellular antioxidant defense against oxidative stress-induced cell damage. In this study, we have determined the combined effect of four PAHs (benzo[a]pyrene (B[a]P), naphthalene (Nap), phenanthrene (Phe) and pyrene (Pyr)) and three heavy metal/loid(s) (arsenic (As), cadmium (Cd) and lead (Pb)) on the Nrf2 antioxidant pathway using the ARE reporter-HepG2 cell line. The mixture study was carried out for binary, ternary, quaternary and seven-component combinations of PAHs and heavy metal/loid(s). Initially, individual dose responses for the PAHs (B[a]P, Nap, Phe and Pyr) and heavy metal/loid(s) (As, Cd and Pb), as well as their respective concentrations that induced an induction ratio of 1.5 (ECIR1.5), were determined. The luciferase assay system was used to quantify the induction of the Nrf2 antioxidant pathway. The individual dose response study showed that both PAHs and heavy metal/loid(s) activated the Nrf2 antioxidant pathway in ARE reporter-HepG2 cells. Among these chemicals, Cd was the most potent inducer, followed by B[a]P and As. Based on the individual dose response findings, PAHs and heavy metal/loid(s) were mixed at equipotent ratios using a fixed concentration ratio, and the effects of the mixtures of PAHs and heavy metal/loid(s) (binary to seven-component) on the Nrf2 antioxidant pathway were determined. The mixture effects were predicted by using the concentration addition (CA) model. Overall, the results showed that the multi-component mixtures of PAHs and heavy metal/loid(s) induced an oxidative stress response in ARE reporter-HepG2 cells, and that the CA model is an appropriate model to predict the interaction effect of these selected mixtures. A human cell line-based reporter gene assay system was successfully used to determine the mixture effects of two groups of common contaminants on oxidative stress response pathway.

Mechanism Profiling of Hepatotoxicity Caused by Oxidative Stress Using Antioxidant Response Element Reporter Gene Assay Models and Big Data

BACKGROUND

Hepatotoxicity accounts for a substantial number of drugs being withdrawn from the market. Using traditional animal models to detect hepatotoxicity is expensive and time-consuming. Alternative in vitro methods, in particular cell-based high-throughput screening (HTS) studies, have provided the research community with a large amount of data from toxicity assays. Among the various assays used to screen potential toxicants is the antioxidant response element beta lactamase reporter gene assay (ARE-bla), which identifies chemicals that have the potential to induce oxidative stress and was used to test > 10,000 compounds from the Tox21 program.

OBJECTIVE

The ARE-bla computational model and HTS data from a big data source (PubChem) were used to profile environmental and pharmaceutical compounds with hepatotoxicity data.

METHODS

Quantitative structure-activity relationship (QSAR) models were developed based on ARE-bla data. The models predicted the potential oxidative stress response for known liver toxicants when no ARE-bla data were available. Liver toxicants were used as probe compounds to search PubChem Bioassay and generate a response profile, which contained thousands of bioassays (> 10 million data points). By ranking the in vitro-in vivo correlations (IVIVCs), the most relevant bioassay(s) related to hepatotoxicity were identified.

RESULTS

The liver toxicants profile contained the ARE-bla and relevant PubChem assays. Potential toxicophores for well-known toxicants were created by identifying chemical features that existed only in compounds with high IVIVCs.

CONCLUSION

Profiling chemical IVIVCs created an opportunity to fully explore the source-to-outcome continuum of modern experimental toxicology using cheminformatics approaches and big data sources.

CITATION

Kim MT, Huang R, Sedykh A, Wang W, Xia M, Zhu H. 2016. Mechanism profiling of hepatotoxicity caused by oxidative stress using antioxidant response element reporter gene assay models and big data. Environ Health Perspect 124:634-641; http://dx.doi.org/10.1289/ehp.1509763.

High-Throughput Phenotypic Screening of Human Astrocytes to Identify Compounds That Protect Against Oxidative Stress

Astrocytes are the predominant cell type in the nervous system and play a significant role in maintaining neuronal health and homeostasis. Recently, astrocyte dysfunction has been implicated in the pathogenesis of many neurodegenerative diseases, including Alzheimer's disease, Parkinson's disease, Huntington's disease, and amyotrophic lateral sclerosis. Astrocytes are thus an attractive new target for drug discovery for neurological disorders. Using astrocytes differentiated from human embryonic stem cells, we have developed an assay to identify compounds that protect against oxidative stress, a condition associated with many neurodegenerative diseases. This phenotypic oxidative stress assay has been optimized for high-throughput screening in a 1,536-well plate format. From a screen of approximately 4,100 bioactive tool compounds and approved drugs, we identified a set of 22 that acutely protect human astrocytes from the consequences of hydrogen peroxide-induced oxidative stress. Nine of these compounds were also found to be protective of induced pluripotent stem cell-differentiated astrocytes in a related assay. These compounds are thought to confer protection through hormesis, activating stress-response pathways and preconditioning astrocytes to handle subsequent exposure to hydrogen peroxide. In fact, four of these compounds were found to activate the antioxidant response element/nuclear factor-E2-related factor 2 pathway, a protective pathway induced by toxic insults. Our results demonstrate the relevancy and utility of using astrocytes differentiated from human stem cells as a disease model for drug discovery and development.

SIGNIFICANCE

Astrocytes play a key role in neurological diseases. Drug discovery efforts that target astrocytes can identify novel therapeutics. Human astrocytes are difficult to obtain and thus are challenging to use for high-throughput screening, which requires large numbers of cells. Using human embryonic stem cell-derived astrocytes and an optimized astrocyte differentiation protocol, it was possible to screen approximately 4,100 compounds in titration to identify 22 that are cytoprotective of astrocytes. This study is the largest-scale high-throughput screen conducted using human astrocytes, with a total of 17,536 data points collected in the primary screen. The results demonstrate the relevancy and utility of using astrocytes differentiated from human stem cells as a disease model for drug discovery and development.

Fingerprinting the reactive toxicity pathways of 50 drinking water disinfection by-products

A set of nine in vitro cellular bioassays indicative of different stages of the cellular toxicity pathway was applied to 50 disinfection by-products (DBPs) to obtain a better understanding of the commonalities and differences in the molecular mechanisms of reactive toxicity of DBPs. An Eschericia coli test battery revealed reactivity towards proteins/peptides for 64% of the compounds. 98% activated the NRf2-mediated oxidative stress response and 68% induced an adaptive stress response to genotoxic effects as indicated by the activation of the tumor suppressor protein p53. All DBPs reactive towards DNA in the E. coli assay and activating p53 also induced oxidative stress, confirming earlier studies that the latter could trigger DBP's carcinogenicity. The energy of the lowest unoccupied molecular orbital ELUMO as reactivity descriptor was linearly correlated with oxidative stress induction for trihalomethanes (r(2)=0.98) and haloacetamides (r(2)=0.58), indicating that potency of these DBPs is connected to electrophilicity. However, the descriptive power was poor for haloacetic acids (HAAs) and haloacetonitriles (r(2) (<) 0.06). For HAAs, we additionally accounted for speciation by including the acidity constant with ELUMO in a two-parameter multiple linear regression model. This increased r(2) to >0.80, indicating that HAAs' potency is connected to both, electrophilicity and speciation. Based on the activation of oxidative stress response and the soft electrophilic character of most tested DBPs we hypothesize that indirect genotoxicity-e.g., through oxidative stress induction and/or enzyme inhibition-is more plausible than direct DNA damage for most investigated DBPs. The results provide not only a mechanistic understanding of the cellular effects of DBPs but the effect concentrations may also serve to evaluate mixture effects of DBPs in water samples.

Modelling the Tox21 10 K chemical profiles for in vivo toxicity prediction and mechanism characterization

Target-specific, mechanism-oriented in vitro assays post a promising alternative to traditional animal toxicology studies. Here we report the first comprehensive analysis of the Tox21 effort, a large-scale in vitro toxicity screening of chemicals. We test ∼ 10,000 chemicals in triplicates at 15 concentrations against a panel of nuclear receptor and stress response pathway assays, producing more than 50 million data points. Compound clustering by structure similarity and activity profile similarity across the assays reveals structure-activity relationships that are useful for the generation of mechanistic hypotheses. We apply structural information and activity data to build predictive models for 72 in vivo toxicity end points using a cluster-based approach. Models based on in vitro assay data perform better in predicting human toxicity end points than animal toxicity, while a combination of structural and activity data results in better models than using structure or activity data alone. Our results suggest that in vitro activity profiles can be applied as signatures of compound mechanism of toxicity and used in prioritization for more in-depth toxicological testing.

Cell-Based Assay for Identifying the Modulators of Antioxidant Response Element Signaling Pathway

The antioxidant response element (ARE) signaling pathway plays an important role in the amelioration of cellular oxidative stress. Thus, assays that detect this pathway can be useful for identifying chemicals that induce or inhibit oxidative stress signaling. The focus of this chapter is to describe a cell-based ARE assay in a quantitative high-throughput screening (qHTS) format to test a large collection of compounds that induce nuclear factor erythroid 2-related factor (Nrf2)/ARE signaling. The assay is described through cell handling, assay preparation, and instrument usage.

Adaptive stress response pathways induced by environmental mixtures of bioaccumulative chemicals in dugongs

To address the poorly understood mixture effects of chemicals in the marine mammal dugong, we coupled equilibrium-based passive sampling in blubber to a range of in vitro bioassays for screening mixtures of bioaccumulative chemicals. The modes of action included early effect indicators along important toxicity pathways, such as induction of xenobiotic metabolism, and some integrative indicators downstream of the molecular initiating event, such as adaptive stress responses. Activation of aryl hydrocarbon receptor (AhR) and Nrf2-mediated oxidative stress response were found to be the most prominent effects, while the p53-mediated DNA damage response and NF-κB-mediated response to inflammation were not significantly affected. Although polychlorinated dibenzo-p-dioxins (PCDDs) quantified in the samples accounted for the majority of AhR-mediated activity, PCDDs explained less than 5% of the total oxidative stress response, despite their known ability to activate this pathway. Altered oxidative stress response was observed with both individual chemicals and blubber extracts subject to metabolic activation by rat liver S9 fraction. Metabolic activation resulted in both enhanced and reduced toxicity, suggesting the relevance and utility of incorporating metabolic enzymes into in vitro bioassays. Our approach provides a first insight into the burden of toxicologically relevant bioaccumulative chemical mixtures in dugongs and can be applied to lipid tissue of other wildlife species.

2013 Lush Science Prize

The annual Lush Science Prize is designed to reward outstanding contributions to 21st Century Toxicology Research. A Background Paper is prepared each year prior to the judging process, in order to provide the judging panel with a brief overview of current developments in the field of Replacement alternatives, particularly those relevant to the concept of toxicity pathways. The Background Paper includes information on some key institutional developments in the area--such as the OECD's Adverse Outcome Pathway Project, the Hamner Institute's work, and the Human Toxome Project, and on the phenomenon of collaborative computer systems relevant to the field. From the literature review that was also performed as part of the background research, the two papers receiving the highest score were recommended for consideration by the judges for the 2013 Science Prize.

Free radical generation from an aniline derivative in HepG2 cells: a possible captodative effect

Xenobiotic metabolism can induce the generation of protein radicals, which are believed to play an important role in the toxicity of chemicals and drugs. It is therefore important to identify chemical structures capable of inducing macromolecular free radical formation in living cells. In this study, we evaluated the ability of four structurally related environmental chemicals, aniline, nitrosobenzene, N,N-dimethylaniline, and N,N-dimethyl-4-nitrosoaniline (DMNA), to induce free radicals and cellular damage in the hepatoma cell line HepG2. Cytotoxicity was assessed using lactate dehydrogenase assays, and morphological changes were observed using phase contrast microscopy. Protein free radicals were detected by immuno-spin trapping using in-cell western experiments and confocal microscopy to determine the subcellular locale of free radical generation. DMNA induced free radical generation, lactate dehydrogenase release, and morphological changes in HepG2 cells, whereas aniline, nitrosobenzene, N,N-dimethylaniline did not. Confocal microscopy showed that DMNA induced free radical generation mainly in the cytosol. Preincubation of HepG2 cells with N-acetylcysteine and 2,2'-dipyridyl significantly prevented free radical generation on subsequent incubation with DMNA, whereas preincubation with apocynin and dimethyl sulfoxide had no effect. These results suggest that DMNA is metabolized to reactive free radicals capable of generating protein radicals which may play a critical role in DMNA toxicity. We propose that the captodative effect, the combined action of the electron-releasing dimethylamine substituent, and the electron-withdrawing nitroso substituent, leads to a thermodynamically stabilized radical, facilitating enhanced protein radical formation by DMNA.

Improving the human hazard characterization of chemicals: a Tox21 update

BACKGROUND

In 2008, the National Institute of Environmental Health Sciences/National Toxicology Program, the U.S. Environmental Protection Agency's National Center for Computational Toxicology, and the National Human Genome Research Institute/National Institutes of Health Chemical Genomics Center entered into an agreement on "high throughput screening, toxicity pathway profiling, and biological interpretation of findings." In 2010, the U.S. Food and Drug Administration (FDA) joined the collaboration, known informally as Tox21.

OBJECTIVES

The Tox21 partners agreed to develop a vision and devise an implementation strategy to shift the assessment of chemical hazards away from traditional experimental animal toxicology studies to one based on target-specific, mechanism-based, biological observations largely obtained using in vitro assays.

DISCUSSION

Here we outline the efforts of the Tox21 partners up to the time the FDA joined the collaboration, describe the approaches taken to develop the science and technologies that are currently being used, assess the current status, and identify problems that could impede further progress as well as suggest approaches to address those problems.

CONCLUSION

Tox21 faces some very difficult issues. However, we are making progress in integrating data from diverse technologies and end points into what is effectively a systems-biology approach to toxicology. This can be accomplished only when comprehensive knowledge is obtained with broad coverage of chemical and biological/toxicological space. The efforts thus far reflect the initial stage of an exceedingly complicated program, one that will likely take decades to fully achieve its goals. However, even at this stage, the information obtained has attracted the attention of the international scientific community, and we believe these efforts foretell the future of toxicology.

Human cell toxicogenomic analysis linking reactive oxygen species to the toxicity of monohaloacetic acid drinking water disinfection byproducts

Chronic exposure to drinking water disinfection byproducts has been linked to adverse health risks. The monohaloacetic acids (monoHAAs) are generated as byproducts during the disinfection of drinking water and are cytotoxic, genotoxic, mutagenic, and teratogenic. Iodoacetic acid toxicity was mitigated by antioxidants, suggesting the involvement of oxidative stress. Other monoHAAs may share a similar mode of action. Each monoHAA generated a significant concentration-response increase in the expression of a β-lactamase reporter under the control of the antioxidant response element (ARE). The monoHAAs generated oxidative stress with a rank order of iodoacetic acid (IAA) > bromoacetic acid (BAA) ≫ chloroacetic acid (CAA); this rank order was observed with other toxicological end points. Toxicogenomic analysis was conducted with a nontransformed human intestinal epithelial cell line (FHs 74 Int). Exposure to the monoHAAs altered the transcription levels of multiple oxidative stress responsive genes, indicating that each exposure generated oxidative stress. The transcriptome profiles showed an increase in thioredoxin reductase 1 (TXNRD1) and sulfiredoxin (SRXN1), suggesting peroxiredoxin proteins had been oxidized during monoHAA exposures. Three possible sources of reactive oxygen species were identified, the hypohalous acid generating peroxidase enzymes lactoperoxidase (LPO) and myeloperoxidase (MPO), nicotinamide adenine dinucleotide phosphate (NADPH)-dependent oxidase 5 (NOX5), and PTGS2 (COX-2) mediated arachidonic acid metabolism. Each monoHAA exposure caused an increase in COX-2 mRNA levels. These data provide a functional association between monoHAA exposure and adverse health outcomes such as oxidative stress, inflammation, and cancer.

Water quality assessment using the AREc32 reporter gene assay indicative of the oxidative stress response pathway

The reporter gene assay AREc32 is based on the induction of the Nrf2 mediated oxidative stress response pathway in the human breast cancer cell line MCF7, where eight copies of the antioxidant response element (ARE) are linked to a reporter gene encoding for luciferase. The Nrf2-ARE pathway is responsive to many chemicals that cause oxidative stress, among them a large number of pesticides and skin irritants. We adopted and validated the AREc32 bioassay for water quality testing. tert-Butylhydroquinone served as the positive control, phenol as the negative control and other reactive chemicals were assessed for their specificity. An environmentally relevant reference chemical, benzo(a)pyrene was the most potent inducer of all tested chemicals. The concentration causing an induction ratio (IR) of 1.5 (EC(IR1.5)) was chosen as the effect benchmark value. The assay was applied to 21 water samples ranging from sewage to drinking water, including secondary treatment and various tertiary treatment options (ozonation, biologically activated carbon filtration, membrane filtration, reverse osmosis, advanced oxidation, chlorination, chloramination). The samples were enriched by solid phase extraction. In most samples the oxidative stress response was far more sensitive than cytotoxicity. The primary and secondary treated effluent exceeded the effect threshold IR 1.5 at a relative enrichment factor (REF) of 1, i.e., the native samples were active. All tertiary treated samples were less potent and their EC(IR1.5) lay between REF 1 and 10. The Nrf2 pathway was induced at a REF of approximately 10 for surface waters and drinking water, and above this enrichment cytotoxicity took over in most samples and quenched the induction. The blank (ultrapure water run through the sample enrichment process) was cytotoxic at an REF of 100, which is the limit of concentrations range that can be evaluated. Treatment typically decreased both the cytotoxicity and oxidative stress response apart from drinking water treatment where chlorination caused an increase in oxidative stress response, presumably due to the formation of disinfection by-products. This study demonstrates the relevance and applicability of the oxidative stress response pathway for water quality monitoring.