Building predictive in vitro pulmonary toxicity assays using high-throughput imaging and artificial intelligence

Applying new approach methodologies to assess next-generation tobacco and nicotine products

In vitro toxicology research has accelerated with the use of in silico, computational approaches and human in vitro tissue systems, facilitating major improvements evaluating the safety and health risks of novel consumer products. Innovation in molecular and cellular biology has shifted testing paradigms, with less reliance on low-throughput animal data and greater use of medium- and high-throughput in vitro cellular screening approaches. These new approach methodologies (NAMs) are being implemented in other industry sectors for chemical testing, screening candidate drugs and prototype consumer products, driven by the need for reliable, human-relevant approaches. Routine toxicological methods are largely unchanged since development over 50 years ago, using high-doses and often employing in vivo testing. Several disadvantages are encountered conducting or extrapolating data from animal studies due to differences in metabolism or exposure. The last decade saw considerable advancement in the development of in vitro tools and capabilities, and the challenges of the next decade will be integrating these platforms into applied product testing and acceptance by regulatory bodies. Governmental and validation agencies have launched and applied frameworks and "roadmaps" to support agile validation and acceptance of NAMs. Next-generation tobacco and nicotine products (NGPs) have the potential to offer reduced risks to smokers compared to cigarettes. These include heated tobacco products (HTPs) that heat but do not burn tobacco; vapor products also termed electronic nicotine delivery systems (ENDS), that heat an e-liquid to produce an inhalable aerosol; oral smokeless tobacco products (e.g., Swedish-style snus) and tobacco-free oral nicotine pouches. With the increased availability of NGPs and the requirement of scientific studies to support regulatory approval, NAMs approaches can supplement the assessment of NGPs. This review explores how NAMs can be applied to assess NGPs, highlighting key considerations, including the use of appropriate in vitro model systems, deploying screening approaches for hazard identification, and the importance of test article characterization. The importance and opportunity for fit-for-purpose testing and method standardization are discussed, highlighting the value of industry and cross-industry collaborations. Supporting the development of methods that are accepted by regulatory bodies could lead to the implementation of NAMs for tobacco and nicotine NGP testing.

Circumvention of Gefitinib Resistance by Repurposing Flunarizine via Histone Deacetylase Inhibition

Gefitinib is an epidermal growth factor receptor tyrosine kinase inhibitor (EGFR TKI) for treating advanced non-small cell lung cancer (NSCLC). However, drug resistance seriously impedes the clinical efficacy of gefitinib. This study investigated the repositioning of the non-oncology drug capable of inhibiting histone deacetylases (HDACs) to overcome gefitinib resistance. A few drug candidates were identified using the in silico repurposing tool "DRUGSURV" and tested for HDAC inhibition. Flunarizine, originally indicated for migraine prophylaxis and vertigo treatment, was selected for detailed investigation in NSCLC cell lines harboring a range of different gefitinib resistance mechanisms (EGFR T790M, KRAS G12S, MET amplification, or PTEN loss). The circumvention of gefitinib resistance by flunarizine was further demonstrated in an EGFR TKI (erlotinib)-refractory patient-derived tumor xenograft (PDX) model in vivo. The acetylation level of cellular histone protein was increased by flunarizine in a concentration- and time-dependent manner. Among the NSCLC cell lines evaluated, the extent of gefitinib resistance circumvention by flunarizine was found to be the most pronounced in EGFR T790M-bearing H1975 cells. The gefitinib-flunarizine combination was shown to induce the apoptotic protein Bim but reduce the antiapoptotic protein Bcl-2, which apparently circumvented gefitinib resistance. The induction of Bim by flunarizine was accompanied by an increase in the histone acetylation and E2F1 interaction with the BIM gene promoter. Flunarizine was also found to upregulate E-cadherin but downregulate the vimentin expression, which subsequently inhibited cancer cell migration and invasion. Importantly, flunarizine was also shown to significantly potentiate the tumor growth suppressive effect of gefitinib in EGFR TKI-refractory PDX in vivo. The findings advocate for the translational application of flunarizine to circumvent gefitinib resistance in the clinic.

