Integrating QSAR models predicting acute contact toxicity and mode of action profiling in honey bees (A. mellifera): Data curation using open source databases, performance testing and validation

A QSAR study for predicting malformation in zebrafish embryo

BACKGROUND

Developmental toxicity tests are extremely expensive, require a large number of animals, and are time-consuming. It is necessary to develop a new approach to simplify the analysis of developmental endpoints. One of these endpoints is malformation, and one group of ongoing methods for simplifying is in silico models. In this study, we aim to develop a quantitative structure-activity relationship (QSAR) model and identify the best algorithm for predicting malformations, as well as the most important and effective physicochemical properties associated with malformation.

METHODS

The dataset was extracted from a reliable database called COMPTOX. Physicochemical properties (descriptors) were calculated using Mordred and RDKit chemoinformatics software. The data were cleaned, preprocessed, and then split into training and testing sets. Machine learning algorithms, such as gradient boosting model (GBM) and logistic regression (LR), as well as deep learning models, including multilayer perceptron (MLP) and neural networks (NNs) trained with train set data and different sets of descriptors. The models were then validated with test set and various statistical parameters, such as Matthew's correlation coefficient (MCC) and balanced accuracy (BAC) score, were used to compare the models.

RESULTS

A set of descriptors containing with 78% AUC was identified as the best set of descriptors. Gradient boosting was determined to be the best algorithm with 78% predictive power.

CONCLUSIONS

The descriptors that were the most effective for developing models directly impact the mechanism of malformation, and GBM is the best model due to its MCC and BAC.

Evaluation of the ecological risk of pesticide residues from the European LUCAS Soil monitoring 2018 survey

The 2018 LUCAS (Land Use and Coverage Area frame Survey) Soil Pesticides survey provides a European Union (EU)-scale assessment of 118 pesticide residues in more than 3473 soil sites. This study responds to the policy need to develop risk-based indicators for pesticides in the environment. Two mixture risk indicators are presented for soil based, respectively, on the lowest and the median of available No Observed Effect Concentration (NOECsoil,min and NOECsoil,50) from publicly available toxicity datasets. Two further indicators were developed based on the corresponding equilibrium concentration in the aqueous phase and aquatic toxicity data, which are available as species sensitivity distributions. Pesticides were quantified in 74.5% of the sites. The mixture risk indicator based on the NOECsoil,min exceeds 1 in 14% of the sites and 0.1 in 23%. The insecticides imidacloprid and chlorpyrifos and the fungicide epoxiconazole are the largest contributors to the overall risk. At each site, one or a few substances drive mixture risk. Modes of actions most likely associated with mixture effects include modulation of acetylcholine metabolism (neonicotinoids and organophosphate substances) and sterol biosynthesis inhibition (triazole fungicides). Several pesticides driving the risk have been phased out since 2018. Following LUCAS surveys will determine the effectiveness of substance-specific risk management and the overall progress toward risk reduction targets established by EU and UN policies. Newly generated data and knowledge will stimulate needed future research on pesticides, soil health, and biodiversity protection. Integr Environ Assess Manag 2024;20:1639-1653. © 2024 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Advancing ecotoxicity assessment: Leveraging pre-trained model for bee toxicity and compound degradability prediction

The prediction of ecological toxicity plays an increasingly important role in modern society. However, the existing models often suffer from poor performance and limited predictive capabilities. In this study, we propose a novel approach for ecological toxicity assessment based on pre-trained models. By leveraging pre-training techniques and graph neural network models, we establish a highperformance predictive model. Furthermore, we incorporate a variational autoencoder to optimize the model, enabling simultaneous discrimination of toxicity to bees and molecular degradability. Additionally, despite the low similarity between the endogenous hormones in bees and the compounds in our dataset, our model confidently predicts that these hormones are non-toxic to bees, which further strengthens the credibility and accuracy of our model. We also discovered the negative correlation between the degradation and bee toxicity of compounds. In summary, this study presents an ecological toxicity assessment model with outstanding performance. The proposed model accurately predicts the toxicity of chemicals to bees and their degradability capabilities, offering valuable technical support to relevant fields.

