Using membrane-water partition coefficients in a critical membrane burden approach to aid the identification of neutral and ionizable chemicals that induce acute toxicity below narcosis levels.

The risk assessment of thousands of chemicals used in our society benefits from adequate grouping of chemicals based on the mode and mechanism of toxic action (MoA). We measure the phospholipid membrane-water distribution ratio (DMLW) using a chromatographic assay (IAM-HPLC) for 121 neutral and ionized organic chemicals and screen other methods to derive DMLW. We use IAM-HPLC based DMLW as a chemical property to distinguish between baseline narcosis and specific MoA, for reported acute toxicity endpoints on two separate sets of chemicals. The first set comprised 94 chemicals of US EPA's acute fish toxicity database: 47 categorized as narcosis MoA, 27 with specific MoA, and 20 predominantly ionic chemicals with mostly unknown MoA. The narcosis MoA chemicals clustered around the median narcosis critical membrane burden (CMBnarc) of 140 mmol kg-1 lipid, with a lower limit of 14 mmol kg-1 lipid, including all chemicals labelled Narcosis_I and Narcosis_II. This maximum 'toxic ratio' (TR) between CMBnarc and the lower limit narcosis endpoint is thus 10. For 23/28 specific MoA chemicals a TR >10 was derived, indicative of a specific adverse effect pathway related to acute toxicity. For 10/12 cations categorized as "unsure amines", the TR <10 suggests that these affect fish via narcosis MoA. The second set comprised 29 herbicides, including 17 dissociated acids, and evaluated the TR for acute toxic effect concentrations to likely sensitive aquatic plant species (green algae and macrophytes Lemna and Myriophyllum), and non-target animal species (invertebrates and fish). For 21/29 herbicides, a TR >10 indicated a specific toxic mode of action other than narcosis for at least one of these aquatic primary producers. Fish and invertebrate TRs were mostly <10, particularly for neutral herbicides, but for acidic herbicides a TR >10 indicated specific adverse effects in non-target animals. The established critical membrane approach to derive the TR provides for useful contribution to the weight of evidence to bin a chemical as having a narcosis MoA or less likely to have acute toxicity caused by a more specific adverse effect pathway. After proper calibration, the chromatographic assay provides consistent and efficient experimental input for both neutral and ionizable chemicals to this approach.

Construction of an In Silico Structural Profiling Tool Facilitating Mechanistically Grounded Classification of Aquatic Toxicants.

The performance of chemical safety assessment within the domain of environmental toxicology is often impeded by a shortfall of appropriate experimental data describing potential hazards across the many compounds in regular industrial use. In silico schemes for assigning aquatic-relevant modes or mechanisms of toxic action to substances, based solely on consideration of chemical structure, have seen widespread employment─including those of Verhaar, Russom, and later Bauer (MechoA). Recently, development of a further system was reported by Sapounidou, which, in common with MechoA, seeks to ground its classifications in understanding and appreciation of molecular initiating events. Until now, this Sapounidou scheme has not seen implementation as a tool for practical screening use. Accordingly, the primary purpose of this study was to create such a resource─in the form of a computational workflow. This exercise was facilitated through the formulation of 183 structural alerts/rules describing molecular features associated with narcosis, chemical reactivity, and specific mechanisms of action. Output was subsequently compared relative to that of the three aforementioned alternative systems to identify strengths and shortcomings as regards coverage of chemical space.

A scheme to evaluate structural alerts to predict toxicity - Assessing confidence by characterising uncertainties.

Structure-activity relationships (SARs) in toxicology have enabled the formation of structural rules which, when coded as structural alerts, are essential tools in in silico toxicology. Whilst other in silico methods have approaches for their evaluation, there is no formal process to assess the confidence that may be associated with a structural alert. This investigation proposes twelve criteria to assess the uncertainty associated with structural alerts, allowing for an assessment of confidence. The criteria are based around the stated purpose, description of the chemistry, toxicology and mechanism, performance and coverage, as well as corroborating and supporting evidence of the alert. Alerts can be given a confidence assessment and score, enabling the identification of areas where more information may be beneficial. The scheme to evaluate structural alerts was placed in the context of various use cases for industrial and regulatory applications. The analysis of alerts, and consideration of the evaluation scheme, identifies the different characteristics an alert may have, such as being highly specific or generic. These characteristics may determine when an alert can be used for specific uses such as identification of analogues for read-across or hazard identification.