Repurposing of triamterene as a histone deacetylase inhibitor to overcome cisplatin resistance in lung cancer treatment

PURPOSE

Cisplatin is the core chemotherapeutic drug used for first-line treatment of advanced non-small cell lung cancer (NSCLC). However, drug resistance is severely hindering its clinical efficacy. This study investigated the circumvention of cisplatin resistance by repurposing non-oncology drugs with putative histone deacetylase (HDAC) inhibitory effect.

METHODS

A few clinically approved drugs were identified by a computational drug repurposing tool called "DRUGSURV" and evaluated for HDAC inhibition. Triamterene, originally indicated as a diuretic, was chosen for further investigation in pairs of parental and cisplatin-resistant NSCLC cell lines. Sulforhodamine B assay was used to evaluate cell proliferation. Western blot analysis was performed to examine histone acetylation. Flow cytometry was used to examine apoptosis and cell cycle effects. Chromatin immunoprecipitation was conducted to investigate the interaction of transcription factors to the promoter of genes regulating cisplatin uptake and cell cycle progression. The circumvention of cisplatin resistance by triamterene was further verified in a patient-derived tumor xenograft (PDX) from a cisplatin-refractory NSCLC patient.

RESULTS

Triamterene was found to inhibit HDACs. It was shown to enhance cellular cisplatin accumulation and potentiate cisplatin-induced cell cycle arrest, DNA damage, and apoptosis. Mechanistically, triamterene was found to induce histone acetylation in chromatin, thereby reducing the association of HDAC1 but promoting the interaction of Sp1 with the gene promoter of hCTR1 and p21. Triamterene was further shown to potentiate the anti-cancer effect of cisplatin in cisplatin-resistant PDX in vivo.

CONCLUSION

The findings advocate further clinical evaluation of the repurposing use of triamterene to overcome cisplatin resistance.

Risk assessment based on dose-responsive and time-responsive genes to build PLS-DA models for exogenously induced lung injury

Xenobiotics can easily harm human lungs owing to the openness of the respiratory system. Identifying pulmonary toxicity remains challenging owing to several reasons: 1) no biomarkers for pulmonary toxicity are available that might help to detect lung injury; 2) traditional animal experiments are time-consuming; 3) traditional detection methods solely focus on poisoning accidents; 4) analytical chemistry methods hardly achieve universal detection. An in vitro testing system able to identify the pulmonary toxicity of contaminants from food, the environment, and drugs is urgently needed. Compounds are virtually infinite, whereas toxicological mechanisms are countable. Therefore, universal methods to identify and predict the risks of contaminants can be designed based on these well-known toxicity mechanisms. In this study, we established a dataset based on transcriptome sequencing of A549 cells upon treatment with different compounds. The representativeness of our dataset was analyzed using bioinformatics methods. Artificial intelligence methods, namely partial least squares discriminant analysis (PLS-DA) models, were employed for toxicity prediction and toxicant identification. The developed model predicted the pulmonary toxicity of compounds with a 92 % accuracy. These models were submitted to an external validation using highly heterogeneous compounds, which supported the accuracy and robustness of our developed methodology. This assay exhibits universal potential applications for water quality monitoring, crop pollution detection, food and drug safety evaluation, as well as chemical warfare agent detection.

Reference compounds for characterizing cellular injury in high-content cellular morphology assays

Robust, generalizable approaches to identify compounds efficiently with undesirable mechanisms of action in complex cellular assays remain elusive. Such a process would be useful for hit triage during high-throughput screening and, ultimately, predictive toxicology during drug development. Here we generate cell painting and cellular health profiles for 218 prototypical cytotoxic and nuisance compounds in U-2 OS cells in a concentration-response format. A diversity of compounds that cause cellular damage produces bioactive cell painting morphologies, including cytoskeletal poisons, genotoxins, nonspecific electrophiles, and redox-active compounds. Further, we show that lower quality lysine acetyltransferase inhibitors and nonspecific electrophiles can be distinguished from more selective counterparts. We propose that the purposeful inclusion of cytotoxic and nuisance reference compounds such as those profiled in this resource will help with assay optimization and compound prioritization in complex cellular assays like cell painting.