Multi-Endpoint Acute Toxicity Assessment of Organic Compounds Using Large-Scale Machine Learning Modeling

In recent years, alternative animal testing methods such as computational and machine learning approaches have become increasingly crucial for toxicity testing. However, the complexity and scarcity of available biomedical data challenge the development of predictive models. Combining nonlinear machine learning together with multicondition descriptors offers a solution for using data from various assays to create a robust model. This work applies multicondition descriptors (MCDs) to develop a QSTR (Quantitative Structure-Toxicity Relationship) model based on a large toxicity data set comprising more than 80,000 compounds and 59 different end points (122,572 data points). The prediction capabilities of developed single-task multi-end point machine learning models as well as a novel data analysis approach with the use of Convolutional Neural Networks (CNN) are discussed. The results show that using MCDs significantly improves the model and using them with CNN-1D yields the best result (R2train = 0.93, R2ext = 0.70). Several structural features showed a high level of contribution to the toxicity, including van der Waals surface area (VSA), number of nitrogen-containing fragments (nN+), presence of S-P fragments, ionization potential, and presence of C-N fragments. The developed models can be very useful tools to predict the toxicity of various compounds under different conditions, enabling quick toxicity assessment of new compounds.

Systematic approaches to machine learning models for predicting pesticide toxicity

Pesticides play an important role in modern agriculture by protecting crops from pests and diseases. However, the negative consequences of pesticides, such as environmental contamination and adverse effects on human and ecological health, underscore the importance of accurate toxicity predictions. To address this issue, artificial intelligence models have emerged as valuable methods for predicting the toxicity of organic compounds. In this review article, we explore the application of machine learning (ML) for pesticide toxicity prediction. This review provides a detailed summary of recent developments, prediction models, and datasets used for pesticide toxicity prediction. In this analysis, we compared the results of several algorithms that predict the harmfulness of various classes of pesticides. Furthermore, this review article identified emerging trends and areas for future direction, showcasing the transformative potential of machine learning in promoting safer pesticide usage and sustainable agriculture.

TKPlate 1.0: An Open-access platform for toxicokinetic and toxicodynamic modelling of chemicals to implement new approach methodologies in chemical risk assessment

This publication is linked to the following EFSA Supporting Publications articles: http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8441/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8440/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8437/full.

Complement System Inhibitory Drugs in a Zebrafish (Danio rerio) Model: Computational Modeling

The dysregulation of complement system activation usually results in acute or chronic inflammation and can contribute to the development of various diseases. Although the activation of complement pathways is essential for innate defense, exacerbated activity of this system may be harmful to the host. Thus, drugs with the potential to inhibit the activation of the complement system may be important tools in therapy for diseases associated with complement system activation. The synthetic peptides Cp40 and PMX205 can be highlighted in this regard, given that they selectively inhibit the C3 and block the C5a receptor (C5aR1), respectively. The zebrafish (Danio rerio) is a robust model for studying the complement system. The aim of the present study was to use in silico computational modeling to investigate the hypothesis that these complement system inhibitor peptides interact with their target molecules in zebrafish, for subsequent in vivo validation. For this, we analyzed molecular docking interactions between peptides and target molecules. Our study demonstrated that Cp40 and the cyclic peptide PMX205 have positive interactions with their respective zebrafish targets, thus suggesting that zebrafish can be used as an animal model for therapeutic studies on these inhibitors.

Machine learning-assisted data filtering and QSAR models for prediction of chemical acute toxicity on rat and mouse

