Comparative Toxicokinetics of Organic Micropollutants in

Jun 11, 2013 - Exposure and depuration experiments for Gammarus pulex and Daphnia magna were conducted to quantitatively analyze biotransformation ...
4 downloads 0 Views 1MB Size
Subscriber access provided by Library, Univ of Limerick | Supported by IReL

Article

Comparative toxicokinetics of organic micropollutants in freshwater crustaceans Junho Jeon, Denise Kurth, Roman Ashauer, and Juliane Hollender Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es400833g • Publication Date (Web): 11 Jun 2013 Downloaded from http://pubs.acs.org on June 17, 2013

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

Comparative toxicokinetics of organic micropollutants in

2

freshwater crustaceans

3 4

Junho Jeon†*, Denise Kurth†‡, Roman Ashauer†§, and Juliane Hollender†‡

5 6



Eawag, Swiss Federal Institute of Aquatic Science and Technology, 8600 Dübendorf,

7 8

Switzerland ‡

Institute of Biogeochemistry and Pollutant Dynamics, ETH Zürich, CH-8092, Zürich,

9 10

Switzerland §

Environment Department, University of York, Heslington, York, United Kingdom

11 12 13

* Corresponding author

14 15

Tel: +41 58 765 5042

16

Fax: +41 58 765 5893

17

E-mail: [email protected]

18 19

Number of words: 5,150

20

Number of figures: 3 (3 large, 1,800 words equivalents)

21

Number of tables: 1 (1 large, 600 words equivalents)

22

7,550 words in total

ACS Paragon Plus Environment

Environmental Science & Technology

23

Abstract

24

Exposure and depuration experiments for Gammarus pulex and Daphnia magna were

25

conducted to quantitatively analyze biotransformation products (BTPs) of organic

26

micropollutants (tramadol, irgarol, and terbutryn). Quantification for BTPs without

27

available standards was performed using an estimation method based on physico-

28

chemical properties. Time-series of internal concentrations of micropollutants and BTPs

29

were used to estimate the toxicokinetic rates describing uptake, elimination, and

30

biotransformation processes. Bioaccumulation factors (BAF) for the parents and retention

31

potential factors (RPF), representing the ratio of the internal amount of BTPs to the

32

parent at steady state, were calculated. Nonlinear correlation of excretion rates with

33

hydrophobicity indicates that BTPs with lower hydrophobicity are not always excreted

34

faster than the parent compound. For irgarol, G.pulex showed comparable elimination,

35

but greater uptake and BAF/RPF values than D.magna. Further, G. pulex had a whole set

36

of secondary transformations that D. magna lacked. Tramadol was transformed more and

37

faster than irgarol and there were large differences in toxicokinetic rates for the

38

structurally similar compounds, irgarol and terbutryn. Thus, predictability of

39

toxicokinetics across species and compounds needs to consider biotransformation and

40

may be more challenging than previously thought because we found large differences in

41

closely related species and similar chemical structures.

ACS Paragon Plus Environment

Page 2 of 37

Page 3 of 37

42 43

Environmental Science & Technology

Introduction Xenobiotic compounds, such as pesticides, personal care products, pharmaceuticals, and

44

industrial chemicals are detected in freshwater environments worldwide. Freshwater

45

organisms accumulate contaminants via active or passive uptake routes and may suffer

46

from adverse effects if concentrations in the organisms reach critical levels.1, 2 The extent

47

of bioaccumulation, quantified by the bioaccumulation factor (BAF), is a key index in

48

assessing chemical risks to the environment.3, 4 Organisms eliminate accumulated

49

chemicals to reduce their internal concentrations in a number of ways e.g., excretion via

50

gills, dermal tissue, urine and feces, and biotransformation. Biotransformation is a

51

biological process aimed at modifying xenobiotics to derivatives with generally increased

52

hydrophilicity to promote the excretion from the body and thus reduce bioaccumulation

53

and potential toxic effects. Since the modifications are mediated by enzymes in most

54

cases and different organisms have different kinds and amounts of enzymes,

55

biotransformation pathways and kinetics vary considerably across species5-7 as well as

56

compounds.8, 9 The processes of uptake, internal distribution, biotransformation and

57

elimination of xenobiotic substances in organisms are termed toxicokinetics (TK) and

58

each process can be quantified by rate constants.

59

Cytochrome P450 isoenzymes (P450) are the heme-thiolate superfamily of enzymes and

60

perform extensive oxidative transformations to form hydroxylation, oxidation,

61

deamination, dealkylation, and dehalogenation products.10 These transformations play an

62

important role in regulating the overall bioaccumulation (or body burden) and potential

63

toxicity of a substance.11-14

ACS Paragon Plus Environment

Environmental Science & Technology

64

There have been many studies on TK of pharmaceuticals and biocides in their target

65

organisms (e.g., mammalians for pharmaceuticals,15-18 algae for biocides19) but limited

66

information is available on non-target organisms such as the freshwater crustaceans,

67

Gammarus pulex and Daphnia magna. G.pulex is a key species in flowing freshwater

68

bodies, where it acts as a shredder of detritus and is consumed by fish. D.magna feeds on

69

algae in lakes and is an important food source for fish and is thus a key species in lake

70

ecosystems. Information on qualitative and quantitative analysis of biotransformation

71

products (BTPs or metabolites) in those organisms is also lacking because conventional

72

analytical techniques such as liquid chromatography coupled with mass spectrometry

73

(LC/MS) or radioactivity detection (LC/RAD), were not sensitive nor selective enough to

74

identify BTPs present at trace concentrations in such small organisms. Recently, liquid

75

chromatography high resolution mass spectrometry (LC-HRMS, e.g., Time-of-flight and

76

Orbitrap) has become a promising technique for metabolite profiling studies due to its

77

sensitivity, selectivity and speed of analysis.20-23 The major obstacle in metabolite

78

profiling with mass spectrometric methods is the limited availability of reference

79

standards for metabolites, leading to the infeasibility of absolute quantification. Semi-

80

quantitative evaluation has been established by comparing peak intensity or area of

81

metabolite to those of parent compound, assuming that the ionization efficiency (IE)

82

during the ion formation in electrospray ionization (ESI) and the response of metabolite

83

in the mass detector corresponds to that of the parent compound.24, 25A number of studies

84

clearly showed the strong relationship between IE and chemical properties such as

85

molecular volume and basicity (pKa),26, 27 which can be used to predict the response in

86

the mass detector. Previously, we identified several BTPs of micropollutants and

ACS Paragon Plus Environment

Page 4 of 37

Page 5 of 37

Environmental Science & Technology

87

elucidated their structures using LC-HRMS/MS (LTQ-Orbitrap).28 It was observed that

88

structure similarity resulted in analogous biotransformation pathways for two pairs of

89

target compounds, i.e. irgarol and terbutryn as well as tramadol and venlafaxine.

90

However, it is still unknown whether TK for the analogous compounds are also similar,

91

hampering our ability to predict biotransformation rates as well as concentrations of

92

BTPs.

93

The present study aimed at better understanding the role of biotransformation in

94

toxicokinetics. Specific objectives were quantitative comparisons in toxicokinetic

95

processes and rates (i) between two analogous compounds having an identical molecular

96