Linking chemicals, genes and morphological perturbations to diseases

The progress in image-based high-content screening technology has facilitated high-throughput phenotypic profiling notably the quantification of cell morphology perturbation by chemicals. However, understanding the mechanism of action of a chemical and linking it to cell morphology and phenotypes remains a challenge in drug discovery. In this study, we intended to integrate molecules that induced transcriptomic perturbations and cellular morphological changes into a biological network in order to assess chemical-phenotypic relationships in humans. Such a network was enriched with existing disease information to suggest molecular and cellular profiles leading to phenotypes. Two datasets were used for this study. Firstly, we used the "Cell Painting morphological profiling assay" dataset, composed of 30,000 compounds tested on human osteosarcoma cells (named U2OS). Secondly, we used the "L1000 mRNA profiling assay" dataset, a collection of transcriptional expression data from cultured human cells treated with approximately 20,000 bioactive small molecules from the Library of Integrated Network-based Cellular Signatures (LINCS). Furthermore, pathways, gene ontology terms and disease enrichments were performed on the transcriptomics data. Overall, our study makes it possible to develop a biological network combining chemical-gene-pathway-morphological perturbation and disease relationships. It contains an ensemble of 9989 chemicals, 732 significant morphological features and 12,328 genes. Through diverse examples, we demonstrated that some drugs shared similar genes, pathways and morphological profiles that, taken together, could help in deciphering chemical-phenotype observations.

Chemical Interactions and Cytotoxicity of Terpene and Diluent Vaping Ingredients

Vaping devices have risen in popularity since their inception in 2007. The practice involves using a variety of commercially available devices. Internal heating systems in devices aerosolize e-liquid formulations of complex mixtures including an active ingredient (e.g., THC, CBD, and nicotine), diluents (or cutting agents), solvents, and flavoring agents (e.g., terpenes and aldehydes). The vaping toxicology literature consists of cytotoxicity studies of individual chemicals and commercial formulas. Because of the variation of e-liquid composition, there is a limited understanding of the toxicity of ingredient combinations. This study analyzed the cytotoxic effects after exposure to individual and binary mixtures of a representative terpene (+-R-limonene) and diluent (triethyl citrate) on human lung cell models. Data were analyzed to determine the effects of 97:3 and 80:20% v/v (triethyl citrate/limonene) binary mixtures. BEAS-2B cells, a bronchial epithelial cell, and A549 cells, a type II alveolar epithelial cell, served as models for comparison. LC₅₀ values were calculated and isobolograms were used to assess chemical interactions. Results show that limonene was more cytotoxic than triethyl citrate. Isobolographic analyses confirmed that the 97:3% v/v mixture resulted in an antagonistic chemical interaction. The 80:20% v/v mixture resulted in a similar result. Further testing of different ratios of binary mixtures is needed for chemical interaction screening to inform safety assessments.

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.

Group VIII Metal Carbonyl Cluster-Boronic Acid Conjugates: Cytotoxicity and Mode of Action Studies

A set of metal carbonyl cluster-boronic acid conjugates of the group VIII metals (Fe, Ru, and Os) were synthesized and their antiproliferative effects measured against two breast cancer cell lines (MCF-7 and MDA-MB-231) and a noncancerous breast epithelial (MCF-10A) cell line. The cytotoxicity followed the order Ru > Os > Fe for the MDA-MB-231 cells, although the latter two exhibited similar cytotoxicity against MCF-7 and MCF-10A cells. The osmium species {Os₃(CO)₁₀(μ-H)[μ-SC₆H₄-p-B(OH)₂]} (2) could be chemically oxidized to its hydroxy analogue [Os₃(CO)₁₀(μ-H)(μ-SC₆H₄ -p-OH)] (2-OH), which showed comparable cytotoxicity. Mode of action studies pointed to an apoptotic pathway for cell death.

