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