A problem formulation framework for the application of in silico toxicology methods in chemical risk assessment.

The first step in the hazard or risk assessment of chemicals should be to formulate the problem through a systematic and iterative process aimed at identifying and defining factors critical to the assessment. However, no general agreement exists on what components an in silico toxicology problem formulation (PF) should include. The present work aims to develop a PF framework relevant to the application of in silico models for chemical toxicity prediction. We modified and applied a PF framework from the general risk assessment literature to peer reviewed papers describing PFs associated with in silico toxicology models. Important gaps between the general risk assessment literature and the analyzed PF literature associated with in silico toxicology methods were identified. While the former emphasizes the need for PFs to address higher-level conceptual questions, the latter does not. There is also little consistency in the latter regarding the PF components addressed, reinforcing the need for a PF framework that enable users of in silico toxicology models to answer the central conceptual questions aimed at defining components critical to the model application. Using the developed framework, we highlight potential areas of uncertainty manifestation in in silico toxicology PF in instances where particular components are missing or implicitly described. The framework represents the next step in standardizing in silico toxicology PF component. The framework can also be used to improve the understanding of how uncertainty is apparent in an in silico toxicology PF, thus facilitating ways to address uncertainty.

The evaluation of skin sensitization potential of the UVCB substance diisopentyl phthalate by in silico and in vitro methods.

Diisopentyl phthalate (DiPeP) is primarily used as a plasticizer or additive within the production of polyvinyl chloride (PVC), and has many additional industrial applications. Its metabolites were recently found in urinary samples of pregnant women; thus, this substance is of concern as relates to human exposure. Depending upon the nature of the alcohol used in its synthesis, DiPeP may exist either as a mixture consisting of several branched positional isomers, or as a single defined structure. This article investigates the skin sensitization potential and immunomodulatory effects of DiPeP CAS No. 84777-06-0, which is currently marketed and classified as a UVCB substance, by in silico and in vitro methods. Our findings showed an immunomodulatory effect for DiPeP in LPS-induced THP-1 activation assay (increased CD54 expression). In silico predictions using QSAR TOOLBOX 4.5, ToxTree, and VEGA did not identify DiPeP, in the form of a discrete compound, as a skin sensitizer. The keratinocyte activation (Key Event 2 (KE2) of the adverse outcome pathway (AOP) for skin sensitization) was evaluated by two different test methods (HaCaT assay and RHE assay), and results were discordant. While the HaCaT assay showed that DiPeP can activate keratinocytes (increased levels of IL-6, IL-8, IL-1α, and ILA gene expression), in the RHE assay, DiPeP slightly increased IL-6 release. Although inconclusive for KE2, the role of DiPeP in KE3 (dendritic cell activation) was demonstrated by the increased levels of CD54 and IL-8 and TNF-α in THP-1 cells (THP-1 activation assay). Altogether, findings were inconclusive regarding the skin sensitization potential of the UVCB DiPeP-disagreeing with the results of DiPeP in the form of discrete compound (skin sensitizer by the LLNA assay). Additional studies are needed to elucidate the differences between DiPeP isomer forms, and to better understand the applicability domains of non-animal methods in identifying skin sensitization hazards of UVCB substances.

Making in silico predictive models for toxicology FAIR.

In silico predictive models for toxicology include quantitative structure-activity relationship (QSAR) and physiologically based kinetic (PBK) approaches to predict physico-chemical and ADME properties, toxicological effects and internal exposure. Such models are used to fill data gaps as part of chemical risk assessment. There is a growing need to ensure in silico predictive models for toxicology are available for use and that they are reproducible. This paper describes how the FAIR (Findable, Accessible, Interoperable, Reusable) principles, developed for data sharing, have been applied to in silico predictive models. In particular, this investigation has focussed on how the FAIR principles could be applied to improved regulatory acceptance of predictions from such models. Eighteen principles have been developed that cover all aspects of FAIR. It is intended that FAIRification of in silico predictive models for toxicology will increase their use and acceptance.

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.

The use of Bayesian methodology in the development and validation of a tiered assessment approach towards prediction of rat acute oral toxicity.

