Search

Abstract The rapid increase of publicly available chemical structures and associated experimental data presents a valuable opportunity to build robust QSAR models for applications in different fields. However, the common concern is the quality of both the chemical structure information and associated experimental data. This is especially true when those data are collected from multiple sources as chemical substance mappings can contain many duplicate structures and molecular inconsistencies. Such issues can impact the resulting molecular descriptors and their mappings to experimental data and, subsequently, the quality of the derived models in terms of accuracy, repeatability, and reliability. Herein we describe the development of an automated workflow to standardize chemical structures according to a set of standard rules and generate two and/or three-dimensional “QSAR-ready” forms prior to the calculation of molecular descriptors. The workflow was designed in the KNIME workflow environment and consists of three high-level steps. First, a structure encoding is read, and then the resulting in-memory representation is cross-referenced with any existing identifiers for consistency. Finally, the structure is standardized using a series of operations including desalting, stripping of stereochemistry (for two-dimensional structures), standardization of tautomers and nitro groups, valence correction, neutralization when possible, and then removal of duplicates. This workflow was initially developed to support collaborative modeling QSAR projects to ensure consistency of the results from the different participants. It was then updated and generalized for other modeling applications. This included modification of the “QSAR-ready” workflow to generate “MS-ready structures” to support the generation of substance mappings and searches for software applications related to non-targeted analysis mass spectrometry. Both QSAR and MS-ready workflows are freely available in KNIME, via standalone versions on GitHub, and as docker container resources for the scientific community. Scientific contribution: This work pioneers an automated workflow in KNIME, systematically standardizing chemical structures to ensure their readiness for QSAR modeling and broader scientific applications. By addressing data quality concerns through desalting, stereochemistry stripping, and normalization, it optimizes molecular descriptors' accuracy and reliability. The freely available resources in KNIME, GitHub, and docker containers democratize access, benefiting collaborative research and advancing diverse modeling endeavors in chemistry and mass spectrometry.

