A review of the electrophilic reaction chemistry involved in covalent protein binding relevant to toxicity

Clustering a Chemical Inventory for Safety Assessment of Fragrance Ingredients: Identifying Read-Across Analogs to Address Data Gaps

A valuable approach to chemical safety assessment is the use of read-across chemicals to provide safety data to support the assessment of structurally similar chemicals. An inventory of over 6000 discrete organic chemicals used as fragrance materials in consumer products has been clustered into chemical class-based groups for efficient search of read-across sources. We developed a robust, tiered system for chemical classification based on (1) organic functional group, (2) structural similarity and reactivity features of the hydrocarbon skeletons, (3) predicted or experimentally verified Phase I and Phase II metabolism, and (4) expert pruning to consider these variables in the context of specific toxicity end points. The systematic combination of these data yielded clusters, which may be visualized as a top-down hierarchical clustering tree. In this tree, chemical classes are formed at the highest level according to organic functional groups. Each subsequent subcluster stemming from classes in this hierarchy of the cluster is a chemical cluster defined by common organic functional groups and close similarity in the hydrocarbon skeleton. By examining the available experimental data for a toxicological endpoint within each cluster, users can better identify potential read-across chemicals to support safety assessments.

Quantification of PET cyclic and linear oligomers in teabags by a validated LC-MS method - In silico toxicity assessment and consumer's exposure

Determination of polyethylene terephthalate (PET) dimer up to heptamer 1st series cyclic oligomers, applying an LC-qTOF-MS method, has been developed and validated. Recoveries ranged between 80 and 112% with RSDs lower than 15%. An innovative semi-quantitative approach has been applied for 2nd and 3rd series cyclic oligomers, using the closest structural-similar 1st series cyclic oligomer standard as analytical reference. Oligomers from the three series were quantified in PET teabags after migration experiments with water and food simulants C (20% v/v ethanol in water) and D1 (50% v/v ethanol in water). No legal migration limits exist currently for these substances. In silico genotoxicity assessment of all identified oligomers has been performed and showed no genotoxicity alert for linear or cyclic molecules. Exposure assessment was performed using EFSA's approach on the total sum of migrating oligomers and on toxicological threshold-of-concern. Amounts found in water were in some cases significantly higher than the respective limits, especially in the worst-case scenario of multiple consumption.

Evaluation of two in silico programs for predicting mutagenicity and carcinogenicity potential for 4-methylimidazole (4-MeI) and known metabolites

4-Methylimidazole (4-MeI) is a nitrogen-containing heterocyclic compound that is used in the manufacture of chemicals, dyes and pharmaceuticals and may be found in a variety of foods following formation during heating. The purpose of this study was to use two different in silico programs, CASE Ultra and Toxtree, to investigate potential structure-activity relationships in 4-MeI and its metabolites for mutagenicity and carcinogenicity, and combine that information with the available literature to draw conclusions regarding the strength of the predictions observed. Neither CASE Ultra nor Toxtree identified any structural alerts that were associated with mutagenic activity. Data for 4-MeI from a single study were used in the development of the CASE Ultra mouse and rat carcinogenicity models, but no additional similar structures were identified in the carcinogenicity model training set. One metabolite, 5-methylhydantoin, was predicted to be positive in the CASE Ultra carcinogenicity male and female mouse models; positive predictivity percentages of 60.9% and 73.7%, respectively. However, low structural similarity between 5-methylhydantoin and the compounds identified in the training set (<25%) decreases confidence in the positive prediction. Three metabolites were predicted to be positive in the CASE Ultra mouse micronucleus model, but again suffered from low structural similarity. Both limited structural similarity and inconsistent responses among the other clastogenicity models suggest that additional structurally similar compounds are needed to assess the predictive capacity of these alerts for biological activity of these compounds.

Toward Less Hazardous Industrial Compounds: Coupling Quantum Mechanical Computations, Biomarker Responses, and Behavioral Profiles To Identify Bioactivity of SN2 Electrophiles in Alternative Vertebrate Models

Sustainable molecular design of less hazardous chemicals promises to reduce risks to public health and the environment. Computational chemistry modeling coupled with alternative toxicology models (e.g., larval fish) present unique high-throughput opportunities to understand structural characteristics eliciting adverse outcomes. Numerous environmental contaminants with reactive properties can elicit oxidative stress, an important toxicological response associated with diverse adverse outcomes (i.e., cancer, diabetes, neurodegenerative disorders, etc.). We examined a common chemical mechanism (bimolecular nucleophilic substitution (S_N2)) associated with oxidative stress using property-based computational modeling coupled with acute (mortality) and sublethal (glutathione, photomotor behavior) responses in the zebrafish (Danio rerio) and the fathead minnow (Pimephales promelas) models to identify whether relationships exist among biological responses and molecular attributes of industrial chemicals. Following standardized methods, embryonic zebrafish and larval fathead minnows were exposed separately to eight different S_N2 compounds for 96 h. Acute and sublethal responses were compared to computationally derived in silico chemical descriptors. Specifically, frontier molecular orbital energies were significantly related to acute LC₅₀ values and photomotor response (PMR) no observed effect concentrations (NOECs) in both fathead minnow and zebrafish. This reactivity index, LC₅₀ values, and PMR NOECs were also significantly related to whole body glutathione (GSH) levels, suggesting that acute and chronic toxicity results from protein adduct formation for S_N2 electrophiles. Shared refractory locomotor response patterns among study compounds and two alternative vertebrate models appear informative of electrophilic properties associated with oxidative stress for S_N2 chemicals. Electrophilic parameters derived from frontier molecular orbitals were predictive of experimental in vivo acute and sublethal toxicity. These observations provide important implications for identifying and designing less hazardous industrial chemicals with reduced potential to elicit oxidative stress through bimolecular nucleophilic substitution.