Artificial Intelligence in Chemistry: Current Trends and Future Directions

The application of artificial intelligence (AI) to chemistry has grown tremendously in recent years. In this Review, we studied the growth and distribution of AI-related chemistry publications in the last two decades using the CAS Content Collection. The volume of both journal and patent publications have increased dramatically, especially since 2015. Study of the distribution of publications over various chemistry research areas revealed that analytical chemistry and biochemistry are integrating AI to the greatest extent and with the highest growth rates. We also investigated trends in interdisciplinary research and identified frequently occurring combinations of research areas in publications. Furthermore, topic analyses were conducted for journal and patent publications to illustrate emerging associations of AI with certain chemistry research topics. Notable publications in various chemistry disciplines were then evaluated and presented to highlight emerging use cases. Finally, the occurrence of different classes of substances and their roles in AI-related chemistry research were quantified, further detailing the popularity of AI adoption in the life sciences and analytical chemistry. In summary, this Review offers a broad overview of how AI has progressed in various fields of chemistry and aims to provide an understanding of its future directions.

Structure-based virtual screening of CYP1A1 inhibitors: towards rapid tier-one assessment of potential developmental toxicants

Cytochrome P450 1A1 (CYP1A1) metabolizes estrogens, melatonin, and other key endogenous signaling molecules critical for embryonic/fetal development. The enzyme has increasing expression during pregnancy, and its inhibition or knockout increases embryonic/fetal lethality and/or developmental problems. Here, we present a virtual screening model for CYP1A1 inhibitors based on the orthosteric and predicted allosteric sites of the enzyme. Using 1001 reference compounds with CYP1A1 activity data, we optimized the decision thresholds of our model and classified the training compounds with 68.3% balanced accuracy (91.0% sensitivity and 45.7% specificity). We applied our final model to 11 known CYP1A1 orthosteric binders and related compounds, and found that our ranking of the known orthosteric binders generally agrees with the relative activity of CYP1A1 in metabolizing these compounds. We also applied the model to 22 new test compounds with unknown/unclear CYP1A1 inhibitory activity, and predicted 16 of them are CYP1A1 inhibitors. The CYP1A1 potency and modes of inhibition of these 22 compounds were experimentally determined. We confirmed that most predicted inhibitors, including drugs contraindicated during pregnancy (amiodarone, bicalutamide, cyproterone acetate, ketoconazole, and tamoxifen) and environmental agents suspected to be endocrine disruptors (bisphenol A, diethyl and dibutyl phthalates, and zearalenone), are indeed potent inhibitors of CYP1A1. Our results suggest that virtual screening may be used as a rapid tier-one method to screen for potential CYP1A1 inhibitors, and flag them out for further experimental evaluations.

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.

Image-based profiling for drug discovery: due for a machine-learning upgrade?

Image-based profiling is a maturing strategy by which the rich information present in biological images is reduced to a multidimensional profile, a collection of extracted image-based features. These profiles can be mined for relevant patterns, revealing unexpected biological activity that is useful for many steps in the drug discovery process. Such applications include identifying disease-associated screenable phenotypes, understanding disease mechanisms and predicting a drug's activity, toxicity or mechanism of action. Several of these applications have been recently validated and have moved into production mode within academia and the pharmaceutical industry. Some of these have yielded disappointing results in practice but are now of renewed interest due to improved machine-learning strategies that better leverage image-based information. Although challenges remain, novel computational technologies such as deep learning and single-cell methods that better capture the biological information in images hold promise for accelerating drug discovery.

Virtual screening of potentially endocrine-disrupting chemicals against nuclear receptors and its application to identify PPARγ-bound fatty acids

