The VERA software: Implementation of the acute fish toxicity endpoint and its application to pharmaceutical compounds.

The Virtual Extensive Read-Across software (VERA) is a new tool for read-across using a global similarity score, molecular groups, and structural alerts to find clusters of similar substances; these clusters are then used to identify suitable similar substances and make an assessment for the target substance. A beta version of VERA GUI is free and available at vegahub.eu; the source code of the VERA algorithm is available on GitHub. In the past we described its use to assess carcinogenicity, a classification endpoint. The aim here is to extend the automated read-across approach to assess continuous endpoints as well. We addressed acute fish toxicity. VERA evaluation on the acute fish toxicity endpoint was done on a dataset containing general substances (pesticides, industrial products, biocides, etc.), obtaining an overall R2 of 0.68. We employed the VERA algorithm also on active pharmaceutical ingredients (APIs). We included a portion of the APIs in the training dataset to predict APIs, successfully achieving an overall R2 of 0.63. VERA evaluates the assessment's reliability, and we reached an R2 of 0.78 and Root Mean Square Error (RMSE) of 0.44 for predictions with high reliability.

New approach methodologies to facilitate and improve the hazard assessment of non-genotoxic carcinogens-a PARC project.

Carcinogenic chemicals, or their metabolites, can be classified as genotoxic or non-genotoxic carcinogens (NGTxCs). Genotoxic compounds induce DNA damage, which can be detected by an established in vitro and in vivo battery of genotoxicity assays. For NGTxCs, DNA is not the primary target, and the possible modes of action (MoA) of NGTxCs are much more diverse than those of genotoxic compounds, and there is no specific in vitro assay for detecting NGTxCs. Therefore, the evaluation of the carcinogenic potential is still dependent on long-term studies in rodents. This 2-year bioassay, mainly applied for testing agrochemicals and pharmaceuticals, is time-consuming, costly and requires very high numbers of animals. More importantly, its relevance for human risk assessment is questionable due to the limited predictivity for human cancer risk, especially with regard to NGTxCs. Thus, there is an urgent need for a transition to new approach methodologies (NAMs), integrating human-relevant in vitro assays and in silico tools that better exploit the current knowledge of the multiple processes involved in carcinogenesis into a modern safety assessment toolbox. Here, we describe an integrative project that aims to use a variety of novel approaches to detect the carcinogenic potential of NGTxCs based on different mechanisms and pathways involved in carcinogenesis. The aim of this project is to contribute suitable assays for the safety assessment toolbox for an efficient and improved, internationally recognized hazard assessment of NGTxCs, and ultimately to contribute to reliable mechanism-based next-generation risk assessment for chemical carcinogens.

Virtual Extensive Read-Across: A New Open-Access Software for Chemical Read-Across and Its Application to the Carcinogenicity Assessment of Botanicals.

Read-across applies the principle of similarity to identify the most similar substances to represent a given target substance in data-poor situations. However, differences between the target and the source substances exist. The present study aims to screen and assess the effect of the key components in a molecule which may escape the evaluation for read-across based only on the most similar substance(s) using a new open-access software: Virtual Extensive Read-Across (VERA). VERA provides a means to assess similarity between chemicals using structural alerts specific to the property, pre-defined molecular groups and structural similarity. The software finds the most similar compounds with a certain feature, e.g., structural alerts and molecular groups, and provides clusters of similar substances while comparing these similar substances within different clusters. Carcinogenicity is a complex endpoint with several mechanisms, requiring resource intensive experimental bioassays and a large number of animals; as such, the use of read-across as part of new approach methodologies would support carcinogenicity assessment. To test the VERA software, carcinogenicity was selected as the endpoint of interest for a range of botanicals. VERA correctly labelled 70% of the botanicals, indicating the most similar substances and the main features associated with carcinogenicity.

In Silico Methods for Carcinogenicity Assessment.

Screening compounds for potential carcinogenicity is of major importance for prevention of environmentally induced cancers. A large sequence of predictive models, ranging from short-term biological assays (e.g., mutagenicity tests) to theoretical models, has been attempted in this field. Theoretical approaches such as (Q)SAR are highly desirable for identifying carcinogens, since they actively promote the replacement, reduction, and refinement of animal tests. This chapter reports and describes some of the most noted (Q)SAR models based on human expert knowledge and statistical approaches, aiming at predicting the carcinogenicity of chemicals. Additionally, the performance of the selected models has been evaluated, and the results are interpreted in details by applying these predictive models to some pharmaceutical molecules.

Exploration of structural requirements for azole chemicals towards human aromatase CYP19A1 activity: Classification modeling, structure-activity relationships and read-across study.