Data-driven learning of narcosis mode of action identifies a CNS transcriptional signature shared between whole organism Caenorhabditis elegans and a fish gill cell line.

With the large numbers of man-made chemicals produced and released in the environment, there is a need to provide assessments on their potential effects on environmental safety and human health. Current regulatory frameworks rely on a mix of both hazard and risk-based approaches to make safety decisions, but the large number of chemicals in commerce combined with an increased need to conduct assessments in the absence of animal testing makes this increasingly challenging. This challenge is catalysing the use of more mechanistic knowledge in safety assessment from both in silico and in vitro approaches in the hope that this will increase confidence in being able to identify modes of action (MoA) for the chemicals in question. Here we approach this challenge by testing whether a functional genomics approach in C. elegans and in a fish cell line can identify molecular mechanisms underlying the effects of narcotics, and the effects of more specific acting toxicants. We show that narcosis affects the expression of neuronal genes associated with CNS function in C. elegans and in a fish cell line. Overall, we believe that our study provides an important step in developing mechanistically relevant biomarkers which can be used to screen for hazards, and which prevent the need for repeated animal or cross-species comparisons for each new chemical.

Development of an Enhanced Mechanistically Driven Mode of Action Classification Scheme for Adverse Effects on Environmental Species.

This study developed a novel classification scheme to assign chemicals to a verifiable mechanism of (eco-)toxicological action to allow for grouping, read-across, and in silico model generation. The new classification scheme unifies and extends existing schemes and has, at its heart, direct reference to molecular initiating events (MIEs) promoting adverse outcomes. The scheme is based on three broad domains of toxic action representing nonspecific toxicity (e.g., narcosis), reactive mechanisms (e.g., electrophilicity and free radical action), and specific mechanisms (e.g., associated with enzyme inhibition). The scheme is organized at three further levels of detail beyond broad domains to separate out the mechanistic group, specific mechanism, and the MIEs responsible. The novelty of this approach comes from the reference to taxonomic diversity within the classification, transparency, quality of supporting evidence relating to MIEs, and that it can be updated readily.

Confidence in Inactive and Active Predictions from Structural Alerts.

Having a measure of confidence in computational predictions of biological activity from in silico tools is vital when making predictions for new chemicals, for example, in chemical risk assessment. Where predictions of biological activity are used as an indicator of a potential hazard, false-negative predictions are the most concerning prediction; however, assigning confidence in inactive predictions is particularly challenging. How can one confidently identify the absence of activating features? In this study, we present methods for assigning confidence to both active and inactive predictions from structural alerts for protein-binding molecular initiating events (MIEs). Structural alerts were derived through an iterative statistical method. Confidence in the activity predictions is assigned by measuring the Tanimoto similarity between Morgan fingerprints of chemicals in the test set to relevant chemicals in the training set, and suitable cutoff values have been defined to give different confidence categories. To avoid a potential compound series bias in the test set and hence overestimate the performance of the method, we measured the biological activity of 27 compounds with 24 proteins, which gave us an additional 648 experimental measurements; many of the measurements are currently nonexistent in the literature and databases. This data set was complemented with newly measured biological activities published in ChEMBL25 and formed a combined independent validation data set. Applying the confidence categories to the computational predictions for the new data leads to the identification of chemicals for which one should be confident of either an inactive or active prediction, allowing model predictions to be used responsibly.

In Silico Guidance for In Vitro Androgen and Glucocorticoid Receptor ToxCast Assays.