Comparing and contrasting the coverage of publicly available structural alerts for protein binding

The molecular initiating event for many mechanisms of toxicological action comprise the reactive, covalent binding between an exogenous electrophile and an endogenous nucleophile. The target sites for electrophiles are typically peptides, proteins, enzymes or DNA. Of these, the formation of covalent adducts with proteins and DNA are perhaps the most established as they are most closely associated with skin sensitisation and genotoxicity endpoints. As such, being able to identify electrophilic features within a chemical structure provides a starting point to characterise its reactivity profile. There are a number of software tools that have been developed to help identify structural features indicative of electrophilic reactive potential to address various purposes, including: 1) to facilitate category formation for read-across of toxicity effects such as skin sensitisation potential, as well as 2) to profile substances to identify potential confounding factors to rationalise their activity in high-throughput screening (HTS) assays. Here, three such schemes that have been published in the literature as collections of SMARTS patterns and their associated chemical-biological reaction domains have been compared. The goals are 1) to better understand their scope and coverage, and 2) to assess their performance relative to a published skin sensitisation dataset where manual annotations to assign likely mechanistic domains based on expert judgement were already available. The 3 schemes were then applied to the Tox21 library and the consensus outcome was reported to highlight the proportion of chemicals likely to exhibit a reactivity response, specific to a mechanistic reaction domain, but non-specific with respect to target-tissue based activity. ToxPrint fingerprints were computed and activity enrichments computed to compare the structural features identified for the skin sensitisation dataset and Tox21 chemicals for each 'consensus' reaction domain. Enriched ToxPrints were also used to identify ToxCast assays potentially informative for reactivity.

Synergistic degradation on phenolic compounds of coal pyrolysis wastewater (CPW) by lignite activated coke-active sludge (LAC-AS) process: Insights into succession of microbial community under selective pressure

This study illustrated synergistic degradation of phenolic compounds by LAC-AS process via the insight into succession of microbial community under selective pressure. The results demonstrated that high phenols exhibited toxicity pressure to single AS process by eliminating non-tolerate bacteria, inducing vicious circulation by intermediates (catechol, nitrate, etc.) accumulation. However, LAC exerted another selective pressure and facilitated positive bio-community succession of moving biological bed reactor (MBBR). Firstly, it created rich microenvironments for diverse bacteria and promoted resilient adsorption for phenols with the assistance of biodegradation. Secondly, LAC enriched facultative bacteria, which developed multiple degradation paths on phenols and nitrogen based on multifunctional genes, counteracting the toxicity pressure. Specifically, phenols were degraded by the combination of anaerobic hydrolysis and oxidation, while conventional and shortcut nitrification-denitrification (SND) and nitrogen fixation all participated in nitrogen removal, achieving high removal of COD (93.49%), Tph (93.74%), TN (92.20%) and NH₄⁺-N (93.20%) under the highest phenols.

Theoretical consideration on the prediction of in vivo toxicity from in vitro toxicity: Effect of bio-uptake equilibrium, kinetics and mode of action

Although in vitro assay is an ideal alternative method for the in vivo toxicity prediction, different in vivo-in vitro correlations have been observed for the toxicity endpoints obtained from different levels of species. In this paper, theoretical in vivo-in vitro toxicity correlations have been developed for cytotoxicity versus human, mammalian and fish toxicity, respectively. These theoretical models were then used to investigate the correlations and the influencing factors between in vivo and in vitro toxicity. Bio-uptake equilibrium theory can well explain why there is a significant correlation between fish and cell toxicity (R² = 0.70); why human toxicity is very close to fish toxicity; and why hydrophobic compounds exhibit relatively greater toxicity than reactive or specifically-acting compounds to human and fish as compared to cells. The kinetic theory can well explain why there is a very poor relationship between mammal and cell toxicity (R² = 0.44). This paper reveals that polar and ionized compounds can more easily pass through cell membrane and have greater bioconcentration potential. Increasing of hydrophobicity and ionization can increase the cytotoxicity. Inclusion of descriptors representing hydrophobicity, ionization, acidity and absorption into the correlation equations can significantly improve the correlations of cytotoxicity with human and fish toxicity (R² > 0.8), but not with mammal toxicity (R² = 0.49). These descriptors reflect the differences of the toxicodynamics and toxicokinetics between cells and organisms.

5-Arylidene(chromenyl-methylene)-thiazolidinediones: Potential New Agents against Mutant Oncoproteins K-Ras, N-Ras and B-Raf in Colorectal Cancer and Melanoma

