Side Chains of Parabens Modulate Antiandrogenic Activity: In Vitro

May 3, 2017 - The antiandrogenic activity of parabens with alkyloxyl side chains decreases as the side chain length ... Pest Management Science 2017 2...
0 downloads 3 Views 1MB Size
Subscriber access provided by Binghamton University | Libraries

Article

Side Chains of Parabens Modulate Antiandrogenic Activity: In Vitro and Molecular Docking Studies Keke Ding, Xiaotian Kong, Jingpeng Wang, Liping Lu, Wenfang Zhou, Tingjie Zhan, Chunlong Zhang, and Shulin Zhuang Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 03 May 2017 Downloaded from http://pubs.acs.org on May 4, 2017

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 free 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 accessible to all readers and 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.

Environmental Science & Technology 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 37

Environmental Science & Technology

1

Side Chains of Parabens Modulate Antiandrogenic Activity: In

2

Vitro and Molecular Docking Studies

3 4

Keke Ding†, Xiaotian Kong‡, Jingpeng Wang†, Liping Lu†, Wenfang Zhou§, Tingjie

5

Zhan†, Chunlong Zhang⊥, Shulin Zhuang†, *

6 7



8

310058, China

9



College of Environmental and Resource Sciences, Zhejiang University, Hangzhou

Institute of Functional Nano and Soft Materials (FUNSOM), Soochow University,

10

Suzhou, Jiangsu 215123, P. R. China.

11

§

College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China



Department of Biological and Environmental Sciences, University of Houston-Clear

12 13

Lake, 2700 Bay Area Blvd., Houston, Texas 77058, USA.

14

15 16 17 18

* Corresponding author: College of Environmental and Resource Sciences, Zhejiang

19

University, Hangzhou 310058, China.

20

Email address: [email protected] (S Zhuang)

1

ACS Paragon Plus Environment

Environmental Science & Technology

21

ABSTRACT: Parabens have been widely used in packaged foods, pharmaceuticals

22

and personal-care products. Considering their potential hydrolysis, we herein

23

investigated structural features leading to the disruption of human androgen receptor

24

(AR) and whether hydrolysis could alleviate such effects using the recombinant yeast

25

two-hybrid assay. Parabens with an aryloxy side chain such as benzyl paraben and

26

phenyl paraben have the strongest antiandrogenic activity. The antiandrogenic activity

27

of parabens with alkyloxyl side chains decreases as the side chain length increases

28

from 1 to 4, and no antiandrogenic effect occurred for heptyl, octyl and dodecyl

29

parabens with the number of alkoxyl carbon atoms longer than 7. The antiandrogenic

30

activity of parabens correlates significantly with their binding energies (R2 = 0.84, p =

31

0.01) and were completely diminished after the hydrolysis, particularly for parabens

32

with aryloxy side chains. The Km for the hydrolysis of parabens with aromatic moiety

33

side chain is one order of magnitude higher than that of the parabens with alkyl side

34

chains. Both in vitro and in silico data, for the first time, suggest parabens with

35

aromatic side chains are less prone to hydrolysis. Our results provide an insight into

36

risk of various paraben and considerations for design of new paraben-related

37

substitutes.

38 39 40 41 42

2

ACS Paragon Plus Environment

Page 2 of 37

Page 3 of 37

Environmental Science & Technology

43

 INTRODUCTION

44

Parabens (PBs) are 4-hydroxybenzoic acid (HbA) esters with aryloxy or alkoxyl

45

moieties in different lengths as side chains (Table 1). They are widely used as

46

anti-microbial preservatives in packaged foods, pharmaceuticals and personal-care

47

products.1-3 The concentrations of total parabens in 215 cosmetic products were as

48

high as 0.01% - 0.87% (w/w),4 with the concentration of up to 0.4% for a single

49

paraben.5 Parabens are commonplace in our daily exposure and they can be absorbed

50

through human skin (the dermal intake of total parabens is estimated to be 31.0

51

µg/kg-bw/day for adult female in the U.S.)2 and gastrointestinal tract.6 Parabens can

52

also be bioaccumulated and biomagnified through food chains due to the relatively

53

high lipophilicity.7 The extensive use of parabens frequently leads to the detection of

54

residues in multiple environmental matrices, including surface waters, soils,

55

sediments, indoor air and dust. 8-10 Parabens were also frequently detected in breast

56

milk, urine and serum with concentrations at µg/L levels.6, 11-14

57

The ubiquitous presence of parabens in the environment has raised public concern

58

about whether they potentially cause adverse reproductive effects and their safety has

59

been widely discussed.15-17 Exposure to parabens was reported to cause estrogenic

60

disrupting effects.18-23 Ethyl- and butyl parabens caused a significant increase in

61

progesterone formation in the adrenal H295R steroidogenesis assay and interfered

62

with the transport of cholesterol to the mitochondria.24 There has been a growing

63

concern in recent years about the health risk of parabens to the androgen hormone

64

system.25-27 Parabens induced an adverse effect on the male reproductive system

3

ACS Paragon Plus Environment

Environmental Science & Technology

65

through reduced testosterone concentrations and sperm production in rats,28

66

potentially by influencing the biosynthesis, transport, and metabolism of steroid

67

hormone.24 Evaluation of the androgen receptor (AR) disruption is thus of

68

significance for better risk assessment of various paraben analogues.

69

There are two major types of side chains of parabens including the alkyl and

70

aromatic chains (Table 1). Whether the length of paraben side chains will affect their

71

potential health risk to androgen system remains largely unexplored. Structural

72

analogues of many environmental contaminants have been examined to determine

73

how a minor change in structure governs their human health risk.24, 29-33 For example,

74

minor changes in the bridged carbon of several bisphenol A substitutes correlated well

75

with experimentally measured toxicities and hormonal activities.29, 34, 35 The binding

76

affinity of twelve perfluorinated carboxylic acids to human liver fatty acid binding

77

protein increased significantly with their carbon number from 4 to 11, and decreased

78

slightly when the number was over 11.36 To date studies relevant to the risk of a