Molecular initiating events (MIEs) are key events in adverse outcome pathways that link molecular chemistry to target biology. As they are based on chemistry, these interactions are excellent targets for computational chemistry approaches to in silico modeling. In this work, we aim to link ligand chemical structures to MIEs for androgen receptor (AR) and glucocorticoid receptor (GR) binding using ToxCast data. This has been done using an automated computational algorithm to perform maximal common substructure searches on chemical binders for each target from the ToxCast dataset. The models developed show a high level of accuracy, correctly assigning 87.20% of AR binders and 96.81% of GR binders in a 25% test set using holdout cross-validation. The 2D structural alerts developed can be used as in silico models to predict these MIEs and as guidance for in vitro ToxCast assays to confirm hits. These models can target such experimental work, reducing the number of assays to be performed to gain required toxicological insight. Development of these models has also allowed some structural alerts to be identified as predictors for agonist or antagonist behavior at the receptor target. This work represents a first step in using computational methods to guide and target experimental approaches.

Neural network activation similarity: a new measure to assist decision making in chemical toxicology.

Deep learning neural networks, constructed for the prediction of chemical binding at 79 pharmacologically important human biological targets, show extremely high performance on test data (accuracy 92.2 ± 4.2%, MCC 0.814 ± 0.093 and ROC-AUC 0.96 ± 0.04). A new molecular similarity measure, Neural Network Activation Similarity, has been developed, based on signal propagation through the network. This is complementary to standard Tanimoto similarity, and the combined use increases confidence in the computer's prediction of activity for new chemicals by providing a greater understanding of the underlying justification. The in silico prediction of these human molecular initiating events is central to the future of chemical safety risk assessment and improves the efficiency of safety decision making.

Structural Alerts and Random Forest Models in a Consensus Approach for Receptor Binding Molecular Initiating Events.

A molecular initiating event (MIE) is the gateway to an adverse outcome pathway (AOP), a sequence of events ending in an adverse effect. In silico predictions of MIEs are a vital tool in a modern, mechanism-focused approach to chemical risk assessment. For 90 biological targets representing important human MIEs, structural alert-based models have been constructed with an automated procedure that uses Bayesian statistics to iteratively select substructures. These models give impressive average performance statistics (an average of 92% correct predictions across targets), significantly improving on previous models. Random Forest models have been constructed from physicochemical features for the same targets, giving similarly impressive performance statistics (93% correct predictions). A key difference between the models is interpretation of predictions-the structural alert models are transparent and easy to interpret, while Random Forest models can only identify the most important physicochemical features for making predictions. The two complementary models have been combined in a consensus model, improving performance compared to each individual model (94% correct predictions) and increasing confidence in predictions. Variation in model performance has been explained by calculating a modelability index (MODI), using Tanimoto coefficient between Morgan fingerprints to identify nearest neighbor chemicals. This work is an important step toward building confidence in the use of in silico tools for assessment of toxicity.

Quantitative Predictions for Molecular Initiating Events Using Three-Dimensional Quantitative Structure-Activity Relationships.

The aim of human toxicity risk assessment is to determine a safe dose or exposure to a chemical for humans. This requires an understanding of the exposure of a person to a chemical and how much of the chemical is required to cause an adverse effect. To do this computationally, we need to understand how much of a chemical is required to perturb normal biological function in an adverse outcome pathway (AOP). The molecular initiating event (MIE) is the first step in an adverse outcome pathway and can be considered as a chemical interaction between a chemical toxicant and a biological molecule. Key chemical characteristics can be identified and used to model the chemistry of these MIEs. In this study, we do just this by using chemical substructures to categorize chemicals and 3D quantitative structure-activity relationships (QSARs) based on comparative molecular field analysis (CoMFA) to calculate molecular activity. Models have been constructed across a variety of human biological targets, the glucocorticoid receptor, mu opioid receptor, cyclooxygenase-2 enzyme, human ether-à-go-go related gene channel, and dopamine transporter. These models tend to provide molecular activity estimation well within one log unit and electronic and steric fields that can be visualized to better understand the MIE and biological target of interest. The outputs of these fields can be used to identify key aspects of a chemical's chemistry which can be changed to reduce its ability to activate a given MIE. With this methodology, the quantitative chemical activity can be predicted for a wide variety of MIEs, which can feed into AOP-based chemical risk assessments, and understanding of the chemistry behind the MIE can be gained.