There exists consensus that the traditional means by which safety of chemicals is assessed-namely through reliance upon apical outcomes obtained following in vivo testing-is increasingly unfit for purpose. Whilst efforts in development of suitable alternatives continue, few have achieved levels of robustness required for regulatory acceptance. An array of "new approach methodologies" (NAM) for determining toxic effect, spanning in vitro and in silico spheres, have by now emerged. It has been suggested, intuitively, that combining data obtained from across these sources might serve to enhance overall confidence in derived judgment. This concept may be formalised in the "tiered assessment" approach, whereby evidence gathered through a sequential NAM testing strategy is exploited so to infer the properties of a compound of interest. Our intention has been to provide an illustration of how such a scheme might be developed and applied within a practical setting-adopting for this purpose the endpoint of rat acute oral lethality. Bayesian statistical inference is drawn upon to enable quantification of degree of confidence that a substance might ultimately belong to one of five LD50-associated toxicity categories. Informing this is evidence acquired both from existing in silico and in vitro resources, alongside a purposely-constructed random forest model and structural alert set. Results indicate that the combination of in silico methodologies provides moderately conservative estimations of hazard, conducive for application in safety assessment, and for which levels of certainty are defined. Accordingly, scope for potential extension of approach to further toxicological endpoints is demonstrated.

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.

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.

Exploring the Potential of ToxCast Data in Supporting Read-Across for Evaluation of Food Chemical Safety.

The intention of this study was to determine the utility of high-throughput screening (HTS) data, as exemplified by ToxCast and Tox21, for application in toxicological read-across in food-relevant chemicals. Key questions were addressed on the extent to which the HTS data could provide information enabling (1) the elucidation of underlying bioactivities associated with apical toxicological outcomes, (2) the closing of existing toxicological data gaps, and (3) the definition of the boundaries of chemical space across which bioactivity could reliably be extrapolated. Results revealed that many biological targets apparently activated within the chemical groupings lack, at this time, validated toxicity pathway associations. Therefore, as means of providing proof-of-principle, a comparatively well-characterized end point-estrogenicity-was selected for evaluation. This was facilitated through the preparation of two exploratory case studies, focusing upon groupings of paraben-gallates and pyranone-type compounds (notably flavonoids). Within both, the HTS data were seen to reflect estrogenic potencies in a manner which broadly corresponded to established structure-activity group relationships, with parabens and flavonoids displaying greater estrogen receptor affinity than benzoate esters and alternative pyranone-containing molecules, respectively. As such, utility in the identification of out-of-domain compounds was demonstrated, indicating potential for application in addressing point (3) as detailed above.

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.

Incorporating lines of evidence from New Approach Methodologies (NAMs) to reduce uncertainties in a category based read-across: A case study for repeated dose toxicity.

A group of triazole compounds was selected to investigate the confidence that may be associated with read-across of a complex data gap: repeated dose toxicity. The read-across was evaluated using Assessment Elements (AEs) from the European Chemicals Agency's (ECHA's) Read-Across Assessment Framework (RAAF), alongside appraisal of associated uncertainties. Following an initial read-across based on chemical structure and properties, uncertainties were reduced by the integration of data streams such as those from New Approach Methodologies (NAM) and other existing data. In addition, addressing the findings of the ECHA RAAF framework, complemented with specific questions concerning uncertainties, increased the confidence that can be placed in read-across. Although a data rich group of compounds with a strong mechanistic basis was analysed, it was clearly demonstrated that NAM data available from publicly available resources could be applied to support read-across. It is acknowledged that most read-across studies will not be so data rich or mechanistically robust, therefore some targeted experimentation may be required to fill the data gaps. In this sense, NAMs should constitute new experimental tests performed with the specific goal of reducing the uncertainties and demonstrating the read-across hypothesis.

Determination of "fitness-for-purpose" of quantitative structure-activity relationship (QSAR) models to predict (eco-)toxicological endpoints for regulatory use.