back-bone structure (irgarol and terbutryn), (ii) two dissimilar compounds with different

97

functionalities (irgarol and tramadol) and (iii) two different species (G. pulex and D.

98

magna). In addition, the role of P450 activity in the TK of D. magna was investigated

99

through experiments with and without the presence of a P450 inhibitor. The implications

100

for species sensitivity differences and predictability of TK in chemical risk assessment

101

are also discussed.

ACS Paragon Plus Environment

Environmental Science & Technology

102

Materials and Methods

103 104

Chemicals and Test organisms

105

Detailed information on chemical standards and the maintenance and culture of test

106

organisms, i.e. Gammarus pulex and Daphnia magna are provided in the supporting

107

information (SI1, SI2 and Table S1). Physico-chemical properties of the target

108

compounds irgarol, terbutryn, and tramadol, and their BTPs identified in our previous

109

study28 are described in Table S2. The BTPs were named with a scheme described as

110

follows; e.g., MIR327A, ‘M’ stands for metabolite, the following two alphabet letters (IR)

111

represent the first two letters of the name of parent compound (irgarol), the following

112

three-figures (327) represent [M+H]+ of BTP, A or B at the end is a differentiator for

113

BTPs having the same exact mass.

114 115 116

Kinetic experiments The TK study closely followed our previous study28 and Ashauer et al.9 Briefly, G.

117

pulex were exposed to the test compounds for 24 h and then transferred into clean media

118

for the 4 day-depuration phase. D. magna were exposed for 20 hours to the test

119

compounds and then transferred to clean media for a 24 h depuration phase. Samples of

120

whole G. pulex and D. magna for tissue residue analysis were taken frequently

121

throughout the experiment. The specific conditions of the experiments are described in

122

SI3. On each sampling date the sampled organisms were sieved, blotted dry on tissue,

123

placed in a 2 ml microcentrifuge tube, weighed and shock frozen in liquid nitrogen. The

ACS Paragon Plus Environment

Page 6 of 37

Page 7 of 37

Environmental Science & Technology

124

organism and exposure solution samples were kept at -80 °C until extraction and

125

chemical analysis. The chemical extraction procedure with Fastprep is described in SI4.

126

In order to examine biotransformation catalyzed by P450, a group of D.magna was pre-

127

exposed to 500 µg/L of the P450 inhibitor piperonyl-butoxide (PBO), 3 hours before

128

adding irgarol to the solution. The remaining procedure was identical to those described

129

above, without the depuration phase. Two controls, one without target chemicals

130

(chemical-negative) and one without test organisms (organism-negative), were included

131

in the experiments. Additionally, a food positive-organism negative control was

132

conducted for G. pulex to quantify the formation of transformation products by microbes

133

on the food source (inoculated horse chestnut leaves).

134 135 136

LC- High resolution mass analysis Organism extracts from the kinetic experiments were cleaned up with online solid phase

137

extraction and then analyzed with liquid chromatography electrospray ionization -

138

tandem mass spectrometry (Q Exactive Orbitrap LC-ESI-MS/MS, Thermo Fisher

139

Scientific Inc.). Aside from the mass detector, the analytical procedure is almost identical

140

to the previous study to identify the BTPs.28 Detailed descriptions are available in SI5.

141 142 143

