Using Molecular Initiating Events to Develop a Structural Alert Based Screening Workflow for Nuclear Receptor Ligands Associated with Hepatic Steatosis

Computational Strategies for Assessing Adverse Outcome Pathways: Hepatic Steatosis as a Case Study

The evolving landscape of chemical risk assessment is increasingly focused on developing tiered, mechanistically driven approaches that avoid the use of animal experiments. In this context, adverse outcome pathways have gained importance for evaluating various types of chemical-induced toxicity. Using hepatic steatosis as a case study, this review explores the use of diverse computational techniques, such as structure-activity relationship models, quantitative structure-activity relationship models, read-across methods, omics data analysis, and structure-based approaches to fill data gaps within adverse outcome pathway networks. Emphasizing the regulatory acceptance of each technique, we examine how these methodologies can be integrated to provide a comprehensive understanding of chemical toxicity. This review highlights the transformative impact of in silico techniques in toxicology, proposing guidelines for their application in evidence gathering for developing and filling data gaps in adverse outcome pathway networks. These guidelines can be applied to other cases, advancing the field of toxicological risk assessment.

Combining graph neural networks and transformers for few-shot nuclear receptor binding activity prediction

Nuclear receptors (NRs) play a crucial role as biological targets in drug discovery. However, determining which compounds can act as endocrine disruptors and modulate the function of NRs with a reduced amount of candidate drugs is a challenging task. Moreover, the computational methods for NR-binding activity prediction mostly focus on a single receptor at a time, which may limit their effectiveness. Hence, the transfer of learned knowledge among multiple NRs can improve the performance of molecular predictors and lead to the development of more effective drugs. In this research, we integrate graph neural networks (GNNs) and Transformers to introduce a few-shot GNN-Transformer, Meta-GTNRP to predict the binding activity of compounds using the combined information of different NRs and identify potential NR-modulators with limited data. The Meta-GTNRP model captures the local information in graph-structured data and preserves the global-semantic structure of molecular graph embeddings for NR-binding activity prediction. Furthermore, a few-shot meta-learning approach is proposed to optimize model parameters for different NR-binding tasks and leverage the complementarity among multiple NR-specific tasks to predict binding activity of compounds for each NR with just a few labeled molecules. Experiments with a compound database containing annotations on the binding activity for 11 NRs shows that Meta-GTNRP outperforms other graph-based approaches. The data and code are available at: https://github.com/ltorres97/Meta-GTNRP .Scientific contributionThe proposed few-shot GNN-Transformer model, Meta-GTNRP captures the local structure of molecular graphs and preserves the global-semantic information of graph embeddings to predict the NR-binding activity of compounds with limited available data; A few-shot meta-learning framework adapts model parameters across NR-specific tasks for different NRs in a joint learning procedure to predict the binding activity of compounds for each NR with just a few labeled molecules in highly imbalanced data scenarios; Meta-GTNRP is a data-efficient approach that combines the strengths of GNNs and Transformers to predict the NR-binding properties of compounds through an optimized meta-learning procedure and deliver robust results valuable to identify potential NR-based drug candidates.

Identifying chemicals based on receptor binding/bioactivation/mechanistic explanation associated with potential to elicit hepatotoxicity and to support structure activity relationship-based read-across

The liver is the most common target organ in toxicology studies. The development of chemical structural alerts for identifying hepatotoxicity will play an important role in in silico model prediction and help strengthen the identification of analogs used in structure activity relationship (SAR)- based read-across. The aim of the current study is development of an SAR-based expert-system decision tree for screening of hepatotoxicants across a wide range of chemistry space and proposed modes of action for clustering of chemicals using defined core chemical categories based on receptor-binding or bioactivation. The decision tree is based on ∼ 1180 different chemicals that were reviewed for hepatotoxicity information. Knowledge of chemical receptor binding, metabolism and mechanistic information were used to group these chemicals into 16 different categories and 102 subcategories: four categories describe binders to 9 different receptors, 11 categories are associated with possible reactive metabolites (RMs) and there is one miscellaneous category. Each chemical subcategory has been associated with possible modes of action (MOAs) or similar key structural features. This decision tree can help to screen potential liver toxicants associated with core structural alerts of receptor binding and/or RMs and be used as a component of weight of evidence decisions based on SAR read-across, and to fill data gaps.