Machine learning (ML) methods provide a new opportunity to build quantitative structure-activity relationship (QSAR) models for predicting chemicals' toxicity based on large toxicity data sets, but they are limited in insufficient model robustness due to poor data set quality for chemicals with certain structures. To address this issue and improve model robustness, we built a large data set on rat oral acute toxicity for thousands of chemicals, then used ML to filter chemicals favorable for regression models (CFRM). In comparison to chemicals not favorable for regression models (CNRM), CFRM accounted for 67% of chemicals in the original data set, and had a higher structural similarity and a smaller toxicity distribution in 2-4 log₁₀ (mg/kg). The performance of established regression models for CFRM was greatly improved, with root-mean-square deviations (RMSE) in the range of 0.45-0.48 log₁₀ (mg/kg). Classification models were built for CNRM using all chemicals in the original data set, and the area under receiver operating characteristic (AUROC) reached 0.75-0.76. The proposed strategy was successfully applied to a mouse oral acute data set, yielding RMSE and AUROC in the range of 0.36-0.38 log₁₀ (mg/kg) and 0.79, respectively.

Risk Assessment of Combined Exposure to Multiple Chemicals at the European Food Safety Authority: Principles, Guidance Documents, Applications and Future Challenges

Human health and animal health risk assessment of combined exposure to multiple chemicals use the same steps as single-substance risk assessment, namely problem formulation, exposure assessment, hazard assessment and risk characterisation. The main unique feature of combined RA is the assessment of combined exposure, toxicity and risk. Recently, the Scientific Committee of the European Food Safety Authority (EFSA) published two relevant guidance documents. The first one "Harmonised methodologies for the human health, animal health and ecological risk assessment of combined exposure to multiple chemicals" provides principles and explores methodologies for all steps of risk assessment together with a reporting table. This guidance supports also the default assumption that dose addition is applied for combined toxicity of the chemicals unless evidence for response addition or interactions (antagonism or synergism) is available. The second guidance document provides an account of the scientific criteria to group chemicals in assessment groups using hazard-driven criteria and prioritisation methods, i.e., exposure-driven and risk-based approaches. This manuscript describes such principles, provides a brief description of EFSA's guidance documents, examples of applications in the human health and animal health area and concludes with a discussion on future challenges in this field.

Toxicity Rank Order (TRO) As a New Approach for Toxicity Prediction by QSAR Models

Quantitative Structure-Activity Relationship (QSAR) models are commonly used for risk assessment of emerging contaminants. The objective of this study was to use a toxicity rank order (TRO) as an integrating parameter to improve the toxicity prediction by QSAR models. TRO for each contaminant was calculated from collected toxicity data including acute toxicity concentration and no observed effect concentration. TRO values associated with toxicity mechanisms were used to classify pollutants into three modes of action consisting of narcosis, transition and reactivity. The selection principle of parameters for QSAR models was established and verified. It showed a reasonable prediction of toxicities caused by organophosphates and benzene derivatives, especially. Compared with traditional procedures, incorporating TRO showed an improved correlation coefficient of QSAR models by approximately 10%. Our study indicated that the proposed procedure can be used for screening modeling parameter data and improve the toxicity prediction by QSAR models, and this could facilitate prediction and evaluation of environmental contaminant toxicity.

Lethal, sublethal, and combined effects of pesticides on bees: A meta-analysis and new risk assessment tools