Organophosphorus flame retardants (OPFRs) are abundant and persistent in the environment but have limited toxicity information. Their similarity in structure to organophosphate pesticides presents great concern for developmental neurotoxicity (DNT). However, current in vivo testing is not suitable to provide DNT information on the amount of OPFRs that lack data. Over the past decade, an in vitro battery was developed to enhance DNT assessment, consisting of assays that evaluate cellular processes in neurodevelopment and function. In this study, behavioral data of small model organisms were also included. To assess if these assays provide sufficient mechanistic coverage to prioritize chemicals for further testing and/or identify hazards, an integrated approach to testing and assessment (IATA) was developed with additional information from the Integrated Chemical Environment (ICE) and the literature. Human biomonitoring and exposure data were identified and physiologically-based toxicokinetic models were applied to relate in vitro toxicity data to human exposure based on maximum plasma concentration. Eight OPFRs were evaluated, including aromatic OPFRs (triphenyl phosphate (TPHP), isopropylated phenyl phosphate (IPP), 2-ethylhexyl diphenyl phosphate (EHDP), tricresyl phosphate (TMPP), isodecyl diphenyl phosphate (IDDP), tert-butylphenyl diphenyl phosphate (BPDP)) and halogenated FRs ((Tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tris(2-chloroethyl) phosphate (TCEP)). Two representative brominated flame retardants (BFRs) (2,2′4,4′-tetrabromodiphenyl ether (BDE-47) and 3,3′,5,5′-tetrabromobisphenol A (TBBPA)) with known DNT potential were selected for toxicity benchmarking. Data from the DNT battery indicate that the aromatic OPFRs have activity at similar concentrations as the BFRs and should therefore be evaluated further. However, these assays provide limited information on the mechanism of the compounds. By integrating information from ICE and the literature, endocrine disruption was identified as a potential mechanism. This IATA case study indicates that human exposure to some OPFRs could lead to a plasma concentration similar to those exerting in vitro activities, indicating potential concern for human health.

Introduction: Skin sensitization, which leads to allergic contact dermatitis, is a key toxicological endpoint with high occupational and consumer prevalence. This study optimized several in vitro assays listed in OECD skin sensitization test guidelines for use on a quantitative high-throughput screening (qHTS) platform and performed in silico model predictions to assess the skin sensitization potential of prioritized compounds from the Tox21 10K compound library.Methods: First, we screened the entire Tox21 10K compound library using a qHTS KeratinoSensTM (KS) assay and built a quantitative structure–activity relationship (QSAR) model based on the KS results. From the qHTS KS screening results, we prioritized 288 compounds to cover a wide range of structural chemotypes and tested them in the solid phase extraction–tandem mass spectrometry (SPE–MS/MS) direct peptide reactivity assay (DPRA), IL-8 homogeneous time-resolved fluorescence (HTRF) assay, CD86 and CD54 surface expression in THP1 cells, and predicted in silico sensitization potential using the OECD QSAR Toolbox (v4.5).Results: Interpreting tiered qHTS datasets using a defined approach showed the effectiveness and efficiency of in vitro methods. We selected structural chemotypes to present this diverse chemical collection and to explore previously unidentified structural contributions to sensitization potential.Discussion: Here, we provide a skin sensitization dataset of unprecedented size, along with associated tools, and analysis designed to support chemical assessments.

All-trans retinoic acid (ATRA) gradients determine skeletal patterning morphogenesis and can be disrupted by diverse genetic or environmental factors during pregnancy, leading to fetal skeleton defects. Adverse Outcome Pathway (AOP) frameworks for ATRA metabolism, signaling, and homeostasis allow for the development of new approach methods (NAMs) for predictive toxicology with less reliance on animal testing. Here, a data-driven model was constructed to identify chemicals associated with both ATRA pathway bioactivity and prenatal skeletal defects. The phenotype data was culled from ToxRefDB prenatal developmental toxicity studies and produced a list of 363 ToxRefDB chemicals with altered skeletal observations. Defects were classified regionally as cranial, post-cranial axial, appendicular, and other (unspecified) features based on ToxRefDB descriptors. To build a multivariate statistical model, high-throughput screening bioactivity data from >8,070 chemicals in ToxCast/Tox21 across 10 in vitro assays relevant to the retinoid signaling system were evaluated and compared to literature-based candidate reference chemicals in the dataset. There were 48 chemicals identified for effects on both in vivo skeletal defects and in vitro ATRA pathway targets for computational modeling. The list included 28 chemicals with prior evidence of skeletal defects linked to retinoid toxicity and 20 chemicals without prior evidence. The combination of thoracic cage defects and DR5 (direct repeats of 5 nucleotides for RAR/RXR transactivation) disruption was the most frequently occurring phenotypic and target disturbance, respectively. This data model provides valuable AOP elucidation and validates current mechanistic understanding. These findings also shed light on potential avenues for new mechanistic discoveries related to ATRA pathway disruption and associated skeletal dysmorphogenesis due to environmental exposures.

Abstract Background The logarithmic acid dissociation constant pKa reflects the ionization of a chemical, which affects lipophilicity, solubility, protein binding, and ability to pass through the plasma membrane. Thus, pKa affects chemical absorption, distribution, metabolism, excretion, and toxicity properties. Multiple proprietary software packages exist for the prediction of pKa, but to the best of our knowledge no free and open-source programs exist for this purpose. Using a freely available data set and three machine learning approaches, we developed open-source models for pKa prediction. Methods The experimental strongest acidic and strongest basic pKa values in water for 7912 chemicals were obtained from DataWarrior, a freely available software package. Chemical structures were curated and standardized for quantitative structure–activity relationship (QSAR) modeling using KNIME, and a subset comprising 79% of the initial set was used for modeling. To evaluate different approaches to modeling, several datasets were constructed based on different processing of chemical structures with acidic and/or basic pKas. Continuous molecular descriptors, binary fingerprints, and fragment counts were generated using PaDEL, and pKa prediction models were created using three machine learning methods, (1) support vector machines (SVM) combined with k-nearest neighbors (kNN), (2) extreme gradient boosting (XGB) and (3) deep neural networks (DNN). Results The three methods delivered comparable performances on the training and test sets with a root-mean-squared error (RMSE) around 1.5 and a coefficient of determination (R2) around 0.80. Two commercial pKa predictors from ACD/Labs and ChemAxon were used to benchmark the three best models developed in this work, and performance of our models compared favorably to the commercial products. Conclusions This work provides multiple QSAR models to predict the strongest acidic and strongest basic pKas of chemicals, built using publicly available data, and provided as free and open-source software on GitHub.

Abstract The median lethal dose for rodent oral acute toxicity (LD50) is a standard piece of information required to categorize chemicals in terms of the potential hazard posed to human health after acute exposure. The exclusive use of in vivo testing is limited by the time and costs required for performing experiments and by the need to sacrifice a number of animals. (Quantitative) structure–activity relationships [(Q)SAR] proved a valid alternative to reduce and assist in vivo assays for assessing acute toxicological hazard. In the framework of a new international collaborative project, the NTP Interagency Center for the Evaluation of Alternative Toxicological Methods and the U.S. Environmental Protection Agency’s National Center for Computational Toxicology compiled a large database of rat acute oral LD50 data, with the aim of supporting the development of new computational models for predicting five regulatory relevant acute toxicity endpoints. In this article, a series of regression and classification computational models were developed by employing different statistical and knowledge-based methodologies. External validation was performed to demonstrate the real-life predictability of models. Integrated modeling was then applied to improve performance of single models. Statistical results confirmed the relevance of developed models in regulatory frameworks, and confirmed the effectiveness of integrated modeling. The best integrated strategies reached RMSEs lower than 0.50 and the best classification models reached balanced accuracies over 0.70 for multi-class and over 0.80 for binary endpoints. Computed predictions will be hosted on the EPA’s Chemistry Dashboard and made freely available to the scientific community.

During the past few decades, the science of toxicology has been undergoing a transformation from observational to predictive science. New approach methodologies (NAMs), including in vitro assays, in silico models, read-across, and in vitro to in vivo extrapolation (IVIVE), are being developed to reduce, refine, or replace whole animal testing, encouraging the judicious use of time and resources. Some of these methods have advanced past the exploratory research stage and are beginning to gain acceptance for the risk assessment of chemicals. A review of the recent literature reveals a burst of IVIVE publications over the past decade. In this review, we propose operational definitions for IVIVE, present literature examples for several common toxicity endpoints, and highlight their implications in decision-making processes across various federal agencies, as well as international organizations, including those in the European Union (EU). The current challenges and future needs are also summarized for IVIVE. In addition to refining and reducing the number of animals in traditional toxicity testing protocols and being used for prioritizing chemical testing, the goal to use IVIVE to facilitate the replacement of animal models can be achieved through their continued evolution and development, including a strategic plan to qualify IVIVE methods for regulatory acceptance.

Assessing skin irritation potential is critical for the safety evaluation of topical drugs and other consumer products such as cosmetics. The use of advanced cellular models, as an alternative to replace animal testing in the safety evaluation for both consumer products and ingredients, is already mandated by law in the European Union (EU) and other countries. However, there has not yet been a large-scale comparison of the effects of topical-use compounds in different cellular skin models. This study assesses the irritation potential of topical-use compounds in different cellular models of the skin that are compatible with high throughput screening (HTS) platforms. A set of 451 topical-use compounds were first tested for cytotoxic effects using two-dimensional (2D) monolayer models of primary neonatal keratinocytes and immortalized human keratinocytes. Forty-six toxic compounds identified from the initial screen with the monolayer culture systems were further tested for skin irritation potential on reconstructed human epidermis (RhE) and full thickness skin (FTS) three-dimensional (3D) tissue model constructs. Skin irritation potential of the compounds was assessed by measuring tissue viability, trans-epithelial electrical resistance (TEER), and secretion of cytokines interleukin 1 alpha (IL-1α) and interleukin 18 (IL-18). Among known irritants, high concentrations of methyl violet and methylrosaniline decreased viability, lowered TEER, and increased IL-1α secretion in both RhE and FTS models, consistent with irritant properties. However, at low concentrations, these two compounds increased IL-18 secretion without affecting levels of secreted IL-1α, and did not reduce tissue viability and TEER, in either RhE or FTS models. This result suggests that at low concentrations, methyl violet and methylrosaniline have an allergic potential without causing irritation. Using both HTS-compatible 2D cellular and 3D tissue skin models, together with irritation relevant activity endpoints, we obtained data to help assess the irritation effects of topical-use compounds and identify potential dermal hazards.

The vast landscape of environmental chemicals has motivated the need for alternative methods to traditional whole-animal bioassays in toxicity testing. Embryonic stem (ES) cells provide an in vitro model of embryonic development and an alternative method for assessing developmental toxicity. Here, we evaluated 309 environmental chemicals, mostly food-use pesticides, from the ToxCast™ chemical library using a mouse ES cell platform. ES cells were cultured in the absence of pluripotency factors to promote spontaneous differentiation and in the presence of DMSO-solubilized chemicals at different concentrations to test the effects of exposure on differentiation and cytotoxicity. Cardiomyocyte differentiation (α,β myosin heavy chain; MYH6/MYH7) and cytotoxicity (DRAQ5™/Sapphire700™) were measured by In-Cell Western™ analysis. Half-maximal activity concentration (AC₅₀) values for differentiation and cytotoxicity endpoints were determined, with 18% of the chemical library showing significant activity on either endpoint. Mining these effects against the ToxCast Phase I assays (∼500) revealed significant associations for a subset of chemicals (26) that perturbed transcription-based activities and impaired ES cell differentiation. Increased transcriptional activity of several critical developmental genes including BMPR2, PAX6 and OCT1 were strongly associated with decreased ES cell differentiation. Multiple genes involved in reactive oxygen species signaling pathways (NRF2, ABCG2, GSTA2, HIF1A) were strongly associated with decreased ES cell differentiation as well. A multivariate model built from these data revealed alterations in ABCG2 transporter was a strong predictor of impaired ES cell differentiation. Taken together, these results provide an initial characterization of metabolic and regulatory pathways by which some environmental chemicals may act to disrupt ES cell growth and differentiation.