In vitro to in vivo extrapolation and high-content imaging for simultaneous characterization of chemically induced liver steatosis and markers of hepatotoxicity

Chemically induced steatosis is characterized by lipid accumulation associated with mitochondrial dysfunction, oxidative stress and nucleus distortion. New approach methods integrating in vitro and in silico models are needed to identify chemicals that may induce these cellular events as potential risk factors for steatosis and associated hepatotoxicity. In this study we used high-content imaging for the simultaneous quantification of four cellular markers as sentinels for hepatotoxicity and steatosis in chemically exposed human liver cells in vitro. Furthermore, we evaluated the results with a computational model for the extrapolation of human oral equivalent doses (OED). First, we tested 16 reference chemicals with known capacities to induce cellular alterations in nuclear morphology, lipid accumulation, mitochondrial membrane potential and oxidative stress. Then, using physiologically based pharmacokinetic modeling and reverse dosimetry, OEDs were extrapolated from data of any stimulated individual sentinel response. The extrapolated OEDs were confirmed to be within biologically relevant exposure ranges for the reference chemicals. Next, we tested 14 chemicals found in food, selected from thousands of putative chemicals on the basis of structure-based prediction for nuclear receptor activation. Amongst these, orotic acid had an extrapolated OED overlapping with realistic exposure ranges. Thus, we were able to characterize known steatosis-inducing chemicals as well as data-scarce food-related chemicals, amongst which we confirmed orotic acid to induce hepatotoxicity. This strategy addresses needs of next generation risk assessment and can be used as a first chemical prioritization hazard screening step in a tiered approach to identify chemical risk factors for steatosis and hepatotoxicity-associated events.

A 10-step framework for use of read-across (RAX) in next generation risk assessment (NGRA) for cosmetics safety assessment

This paper presents a 10-step read-across (RAX) framework for use in cases where a threshold of toxicological concern (TTC) approach to cosmetics safety assessment is not possible. RAX builds on established approaches that have existed for more than two decades using chemical properties and in silico toxicology predictions, by further substantiating hypotheses on toxicological similarity of substances, and integrating new approach methodologies (NAM) in the biological and kinetic domains. NAM include new types of data on biological observations from, for example, in vitro assays, toxicogenomics, metabolomics, receptor binding screens and uses physiologically-based kinetic (PBK) modelling to inform about systemic exposure. NAM data can help to substantiate a mode/mechanism of action (MoA), and if similar chemicals can be shown to work by a similar MoA, a next generation risk assessment (NGRA) may be performed with acceptable confidence for a data-poor target substance with no or inadequate safety data, based on RAX approaches using data-rich analogue(s), and taking account of potency or kinetic/dynamic differences.

COSMOS next generation - A public knowledge base leveraging chemical and biological data to support the regulatory assessment of chemicals

The COSMOS Database (DB) was originally established to provide reliable data for cosmetics-related chemicals within the COSMOS Project funded as part of the SEURAT-1 Research Initiative. The database has subsequently been maintained and developed further into COSMOS Next Generation (NG), a combination of database and in silico tools, essential components of a knowledge base. COSMOS DB provided a cosmetics inventory as well as other regulatory inventories, accompanied by assessment results and in vitro and in vivo toxicity data. In addition to data content curation, much effort was dedicated to data governance - data authorisation, characterisation of quality, documentation of meta information, and control of data use. Through this effort, COSMOS DB was able to merge and fuse data of various types from different sources. Building on the previous effort, the COSMOS Minimum Inclusion (MINIS) criteria for a toxicity database were further expanded to quantify the reliability of studies. COSMOS NG features multiple fingerprints for analysing structure similarity, and new tools to calculate molecular properties and screen chemicals with endpoint-related public profilers, such as DNA and protein binders, liver alerts and genotoxic alerts. The publicly available COSMOS NG enables users to compile information and execute analyses such as category formation and read-across. This paper provides a step-by-step guided workflow for a simple read-across case, starting from a target structure and culminating in an estimation of a NOAEL confidence interval. Given its strong technical foundation, inclusion of quality-reviewed data, and provision of tools designed to facilitate communication between users, COSMOS NG is a first step towards building a toxicological knowledge hub leveraging many public data systems for chemical safety evaluation. We continue to monitor the feedback from the user community at support@mn-am.com.

