Defining molecular initiating events in the adverse outcome pathway framework for risk assessment.

Consumer and environmental safety decisions are based on exposure and hazard data, interpreted using risk assessment approaches. The adverse outcome pathway (AOP) conceptual framework has been presented as a logical sequence of events or processes within biological systems which can be used to understand adverse effects and refine current risk assessment practices in ecotoxicology. This framework can also be applied to human toxicology and is explored on the basis of investigating the molecular initiating events (MIEs) of compounds. The precise definition of the MIE has yet to reach general acceptance. In this work we present a unified MIE definition: an MIE is the initial interaction between a molecule and a biomolecule or biosystem that can be causally linked to an outcome via a pathway. Case studies are presented, and issues with current definitions are addressed. With the development of a unified MIE definition, the field can look toward defining, classifying, and characterizing more MIEs and using knowledge of the chemistry of these processes to aid AOP research and toxicity risk assessment. We also present the role of MIE research in the development of in vitro and in silico toxicology and suggest how, by using a combination of biological and chemical approaches, MIEs can be identified and characterized despite a lack of detailed reports, even for some of the most studied molecules in toxicology.The authors acknowledge the financial support of Unilever.This is the accepted manuscript. The final published version is available from ACS at http://dx.doi.org/10.1021/tx500345

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

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.Unileve

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.Unileve

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.Unilever (Post-Doctoral Funding to T.E.H.A) and Girton College, Cambridge (Research Fellowship to M.N.G.) for financial support. Part of this work was performed using the Darwin Supercomputer of the University of Cambridge High Performance Computing Service (http://www.hpc.cam.ac.uk/), provided by Dell Inc. using Strategic Research Infrastructure Funding from the Higher Education Funding Council for England and funding from the Science and Technology Facilities Council

Improved Algal Toxicity Test System for Robust <i>Omics</i>-Driven Mode-of-Action Discovery in <i>Chlamydomonas reinhardtii</i>

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

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.Unileve

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.Unilever Research Participation Program of the U.S. Environmental Protection Agency, Office of Research and Development, administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the U.S. Department of Energy and the U.S. EP

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