Vision of a near future: Bridging the human health-environment divide. Toward an integrated strategy to understand mechanisms across species for chemical safety assessment.

There is a growing recognition that application of mechanistic approaches to understand cross-species shared molecular targets and pathway conservation in the context of hazard characterization, provide significant opportunities in risk assessment (RA) for both human health and environmental safety. Specifically, it has been recognized that a more comprehensive and reliable understanding of similarities and differences in biological pathways across a variety of species will better enable cross-species extrapolation of potential adverse toxicological effects. Ultimately, this would also advance the generation and use of mechanistic data for both human health and environmental RA. A workshop brought together representatives from industry, academia and government to discuss how to improve the use of existing data, and to generate new NAMs data to derive better mechanistic understanding between humans and environmentally-relevant species, ultimately resulting in holistic chemical safety decisions. Thanks to a thorough dialogue among all participants, key challenges, current gaps and research needs were identified, and potential solutions proposed. This discussion highlighted the common objective to progress toward more predictive, mechanistically based, data-driven and animal-free chemical safety assessments. Overall, the participants recognized that there is no single approach which would provide all the answers for bridging the gap between mechanism-based human health and environmental RA, but acknowledged we now have the incentive, tools and data availability to address this concept, maximizing the potential for improvements in both human health and environmental RA.

Improved Algal Toxicity Test System for Robust Omics-Driven Mode-of-Action Discovery in Chlamydomonas reinhardtii.

Algae are key components of aquatic food chains. Consequently, they are internationally recognised test species for the environmental safety assessment of chemicals. However, existing algal toxicity test guidelines are not yet optimized to discover molecular modes of action, which require highly-replicated and carefully controlled experiments. Here, we set out to develop a robust, miniaturised and scalable Chlamydomonas reinhardtii toxicity testing approach tailored to meet these demands. We primarily investigated the benefits of synchronised cultures for molecular studies, and of exposure designs that restrict chemical volatilisation yet yield sufficient algal biomass for omics analyses. Flow cytometry and direct-infusion mass spectrometry metabolomics revealed significant and time-resolved changes in sample composition of synchronised cultures. Synchronised cultures in sealed glass vials achieved adequate growth rates at previously unachievably-high inoculation cell densities, with minimal pH drift and negligible chemical loss over 24-h exposures. Algal exposures to a volatile test compound (chlorobenzene) yielded relatively high reproducibility of metabolic phenotypes over experimental repeats. This experimental test system extends existing toxicity testing formats to allow highly-replicated, omics-driven, mode-of-action discovery.

Chemical hazard prediction and hypothesis testing using quantitative adverse outcome pathways.

Current efforts in chemical safety are focused on utilizing human in vitro or alternative animal data in biological pathway context. However, it remains unclear how biological pathways, and toxicology data developed in that context, can be used to quantitatively facilitate decision-making.  The objective of this work is to determine if hypothesis testing using Adverse Outcome Pathways (AOPs) can provide quantitative chemical hazard predictions.  Current methods for predicting hazards of chemicals in a biological pathway context were extensively reviewed, specific case studies examined and computational modeling used to demonstrate quantitative hazard prediction based on an AOP. Since AOPs are chemically agnostic, we propose that AOPs function as hypotheses for how specific chemicals may cause adverse effects via specific pathways. Three broad approaches were identified for testing the hypothesis with AOPs, semi-quantitative weight of evidence, probabilistic, and mechanistic modeling. We then demonstrate how these approaches could be used to test hypotheses using high throughput in vitro data and alternative animal data. Finally, we discuss standards in development and documentation that would facilitate use in a regulatory context. We conclude that quantitative AOPs provide a flexible hypothesis framework for predicting chemical hazards. It accommodates a wide range of approaches that are useful at many stages and build upon one another to become increasingly quantitative.