Novel QSAR Models for Molecular Initiating Event Modeling in Two Intersecting Adverse Outcome Pathways Based Pulmonary Fibrosis Prediction for Biocidal Mixtures

The adverse outcome pathway (AOP) was introduced as an alternative method to avoid unnecessary animal tests. Under the AOP framework, an in silico methods, molecular initiating event (MIE) modeling is used based on the ligand-receptor interaction. Recently, the intersecting AOPs (AOP 347), including two MIEs, namely peroxisome proliferator-activated receptor-gamma (PPAR-γ) and toll-like receptor 4 (TLR4), associated with pulmonary fibrosis was proposed. Based on the AOP 347, this study developed two novel quantitative structure-activity relationship (QSAR) models for the two MIEs. The prediction performances of different MIE modeling methods (e.g., molecular dynamics, pharmacophore model, and QSAR) were compared and validated with in vitro test data. Results showed that the QSAR method had high accuracy compared with other modeling methods, and the QSAR method is suitable for the MIE modeling in the AOP 347. Therefore, the two QSAR models based on the AOP 347 can be powerful models to screen biocidal mixture related to pulmonary fibrosis.

Development of a Battery of In Silico Prediction Tools for Drug-Induced Liver Injury from the Vantage Point of Translational Safety Assessment

Drug-induced liver injury (DILI) remains a challenge when translating knowledge from the preclinical stage to human use cases. Attempts to model human DILI directly based on the information from drug labels have had some success; however, the approach falls short of providing insights or addressing uncertainty due to the difficulty of decoupling the idiosyncratic nature of human DILI outcomes. Our approach in this comparative analysis is to leverage existing preclinical and clinical data as well as information on metabolism to better translate mammalian to human DILI. The human DILI knowledge base from the United States Food and Drug Administration (U.S. FDA) National Center for Toxicology Research contains 1036 pharmaceuticals from diverse therapeutic categories. A human DILI training set of 305 oral marketed drugs was prepared and a binary classification scheme applied. The second knowledge base consists of mammalian repeated dose toxicity with liver toxicity data from various regulatory sources. Within this knowledge base, we identified 278 pharmaceuticals containing 198 marketed or withdrawn oral drugs with data from the U.S. FDA new drug application and 98 active pharmaceutical ingredients from ToxCast. From this collection, a set of 225 oral drugs was prepared as the mammalian hepatotoxicity training set with particular end points of pathology findings in the liver and bile duct. Both human and mammalian data sets were processed using various learning algorithms, including artificial intelligence approaches. The external validations for both models were comparable to the training statistics. These data sets were also used to extract species-differentiating chemotypes that differentiate DILI effects on humans from mammals. A systematic workflow was devised to predict human DILI and provide mechanistic insights. For a given query molecule, both human and mammalian models are run. If the predictions are discordant, both metabolites and parents are investigated for quantitative structure-activity relationship and species-differentiating chemotypes. Their results are combined using the Dempster-Shafer decision theory to yield a final outcome prediction for human DILI with estimated uncertainty. Finally, these tools are implementable within an in silico platform for systematic evaluation.

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.

A Robust, Mechanistically Based In Silico Structural Profiler for Hepatic Cholestasis