In silico models are used to predict toxicity and molecular properties in chemical safety assessment, gaining widespread regulatory use under a number of legislations globally. This study has rationalised previously published criteria to evaluate quantitative structure-activity relationships (QSARs) in terms of their uncertainty, variability and potential areas of bias, into ten assessment components, or higher level groupings. The components have been mapped onto specific regulatory uses (i.e. data gap filling for risk assessment, classification and labelling, and screening and prioritisation) identifying different levels of uncertainty that may be acceptable for each. Twelve published QSARs were evaluated using the components, such that their potential use could be identified. High uncertainty was commonly observed with the presentation of data, mechanistic interpretability, incorporation of toxicokinetics and the relevance of the data for regulatory purposes. The assessment components help to guide strategies that can be implemented to improve acceptability of QSARs through the reduction of uncertainties. It is anticipated that model developers could apply the assessment components from the model design phase (e.g. through problem formulation) through to their documentation and use. The application of the components provides the possibility to assess QSARs in a meaningful manner and demonstrate their fitness-for-purpose against pre-defined criteria.

A Systematic Review of Published Physiologically-based Kinetic Models and an Assessment of their Chemical Space Coverage.

Across multiple sectors, including food, cosmetics and pharmaceutical industries, there is a need to predict the potential effects of xenobiotics. These effects are determined by the intrinsic ability of the substance, or its derivatives, to interact with the biological system, and its concentration-time profile at the target site. Physiologically-based kinetic (PBK) models can predict organ-level concentration-time profiles, however, the models are time and resource intensive to generate de novo. Read-across is an approach used to reduce or replace animal testing, wherein information from a data-rich chemical is used to make predictions for a data-poor chemical. The recent increase in published PBK models presents the opportunity to use a read-across approach for PBK modelling, that is, to use PBK model information from one chemical to inform the development or evaluation of a PBK model for a similar chemical. Essential to this process, is identifying the chemicals for which a PBK model already exists. Herein, the results of a systematic review of existing PBK models, compliant with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) format, are presented. Model information, including species, sex, life-stage, route of administration, software platform used and the availability of model equations, was captured for 7541 PBK models. Chemical information (identifiers and physico-chemical properties) has also been recorded for 1150 unique chemicals associated with these models. This PBK model data set has been made readily accessible, as a Microsoft Excel® spreadsheet, providing a valuable resource for those developing, using or evaluating PBK models in industry, academia and the regulatory sectors.

Guidance for good practice in the application of machine learning in development of toxicological quantitative structure-activity relationships (QSARs).

Recent years have seen a substantial growth in the adoption of machine learning approaches for the purposes of quantitative structure-activity relationship (QSAR) development. Such a trend has coincided with desire to see a shifting in the focus of methodology employed within chemical safety assessment: away from traditional reliance upon animal-intensive in vivo protocols, and towards increased application of in silico (or computational) predictive toxicology. With QSAR central amongst techniques applied in this area, the emergence of algorithms trained through machine learning with the objective of toxicity estimation has, quite naturally, arisen. On account of the pattern-recognition capabilities of the underlying methods, the statistical power of the ensuing models is potentially considerable-appropriate for the handling even of vast, heterogeneous datasets. However, such potency comes at a price: this manifesting as the general practical deficits observed with respect to the reproducibility, interpretability and generalisability of the resulting tools. Unsurprisingly, these elements have served to hinder broader uptake (most notably within a regulatory setting). Areas of uncertainty liable to accompany (and hence detract from applicability of) toxicological QSAR have previously been highlighted, accompanied by the forwarding of suggestions for "best practice" aimed at mitigation of their influence. However, the scope of such exercises has remained limited to "classical" QSAR-that conducted through use of linear regression and related techniques, with the adoption of comparatively few features or descriptors. Accordingly, the intention of this study has been to extend the remit of best practice guidance, so as to address concerns specific to employment of machine learning within the field. In doing so, the impact of strategies aimed at enhancing the transparency (feature importance, feature reduction), generalisability (cross-validation) and predictive power (hyperparameter optimisation) of algorithms, trained upon real toxicity data through six common learning approaches, is evaluated.

