High-Throughput Screening and Quantitative Chemical Ranking for

1 day ago - Thyroid uptake of iodide via the sodium-iodide symporter (NIS) is the first step in the biosynthesis of thyroid hormones that are critical...
7 downloads 10 Views 1MB Size
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