Quantification of target molecules To quantify compounds with a reference standard available (i.e. the parent compounds

144

and the BTP (2-methylthio-4-tert-butyl amino-s-triazine), calibration curves were

145

established on the basis of the area response ratio of the reference standard to the

146

respective deuterated internal standard (IS). For the remaining BTPs where reference

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 37

147

standards were not commercially available, we used an adjusted response area ratio. The

148

ionization efficiency (IE) ratio of BTP to the molar concentration of the parent was

149

estimated using physico-chemical properties as described elsewhere26, 27 and the

150

concentration of the BTP was calculated based on the calibration curve established for the

151

parent. The change in the IE caused by the chromatographic retention time difference of

152

parent and BTP was also considered. The concentration was considered as 0 if the signal

153

to noise ratio (S/N) of the target peak in chromatogram was below the limit of detection

154

(LOD). When the S/N was between LOD and the limit of quantification (LOQ), the

155

concentrations were set to (LOD+LOQ)/2. Detailed information on the estimation of BTP

156

concentration and LOD/LOQ are provided in the SI6 and Table S3, respectively.

157 158

Evaluation of biotransformation kinetics

159

A first-order model was applied for the estimation of uptake, biotransformation and

160

elimination rate constants as described below. Ashauer et al.9 reported that the percentage

161

of dietary uptake in G. pulex for organic compounds (log KOW range from 1.1 to 5.18)

162

was generally less than 1 % of total uptake. We included one combined rate constant to

163

estimate the contribution of uptake via food and water in the BAF estimation for G. pulex.

164

As the contribution of uptake via food is likely very small the resulting BAF will be close

165

to the bioconcentration factor (BCF). The TK model is described by:

166

dCorg _ p (t ) 167 168

dt

= C water (t ) × ku − Corg _ p (t ) × k e (1)

And at steady-state:

ACS Paragon Plus Environment

Page 9 of 37

169

Environmental Science & Technology

BAF (or BCF) =

C org_p (t ) C water (t )

=

ku ku = k e k ex + k bt

(2)

170 171

Where Corg_p(t) is the time course of the concentration of the parent compound in the

172

organism [nmol kgwet weight−1], Cwater(t) is the time course of the concentration of the

173

parent compound in water [nmol L−1], ku is the uptake rate constant [L kgwet weight−1 d−1],

174

ke is the overall elimination rate constant, corresponding to the sum of the direct excretion

175

(kex) and the total biotransformation rate constant (kbt) of the parent compound [d−1].

176

According to the proposed pathways in our previous study, the biotransformation rate

177

depends on the internal concentration of either the parent compound or the respective

178

BTP formed upstream in the biotransformation pathway (Figure S1). For the rate

179

estimation, BTPs are classified into two groups: primary BTPs formed directly from

180

parent compounds and secondary BTPs formed from a preceding BTP. The TK model

181

was refined accordingly as follows, assuming that biotransformation can be modeled as a

182

first-order process.

183 184

For the primary BTPs:

dCorg_1stBTn(t) 185 186

dt

= Corg_p(t) × k1stBTn − Corg_1stBTn(t) × k2 ndBTn − Corg_1stBTn(t) × kex_1stBTn (3)

187 188

For the secondary BTPs

ACS Paragon Plus Environment

Environmental Science & Technology

dC org_ 2 ndBTn (t) 189 190

dt

= Corg_ 1stBTn (t) × k 2 ndBTn − Corg_ 2 ndBTn (t) × k ex_ 2 ndBTn

(4) 191 192

where Corg_1stBTn(t) and Corg_2ndBTn(t) are the time courses of the concentration [nmol kgwet

193

weight

194

k2ndBTn are the biotransformation rate constants [L kgwet weight−1 d−1] of the primary and

195

secondary BTPs, respectively, and kex_1stBTn and kex_2ndBTn are the elimination rate

196

constants [d−1] of the primary and secondary BTPs, respectively.

−1

] of the primary and secondary BTPs in the organism, respectively, k1stBTn and

197 198

Model calibration

199

The rate constants were estimated for each BTP by first approximating the parameters of

200

the primary BTPs and then simultaneously fitting all toxicokinetic parameters (model

201

calibration). The model calibration minimized the differences between simulated internal

202

concentrations for all parent compounds and BTPs and measured concentrations of all

203

parent compounds and BTPs at all sampling dates. The model calibration assumed

204

normally distributed, independent errors in the internal concentration data and thus

205

employed a least squares minimization using the Levenberg-Marquard algorithm.

206

Simulations and parameter estimation were carried out in ModelMaker (v4.0, Cherwell

207

Scientific Ltd., Oxford, UK). Further information on model calibration is available in SI7.

208 209 210 211

The Retention Potential Factor We propose the retention potential factor (RPF) as a measure of the relative significance of BTP in bioaccumulation. RPF is the ratio of the concentration of a BTP to the

ACS Paragon Plus Environment

Page 10 of 37

Page 11 of 37

Environmental Science & Technology

212

concentration of the parent in the organism at steady state, corresponding to the ratio of

213

the biotransformation rate constant to the elimination rate constant as described below.

214 215

For the primary BTPs,

216

RPF1stBTn =

Corg_1stBTn (t ) Corg_p (t )

k1stBTn kex_1stBTn + k2ndBTn

=

(6)

217 218

For the secondary BTPs,

RPF2ndBTn = 219

=

C org_2ndBTn (t ) C org_p (t )

=

C org_1stBTn (t ) C org_p (t )

×

C org_2ndBTn (t ) C org_1stBTn (t )

k1stBTn k × 2ndBTn k ex_1stBTn + k 2ndBTn kex_2ndBTn

220

(7)

221 222

Thus, the relative internal concentration of BTPs to their parent compound at steady state

223

can be calculated on the basis of RPF values.

224

ACS Paragon Plus Environment

Environmental Science & Technology

225

Results and discussions

226 227

Validation of BTP (semi)quantification method

228

BTPs without reference standards were quantified using estimated relative ionization

229

efficiencies of BTPs to the parent. To verify the quantification method, the estimated

230

concentration of the BTP MIR214 (MTE214) was compared with the concentration

231

measured using the available reference standard. Overall, on average, the estimated

232

concentration was 81.2 % of the measured concentration (standard deviation (sd), ±23.5).

233

This compares favorably with another common (semi)quantification method, which

234

simply applies the area ratio of BTP to IS of the parent compound directly to the

235

calibration curve of the parent without any conversion factor.24 Applying this method to

236

MIR214 significantly underestimates the concentration (average ± sd, 26.4 % ± 6.6). This

237

comparison clearly shows that our quantification method can improve the accuracy of the

238

quantification at least for structurally similar BTPs. The detailed information on the

239

validation of quantification method is available in SI8.

240 241

Kinetic rate estimation

242

The bioaccumulation of the parent compounds in G. pulex and D. magna reached a

243

steady state within 24h and 20h exposure, respectively (Figure 1). At the end of the 96h

244

depuration period, the parent compounds were almost completely eliminated from G.

245

pulex (less than 1 % of the highest internal concentration remained), while MIR270A,

246

MIR327A, MIR327B, and MTR266 were still present in the organism. In contrast, the

247

parent compounds in D. magna and all BTPs except MTE214 were almost completely

ACS Paragon Plus Environment

Page 12 of 37

Page 13 of 37

Environmental Science & Technology

248

eliminated after 24h depuration. MTE214 concentrations at the end of the depuration

249

phase in D. magna were still equivalent to about 5 % of the highest internal concentration

250

during the exposure phase.

251

Since all kinetic parameters for a parent compound and its BTPs were simultaneously

252

estimated in the complete TK model, biotransformation and elimination rate constants are

253

mutually related to each other in the model (co-variation). This co-variation is reflected in

254

the parameter standard errors (SEs), which are consequently larger than SEs derived from

255

estimations for single reaction rates (Table S4). Our simultaneous fit integrates the error

256

and uncertainty for all reactions.

257 258

Comparison between species with irgarol

259

Toxicokinetics for irgarol and 4 BTPs were compared between G. pulex and D. magna

260

(Figure 2(a) and (b)) as MIR214, MIR270A, MIR270B, and MIR327B were successfully

261

quantified in both species. The comparison of TK in both species indicates that N-

262

dealkylation and hydroxylation are major transformation mechanisms for irgarol in D.

263

magna, while the products formed by hydroxylation and glutathione conjugation resulting

264

in the cysteine conjugate formation are dominant in G. pulex.

265

G. pulex exhibited a higher uptake (ku) and excretion rate (kex) of irgarol compared with

266

D. magna, but the overall elimination rate, ke (sum of excretion and biotransformation

267

rates, kex+ k1stBT1+ k1stBT2+ k1stBT3+ k1stBT4) was comparable in both species (Table S4),

268

resulting in a higher BAF in G. pulex (Table 1). The extremely low kex value (0.0002) and

269

relatively high discrepancy between ke and kex (ke- kex=0.35) in D. magna indicate that

270

most of the accumulated parent molecules are transformed and thus only few are excreted

ACS Paragon Plus Environment

Environmental Science & Technology

271

in the unchanged form (Figure 2(b)). Contrarily, the relatively high kex and low sum of

272

biotransformation rates for G. pulex indicate that the significance of biotransformation for

273

the elimination of irgarol in G. pulex is less than that in D. magna. For the primary BTPs,

274

MIR270A, MIR270B, and MIR327B, which all produce the secondary BTPs in G. pulex,

275

the overall elimination is regulated by both the transformation to the secondary BTP and

276

the direct excretion to surrounding medium. The direct excretion of MIR270A and

277

MIR270B was the dominant elimination process in G. pulex, making up about 67 % and

278

98 %, respectively, of the overall elimination, whereas almost all MIR327B was

279

transformed to MIR327A (Figure 2(a)). On the other hand, no secondary BTPs were

280

identified in D. magna, due either to a relatively low transformation or high excretion rate.

281

We simulated a maximum internal concentration of a secondary BTP, MIR286 in D.

282

magna, by assuming that k2ndBT2 and k2ndBT3 were identical to k1stBT3 and k1stBT2 (0.067 and

283

0.259 h-1), respectively, and the excretion rate was identical to kex_2ndBT2 estimated for G.

284

pulex (0.346 h-1) (see Table 1). The simulated concentration of approximately 350

285

nmol/kg is above the LOQ (about 6 nmol/kg) and indicates that the non-detection of the

286

BTP in D. magna is mainly due to faster excretion than that in G. pulex. If this hypothesis

287

is true, the calculated excretion rate of MIR286 in D. magna necessary to reduce internal

288

concentrations below the LOQ should be greater than 4 h-1, which is about 12 times

289

greater than in G. pulex. This faster excretion is consistent with the comparison of the

290

elimination kinetics between the two species: the total sum of elimination (or excretion)

291

rate for the primary BTPs in D. magna is about 6.6 (or 7) times greater than that in G.

292

pulex.

ACS Paragon Plus Environment

Page 14 of 37

Page 15 of 37

Environmental Science & Technology

293

Biotransformation rates (k1stBT1, k1stBT2, and k1stBT3) for oxidation and hydroxylation in D.

294

magna are greater than those in G. pulex; about 4 times for MIR214, 18 times for

295

MIR270A and 4 times greater for MIR270B. Consequently, the elimination rates,

296

kex_1stBT1, kex_1stBT2, and kex_1stBT3 in D. magna are about 3 times, 100 times and 2 times

297

greater than the overall elimination rate in G. pulex, respectively. Note that these

298

differences in rate constants are larger than what could be explained by the temperature

299

difference in the D.magna and G.pulex studies. For MIR327B we estimated a greater

300

biotransformation and a lower overall elimination rate in G. pulex, resulting in a RPF

301

about 13 times higher than for D. magna (Table 1). In addition to MIR327B, the RPF of

302

MIR270A in G. pulex is greater than that in D. magna, whereas we found a greater RPF

303

for MIR214 and MIR270B in D. magna. Overall, the total sum of RPFs representing the

304

ratio of the total amount of BTPs to the parent at steady state was 0.56 and 0.17 for G.

305

pulex and D. magna, respectively. The greatest RPF in G. pulex was estimated for MIR

306

327B, comparable to MIR270B, followed by >>MIR214 ≥ MIR270B whereas the RPFs

307

in D. magna were all similar, although declining in the order MIR214≥ MIR270B ≥

308

MIR270A>MIR327B (Table 1).

309

The greater RPFs for G. pulex are attributable to the lower elimination rates of the BTPs.

310

However, this does not mean that the overall toxicity of irgarol in G. pulex is less than in

311

D. magna since the absolute amount of irgarol accumulated was higher in G. pulex and

312

the relative toxicity potential of each of the BTPs are not known. In this study, no adverse

313

effects were observed for either organism during the exposure at 200 µg/L of irgarol, so

314

that our insight on the overall toxicity and in particular the contribution of BTPs is

315

limited. Nevertheless, the different amounts and compositions of BTPs in different

ACS Paragon Plus Environment

Environmental Science & Technology

316

species as indicated by the RPF values are likely related to differences in species

317

sensitivities to toxicants, which should be further explored.

318 319

Nonlinearity of excretion rate over hydrophobicity

320

In general, the ability to biotransform xenobiotics contributes to organisms’ defenses

321

against toxicants by generating BTPs that lack the toxicophore, or by generating BTPs

322

that are less hydrophobic and thus excreted faster. The faster excretion of BTPs leads to

323

lower internal concentrations and thus to a reduction of toxicity which is a function of the

324

number of accumulated molecules in target sites. We observed that biotransformation in

325

G. pulex and D. magna generally reduces hydrophobicity (see log Dow values, Table 1),

326

but that the excretion of some BTPs is not faster than that of the parents (Table 1 and

327

Table S4). The common assumption that reduced hydrophobicity results in faster

328

excretion of BTPs in aquatic organisms has rarely been tested empirically, because

329

studies with in depth TK characterization such as the one presented here are challenging.

330

To extensively establish the correlation of the excretion rate with hydrophobicity, we

331

combined datasets from both the present study and Ashauer et al. (only for G. pulex).9

332

Interestingly, the hydrophobicity correlates non-linearly with the excretion rates for both

333

organisms (Figure 3). The excretion rate reaches the maximum between log Dow = 1 - 2

334

and decreased with increasing/decreasing hydrophobicity of molecules which is

335

consistent with an observation in fish29. As the hydrophobicity of tramadol BTPs (except

336

MTR280B) in G. pulex, is reduced, the excretion rates decreased, resulting in the

337

increased elimination half-lives (Table 1). The low excretion for hydrophilic (ionized)

338

molecules can be explained by the lower permeability across cell membranes.30 In the

ACS Paragon Plus Environment

Page 16 of 37

Page 17 of 37

Environmental Science & Technology

339

membrane permeation concept based on fugacity theory,30-32 more hydrophobic

340

compounds (e.g. a parent) can more easily gain access into a cell membrane by

341

hydrophobic interaction. Once the compound permeates the membrane and is

342

biotransformed in the cell, the BTPs have less potential to overcome the membrane due to

343

the reduced hydrophobicity leading to the retarded excretion. Hendriks et al.33 also

344

reported declining elimination rates with decreasing log Kow below 3, while the increase

345

in elimination rates with decreasing log Kow from 7 to 3 was also shown in other

346

studies.33-36 However the authors could not fully confirm declining elimination rates for

347

substances with lower log Kow value due to the limited number of datasets on hydrophilic

348

substances. The present study and Ashauer et al.9 support the nonlinear correlation of

349

excretion rate with hydrophobicity by providing datasets covering hydrophilic

350

compounds and their BTPs in the range between -1 and 5.2 of log Dow. This observation

351

implies that biotransformation of polar compounds does not always induce faster

352

excretion, but can result in increased BTP retention compared with the parent where the

353

retention potential factor (RPF) is greater than 1 (Table 1).

354 355

Comparison between compounds: Irgarol vs. Tramadol

356

Here, the toxicokinetics in G. pulex are compared for substances with differing physico-

357

chemical properties. The BAF for tramadol is much smaller than irgarol, resulting mainly

358

from the significantly lower ku for tramadol due to the slower uptake of cations. In

359

addition, tramadol transformations to primary BTPs (i.e. MTR250, MTR280A, and

360

MTR280B) comprised about 50 % of the total elimination of the parent, while only 11 %

361

of irgarol was transformed. This higher biotransformation rate for tramadol apparently

ACS Paragon Plus Environment

Environmental Science & Technology

362

contributed to the lower BAF value. In general, the bioaccumulation potential of organic

363

xenobiotic compounds in aquatic organisms is predicted by quantitative structure-activity

364

relationships (QSARs) considering only specific chemical properties such as

365

hydrophobicity.33, 37 However, since the extent of biotransformation is independent of the

366

molecule’s hydrophobicity, ignoring biotransformation seriously limits the predictive

367

power of such models. Thus, it is desirable to develop novel relationships for QSARs

368

which include the extent of biotransformation, for example, based on the presence of

369

certain functional groups in the molecular structure of xenobiotics.

370

Hydroxylation was the most rapid biotransformation mechanism for irgarol, whereas O-

371

demethylation and N-oxidation were key mechanisms for tramadol in G. pulex. The

372

higher sum of RPFs for tramadol BTPs than for irgarol indicates that a relatively large

373

proportion of tramadol is transformed and the BTPs are retained longer in the body than

374

the parent (Table 1). Tramadol is designed to be easily oxidized by P450 in humans to

375

form the active metabolite, O-desmethyltramadol (MTR250) which has a 200 times

376

greater affinity to the target site than the parent.38 Obviously, this fast biotransformation

377

kinetics of tramadol is also similar in freshwater organisms. Except for the glutathione

378

conjugation products, all BTPs identified for both compounds were formed through

379

oxidative reactions mainly mediated by P450 isoenzymes. Hence, the different

380

transformation kinetics for oxidation products can be attributed to the dissimilarity either

381

in physico-chemical properties affecting the affinity of the molecule to the activation site

382

of enzymes or in the amount of activated enzyme affecting the probability for the binding.

383 384

Comparison between compounds: Irgarol vs. Terbutryn

ACS Paragon Plus Environment

Page 18 of 37

Page 19 of 37

Environmental Science & Technology

385

Irgarol and terbutryn are both triazine compounds and used as biocides. Both have a

386

similar hydrophobicity, which is reflected in their similar BAF values. The difference in

387

the structure is a branch moiety bound to an amine: the cyclopropyl in irgarol and the

388

ethyl in terbutryn. Interestingly, this different structural moiety led to a dissimilarity in

389

the biotransformation kinetics. N-dealkylation resulted in an identical BTP (MIR214 and

390

MTE214), but different leaving groups led to a higher kinetic rate in terbutryn. This

391

means that the cleavage of the amine-ethyl bond is faster than that of amine-cyclopropyl.

392

Supposing the enzymes involved in the reaction and the quantity of the enzyme are

393

identical, the discrepancy in the formation rate is likely due to different steric hindrance

394

of the moiety. The ethyl group is smaller and more flexible and thus activated enzymes

395

can easily access the amine-ethyl bond. In contrast, the hydroxylation (k1stBT3) at

396

cyclopropyl and ethyl and the glutathione conjugation rate (k1stBT4) for both compounds

397

were almost the same.

398 399

Significance of P450

400

Internal concentrations of irgarol and terbutryn in D. magna concurrently exposed with

401

the P450 inhibitor PBO were greater than those in PBO-free D. magna. As shown in

402

Figure S2, the time courses of the concentration of the parents and their oxidation

403

products (MIR270A, MIR270B, MTE258A and MTE258B) were affected by PBO

404

leading to the higher bioconcentration of the parent at the end of the exposure phase.

405

About 2-fold (irgarol) and 1.5-fold (terbutryn) greater BAF for PBO-exposed D. magna

406

indicates that more than half of the irgarol and one third of the terbutryn accumulated in

407

D. magna were oxidatively biotransformed by P450. No significant changes for MIR214,

ACS Paragon Plus Environment

Environmental Science & Technology

408

MIR327B, and MTE315B were observed. MIR327B and MTE315B are cysteine

409

conjugates formed by removal of glycyl and glutamyl moieties catalyzed by peptidases

410

and gamma-glutamyltransferase after their predecessors, GSH conjugates are formed by

411

glutathione-s-transferases.39, 40 Thus, the unaffected kinetics of cysteine conjugates prove

412

that all enzymes involved in the formation of cysteine conjugates are activated

413

independently even when the activity of P450 is inhibited. In the case of MIR214, we

414

expected that the increase of internal concentration with time would be retarded in the

415

presence of PBO as observed for MTE214, because general dealkylation reactions are

416

also mediated by CYP450s. However, no PBO effect was observed for MIR214. A

417

possible explanation is that the increasing internal concentration is achieved not only by

418

the transformation of irgarol, but also by the aqueous uptake of MIR214 originated from

419

impurity of the irgarol standard. We could not detect MIR214 in the exposure medium at

420

the beginning of exposure, but observed a small peak of MIR214 in the measurement of

421

the irgarol standard, corresponding to less than 1 % of irgarol. Therefore, the estimated

422

kinetic rates for MIR214 should be interpreted with caution because the rate constants

423

were calculated assuming the absence of MIR214 in the medium.

424 425

Comparison with other species and substances

426

Information on the toxicokinetics of irgarol in different species is currently limited, but

427

some knowledge about the mammalian metabolism of terbutryn and other triazines (e.g.,

428

atrazine) is available. Interestingly, major metabolites of triazines identified in previous

429

studies are inconsistent even when identical species were used in assays.41-43 This might

430

be due to the different exposure schemes, sample treatments and analytical methods.

ACS Paragon Plus Environment

Page 20 of 37

Page 21 of 37

Environmental Science & Technology

431

Despite this, dealkylation products (i.e., desethyl-, desisopropyl, and didealkylatrazine)

432

and conjugates (e.g., mercapturate) were the most frequently identified metabolites in

433

rodent and human samples at moderate levels, whereas hydroxylation products were

434

measured at trace levels or undetectable.42-45 This is different to our observations

435

revealing that hydroxylation was the most significant transformation mechanism for

436

irgarol and terbutryn in both crustaceans tested. In addition to triazines, hydroxylation

437

and oxidation products of tramadol were also predominantly formed in G. pulex. These

438

metabolites were considered minor products in humans and other mammal samples,

439

where O- and N-desmethytramadol were identified as the major metabolites.46, 47 Thus, it

440

is likely that the significance of hydroxylation (or oxidation) in detoxification processes

441

in freshwater crustaceans is greater than in mammals. This characteristic might have

442

evolved because it may be beneficial for responding to naturally occurring chemical

443

stressors such as algal toxins. For example, a number of studies reported that the

444

formation of oxidation products of microcystins by ozonation48 and ultrasonication49, 50

445

apparently decreased toxicity. On the other hand, Ufelmann et al. 51 demonstrated that

446

demethylated derivatives of cyanobacterial toxins, microcystins and nodularins exhibited

447

higher hepatocyte toxicity than their parent molecules. This supports the hypothesis that

448

hydroxylation or oxidation is more important for freshwater crustaceans than

449

demethylation.

450

Comparative studies on toxicokinetics of organophosphates in both test organisms have

451

been conducted. Kretschmann et al.52 showed that biotranformation rates of diazinon,

452

forming diazoxon and dearyl products, were more than 2 time and 6 times greater,

453

respectively, in D. magna than in G. pulex, which is consistent with our observations.

ACS Paragon Plus Environment

Environmental Science & Technology

454

Meanwhile, Rubach et al.5 compared uptake and elimination rates and bioconcentration

455

factors for chlorpyrifos and reported about 2 times greater uptake rates but slightly

456

smaller elimination rates in D. magna compared to G. pulex. Thus, for chlorpyrifos, the

457

pattern seems to deviate from that observed in our study, but the underlying causes

458

remain unclear because the biotransformation kinetics were not quantified in the study.

459

The resulting BAFs of chlorpyrifos for G. pulex and D. magna were two to three orders

460

of magnitude higher than for any of the compounds tested in this study, mainly due to the

461

greater hydrophobicity of chlorpyrifos.

462 463

Implications for the study of species sensitivity differences

464

Internal exposure is the biologically relevant dose and was heavily modified by

465

biotransformation kinetics in both species. The extent of biotransformation differed

466

according to species and chemicals, leading to the dissimilarity in the profile of internal

467

BTPs. Thus, we should expect large differences in biotransformation even amongst fairly

468

closely related species and extrapolation of toxic effects from one species to others must

469

account for biotransformation differences, for example in risk assessment for chemicals.

470

This is currently challenging as the causes of these biotransformation differences are not

471

well understood. In general, the higher the uptake rate, the higher the sensitivity (e.g.,

472

lower EC50),5 as the internal dose at the target site increases faster. However, Rubach et

473

al.5 observed that BAF was not significantly correlated with arthropod sensitivity

474

(R2=0.054), implying that the ability to perform biotransformation and the susceptibility

475

of target sites to the parent and its BTPs may vary across species. This indicates that a

476

simple quantitative analysis for the parent compound (i.e. BAF) is insufficient to fully

ACS Paragon Plus Environment

Page 22 of 37

Page 23 of 37

Environmental Science & Technology

477

understand the observed toxicity and species sensitivity. Thus, the quantitative estimates

478

for biotransformation and elimination rate of BTPs are essential information for the better

479

mechanistic elucidation of toxicity and species sensitivity.

480

ACS Paragon Plus Environment

Environmental Science & Technology

481

Acknowledgement

482

We thank to Heinz Singer, Matthias Ruff, and Philipp Longree for Q Exactive Orbitrap

483

analysis, to Anita Hintermeister for collection and maintenance of G. pulex, to Jelena

484

Simovic for providing data on responses of molecules in mass spectrometry, and to

485

Emma Schymanski and Michael Stravs for comments on the manuscript. This study was

486

supported by funding from the Swiss National Science foundation (Grant No. 315230-

487

141190).

488 489

Supporting information available

490

Details on chemicals and solutions, maintenance and culture of test organisms, kinetic

491

experiments, procedure of extract sample (Fastprep method), LC-HRMS/MS analysis,

492

(semi)quantification of BTPs, model calibration, and validation of quantification method.

493

This material is available free of charge via the Internet at http://pubs.acs.org.

ACS Paragon Plus Environment

Page 24 of 37

Page 25 of 37

Environmental Science & Technology

Table 1. Bioaccumulation factor (BAF), retention potential factor (RPF), and elimination half-life (t1/2) with standard error (SE). Parent log Dow compounds and at pH 7.9* their BTP Irgarol MIR214 MIR270A MIR270B MIR284 MIR286 MIR327A MIR327B Tramadol MTR236 MTR250 MTR266 MTR280A MTR280B Terbutryn MTE214 MTE258A MTE258B MTE315B

2.94 2.21 1.94 2.07 -1.00 1.08 -0.53 -0.68 1.10 0.03 0.95 1.16 0.00 1.33 2.79 2.21 1.83 1.83 -0.79

BAF(L/kg) or RPF(mol/mol) , ±SE D. magna G. pulex without PBO with PBO 19.20±1.89 12.84±25.28 29.20±40.53 0.03±0.05 0.06±0.08 0.03±0.03 0.27±0.40 0.05±0.37 0.02±0.12 0.02±0.09 0.05±0.14 0.02±0.05 0.01±0.13 0.01±0.07 0.03±0.18 0.18±1.03 0.01±0.06 0.01±0.03 0.55±0.09 0.65±0.82 1.12±0.48 4.58±1.57 0.42±0.38 0.32±0.06 12.17±11.31 17.68±11.81 0.58±0.08 0.33±0.04 0.03±0.19 0.00±0.02 0.02±0.30 0.00±0.01 0.02±0.07 0.01±0.05

t1/2 (h)** , ±SE D. magna G. pulex (without PBO) 1.8±0.4 2.0±1.7 3.9±2.1 1.5±0.7 13.1±9.0 0.1±0.1 0.8±0.6 0.5±0.4 32.6±27.7 2.0±0.4 3.4±2.9 20.1±19.6 3.6±2.7 1.1±0.3 7.6±2.4 4.4±1.6 18.2±1.7 12.6±5.3 2.1±0.7 1.6±1.2 2.5±0.2 0.5±0.4 0.2±0.2 4.9±3.6

* Octanol-water partitioning coefficient at pH 7.9 predicted by MarvinSketch Ver. 5.10, **Calculated by t1/2=ln(2)/ke

ACS Paragon Plus Environment

Environmental Science & Technology

Figure 1. Uptake, biotransformation, elimination kinetics of a) irgarol and b) tramadol in G. pulex and c) irgarol and d) terbutryn in D. magna. The exposure and depuration time for G. pulex are 24h and 96h, respectively, whereas 20h and 24h for D. magna. The kinetics of the low concentrated BTPs are depicted in the panel cut-outs in more detail. a) 2500