Owing to the primary role which it holds within metabolism of xenobiotics, the liver stands at heightened risk of exposure to, and injury from, potentially hazardous substances. A principal manifestation of liver dysfunction is cholestasis-the impairment of physiological bile circulation from its point of origin within the organ to the site of action in the small intestine. The capacity for early identification of compounds liable to exert cholestatic effects is of particular utility within the field of pharmaceutical development, where contribution toward candidate attrition is great. Shortcomings associated with the present in vitro methodologies forecasting cholestasis render their predictivity questionable, permitting scope for the adoption of computational toxicology techniques. As such, the intention of this study has been to construct an in silico profiler, founded upon clinical data, highlighting structural motifs most reliably associated with the end point. Drawing upon a list of >1500 small molecular drugs, compiled and annotated by Kotsampasakou, E. and Ecker, G. F. (J. Chem. Inf. Model. 2017, 57, 608-615), we have formulated a series of 15 structural alerts. These describe fragments intrinsic within distinct pharmaceutical classes including psychoactive tricyclics, β-lactam antimicrobials, and estrogenic/androgenic steroids. Description of the coverage and selectivity of each are provided, alongside consideration of the underlying reactive mechanisms and relevant structure-activity concerns. Provision of mechanistic anchoring ensures that potential exists for framing within the adverse outcome pathway paradigm-the chemistry conveyed through the alert, in particular enabling rationalization at the level of the molecular initiating event.

In Silico Identification of Chemicals Capable of Binding to the Ecdysone Receptor

The process of molting, known alternatively as ecdysis, is a feature integral in the life cycles of species across the arthropod phylum. Regulation occurs as a function of the interaction of ecdysteroid hormones with the arthropod nuclear ecdysone receptor-a process preceding the triggering of a series of downstream events constituting an endocrine signaling pathway highly conserved throughout environmentally prevalent insect, crustacean, and myriapod organisms. Inappropriate ecdysone receptor binding and activation forms the essential molecular initiating event within possible adverse outcome pathways relating abnormal molting to mortality in arthropods. Definition of the characteristics of chemicals liable to stimulate such activity has the potential to be of great utility in mitigation of hazards posed toward vulnerable species. Thus the aim of the present study was to develop a series of rule-sets, derived from the key structural and physicochemical features associated with identified ecdysone receptor ligands, enabling construction of Konstanz Information Miner (KNIME) workflows permitting the flagging of compounds predisposed to binding at the site. Data describing the activities of 555 distinct chemicals were recovered from a variety of assays across 10 insect species, allowing for formulation of KNIME screens for potential binding activity at the molecular initiating event and adverse outcome level of biological organization. Environ Toxicol Chem 2020;39:1438-1450. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

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.

The Se…S/N interactions as a possible mechanism of δ-aminolevulinic acid dehydratase enzyme inhibition by organoselenium compounds: A computational study

•

DPDS and PSA interacts with cysteine residues from AlaD active site.

•

The Se…S interactions could be involved in the δ-AlaD inhibition.

•

δ-AlaD from Cucumis sativus does not present cysteine residues in the active site.

•

Se…N interactions could be involved in the organoselenium action.

Organoselenium compounds present many pharmacological properties and are promising drugs. However, toxicological effects associated with inhibition of thiol-containing enzymes, such as the δ-aminolevulinic acid dehydratase (δ-AlaD), have been described. The molecular mechanism(s) by which they inhibit thiol-containing enzymes at the atomic level, is still not well known. The use of computational methods to understand the physical–chemical properties and biological activity of chemicals is essential to the rational design of new drugs. In this work, we propose an in silico study to understand the δ-AlaD inhibition mechanism by diphenyl diselenide (DPDS) and its putative metabolite, phenylseleninic acid (PSA), using δ-AlaD enzymes from Homo sapiens (Hsδ-AlaD), Drosophila melanogaster (Dmδ-AlaD) and Cucumis sativus (Csδ-AlaD). Protein modeling homology, molecular docking, and DFT calculations are combined in this study. According to the molecular docking, DPDS and PSA might bind in the Hsδ-AlaD and Dmδ-AlaD active sites interacting with the cysteine residues by Se^…S interactions. On the other hand, the DPDS does not access the active site of the Csδ-AlaD (a non-thiol protein), while the PSA interacts with the amino acids residues from the active site, such as the Lys291. These interactions might lead to the formation of a covalent bond, and consequently, to the enzyme inhibition. In fact, DFT calculations (mPW1PW91/def2TZVP) demonstrated that the selenylamide bond formation is energetically favored. The in silico data showed here are in accordance with previous experimental studies, and help us to understand the reactivity and biological activity of organoselenium compounds.