Nuclear receptors (NRs) are key regulators of energy homeostasis, body development, and sexual reproduction. Xenobiotics binding to NRs may disrupt natural hormonal systems and induce undesired adverse effects in the body. However, many chemicals of concerns have limited or no experimental data on their potential or lack-of-potential endocrine-disrupting effects. Here, we propose a virtual screening method based on molecular docking for predicting potential endocrine-disrupting chemicals (EDCs) that bind to NRs. For 12 NRs, we systematically analyzed how multiple crystal structures can be used to distinguish actives and inactives found in previous high-throughput experiments. Our method is based on (i) consensus docking scores from multiple structures at a single functional state (agonist-bound or antagonist-bound), (ii) multiple functional states (agonist-bound and antagonist-bound), and (iii) multiple pockets (orthosteric site and alternative sites) of these NRs. We found that the consensus enrichment from multiple structures is better than or comparable to the best enrichment from a single structure. The discriminating power of this consensus strategy was further enhanced by a chemical similarity-weighted scoring scheme, yielding better or comparable enrichment for all studied NRs. Applying this optimized method, we screened 252 fatty acids against peroxisome proliferator-activated receptor gamma (PPARγ) and successfully identified 3 previously unknown fatty acids with Kd = 100-250 μM including two furan fatty acids: furannonanoic acid (FNA) and furanundecanoic acid (FUA), and one cyclopropane fatty acid: phytomonic acid (PTA). These results suggested that the proposed method can be used to rapidly screen and prioritize potential EDCs for further experimental evaluations.

Optimum concentration-response curve metrics for supervised selection of discriminative cellular phenotypic endpoints for chemical hazard assessment

High-content imaging (HCI) provides quantitative and information-rich measurements of chemical effects on human in vitro cell models. Identification of discriminative phenotypic endpoints from cellular features obtained from HCI is required for accurate assessments of potential chemical hazards. However, the use of suboptimal metrics to quantify the concentration-response curves (CRC) of chemicals based on these features may obscure discriminative features, and lead to non-predictive endpoints and poor chemical classifications or hazard assessments. Here, we present a systematic and data-driven study on the performances of different CRC metrics in identifying image-based phenotypic features that can accurately classify the effects of reference chemicals with known in vivo toxicities. We studied four previous HCI in vitro nephro- or pulmono-toxicity datasets, which contain phenotypic feature measurements from different cell and feature types. Within a feature type, we found that efficacy metrics at higher chemical concentrations tend to give higher classification accuracy, whereas potency metrics do not have obvious trends across different response levels. Across different cell and feature types, efficacy metrics generally gave higher classification accuracy than potency metrics and area under the curve (AUC). Our results suggest that efficacy metrics, especially at higher concentrations, are more likely to help us to identify discriminative phenotypic endpoints. Therefore, HCI experiments for toxicological applications should include measurements at sufficiently high chemical concentrations, and efficacy metrics should always be analyzed. The identified features may be used as specific toxicity endpoints for further chemical hazard assessment.

Predicting direct hepatocyte toxicity in humans by combining high-throughput imaging of HepaRG cells and machine learning-based phenotypic profiling

Accurate prediction of drug- and chemical-induced hepatotoxicity remains to be a problem for pharmaceutical companies as well as other industries and regulators. The goal of the current study was to develop an in vitro/in silico method for the rapid and accurate prediction of drug- and chemical-induced hepatocyte injury in humans. HepaRG cells were employed for high-throughput imaging in combination with phenotypic profiling. A reference set of 69 drugs and chemicals was screened at a range of 7 concentrations, and the cellular response values were used for training a supervised classifier and for determining assay performance by using tenfold cross-validation. The results showed that the best performing phenotypic features were related to nuclear translocation of RELA (RELA proto-oncogene, NF-kB subunit; also known as NF-kappa B p65), DNA organization, and the F-actin cytoskeleton. Using a subset of 30 phenotypic features, direct hepatocyte toxicity in humans could be predicted with a test sensitivity, specificity and balanced accuracy of 73%, 92%, and 83%, respectively. The method was applied to another set of 26 drugs and chemicals with unclear annotation and their hepatocyte toxicity in humans was predicted. The results also revealed that the identified discriminative phenotypic changes were related to cell death and cellular senescence. Whereas cell death-related endpoints are widely applied in in vitro toxicology, cellular senescence-related endpoints are not, although cellular senescence can be induced by various drugs and other small molecule compounds and plays an important role in liver injury and disease. These findings show how phenotypic profiling can reveal unexpected chemical-induced mechanisms in toxicology.