Background and objectives: Cancer represents the miscommunication between and within the body cells. The mutations of the oncogenes encoding the MAPK pathways play an important role in the development of tumoral diseases. The mutations of KRAS and BRAF oncogenes are involved in colorectal cancer and melanoma, while the NRAS mutations are associated with melanoma. Thiazolidine-2,4-dione is a versatile scaffold in medicinal chemistry and a useful tool in the development of new antitumoral compounds. The aim of our study was to predict the pharmacokinetic/pharmacodynamic properties, the drug-likeness and lead-likeness of two series of synthetic 5-arylidene(chromenyl-methylene)-thiazolidinediones, the molecular docking on the oncoproteins K-Ras, N-Ras and B-Raf, and to investigate the cytotoxicity of the compounds, in order to select the best structural profile for potential anticancer agents. Materials and Methods: In our paper we studied the cytotoxicity of two series of thiazolidine-2,4-dione derivatives, their ADME-Tox properties and the molecular docking on a mutant protein of K-Ras, two isoforms of N-Ras and an isoform of B-Raf with 16 mutations. Results: The heterocyclic compounds strongly interact with K-Ras and N-Ras right after their posttranslational processing and/or compete with GDP for the nucleotide-binding site of the two GTPases. They are less active against the GDP-bound states of the two targets. All derivatives have a similar binding pattern in the active site of B-Raf. Conclusions: The data obtained encourage the further investigation of the 5-arylidene(chromenyl-methylene)-thiazolidinediones as potential new agents against the oncoproteins K-Ras, N-Ras and B-Raf.

Protein Adductomics: Methodologies for Untargeted Screening of Adducts to Serum Albumin and Hemoglobin in Human Blood Samples

The reaction products of electrophiles in vivo can be measured as adducts to the abundant proteins, hemoglobin (Hb), and human serum albumin (HSA), in human blood samples. During the last decade, methods for untargeted screening of such adducts, called "adductomics", have used liquid chromatography-mass spectrometry to detect large numbers of previously unknown Hb and HSA adducts. This review presents methodologies that were developed and used in our laboratories for Hb and HSA adductomics, respectively. We discuss critical aspects regarding choice of target protein, sample preparation, mass spectrometry, data evaluation, and strategies for identification of detected unknown adducts. With this review we give an overview of these two methodologies used for protein adductomics and the precursor electrophiles that have been elucidated from the adducts.

Toxicity of some prevalent organic chemicals to tadpoles and comparison with toxicity to fish based on mode of toxic action

Although mode of action (MOA) plays a key role in the understanding of the toxic mechanism of chemicals, the MOAs of class-based compounds to tadpoles have not been investigated. To explore the MOAs, acute toxicity (expressed as log 1/LC₅₀) to Rana chensinensis tadpoles were determined and molecular descriptors were calculated. Quantitative structure-activity relationship (QSAR) showed that toxicity to tadpoles is closely related to the chemical octanol/water partition coefficient (log K_OW), energy of the lowest unoccupied molecular orbital (E_LUMO), and number of hydrogen bond donors and acceptors (NH_DA), representing the bio-uptake potential in tadpoles, the electrophilicity and hydrogen bonding capacity with target site(s), respectively. Comparison of the toxicity values between tadpoles and fish revealed that there were no significant differences for the overlapping compounds (average residual = 0.29 between tadpole and fish toxicity) with P values of interspecies correlation substantially less than 0.001. Classification of MOAs for the class-based compounds based on the excess toxicity calculated from toxicity ratio suggested that baseline, less inert compounds and some reactive or specifically-acting compounds share same MOAs between tadpoles and fish. Fish and tadpoles can serve as surrogates for each other in the safety evaluation of organic pollutants.

Structural alerts for the identification of bioaccumulative compounds

Legislators have included bioaccumulation in the evaluation of chemicals in the framework of the European Registration, Evaluation, Authorization, and Restriction of Chemicals (REACH) regulation. REACH requires information on the bioconcentration factor (BCF), which is a parameter for assessing bioaccumulation and encourages the use of a weight-of-evidence approach, including predictions from quantitative structure-activity relationships (QSARs). This study presents a novel approach, based on structural alerts, to be used as a decision-support system for the identification of substances with bioaccumulation potential. In a regulatory framework, these alerts can be integrated with other sources of information, such as experimental and in silico data, to reduce the uncertainty of the assessment, thereby supporting a weight-of-evidence approach. Moreover, the identified alerts have a direct connection with relevant structural features, thus fostering the applicability and interpretability of the approach. The structural alerts were identified on 779 chemicals annotated for their fish BCF, and the approach was then validated on 278 external molecules. The developed decision-support system allowed identification of 77% of bioaccumulative chemicals and was competitive with more complex QSAR models used in regulatory assessments. The approach is implemented in an easy-to-use workflow, provided free of charge. Integr Environ Assess Manag 2019;15:19-28. © 2018 SETAC.

Discrimination of typical cyclic compounds and selection of toxicity evaluation bioassays for coal gasification wastewater (CGW) based on toxicity mechanism of actions (MOAs)

This paper originally investigated toxicity discrimination of typical cyclic compounds and bioassays selection on toxicity evaluation for coal gasification wastewater (CGW) effluent with mechanism-oriented investigation. Initially, representative cyclic toxicants were selected and classified with quantitative structure-toxicity relationship (QSTR). Nitrogen heterocyclic compounds (NHCs) and polycyclic aromatic hydrocarbons (PAHs) were basically discriminated as nonpolar narcotics with significant correlation to hydrophobicity (p < 0.05, R² = 0.8668-0.9635), while phenols were regarded as polar narcotics and reactive compounds due to slight correlation to hydrophobicity (p > 0.05, R² < 0.5). Furthermore, specific mechanism of actions (MOAs) to various organisms revealed that phenols were discriminated as critical source of acute toxicity in CGW, with short-term visible and irreversible damage. However, NHCs and PAHs, which exerted accumulation toxicity rather than acute toxicity, might result in potential mutagenicity and unpredictable risk along the food chain. Afterwards, based on species sensitivity to typical toxicants and application in real CGW effluent, non-applicability of Chlorella vulgaris (C. vulgaris) was validated in toxicity evaluation. While Daphnia magna (D. magna) was suggested as a toxicity bioassay in entire effluent due to the highest sensitivity and applicability. Tetrahymena thermophile (T. pyriformis) might be applicable in effluent with low biodegradability due to similar evaluation results (TU = 8.90) to D. magna (TU = 6.67) in aerobic effluent. Finally, the relationship between toxicity and bioavailability based on typical pollutants and model species illustrated necessity for dualism toxicity-biodegradability investigation on CGW.

