Y. G.
Chushak
*a,
H. W.
Shows
b,
J. M.
Gearhart
a and
H. A.
Pangburn
c
aHenry M Jackson Foundation for the Advancement of Military Medicine, Wright-Patterson AFB, Ohio 45433, USA. E-mail: Chushak@yahoo.com
bBiological Sciences Department, Wright State University, Dayton, Ohio 45435, USA
cUnited States Air Force School of Aerospace Medicine, Aeromedical Research Department, Force Health Protection, Wright-Patterson AFB, Ohio 45433, USA
First published on 12th March 2018
There are many mechanisms of neurotoxicity that are initiated by the interaction of chemicals with different neurological targets. Under the U.S. Environmental Protection Agency's ToxCast program, the biological activity of thousands of chemicals was screened in biochemical and cell-based assays in a high-throughput manner. Two hundred sixteen assays in the ToxCast screening database were identified as targeting a total of 123 proteins having neurological functions according to the Gene Ontology database. Data from these assays were imported into the Organization for Economic Co-operation and Development QSAR Toolbox and used to predict neurological targets for chemical neurotoxins. Two sets of data were generated: one set was used to classify compounds as active or inactive and another set, composed of AC50s for only active compounds, was used to predict AC50 values for unknown chemicals. Chemical grouping and read-across within the QSAR Toolbox were used to identify neurologic targets and predict interactions for pyrethroids, a class of compounds known to elicit neurotoxic effects in humans. The classification prediction results showed 79% accuracy while AC50 predictions demonstrated mixed accuracy compared with the ToxCast screening data.
With the recent advances made in the field of in vitro high-throughput screening (HTS), it is now possible to screen the biological activity of large chemical libraries in a cost efficient and timely manner. In 2006, the U.S. Environmental Protection Agency (EPA) initiated the ToxCast program to develop and evaluate in vitro biochemical and cell-based assays for screening thousands of chemicals at multiple concentrations in the high-throughput mode.3 In Phase II of the program, approximately 1800 compounds were tested in ∼900 HTS assays. The chemicals in the library included pesticides, commercial compounds, and some failed pharmaceuticals. In 2008, the ToxCast program was merged with a large multiagency Tox21 collaboration. Under this new program, ∼8400 chemicals were screened in ∼70 HTS assays.4 These in vitro screenings generated an enormous volume of data, which are publicly available at https://www.epa.gov/chemical-research/toxicity-forecaster-toxcasttm-data. Although ToxCast and other in vitro screening programs provide a significant amount of information about the biological activities of thousands of chemicals, a great deal of information for millions of chemicals remains missing. Computational methods together with the HTS data offer a great opportunity to partially address this data gap and to identify molecular targets and other endpoints for chemical toxins of interest.
According to the European Chemical Agency (ECHA) guidance on information requirements and chemical safety assessment,5 two computational methods – (quantitative) structure–activity relationship [(Q)SAR] and grouping of chemicals with read-across – can be used for evaluating the intrinsic properties of chemicals. Both methods are based on the similarity principle, i.e., similar molecules have similar properties and biological activities are defined by the molecular structure.6 (Q)SAR methods are statistical in nature, as they try to correlate the molecular descriptors of chemicals with their properties. Furthermore, these methods are global in their scope, as they build models for all chemicals in the training dataset and make predictions for a wide range of chemicals within the applicability domain. QSAR modelling was applied to develop predictive models based on ToxCast HTS data. Some of the models were successful,7–9 while other QSAR models yielded low predictive performance.10,11
Grouping of chemicals into a category and read-across is another important technique for data gap filling in chemical hazard assessment. This approach is local in scope, as its predictions are based on the properties of a small set of similar chemicals. The OECD Guidance on Grouping of Chemicals defines a chemical category as a group of chemicals whose physico-chemical and toxicological properties are similar or follow a regular pattern as a result of structural similarity.12 The similarities may be based on common functional groups, common modes or mechanisms of action, common constituents or chemical classes, etc. Read-across is a technique to predict the unknown properties of chemicals of interest based on the known properties of chemicals in the same chemical group.5 Grouping of chemicals and the read-across technique are implemented in several freely available tools such as QSAR Toolbox,13 Toxmatch,14 and ToxRead.15 QSAR Toolbox is a software platform developed by the Organisation for Economic Co-operation and Development (OECD) in collaboration with the European Chemical Agency, intended to be used to fill data gaps in the hazard assessment of chemicals. QSAR Toolbox v.3.5 has a database with about 200000 chemicals provided by governmental and commercial institutions. Furthermore, it allows access to import custom databases and use data for hazard assessment. The main aim of the present study was to explore the application of the QSAR Toolbox and data from ToxCast HTS assays to identify and predict the molecular interactions of chemical neurotoxins with their targets.
Recently, activities of 86 compounds from the ToxCast library were tested in neuronal cultures on multi-well microelectrode arrays (MEAs).16 Activities of these compounds on MEAs were compared with their activities on 20 ToxCast binding assays that measured the interaction of chemicals with 8 different ion channels. In our approach, we identified 123 proteins from ToxCast HTS assays that are related to neurological functions. This set of proteins includes ion channels, G-protein coupled receptors, nuclear receptors, transporters, and enzymes as potential neurological targets. The developed approach was evaluated by predicting neurological targets for pyrethroids and comparing the predicted results with ToxCast screening data.
The DB database combines detailed drug data with comprehensive drug–target information18 containing 8261 drugs and 4338 non-redundant proteins that are linked to these drug entries. Twenty-nine chemicals from the DB database were screened in selected ToxCast assays and were used to evaluate the activity of neurological proteins in ToxCast screening.
The PDSP Ki database, which is funded by the U.S. National Institute of Mental Health Psychoactive Drug Screening Program, serves as a data warehouse for published and internally derived Ki, or affinity, values for a large number of drugs and drug candidates at an expanding number of G-protein coupled receptors, ion channels, transporters, and enzymes.19 Currently, it has ∼60000 Ki values. Seventeen chemicals from this database were tested on their targets in selected ToxCast HTS assays and were used to evaluate the performance of ToxCast assays.
Another database, the Toxin and Toxin Target Database (http://www.t3db.ca/), also provides the mechanisms of toxicity and target proteins for toxins. However, ToxCast HTS data are already included in this database for chemical–protein associations. Therefore, this database information was not used for evaluation of ToxCast screening assays to avoid bias.
Fig. 1 Distribution of the number of active chemicals for neurological proteins. The inset shows 22 proteins with more than 100 active chemicals. |
As previously mentioned, neurological proteins were screened in 263 ToxCast high-throughput assays. In the ToxCast screening protocol, initially all chemicals were screened in each assay at a single concentration of 25 μM (10 μM for CYP assays) to identify active chemical-assay combinations in which the mean assay signal differed by at least 30% from the DMSO control signal.28 Next, for active chemical-assay pairs, each chemical was run in eight point serial dilutions starting from a top concentration of 50 μM down to 0.023 μM to estimate the AC50s. We generated two sets of data for neurological proteins from the ToxCast dataset: one set for classification of compounds as active or inactive and another set for prediction of AC50 values which were subsequently imported into the QSAR Toolbox to identify and predict the molecular interactions of chemical neurotoxins with their targets.
We evaluated the performance of the ToxCast HTS assays on neurological proteins by using chemical neurotoxins with known mechanisms of action from the two databases: DrugBank and Ki databases from the Psychoactive Drug Screening Program. Twenty-nine chemicals from the DB database were screened in ToxCast assays on neurological protein targets. The results of the comparison of drug–target pairs from DB and ToxCast screening are shown in Table 1. Only proteins that were screened in ToxCast assays are presented in this table. Chemicals in green rows were found to have all of their targets from DB represented in the ToxCast database. For chemicals in the yellow rows, only some DB chemical–target pairs were reproduced in ToxCast screening, and for chemicals in the red rows, their targets were not present in the ToxCast screening. Overall, there are 46 chemical–protein interactions between the selected compounds and neurological proteins in the DB database. Thirty-three of these interactions, corresponding to ∼72% of the interactions, were reproduced in the ToxCast HTS screening. ESR1 is a target for seven of the screened drugs in the DB database, and it was identified as a target for all selected chemicals in the ToxCast screening. The androgen receptor (AR) interacts with six drugs in DB, and the progesterone receptor (PGR) interacts with five drugs; all these interactions were also identified in ToxCast screening. Among the 13 missing chemical–protein pairs, 4 of them involve the interactions of prostaglandin-endoperoxide synthase PTGS2 with common analgesics such as aspirin or acetaminophen. Only one interaction of celecoxib with PTGS2 was identified in ToxCast screening.
Seventeen chemicals from the Ki database were screened in the ToxCast assays, and their neurological targets are shown in Table 2. The colouring schema in Table 2 is the same as in Table 1. Although only 17 chemicals from the Ki database were screened in the ToxCast assays, there are 74 chemical–protein interactions. The majority of these chemical–protein pairs represent the interactions of two drugs – haloperidol and volinanserin – with their targets. Haloperidol, which is used to treat schizophrenia and other psychoses, has 5 protein targets in the DB database (only 3 of them were used in the ToxCast screening), but it has 27 protein targets in the PDSP database. In the ToxCast screening, haloperidol was active on 37 neurological protein targets, which indicates the very high activity of this drug having potential side effects. Volinanserin was tested in clinical trials as a potential antipsychotic drug, but it was never advanced to the market place. It interacts with 11 protein targets in the PDSP database and has 26 interactions in the ToxCast screening assays. Among the chemicals whose protein interactions were not identified in the ToxCast screening, celecoxib and reserpine are also present in the DB database. The anti-inflammatory drug celecoxib interacts with prostaglandin synthase PTGS2 in DB, and this interaction was captured in ToxCast screening. However, its interaction with the sodium-dependent dopamine transporter SLC6A3 from the Ki database was not observed. The antipsychotic drug reserpine has the synaptic vesicular amine transporter SLC18A2 listed as a target in DB, which was identified in ToxCast, while in PDSP it interacts with serotonin and dopamine receptors that were not observed in the ToxCast screening. Overall, out of 73 chemical–protein interactions between the selected compounds and neurological proteins in the Ki database, 58 chemical–protein pairs that correspond to 80% of interactions were identified in the ToxCast screening assays.
The neuronal activities of 86 compounds from the ToxCast library were recently examined in primary cortical cultures with MEAs.16 The changes in the weighted mean firing rate were used as indicators of chemical activities on neuronal networks. The effects of tested chemicals on spontaneous neuronal activity in MEAs were compared with their activities on 20 ToxCast assays that measured the interaction of chemicals with ion channels. Pyrethroids were among the classes of chemicals that affected neuronal activity in MEAs. The primary mechanism of pyrethroid neurotoxicity is proposed to be via their effects on the sodium channels of nerve cells.29 However, only one of the six studied pyrethroid compounds – allethrin – showed activity in the ToxCast ion channel assays, while all but permethrin were active in MEA measurements.16 This indicates the presence of different mechanisms of toxicity for pyrethroids. We selected chemicals from this group to evaluate the predictive capability of the proposed approach.
The first step in the read-across approach is to identify a group of chemicals from the ToxCast screening database that belong to the same class as a studied compound. There are several ways to perform a similarity search in the QSAR Toolbox (“Category Definition”). We used two terms, “pyrethroids” and “esters”, to identify pyrethroids in the ToxCast database. The search recognized 12 compounds that are shown in Fig. 2. Two different types of pyrethroids are recognized based on the differences in their structure and the symptoms of poisoning.29 Structures of four of the selected compounds include a cyano group, and these chemicals belong to the Type II pyrethroid group, while the remaining chemicals, which lack a cyano group, belong to the Type I pyrethroid group.
The majority of pyrethroids are active in 10 to 20 ToxCast HTS assays with neurological proteins (Fig. 3). However, allethrin shows very high activity as it is active in 26 assays, while bifenthrin is only active in 4 assays and tefluthrin in 6 assays. This is a clear indication that some compounds in the pyrethroid group have unique structural features that are reflected in their distinctive activities. As for neuronal targets, the ATG_DR5_CIS assay that targeted retinoic acid receptors RARA, RARB, and RARG exhibited activity for 10 of the 12 pyrethroids examined. Another assay that also targeted RARA – ATG_RARa_TRANS – has seven active compounds, indicating that retinoic acid receptors are targeted by pyrethroids.
Fig. 3 Number of active assays with neurological proteins for pyrethroids from the database of chemicals. |
Another potential target of pyrethroid compounds is the dopamine transporter SLC63. Two assays that targeted human and guinea pig transporters – NVS_TR_hDAT and NVS_NVS_gDAT – have nine and five active pyrethroids, respectively. Among the neuronal proteins with only one or two active chemicals are GABA receptor GABRA1, which is affected only by allethrin, and muscarinic cholinergic receptors CHRM1 and CHMR3, which are affected by prallethrin. Since we can only make predictions for protein targets that are affected by a group of chemicals and not just by a single chemical, we selected a set of 14 ToxCast assays with 5 or more active pyrethroids to assess the predictive capability of the proposed approach. The detailed information about these assays is presented in ESI Table S2.† Eight of the studied assays are cell-based assays and the effect of chemicals was investigated by measuring the activity of proteins or changes in the protein expression level. Five of the biochemical assays measured protein–ligand binding, while in one assay enzymatic activity was measured to evaluate the activity of chemicals.
The predictions for pyrethroids were based on the activity of five nearest neighbour chemicals, and data from target chemicals were not used in the prediction. The detailed classification prediction results are presented in ESI Table S3.† For four compounds – allethrin, S-bioallethrin, prallethrin, and fenpropathrin – activity in only one of the studied assays was incorrectly predicted. On the other hand, activities for tetramethrin and for bifenthrin in 7 and 5 assays, respectively, were incorrectly predicted. Tetramethrin is the only compound among the studied pyrethroids that has an imide group that potentially explains the high number of incorrect predictions. Bifenthrin, as mentioned before, showed extremely low activity in the ToxCast screening, as it was active with neurological proteins only in four assays (see Fig. 3). The summary of classification predictions for pyrethroid compounds is presented in Table 3. Overall, there are 168 chemical–protein pairs; 89 are active and 79 are inactive. Eighty three percent of active interactions, 73 of 89 interactions in the ToxCast database, were correctly predicted using our approach. We also correctly predicted 60 inactive chemical–protein pairs, which corresponds to a 76% hit rate.
ToxCast | Predicted | ||
---|---|---|---|
Active | Inactive | Hit rate | |
Active | 73 | 16 | 83% |
Inactive | 19 | 60 | 76% |
Accuracy | 79% | 79% |
Table 4 shows the predicted AC50 values (in mg L−1) for a set of five pyrethroid compounds in comparison with the ToxCast screening results. The AC50 prediction results for all 12 pyrethroids are presented in ESI Table S4.† For cyfluthrin, we correctly predicted activity in 11 out of 14 assays. The predicted AC50 values for cyfluthrin are in good correlation with those observed experimentally. Activity predictions for permethrin yielded incorrect results for two assays while only a few predicted AC50 values were in agreement with those measured in ToxCast assays. In several in vitro assays, permethrin shows very low activity, with AC50 values in the range of 8.6–16.0 mg L−1. This is in agreement with the low potency of permethrin, which showed no effect on the neuronal network activity measured using the MEA technique.16 The predicted values are based on the values measured for a group of similar chemicals. As a result, the predicted AC50 values for permethrin in these assays are lower compared with the observed. For cypermethrin, only the BSK_3C_HLADR assay activity was incorrectly predicted, while the predicted AC50 values show, in general, a good correlation with the measured ones. Although tetramethrin is active in 16 neurological assays, it is active only in 5 of the assays selected for the evaluation. As mentioned before, classification incorrectly predicted the activity of tetramethrin in seven of the ToxCast assays. The predicted AC50 values also showed mixed accuracy compared with the ToxCast data. These results indicate that it is currently difficult to predict the activity of a compound that has different structural and activity properties compared with the other chemicals in the same group.
Grouping of chemicals. The first and most critical part in the read-across prediction is finding similar chemicals in the database. There is no unique way for grouping chemicals into a category. Therefore, the same chemical can belong to different groups based on the similarity criteria.30 We used two terms to identify pyrethroids in the ToxCast database, “pyrethroids” and “esters”. However, 12 pyrethroid compounds that were identified showed very different activities in the ToxCast assays. For example, allethrin is active in 26 assays, while bifenthrin is active only in 4 assays and tefluthrin is active in 6 assays. Furthermore, this set of pyrethroids contains eight Type I and four Type II pyrethroids classified according to eliciting different symptomology. Clearly, the difference in the chemical structure of these chemicals resulted in their different activities. One way to improve the similarity is to subcategorize the list of chemicals based on the structural similarity on a chemical-by-chemical basis.
Low activity of chemicals. In some cases, a chemical showed low activity in a ToxCast assay with multiple points above the baseline, but the AC50 value could not be calculated from the dose–response data.28 In such cases, chemicals were marked as “inactive” in the ToxCast database, making it difficult to identify the active/inactive boundary for chemicals with low activity.
The number of active/inactive compounds in the classification dataset should be similar. The ATG_DR5_CIS assay has 10 active pyrethroids out of 12 and several assays have 9 active pyrethroids. It is impossible to predict inactive compounds for ATG_DR5_CIS (like tetramethrin and bifenthrin) when 10 neighbours are active in the assay. The same is true also for assays that have just 1–3 active compounds. Therefore, these assays were not considered in the current approach. In this case, a different grouping schema needs to be applied that will identify a different set of similar chemicals.
Descriptors and applicability domain. The QSAR Toolbox allows the use of a different set of molecular descriptors for read-across prediction. Special attention is needed when selecting descriptors to be sure that they show correlation with the AC50 values and the tested chemical is within the applicability domain for the selected descriptors. Another approach is to use trend analysis, which can improve the predicted results when the chemical is outside the applicability domain for the selected descriptor. However, that approach was not used in the current study.
Our results demonstrate that ToxCast in vitro screening data and the QSAR Toolbox can be used to identify neurological targets and predict AC50 values for chemical neurotoxins that have similar chemicals in the ToxCast database. This approach can be combined with MEA measurements to link molecular interactions to cellular response. Very recently, MEA techniques were used to evaluate the effect of a set of pyrethroids on the spontaneous activity of cortical neuronal networks.31 Besides the weighted mean firing rate, authors also measured other burst parameters such as the mean burst duration, the mean interspike interval in burst, etc. In the present study, we show that pyrethroids affect multiple protein targets that can result in complex neurological responses. Therefore, it is of interest to investigate the correlation of molecular interactions with the multiple endpoints in MEA measurements. Such a combined in vitro/in silico approach would allow us to connect the MIE with the cellular response of chemical neurotoxins.
Footnote |
† Electronic supplementary information (ESI) available: Tables S1–S4. See DOI: 10.1039/c7tx00268h |
This journal is © The Royal Society of Chemistry 2018 |