Quantitative adverse outcome pathway (qAOP) models for toxicity prediction

The quantitative adverse outcome pathway (qAOP) concept is gaining interest due to its potential regulatory applications in chemical risk assessment. Even though an increasing number of qAOP models are being proposed as computational predictive tools, there is no framework to guide their development and assessment. As such, the objectives of this review were to: (i) analyse the definitions of qAOPs published in the scientific literature, (ii) define a set of common features of existing qAOP models derived from the published definitions, and (iii) identify and assess the existing published qAOP models and associated software tools. As a result, five probabilistic qAOPs and ten mechanistic qAOPs were evaluated against the common features. The review offers an overview of how the qAOP concept has advanced and how it can aid toxicity assessment in the future. Further efforts are required to achieve validation, harmonisation and regulatory acceptance of qAOP models.

In silico toxicology: From structure–activity relationships towards deep learning and adverse outcome pathways

In silico toxicology is an emerging field. It gains increasing importance as research is aiming to decrease the use of animal experiments as suggested in the 3R principles by Russell and Burch. In silico toxicology is a means to identify hazards of compounds before synthesis, and thus in very early stages of drug development. For chemical industries, as well as regulatory agencies it can aid in gap‐filling and guide risk minimization strategies. Techniques such as structural alerts, read‐across, quantitative structure–activity relationship, machine learning, and deep learning allow to use in silico toxicology in many cases, some even when data is scarce. Especially the concept of adverse outcome pathways puts all techniques into a broader context and can elucidate predictions by mechanistic insights.

This article is categorized under:

Structure and Mechanism > Computational Biochemistry and Biophysics

Data Science > Chemoinformatics

Going All In: A Strategic Investment in In Silico Toxicology

As vast numbers of new chemicals are introduced to market annually, we are faced with the grand challenge of protecting humans and the environment while minimizing economically and ethically costly animal testing. In silico models promise to be the solution we seek, but we find ourselves at crossroads of future development efforts that would ensure standalone applicability and reliability of these tools. A conscientious effort that prioritizes experimental testing to support the needs of in silico models (versus regulatory needs) is called for to achieve this goal. Using economic analogy in the title of this work, we argue that a prudent investment is to go all-in to support in silico model development, rather than gamble our future by keeping the status quo of a "balanced portfolio" of testing approaches. We discuss two paths to future in silico toxicology-one based on big-data statistics ("broadsword"), and the other based on direct modeling of molecular interactions ("scalpel")-and offer rationale that the latter approach is more transparent, is better aligned with our quest for fundamental knowledge, and has a greater potential to succeed if we are willing to transform our toxicity-testing paradigm.

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.

In silico prediction of toxicity and its applications for chemicals at work

OBJECTIVE AND METHODS

This study reviewed the concept of in silico prediction of chemical toxicity for prevention of occupational cancer and future prospects in workers' health. In this review, a new approach to determine the credibility of in silico predictions with raw data is explored, and the method of determining the confidence level of evaluation based on the credibility of data is discussed. I searched various papers and books related to the in silico prediction of chemical toxicity and carcinogenicity. The intention was to utilize the most recent reports after 2015 regarding in silico prediction.

RESULTS AND CONCLUSION

The application of in silico methods is increasing with the prediction of toxic risks to human and the environment. The various toxic effects of industrial chemicals have triggered the recognition of the importance of using a combination of in silico models in the risk assessments. In silico occupational exposure models, industrial accidents, and occupational cancers are effectively managed and chemicals evaluated. It is important to identify and manage hazardous substances proactively through the rigorous evaluation of chemicals.

Advancing Computational Toxicology in the Big Data Era by Artificial Intelligence: Data-Driven and Mechanism-Driven Modeling for Chemical Toxicity