Multiple stressors threaten bee health, a major one being pesticides. Bees are simultaneously exposed to multiple pesticides that can cause both lethal and sublethal effects. Risk assessment and most research on bee health, however, focus on lethal individual effects. Here, we performed a systematic literature review and meta-analysis that summarizes and re-interprets the available qualitative and quantitative information on the lethal, sublethal, and combined toxicity of a comprehensive range of pesticides on bees. We provide results (1970-2019) for multiple bee species (Bombus, Osmia, Megachile, Melipona, Partamona, Scaptotrigona), although most works focused on Apis mellifera L. (78 %). Our harmonised results document the lethal toxicity of pesticides in bees (n = 377 pesticides) and the types of sublethal testing methods and related effects that cause a sublethal effect (n = 375 sublethal experiments). We identified the most common combinations of pesticides and mode of actions tested, and summarize the experimental methods, magnitude of the interactions, and robustness of available data (n = 361 experiments). We provide open access searchable, comprehensive, and integrated list of pesticides and their levels causing lethal, sublethal, and combined effects. We report major data gaps related to pesticide's sublethal (71 %) and combined (e.g., ~99 %) toxicity. We identified pesticides and mode of actions of greatest concern in terms of sublethal (chlorothalonil, pymetrozine, glyphosate; neonicotinoids) and combined (tau-fluvalinate combinations; acetylcholinesterase inhibitors and neonicotinoids) effects. Although certain pesticides have faced regulatory restrictions in specific countries (chlorothalonil, pymetrozine, neonicotinoids), most are still widely used worldwide (e.g., glyphosate). This work aims at facilitating the implementation of more comprehensive and harmonised research and risk assessments, considering sublethal and combined effects. To ensure safeguarding pollinators and the environment, we advocate for a more refined and holistic assessment that do not only focus on lethality but uses harmonised methods to test sublethal and relevant combinations.

A regression-based QSAR-model to predict acute toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica): Calibration, validation, and future developments to support risk assessment of chemicals in amphibians

Amphibian populations are undergoing a global decline worldwide. Such decline has been attributed to their unique physiology, ecology, and exposure to multiple stressors including chemicals, temperature, and biological hazards such as fungi of the Batrachochytrium genus, viruses such as Ranavirus, and habitat reduction. There are limited toxicity data for chemicals available for amphibians and few quantitative structure-activity relationship (QSAR) models have been developed and are publicly available. Such QSARs provide important tools to assess the toxicity of chemicals particularly in a data poor context. QSARs provide important tools to assess the toxicity of chemicals particularly when no toxicological data are available. This manuscript provides a description and validation of a regression-based QSAR model to predict, in a quantitative manner, acute lethal toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica). QSAR models for acute median lethal molar concentrations (LC50-12 h) of waterborne chemicals using the Monte Carlo method were developed. The statistical characteristics of the QSARs were described as average values obtained from five random distributions into training and validation sets. Predictions from the model gave satisfactory results for the overall training set (R² = 0.72 and RMSE = 0.33) and were even more robust for the validation set (R² = 0.96 and RMSE = 0.11). Further development of QSAR models in amphibians, particularly for other life stages and species, are discussed.

A critical review of effect modeling for ecological risk assessment of plant protection products

A wide diversity of plant protection products (PPP) is used for crop protection leading to the contamination of soil, water, and air, which can have ecotoxicological impacts on living organisms. It is inconceivable to study the effects of each compound on each species from each compartment, experimental studies being time consuming and cost prohibitive, and animal testing having to be avoided. Therefore, numerous models are developed to assess PPP ecotoxicological effects. Our objective was to provide an overview of the modeling approaches enabling the assessment of PPP effects (including biopesticides) on the biota. Six categories of models were inventoried: (Q)SAR, DR and TKTD, population, multi-species, landscape, and mixture models. They were developed for various species (terrestrial and aquatic vertebrates and invertebrates, primary producers, micro-organisms) belonging to diverse environmental compartments, to address different goals (e.g., species sensitivity or PPP bioaccumulation assessment, ecosystem services protection). Among them, mechanistic models are increasingly recognized by EFSA for PPP regulatory risk assessment but, to date, remain not considered in notified guidance documents. The strengths and limits of the reviewed models are discussed together with improvement avenues (multigenerational effects, multiple biotic and abiotic stressors). This review also underlines a lack of model testing by means of field data and of sensitivity and uncertainty analyses. Accurate and robust modeling of PPP effects and other stressors on living organisms, from their application in the field to their functional consequences on the ecosystems at different scales of time and space, would help going toward a more sustainable management of the environment. Graphical Abstract Combination of the keyword lists composing the first bibliographic query. Columns were joined together with the logical operator AND. All keyword lists are available in Supplementary Information at https://doi.org/10.5281/zenodo.5775038 (Larras et al. 2021).