Rapid hazard characterization of environmental chemicals using a compendium of human cell lines from different organs

The lack of adequate toxicity data for the vast majority of chemicals in the environment has spurred the development of new approach methodologies (NAMs). This study aimed to develop a practical high-throughput in vitro model for rapidly evaluating potential hazards of chemicals using a small number of human cells. Forty-two compounds were tested using human induced pluripotent stem cell (iPSC)-derived cells (hepatocytes, neurons, cardiomyocytes and endothelial cells), and a primary endothelial cell line. Both functional and cytotoxicity endpoints were evaluated using high-content imaging. Concentration-response was used to derive points-of-departure (POD). PODs were integrated with ToxPi and used as surrogate NAM-based PODs for risk characterization. We found chemical class-specific similarity among the chemicals tested; metal salts exhibited the highest overall bioactivity. We also observed cell type-specific patterns among classes of chemicals, indicating the ability of the proposed in vitro model to recognize effects on different cell types. Compared to available NAM datasets, such as ToxCast/Tox21 and chemical structure-based descriptors, we found that the data from the five-cell-type model was as good or even better in assigning compounds to chemical classes. Additionally, the PODs from this model performed well as a conservative surrogate for regulatory in vivo PODs and were less likely to underestimate in vivo potency and potential risk compared to other NAM-based PODs. In summary, we demonstrate the potential of this in vitro screening model to inform rapid risk-based decision-making through ranking, clustering, and assessment of both hazard and risks of diverse environmental chemicals.

Assessing Toxicity with Human Cell-Based In Vitro Methods

In toxicology, there is a strong push towards replacing animal experiments with alternative methods, which include cell-based in vitro methods for the assessment of adverse health effects in humans. High-throughput methods are of central interest due to the large and steadily growing numbers of compounds that require assessment. Tremendous progress has been made during the last decade in developing and applying such methods. Innovative technologies for addressing complex biological interactions include induced pluripotent stem cell- and organoid-based approaches, organotypic coculture systems, and microfluidic 'multiorgan' chips. Combining in vitro methods with bioinformatics and in silico modeling generates new powerful tools for toxicity assessment, and the rapid progress in the field is expected to continue.

A Case Study with Triazole Fungicides to Explore Practical Application of Next-Generation Hazard Assessment Methods for Human Health

The ongoing developments in chemical risk assessment have led to new concepts building on integration of sophisticated nonanimal models for hazard characterization. Here we explore a pragmatic approach for implementing such concepts, using a case study of three triazole fungicides, namely, flusilazole, propiconazole, and cyproconazole. The strategy applied starts with evaluating the overall level of concern by comparing exposure estimates to toxicological potential, followed by a combination of in silico tools and literature-derived high-throughput screening assays and computational elaborations to obtain insight into potential toxicological mechanisms and targets in the organism. Additionally, some targeted in vitro tests were evaluated for their utility to confirm suspected mechanisms of toxicity and to generate points of departure. Toxicological mechanisms instead of the current "end point-by-end point" approach should guide the selection of methods and assays that constitute a toolbox for next-generation risk assessment. Comparison of the obtained in silico and in vitro results with data from traditional in vivo testing revealed that, overall, nonanimal methods for hazard identification can produce adequate qualitative hazard information for risk assessment. Follow-up studies are needed to further refine the proposed approach, including the composition of the toolbox, toxicokinetics models, and models for exposure assessment.