Comparison of modes of action between fish and zebrafish embryo toxicity for baseline, less inert, reactive and specifically-acting compounds

The mode of action (MOA) plays a key role in the risk assessment of pollutants in water. Although fish is a key model organism used in the risk assessment of pollutants in water, the MOAs have not been compared between fish and embryo toxicity for classified compounds. In this paper, regression analysis was carried out for fish and embryo toxicities against the calculated molecular descriptors and MOAs were evaluated from toxicity ratio. The toxicity significantly related with the chemical hydrophobicity for baseline and less inert compounds, respectively, indicates that these two classes of compounds share the same MOAs between fish and embryos. Comparison of the toxicity ratios shows that reactive compounds exhibit excess toxicity to both fish and embryos. These compounds can react covalently with biologically target molecules through nucleophilic addition reactions, Michael addition oxidation, or amination. Comparing with baseline, less inert and reactive compounds, many specifically-acting compounds have strong docking capacity with protein molecules. Some specifically-acting compounds, such as fungicides, have very similar toxic effect to both fish and embryos. However, insecticides are more toxic to fish than embryos; herbicides and medications are more toxic to embryos than fish. Differences in the interactions of chemicals with target molecules or bioconcentration potentials between fish and embryos may result in the differences in toxic effects. There are some factors that influence the identification of MOAs, such as quality of toxicity data, bioavailability and ionization. These factors should be considered in the identification of MOAs in the risk assessment of organic pollutants.

Discovery of Novel N-(4-Hydroxybenzyl)valine Hemoglobin Adducts in Human Blood

Humans are exposed to a wide range of electrophilic compounds present in our diet and environment or formed endogenously as part of normal physiological processes. These electrophiles can modify nucleophilic sites of proteins and DNA to form covalent adducts. Recently, powerful untargeted adductomic approaches have been developed for systematic screening of these adducts in human blood. Our earlier untargeted adductomics study detected 19 unknown adducts to N-terminal valine in hemoglobin (Hb) in human blood. We now describe a full characterization of one of these adducts, which corresponds to the addition of a 4-hydroxybenzyl (4-OHBn) group to N-terminal valine in Hb to form N(4-hydroxybenzyl)valine (4-OHBn-Val). The adduct structure was determined by comparison of its accurate mass, HPLC retention time, and MS/MS fragmentation to that of authentic standards prepared by chemical synthesis. Average 4-OHBn-Val adduct concentrations in 12 human blood samples were estimated to 380 ± 160 pmol/g Hb. Two possible routes of 4-OHBnVal adduct formation are proposed using two different precursor electrophiles: 4-quinone methide (4-QM) and 4-hydroxybenzaldehyde (4-OHBA). We found that 4-QM reacts rapidly with valine to form the 4-OHBn-Val adduct; however, the quinone methide is unstable under physiological conditions due to hydrolysis. It was shown that 4-OHBA forms reversible Schiff base adducts with valine, which can be stabilized via reduction in blood generating the 4-OHBn-Val adduct. In addition, trace amounts of isomeric 2-hydroxybenzyl-valine (2-OHBn-Val) adducts were detected in 12 human blood samples (estimated mean adduct level, 5.0 ± 1.4 pmol/g Hb). Further studies are needed to quantify the contributions from identified possible precursor electrophiles to the observed hydroxybenzyl adducts in humans.

A critical review of adverse effects to the kidney: mechanisms, data sources, and in silico tools to assist prediction

INTRODUCTION

The kidney is a major target for toxicity elicited by pharmaceuticals and environmental pollutants. Standard testing which often does not investigate underlying mechanisms has proven not to be an adequate hazard assessment approach. As such, there is an opportunity for the application of computational approaches that utilize multiscale data based on the Adverse Outcome Pathway (AOP) paradigm, coupled with an understanding of the chemistry underpinning the molecular initiating event (MIE) to provide a deep understanding of how structural fragments of molecules relate to specific mechanisms of nephrotoxicity. Aims covered: The aim of this investigation was to review the current scientific landscape related to computational methods, including mechanistic data, AOPs, publicly available knowledge bases and current in silico models, for the assessment of pharmaceuticals and other chemicals with regard to their potential to elicit nephrotoxicity. A list of over 250 nephrotoxicants enriched with, where possible, mechanistic and AOP-derived understanding was compiled. Expert opinion: Whilst little mechanistic evidence has been translated into AOPs, this review identified a number of data sources of in vitro, in vivo, and human data that may assist in the development of in silico models which in turn may shed light on the interrelationships between nephrotoxicity mechanisms.

Design, synthesis, structural characterization and in vitro evaluation of new 1,4-disubstituted-1,2,3-triazole derivatives against glioblastoma cells