Using 2D Structural Alerts to Define Chemical Categories for Molecular Initiating Events.

Molecular initiating events (MIEs) are important concepts for in silico predictions. They can be used to link chemical characteristics to biological activity through an adverse outcome pathway (AOP). In this work, we capture chemical characteristics in 2D structural alerts, which are then used as models to predict MIEs. An automated procedure has been used to identify these alerts, and the chemical categories they define have been used to provide quantitative predictions for the activity of molecules that contain them. This has been done across a diverse group of 39 important pharmacological human targets using open source data. The alerts for each target combine into a model for that target, and these models are joined into a tool for MIE prediction with high average model performance (sensitivity = 82%, specificity = 93%, overall quality = 93%, Matthews correlation coefficient = 0.57). The result is substantially improved from our previous study (Allen, T. E. H., Goodman, J. M., Gutsell, S., and Russell, P. J. 2016. A history of the molecular initiating event. Chem. Res. Toxicol. 29, 2060-2070) for which the mean sensitivity for each target was only 58%. This tool provides the first step in an AOP-based risk assessment, linking chemical structure to toxicity endpoint.

Using Transition State Modeling To Predict Mutagenicity for Michael Acceptors.

The Ames mutagenicity assay is a long established in vitro test to measure the mutagenicity potential of a new chemical used in regulatory testing globally. One of the key computational approaches to modeling of the Ames assay relies on the formation of chemical categories based on the different electrophilic compounds that are able to react directly with DNA and form a covalent bond. Such approaches sometimes predict false positives, as not all Michael acceptors are found to be Ames-positive. The formation of such covalent bonds can be explored computationally using density functional theory transition state modeling. We have applied this approach to mutagenicity, allowing us to calculate the activation energy required for α,ß-unsaturated carbonyls to react with a model system for the guanine nucleobase of DNA. These calculations have allowed us to identify that chemical compounds with activation energies greater than or equal to 25.7 kcal/mol are not able to bind directly to DNA. This allows us to reduce the false positive rate for computationally predicted mutagenicity assays. This methodology can be used to investigate other covalent-bond-forming reactions that can lead to toxicological outcomes and learn more about experimental results.

In silico toxicology protocols.

The present publication surveys several applications of in silico (i.e., computational) toxicology approaches across different industries and institutions. It highlights the need to develop standardized protocols when conducting toxicity-related predictions. This contribution articulates the information needed for protocols to support in silico predictions for major toxicological endpoints of concern (e.g., genetic toxicity, carcinogenicity, acute toxicity, reproductive toxicity, developmental toxicity) across several industries and regulatory bodies. Such novel in silico toxicology (IST) protocols, when fully developed and implemented, will ensure in silico toxicological assessments are performed and evaluated in a consistent, reproducible, and well-documented manner across industries and regulatory bodies to support wider uptake and acceptance of the approaches. The development of IST protocols is an initiative developed through a collaboration among an international consortium to reflect the state-of-the-art in in silico toxicology for hazard identification and characterization. A general outline for describing the development of such protocols is included and it is based on in silico predictions and/or available experimental data for a defined series of relevant toxicological effects or mechanisms. The publication presents a novel approach for determining the reliability of in silico predictions alongside experimental data. In addition, we discuss how to determine the level of confidence in the assessment based on the relevance and reliability of the information.

The Role of Omics in the Application of Adverse Outcome Pathways for Chemical Risk Assessment.