In 2016, the Frank R. Lautenberg Chemical Safety for the 21st Century Act became the first US legislation to advance chemical safety evaluations by utilizing novel testing approaches that reduce the testing of vertebrate animals. Central to this mission is the advancement of computational toxicology and artificial intelligence approaches to implementing innovative testing methods. In the current big data era, the terms volume (amount of data), velocity (growth of data), and variety (the diversity of sources) have been used to characterize the currently available chemical, in vitro, and in vivo data for toxicity modeling purposes. Furthermore, as suggested by various scientists, the variability (internal consistency or lack thereof) of publicly available data pools, such as PubChem, also presents significant computational challenges. The development of novel artificial intelligence approaches based on public massive toxicity data is urgently needed to generate new predictive models for chemical toxicity evaluations and make the developed models applicable as alternatives for evaluating untested compounds. In this procedure, traditional approaches (e.g., QSAR) purely based on chemical structures have been replaced by newly designed data-driven and mechanism-driven modeling. The resulting models realize the concept of adverse outcome pathway (AOP), which can not only directly evaluate toxicity potentials of new compounds, but also illustrate relevant toxicity mechanisms. The recent advancement of computational toxicology in the big data era has paved the road to future toxicity testing, which will significantly impact on the public health.

Nuclear Receptors Database Including Negative Data (NR-DBIND): A Database Dedicated to Nuclear Receptors Binding Data Including Negative Data and Pharmacological Profile

Nuclear receptors (NRs) are transcription factors that regulate gene expression in various physiological processes through their interactions with small hydrophobic molecules. They constitute an important class of targets for drugs and endocrine disruptors and are widely studied for both health and environment concerns. Since the integration of negative data can be critical for accurate modeling of ligand activity profiles, we manually collected and annotated NRs interaction data (positive and negative) through a sharp review of the corresponding literature. 15 116 positive and negative interactions data are provided for 28 NRs together with 593 PDB structures in the freely available Nuclear Receptors Database Including Negative Data ( http://nr-dbind.drugdesign.fr ). The NR-DBIND contains the most extensive information about interaction data on NRs, which should bring valuable information to chemists, biologists, pharmacologists and toxicologists.

A Mechanistic Framework for Integrating Chemical Structure and High-Throughput Screening Results to Improve Toxicity Predictions

Adverse Outcome Pathways (AOPs) establish a connection between a molecular initiating event (MIE) and an adverse outcome. Detailed understanding of the MIE provides the ideal data for determining chemical properties required to elicit the MIE. This study utilized high-throughput screening data from the ToxCast program, coupled with chemical structural information, to generate chemical clusters using three similarity methods pertaining to nine MIEs within an AOP network for hepatic steatosis. Three case studies demonstrate the utility of the mechanistic information held by the MIE for integrating biological and chemical data. Evaluation of the chemical clusters activating the glucocorticoid receptor identified activity differences in chemicals within a cluster. Comparison of the estrogen receptor results with previous work showed that bioactivity data and structural alerts can be combined to improve predictions in a customizable way where bioactivity data are limited. The aryl hydrocarbon receptor (AHR) highlighted that while structural data can be used to offset limited data for new screening efforts, not all ToxCast targets have sufficient data to define robust chemical clusters. In this context, an alternative to additional receptor assays is proposed where assays for proximal key events downstream of AHR activation could be used to enhance confidence in active calls. These case studies illustrate how the AOP framework can support an iterative process whereby in vitro toxicity testing and chemical structure can be combined to improve toxicity predictions. In vitro assays can inform the development of structural alerts linking chemical structure to toxicity. Consequently, structurally related chemical groups can facilitate identification of assays that would be informative for a specific MIE. Together, these activities form a virtuous cycle where the mechanistic basis for the in vitro results and the breadth of the structural alerts continually improve over time to better predict activity of chemicals for which limited toxicity data exist.

A Workflow for Identifying Metabolically Active Chemicals to Complement in vitro Toxicity Screening

The new paradigm of toxicity testing approaches involves rapid screening of thousands of chemicals across hundreds of biological targets through use of in vitro assays. Such assays may lead to false negatives when the complex metabolic processes that render a chemical bioactive in a living system are unable to be replicated in an in vitro environment. In the current study, a workflow is presented for complementing in vitro testing results with in silico and in vitro techniques to identify inactive parents that may produce active metabolites. A case study applying this workflow involved investigating the influence of metabolism for over 1,400 chemicals considered inactive across18 in vitro assays related to the estrogen receptor (ER) pathway. Over 7,500 first-generation and second-generation metabolites were generated for these in vitro inactive chemicals using an in silico software program. Next, a consensus model comprised of four individual quantitative structure activity relationship (QSAR) models was used to predict ER-binding activity for each of the metabolites. Binding activity was predicted for ~8-10% of metabolites in each generation, with these metabolites linked to 259 in vitro inactive parent chemicals. Metabolites were enriched in substructures consisting of alcohol, aromatic, and phenol bonds relative to their inactive parent chemicals, suggesting these features are potentially favorable for ER-binding. The workflow presented here can be used to identify parent chemicals that can be potentially bioactive, to aid confidence in high throughput risk screening.

Comparative Proteomic Analysis of Liver Steatosis and Fibrosis after Oral Hepatotoxicant Administration in Sprague-Dawley Rats

The past decade has seen an increase in the development and clinical use of biomarkers associated with histological features of liver disease. Here, we conduct a comparative histological and global proteomics analysis to identify coregulated modules of proteins in the progression of hepatic steatosis or fibrosis. We orally administered the reference chemicals bromobenzene (BB) or 4,4'-methylenedianiline (4,4'-MDA) to male Sprague-Dawley rats for either 1 single administration or 5 consecutive daily doses. Livers were preserved for histopathology and global proteomics assessment. Analysis of liver sections confirmed a dose- and time-dependent increase in frequency and severity of histopathological features indicative of lipid accumulation after BB or fibrosis after 4,4'-MDA. BB administration resulted in a dose-dependent increase in the frequency and severity of inflammation and vacuolation. 4,4'-MDA administration resulted in a dose-dependent increase in the frequency and severity of periportal collagen accumulation and inflammation. Pathway analysis identified a time-dependent enrichment of biological processes associated with steatogenic or fibrogenic initiating events, cellular functions, and toxicological states. Differentially expressed protein modules were consistent with the observed histology, placing physiologically linked protein networks into context of the disease process. This study demonstrates the potential for protein modules to provide mechanistic links between initiating events and histopathological outcomes.