Update of the Cancer Potency Database (CPDB) to enable derivations of Thresholds of Toxicological Concern (TTC) for cancer potency.

The purpose of this study was to update the existing Cancer Potency Database (CPDB) in order to support the development of a dataset of compounds, with associated points of departure (PoDs), to enable a review and update of currently applied values for the Threshold of Toxicological Concern (TTC) for cancer endpoints. This update of the current CPDB, last reviewed in 2012, includes the addition of new data (44 compounds and 158 studies leading to additional 359 dose-response curves). Strict inclusion criteria were established and applied to select compounds and studies with relevant cancer potency data. PoDs were calculated from dose-response modeling, including the benchmark dose (BMD) and the lower 90% confidence limits (BMDL) at a specified benchmark response (BMR) of 10%. The updated full CPDB database resulted in a total of 421 chemicals which had dose-response data that could be used to calculate PoDs. This candidate dataset for cancer TTC is provided in a transparent and adaptable format for further analysis of TTC to derive cancer potency thresholds.

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.

In Silico Toxicology Data Resources to Support Read-Across and (Q)SAR.

A plethora of databases exist online that can assist in in silico chemical or drug safety assessment. However, a systematic review and grouping of databases, based on purpose and information content, consolidated in a single source, has been lacking. To resolve this issue, this review provides a comprehensive listing of the key in silico data resources relevant to: chemical identity and properties, drug action, toxicology (including nano-material toxicity), exposure, omics, pathways, Absorption, Distribution, Metabolism and Elimination (ADME) properties, clinical trials, pharmacovigilance, patents-related databases, biological (genes, enzymes, proteins, other macromolecules etc.) databases, protein-protein interactions (PPIs), environmental exposure related, and finally databases relating to animal alternatives in support of 3Rs policies. More than nine hundred databases were identified and reviewed against criteria relating to accessibility, data coverage, interoperability or application programming interface (API), appropriate identifiers, types of in vitro, in vivo,-clinical or other data recorded and suitability for modelling, read-across, or similarity searching. This review also specifically addresses the need for solutions for mapping and integration of databases into a common platform for better translatability of preclinical data to clinical data.

In vitro Liver Zonation of Primary Rat Hepatocytes.

The ability of the liver to simultaneously carry out multiple functions is dependent on the metabolic heterogeneity of hepatocytes spatially located within a liver lobule spanning from the portal triad to the central vein. This complex zonal architecture of the liver, however, makes accurate in vitro modeling a challenge and often standard culture systems assume a homogenous model which may lead to inaccurate translatability of results. Here, we use a combination of mathematical modeling and experimental data to demonstrate a readily constructible in vitro flow system capable of liver zonation in primary rat hepatocytes. We show the differential expression of zonation markers, enhanced functionality when compared to standard static cultures and zone-specific metabolism and cell damage in the presence of paracetamol, a known zone-specific toxin. This type of advanced system provides a more in-depth and essential understanding of liver physiology and pathophysiology as well as the accurate evaluation of pharmacological interventions at a zone-specific level.

Chemoinformatic Consideration of Novel Psychoactive Substances: Compilation and Preliminary Analysis of a Categorised Dataset.

Recent years have seen the emergence into circulation of a growing array of novel psychoactive substances (NPS). Knowledge of the pharmacological profiles and risk liability of these compounds is typically very scarce. Development of chemoinformatic tools enabling prediction of properties within uncharacterised analogues has potential be of particular use. In order to facilitate this, compilation of a chemical inventory comprising known NPS is a necessity. Sourcing a variety of published governmental and analytical reports, a dataset composed of 690 distinct acknowledged NPS, complete with defined chemical structures, has been constructed. This is supplemented by a complementary series of 155 established psychoactive drugs of abuse (EPDA). Classification was performed in accordance with their key molecular structural features, subjective effect profiles and pharmacological mechanisms of action. In excess of forty chemical groupings, spanning seven subjective effect categories and six broad mechanisms of pharmacological action, were identified. Co-occurrence of NPS and EPDA within specific classes was common, showcasing inherent scope both for chemical read-across and for the derivation of structural alerts.

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

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