In conjunction with the second International Environmental Omics Symposium (iEOS) conference, held at the University of Liverpool (United Kingdom) in September 2014, a workshop was held to bring together experts in toxicology and regulatory science from academia, government and industry. The purpose of the workshop was to review the specific roles that high-content omics datasets (eg, transcriptomics, metabolomics, lipidomics, and proteomics) can hold within the adverse outcome pathway (AOP) framework for supporting ecological and human health risk assessments. In light of the growing number of examples of the application of omics data in the context of ecological risk assessment, we considered how omics datasets might continue to support the AOP framework. In particular, the role of omics in identifying potential AOP molecular initiating events and providing supportive evidence of key events at different levels of biological organization and across taxonomic groups was discussed. Areas with potential for short and medium-term breakthroughs were also discussed, such as providing mechanistic evidence to support chemical read-across, providing weight of evidence information for mode of action assignment, understanding biological networks, and developing robust extrapolations of species-sensitivity. Key challenges that need to be addressed were considered, including the need for a cohesive approach towards experimental design, the lack of a mutually agreed framework to quantitatively link genes and pathways to key events, and the need for better interpretation of chemically induced changes at the molecular level. This article was developed to provide an overview of ecological risk assessment process and a perspective on how high content molecular-level datasets can support the future of assessment procedures through the AOP framework.

A History of the Molecular Initiating Event.

The adverse outcome pathway (AOP) framework provides an alternative to traditional in vivo experiments for the risk assessment of chemicals. AOPs consist of a number of key events (KEs) linked by key event relationships across a range of biological organization backed by scientific evidence. The first KE in the pathway is the molecular initiating event (MIE)-the initial chemical trigger that starts an AOP. Over the past 3 years the AOP conceptual framework has gained a large amount of momentum in toxicology as an alternative to animal methods, and so the MIE has come into the spotlight. What is an MIE? How can MIEs be measured or predicted? What research is currently contributing to our understanding of MIEs? In this Perspective we outline answers to these key questions.

Using Molecular Initiating Events To Generate 2D Structure-Activity Relationships for Toxicity Screening.

Molecular initiating events (MIEs) can be boiled down to chemical interactions. Chemicals that interact must have intrinsic properties that allow them to exhibit this behavior, be these properties stereochemical, electronic, or otherwise. In an attempt to discover some of these chemical characteristics, we have constructed structural alert-style structure-activity relationships (SARs) to computationally predict MIEs. This work utilizes chemical informatics approaches, searching the ChEMBL database for molecules that bind to a number of pharmacologically important human toxicology targets, including G-protein coupled receptors, enzymes, ion channels, nuclear receptors, and transporters. By screening these compounds to find common 2D fragments and combining this approach with a good understanding of the literature, bespoke 2D structural alerts have been written. These SARs form the beginning of a tool for screening novel chemicals to establish the kind of interactions that they may be able to make in humans. These SARs have been run through an internal validation to test their quality, and the results of this are also discussed. MIEs have proven to be difficult to find and characterize, but we believe we have taken a key first step with this work.

Adverse Outcome Pathways for Regulatory Applications: Examination of Four Case Studies With Different Degrees of Completeness and Scientific Confidence.

Adverse outcome pathways (AOPs) offer a pathway-based toxicological framework to support hazard assessment and regulatory decision-making. However, little has been discussed about the scientific confidence needed, or how complete a pathway should be, before use in a specific regulatory application. Here we review four case studies to explore the degree of scientific confidence and extent of completeness (in terms of causal events) that is required for an AOP to be useful for a specific purpose in a regulatory application: (i) Membrane disruption (Narcosis) leading to respiratory failure (low confidence), (ii) Hepatocellular proliferation leading to cancer (partial pathway, moderate confidence), (iii) Covalent binding to proteins leading to skin sensitization (high confidence), and (iv) Aromatase inhibition leading to reproductive dysfunction in fish (high confidence). Partially complete AOPs with unknown molecular initiating events, such as 'Hepatocellular proliferation leading to cancer', were found to be valuable. We demonstrate that scientific confidence in these pathways can be increased though the use of unconventional information (eg, computational identification of potential initiators). AOPs at all levels of confidence can contribute to specific uses. A significant statistical or quantitative relationship between events and/or the adverse outcome relationships is a common characteristic of AOPs, both incomplete and complete, that have specific regulatory uses. For AOPs to be useful in a regulatory context they must be at least as useful as the tools that regulators currently possess, or the techniques currently employed by regulators.