An automated framework for QSAR model building

Background

In-silico quantitative structure–activity relationship (QSAR) models based tools are widely used to screen huge databases of compounds in order to determine the biological properties of chemical molecules based on their chemical structure. With the passage of time, the exponentially growing amount of synthesized and known chemicals data demands computationally efficient automated QSAR modeling tools, available to researchers that may lack extensive knowledge of machine learning modeling. Thus, a fully automated and advanced modeling platform can be an important addition to the QSAR community.

Results

In the presented workflow the process from data preparation to model building and validation has been completely automated. The most critical modeling tasks (data curation, data set characteristics evaluation, variable selection and validation) that largely influence the performance of QSAR models were focused. It is also included the ability to quickly evaluate the feasibility of a given data set to be modeled. The developed framework is tested on data sets of thirty different problems. The best-optimized feature selection methodology in the developed workflow is able to remove 62–99% of all redundant data. On average, about 19% of the prediction error was reduced by using feature selection producing an increase of 49% in the percentage of variance explained (PVE) compared to models without feature selection. Selecting only the models with a modelability score above 0.6, average PVE scores were 0.71. A strong correlation was verified between the modelability scores and the PVE of the models produced with variable selection.

Conclusions

We developed an extendable and highly customizable fully automated QSAR modeling framework. This designed workflow does not require any advanced parameterization nor depends on users decisions or expertise in machine learning/programming. With just a given target or problem, the workflow follows an unbiased standard protocol to develop reliable QSAR models by directly accessing online manually curated databases or by using private data sets. The other distinctive features of the workflow include prior estimation of data modelability to avoid time-consuming modeling trials for non modelable data sets, an efficient variable selection procedure and the facility of output availability at each modeling task for the diverse application and reproduction of historical predictions. The results reached on a selection of thirty QSAR problems suggest that the approach is capable of building reliable models even for challenging problems.

Electronic supplementary material

The online version of this article (10.1186/s13321-017-0256-5) contains supplementary material, which is available to authorized users.

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.

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.