79

homolog series of parabens are limited. 37, 38 The side chains of parabens can be

80

removed via hydrolysis in human tissues and multiple matrices (air, water, soil, or

81

sediment), possibly affecting their environmental risks. 8, 39, 40 Till now, much still

82

remains unknown as to whether and how the hydrolysis of side chains of parabens

83

disrupts human AR. The influence of the sizes of side chains on their hydrolysis and

84

the associated kinetic parameters are essential to the risk assessment of parabens.

85

Such information on the effect of different side chains before and after the hydrolysis

86

is of significance for better health risk assessment of parabens.

4

ACS Paragon Plus Environment

Page 4 of 37

Page 5 of 37

Environmental Science & Technology

87

In the present study, we investigated the effect of parabens with different sizes of

88

alkyl and aromatic side chains on their induction of AR disruption and the impact of

89

hydrolysis on this effect through an in vitro removal of side chains via

90

esterase-mediated hydrolysis. Molecular docking was conducted to elucidate the

91

contribution of different side chains of parabens to their interactions with AR ligand

92

binding domain (LBD). We revealed the significant role of aromatic moiety toward

93

the AR disrupting effect of parabens. Our results provided a structure activity

94

relationship (SAR) to decipher how hydrolytic metabolism and side chain change of

95

parabens can exert considerable difference in activity and potential risk. Such

96

information is crucial to assess the risk of various paraben-containing consumer

97

products and design environmentally benign paraben substitutes by considering the

98

multiple adverse effects induced by various local molecular moieties.

99

 MATERIALS AND METHODS

100

Chemicals. Parabens (purity ≥ 98%) with different side chain moieties (Table 1)

101

were purchased from J&K Chemicals Corporation (Shanghai, China). Dimethyl

102

sulfoxide (DMSO, 99.5%) and dihydrotestosterone (DHT, 98%) were purchased from

103

Sigma Chemical (St. Louis, MO, USA). Porcine liver esterase (PLE, Catalogue:

104

E3019-20KU) was also purchased from Sigma Chemical (St. Louis, MO, USA).

105

SD/-Leu/-Trp medium was purchased from Mobitec Company (Catalogue: 4823-6).

106

The constituents in mobile phase including water, methanol and acetic acid were

107

chromatographic grade. All other chemicals were analytical grade. All parabens were

108

dissolved in DMSO and diluted using deionized water (18.2 MΩ, Millipore, Bedford, 5

ACS Paragon Plus Environment

Environmental Science & Technology

109

MA) for activity evaluation. All prepared solutions were stored at 4 oC in centrifuge

110

tubes. Recombinant Yeast Two-Hybrid Assay. The recombinant two-hybrid yeast

111 112

system employed yeast cells transformed with the human AR plasmid, coactive

113

plasmid, and the reporter gene expressing β-galactosidase. pGBT9 AR plasmid was

114

provided by Dr. Erik Jan Dubbink (Erasmus University, Holland) and the yeast

115

transformed with the pGBT9AR plasmid was provided by Prof. Zijian Wang

116

(Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences,

117

Beijing, China). The transfection efficiencies of the yeast and positive control samples

118

were greater than 4.5 × 104 cfu / µg. The blank control was 0, and there was a

119