Human aromatase, also called CYP19A1, plays a major role in the conversion of androgens into estrogens. Inhibition of aromatase is an important target for estrogen receptor (ER)-responsive breast cancer therapy. Use of azole compounds as aromatase inhibitors is widespread despite their low selectivity. A toxicological evaluation of commonly used azole-based drugs and agrochemicals with respect to CYP19A1 is currently requested by the European Union- Registration, Evaluation, Authorization and Restriction of Chemicals (EU-REACH) regulations due to their potential as endocrine disruptors. In this connection, identification of structural alerts (SAs) is an effective strategy for the toxicological assessment and safe drug design. The present study describes the identification of SAs of azole-based chemicals as guiding experts to predict the aromatase activity. Total 21 SAs associated with aromatase activity were extracted from dataset of 326 azole-based drugs/chemicals obtained from Tox21 library. A cross-validated classification model having high accuracy (error rate 5%) was proposed which can precisely classify azole chemicals into active/inactive toward aromatase. In addition, mechanistic details and toxicological properties (agonism/antagonism) of azoles with respect to aromatase were explored by comparing active and inactive chemicals using structure-activity relationships (SAR). Lastly, few structural alerts were applied to form chemical categories for read-across applications.

Structures of Endocrine-Disrupting Chemicals Correlate with the Activation of 12 Classic Nuclear Receptors.

