Subscriber access provided by UNIV OF DURHAM
Ecotoxicology and Human Environmental Health
High-Throughput Screening and Quantitative Chemical Ranking for Sodium Iodide Symporter (NIS) Inhibitors in ToxCast Phase I Chemical Library Jun Wang, Daniel Hallinger, Ashley Murr, Angela Buckalew, Steven O'Neal Simmons, Susan C. Laws, and Tammy Stoker Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b06145 • Publication Date (Web): 03 Apr 2018 Downloaded from http://pubs.acs.org on April 3, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 34
Environmental Science & Technology
1
High-Throughput Screening and Quantitative Chemical Ranking for
2
Sodium Iodide Symporter (NIS) Inhibitors in ToxCast Phase I Chemical
3
Library
4
Jun Wang1,2, Daniel R. Hallinger2, Ashley S. Murr2, Angela R. Buckalew2, Steven O. Simmons3,
5
Susan C. Laws2,*, Tammy E. Stoker2,*
6
1
7
37831, USA; 2Endocrine Toxicology Branch, Toxicity Assessment Division, National Health
8
and Environmental Effects Research Laboratory, Office of Research and Development, U.S.
9
Environmental Protection Agency, Research Triangle Park, NC, 27711, USA; 3National Center
10
for Computational Toxicology, Office of Research and Development, U.S. Environmental
11
Protection Agency, Research Triangle Park, NC, 27711, USA
12
*Corresponding authors:
[email protected] (Phone: 919-541-0173 Fax: 919-541-5138) and
13
[email protected] (Phone: 919-541-2783 Fax: 919-541-5138)
Oak Ridge Institute for Science and Education, U.S. Department of Energy, Oak Ridge, TN
1 ACS Paragon Plus Environment
Environmental Science & Technology
14
Abstract
15
Thyroid uptake of iodide via the sodium-iodide symporter (NIS) is the first step in the
16
biosynthesis of thyroid hormones that are critical for health and development in humans and
17
wildlife. Despite having long been a known target of endocrine disrupting chemicals such as
18
perchlorate, information regarding NIS inhibition activity is still unavailable for the vast majority
19
of environmental chemicals. This study applied a previously validated high-throughput approach
20
to screen for NIS inhibitors in the ToxCast phase I library, representing 293 important
21
environmental chemicals. Here 310 blinded samples were screened in a tiered-approach by an
22
initial single-concentration (100µM) radioactive-iodide uptake (RAIU) assay, followed with 169
23
samples further evaluated in multi-concentration (0.001µM-100µM) testing in parallel RAIU and
24
cell viability assays. A novel chemical ranking system that incorporates multi-concentration
25
RAIU and cytotoxicity responses was also developed as a standardized method for chemical
26
prioritization in current and future screenings. Representative chemical responses and thyroid
27
effects of high-ranking chemicals are further discussed. This study significantly expands current
28
knowledge of NIS inhibition potentials in environmental chemicals, and provides critical support
29
to U.S.EPA’s Endocrine Disruptor Screening Program (EDSP) initiative to expand coverage of
30
thyroid molecular targets as well as the development of thyroid adverse outcome pathways
31
(AOPs).
2 ACS Paragon Plus Environment
Page 2 of 34
Page 3 of 34
Environmental Science & Technology
32
Introduction
33
The presence of endocrine disrupting chemicals (EDCs) in the environment continues to be a top
34
public health concern not only for effects on the regulation of androgen and estrogen pathways,
35
but also thyroid hormone homeostasis. Thyroid hormones [TH; i.e., thyroxine (T4) and
36
triiodothyronine (T3)] regulate an array of physiological processes that are essential for
37
metabolism, cardiovascular function, bone maintenance, as well as fetal and post-natal
38
neurodevelopment.1-3 Over the past two decades, a number of structurally diverse xenobiotics
39
have been shown to interfere with TH homeostasis and result in physiological and morphological
40
perturbations.4-7 Further work in this area has identified multiple molecular targets of chemical-
41
mediated thyroid disruption including the regulation of circulating TH through feedback
42
mechanisms within the hypothalamic-pituitary-thyroid (HPT) axis, TH synthesis and secretion,
43
TH distribution and transport, TH metabolism, and TH receptor binding and action.8-13 These
44
studies demonstrate a need for a better understanding of the structural characteristics of
45
chemicals with thyroid disrupting activity, the potential for environmental exposures, and the
46
possible health risks to humans and wildlife. To address these concerns, the U.S. EPA’s
47
Endocrine Disruptor Screening Program (EDSP21, https://www.epa.gov/endocrine-disruption) in
48
concert with the U.S. EPA’s Office of Research and Development (ORD) recently expanded the
49
coverage of molecular targets for thyroid disruption by developing and implementing the use of
50
high-throughput screening (HTS) assays to identify inhibitors of TH synthesis (e.g.,
51
sodium/iodide symporter (NIS) and thyroid peroxidase) and T4 metabolism (e.g., deiodinases).
52
Methods for a thyroid peroxidase assay were previously described and used to screen chemicals
53
in ToxCast chemical libraries.14, 15 In addition, our laboratory recently demonstrated the utility
3 ACS Paragon Plus Environment
Environmental Science & Technology
Page 4 of 34
54
of another HTS assay to detect chemicals that disrupt NIS-mediated transport of extracellular
55
iodide across the cellular membrane.16
56
NIS is a glycoprotein with 13 transmembrane helices that actively transports iodide (I-) into
57
thyroid follicular cells.17-19 This activity relies on the Na+ electrochemical gradient maintained by
58
Na+/K+ ATPases with an electrogenic stoichiometry of 2 Na+ per I-.20,
59
physiological conditions, NIS can concentrate iodide into the thyroid gland 20- to 40-fold greater
60
than serum levels.22 Malfunction of the NIS protein is known to disrupt TH homeostasis; patients
61
with various gene mutations leading to aberrant NIS protein have been diagnosed with
62
hypothyroidism.23 NIS activity can also be disrupted by xenobiotic chemicals. A well-known
63
example is perchlorate (ClO4−), which is a widespread environmental contaminant that has been
64
detected in soil, ground/surface water, food, and consumer products across the world24-28.
65
Perchlorate, along with several other environmental anions including thiocyanate (SCN−) and
66
nitrate (NO3−), have been demonstrated to competitively inhibit I- uptake by NIS and disrupt TH
67
synthesis in humans
68
known target of endocrine disruption, studies on NIS inhibition activity have been largely
69
restricted to only a few environmental chemicals. However, two studies have recently identified
70
triclosan, triclocarban, BDE-47, bisphenol A34, and several small drug-like organic molecules as
71
NIS-mediated iodide transport inhibitors35, suggesting that NIS inhibition activity extends
72
beyond traditionally known anions to the more complex organic compounds. Considering the
73
pivotal role of NIS in the thyroid hormone system, there is an urgent need to expand the
74
knowledge of NIS inhibition potentials to a broader range of environmental chemicals.
29-31
21
Under normal
and multiple species of vertebrates32, 33. Despite having long been a
4 ACS Paragon Plus Environment
Page 5 of 34
Environmental Science & Technology
75
In this study, we applied a previously validated HTS approach to screen the 293 ToxCast
76
(https://www.epa.gov/chemical-research/toxicity-forecasting) phase I chemicals and developed a
77
specific chemical ranking approach to assist with chemical prioritization in current and future
78
screenings. ToxCast phase I library contains a diverse collection of environmental chemicals
79
(mostly pesticides and antimicrobials) that are currently under the U.S. EPA’s regulatory
80
purview36, as well as a subset of chemicals that have been extensively tested in the Agency’s
81
EDSP Tier I Screening Battery.37 This study significantly expands our current knowledge of
82
environmentally-relevant chemicals with NIS inhibition potential and contributes to a broader
83
understanding of thyroid disruptive mechanisms. In addition, the results support the development
84
of adverse outcome pathways (https://aopwiki.org/) related to thyroid hormone disruption and
85
Organization for Economic Co-operation and Development (OECD) test guidelines38 to assess
86
thyroid disruptive chemicals. With an accompanied R package made publicly available, the
87
newly developed chemical ranking approach could also be extended as a standardized method to
88
prioritize chemicals for other HTS studies.
5 ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 34
89
Materials and Methods
90
Chemicals
91
Chemical names, CAS numbers, and maximum concentrations tested are shown in Table S1 (test
92
chemicals) and Table S2 (assay controls). All control chemicals were initially solubilized in
93
DMSO (EMD Millipore Corp., Darmstadt, Germany) at 20mM and included sodium perchlorate
94
(NaClO4; RAIU assay positive control), sodium nitrate (NaNO3; RAIU assay EC80 control),
95
sodium thiocyanate (NaSCN; RAIU assay EC20 control), 2,4-dichlorophenoxyacetic acid (2,4-D;
96
RAIU assay negative control), and 2,3-dichloro-1,4-napthoquinone (DCNQ; cell viability assay
97
positive control) (Sigma Aldrich, St. Louis, MO).
98
Chemicals of the ToxCast phase I_v2 library were obtained from the National Center of
99
Computational Toxicology (NCCT), US EPA, Research Triangle Park, NC, USA.36
Each
100
chemical was solubilized in DMSO (≤ 20mM, Table S1) and provided as 310 blinded samples in
101
five 96-well plates (62 samples per plate) (Evotec Inc., South San Francisco, CA). Stock
102
chemical plates were visually examined under microscope to check for solubility/precipitate in
103
each well (Table S1). Of the 310 samples, there were 293 unique chemicals; the remaining 17
104
served as internal quality control replicates of 12 chemicals randomly distributed among the five
105
plates. All chemicals were transferred to bioassay plates using a BioMek FX Automated
106
Laboratory Workstation (Beckman Coulter, Indianapolis, IN) equipped with a stainless steel high
107
density replicating (HDR) tool (96-pin) to transport 0.35µL per well for the RAIU and cell
108
viability assays.
6 ACS Paragon Plus Environment
Page 7 of 34
Environmental Science & Technology
109
Radioactive Iodide Uptake (RAIU) and Cell Viability Assays
110
RAIU and cell viability assays were conducted with low passage hNIS-HEK293T-EPA cells (
50%, absolute EC50
163
(absEC50) were also reported. absEC50 was determined as the log concentration where the
164
modeled activity equals 50% of control activity. absEC80 was determined as the log
165
concentration where the modeled activity equals 80% of control activity.
166
For chemicals that demonstrated significant cytotoxicity (exceeding 3bMAD of cell viability
167
assay, 17.7%), the concentration where a significant reduction in cell viability for each chemical
168
was determined and referred to as the cytotox-point. The cytotox-point is the equivalent of
169
absEC82.3, the log concentration where the modeled activity equals the cutoff value for
170
significant cytotoxicity (82.3% of control viability).
( ) |
, where σpos and µpos are the standard deviation and
9 ACS Paragon Plus Environment
Environmental Science & Technology
171
Chemical ranking score
172
To prioritize the chemicals for potential NIS inhibition activity and further evaluation, a new
173
scoring system was developed based on two metrics that take into account the confounding
174
impact of cytotoxicity on identifying RAIU inhibition activity: 1) toxicity-adjusted area (TAA)
175
and 2) the difference of median responses of RAIU and cell viability at maximum tested
176
concentration (Median-Difference) (Figure1). Ranking analysis was only performed if a
177
chemical produced significant RAIU inhibition in multi-concentration screening. To obtain TAA
178
(gray stripe area illustrated in Figure 1), RAIU inhibition area and cytotoxicity area were first
179
calculated. RAIU inhibition area was defined by the RAIU 3bMAD (23.8%) significant
180
threshold horizontal line (top border), maximum concentration vertical line (right border), and
181
the RAIU dose-response curve. Cytotoxicity area was defined with the same top and right
182
borders and the cell viability dose-response curve. TAA was then obtained by subtracting the
183
cytotoxicity area from the RAIU inhibition area. Therefore, the numeric value of TAA is
184
penalized when a chemical demonstrates strong cytotoxicity. Median-Difference was calculated
185
using the median of cell viability responses minus the median of RAIU responses at the
186
maximum tested concentration (usually 100 µM). Larger Median-Difference values represent
187
larger separations between RAIU and cell viability.
188
To rank test chemicals, the well-documented NIS inhibitor, NaClO4, was chosen as the reference
189
chemical to normalize the TAA and Median-Difference of each test chemical. Specifically, the
190
TAA and Median-Difference values of the NaClO4 positive control included on each of the 54
191
multi-concentration testing plates were first calculated to obtain the median of NaClO4 TAA and
192
Median-Difference (150.03 and 95.67 respectively). The TAA and Median-Difference of test
10 ACS Paragon Plus Environment
Page 10 of 34
Page 11 of 34
Environmental Science & Technology
193
chemicals were normalized as the percentage of the median NaClO4 TAA and Median-
194
Difference separately and then summed to obtain a chemical ranking score. The ranking score of
195
200 represents the potency of the reference NaClO4 (Figure 1A). Figure 1B and 1C show the
196
dose-response of two test chemicals with lower level ranking scores due to less RAIU inhibition
197
and increased cytotoxicity level.
198
Functions for the calculation of the ranking score, TAA, and Median-Difference, along with
199
dose-response modeling and visualization are made available in the R package ‘toxplot’
200
(https://cran.r-project.org/package=toxplot).
11 ACS Paragon Plus Environment
Environmental Science & Technology
201
Results and Discussion
202
Assay performance and quality control
203
Performance of the RAIU and cell viability assays was monitored for each assay plate through
204
CV of DMSO vehicle control, Z’ and AC50 for the positive control chemicals NaClO4/DCNQ
205
(Table S3). Fifteen plates in single-concentration RAIU assay and 54 plates were used in multi-
206
concentration RAIU and cell viability assays, respectively. Both RAIU and cell viability assays
207
performed well with excellent dynamic range, reproducibility, and reliability during both single
208
and multi-concentration screening. CV of DMSO were ≤ 11.5% and the standard deviation of
209
control AC50 were ≤0.13 (logM). The Z’ scores were also consistently above 0.64 in RAIU
210
screening. Additional controls including NaNO3, NaSCN, and 2,4-D (RAIU assay EC80, EC20
211
and negative controls) also showed great consistency with expected responses in the RAIU and
212
cell viability assays with small variances across the entire screening (Table S4).
213
The phase I_v2 chemical library consisted of 310 blinded samples that included 293 unique
214
chemicals. Twelve of these chemicals were internally replicated (7 chemicals replicated twice
215
and 5 chemicals replicated three times) to assess assay reproducibility. The robustness of the
216
RAIU assay is shown in Figure S2a, where all 12 replicated chemicals, excluding bisphenol A,
217
produced highly reproducible results. Of the 12 replicated chemicals, 10 exceeded the 20%
218
inhibition threshold and were subsequently tested in multi-concentration screening. The AC50 for
219
each replicate of the 10 chemicals was calculated to assess the reproducibility of the RAIU assay
220
in the multi-concentration setting (Figure S2b, Table S5). These 10 internal replicate samples
12 ACS Paragon Plus Environment
Page 12 of 34
Page 13 of 34
Environmental Science & Technology
221
had highly reproducible AC50 (logM) values with the maximum variation range less than 0.25
222
(logM).
223
Single-concentration screening
224
To facilitate the screening process, the 310 blinded samples were first tested in the RAIU assay
225
at their maximum permissible concentration (typically 100µM, Table S1) to select potentially
226
active compounds for multi-concentration evaluation. Single-concentration RAIU screening
227
results for all samples were ordered by increasing inhibitory median responses and plotted as
228
median and maximum/minimum responses for each sample (Figure S3). These samples
229
demonstrated a wide range of iodide uptake inhibition activity in the assay with median
230
responses ranging from 2.8% to 116.6% of maximum iodide uptake relative to the DMSO
231
control. Of the 310 samples, 169 samples (54.5%) produced over the 20% inhibition threshold
232
and were subjected to multi-concentration testing (Table S1).
233
All chemical samples were also tested with Sandell-Kolthoff reactions41 to determine the
234
presence of any contaminant iodide that may lead to potential false positive RAIU inhibition (see
235
SI Part I). Only two chemicals (iodosulfuron-methyl-sodium and 3-iodo-2-propynyl-N-
236
butylcarbamate) produced positive SK reactions, suggesting that overall iodide contamination
237
was not introducing false positive results in the RAIU assay.
238
Multi-concentration testing
239
The 169 blinded chemical samples with ≥20% inhibition in single-concentration RAIU screening
240
were further tested at 6 concentrations (0.001µM – 100µM) in parallel RAIU and cell viability
13 ACS Paragon Plus Environment
Environmental Science & Technology
241
assays. Dose-response curves for test chemicals and positive controls (NaClO4 and DCNQ) are
242
available in SI. The significant activity threshold (3bMAD, calculated across all 54 assay plates
243
in multi-concentration screening) for the cell viability and RAIU assays were 17.7% and 23.8%,
244
respectively. Among the 169 tested samples, 137 showed significant RAIU inhibition activity
245
(>23.8% inhibition). Within the 137 samples with significant RAIU inhibition activity, 26
246
samples exhibited no change in cell viability, while 111 displayed a cytotoxic response at one or
247
more concentrations tested.
248
In reporting the RAIU and cell viability responses of each test chemical, several activity metrics,
249
including AC50, absEC50, and cytotox-point (Table 1 and Table S1), were used together to aid in
250
the interpretation of the assay results. AC50 and absEC50 indicate the chemical potency from two
251
different perspectives, as they differ in the inhibition activity level each represents. absEC50 is
252
the concentration that causes 50% inhibition, while AC50 is the concentration that triggers half-
253
maximal inhibition, which is 100, 17 had ranking
280
scores between 50 and 100, and 116 had ranking scores KClO4>NaClO4>NaBF4>NaSCN) is also in full
286
agreement with previously reported potencies.16
287
Cytotoxicity information is critical for interpreting the observed inhibition in the RAIU assay. A
288
traditional toxicological approach would eliminate testing results for all concentrations with
289
significant cytotoxicity. If using this approach, 46 chemicals in this study produced significant
290
RAIU inhibition without significant cytotoxicity at one or two concentrations (Table 1, S1 and
291
Figure 2). Among the 46 chemicals, 5 of them produced >50% RAIU inhibition at concentrations
292
where no cytotoxicity was observed. However, simply excluding RAIU data points where
293
significant cytotoxicity is observed may also cause high false negative rate that is undesirable in
294
HTS, especially for chemicals that only produced borderline cytotoxicity. Although it is
295
impossible to completely dissociate RAIU and cytotoxic activities, a chemical with stronger
296
cytotoxicity at lower concentrations is less likely to be an RAIU inhibitor than one with less or
297
no cytotoxicity, given the same level of observed RAIU inhibition. The new ranking system was
298
designed with this rationale to provide a continuous metric that adjusts cytotoxic chemicals
299
without completely excluding them, thereby reducing the chance of false negatives. For example,
300
a chemical such as triphenyltin hydroxide (Figure 3a) that had significant and strong RAIU
301
inhibition but with low cytotoxicity (at 1E-5M and 1E-4M) will not only be included, but also
302
receive a high ranking score.
303
As described previously, the ranking score incorporates cytotoxicity responses when calculating
304
the two underlying metrics, TAA and Median-Difference. Unlike other single-point metrics such
305
as AC50, the TAA incorporates levels of RAIU inhibition and cytotoxicity at all tested
306
concentrations and the shape of the dose-response curves. The addition of Median-Difference as 16 ACS Paragon Plus Environment
Page 16 of 34
Page 17 of 34
Environmental Science & Technology
307
a component in the ranking score further promotes chemicals with high RAIU inhibition and low
308
cytotoxicity at the highest concentration so that chemicals like PFOS, oxyfluorfen, and
309
cypronidil (Figure 3a) were ranked higher than using TAA alone. These features of the ranking
310
system enhance the quality of chemical prioritization and it produced a vastly different ranking
311
compared to using the AC50 or absEC50 of RAIU response, as ranking scores had only moderate
312
correlation with AC50 or absEC50 for the 137 chemicals (Spearman’s rank correlation, r = 0.46
313
and 0.41, p 100, 17 having ranking score between 50 and 100, and 116 having ranking score
606