A new series of 1,4-disubstituted-1,2,3-triazole derivatives were synthesized through the copper-catalyzed azide-alkyne 1,3-dipolar cycloaddition (Click chemistry) and their inhibitory activities were evaluated against different human glioblastoma (GBM) cell lines, including highly drug-resistant human cell lines GBM02, GBM95. The most effective compounds were 9d, containing the methylenoxy moiety linked to triazole and the tosyl-hydrazone group, and the symmetrical bis-triazole 10a, also containing methylenoxy moiety linked to triazole. Single crystal X-ray diffraction analysis was employed for structural elucidation of compound 9d. In silico analyses of physicochemical, pharmacokinetic, and toxicological properties suggest that compounds 8a, 8b, 8c, 9d, and 10a are potential candidates for central nervous system-acting drugs.

Toxicological assessment of epinephrine and norepinephrine by analog approach

Epinephrine and norepinephrine have been used in the management of anaphylactic reactions and cardiac resuscitation, along with treatment of asthma and glaucoma extensively, but their toxicological profiles are not yet completed. Based on this circumstance, various toxicological endpoints of epinephrine and norepinephrine were explored. Since there is a paucity of some endpoints' data, readacross was applied to fill the data gaps using analog approach. Along with structural similarity, biological and mechanistic plausibility were also considered in analog selection. The similarity justification and supporting experimental data were provided for uncertainty evaluation. Short term repeated dose toxicity values as NOAEL and LOAEL belonging to epinephrine were used to estimate the repeated dose toxicity of norepinephrine. The in vivo and in vitro mutagenicity tests were considered representative of genotoxicity. Both chemicals are showed to be non-genotoxic. They are experimentally reported to cause developmental and reproductive toxicity. For the carcinogenicity endpoint, a conclusion could not be reached because similar compounds were seen to show conflicting results.

Development of thresholds of excess toxicity for environmental species and their application to identification of modes of acute toxic action

The acute toxicity of organic pollutants to fish, Daphnia magna, Tetrahymena pyriformis, and Vibrio fischeri was investigated. The results indicated that the Toxicity Ratio (TR) threshold of log TR =1, which has been based on the distribution of toxicity data to fish, can also be used to discriminate reactive or specifically acting compounds from baseline narcotics for Daphnia magna and Vibrio fischeri. A log TR=0.84 is proposed for Tetrahymena pyriformis following investigation of the relationships between the species sensitivity and the absolute averaged residuals (AAR) between the predicted baseline toxicity and the experimental toxicity. Less inert compounds exhibit relatively higher toxicity to the lower species (Tetrahymena pyriformis and Vibrio fischeri) than the higher species (fish and Daphnia magna). A greater number of less inert compounds with log TR greater than the thresholds was observed for Tetrahymena pyriformis and Vibrio fischeri. This may be attributed to the hydrophilic compounds which may pass more easily through cell membranes than the skin or exoskeleton of organisms and have higher bioconcentration factors in the lower species, leading to higher toxicity. Most of classes of chemical associated with excess toxicity to one species also exhibited excess toxicity to other species, however, a few classes with excess toxicity to one species exhibiting narcotic toxicity to other species and thus may have different MOAs between species. Some ionizable compounds have log TR much lower than one because of the over-estimated log K_OW. The factors that influence the toxicity ratio calculated from baseline level are discussed in this paper.

Ab initio chemical safety assessment: A workflow based on exposure considerations and non-animal methods

We describe and illustrate a workflow for chemical safety assessment that completely avoids animal testing. The workflow, which was developed within the SEURAT-1 initiative, is designed to be applicable to cosmetic ingredients as well as to other types of chemicals, e.g. active ingredients in plant protection products, biocides or pharmaceuticals. The aim of this work was to develop a workflow to assess chemical safety without relying on any animal testing, but instead constructing a hypothesis based on existing data, in silico modelling, biokinetic considerations and then by targeted non-animal testing. For illustrative purposes, we consider a hypothetical new ingredient x as a new component in a body lotion formulation. The workflow is divided into tiers in which points of departure are established through in vitro testing and in silico prediction, as the basis for estimating a safe external dose in a repeated use scenario. The workflow includes a series of possible exit (decision) points, with increasing levels of confidence, based on the sequential application of the Threshold of Toxicological (TTC) approach, read-across, followed by an "ab initio" assessment, in which chemical safety is determined entirely by new in vitro testing and in vitro to in vivo extrapolation by means of mathematical modelling. We believe that this workflow could be applied as a tool to inform targeted and toxicologically relevant in vitro testing, where necessary, and to gain confidence in safety decision making without the need for animal testing.

Relationship between acute and chronic toxicity for prevalent organic pollutants in Vibrio fischeri based upon chemical mode of action

Chemicals show diverse modes of action (MOAs) in aquatic organisms depending upon acute and chronic toxicity evaluations. Here, toxicity data for Vibrio fischeri involving 52 compounds for acute and chronic toxicity were used to determine the congruence of acute and chronic toxicity for assessing MOAs. Using toxic ratios, most of the compounds categorized into MOAs that included baseline, less inert or reactive compounds with acute toxicity were also categorized as baseline, less inert or reactive compounds with chronic toxicity. However, significantly different toxic effects were observed with acute and chronic toxicity for the reactive and specific-acting compounds. The acute-chronic toxic ratios were smaller and less variable for the baseline and less inert compounds, but were greater and more variable for the reactive and specific-acting compounds. Baseline and less inert compounds share same MOAs, but reactive and specific-acting compounds have different MOAs between acute and chronic toxicity. Bioconcentration processes cannot reach an equilibrium for highly hydrophilic and ionized compounds with short-term exposure, resulting in lower toxicity compared to long-term exposure. Pronounced differences for the antibiotics were not only due to the difference in bioconcentration, but also due to a predicted difference in MOAs during acute and chronic exposures.