Irgarol MIR214 MIR270A MIR270B MIR284 MIR286 MIR327A MIR327B

2000 16000

nmol/kg

14000 12000

1500

1000

nmol/kg

10000 8000

500

6000 0

4000

0

20

40

60

80

100

time (h)

2000 0 0

20

40

60

80

100

120

time (h)

ACS Paragon Plus Environment

120

Page 26 of 37

Page 27 of 37

Environmental Science & Technology

b) 1800 1600 Tramadol MTR236 MTR250 MTR266 MTR280A MTR280B

1400

nmol/kg

1200 1000 800 600 400 200 0 0

20

40

60

80

100

120

time (h)

ACS Paragon Plus Environment

Environmental Science & Technology

Page 28 of 37

c) 500

400

nmol/kg

10000

nmol/kg

8000

Irgarol MIR214 MIR270A MIR270B MIR327B

300

200

6000

100

4000

0 0

10

20

30

time (h)

2000

0 0

10

20

30

40

time (h)

ACS Paragon Plus Environment

40

Page 29 of 37

Environmental Science & Technology

d) 300

250 Terbutryn MTE214 MTE258A MTE258B MTE315B

200

nmol/kg

12000

10000

150

100

nmol/kg

8000 50 6000 0 0

10

20

30

4000

time (h) 2000