significant difference between the positive and negative controls (p PnP > MeP > EtP > PrP > BuP >

214

HbA. This agrees with the previous report that MeP, PrP and BuP inhibited

215

testosterone-induced transcriptional activity by 40%, 33% and 19%, respectively, at a

216

concentration of 10 µM using human embryonic kidney (HEK) 293 cells that lack

217

critical steroid metabolizing enzymes.47 Although parabens have no predominant

218

antiandrogenic effect at the low concentration, their bioaccumulation and

10

ACS Paragon Plus Environment

Page 10 of 37

Page 11 of 37

Environmental Science & Technology

219

biomagnification through food chains may lead to the higher concentration in humans,

220

thereby inducing potentially higher adverse health risks.

221

With the 10 analogues of parabens tested, some structural dependency of the

222

antiandrogenic activity could be inferred. The antiandrogenic activity potency could

223

be associated with the length of the ester alkyl side chain of parabens. The AR

224

disrupting potency decreased with an increased alkyl side chain length from 1 to 4

225

(MeP > EtP > PrP > BuP), whereas the lipophilicity and their penetration through

226

human epidermis-dermis layers are in the reverse order (MeP < EtP < PrP < BuP).48

227

No significant AR disrupting potency was found for HeP, OcP and DdP with the alkyl

228

side chain length equals to 7, 8 and 12, respectively. This is an interesting observation

229

for parabens with alkyl side chains, which implies that the adverse effects of MeP

230

with its highest antiandrogenic activity can be reduced to the greatest extent by the

231

nature of its low lipophilicity and dermal penetration. The parabens with aromatic

232

side chain have distinct disrupting effect in comparison with parabens with alkyl side

233

chain. Side chains with an aromatic moiety showed a stronger AR disrupting potency

234

compared with the alkyl side chains. These parabens with the aromatic side chain are

235

less abundant (e.g., 15% in baby care products) 2 than MeP, EtP, and PrP but their

236

higher logKow values might potentially aggravate their adverse health hazards. This

237

observed structural dependence agrees with an acute aquatic toxicity assessment of

238

parabens toward Vibrio fischeri, in which the acute toxicity value (EC50) of benzyl

239

paraben was 3.8 µg/L (0.017 µM), 2 to 3 orders of magnitude lower than that of

240

parabens with alkyl side chains.49 Therefore, the different side chains of parabens

11

ACS Paragon Plus Environment

Environmental Science & Technology

241

have significant influence on its androgenic disrupting effect and the aromatic moiety

242

of the side chain contributes predominantly to the observed toxicity.

243

Side Chains Affect the Molecular Recognition toward AR LBD. To elucidate

244

why the side chain of different sizes affects the AR disruption effect, the molecular

245

docking was applied to examine their interactions with AR LBD at the atomic level.

246

The ligand in the crystal structure (PDB entry: 1Z95) was extracted and re-docked

247

into the binding pocket. The RMSD between the predicted binding pose and the

248

respective experimental structure was calculated.50-54 The low RMSD (0.49 Å)

249

suggests Glide docking can accurately reproduce the experimental binding poses. A

250

series of parabens and DHT were then all automatically docked into the binding

251

cavity of AR LBD. The docking free energy and three decomposed energy

252

components including hydrophobic interactions (∆E ), electrostatic interactions

253

(∆E

 ) and van der Waals interactions (∆E ) are listed in Table S1. DHT has the

254

smallest docking score compared with the parabens, indicating its strongest binding

255

affinities with AR LBD. The energy components of DHT to AR LBD including

256

∆E , ∆E

 and ∆E also have the most negative values. Besides DHT, the

257

binding potency of parabens to AR LBD was in the order of BzP > PnP > MeP > EtP >

258

PrP > BuP > HbA > HeP > OcP. The docking scores appear to be relevant to the

259

antiandrogen of parabens corresponding to the structural variations of different side

260

chains. Consequently, the correlation between Glide docking energy and the

261

maximum antiandrogenic activity of parabens were performed. As shown in Figure 2,

262

there exists a significant linear correlation between the docking energy and the

12

ACS Paragon Plus Environment

Page 12 of 37

Page 13 of 37

Environmental Science & Technology

263

antiandrogenic activity (R2 = 0.84, p = 0.01). The parabens with the smaller docking

264

score have the stronger antiandrogenic activities, further indicating the agreement of

265

our in vitro and in silico data based on yeast assay and molecular docking.

266

We examined the surrounding amino acids and binding sites of parabens to AR

267

LBD. Before we conducted the molecular docking of parabens to AR LBD, the

268

natural androgen hormone DHT was first docked into the binding cavity of AR LBD.

269

Our results suggested that two amino acids Thr877 and Asn705 were involved in

270

forming hydrogen bonds between DHT and AR LBD (Figure S2). This agrees with a

271

previous study on the interactions of DHT with AR LBD,55 hence validating our

272

molecular docking approach. As shown in Figure 3, benzyl-, phenyl paraben and the

273

alkoxyl parabens with the number of carbon atom equals to 1 to 4, 7 and 8 were

274

docked into the AR LBD, while the larger side chain (e.g., dodecyl paraben) could not

275

be docked into the AR LBD (data not shown). The reason why dodecyl paraben could

276

not be docked into the AR LBD attributes to the relatively bulkier volume of dodecyl

277

paraben which is too large to fit the binding cavity of AR LBD. The hydrogen bonds

278

were also formed between the hydroxyl oxygen of parabens and two amino acids

279

(Thr877 and Asn705), similar to the binding mode of DHT to AR LBD. The HbA

280

formed the hydrogen bonds with the amino acids Gln711 and Arg752 of AR LBD and

281

did not show any antiandrogenic effect. These results clearly suggested that the

282

parabens capable of binding Thr877 and Asn705 of AR LBD can compete with DHT,

283

thereby causing the antiandrogenic activity. The HbA and OcP formed hydrogen

284

bonds with Arg752 and Gln711, but did not show any antiandrogenic activity. This

13

ACS Paragon Plus Environment

Environmental Science & Technology

285

further suggests that Thr877 and Asn705 are the amino acids vital to the induction of

286

the antiandrogenic activity. It is worth noting that heptyl paraben did not show any

287

antiandrogenic activity. Nevertheless, heptyl paraben was docked into the AR LBD

288

and formed the hydrogen bond with Thr877 and Asn705. This could be explained by

289

the fact that the yeast two-hybrid system is not sufficiently sensitive to respond to the

290

very weak AR disrupting effect of heptyl paraben.

291

The Hydrolysis Reduced Antiandrogenic Activities. It was reported that the

292

metabolism of parabens by carboxylesterases in rat liver microsomes markedly

293

reduced the estrogenic activity.56 As important Phase I metabolizing enzymes, the

294

esterase enzymes play an important role in the hydrolysis of various exogenous

295

contaminants. The porcine liver esterase was chosen in this study to determine

296

whether the hydrolysis of parabens will influence the potential for androgen

297

disruption in vitro. Six parent parabens with obvious AR disrupting effects were

298

subjected to hydrolysis, including MeP, EtP, PrP, BuP, PnP and BzP. In the absence of

299

porcine liver esterase, the antiandrogenic activity remained unchanged when

300

incubating parent parabens with the incubation solutions, demonstrating the

301

incubation solutions did not influence AR disruption. In the presence of porcine liver

302

esterase, antiandrogenic activities were tested following the hydrolysis of six parabens

303

at 5, 25 and 50 µM by comparing antiandrogenic activities before and after the

304

hydrolysis. As shown in Figure 4, the antiandrogenic activities of parent parabens

305

(benzyl, phenyl, methyl, ethyl, propyl and butyl paraben) at three concentrations of 5,

306

25 and 50 µM were all significantly decreased after the incubation for 2 h. The

14

ACS Paragon Plus Environment

Page 14 of 37

Page 15 of 37

Environmental Science & Technology

307

metabolic products from the hydrolysis of six parent parabens did not show any

308

noticeable AR inhibition, indicating a very weak antiandrogenic activity and a

309

significant hydrolysis deactivation. A similar metabolic deactivation of paraben

310

toward ERα and ERβ was confirmed by incubating butyl paraben with rat liver

311

microsomes.56 Besides hydrolysis, parabens can also be metabolized via other

312

pathways such as photochemical transformation, phase I and phase II metabolism.

313

Parabens were reportedly detected in the conjugated form

314

(glucuronidation/sulfonation) in vivo and their activities were markedly decreased.12

315

The androgenic disrupting effect of other forms of paraben metabolites should be

316

considered in future study.

317

The Hydrolysis Products Identified by HPLC-UV and UPLC-MS/MS. To

318

further identify the metabolites of parabens after hydrolysis, the products of parabens

319

were firstly qualified using HPLC-UV and UPLC-MS/MS. Representative

320

chromatograms of parent parabens, including MeP, EtP, PrP, BuP, BzP and PnP, are

321

shown in Figure S3 (A) - (F), and the retention times of the corresponding parabens

322

are 10.08, 11.35, 13.89, 18.64, 6.75 and 8.22 min, respectively. After a 2 h incubation

323

with four individual parabens (MeP, EtP, PrP and BuP) and porcine liver esterase at 25

324

o

325

peak at the retention time of around 9.00 min, similar to the chromatogram of a

326

standard containing HbA under the same chromatographic condition. The

327

corresponding UV-Vis absorbance spectrum is characterized by a maximum

328

absorbance wavelength at 254.6 nm, also similar to the UV-Vis spectrum of the

C, the peak of each parent paraben completely disappeared and all presented a new

15

ACS Paragon Plus Environment

Environmental Science & Technology

329

standard HbA. For PnP and BzP, however, after 2 h incubation with porcine liver

330

esterase, the enzymatic hydrolysis products peaks can only be resolved by changing

331

the mobile phase to 65% methanol and 35% acetic acid (0.01%, V/V). As shown in

332

Figure S3 (e) and (f), the peak of hydrolysis products in the presence of porcine liver

333

esterase was observed at the retention time of 8.98 and 9.00 min for the incubation of

334

phenyl and benzyl paraben, respectively. The corresponding UV-Vis spectrum of

335

hydrolysis product also has a maximum absorbance wavelength at 254.6 nm, similar

336

to the UV-Vis spectrum of standard HbA. The chromatograms of six parabens

337

obtained from UPLC-MS/MS were the same to the above HPLC-UV analysis with the

338

metabolite peak at the retention time of around 9.00 min, validating the

339

chromatographic separation of the products. The precursor ion of metabolite was [M -

340

H]- with the m/z around 137 Da (Figure 5), one Da lower than HbA with a molecular

341

weight of 138 Da. The m/z values of corresponding daughter ion were around 93 Da

342

([M - H - 44]-), in agreement with the loss of a 44 Da (CO2). They matched the

343

retention time and MS2 spectra with the standard HbA, indicating that only the ester

344

bond was cleaved for these parabens. The peak at the retention time of about 9.00 min

345

from the hydrolysis of six parent parabens can be identified as HbA.

346

The hydrolysis products of six parabens did not display the antiandrogenic

347

activity at the test concentrations, in line with the antiandrogenic activity of HbA

348

standard (Figure 1). Our multiple lines of evidence indicate the formation of HbA

349

product after the hydrolysis of BzP, PnP, MeP, EtP, PrP and BuP. These results were

350

supported by the reported pathways of paraben hydrolysis using skin microsomes and

16

ACS Paragon Plus Environment

Page 16 of 37

Page 17 of 37

Environmental Science & Technology

351

cytosol from human and minipigs.57 The formation of HbA induced negligible

352

antiandrogenic activities of paraben. The dissipated antiandrogenic activities after the

353

hydrolysis mediated by porcine liver esterase are consistent with the observed

354

cleavage of the side chain of parabens.

355

The Hydrolysis Rate was Side Chain-Dependent. To determine the effects of

356

parabens with different side chains on the rate of hydrolysis mediated by porcine liver

357

esterase, the commonly used Lineweaver-Burk plot between 1/V0 and 1/S was applied

358

to determine the maximum rate of reaction (Vmax) and the specificity of enzyme

359

porcine liver esterase for parabens with different side chains. The Lineweaver-Burk

360

plot displayed a good linear correlation based on the R2 values ranging from 0.9958 to

361

0.9997 (Figure S4). The Km and Vmax calculated from the slopes and intercepts are

362

shown in Table 2. The Vmax values were 0.09 - 0.19 mM/min for alkyl parabens and

363

3.43 - 4.91 mM/min for parabens with aromatic side chains, suggesting the very quick

364

hydrolysis rate mediated by porcine liver esterase for alkyl parabens. This high rate

365

suggests that the rapid metabolism may prevent these parabens from accumulating in

366

biological tissues. Hydrolysis rate in porcine liver esterase may differ from that in

367

human plasma likely due to the different enzymatic concentration and enzyme

368

specificity. The Km values for MeP, EtP, PrP, BuP, PnP and BzP incubated with

369

porcine liver esterase were 0.40 ± 0.01, 0.33 ± 0.02, 0.27 ± 0.02, 0.17 ± 0.04, 6.00 ±

370

0.28 and 8.23 ± 0.15 mM, respectively (Table 2). It was suggested that the Km values

371

decreased with the increasing alkyl side chain length of parabens, indicating the

372

parabens with the longer alkyl side chains have the lower hydrolysis rates. The Km

17

ACS Paragon Plus Environment

Environmental Science & Technology

373

values of parabens with aromatic moiety side chain have over one order of magnitude

374

higher than that of the parabens with alkyl side chains, demonstrating that parabens

375

with alkyl side chains are the preferred substrate compared with the parabens with

376

aromatic side chains for porcine liver esterase. In contrast to parabens with aromatic

377

side chain moieties, parabens with alkyl side chains are more susceptible to hydrolysis

378

by porcine liver esterase. These results strongly indicated that the parabens with an

379

aromatic moiety in their side chains tend to be bioaccumulated more readily and

380

hence pose a greater threat to organisms.

381

 ENVIRONMENTAL IMPLICATIONS

382

The present study for the first time provides both the in vitro and in silico evidence on

383

the structure dependent human AR disrupting effects of parabens. It becomes

384

increasingly clear that certain parabens can exert adverse biological activities and

385

their environmental and human health risk should be of concern. Although direct

386

extrapolation from such in vitro and in silico data should be cautious, the effect

387

concentrations observed in this study should help further delineation of human health

388

exposure. The parabens with aromatic side chain moieties have a stronger AR

389

disrupting potency compared with the alkyl side chains, indicative of the aromatic

390

moiety of parabens as the toxicophore toward AR disruption. Parabens with aromatic

391

moieties in the side chains seem to have slower hydrolysis rates and could lead to a

392

higher bioaccumulation and hence a greater threat to organisms. Insight into the

393

structure dependence of antiandrogenic activities and other biological activities of

394

parabens is valuable to formulate and design safe parabens and their substitutes for 18

ACS Paragon Plus Environment

Page 18 of 37

Page 19 of 37

Environmental Science & Technology

395

use in consumer products. Parabens-containing consumer products can be better

396

formulated by modulating the compositions of the primary parabens (i.e., methyl-,

397

ethyl- and propyl parabens) to ensure their safe uses. Our results will be beneficial to

398

understand the environmental transformation and fate of paraben esters and provide

399

the guideline for the production of paraben-containing products with low human

400

toxicity and high efficiency in preservation as antibacterial agents.

401

 ASSOCIATED CONTENT

402

Supporting Information

403

Table S1: The decomposition of binding free energies. Table S2: The amino acids

404

involved with hydrogen bonds and its binding energy between parabens and AR LBD.

405

Figure S1: Androgenic activity of ten parabens evaluated by two-hybrid yeast assay.

406

Figure S2: The binding mode of DHT to AR LBD. Figure S3: High performance

407

liquid chromatogram of parent parabens and products from hydrolysis mediated by

408

porcine liver esterase. Figure S4: The Lineweaver-Burk plots in describing the

409

kinetics of hydrolysis of parabens incubated with porcine liver esterase. The

410

Supporting Information is available free of charge on the ACS Publications website.

411

 AUTHOR INFORMATION

412

Corresponding Author

413

* E-mail: [email protected]. Phone: +86-571-8898-2344.

414

ORCID

415

Shulin Zhuang: 0000-0002-7774-7239

19

ACS Paragon Plus Environment

Environmental Science & Technology

416

Notes

417

The authors declare no competing financial interest.

418

 ACKNOWLEDGEMENTS

419

The authors thank the National Natural Science Foundation of China (No. 21477113)

420

and the Major Research Plan of the National Natural Science Foundation of China

421

(No. 91643107). CZ acknowledged the Welch Foundation for partial support.

422 423

 REFERENCES

424

1.

425

foodstuffs from the United States. Environ Sci Technol 2013, 47, (8), 3918-3925.

426

2.

427

products from the United States and its implications for human exposure. Environ Sci

428

Technol 2013, 47, (24), 14442-14449.

429

3.

430

toxicity, absorption, esterase and human exposure, and discussion of potential human

431

health risks. J Appl Toxicol 2008, 28, (5), 561-578.

432

4.

433

ethyl-, propyl-, butyl- and benzylparaben in cosmetic products. Contact Dermatitis

434

1995, 32, (1), 28-30.

435

5.

436

ethylparaben, propylparaben, isopropylparaben, butylparaben, isobutylparaben, and

437

benzylparaben as used in cosmetic products. Int J Toxicol 2008, 27, 1-82.

Liao, C.; Liu, F.; Kannan, K., Occurrence of and dietary exposure to parabens in

Guo, Y.; Kannan, K., A survey of phthalates and parabens in personal care

Darbre, P. D.; Harvey, P. W., Paraben esters: review of recent studies of endocrine

Rastogi, S. C.; Schouten, A.; de Kruijf, N.; Weijland, J. W., Contents of methyl-,

Andersen, F. A., Final amended report on the safety assessment of methylparaben,

20

ACS Paragon Plus Environment

Page 20 of 37

Page 21 of 37

Environmental Science & Technology

438

6.

Frederiksen, H.; Jorgensen, N.; Andersson, A. M., Parabens in urine, serum and

439

seminal

440

chromatography-tandem mass spectrometry (LC-MS/MS). J Expo Sci Environ

441

Epidemiol 2011, 21, (3), 262-271.

442

7.

443

Accumulation of Parabens and their Metabolites in Marine Mammals from the United

444

States Coastal Waters. Environ Sci Technol 2015, 49, (20), 12071-12079.

445

8.

446

K., Occurrence and human exposure of p-hydroxybenzoic acid esters (parabens),

447

bisphenol A diglycidyl ether (BADGE), and their hydrolysis products in indoor dust

448

from the United States and three East Asian countries. Environ Sci Technol 2012, 46,

449

(21), 11584-11593.

450

9.

451

studies to human health. Environ Int 2014, 67, 27-42.

452

10. Liao, C.; Lee, S.; Moon, H. B.; Yamashita, N.; Kannan, K., Parabens in sediment

453

and sewage sludge from the United States, Japan, and Korea: spatial distribution and

454

temporal trends. Environ Sci Technol 2013, 47, (19), 10895-10902.

455

11. Calafat, A. M.; Ye, X.; Wong, L. Y.; Bishop, A. M.; Needham, L. L., Urinary

456

concentrations of four parabens in the U.S. population: NHANES 2005-2006. Environ

457

Health Perspect 2010, 118, (5), 679-685.

458

12. Ye, X. Y.; Bishop, A. M.; Reidy, J. A.; Needham, L. L.; Calafat, A. M., Parabens

459

as urinary biomarkers of exposure in humans. Environ Health Persp 2006, 114, (12),

plasma

from

healthy

Danish

men

determined

by

liquid

Xue, J.; Sasaki, N.; Elangovan, M.; Diamond, G.; Kannan, K., Elevated

Wang, L.; Liao, C.; Liu, F.; Wu, Q.; Guo, Y.; Moon, H. B.; Nakata, H.; Kannan,

Bledzka, D.; Gromadzinska, J.; Wasowicz, W., Parabens. From environmental

21

ACS Paragon Plus Environment

Environmental Science & Technology

460

1843-1846.

461

13. Hines, E. P.; Mendola, P.; von Ehrenstein, O. S.; Ye, X.; Calafat, A. M.; Fenton, S.

462

E., Concentrations of environmental phenols and parabens in milk, urine and serum of

463

lactating North Carolina women. Reprod Toxicol 2015, 54, 120-128.

464

14. Barr, L.; Metaxas, G.; Harbach, C. A. J.; Savoy, L. A.; Darbre, P. D.,

465

Measurement of paraben concentrations in human breast tissue at serial locations

466

across the breast from axilla to sternum. J Appl Toxicol 2012, 32, (3), 219-232.

467

15. Soni, M. G.; Carabin, I. G.; Burdock, G. A., Safety assessment of esters of

468

p-hydroxybenzoic acid (parabens). Food Chem Toxicol 2005, 43, (7), 985-1015.

469

16. Soni, M. G.; Burdock, G. A.; Taylor, S. L.; Greenberg, N. A., Safety assessment

470

of propyl paraben: a review of the published literature. Food Chem Toxicol 2001, 39,

471

(6), 513-532.

472

17. Soni, M. G.; Taylor, S. L.; Greenberg, N. A.; Burdock, G. A., Evaluation of the

473

health aspects of methyl paraben: a review of the published literature. Food Chem

474

Toxicol 2002, 40, (10), 1335-1373.

475

18. Alslev, B.; Korsgaard, B.; Bjerregaard, P., Estrogenicity of butylparaben in

476

rainbow trout Oncorhynchus mykiss exposed via food and water. Aquat Toxicol 2005,

477

72, (4), 295-304.

478

19. Pedersen, K. L.; Pedersen, S. N.; Christiansen, L. B.; Korsgaard, B.; Bjerregaard,

479

P., The preservatives ethyl-, propyl- and butylparaben are oestrogenic in an in vivo

480

fish assay. Pharmacol Toxicol 2000, 86, (3), 110-113.

481

20. Hu, Y.; Zhang, Z. B.; Sun, L. B.; Zhu, D. S.; Liu, Q. C.; Jiao, J.; Li, J.; Qi, M. W.,

22

ACS Paragon Plus Environment

Page 22 of 37

Page 23 of 37

Environmental Science & Technology

482

The estrogenic effects of benzylparaben at low doses based on uterotrophic assay in

483

immature SD rats. Food Chem Toxicol 2013, 53, 69-74.

484

21. Sun, L. B.; Yu, T.; Guo, J. L.; Zhang, Z. B.; Hu, Y.; Xiao, X.; Sun, Y. L.; Xiao, H.;

485

Li, J. Y.; Zhu, D. S.; Sai, L. L.; Li, J., The estrogenicity of methylparaben and

486

ethylparaben at doses close to the acceptable daily intake in immature

487

Sprague-Dawley rats. Sci Rep 2016, 6, 25173.

488

22. Zhang, Z. B.; Sun, L. B.; Hu, Y.; Jiao, J.; Hu, J. Y., Inverse antagonist activities of

489

parabens on human oestrogen-related receptor gamma (ERR gamma): In vitro and in

490

silico studies. Toxicol Appl Pharm 2013, 270, (1), 16-22.

491

23. Darbre, P. D.; Byford, J. R.; Shaw, L. E.; Hall, S.; Coldham, N. G.; Pope, G. S.;

492

Sauer, M. J., Oestrogenic activity of benzylparaben. J Appl Toxicol 2003, 23, (1),

493

43-51.

494

24. Taxvig, C.; Vinggaard, A. M.; Hass, U.; Axelstad, M.; Boberg, J.; Hansen, P. R.;

495

Frederiksen, H.; Nellemann, C., Do parabens have the ability to interfere with

496

steroidogenesis? Toxicol Sci 2008, 106, (1), 206-213.

497

25. Ma, D.; Chen, L.; Zhu, X.; Li, F.; Liu, C.; Liu, R., Assessment of combined

498

antiandrogenic effects of binary parabens mixtures in a yeast-based reporter assay.

499

Environ Sci Pollut Res Int 2014, 21, (10), 6482-6494.

500

26. Kolsek, K.; Gobec, M.; Mlinaric Rascan, I.; Sollner Dolenc, M., Screening of

501

bisphenol A, triclosan and paraben analogues as modulators of the glucocorticoid and

502

androgen receptor activities. Toxicol In Vitro 2015, 29, (1), 8-15.

503

27. Satoh, K.; Nonaka, R.; Ohyama, K.; Nagai, F., Androgenic and antiandrogenic

23

ACS Paragon Plus Environment

Environmental Science & Technology

504

effects of alkylphenols and parabens assessed using the reporter gene assay with

505

stably transfected CHO-K1 cells (AR-EcoScreen System). J Health Sci 2005, 51, (5),

506

557-568.

507

28. Oishi, S., Effects of butyl paraben on the male reproductive system in mice. Arch

508

Toxicol 2002, 76, (7), 423-429.

509

29. Delfosse, V.; Grimaldi, M.; Pons, J. L.; Boulahtouf, A.; le Maire, A.; Cavailles, V.;

510

Labesse, G.; Bourguet, W.; Balaguer, P., Structural and mechanistic insights into

511

bisphenols action provide guidelines for risk assessment and discovery of bisphenol A

512

substitutes. P Natl Acad Sci USA 2012, 109, (37), 14930-14935.

513

30. Zhuang, S. L.; Wang, H. F.; Ding, K. K.; Wang, J. Y.; Pan, L. M.; Lu, Y. L.; Liu,

514

Q. J.; Zhang, C. L., Interactions of benzotriazole UV stabilizers with human serum

515

albumin: Atomic insights revealed by biosensors, spectroscopies and molecular

516

dynamics simulations. Chemosphere 2016, 144, 1050-1059.

517

31. Yuan, B.; Fu, J.; Wang, Y.; Jiang, G., Short-chain chlorinated paraffins in soil,

518

paddy seeds (Oryza sativa) and snails (Ampullariidae) in an e-waste dismantling area

519

in China: Homologue group pattern, spatial distribution and risk assessment. Environ

520

Pollut 2017, 220, (Pt A), 608-615.

521

32. Zhang, F.; Zhang, J.; Tong, C.; Chen, Y.; Zhuang, S.; Liu, W., Molecular

522

interactions of benzophenone UV filters with human serum albumin revealed by

523

spectroscopic techniques and molecular modeling. J Hazard Mater 2013, 263 Pt 2,

524

618-26.

525

33. Xiao, R. Y.; Ye, T. T.; Wei, Z. S.; Luo, S.; Yang, Z. H.; Spinney, R., Quantitative

24

ACS Paragon Plus Environment

Page 24 of 37

Page 25 of 37

Environmental Science & Technology

526

Structure-Activity Relationship (QSAR) for the Oxidation of Trace Organic

527

Contaminants by Sulfate Radical. Environ Sci Technol 2015, 49, (22), 13394-13402.

528

34. Zhuang, S.; Zhang, C.; Liu, W., Atomic insights into distinct hormonal activities

529

of Bisphenol A analogues toward PPARgamma and ERalpha receptors. Chem Res

530

Toxicol 2014, 27, (10), 1769-1779.

531

35. Rochester, J. R.; Bolden, A. L., Bisphenol S and F: A Systematic Review and

532

Comparison of the Hormonal Activity of Bisphenol A Substitutes. Environ Health

533

Perspect 2015, 123, (7), 643-650.

534

36. Zhang, L.; Ren, X. M.; Guo, L. H., Structure-based investigation on the

535

interaction of perfluorinated compounds with human liver fatty acid binding protein.

536

Environ Sci Technol 2013, 47, (19), 11293-11301.

537

37. Gao, Y. P.; Ji, Y. M.; Li, G. Y.; An, T. C., Theoretical investigation on the kinetics

538

and mechanisms of hydroxyl radical-induced transformation of parabens and its

539

consequences for toxicity: Influence of alkyl-chain length. Water Res 2016, 91, 77-85.

540

38. Okubo, T.; Yokoyama, Y.; Kano, K.; Kano, I., ER-dependent estrogenic activity

541

of parabens assessed by proliferation of human breast cancer MCF-7 cells and

542

expression of ER alpha and PR. Food and Chem Toxicol 2001, 39, (12), 1225-1232.

543

39. Wang, W.; Kannan, K., Fate of Parabens and Their Metabolites in Two

544

Wastewater Treatment Plants in New York State, United States. Environ Sci Technol

545

2016, 50, (3), 1174-1181.

546

40. Pugazhendhi, D.; Pope, G. S.; Darbre, P. D., Oestrogenic activity of

547

p-hydroxybenzoic acid (common metabolite of paraben esters) and methylparaben in

25

ACS Paragon Plus Environment

Environmental Science & Technology

548

human breast cancer cell lines. J Appl Toxicol 2005, 25, (4), 301-309.

549

41. Li, J.; Li, N.; Ma, M.; Giesy, J. P.; Wang, Z. J., In vitro profiling of the endocrine

550

disrupting potency of organochlorine pesticides. Toxicoly Letts 2008, 183, (1-3),

551

65-71.

552

42. Doesburg, P.; Kuil, C. W.; Berrevoets, C. A.; Steketee, K.; Faber, P. W.; Mulder,

553

E.; Brinkmann, A. O.; Trapman, J., Functional in vivo interaction between the

554

amino-terminal, transactivation domain and the ligand binding domain of the

555

androgen receptor. Biochemistry 1997, 36, (5), 1052-1064.

556

43. Brausch, J. M.; Rand, G. M., A review of personal care products in the aquatic

557

environment: Environmental concentrations and toxicity. Chemosphere 2011, 82, (11),

558

1518-1532.

559

44. Li, J.; Ma, M.; Wang, Z., A two-hybrid yeast assay to quantify the effects of

560

xenobiotics on thyroid hormone-mediated gene expression. Environ Toxicol Chem

561

2008, 27, (1), 159-167.

562

45. Zhuang, S.; Lv, X.; Pan, L.; Lu, L.; Ge, Z.; Wang, J.; Wang, J.; Liu, J.; Liu, W.;

563

Zhang, C., Benzotriazole UV 328 and UV-P showed distinct antiandrogenic activity

564

upon human CYP3A4-mediated biotransformation. Environ Pollut 2017, 220, (Pt A),

565

616-624.

566

46. Lv, X.; Pan, L.; Wang, J.; Lu, L.; Yan, W.; Zhu, Y.; Xu, Y.; Guo, M.; Zhuang, S.,

567

Effects of triazole fungicides on androgenic disruption and CYP3A4 enzyme activity.

568

Environ Pollut 2017, 222, 504-512.

569

47. Chen, J.; Ahn, K. C.; Gee, N. A.; Gee, S. J.; Hammock, B. D.; Lasley, B. L.,

26

ACS Paragon Plus Environment

Page 26 of 37

Page 27 of 37

Environmental Science & Technology

570

Antiandrogenic properties of parabens and other phenolic containing small molecules

571

in personal care products. Toxicol Appl Pharmacol 2007, 221, (3), 278-284.

572

48. El Hussein, S.; Muret, P.; Berard, M.; Makki, S.; Humbert, P., Assessment of

573

principal parabens used in cosmetics after their passage through human

574

epidermis-dermis layers (ex-vivo study). Exp Dermatol 2007, 16, (10), 830-836.

575

49. Terasaki, M.; Makino, M.; Tatarazako, N., Acute toxicity of parabens and their

576

chlorinated by-products with Daphnia magna and Vibrio fischeri bioassays. J Appl

577

Toxicol 2009, 29, (3), 242-247.

578

50. Tian, S.; Sun, H.; Li, Y.; Pan, P.; Li, D.; Hou, T., Development and evaluation of

579

an integrated virtual screening strategy by combining molecular docking and

580

pharmacophore searching based on multiple protein structures. J Chem Inf Model

581

2013, 53, (10), 2743-2756.

582

51. Zhou, S.; Li, Y.; Hou, T., Feasibility of using molecular docking-based virtual

583

screening for searching dual target kinase inhibitors. J Chem Inf Model 2013, 53, (4),

584

982-996.

585

52. Shen, M.; Zhou, S.; Li, Y.; Pan, P.; Zhang, L.; Hou, T., Discovery and

586

optimization of triazine derivatives as ROCK1 inhibitors: molecular docking,

587

molecular dynamics simulations and free energy calculations. Mol Biosyst 2013, 9, (3),

588

361-374.

589

53. Shen, M.; Zhou, S.; Li, Y.; Li, D.; Hou, T., Theoretical study on the interaction of

590

pyrrolopyrimidine derivatives as LIMK2 inhibitors: insight into structure-based

591

inhibitor design. Mol Biosyst 2013, 9, (10), 2435-2446.

27

ACS Paragon Plus Environment

Environmental Science & Technology

592

54. Shen, M.; Yu, H.; Li, Y.; Li, P.; Pan, P.; Zhou, S.; Zhang, L.; Li, S.; Lee, S. M.;

593

Hou, T., Discovery of Rho-kinase inhibitors by docking-based virtual screening. Mol

594

Biosyst 2013, 9, (6), 1511-1521.

595

55. Bohl, C. E.; Gao, W.; Miller, D. D.; Bell, C. E.; Dalton, J. T., Structural basis for

596

antagonism and resistance of bicalutamide in prostate cancer. Proc Natl Acad Sci U S

597

A 2005, 102, (17), 6201-6206.

598

56. Watanabe, Y.; Kojima, H.; Takeuchi, S.; Uramaru, N.; Ohta, S.; Kitamura, S.,

599

Comparative study on transcriptional activity of 17 parabens mediated by estrogen

600

receptor alpha and beta and androgen receptor. Food Chem Toxicol 2013, 57,

601

227-234.

602

57. Jewell, C.; Prusakiewicz, J. J.; Ackermann, C.; Payne, N. A.; Fate, G.; Voorman,

603

R.; Williams, F. M., Hydrolysis of a series of parabens by skin microsomes and

604

cytosol from human and minipigs and in whole skin in short-term culture. Toxicol

605

Appl Pharmacol 2007, 225, (2), 221-228.

606 607 608 609 610 611 612 613

28

ACS Paragon Plus Environment

Page 28 of 37

Page 29 of 37

Environmental Science & Technology

614

615

Table 1. The CAS numbers, structures and Log Kow of the test parabens.

a

Structure

Log Kow a

Compound

CAS

4-hydroxybenzoic acid (HbA)

99-96-7

1.58

methyl paraben (MeP)

99-76-3

1.96

ethyl paraben (EtP)

120-47-8

2.47

propyl paraben (PrP)

94-13-3

3.04

butyl paraben (BuP)

94-26-8

3.57

heptyl paraben (HeP)

1085-12-7

4.83

octyl paraben (OcP)

1219-38-1

5.43

dodecyl paraben (DdP)

2664-60-0

7.40

phenyl paraben (PnP)

17696-62-7

3.21

benzyl paraben (BzP)

94-18-8

3.56

The values of Log Kow were calculated by the U.S. EPA EPI Suite version 4.1.

616 617 618 619 29

ACS Paragon Plus Environment

Environmental Science & Technology

Page 30 of 37

620

Table 2. Parameters (Km, Vmax) of the Lineweaver-Burk plots for the hydrolysis

621

kinetics of parabens. Km (mM)

Vmax (mM/min)

R2

methyl paraben

0.40 ± 0.01

0.19 ± 0.02

0.9976

ethyl paraben

0.33 ± 0.02

0.17 ± 0.01

0.9981

propyl paraben

0.27 ± 0.04

0.13 ± 0.02

0.9958

butyl paraben

0.17 ± 0.04

0.09 ± 0.01

0.9675

phenyl paraben

6.00 ± 0.28

3.43 ± 0.10

0.9939

benzyl paraben

8.23 ± 0.15

4.91 ± 0.08

0.9997

Compounds

622 623 624 625 626 627 628 629 630 631 632 633 634 635 636

30

ACS Paragon Plus Environment

Page 31 of 37

Environmental Science & Technology

637

Figure Legends

638

Figure 1. The antiandrogenic activity of paraben analogues measured by the yeast

639

two-hybrid assay. The antiandrogenic activity was expressed as the values of

640

β-galactosidase activity (U) for positive control substrate of the samples divided by

641

the values of U for the positive control. Values were measured from three independent

642

experiments and were presented as the mean ± standard error. The symbol * stands for

643

the significant difference in comparison with the control (p < 0.05).

644

Figure 2. The linear regression curves for the maximum antiandrogenic activity of

645

parabens vs. the calculated binding energy.

646

Figure 3. The binding mode of parabens to AR-LBD. The protein is shown in cartoon

647

form. Parabens and surrounding residues are represented in sticks form. The dashed

648

lines indicate the formation of hydrogen bonds.

649

Figure 4. The antiandrogenic activity of parabens before and after hydrolysis

650

mediated by porcine liver esterase. The symbol * stands for the significant difference

651

in comparison without the incubation with porcine liver esterase (p < 0.05).

652

Figure 5. Mass spectrometric fragmentation of metabolite formed after the hydrolysis

653

of parent parabens. (A) and (B) are the mass spectra in full-scan and MS2 modes,

654

respectively.

655 656 657

31

ACS Paragon Plus Environment

Environmental Science & Technology

Page 32 of 37

658

Figure 1

659 660 661 662 663 664 665 666 667 668 669 670 671 672 673

32

ACS Paragon Plus Environment

Page 33 of 37

Environmental Science & Technology

674

Figure 2

675 676 677 678 679 680 681 682 683 684 685 686 687 688

33

ACS Paragon Plus Environment

Environmental Science & Technology

Page 34 of 37

689

690 691

Figure 3

692 693 694 695 696 697 698 699 700 701 702 703

34

ACS Paragon Plus Environment

Page 35 of 37

Environmental Science & Technology

704

705

Figure 4

706 707 708 709 710

35

ACS Paragon Plus Environment

Environmental Science & Technology

Page 36 of 37

711

Figure 5

712 713 714 715 716 717 718 719 720 721 722 723 36

ACS Paragon Plus Environment

Page 37 of 37

Environmental Science & Technology

724

Table of Contents Graphic

725

726

37

ACS Paragon Plus Environment