Chemometric modeling of plant protection products (PPPs) for the prediction of acute contact toxicity against honey bees (A. mellifera): A 2D-QSAR approach

Honey bees (Apis mellifera) are vital for economic, viable agriculture and for food safety. Although Plant Protection Products (PPPs) are of undeniable importance in the global agricultural system, these have become potential threats for non-target organisms like pollinators (e.g., honey bees etc.), resulting in the disruption of the ecological balance. In the current work, we have used the 113 PPP analogs to develop a 2D-QSAR model and explored the structural features modulating the toxic effects on honey bees, following the Organization for Economic Co-operation and Development (OECD) guidelines. The extensive validation of the developed model has been performed using internal and external validation metrics to make sure that the model is statistically sound and interpretable enough to be acceptable. The obtained results (R² = 0.666, Q² = 0.594, Q²_F1 = 0.647 and Q²_F2 = 0.646) determine the predictability and reliability of the developed model. This model should be useful for the predictions (acute contact toxicity (LD₅₀)) of the new and untested compounds located inside the applicability domain of the developed model. Moreover, we have performed the in-silico prediction of toxicity against honey bees of a total of 709 compounds obtained from the pesticide properties database (PPDB) using the developed model.

Ecological risk assessment for difenoconazole in aquatic ecosystems using a web-based interspecies correlation estimation (ICE)-species sensitivity distribution (SSD) model

Difenoconazole is a typical triazole fungicide that can inhibit demethylation during ergosterol synthesis. Due to its wide use, difenoconazole is frequently detected in surface water, paddy water, agricultural water, and other aquatic environments. Presently, an assessment of the ecological risk posed by difenoconazole in aquatic ecosystems is lacking. Here, a web-based interspecies correlation estimation (ICE)-species sensitivity distribution (SSD) model was first applied to assess the ecological risk of difenoconazole in aquatic environments. Meanwhile, maximum acceptable concentration (MAC), maximum risk-free concentration (MRFC), and risk quotient (RQ) values were used to evaluate the potential risk of difenoconazole to aquatic organisms. Our results showed that an aquatic MAC value of 0.31 μg/L was acceptable for difenoconazole in aquatic environments. Further, the detected concentration of difenoconazole was lower than the MRFC value of 0.09 μg/L indicating no risk to aquatic organisms. Assessment data suggested that difenoconazole exhibited potential risks to eight studied aquatic ecosystems (including surface water, paddy water, and agricultural water) in different countries (RQ > 1), indicating that difenoconazole overuse could cause adverse effects to aquatic organisms in these aquatic ecosystems. Thus, restricted use and rational use of difenoconazole are recommended.

In Silico Methods for Environmental Risk Assessment: Principles, Tiered Approaches, Applications, and Future Perspectives

This chapter aims to introduce the reader to the basic principles of environmental risk assessment of chemicals and highlights the usefulness of tiered approaches within weight of evidence approaches in relation to problem formulation i.e., data availability, time and resource availability. In silico models are then introduced and include quantitative structure-activity relationship (QSAR) models, which support filling data gaps when no chemical property or ecotoxicological data are available. In addition, biologically-based models can be applied in more data rich situations and these include generic or species-specific models such as toxicokinetic-toxicodynamic models, dynamic energy budget models, physiologically based models, and models for ecosystem hazard assessment i.e. species sensitivity distributions and ultimately for landscape assessment i.e. landscape-based modeling approaches. Throughout this chapter, particular attention is given to provide practical examples supporting the application of such in silico models in real-world settings. Future perspectives are discussed to address environmental risk assessment in a more holistic manner particularly for relevant complex questions, such as the risk assessment of multiple stressors and the development of harmonized approaches to ultimately quantify the relative contribution and impact of single chemicals, multiple chemicals and multiple stressors on living organisms.