0 0

10

20

30

40

time (h)

ACS Paragon Plus Environment

40

Environmental Science & Technology

Figure 2. Comparison in toxicokinetics of irgarol between a) G. pulex and b) D. magna. The number of cubes corresponds to the relative internal concentration of the metabolite to the parent whereas the red/blue indicators represent the magnitude of biotransformation/excretion rate, respectively. a)

ACS Paragon Plus Environment

Page 30 of 37

Page 31 of 37

Environmental Science & Technology

b)

ACS Paragon Plus Environment

Environmental Science & Technology

Page 32 of 37

Figure 3. Relationship between hydrophobicity (log Dow) and excretion rate (log kex). Datasets shown here are collected from results of both the present study and Ashauer et al.9 (only for G.pulex). For logarithm conversion, the excretion rates estimated as zero in Ashauer et al. were considered as 0.0001 to avoid negative rate constants in the model estimation. Non-linear curve fitting was performed by adopting an equation described in Hendriks et al.33 The equation is the following. ݇௘௫ ൌ ௉

ଵ ಴ಹమ ൉ሺ஽೚ೢ

ൈ ିଵሻାଵ

ௐ షೖ ഐ భ ఘಹమೀ ା ಴ಹమ ା ವ೚ೢ



Where w = weight of organism (kg), k = rate exponent, PCH2 = lipid fraction of organism, ρH2O = water layer diffusion resistance (d·kgk ), ρCH2 = lipid layer permeation resistance (d·kg-k), γ=water absorption-excretion coefficient(kgk ·d-1). Parameter values used are w=0.00003 and 0.000001 for G. pulex and D. magna, respectively, k=0.25, PCH2=0.01·w-0.04 and γ = 200 kgk ·d-1 for both organisms. The remaining parameters, ρH2O and ρCH2 are optimized to obtain the best fitting curves. 6