In Silico Prediction of Organ Level Toxicity: Linking Chemistry to Adverse Effects

In silico methods to predict toxicity include the use of (Quantitative) Structure-Activity Relationships ((Q)SARs) as well as grouping (category formation) allowing for read-across. A challenging area for in silico modelling is the prediction of chronic toxicity and the No Observed (Adverse) Effect Level (NO(A)EL) in particular. A proposed solution to the prediction of chronic toxicity is to consider organ level effects, as opposed to modelling the NO(A)EL itself. This review has focussed on the use of structural alerts to identify potential liver toxicants. In silico profilers, or groups of structural alerts, have been developed based on mechanisms of action and informed by current knowledge of Adverse Outcome Pathways. These profilers are robust and can be coded computationally to allow for prediction. However, they do not cover all mechanisms or modes of liver toxicity and recommendations for the improvement of these approaches are given.

Xanthine oxidase inhibitors beyond allopurinol and febuxostat; an overview and selection of potential leads based on in silico calculated physico-chemical properties, predicted pharmacokinetics and toxicity

Xanthine oxidase (XO), a versatile metalloflavoprotein enzyme, catalyzes the oxidative hydroxylation of hypoxanthine and xanthine to uric acid in purine catabolism while simultaneously producing reactive oxygen species. Both lead to the gout-causing hyperuricemia and oxidative damage of the tissues where overactivity of XO is present. Over the past years, significant progress and efforts towards the discovery and development of new XO inhibitors have been made and we believe that not only experts in the field, but also general readership would benefit from a review that addresses this topic. Accordingly, the aim of this article was to overview and select the most potent recently reported XO inhibitors and to compare their structures, mechanisms of action, potency and effectiveness of their inhibitory activity, in silico calculated physico-chemical properties as well as predicted pharmacokinetics and toxicity. Derivatives of imidazole, 1,3-thiazole and pyrimidine proved to be more potent than febuxostat while also displaying/possessing favorable predicted physico-chemical, pharmacokinetic and toxicological properties. Although being structurally similar to febuxostat, these optimized inhibitors bear some structural freshness and could be adopted as hits for hit-to-lead development and further evaluation by in vivo studies towards novel drug candidates, and represent valuable model structures for design of novel XO inhibitors.

General baseline toxicity QSAR for nonpolar, polar and ionisable chemicals and their mixtures in the bioluminescence inhibition assay with Aliivibrio fischeri

The Microtox assay, a bioluminescence inhibition assay with the marine bacterium Aliivibrio fischeri, is one of the most popular bioassays for assessing the cytotoxicity of organic chemicals, mixtures and environmental samples. Most environmental chemicals act as baseline toxicants in this short-term screening assay, which is typically run with only 30 min of exposure duration. Numerous Quantitative Structure-Activity Relationships (QSARs) exist for the Microtox assay for nonpolar and polar narcosis. However, typical water pollutants, which have highly diverse structures covering a wide range of hydrophobicity and speciation from neutral to anionic and cationic, are often outside the applicability domain of these QSARs. To include all types of environmentally relevant organic pollutants we developed a general baseline toxicity QSAR using liposome-water distribution ratios as descriptors. Previous limitations in availability of experimental liposome-water partition constants were overcome by reliable prediction models based on polyparameter linear free energy relationships for neutral chemicals and the COSMOmic model for charged chemicals. With this QSAR and targeted mixture experiments we could demonstrate that ionisable chemicals fall in the applicability domain. Most investigated water pollutants acted as baseline toxicants in this bioassay, with the few outliers identified as uncouplers or reactive toxicants. The main limitation of the Microtox assay is that chemicals with a high melting point and/or high hydrophobicity were outside of the applicability domain because of their low water solubility. We quantitatively derived a solubility cut-off but also demonstrated with mixture experiments that chemicals inactive on their own can contribute to mixture toxicity, which is highly relevant for complex environmental mixtures, where these chemicals may be present at concentrations below the solubility cut-off.

(Q)SARs to predict environmental toxicities: current status and future needs

The current state of the art of (Quantitative) Structure-Activity Relationships ((Q)SARs) to predict environmental toxicity is assessed along with recommendations to develop these models further. The acute toxicity of compounds acting by the non-polar narcotic mechanism of action can be well predicted, however other approaches, including read-across, may be required for compounds acting by specific mechanisms of action. The chronic toxicity of compounds to environmental species is more difficult to predict from (Q)SARs, with robust data sets and more mechanistic information required. In addition, the toxicity of mixtures is little addressed by (Q)SAR approaches. Developments in environmental toxicology including Adverse Outcome Pathways (AOPs) and omics responses should be utilised to develop better, more mechanistically relevant, (Q)SAR models.

The adverse outcome pathway for skin sensitisation: Moving closer to replacing animal testing

This article outlines the work of the Organisation for Economic Co-operation and Development (OECD) that led to being jointly awarded the 2015 Lush Black Box Prize. The award-winning work centred on the development of 'The Adverse Outcome Pathway for Skin Sensitisation Initiated by Covalent Binding to Proteins'. This Adverse Outcome Pathway (AOP) has provided the mechanistic basis for the integration of skin sensitisation-related information. Recent developments in integrated approaches to testing and assessment, based on the AOP, are summarised. The impact of the AOP on regulatory policy and on the Three Rs are discussed. An overview of the next generation of the skin sensitisation AOP module in the OECD QSAR Toolbox, based on more-recent work at the Laboratory of Mathematical Chemistry, is also presented.