Binary and ternary toxicological interactions of clothianidin and eight commonly used pesticides on honey bees (Apis mellifera)

Although many toxicological evaluations have been conducted for honey bees (Apis mellifera), most of these studies have only focused on the effects of individual chemicals. However, honey bees are usually exposed to pesticide mixtures under field conditions. In this study, we examined the effects of individual pesticides and mixtures of clothianidin (CLO) with eight other pesticides [carbaryl (CAR), thiodicarb (THI), chlorpyrifos (CHL), beta-cyfluthrin (BCY), gamma-cyhalothrin (GCY), tetraconazole (TET), spinosad (SPI) and indoxacarb (IND)] on honey bees using a feeding method. Toxicity tests of a 4-day exposure to individual pesticides revealed that CLO had the highest toxicity to A. mellifera, with an LC₅₀ value of 0.24 μg a.i. mL^-1, followed by IND and CHL with LC₅₀ values of 3.40 and 3.56 μg a.i. mL^-1, respectively. SPI and CAR had relatively low toxicities, with LC₅₀ values of 7.19 and 8.42 μg a.i. mL^-1, respectively. In contrast, TET exhibited the least toxicity, with an LC₅₀ value of 258.7 μg a.i. mL^-1. Most binary mixtures of CLO with other pesticides exerted additive and antagonistic effects. However, all the ternary mixtures containing CLO and TET (except for CLO+TET+THD) elicited synergistic responses to bees. Either increased numbers of components in the mixture or/and a unique mode of action appeared to be responsible for the higher toxicity of mixtures. Our findings emphasized the need for risk assessment of pesticide mixtures rather than the individual chemicals. Our data also provided information that might help growers avoid increased toxicity and unnecessary injury to pollinators.

In Silico Resources of Drug-Likeness as a Mirror: What Are We Lacking in Pesticide-Likeness?

Unfavorable bioavailability is an important aspect underlying the failure of drug candidates. Computational approaches for evaluating drug-likeness can minimize these risks. Over the past decades, computational approaches for evaluating drug-likeness have sped up the process of drug development and were also quickly derived to pesticide-likeness. As a result of many critical differences between drugs and pesticides, many kinds of methods for drug-likeness cannot be used for pesticide-likeness. Therefore, it is crucial to comprehensively compare and analyze the differences between drug-likeness and pesticide-likeness, which may provide a basis for solving the problems encountered during the evaluation of pesticide-likeness. Here, we systematically collected the recent advances of drug-likeness and pesticide-likeness and compared their characteristics. We also evaluated the current lack of studies on pesticide-likeness, the molecular descriptors and parameters adopted, the pesticide-likeness model on pesticide target organisms, and comprehensive analysis tools. This work may guide researchers to use appropriate methods for developing pesticide-likeness models. It may also aid non-specialists to understand some important concepts in drug-likeness and pesticide-likeness.

Whole-Body Acute Contact Toxicity of Formulated Insecticide Mixtures to Blue Orchard Bees (Osmia lignaria)

Blue orchard bees, [Osmia lignaria (Say) (Hymenoptera: Megachilidae)], have been developed as an important pollinator for orchard crops in North America over the last 40 years. The toxicity of several pesticides to O. lignaria and other Osmia species has been previously reported. However, the field-realistic toxicity of formulated premix insecticides comprised of multiple active ingredients (each with a different mode of action) to O. lignaria has not been assessed. Here, we use a customized spray tower in a laboratory setting to assess adult male and female whole-body direct contact exposure to four formulated pesticide mixtures: thiamethoxam + lambda-cyhalothrin (TLC), imidacloprid + beta-cyfluthrin (IBC), chlorantraniliprole + lambda-cyhalothrin (CLC) and methoxyfenozide + spinetoram (MS) by directly spraying anesthetized bees in Petri dishes. Separately, adult male and female whole-body direct contact exposure to formulated imidacloprid (I), beta-cyfluthrin (BC) and their 1:1 binary combination (IBC) was assessed using the same experimental method. Resulting mortality in each study was screened up to 96 h post-treatment to determine acute whole-body contact toxicity. In the first study, TLC and IBC resulted in statistically higher mortality at 24 and 48 h than the two other insecticide combinations tested. The CLC and MS combinations were slower acting and the highest mortality for O. lignaria exposed to these mixtures was recorded at 96 h. We did observe significant differences in toxicity between CLC and MS. In the second study, exposure to the 1:1 binary combination of IBC caused overall significantly higher mortality than exposure to I or BC alone. Both active ingredients alone, however, demonstrated equivalent levels of mortality to the 1:1 binary combination treatment at the 96 h observation reading, indicating increased speed of kill, but not necessarily increased toxicity. Significant differences in the onset of mortality following acute contact whole-body exposure to the formulated insecticide mixtures and individual active ingredients tested were consistently observed across all experiments in both studies.