4

G.pulex D.magna Reg. for G.pulex Reg. for D.magna

log kex

2

0

-2

-4

-6 0

2

4

6

log Dow

ACS Paragon Plus Environment

Page 33 of 37

Environmental Science & Technology

References (1)

(2)

(3) (4)

(5)

(6) (7)

(8)

(9)

(10) (11)

(12)

(13) (14) (15)

Meador, J.P.; McCarty, L.S.; Escher, B.I.; Adams, W.J., 10th Anniversary Critical Review: The tissue-residue approach for toxicity assessment: Concepts, issues, application, and recommendations. J. Environ. Monit. 2008, 10, (12), 1486-1498. Escher, B.I.; Ashauer, R.; Dyer, S.; Hermens, J.L.M.; Lee, J.-H.; Leslie, H.A.; Mayer, P.; Meador, J.P.; Warne, M.S.J., Crucial role of mechanisms and modes of toxic action for understanding tissue residue toxicity and internal effect concentrations of organic chemicals. Integrated Environ. Assess. Manag. 2010, 7, (1), 28-49. Van Leeuwen, C.J.; Vermeire, T.G., 2nd ed.;, Risk assessment of chemicals - an introduction. 2nd. edition, Springer: Dordrecht, The Netherlands 2007. Howard, P.H.; Muir, D.C.G., Identifying New Persistent and Bioaccumulative Organics Among Chemicals in Commerce II: Pharmaceuticals. Environ. Sci. Technol. 2011, 45, (16), 6938–6946. Rubach, M.N.; Ashauer, R.; Maund, S.J.; Baird, D.J.; Van den Brink, P.J., Toxicokinetic Variation in 15 Freshwater Arthropod Species Exposed to the Insecticide Chlorpyrifos. Environ. Toxicol. Chem. 2010, 29, (10), 2225-2234. Katagi, T., Bioconcentration, bioaccumulation, and metabolism of pesticides in aquatic organisms. Rev. Environ. Contam. Toxicol. 2010, 204, 1-132. Koenig, S.; Fernandez, P.; Sole, M., Differences in cytochrome P450 enzyme activities between fish and crustacea: Relationship with the bioaccumulation patterns of polychlorobiphenyls (PCBs). Aquat. Toxicol. 2012, 108, 11-17. Arnot, J.A.; Mackay, D.; Parkerton, T.F.; Bonnell, M., A Database of Fish Biotransformation Rates for Organic Chemicals. Environ. Toxicol. Chem. 2008, 27, (11), 2263-2270. Ashauer, R.; Hintermeister, A.; O'Connor, I.; Elumelu, M.; Hollender, J.; Escher, B.I., Significance of Xenobiotic Metabolism for Bioaccumulation Kinetics of Organic Chemicals in Gammarus pulex. Environ. Sci. Technol. 2012, 46, (6), 3498-3508. Snyder, M.J., Cytochrome P450 enzymes in aquatic invertebrates: recent advances and future directions. Aquat. Toxicol. 2000, 48, (4), 529-547. Ankley, G.T.; Collyard, S.A., Influence of Piperonyl Butoxide on the Toxicity of Organophosphate Insecticides to 3 Species of Fresh-Water Benthic Invertebrates. Comp. Biochem. Physiol., C: Pharmacol. Toxicol. Endocrinol. 1995, 110, (2), 149-155. Baldwin, W.S.; Leblanc, G.A., In-Vivo Biotransformation of Testosterone by Phase-I and Phase-Ii Detoxication Enzymes and Their Modulation by 20Hydroxyecdysone in Daphnia-Magna. Aquat. Toxicol. 1994, 29, (1-2), 103-117. Ikenaka, Y.; Eun, H.; Ishizaka, M.; Miyabara, Y., Metabolism of pyrene by aquatic crustacean, Daphnia magna. Aquat. Toxicol. 2006, 80, (2), 158-165. Akkanen, J.; Kukkonen, J.V.K., Biotransformation and bioconcentration of pyrene in Daphnia magna. Aquat. Toxicol. 2003, 64, (1), 53-61. Katsuki, R.; Fujita, T.; Koga, A.; Liu, T.; Nakatsuka, T.; Nakashima, M.; Kumamoto, E., Tramadol, but not its major metabolite (mono-O-demethyl