Endocrine-disrupting chemicals (EDCs) can inadvertently interact with 12 classic nuclear receptors (NRs) that disrupt the endocrine system and cause adverse effects. There is no widely accepted understanding about what structural features make thousands of EDCs able to activate different NRs as well as how these structural features exert their functions and induce different outcomes at the cellular level. This paper applies the hierarchical characteristic fragment methodology and high-throughput screening molecular docking to comprehensively explore the structural and functional features of EDCs for the 12 NRs based on more than 7000 chemicals from curated datasets. EDCs share three levels of key fragments. The primary and secondary fragments are associated with the binding of EDCs to four groups of receptors: steroidal nuclear receptors (SNRs, including androgen, estrogen, glucocorticoid, mineralocorticoid, and progesterone), retinoic acid receptors, thyroid hormone receptors, and vitamin D receptors. The tertiary fragments determine the activity type by interacting with two key locations in the ligand-binding domains of NRs (N-H5-H3-C and N-H7-H11-C for SNRs and N-H5-H5'-H2'-H3-C and N-H6'-H11-C for non-SNRs). The resulting compiled structural fragments of EDCs together with elucidated compound NR binding modes provide a framework for understanding the interactions between EDCs and NRs, facilitating faster and more accurate screening of EDCs for multiple NRs in the future.

Quantitative Structure-Activity Relationship Modeling of the Amplex Ultrared Assay to Predict Thyroperoxidase Inhibitory Activity.

The thyroid system plays a major role in the regulation of several physiological processes. The dysregulation of the thyroid system caused by the interference of xenobiotics and contaminants may bring to pathologies like hyper- and hypothyroidism and it has been recently correlated with adverse outcomes leading to cancer, obesity, diabetes and neurodevelopmental disorders. Thyroid disruption can occur at several levels. For example, the inhibition of thyroperoxidase (TPO) enzyme, which catalyses the synthesis of thyroid hormones, may cause dysfunctions related to hypothyroidism. The inhibition of the TPO enzyme can occur as a consequence of prolonged exposure to chemical compounds, for this reason it is of utmost importance to identify alternative methods to evaluate the large amount of pollutants and other chemicals that may pose a potential hazard to the human health. In this work, quantitative structure-activity relationship (QSAR) models to predict the TPO inhibitory potential of chemicals are presented. Models are developed by means of several machine learning and data selection approaches, and are based on data obtained in vitro with the Amplex UltraRed-thyroperoxidase (AUR-TPO) assay. Balancing methods and feature selection are applied during model development. Models are rigorously evaluated through internal and external validation. Based on validation results, two models based on Balanced Random Forest (BRF) and K-Nearest Neighbours (KNN) algorithms were selected for a further validation phase, that leads predictive performance (BA = 0.76-0.78 on external data) that is comparable with the reported experimental variability of the AUR-TPO assay (BA ∼0.70). Finally, a consensus between the two models was proposed (BA = 0.82). Based on the predictive performance, these models can be considered suitable for toxicity screening of environmental chemicals.

A descriptor-based analysis to highlight the mechanistic rationale of mutagenicity.

Cancer is a main concern for human health and there is a need of alternative methodologies to rapidly screen large quantitative of compounds that may represent a toxicological risk. Here a statistical analyses is performed on a benchmark database of experimental Ames data to identify chemical descriptors discriminating mutagens and non-mutagens. A total of 53 activating and deactivating modulators are identified, that flagged respectively a percentage of mutagen and non-mutagen up to 87%. Modulators are further combined to form synergistic cross-terms, accounting for the effect that combined properties may have on the final toxicity. Exclusion rules are defined as exception to the modulators. Synergistic cross-terms and exclusion rules improve the enrichment of mutagens/non-mutagens with respect of the original abundance in the dataset to values higher than 95%. The external predictivity of modulators and cross-terms reach balanced accuracy up to 0.775 that is analogous to other mutagenicity models from the literature, confirming the suitability of the rules to real-life screening of chemicals. Modulators are discussed for their mechanistic link to mutagenicity. This analysis confirms the key role of some properties (polarizability, shape, mass, presence of reactive functional groups or unsaturated planar systems) as driving elements for the initiation of the mutagenicity.

The index of ideality of correlation improves the predictive potential of models of the antioxidant activity of tripeptides from frog skin (Litoria rubella).

It is usually held that good-quality models for the biological activity of peptides must take into account their 3D architecture and descriptors of quantum mechanics. However, the present study shows that it is possible to build up models without these complex calculations. The structure of tripeptides represented by sequences of one-symbol abbreviations of the corresponding amino acids serves to build up quantitative structure-activity relationships for the antioxidant activity of tripeptides from frog skin. The statistical quality of the best model for the validation set is n = 27, r2 = 0.93, RMSE = 0.15.

Pesticides, cosmetics, drugs: identical and opposite influences of various molecular features as measures of endpoints similarity and dissimilarity.

The similarity is an important category in natural sciences. A measure of similarity for a group of various biochemical endpoints is suggested. The list of examined endpoints contains (1) toxicity of pesticides towards rainbow trout; (2) human skin sensitization; (3) mutagenicity; (4) toxicity of psychotropic drugs; and (5) anti HIV activity. Further applying and evolution of the suggested approach is discussed. In particular, the conception of the similarity (dissimilarity) of endpoints can play the role of a "useful bridge" between quantitative structure property/activity relationships (QSPRs/QSARs) and read-across technique.

Review and priority setting for substances that are listed without a specific migration limit in Table 1 of Annex 1 of Regulation 10/2011 on plastic materials and articles intended to come into contact with food.

The EFSA Panel on Food Contact Materials, Enzymes and Processing Aids (CEP) was requested by the European Commission to review the substances for which a Specific Migration Limit (SML) is not assigned in Regulation (EU) No 10/2011. These substances had been covered by the Generic SML of 60 mg/kg food, but with Regulation (EU) 2016/1416 it was removed, necessitating their re-examination. EFSA was requested to identify those substances requiring an SML to ensure the authorisation is sufficiently protective to health, grouping them in high, medium and low priority to serve as the basis for future re-evaluations of individual substances. The CEP Panel established a stepwise procedure. This took into account existing hazard assessments for each substance on carcinogenicity/mutagenicity/reprotoxicity (CMR), bioaccumulation and endocrine disruptor (ED) properties along with the use of in silico generated predictions on genotoxicity. Molecular weights and boiling points were considered with regard to their effect on potential consumer exposure. This prioritisation procedure was applied to a total of 451 substances, from which 78 substances were eliminated at the outset, as they had previously been evaluated by EFSA as food contact substances. For 89 substances, the Panel concluded that a migration limit should not be needed. These are in the lists 0 and 1 of the Scientific Committee for Food (SCF), defined as substances for which an Acceptable Daily Intake (ADI) does not need to be established, along with substances that are controlled by existing restrictions and/or generic limits. Of the remaining 284 substances, 179 were placed into the low priority group, 102 were placed into the medium priority group and 3 were placed into the high priority group, i.e. salicylic acid (FCM No 121), styrene (FCM No 193) and lauric acid, vinyl ester (FCM No 436).

QSAR model for pesticides toxicity to Rainbow Trout based on "ideal correlations".

Pesticides have an impact on the aquatic environment, with ecological effects. The regulation of this impact is of key importance. One of the components of the planning of agricultural and industrial activities is the development of databases and models in order to identify substances that may cause damage. In this study, a quantitative structure-activity relationship (QSAR) approach was established for the prediction of acute toxicity toward rainbow trout of various pesticides. The so-called index of ideality of correlation is the main component of this approach. The validation of this approach has been carried out with three random splits into the training and validation sets. The range of statistical quality of models obtained here for the validation set is R2 = [0.81-0.86] and RMSE = [0.55-0.65].

Structures of Endocrine-Disrupting Chemicals Determine Binding to and Activation of the Estrogen Receptor α and Androgen Receptor.

Endocrine-disrupting chemicals (EDCs) can interact with nuclear receptors, including estrogen receptor α (ERα) and androgen receptor (AR), to affect the normal endocrine system function, causing severe symptoms. Limited studies queried the EDC mechanisms, focusing on limited chemicals or a set of structurally similar compounds. It remained uncertain how hundreds of diverse EDCs could bind to ERα and AR and cause distinct functional consequences. Here, we employed a series of computational methodologies to investigate the structural features of EDCs that bind to and activate ERα and AR based on more than 4000 compounds. We used molecular docking and molecular dynamics simulations to elucidate the functional consequences and validated structure-function correlations experimentally using a time-resolved fluorescence resonance energy-transfer assay. We found that EDCs share three levels of key fragments. Primary (20 for ERα and 18 for AR) and secondary fragments (38 for ERα and 29 for AR) are responsible for the binding to receptors, and tertiary fragments determine the activity type (agonist, antagonist, or mixed). In summary, our study provides a general mechanism for the EDC function. Discovering the three levels of key fragments may drive fast screening and evaluation of potential EDCs from large sets of commercially used synthetic compounds.

Ecosystem ecology: Models for acute toxicity of pesticides towards Daphnia magna.

Quantitative structure - activity relationships (QSARs) which are obtained with a representation of the molecular architecture via simplified molecular input-line entry system (SMILES) are applied to build up predictive models of acute toxicity of pesticides towards Daphnia magna. The acute toxicity towards Daphnia magna is an adequate measure of the ecological impact of various substances. The Monte Carlo technique is the basis to build up the above QSAR models. The statistical quality of suggested models is good: the best model is characterized by n = 103, R2 = 0.76, RMSE = 0.91 (training set); n = 53, R2 = 0.82, RMSE = 0.87 (validation set). The approach provides the mechanistic interpretation (e.g. aromaticity and branching of carbon skeleton are promoters of increase for toxicity towards Daphnia magna in the case of the examined set of pesticides). The approach is attractive to build up predictive models since instead of a large number of different molecular descriptors the corresponding model is based on solely one optimal descriptor calculated with SMILES and all necessary calculations can be done using the CORAL software available on the Internet (http://ww.insilico.eu/coral).

Homology Modeling of the Human P-glycoprotein (ABCB1) and Insights into Ligand Binding through Molecular Docking Studies.

The ABCB1 transporter also known as P-glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette super-family of transporters; it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping the compounds out of cells. P-gp contributes to a decrease of toxicity and possesses broad substrate specificity. It is involved in the failure of numerous anticancer and antiviral chemotherapies due to the multidrug resistance (MDR) phenomenon, where it removes the chemotherapeutics out of the targeted cells. Understanding the details of the ligand-P-gp interaction is therefore crucial for the development of drugs that might overcome the MRD phenomenon and for obtaining a more effective prediction of the toxicity of certain compounds. In this work, an in silico modeling was performed using homology modeling and molecular docking methods with the aim of better understanding the ligand-P-gp interactions. Based on different mouse P-gp structural templates from the PDB repository, a 3D model of the human P-gp (hP-gp) was constructed by means of protein homology modeling. The homology model was then used to perform molecular docking calculations on a set of thirteen compounds, including some well-known compounds that interact with P-gp as substrates, inhibitors, or both. The sum of ranking differences (SRD) was employed for the comparison of the different scoring functions used in the docking calculations. A consensus-ranking scheme was employed for the selection of the top-ranked pose for each docked ligand. The docking results showed that a high number of π interactions, mainly π-sigma, π-alkyl, and π-π type of interactions, together with the simultaneous presence of hydrogen bond interactions contribute to the stability of the ligand-protein complex in the binding site. It was also observed that some interacting residues in hP-gp are the same when compared to those observed in a co-crystallized ligand (PBDE-100) with mouse P-gp (PDB ID: 4XWK). Our in silico approach is consistent with available experimental results regarding P-gp efflux transport assay; therefore it could be useful in the prediction of the role of new compounds in systemic toxicity.

Value and limitation of structure-based profilers to characterize developmental and reproductive toxicity potential.

The uncertainty regarding the safety of chemicals leaching from food packaging triggers attention. In silico models provide solutions for screening of these chemicals, since many are toxicologically uncharacterized. For hazard assessment, information on developmental and reproductive toxicity (DART) is needed. The possibility to apply in silico toxicology to identify and quantify DART alerts was investigated. Open-source models and profilers were applied to 195 packaging chemicals and analogues. An approach based on DART and estrogen receptor (ER) binding profilers and molecular docking was able to identify all except for one chemical with documented DART properties. Twenty percent of the chemicals in the database known to be negative in experimental studies were classified as positive. The scheme was then applied to 121 untested chemicals. Alerts were identified for sixteen of them, five being packaging substances, the others structural analogues. Read-across was then developed to translate alerts into quantitative toxicological values. They can be used to calculate margins of exposure (MoE), the size of which reflects safety concern. The application of this approach appears valuable for hazard characterization of toxicologically untested packaging migrants. It is an alternative to the use of default uncertainty factor (UF) applied to animal chronic toxicity value to handle absence of DART data in hazard characterization.

CoMPARA: Collaborative Modeling Project for Androgen Receptor Activity.

BACKGROUND: Endocrine disrupting chemicals (EDCs) are xenobiotics that mimic the interaction of natural hormones and alter synthesis, transport, or metabolic pathways. The prospect of EDCs causing adverse health effects in humans and wildlife has led to the development of scientific and regulatory approaches for evaluating bioactivity. This need is being addressed using high-throughput screening (HTS) in vitro approaches and computational modeling. OBJECTIVES: In support of the Endocrine Disruptor Screening Program, the U.S. Environmental Protection Agency (EPA) led two worldwide consortiums to virtually screen chemicals for their potential estrogenic and androgenic activities. Here, we describe the Collaborative Modeling Project for Androgen Receptor Activity (CoMPARA) efforts, which follows the steps of the Collaborative Estrogen Receptor Activity Prediction Project (CERAPP). METHODS: The CoMPARA list of screened chemicals built on CERAPP's list of 32,464 chemicals to include additional chemicals of interest, as well as simulated ToxCast™ metabolites, totaling 55,450 chemical structures. Computational toxicology scientists from 25 international groups contributed 91 predictive models for binding, agonist, and antagonist activity predictions. Models were underpinned by a common training set of 1,746 chemicals compiled from a combined data set of 11 ToxCast™/Tox21 HTS in vitro assays. RESULTS: The resulting models were evaluated using curated literature data extracted from different sources. To overcome the limitations of single-model approaches, CoMPARA predictions were combined into consensus models that provided averaged predictive accuracy of approximately 80% for the evaluation set. DISCUSSION: The strengths and limitations of the consensus predictions were discussed with example chemicals; then, the models were implemented into the free and open-source OPERA application to enable screening of new chemicals with a defined applicability domain and accuracy assessment. This implementation was used to screen the entire EPA DSSTox database of ∼875,000 chemicals, and their predicted AR activities have been made available on the EPA CompTox Chemicals dashboard and National Toxicology Program's Integrated Chemical Environment. https://doi.org/10.1289/EHP5580.

The using of the Index of Ideality of Correlation (IIC) to improve predictive potential of models of water solubility for pesticides.

Models for water solubility of pesticides suggested in this manuscript are important data from point of view of ecologic engineering. The Index of Ideality of Correlation (IIC) of groups of quantitative structure-property relationships (QSPRs) for water solubility of pesticides related to the calibration sets was used to identify good in silico models. This comparison confirmed the high IIC set provides better statistical quality of the model for the validation set. Though there are large databases on solubility, the reliable prediction of the endpoint for new substances which are potential pesticides is an important ecologic task. Unfortunately, predictive models for various endpoints suffer overtraining, and the IIC serves to avoid or at least reduce this. Thus, the approach suggested has both theoretical and economic effects for ecology.

QSAR Models for Human Carcinogenicity: An Assessment Based on Oral and Inhalation Slope Factors.

Carcinogenicity is a crucial endpoint for the safety assessment of chemicals and products. During the last few decades, the development of quantitative structure-activity relationship ((Q)SAR) models has gained importance for regulatory use, in combination with in vitro testing or expert-based reasoning. Several classification models can now predict both human and rat carcinogenicity, but there are few models to quantitatively assess carcinogenicity in humans. To our knowledge, slope factor (SF), a parameter describing carcinogenicity potential used especially for human risk assessment of contaminated sites, has never been modeled for both inhalation and oral exposures. In this study, we developed classification and regression models for inhalation and oral SFs using data from the Risk Assessment Information System (RAIS) and different machine learning approaches. The models performed well in classification, with accuracies for the external set of 0.76 and 0.74 for oral and inhalation exposure, respectively, and r2 values of 0.57 and 0.65 in the regression models for oral and inhalation SFs in external validation. These models might therefore support regulators in (de)prioritizing substances for regulatory action and in weighing evidence in the context of chemical safety assessments. Moreover, these models are implemented on the VEGA platform and are now freely downloadable online.