EFSA's OpenFoodTox: An open source toxicological database on chemicals in food and feed and its future developments

Since its creation in 2002, the European Food Safety Authority (EFSA) has produced risk assessments for over 5000 substances in >2000 Scientific Opinions, Statements and Conclusions through the work of its Scientific Panels, Units and Scientific Committee. OpenFoodTox is an open source toxicological database, available both for download and data visualisation which provides data for all substances evaluated by EFSA including substance characterisation, links to EFSA's outputs, applicable legislations regulations, and a summary of hazard identification and hazard characterisation data for human health, animal health and ecological assessments. The database has been structured using OECD harmonised templates for reporting chemical test summaries (OHTs) to facilitate data sharing with stakeholders with an interest in chemical risk assessment, such as sister agencies, international scientific advisory bodies, and others. This manuscript provides a description of OpenFoodTox including data model, content and tools to download and search the database. Examples of applications of OpenFoodTox in chemical risk assessment are discussed including new quantitative structure-activity relationship (QSAR) models, integration into tools (OECD QSAR Toolbox and AMBIT-2.0), assessment of environmental footprints and testing of threshold of toxicological concern (TTC) values for food related compounds. Finally, future developments for OpenFoodTox 2.0 include the integration of new properties, such as physico-chemical properties, exposure data, toxicokinetic information; and the future integration within in silico modelling platforms such as QSAR models and physiologically-based kinetic models. Such structured in vivo, in vitro and in silico hazard data provide different lines of evidence which can be assembled, weighed and integrated using harmonised Weight of Evidence approaches to support the use of New Approach Methodologies (NAMs) in chemical risk assessment and the reduction of animal testing.

A Review of In Silico Tools as Alternatives to Animal Testing: Principles, Resources and Applications

Across the spectrum of industrial sectors, including pharmaceuticals, chemicals, personal care products, food additives and their associated regulatory agencies, there is a need to develop robust and reliable methods to reduce or replace animal testing. It is generally recognised that no single alternative method will be able to provide a one-to-one replacement for assays based on more complex toxicological endpoints. Hence, information from a combination of techniques is required. A greater understanding of the time and concentration-dependent mechanisms, underlying the interactions between chemicals and biological systems, and the sequence of events that can lead to apical effects, will help to move forward the science of reducing and replacing animal experiments. In silico modelling, in vitro assays, high-throughput screening, organ-on-a-chip technology, omics and mathematical biology, can provide complementary information to develop a complete picture of the potential response of an organism to a chemical stressor. Adverse outcome pathways (AOPs) and systems biology frameworks enable relevant information from diverse sources to be logically integrated. While individual researchers do not need to be experts across all disciplines, it is useful to have a fundamental understanding of what other areas of science have to offer, and how knowledge can be integrated with other disciplines. The purpose of this review is to provide those who are unfamiliar with predictive in silico tools, with a fundamental understanding of the underlying theory. Current applications, software, barriers to acceptance, new developments and the use of integrated approaches are all discussed, with additional resources being signposted for each of the topics.