Integrating Drug's Mode of Action into Quantitative Structure-Activity Relationships for Improved Prediction of Drug-Induced Liver Injury

Drug-induced liver injury (DILI) is complex in mechanism. Different drugs could undergo different mechanisms but result in the same DILI type, while the same drug could lead to different DILI types via different mechanisms. Therefore, predicting a drug's potential for DILI should take its underlying mechanisms into consideration. To achieve that, we constructed a novel approach by incorporating the drug's Mode of Action (MOA) into Quantitative Structure-Activity Relationship (QSAR) modeling. This MOA-DILI approach was examined using a data set of 333 drugs. The drugs were first grouped according to their MOA profiles (positive or negative in each MOA) based on the Tox21 qHTS assays. QSAR models for individual MOA assays were developed and subsequently combined to obtain the MOA-DILI model. A hold-out testing strategy (222 drugs for training and 111 drugs as a test set) was employed, which yielded a predictive accuracy of 0.711. The MOA-DILI model was directly compared with the standard QSAR approach using the same hold-out strategy, and the QSAR model yielded an accuracy of 0.662. To minimize the random chance in splitting training/test sets, the hold-out testing process was repeated 1000 times, and the observed difference in prediction accuracy between MOA-DILI and QSARs was statistically significant (P value <0.0001). Out of 17 MOAs used, four assays (i.e., antioxidant response elements, PPAR-gamma, estrogen receptor, and thyroid receptor assays) contributed most to the improved prediction of the MOA-DILI model over QSARs. In conclusion, the MOA-DILI approach has the potential to significantly improve predictive outcomes and to reveal complex relationships between MOAs and DILI, all of which would be helpful in developing DILI predictive models in drug screening and for risk assessment of industrial chemicals.

Systems Toxicology: Real World Applications and Opportunities

Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized from empirical end points to describing modes of action as adverse outcome pathways and perturbed networks. Toward this aim, Systems Toxicology entails the integration of in vitro and in vivo toxicity data with computational modeling. This evolving approach depends critically on data reliability and relevance, which in turn depends on the quality of experimental models and bioanalysis techniques used to generate toxicological data. Systems Toxicology involves the use of large-scale data streams ("big data"), such as those derived from omics measurements that require computational means for obtaining informative results. Thus, integrative analysis of multiple molecular measurements, particularly acquired by omics strategies, is a key approach in Systems Toxicology. In recent years, there have been significant advances centered on in vitro test systems and bioanalytical strategies, yet a frontier challenge concerns linking observed network perturbations to phenotypes, which will require understanding pathways and networks that give rise to adverse responses. This summary perspective from a 2016 Systems Toxicology meeting, an international conference held in the Alps of Switzerland, describes the limitations and opportunities of selected emerging applications in this rapidly advancing field. Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized, from empirical end points to pathways of toxicity. This requires the integration of in vitro and in vivo data with computational modeling. Test systems and bioanalytical technologies have made significant advances, but ensuring data reliability and relevance is an ongoing concern. The major challenge facing the new pathway approach is determining how to link observed network perturbations to phenotypic toxicity.

Generic and Automated Data Evaluation in Analytical Measurement

In the past year, automation has become more and more important in the field of elemental and structural chemical analysis to reduce the high degree of manual operation and processing time as well as human errors. Thus, a high number of data points are generated, which requires fast and automated data evaluation. To handle the preprocessed export data from different analytical devices with software from various vendors offering a standardized solution without any programming knowledge should be preferred. In modern laboratories, multiple users will use this software on multiple personal computers with different operating systems (e.g., Windows, Macintosh, Linux). Also, mobile devices such as smartphones and tablets have gained growing importance. The developed software, Project Analytical Data Evaluation (ADE), is implemented as a web application. To transmit the preevaluated data from the device software to the Project ADE, the exported XML report files are detected and the included data are imported into the entities database using the Data Upload software. Different calculation types of a sample within one measurement series (e.g., method validation) are identified using information tags inside the sample name. The results are presented in tables and diagrams on different information levels (general, detailed for one analyte or sample).