ACS Paragon Plus Environment

Environmental Science & Technology

(16)

(17)

(18) (19)

(20)

(21)

(22)

(23) (24)

(25)

(26) (27)

(28)

tramadol) depresses compound action potentials in frog sciatic nerves. Br. J. Pharmacol. 2006, 149, (3), 319-327. Claahsen-van der Grinten, H.L.; Verbruggen, I.; van den Berg, P.P.; Sporken, J.M.J.; Kollee, L.A.A., Different pharmacokinetics of tramadol in mothers treated for labour pain and in their neonates. Eur. J. Clin. Pharmacol. 2005, 61, (7), 523529. Wu, W.N.; McKown, L.A.; Liao, S., Metabolism of the analgesic drug ULTRAM (R) (tramadol hydrochloride) in humans: API-MS and MS/MS characterization of metabolites. Xenobiotica 2002, 32, (5), 411-425. Dayer, P.; Desmeules, J.; Collart, L., The pharmacology of tramadol. Drugs 1997, 53, 18-24. Mohr, S.; Berghahn, R.; Mailahn, W.; Schmiediche, R.; Feibicke, M.; Schmidt, R., Toxic and Accumulative Potential of the Antifouling Biocide and TBT Successor Irgarol on Freshwater Macrophytes: A Pond Mesocosm Study. Environ. Sci. Technol. 2009, 43, (17), 6838-6843. Yang, Y.O.; Grubb, M.F.; Luk, C.E.; Humphreys, W.G.; Josephs, J.L., Quantitative estimation of circulating metabolites without synthetic standards by ultra-high-performance liquid chromatography/high resolution accurate mass spectrometry in combination with UV correction. Rapid Commun. Mass Spectrom. 2011, 25, (21), 3245-3251. Lu, W.Y.; Clasquin, M.F.; Melamud, E.; Amador-Noguez, D.; Caudy, A.A.; Rabinowitz, J.D., Metabolomic Analysis via Reversed-Phase Ion-Pairing Liquid Chromatography Coupled to a Stand Alone Orbitrap Mass Spectrometer. Anal. Chem. 2010, 82, (8), 3212-3221. Yang, Y.O.; Luk, C.E.; Grubb, M.; Zhang, H.Y.; Humphreys, W.G.; Josephs, J., A new strategy based on UHPLC-LTQ/Orbitrap for simultaneous metabolite identification and quantitation of in vivo PK samples. Drug Metab. Rev. 2009, 41, 134-135. Makarov, A.; Scigelova, M., Coupling liquid chromatography to Orbitrap mass spectrometry. J. Chromatogr., A 2010, 1217, (25), 3938-3945. Du, F.Y.; Liu, T.; Shen, T.; Zhu, F.P.; Xing, J., Qualitative-(semi)quantitative data acquisition of artemisinin and its metabolites in rat plasma using an LTQ/Orbitrap mass spectrometer. J. Mass Spectrom. 2012, 47, (2), 246-252. Bateman, K.P.; Kellmann, M.; Muenster, H.; Papp, R.; Taylor, L., QuantitativeQualitative Data Acquisition Using a Benchtop Orbitrap Mass Spectrometer. J. Am. Soc. Mass Spectrom. 2009, 20, (8), 1441-1450. Oss, M.; Kruve, A.; Herodes, K.; Leito, I., Electrospray Ionization Efficiency Scale of Organic Compounds. Anal. Chem. 2010, 82, (7), 2865-2872. Leito, I.; Herodes, K.; Huopolainen, M.; Virro, K.; Kunnapas, A.; Kruve, A.; Tanner, R., Towards the electrospray ionization mass spectrometry ionization efficiency scale of organic compounds. Rapid Commun. Mass Spectrom. 2008, 22, (3), 379-384. Jeon, J.; Kurth, D.; Hollender, J., Biotransformation Pathways of Biocides and Pharmaceuticals in Freshwater Crustaceans Based on Structure Elucidation of Metabolites Using High Resolution Mass Spectrometry. Chem. Res. Toxicol. 2013, 26, (3), 313-324.

ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37