Kinetic-Based Reactivity for Michael Acceptors: Structural Activity Relationships and Its Relationship to Excess Acute Fish Toxicity

Acute aquatic toxicity is divided into the "physical" mode governed by weak, non-covalent interactions and the "chemical" mode governed by covalent reactions. The potency of chemical interactions is typically expected to be greater than that for physical ones. This enhanced potency is called "excess" toxicity. As databases have become complex, substances thought to elicit a chemical mode reveal a lack of excess toxicity. One mechanism where the latter is prevalent is Michael-type addition. A series of α-β-unsaturated substances were evaluated for reactivity. Second order rate constants (k') were calculated (M^-1 s^-1) and found to vary from >4000 to <0.0003. The electron-withdrawing capacity of the polar group impacts k' values; the sequence is nitro > carbonyl or sulfone ≫ sulfoxide, nitrile or amide. When the α-carbon and/or the β-carbon of the π-system are substituted, the k' value is sharply reduced. Excess toxicity is associated with k' values >0.01 (M^-1 s^-1).

Adverse Outcome Pathway (AOP) Informed Modeling of Aquatic Toxicology: QSARs, Read-Across, and Interspecies Verification of Modes of Action

Alternative approaches have been promoted to reduce the number of vertebrate and invertebrate animals required for the assessment of the potential of compounds to cause harm to the aquatic environment. A key philosophy in the development of alternatives is a greater understanding of the relevant adverse outcome pathway (AOP). One alternative method is the fish embryo toxicity (FET) assay. Although the trends in potency have been shown to be equivalent in embryo and adult assays, a detailed mechanistic analysis of the toxicity data has yet to be performed; such analysis is vital for a full understanding of the AOP. The research presented herein used an updated implementation of the Verhaar scheme to categorize compounds into AOP-informed categories. These were then used in mechanistic (quantitative) structure-activity relationship ((Q)SAR) analysis to show that the descriptors governing the distinct mechanisms of acute fish toxicity are capable of modeling data from the FET assay. The results show that compounds do appear to exhibit the same mechanisms of toxicity across life stages. Thus, this mechanistic analysis supports the argument that the FET assay is a suitable alternative testing strategy for the specified mechanisms and that understanding the AOPs is useful for toxicity prediction across test systems.

Identification of Small-Molecule Frequent Hitters of Glutathione S-Transferase–Glutathione Interaction

In high-throughput screening (HTS) campaigns, the binding of glutathione S-transferase (GST) to glutathione (GSH) is used for detection of GST-tagged proteins in protein-protein interactions or enzyme assays. However, many false-positives, so-called frequent hitters (FH), arise that either prevent GST/GSH interaction or interfere with assay signal generation or detection. To identify GST-FH compounds, we analyzed the data of five independent AlphaScreen-based screening campaigns to classify compounds that inhibit the GST/GSH interaction. We identified 53 compounds affecting GST/GSH binding but not influencing His-tag/Ni2+-NTA interaction and general AlphaScreen signals. The structures of these 53 experimentally identified GST-FHs were analyzed in chemoinformatic studies to categorize substructural features that promote interference with GST/GSH binding. Here, we confirmed several existing chemoinformatic filters and more importantly extended them as well as added novel filters that specify compounds with anti–GST/GSH activity. Selected compounds were also tested using different antibody-based GST detection technologies and exhibited no interference clearly demonstrating specificity toward their GST/GSH interaction. Thus, these newly described GST-FH will further contribute to the identification of FH compounds containing promiscuous substructures. The developed filters were uploaded to the OCHEM website ( http://ochem.eu ) and are publicly accessible for analysis of future HTS results.

Bioalerts: a python library for the derivation of structural alerts from bioactivity and toxicity data sets

BACKGROUND

Assessing compound toxicity at early stages of the drug discovery process is a crucial task to dismiss drug candidates likely to fail in clinical trials. Screening drug candidates against structural alerts, i.e. chemical fragments associated to a toxicological response prior or after being metabolized (bioactivation), has proved a valuable approach for this task. During the last decades, diverse algorithms have been proposed for the automatic derivation of structural alerts from categorical toxicity data sets.

RESULTS AND CONCLUSIONS

Here, the python library bioalerts is presented, which comprises functionalities for the automatic derivation of structural alerts from categorical (dichotomous), e.g. toxic/non-toxic, and continuous bioactivity data sets, e.g. [Formula: see text] or [Formula: see text] values. The library bioalerts relies on the RDKit implementation of the circular Morgan fingerprint algorithm to compute chemical substructures, which are derived by considering radial atom neighbourhoods of increasing bond radius. In addition to the derivation of structural alerts, bioalerts provides functionalities for the calculation of unhashed (keyed) Morgan fingerprints, which can be used in predictive bioactivity modelling with the advantage of allowing for a chemically meaningful deconvolution of the chemical space. Finally, bioalerts provides functionalities for the easy visualization of the derived structural alerts.

Discrimination of excess toxicity from baseline level for ionizable compounds: Effect of pH

The toxic effect can be affected by pH in water through affecting the degree of ionization of ionizable compounds. Wrong classification of mode of action can be made from the apparent toxicities. In this paper, the toxicity data of 61 compounds to Daphnia magna determined at three pH values were used to investigate the effect of pH on the discrimination of excess toxicity. The results show that the apparent toxicities are significantly less than the baseline level. Analysis on the effect of pH on bioconcentration factor (BCF) shows that the log BCF values are significantly over-estimated for the strongly ionizable compounds, leading to the apparent toxicities greatly less than the baseline toxicities and the toxic ratios greatly less than zero. A theoretical equation between the apparent toxicities and pH has been developed basing on the critical body residue (CBR). The apparent toxicities are non-linearly related to pH, but linearly to fraction of unionized form. The determined apparent toxicities are well fitted with the toxicities predicted by the equation. The toxicities in the unionized form calculated from the equation are close to, or greater than the baseline level for almost all the strongly ionizable compounds, which are very different from the apparent toxicities. The studied ionizable compounds can be either classified as baseline, less inert or reactive compounds in D. magna toxicity. Some ionizable compounds do not exhibit excess toxicity at a certain pH, due not to their poor reactivity with target molecules, but because of the ionization in water.

Comparison of Toxicities to Vibrio fischeri and Fish Based on Discrimination of Excess Toxicity from Baseline Level

Investigations on the relationship of toxicities between species play an important role in the understanding of toxic mechanisms to environmental organisms. In this paper, the toxicity data of 949 chemicals to fish and 1470 chemicals to V. fischeri were used to investigate the modes of action (MOAs) between species. The results show that although there is a positive interspecies correlation, the relationship is poor. Analysis on the excess toxicity calculated from toxic ratios (TR) shows that many chemicals have close toxicities and share the same MOAs between the two species. Linear relationships between the toxicities and octanol/water partition coefficient (log KOW) for baseline and less inert compounds indicate that the internal critical concentrations (CBRs) approach a constant both to fish and V. fischeri for neutral hydrophobic compounds. These compounds share the same toxic mechanisms and bio-uptake processes between species. On the other hand, some hydrophilic compounds exhibit different toxic effects with greatly different log TR values between V. fischeri and fish species. These hydrophilic compounds were identified as reactive MOAs to V. fischeri, but not to fish. The interspecies correlation is improved by adding a hydrophobic descriptor into the correlation equation. This indicates that the differences in the toxic ratios between fish and V. fischeri for these hydrophilic compounds can be partly attributed to the differences of bioconcentration between the two species, rather than the differences of reactivity with the target macromolecules. These hydrophilic compounds may more easily pass through the cell membrane of V. fischeri than the gill and skin of fish, react with the target macromolecules and exhibit excess toxicity. The compounds with log KOW > 7 exhibiting very low toxicity (log TR < -1) to both species indicate that the bioconcentration potential of a chemical plays a very important role in the identification of excess toxicity and MOAs.

Application of in silico and in vitro methods in the development of adverse outcome pathway constructs in wildlife

There is a long history of using both in silico and in vitro methods to predict adverse effects in humans and environmental species where toxicity data are lacking. Currently, there is a great deal of interest in applying these methods to the development of so-called 'adverse outcome pathway' (AOP) constructs. The AOP approach provides a framework for organizing information at the chemical and biological level, allowing evidence from both in silico and in vitro studies to be rationally combined to fill gaps in knowledge concerning toxicological events. Fundamental to this new paradigm is a greater understanding of the mechanisms of toxicity and, in particular, where these mechanisms may be conserved across taxa, such as between model animals and related wild species. This presents an opportunity to make predictions across diverse species, where empirical data are unlikely to become available as is the case for most species of wildlife.

Development of an in Silico Profiler for Mitochondrial Toxicity

This study outlines the analysis of mitochondrial toxicity for a variety of pharmaceutical drugs extracted from Zhang et al. ((2009) Toxicol. In Vitro, 23, 134-140). These chemicals were grouped into categories based upon structural similarity. Subsequently, mechanistic analysis was undertaken for each category to identify the molecular initiating event driving mitochondrial toxicity. The mechanistic information elucidated during the analysis enabled mechanism-based structural alerts to be developed and combined together to form an in silico profiler. This profiler is envisaged to be used to develop chemical categories based upon similar mechanisms as part of the adverse outcome pathway paradigm. Additionally, the profiler could be utilized in screening large data sets in order to identify chemicals with the potential to induce mitochondrial toxicity.

Automated computational screening of the thiol reactivity of substituted alkenes

Electrophilic olefins can react with the S-H moiety of cysteine side chains. The formation of a covalent adduct through this mechanism can result in the inhibition of an enzyme. The reactivity of an olefin towards cysteine depends on its functional groups. In this study, 325 reactions of thiol-Michael-type additions to olefins were modeled using density functional theory. All combinations of ethenes with hydrogen, methyl ester, amide, and cyano substituents were included. An automated workflow was developed to perform the construction, conformation search, minimization, and calculation of molecular properties for the reactant, carbanion intermediate, and thioether products for a model reaction of the addition of methanethiol to the electrophile. Known cysteine-reactive electrophiles present in the database were predicted to react exergonically with methanethiol through a carbanion with a stability in the 30-40 kcal mol(-1) range. 13 other compounds in our database that are also present in the PubChem database have similar properties. Natural bond orbital parameters were computed and regression analysis was used to determine the relationship between properties of the olefin electronic structure and the product and intermediate stability. The stability of the intermediates is very sensitive to electronic effects on the carbon where the anionic charge is centered. The stability of the products is more sensitive to steric factors.