Environmental Science & Technology

(29)

(30)

(31) (32)

(33)

(34)

(35)

(36)

(37)

(38) (39)

(40) (41)

(42)

(43)

(44)

Sijm, D.T.H.M.; Vanderlinde, A., Size-Dependent Bioconcentration Kinetics of Hydrophobic Organic-Chemicals in Fish Based on Diffusive Mass-Transfer and Allometric Relationships. Environ. Sci. Technol. 1995, 29, (11), 2769-2777. Gobas, F.A.P.C.; Opperhuizen, A.; Hutzinger, O., Bioconcentration of Hydrophobic Chemicals in Fish - Relationship with Membrane Permeation. Environ. Toxicol. Chem. 1986, 5, (7), 637-646. Gobas, F.A.P.C.; Mackay, D., Dynamics of Hydrophobic Organic-Chemical Bioconcentration in Fish. Environ. Toxicol. Chem. 1987, 6, (7), 495-504. Gobas, F.A.P.C.; Muir, D.C.G.; Mackay, D., Dynamics of Dietary Bioaccumulation and Fecal Elimination of Hydrophobic Organic-Chemicals in Fish. Chemosphere 1988, 17, (5), 943-962. Hendriks, A.J.; van der Linde, A.; Cornelissen, G.; Sijm, D.T.H.M., The power of size. 1. Rate constants and equilibrium ratios for accumulation of organic substances related to octanol-water partition ratio and species weight. Environ. Toxicol. Chem. 2001, 20, (7), 1399-1420. van der Linde, A.; Hendriks, A.J.; Sijm, D.T.H.M., Estimating biotransformation rate constants of organic chemicals from modeled and measured elimination rates. Chemosphere 2001, 44, (3), 423-435. Legierse, K.C.H.M.; Sijm, D.T.H.M.; van Leeuwen, C.J.; Seinen, W.; Hermens, J.L.M., Bioconcentration kinetics of chlorobenzenes and the organophosphorus pesticide chlorthion in the pond snail Lymnaea stagnalis - a comparison with the guppy Poecilia reticulata. Aquat. Toxicol. 1998, 41, (4), 301-323. Tolls, J.; Sijm, D.T.H.M., A Preliminary Evaluation of the Relationship between Bioconcentration and Hydrophobicity for Surfactants. Environ. Toxicol. Chem. 1995, 14, (10), 1675-1685. Arnot, J.A.; Gobas, F.A.P.C., A food web bioaccumulation model for organic chemicals in aquatic ecosystems. Environ. Toxicol. Chem. 2004, 23, (10), 23432355. Grond, S.; Sablotzki, A., Clinical pharmacology of tramadol. Clin. Pharmacokinet. 2004, 43, (13), 879-923. Lamoureux, G.L.; Stafford, L.E.; Shimabukuro, R.H.; Zaylskie, R.G., Atrazine metabolism in sorghum: catabolism of the glutathione conjugate of atrazine. J. Agric. Food Chem. 1973, 21, (6), 1020-30. Dekant, W.; Vamvakas, S., Glutathione-Dependent Bioactivation of Xenobiotics. Xenobiotica 1993, 23, (8), 873-887. McMullin, T.S.; Brzezicki, J.M.; Cranmer, B.K.; Tessari, J.D.; Andersen, M.E., Pharmacokinetic modeling of disposition and time-course studies with [C14]atrazine. J. Toxicol. Environ. Health, Part A 2003, 66, (10), 941-964. Lang, D.; Criegee, D.; Grothusen, A.; Saalfrank, R.W.; Bocker, R.H., In vitro metabolism of atrazine, terbuthylazine, ametryne, and terbutryne in rats, pigs, and humans. Drug Metab. Dispos. 1996, 24, (8), 859-865. Abel, E.L.; Opp, S.M.; Verlinde, C.L.M.J.; Bammler, T.K.; Eaton, D.L., Characterization of atrazine biotransformation by human and murine glutathione S-transferases. Toxicol. Sci. 2004, 80, (2), 230-238. Ross, M.K.; Filipov, N.M., Determination of atrazine and its metabolites in mouse urine and plasma by LC-MS analysis. Anal. Biochem. 2006, 351, (2), 161-173.

ACS Paragon Plus Environment

Environmental Science & Technology

(45)

(46)

(47)

(48)

(49)

(50)

(51)

(52)

Lucas, A.D.; Jones, A.D.; Goodrow, M.H.; Saiz, S.G.; Blewett, C.; Seiber, J.N.; Hammock, B.D., Determination of Atrazine Metabolites in Human Urine Development of a Biomarker of Exposure. Chem. Res. Toxicol. 1993, 6, (1), 107116. Garcia-Quetglas, E.; Azanza, J.R.; Sadaba, B.; Munoz, M.J.; Gil, I.; Campanero, M.A., Pharmacokinetics of tramadol enantiomers and their respective phase I metabolites in relation to CYP2D6 phenotype. Pharmacol. Res. 2007, 55, (2), 122-130. Ardakani, Y.H.; Rouini, M.R., Pharmacokinetics of tramadol and its three main metabolites in healthy male and female volunteers. Biopharm. Drug Dispos. 2007, 28, (9), 527-534. Miao, H.F.; Qin, F.; Tao, G.J.; Tao, W.Y.; Ruan, W.Q., Detoxification and degradation of microcystin-LR and -RR by ozonation. Chemosphere 2010, 79, (4), 355-361. Song, W.; De La Cruz, A.A.; Rein, K.; O'Shea, K.E., Ultrasonically induced degradation of microcystin-LR and -RR: Identification of products, effect of pH, formation and destruction of peroxides. Environ. Sci. Technol. 2006, 40, (12), 3941-3946. Song, W.H.; Teshiba, T.; Rein, K.; O'Shea, K.E., Ultrasonically induced degradation and detoxification of microcystin-LR (cyanobacterial toxin). Environ. Sci. Technol. 2005, 39, (16), 6300-6305. Ufelmann, H.; Kruger, T.; Luckas, B.; Schrenk, D., Human and rat hepatocyte toxicity and protein phosphatase 1 and 2A inhibitory activity of naturally occurring desmethyl-microcystins and nodularins. Toxicology 2012, 293, (1-3), 59-67. Kretschmann, A.; Ashauer, R.; Hollender, J.; Escher, B.I., Toxicokinetic and toxicodynamic model for diazinon toxicity-mechanistic explanation of differences in the sensitivity of Daphnia magna and Gammarus pulex. Environ. Toxicol. Chem. 2012, 31, (9), 2014-22.

ACS Paragon Plus Environment

Page 36 of 37

Page 37 of 37

Environmental Science & Technology

986x528mm (96 x 96 DPI)

ACS Paragon Plus Environment