Subscriber access provided by HACETTEPE UNIVERSITESI KUTUPHANESI
Article
Mechanistic Pharmacokinetic Modelling of the Bioamplification of Persistent Lipophilic Organic Pollutants in Humans during Weight Loss Li Li, and Frank Wania Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 24 Apr 2017 Downloaded from http://pubs.acs.org on April 24, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 26
Environmental Science & Technology
1
Mechanistic Pharmacokinetic Modelling of the Bioamplification of Persistent Lipophilic Organic
2
Pollutants in Humans during Weight Loss
3
Li Li* and Frank Wania
4
Affiliations: Department of Physical and Environmental Sciences, University of Toronto Scarborough,
5
1265 Military Trail, Toronto, Ontario, Canada, M1C 1A4
6
Correspondence to: L. Li, Department of Physical and Environmental Sciences, University of Toronto
7
Scarborough, 1265 Military Trail, Toronto, Ontario, Canada M1C 1A4. Telephone: (416) 287-5659. Fax:
8
(416) 287-7279. E-mail:
[email protected].
9
TOC Art
10 11
1
ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 26
12
Abstract
13
Bioamplification means the liberation of persistent lipophilic organic pollutants (PLOPs) into blood
14
from their storage in inert adipose tissue during rapid weight loss. Here, using a modified mechanistic
15
pharmacokinetic model, we investigated how chemical properties and anthropometric parameters
16
interact to influence the bioamplification of various PLOPs in humans. The model succeeds in
17
reproducing literature documented weight loss-induced increments in human blood PLOP
18
concentrations. We simulated the degree of bioamplification, as characterized by the bioamplification
19
factor (BAmF), of hypothetical PLOPs with different combinations of partitioning and biotransformation
20
properties at various rates of lipid loss. We also investigated how BAmF evolves with the duration of
21
weight loss. Results show that bioamplification is expected to occur for any chemical with even
22
moderate lipophilicity (logKOW>2 and logKOA>6) as long as the half-life for metabolic elimination is
23
long relative to the time scale of relative lipid loss (e.g. exceeding 104 hours in the case of lipid loss of 3
24
kg month-1 with an initial lipid mass of 40 kg). While BAmF of a chemical is time-variant, whether
25
bioamplification occurs for a chemical or not is independent of the duration of weight loss. The
26
successful application of such a simple model demonstrates that it is the lipid dynamics that
27
predominantly govern the dynamics of PLOPs rather than vice versa.
2
ACS Paragon Plus Environment
Page 3 of 26
Environmental Science & Technology
28
Introduction
29
Humans are exposed to a variety of persistent lipophilic organic pollutants (PLOPs) of concern;
30
prominent examples are polybrominated diphenyl ethers and polychlorinated biphenyls. Within the
31
human body, PLOPs preferentially reside in adipose tissue such as subcutaneous and visceral fats1-3 and
32
only a tiny fraction partitions into plasma/serum lipids (e.g., cholesterol, triglycerides, and phospholipids)
33
and circulates throughout the body.1,
34
toxicological effects in other organs and tissues such as reproductive organs and brain, and leads to slow
35
elimination by limiting biotransformation by xenobiotic metabolizing enzymes and excretion via
36
egestion.1, 5 However, such storage is not necessarily permanent – when an individual loses weight due
37
to gastric surgery, exercise, a low-calorie-diet or disease progression, PLOPs stored in adipose tissue can
38
be mobilized and released into the circulating fluids along with lipolysis.6 If the rate of mobilization
39
exceeds the rate of elimination, the plasma/serum concentration of a lipophilic organic chemical is
40
anticipated to amplify; this phenomenon, in particular when referring to animal studies, is called
41
“bioamplification”.7 The essence of bioamplification is “solvent depletion”;7 namely, as adipose tissue
42
(which can be viewed as solvent) is depleted, the chemical potential of PLOPs (as solutes) therein rises,
43
creating a driving force for redistribution from adipose to other tissues.8
4
Storage in adipose tissue prevents PLOPs from exerting
44
A weight loss-induced increase in plasma/serum chemical concentration can be problematic as it has
45
been linked with elevated adverse health outcomes such as endocrine disruption.9, 10 Bioamplification
46
has even been hypothesized as one of plausible explanations for earlier epidemiological observations of
47
elevated morbidity and mortality among post-weight-loss individuals.11, 12 Moreover, bioamplification
48
during postpartum weight loss may elevate PLOP concentrations in breast milk, which could result in
49
greater postnatal exposure of breastfed infants.13,
50
cost-effective intervention to combat the global obesity epidemic,15 it is of tremendous interest to
51
establish to what extent the negative consequences of bioamplification could counteract the health
52
benefits of weight loss. The above reasons necessitate building a deeper understanding of the mechanism,
53
influencing factors, and adverse outcomes of bioamplification.
14
As weight loss is now recommended as a
3
ACS Paragon Plus Environment
Environmental Science & Technology
54
Page 4 of 26
A number of lab and biomonitoring studies have documented bioamplification of PLOPs in both
55
animals16-22 and humans.9, 10,
12, 23-28
56
inconsistency. For example, Imbeault et al.10 observed statistically significant increases in human plasma
57
concentrations for all compounds but β-hexachlorocyclohexane (HCH) during the same weight loss
58
period. For the same compound, e.g., β-HCH, the increase can be significant in some studies24, 25 but not
59
in others.9, 10 Lab and biomonitoring studies generally provide “snap-shots” of the dynamics of PLOPs at
60
two specific temporal cross-sections before and after a weight loss event. As such, they provide only
61
limited insight into pharmacokinetic processes occurring during a weight loss event, which hinders us to
62
gain comprehensive understanding of how chemical and anthropometric factors interact in determining
63
human bioamplification. To cope with this issue, process-oriented pharmacokinetic (PK) models
64
incorporating parameters that account for chemical and anthropometric properties, can complement the
65
experimental studies. However, while there indeed exist a few mechanistic PK modelling studies on the
66
dynamics of PLOPs in response to periodical (or seasonal) change in body weight of animals like birds29,
67
30
68
model simulations”.7 Furthermore, bioamplification of PLOPs in humans has not yet been subject to PK
69
modelling.
However, their findings are often plagued with variability and
and fish,31 few of them “were focused explicitly on interpreting bioamplification […] resulting from
70
In this work, we use a one-compartment PK model to explore the factors that influence
71
bioamplification of different chemicals in humans and thus may account for the observed variability.
72
Specifically, we first evaluate the model performance by seeking to reproduce bioamplification
73
observations reported in publications. With confidence in its predictive capability established, we
74
investigate bioamplification using diverse combinations of lipid loss rates, chemical partitioning
75
properties and biotransformation half-lives. Furthermore, we investigate how bioamplification changes
76
with time during weight loss.
77
Methods
78
Model Equations 4
ACS Paragon Plus Environment
Page 5 of 26
79
Environmental Science & Technology
The starting point of our PK model is the human module within the ACC-Human model32 (Equation
80
1), which describes the dynamics of a chemical within the human body using fugacity notation:
81
dm d (VH ⋅ Z H ⋅ f H ) = = EOH ∑ ( DuHi ⋅ f UHi ) + DreH f A dt dt i
− ( DmH + EOH DeH + DreH + DperH + DlH + DchH + DurH ) ⋅ f H
(Eq.1)
82
m is the total amount of a chemical within the human body (in mol), which is the product of an
83
“effective” body volume (VH in m3), a fugacity capacity (ZH in mol Pa-1.m-3) and a fugacity (fH in Pa); t is
84
time (in h); DmH, DuHi, DeH, DreH, DperH, DlH, DchH and DurH are D-values (with units of mol Pa-1.h-1)
85
representing chemical transfer via metabolism (or referred to as biotransformation), food uptake,
86
egestion, respiration, percutaneous excretion, lactation, childbirth, and urinary excretion, respectively;
87
EOH is the efficiency of gastrointestinal absorption of the chemical; and fUHi and fA are fugacity of the
88
chemical in foods and air. These terms are calculated according to equations given in ref.32; they are
89
gender-, age- and body-weight-dependent. To make the model geared toward bioamplification
90
simulation, we make two simplifying assumptions:
91
(I) Chemical uptake via diet and inhalation is negligible throughout the period of short-term weight
92
loss, i.e.,
∑(D
uHi
i
⋅ f UHi ) = 0 and DreH f A = 0 . An individual losing weight often has a reduced
93
intake of lipid-rich food12 which is the dominant route for non-occupational PLOP exposure in the
94
general population.33 Meanwhile, a sensitivity analysis (Figure S1) indicates that omitting chemical
95
uptake has a minor (less than a factor of 2) influence on calculation results, irrespective of chemical
96
partitioning properties, when daily intake of a compound is below a certain level (e.g., 10-4
97
mg-chemical kg-1-body-weight d-1 in the case of the absolute rate of lipid loss of 3 kg month-1).
98
Such an intake level exceeds the realistic daily intakes of most PLOPs even for populations that are
99
not dieting (e.g., ref.34, 35);
100
101
(II) Only the fraction of chemical in blood is available for metabolism, which occurs in the liver alone:
ln 2 DmH = HLB′
⋅ Vblood ⋅ Z blood
(Eq. 2)
5
ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 26
102
whereby HLB′ (in h) is the total-blood-volume-adjusted biotransformation half-life (it is an
103
“intrinsic” half-life36 and independent of weight change), Vblood the human blood volume which is
104
height- and body-weight-related,37 and Zblood is the fugacity capacity of blood. In comparison with
105
the original ACC-Human where DmH is related to the whole “effective” body volume (VH) and a
106
whole body biotransformation half-life (HLB), this approach is more appropriate for
107
bioamplification simulation because the smaller change in Vblood than VH during weight loss ensures
108
a non-dramatic change in DmH, i.e., relatively stable hepatic clearance.
109 110
111
Given f H =
m , the assumptions simplify equation 1 to: VH ⋅ Z H
dm m = − ( DmH + EOH DeH + DreH + DperH + DlH + DchH + DurH ) ⋅ dt VH ⋅ Z H
(Eq.3)
Thus, the lipid-normalized chemical concentration C (in molPLOP kg-1lipid) in the human body can be
112
calculated by dividing the total chemical amount by the total lipid mass (Mlipid in kg) (Equation 4):
113
dC d m = dt dt M lipid
D EOH DeH + DreH + DperH + DlH + DchH + DurH dM lipid + = − mH + ⋅C VH ⋅ Z H M lipid ⋅ dt VH ⋅ Z H
(Eq.4)
114
Here, Mlipid includes both the lipids in plasma/serum and the whole body lipid pool (lipids in adipose
115
tissue, liver, muscle, etc.). Since the lipid-normalized chemical concentration C is more often expressed
116
on a mass basis (i.e., ngPLOP kg-1lipid), we do the same conversion throughout this paper.
117 118
The model is coded using the Visual Basic for Applications (VBA) programming language.
Reproducing Observed Bioamplification Reported in the Literature
119
In order to investigate whether Equation 4 is capable of reliably simulating realistic bioamplification
120
situations, we sought to reproduce increments in lipid-normalized plasma/serum concentrations of
121
PLOPs observed in study populations.9, 10, 12, 24-26 The observational data set comprises 91 measurements
122
involving 20 compounds in 9 intentional weight loss events (see details in Supporting Information;
123
Table S1). Here, event refers to the period between two subsequent measurements of a compound of
124
interest. For example, multiple samplings after different durations (e.g., 3, 6 and 9 months after 6
ACS Paragon Plus Environment
Page 7 of 26
Environmental Science & Technology
125
observation starts25) within the same observational cohort were classified as different events;
126
simultaneous observations on male and female subgroups10 were treated as two events.
127
The collected observational data are on the population level (see Supporting Information; Table S1)
128
and do not contain information on the inter-individual variations in anthropometric parameters and
129
chemical concentrations. Here, we used a Monte Carlo simulation to construct an individual-specific
130
dataset with inter-individual random variations in the reported (i) gender, (ii) age, (iii) body weight (both
131
before and post the event), and (iv) the initial lipid-normalized blood concentration of target compounds
132
at the beginning of the event, as information on their variations is available in almost all 91
133
measurements. Combination of variations in the four variables leads to variations in all terms on the
134
right hand of Eq. 4 because the latter are gender-, age- and body-weight-dependent. For each weight loss
135
event, we ran 1000 Monte Carlo iterations, each of which represents a hypothetical individual. For a
136
given hypothetical individual, we randomly assigned its gender according to the published fraction of
137
male subjects (male %; see Supporting Information; Table S1) in the observation population, and
138
sampled its age, body weight and lipid-normalized chemical concentrations based on the respective
139
population-level statistics given in the literature. Given that the population-level statistics are presented
140
in different forms in the literature, we assumed that a variable was either normally distributed if mean
141
and standard deviation were reported, or triangularly distributed if median (or mean) and range (or
142
percentile interval) were reported. Furthermore, we made the following pharmacokinetics-related
143
assumptions:
144
(I) Lipid-normalized PLOP concentrations in plasma/serum are the same as in adipose tissue. This
145
assumption is reasonable for PLOPs because the lipid content is the predominant determinant
146
controlling the partitioning of a lipophilic compound between different tissues,38,
147
residence time within the human body is sufficiently long for a persistent compound to reach
148
equilibrium between the lipids in different tissues. Previous biomonitoring studies show
149
mean/median lipid-normalized PLOP concentrations in the two matrices to be within a factor of
150
4,40-42 which is much smaller than the variability range of the lipid-normalized PLOP concentrations
7
ACS Paragon Plus Environment
39
and the
Environmental Science & Technology
Page 8 of 26
151
themselves. Meanwhile, we do not discriminate between different types of adipose tissues although
152
we recognize that weight loss can affect PLOPs behavior differently in subcutaneous and visceral
153
fats;12
154
(II) Given that either plasma or serum measurements are available in the literature, we do not
155
discriminate whether a concentration refers to that in plasma or serum. This assumption, which has
156
also been adopted in other PK modeling studies such as Verner et al.,43 is justifiable because earlier
157
analytical studies have demonstrated the equivalence of the two concentrations for various PLOPs;44,
158
45
159
(III) Since the total lipid mass Mlipid before and after weight loss is not available in most measurements,
160
we calculate Mlipid from the total body fat mass Mfat by assuming that lipids constitute 83% of total
161
body fat.46 Unless reported in the literature, Mfat is calculated as the product of body weight and
162
body fat percentage (BF%). The BF% is further calculated using an empirical regression47 that
163
correlates BF% with the body mass index (BMI, defined as body weight divided by squared height
164
in kg m-2), age, and sex. Height is age-related (built-in in ACC-Human) and is assumed to be
165
constant during weight loss;
166
(IV) During the weight loss, an individual’s body weight changes linearly.
167
At each Monte Carlo iteration, we simulated the increment in lipid-normalized PLOP concentration
168
(∆Cmodelled) in response to weight loss. The modelled results of the 1000 hypothetical obese individuals
169
were then statistically summarized and compared with the published observed increments.
170
For PLOPs used for the model evaluation, partitioning coefficients at 25 °C were taken from the
171
literature48-51 or calculated using the EPI SuiteTM software. Internal energies (∆UOA, ∆UAW, ∆UOW) for
172
adjusting the cited partitioning coefficients to normal body temperature (37 °C) were calculated based
173
on the method in ref.52 Derivation of HL′ from HLB is detailed in Text S1 of Supporting Information. B
174
The HLB estimates of PLOPs other than β-HCH were derived using the quantitative structure-activity
175
relationship (QSAR) in ref.53 The HLB estimate of β-HCH was calculated from the observed total
176
elimination half-life (HLT)54 using the “1-CoPK model” method in ref.,53 because the QSAR fails to 8
ACS Paragon Plus Environment
Page 9 of 26
Environmental Science & Technology
177
provide isomer-specific estimates to discriminate isomers with slow (β-HCH, HLT=7.2 – 7.6 yrs54) and
178
fast (γ-HCH, HLT=26 h55) elimination.
179
Investigating Factors Influencing Bioamplification
180 181
The “bioamplification factor (BAmF)” is defined as the ratio of lipid-normalized chemical concentrations after (Cafter) and before (Cbefore) weight loss (Equation 5):7
182
BAmF = Cafter / Cbefore,
(Eq. 5)
183
We used the model to calculate BAmF for hypothetical chemicals with octanol-air partitioning
184
coefficients log KOA ranging from 0 to 12 and octanol-water partitioning coefficients log KOW from 0
185
and 12. Internal energies of phase transfer of these hypothetical chemicals were again calculated using
186
the
187
total-blood-volume-adjusted biotransformation half-life HLB′ ranging from 1 h to infinite, and (ii)
188
absolute rate of lipid loss of either 1.5 kg month-1 (representing the range of 0.5 – 1 lb week-1 weight
189
loss recommended for overweight patients15), 3 kg month-1 (representing the range of 1 – 2 lb week-1
190
recommended for obese patients15), or 6 kg month-1 (representing the range of 110 – 220 lb over a period
191
of 6 to 12 months recommended for severely obese patients who receive bariatric surgery15). All other
192
parameters were the same as in our above efforts reproducing literature reported increments. BAmF
193
values after n months of weight loss (BAmFn) are displayed in a chemical partitioning space defined by
194
log KOA and log KOW.
195
Results and discussion
196
Capability of the model to reproduce literature-reported bioamplification
empirical
relationships
in
ref.52
We
tested
different
combinations
of
(i)
the
197
Figure 1 plots the modelled increments in lipid-normalized PLOP concentration (∆Cmodelled) against
198
the increments derived from measured data reported in the literature.9, 10, 12, 24-26 Compound-specific
199
modelled and measured increments are also presented in numeric form in Table S1 in Supporting
200
Information. The modelled mean/median increments are significantly correlated with measurements (r =
9
ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 26
201
0.80, p < 0.05). The root-mean-square deviation (RMSD) between logarithmical measured and modelled
202
pairs, i.e.,
203
RMSD =
1 n 2 ( log ∆Cmeasured − log ∆Cmodelled ) ∑ n i =1
204
is calculated to be 0.31 (6 pairs with non-positive measured increments are excluded because the real
205
logarithm is undefined for non-positive real numbers, see Table S1 in Supporting Information). That
206
RMSD is less than 1 means simulations and measurements are generally within an order of magnitude,
207
i.e., much smaller than the variability associated with the measured data (denoted as the error bars).
208
Recent studies highlighted the complex interactions between the dynamics of PLOPs and lipids
209
during weight loss. For instance, a negative feedback loop has been suggested, whereby a weight
210
loss-induced increase in PLOP concentrations may trigger weight gain, or at least decelerate weight loss,
211
due to the “obesogenic” effect.56, 57 Meanwhile, the induced expression and activation of xenobiotic
212
metabolizing enzymes (e.g., cytochrome P450)27 is believed to enhance the hepatic biotransformation,58
213
which would constitute a positive feedback loop. That our simple one-compartment PK model
214
performed well in explaining the observations, despite ignoring all of these feedbacks, strongly suggests
215
that it is the lipid dynamics that predominantly govern the dynamics of PLOPs rather than vice versa.
216
Since epidemiological studies are often plagued with misinterpretations or misclassifications when
217
involving multiple complicated and variable loops,59 such a simple model helps to recognize the essence
218
of bioamplification, which could be particularly useful for future investigations on the adverse health
219
outcome of PLOP bioamplification.
220
Influences of chemical properties and lipid loss rate on bioamplification
221
Figure 2 displays the calculated BAmF3, i.e., the BAmF at the end of a 3-month weight loss
222
treatment, as a function of partition properties under different scenarios of biotransformation half-lives
223
(horizontal panels in Figure 2) and absolute rates of lipid loss (vertical panels in Figure 2). In order to
224
discuss and understand the results displayed in these plots, it is helpful to introduce a number of
225
additional parameters. Specifically, we can define: 10
ACS Paragon Plus Environment
Page 11 of 26
226 227 228
Environmental Science & Technology
(i) a metabolic elimination rate kelim,metab = DmH / (VHZH), and the corresponding metabolic elimination half time HLelim,metab = ln2 / kelim,metab; (ii) a non-metabolic elimination rate kelim, non-metab =
EOH DeH + DreH + DperH + DlH +DchH + DurH VH ⋅ Z H
, and the
229
corresponding non-metabolic elimination half time: HLelim,non-metab = ln2 / kelim,non-metab. Here, kelim, non-metab
230
is a function of the chemical partitioning properties log KOA and log KOW (see Supporting Information;
231
Figure S2).
232
(iii) a rate describing the mobilization of PLOPs due to lipolysis, which corresponds to the relative
dM lipid
233
rate of lipid loss: kmobil = −
234
mobilization: DLelim,non-metab = ln2/ kmobil.
235 236
M lipid ⋅ dt
, and the doubling time of blood chemical concentration due to
Note that while kelim,metab, kelim,non-metab and kmobil are labeled “rate constants”, they are time variant as lipid is lost (as discussed later). Substituting these rate constants into Equation 4, we get:
237
dC = ( kmobil − kelim,metab − kelim,non-metab ) ⋅ C dt
238
Equation 6 echoes the earlier understanding7 that the change of lipid-normalized chemical
239
concentration depends on two competing processes: the “source”, i.e., the mobilization of PLOPs from
240
adipose tissues into the bloodstream (i.e., kmobil) and the “sink”, i.e., the elimination of PLOPs from
241
blood (i.e., kelim = kelim, metab + kelim, non-metab). Mathematically, bioamplification will occur (BAmF > 1)
242
when kmobil – kelim > 0, or DLmobil < HLelim, which provides a quantitative basis for the statement by Daley
243
et al.21 that “[b]ioamplification can be defined as the special case where the partition capacity (or
244
fugacity capacity) of an animal’s tissues decreases faster than chemical can be eliminated.”
(Eq.6)
245
Under all scenarios displayed in Figure 2, more lipophilic chemicals (i.e., with larger log KOW and
246
log KOA, in the upper right corner) are prone to possess higher BAmF3 values (Figure 2). Specifically,
247
bioamplification occurs (i.e. BAmF3 > 1) (Figure 2b,c,d,f,g,h,k,j,l), for compounds with log KOW > 2 and
248
log KOA > 6, which coincides with the part of the partitioning space in which non-metabolic elimination
249
is inefficient (HLelim, non-metab >104 h; see Supporting Information; Figure S2). This finding agrees well 11
ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 26
250
with a conclusion drawn from animal studies that “for terrestrial food webs, chemicals that possess a log
251
KOW > 2 and a log KOA > 6 should be regarded as having a high bioamplification potential”.7 Soluble
252
chemicals with a log KOW 2 and log KOA > 6 varies a lot among scenarios, in the remaining part of the partitioning space
256
BAmF3 remains almost the same under all scenarios.
257
At the contour lines indicating BAmF3 =1 mobilization and elimination processes are in balance and
258
no change in contaminant concentration occurs. Mobilization prevails in the region BAmF3 >1 (upper
259
right in Figure 2) while elimination prevails in the region BAmF3 5 and
265
log KOA > 6) are almost constant (Figure 2), because their non-metabolic elimination half times (HLelim,
266
non-metab
=104 – 106 h) far exceed mobilization doubling times (DLmobil = 103.5 – 104 h, assuming that an
267
individual with an initial total lipid mass of 40 kg loses 1.5 – 6 kg lipid per month) and thus the
268
mobilization process is the dominant determinant for bioamplification. The finding that PLOPs share
269
similar BAmF values is also borne out in observations in human monitoring studies.10, 24 For instance,
270
Chevrier et al.24 observed similar percent increments between 10% and 20% in plasma concentrations of
271
14 of 17 detected chemicals after weight loss in a 39-individual population. This is also true for animal
272
studies. For example, the depletion of yolk lipids during hatching was observed to amplify PLOP
273
concentrations in hatched Chinook salmon larvae, with the BAmF between fresh eggs and larvae
274
increasing when log KOW exceeded 5.5 “until it reached an asymptote at log KOW values exceeding 7”.20
12
ACS Paragon Plus Environment
Page 13 of 26
Environmental Science & Technology
275
Raising the absolute rate of weight loss (moving from top to bottom in Figure 2) increases the
276
BAmF3 value for all chemicals in the partitioning space because of elevated mobilization. This
277
simulation explains an earlier finding that for most compounds the increment in plasma chemical
278
concentration is significantly positively correlated with the reduction in body weight.12, 24 The increment
279
is more profound for the most lipophilic PLOPs (upper right corner in Figure 2) because of inefficient
280
non-metabolic elimination (see Supporting Information; Figure S2). As the BAmF3 increases, the
281
contour line of BAmF3 =1 slightly moves towards the bottom left, which means that the faster lipid is
282
lost the more chemicals will exhibit bioamplification.
283
As expected, increasing resistance to metabolic elimination (moving from left to right in Figure 2)
284
increases BAmF3 for all chemical partitioning property combinations. The increase in BAmF3 is most
285
prominent for the highly lipophilic chemicals in the upper right of the partitioning space. As indicated by
286
the prevalence of green and yellow shadings in Figures 2a, 2e and 2i, readily degradable chemicals with
287
a HLB′ < 10 hours (equivalent to an HLelim,metab of 2 and log KOA > 6) are quite different depending on whether HLB′ is 1 or 10 hours (equivalent to
290
an HLelim,metab of 103 or 104 hours) (Figures 2a, b, e, f, i, and j), but are fairly similar for HLB′ of 100
291
hours (equivalent to an HLelim,metab of 105 hours) or infinite (Figures 2c, d, g, h, k, and l), which suggests
292
that the BAmF3 of PLOPs is no longer sensitive to HLB′ if HLB′ is sufficiently long. Here is the
293
explanation: Given that their HLelim,non-metab values are around 104 – 106 hours (see Supporting
294
Information; Figure S1), most PLOPs cannot be efficiently eliminated from humans via non-metabolic
295
elimination; thus, biotransformation dominates total elimination if HLelim,metab is shorter than 104 hours,
296
hence rendering the BAmF3 of the lipophilic chemicals highly sensitive to the HLelim,metab. This
297
conclusion implies that an accurate HLB′ value is crucial for quantifying bioamplification for
298
non-persistent compounds.
13
ACS Paragon Plus Environment
Environmental Science & Technology
299
Page 14 of 26
Temporal evolution of the bioamplification factor
300
The BAmF changes with time because both mobilization and elimination are dynamic as weight loss
301
proceeds (Figure 3). As indicated by Eqs 4 and 6, on the one hand, while we assume the absolute rate of
302
lipid loss
dM lipid
dt
to be constant over the period of weight loss, the relative rate
dM lipid
M lipid ⋅ dt
keeps
303
increasing as Mlipid decreases, which means that mobilization (kmobil) accelerates. On the other hand, as
304
lipid is lost, the total effective volume (VH) and the Z-value (ZH) of an individual decrease, which
305
accelerates both metabolic and non-metabolic elimination (kelim). Accelerated elimination in leaner
306
individuals is supported by earlier observations, e.g., the decreasing yearly “decay rate” of
307
polychlorinated dibenzo-p-dioxins and dibenzofurans against increasing BF% in an occupationally
308
exposed population (N = 48).60
309
The deduction in Text S2 in Supporting Information indicates that both the mobilization (kmobil) and
310
elimination (kelim) rates are roughly proportional to the reciprocal total lipid mass (Mlipid), thus the loss of
311
a given amount of lipid will increase the absolute values of both rates by the same factor. Hence, if the
312
absolute values of the two rates are not equal, the one with a larger initial absolute value will increase
313
more; the accumulation of larger increments further enlarges the difference between the two rates. As
314
such, we can predict that: (i) if a compound falls into the elimination-dominated part of the partitioning
315
space and has a BAmF below 1 at the start of lipid loss, the BAmF will keep decreasing throughout the
316
entire weight loss period, i.e., BAmF1 > BAmF3 > BAmF6 (Figure 3a); (ii) if a compound falls into the
317
mobilization-dominated part of the partitioning space and has BAmF above 1 at the start of lipid loss, the
318
BAmF will keep increasing throughout the entire weight loss period, i.e., BAmF1 < BAmF3 < BAmF6
319
(Figure 3c); and (iii) the BAmF will remain almost unchanged if its initial value is close to 1 (Figure 3b).
320
This implies that, while the value of the BAmF is time-variant, whether a chemical is able to bio-amplify
321
in humans or not is independent of the length of weight loss. If bioamplification does occur, the longer
322
the weight loss duration the larger the observed BAmF.
14
ACS Paragon Plus Environment
Page 15 of 26
Environmental Science & Technology
323
From Equation 6, it is clear that the derivative of the BAmF curve, i.e., the change in BAmF per unit
324
of time (∆BAmF/∆t), relies not only on the relative importance of the two rates kmobil and kelim, but also
325
on the time-dependent lipid-normalized chemical concentration C. When elimination dominates (Figure
326
3a), the concentration declines with time, which in turn decelerates the decline in BAmF and BAmF
327
asymptotically approaches 0 (a negative-feedback-loop-like correlation); in contrast, when mobilization
328
dominates (Figure 3c), the concentration keeps increasing as lipid is lost, which further accelerates the
329
increase in BAmF (a positive-feedback-loop-like correlation). Therefore, the time-dependent BAmF
330
curves in both cases are concave and have increasing derivatives.
331
Implications for future studies
332
The success of a modified and well-parameterized one-compartment PK model to simulate
333
bioamplification of PLOPs in humans has profound implications for future studies. For example, it can
334
serve as a simple and fast predictive tool for batch screening chemicals for their bioamplification
335
potential. For specific compounds (e.g., 2,3,7,8-tetrachlorodibenzo-p-dioxin) more sophisticated PK
336
models with additional physiological details (e.g., induction of a hepatic binding protein, and diffusion
337
limitation for movement of chemical between blood and tissue) will be required.58
338
Second, while we focused here on intentional weight loss, the model can be applied to investigate
339
the dynamics of PLOPs in the human body during unintentional, in particular illness-related, weight loss,
340
as well as its association with obesity-related disease. Such a simulation could help to avoid potential
341
confounders or reverse causation in developing exposure-outcome association in epidemiological
342
studies.61 For instance, future modelling studies may clarify whether high concentrations of PLOPs
343
observed in diagnosed diabetics are really epidemiological evidence of the endocrine disruptive effects
344
of the PLOPs, or just due to rapid bioamplification during the pathological weight reduction in the
345
pre-diagnostic phase of pre-existing illness.62
346
Finally, while we deemed ongoing dietary or respiratory exposure negligible during the short-term
347
simulations presented here, it could be included in long-term studies so as to assess its contribution to
348
the chemical burden in human. Such a simulation could help to elucidate the “apparent” 15
ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 26
349
bioamplification of non-persistent chemicals (e.g., bisphenol A, triclosan and phthalates) due to ongoing
350
exposure,63, 64 and to explore differences in “cross-sectional body burden versus BMI trends” of diverse
351
PLOPs due to their distinct emission histories and cumulative external exposures.59, 65
352
Acknowledgement
353
We thank Dr. Jon A. Arnot for QSAR calculation of whole-body biotransformation half-lives in humans.
354
Funding was provided by the Natural Sciences and Engineering Research Council of Canada.
355
Supporting information
356
The Supporting Information is available free of charge on the ACS Publications website at DOI:…
357
Text describing deduction of the total-blood-volume-adjusted biotransformation half-life, and the
358
approximate dependence on total lipid mass of mobilization and elimination rates; tables listing
359
literature-reported bioamplification measurements; and figures depicting results of sensitivity analysis,
360
and the chemical space of non-metabolic elimination half-life.
361
References
362
1.
La Merrill, M.; Emond, C.; Kim, M. J.; Antignac, J.-P.; Le Bizec, B.; Clément, K.; Birnbaum, L. S.;
363
Barouki, R., Toxicological function of adipose tissue: focus on persistent organic pollutants.
364
Environ. Health Perspect. 2013, 121, (2), 162.
365
2.
Malarvannan, G.; Dirinck, E.; Dirtu, A. C.; Pereira-Fernandes, A.; Neels, H.; Jorens, P. G.; Gaal, L.
366
V.; Blust, R.; Covaci, A., Distribution of persistent organic pollutants in two different fat
367
compartments from obese individuals. Environ. Int. 2013, 55, 33-42.
368
3.
pollutants. Physiol. Res. 2007, 56, (4), 375-381.
369 370 371
Mullerova, D.; Kopecky, J., White adipose tissue: Storage and effector site for environmental
4.
Brown, J. F.; Lawton, R. W., Polychlorinated biphenyl (PCB) partitioning between adipose tissue and serum. Bull. Environ. Contam. Toxicol. 1984, 33, (1), 277-280. 16
ACS Paragon Plus Environment
Page 17 of 26
372
5.
373 374
Environmental Science & Technology
Barouki, R., Linking long-term toxicity of xeno-chemicals with short-term biological adaptation. Biochimie 2010, 92, (9), 1222-1226.
6.
Louis, C.; Van den Daelen, C.; Tinant, G.; Bourez, S.; Thomé, J.-P.; Donnay, I.; Larondelle, Y.;
375
Debier, C., Efficient in vitro adipocyte model of long-term lipolysis: A tool to study the behavior of
376
lipophilic compounds. In Vitro Cell. Dev. Biol. - Animal 2014, 50, (6), 507-518.
377
7.
Daley, J. M.; Paterson, G.; Drouillard, K. G., Bioamplification as a bioaccumulation mechanism for
378
persistent organic pollutants (POPs) in wildlife. In Reviews of Environmental Contamination and
379
Toxicology, Whitacre, D. M., Ed. Springer: Switzerland, 2014; Vol. 227, pp 107-155.
380
8.
environment. Environ. Sci. Technol. 2002, 36, 457A-462A.
381 382
Macdonald, R.; Mackay, D.; Hickie, B., Peer reviewed: Contaminant amplification in the
9.
Pelletier, C.; Doucet, E.; Imbeault, P.; Tremblay, A., Associations between weight loss-induced
383
changes in plasma organochlorine concentrations, serum T3 concentration, and resting metabolic
384
rate. Toxicol. Sci. 2002, 67, (1), 46-51.
385
10. Imbeault, P.; Chevrier, J.; Ayotte, P.; Despres, J.; Mauriege, P.; Tremblay, A., Increase in plasma
386
pollutant levels in response to weight loss is associated with the reduction of fasting insulin levels in
387
men but not in women. Metabolism 2002, 51, (4), 482-486.
388
11. Lim, J. S.; Son, H. K.; Park, S. K.; Jacobs, D. R., Jr.; Lee, D. H., Inverse associations between
389
long-term weight change and serum concentrations of persistent organic pollutants. Int. J. Obes.
390
2011, 35, (5), 744-747.
391
12. Dirinck, E.; Dirtu, A.; Jorens, P.; Malarvannan, G.; Covaci, A.; Van Gaal, L., Pivotal role for the
392
visceral fat compartment in the release of persistent organic pollutants during weight loss. J. Clin.
393
Endocrinol. Metab. 2015, 100, (12), 4463-4471.
394
13. Lignell, S.; Winkvist, A.; Bertz, F.; Rasmussen, K. M.; Glynn, A.; Aune, M.; Brekke, H. K.,
395
Environmental organic pollutants in human milk before and after weight loss. Chemosphere 2016,
396
159, 96-102.
17
ACS Paragon Plus Environment
Environmental Science & Technology
Page 18 of 26
397
14. Vaz, R.; Slorach, S. A.; Hofvander, Y., Organochlorine contaminants in Swedish human milk:
398
Studies conducted at the national food administration 1981–1990. Food Addit. Contam. 1993, 10,
399
(4), 407-418.
400
15. US National Heart Lung and Blood Institute Clinical Guidelines on the Identification, Evaluation,
401
and Treatment of Overweight and Obesity in Adults: The Evidence Report (Report No.: 98-4083);
402
1998.
403 404 405 406 407 408 409 410
16. Fitzhugh, O. G.; Nelson, A. A., The chronic oral toxicity of DDT (2, 2-bis (p-chlorophenyl-1, 1, 1-trichloroethane). J. Pharmacol. Exp. Ther. 1947, 89, (1), 18-30. 17. Findlay, G. M.; DeFreitas, A. S. W., DDT Movement from adipocyte to muscle cell during lipid utilization. Nature 1971, 229, 63-65. 18. Deichmann, W. B.; Beasley, A. G.; Cubit, D. A.; Macdonald, W. E., Effects of starvation in rats with elevated DDT and Dieldrin tissue levels. Int. Arch. Arbeitsmed. 1972, 29, (3), 233-252. 19. Ohmiya, Y.; Nakai, K., Effect of starvation on excretion, distribution and metabolism of DDT in mice. Tohoku J. Exp. Med. 1977, 122, (2), 143-153.
411
20. Daley, J. M.; Leadley, T. A.; Pitcher, T. E.; Drouillard, K. G., Bioamplification and the selective
412
depletion of persistent organic pollutants in Chinook Salmon larvae. Environ. Sci. Technol. 2012, 46,
413
(4), 2420-2426.
414
21. Daley, J. M.; Leadley, T. A.; Drouillard, K. G., Evidence for bioamplification of nine
415
polychlorinated biphenyl (PCB) congeners in yellow perch (Perca flavascens) eggs during
416
incubation. Chemosphere 2009, 75, (11), 1500-1505.
417
22. Tartu, S.; Bourgeon, S.; Aars, J.; Andersen, M.; Polder, A.; Thiemann, G. W.; Welker, J. M.; Routti,
418
H., Sea ice-associated decline in body condition leads to increased concentrations of lipophilic
419
pollutants in polar bears (Ursus maritimus) from Svalbard, Norway. Sci. Total Environ. 2017, 576,
420
409-419.
18
ACS Paragon Plus Environment
Page 19 of 26
Environmental Science & Technology
421
23. Backman, L.; Kolmodin-Hedman, B., Concentration of DDT and DDE in plasma and subcutaneous
422
adipose tissue before and after intestinal bypass operation for treatment of obesity. Toxicol. Appl.
423
Pharm. 1978, 46, (3), 663-669.
424
24. Chevrier, J.; Dewailly, E.; Ayotte, P.; Mauriege, P.; Despres, J.; Tremblay, A., Body weight loss
425
increases plasma and adipose tissue concentrations of potentially toxic pollutants in obese
426
individuals. Int. J. Obes. 2000, 24, (10), 1272-1278.
427
25. Dirtu, A. C.; Dirinck, E.; Malarvannan, G.; Neels, H.; Van Gaal, L.; Jorens, P. G.; Covaci, A.,
428
Dynamics of organohalogenated contaminants in human serum from obese individuals during one
429
year of weight loss treatment. Environ. Sci. Technol. 2013, 47, 12441-12449.
430
26. Charlier, C.; Desaive, C.; Plomteux, G., Human exposure to endocrine disrupters: consequences of
431
gastroplasty on plasma concentration of toxic pollutants. Int. J. Obes. 2002, 26, (11), 1465-1468.
432
27. Kim, M.-J.; Marchand, P.; Henegar, C.; Antignac, J.-P.; Alili, R.; Poitou, C.; Bouillot, J.-L.;
433
Basdevant, A.; Le Bizec, B.; Barouki, R., Fate and complex pathogenic effects of dioxins and
434
polychlorinated biphenyls in obese subjects before and after drastic weight loss. Environ. Health
435
Perspect. 2011, 119, (3), 377.
436
28. Hue, O.; Marcotte, J.; Berrigan, F.; Simoneau, M.; Doré, J.; Marceau, P.; Marceau, S.; Tremblay, A.;
437
Teasdale, N., Increased plasma levels of toxic pollutants accompanying weight loss induced by
438
hypocaloric diet or by bariatric surgery. Obes. Surg. 2006, 16, (9), 1145-1154.
439
29. Norstrom, R. J.; Clark, T. P.; Enright, M.; Leung, B.; Drouillard, K. G.; Macdonald, C. R., ABAM, a
440
model for bioaccumulation of POPs in birds: validation for adult herring gulls and their eggs in
441
Lake Ontario. Environ. Sci. Technol. 2007, 41, (12), 4339-4347.
442
30. Drouillard, K. G.; Norstrom, R. J.; Fox, G. A.; Gilman, A.; Peakall, D. B., Development and
443
validation of a herring gull embryo toxicokinetic model for PCBs. Ecotoxicology 2003, 12, (1),
444
55-68.
19
ACS Paragon Plus Environment
Environmental Science & Technology
Page 20 of 26
445
31. Drouillard, K. G.; Paterson, G.; Haffner, G. D., A combined food web toxicokinetic and species
446
bioenergetic model for predicting seasonal PCB elimination by Yellow Perch (Perca flavescens).
447
Environ. Sci. Technol. 2009, 43, (8), 2858-2864.
448 449 450 451
32. Czub, G.; McLachlan, M. S., A food chain model to predict the levels of lipophilic organic contaminants in humans. Environ. Toxicol. Chem. 2004, 23, (10), 2356-2366. 33. Undeman, E.; Brown, T. N.; Wania, F.; McLachlan, M. S., Susceptibility of human populations to environmental exposure to organic contaminants. Environ. Sci. Technol. 2010, 44, (16), 6249-6255.
452
34. Rawn, D. F. K.; Cao, X. L.; Doucet, J.; Davies, D. J.; Sun, W. F.; Dabeka, R. W.; Newsome, W. H.,
453
Canadian Total Diet Study in 1998: Pesticide levels in foods from Whitehorse, Yukon, Canada, and
454
corresponding dietary intake estimates. Food Addit. Contam. 2004, 21, (3), 232-250.
455 456
35. Domingo, J. L., Polybrominated diphenyl ethers in food and human dietary exposure: A review of the recent scientific literature. Food Chem. Toxicol. 2012, 50, (2), 238-249.
457
36. Ritter, R.; Scheringer, M.; MacLeod, M.; Moeckel, C.; Jones, K. C.; Hungerbühler, K., Intrinsic
458
human elimination half-lives of polychlorinated biphenyls derived from the temporal evolution of
459
cross-sectional biomonitoring data from the United Kingdom. Environ. Health Perspect. 2011, 119,
460
(2), 225-231.
461
37. Sprenger, K. B. G.; Huber, K.; Kratz, W.; Henze, E., Nomograms for the prediction of patient's
462
plasma volume in plasma exchange therapy from height, weight, and hematocrit. J. Clin. Apheresis
463
1987, 3, (3), 185-190.
464
38. Haddad, S.; Poulin, P.; Krishnan, K., Relative lipid content as the sole mechanistic determinant of
465
the adipose tissue:blood partition coefficients of highly lipophilic organic chemicals. Chemosphere
466
2000, 40, (8), 839-843.
467
39. Endo, S.; Brown, T. N.; Goss, K.-U., General model for estimating partition coefficients to
468
organisms and their tissues using the biological compositions and polyparameter linear free energy
469
relationships. Environ. Sci. Technol. 2013, 47, 6630-6639.
20
ACS Paragon Plus Environment
Page 21 of 26
Environmental Science & Technology
470
40. Artacho-Cordon, F.; Fernandez-Rodriguez, M.; Garde, C.; Salamanca, E.; Iribarne-Duran, L. M.;
471
Torne, P.; Exposito, J.; Papay-Ramirez, L.; Fernandez, M. F.; Olea, N.; Arrebola, J. P., Serum and
472
adipose tissue as matrices for assessment of exposure to persistent organic pollutants in breast
473
cancer patients. Environ. Res. 2015, 142, 633-643.
474
41. Arrebola, J. P.; Cuellar, M.; Claure, E.; Quevedo, M.; Antelo, S. R.; Mutch, E.; Ramirez, E.;
475
Fernandez, M. F.; Olea, N.; Mercado, L. A., Concentrations of organochlorine pesticides and
476
polychlorinated biphenyls in human serum and adipose tissue from Bolivia. Environ. Res. 2012, 112,
477
40-47.
478
42. Whitcomb, B. W.; Schisterman, E. F.; Buck, G. M.; Weiner, J. M.; Greizerstein, H.; Kostyniak, P. J.,
479
Relative concentrations of organochlorines in adipose tissue and serum among reproductive age
480
women. Environ. Toxicol. Pharmacol. 2005, 19, (2), 203-213.
481
43. Verner, M.-A.; McDougall, R.; Glynn, A.; Andersen, M. E.; Clewell, H. J.; Longnecker, M. P., Is the
482
relationship between prenatal exposure to PCB-153 and decreased birth weight attributable to
483
pharmacokinetics? Environ. Health Perspect. 2013, 121, (10), 1219-1224.
484
44. Ehresman, D. J.; Froehlich, J. W.; Olsen, G. W.; Chang, S.-C.; Butenhoff, J. L., Comparison of
485
human whole blood, plasma, and serum matrices for the determination of perfluorooctanesulfonate
486
(PFOS), perfluorooctanoate (PFOA), and other fluorochemicals. Environ. Res. 2007, 103, (2),
487
176-184.
488
45. Mes, J.; Marchand, L.; Turton, D.; Lau, P. Y.; Ganz, P. R., The determination of polychlorinated
489
biphenyl congeners and other chlorinated hydrocarbon residues in human blood, serum and plasma.
490
A comparative study. Int. J. Environ. Analyt. Chem. 1992, 48, (3-4), 175-186.
491
46. Price, P. S.; Conolly, R. B.; Chaisson, C. F.; Gross, E. A.; Young, J. S.; Mathis, E. T.; Tedder, D. R.,
492
Modeling interindividual variation in physiological factors used in PBPK models of humans. Crit.
493
Rev. Toxicol. 2003, 33, (5), 469-503.
494 495
47. Deurenberg, P.; Yap, M.; van Staveren, W. A., Body mass index and percent body fat: a meta-analysis among different ethnic groups. Int. J. Obes. 1998, 22, (12), 1164-1171.
21
ACS Paragon Plus Environment
Environmental Science & Technology
496 497
Page 22 of 26
48. Mackay, D.; Shiu, W. Y.; Ma, K. C.; Lee, S. C., Handbook of Physical-Chemical Properties and Environmental Fate for Organic Chemicals (Second Edition). CRC Press: 2006.
498
49. Li, N.; Wania, F.; Lei, Y. D.; Daly, G. L., A comprehensive and critical compilation, evaluation, and
499
selection of physical–chemical property data for selected polychlorinated biphenyls. J. Phys. Chem.
500
Ref. Data 2003, 32, (4), 1545-1590.
501 502 503 504 505 506 507 508
50. Shen, L.; Wania, F., Compilation, evaluation, and selection of physical−chemical property data for organochlorine pesticides. J. Chem. Eng. Data 2005, 50, (3), 742-768. 51. Xiao, H.; Li, N.; Wania, F., Compilation, evaluation, and selection of physical-chemical property data for α-, β-, and γ-Hexachlorocyclohexane. J. Chem. Eng. Data 2004, 49, (2), 173-185. 52. MacLeod, M.; Scheringer, M.; Hungerbuhler, K., Estimating enthalpy of vaporization from vapor pressure using Trouton's rule. Environ. Sci. Technol. 2007, 41, (8), 2827-2832. 53. Arnot, J. A.; Brown, T. N.; Wania, F., Estimating screening-level organic chemical half-lives in humans. Environ. Sci. Technol. 2014, 48, 723-730.
509
54. Jung, D.; Becher, H.; Edler, L.; Flesch-Janys, D.; Gurn, P.; Konietzko, J.; Manz, A.; Papke, O.,
510
Elimination of β-hexachlorocyclohexane in occupationally exposed persons. J. Toxicol. Environ.
511
Health 1997, 51, (1), 23-34.
512 513 514 515
55. Feldmann, R. J.; Maibach, H. I., Percutaneous penetration of some pesticides and herbicides in man. Toxicol. Appl. Pharm. 1974, 28, (1), 126-132.
56. Holtcamp, W., Obesogens: an environmental link to obesity. Environ. Health Perspect. 2012, 120, (2), A63-A68.
516
57. Lee, D.-H.; Steffes, M. W.; Sjödin, A.; Jones, R. S.; Needham, L. L.; Jacobs, D. R., Jr., Low dose
517
organochlorine pesticides and polychlorinated biphenyls predict obesity, dyslipidemia, and insulin
518
resistance among people free of diabetes. PLOS ONE 2011, 6, (1), e15977.
519
58. Emond, C.; Birnbaum, L. S.; DeVito, M. J., Use of a physiologically based pharmacokinetic model
520
for rats to study the influence of body fat mass and induction of CYP1A2 on the pharmacokinetics
521
of TCDD. Environ. Health Perspect. 2006, 114, (9), 1394-1400.
22
ACS Paragon Plus Environment
Page 23 of 26
Environmental Science & Technology
522
59. Wolff, M. S.; Anderson, H. A.; Britton, J. A.; Rothman, N., Pharmacokinetic variability and modern
523
epidemiology—the example of dichlorodiphenyltrichloroethane, body mass index, and birth cohort.
524
Cancer Epidemiol. Biomarkers Prev. 2007, 16, (10), 1925-1930.
525 526
60. Flesch-Janys, D., Elimination of polychlorinated dibenzo-p-dioxins and dibenzofurans in occupationally exposed persons. J. Toxicol. Environ. Health 1996, 47, (4), 363-378.
527
61. Flegal, K. M.; Graubard, B. I.; Williamson, D. F.; Cooper, R. S., Reverse causation and
528
illness-related weight loss in observational studies of body weight and mortality. Am. J. Epidemiol.
529
2011, 173, (1), 1-9.
530
62. Kerger, B. D.; Scott, P. K.; Pavuk, M.; Gough, M.; Paustenbach, D. J., Re-analysis of Ranch Hand
531
study supports reverse causation hypothesis between dioxin and diabetes. Crit. Rev. Toxicol. 2012,
532
42, (8), 669-687.
533
63. Dirtu, A. C.; Geens, T.; Dirinck, E.; Malarvannan, G.; Neels, H.; Van Gaal, L.; Jorens, P. G.; Covaci,
534
A., Phthalate metabolites in obese individuals undergoing weight loss: Urinary levels and estimation
535
of the phthalates daily intake. Environ. Int. 2013, 59, 344-353.
536
64. Geens, T.; Dirtu, A. C.; Dirinck, E.; Malarvannan, G.; Van Gaal, L.; Jorens, P. G.; Covaci, A., Daily
537
intake of bisphenol A and triclosan and their association with anthropometric data, thyroid
538
hormones and weight loss in overweight and obese individuals. Environ. Int. 2015, 76, 98-105.
539
65. Wood, S.; Xu, F.; Armitage, J. M.; Wania, F., Unravelling the relationship between body mass index
540
and polychlorinated biphenyl concentrations using a mechanistic model. Environ. Sci. Technol. 2016,
541
50, 10055-10064.
542
23
ACS Paragon Plus Environment
Environmental Science & Technology
Page 24 of 26
543
Figure Legends
544
Figure 1 Comparison between measured and modelled increments in blood concentrations of PLOPs.
545
Dots represent means (if means were reported in the cited literature) or medians (if medians were
546
reported in the cited literature) of the increments, and error bars represent semi-interquartile range
547
(intervals between the 75th and 25th percentiles) of the increments (reported in the cited literature
548
or converted from the data in the literature). Dashed diagonal lines denote agreement between
549
measurements and model within a factor of 10. Results for mirex9,24 cannot be displayed because
550
their measured increments are negligible.
551
552
24
ACS Paragon Plus Environment
Page 25 of 26
Environmental Science & Technology
553
Figure 2 The bioamplification factor (BAmF) of hypothetical chemicals as a function of partitioning
554
coefficients (log KOW and log KOA), absolute rate of lipid loss (in kg month-1), and
555
total-blood-volume-adjusted biotransformation half-life (in h).
556
557
25
ACS Paragon Plus Environment
Environmental Science & Technology
Page 26 of 26
558
Figure 3 Time courses of bioamplification factor (BAmF), elimination (kelim = kelim, metab + kelim, non-metab)
559
and mobilization (kmobil) rates of hypothetical chemicals (a) falling in the elimination-dominated
560
partitioning space (log KOW = log KOA = 4.45), (b) with equal elimination and mobilization rates
561
(log KOW = log KOA = 5.45), and (c) falling in the mobilization-dominated partitioning space (log
562
KOW = log KOA = 6.45). An initial total lipid mass of Mlipid = 40 kg and a lipid loss rate of 3 kg
563
month-1 was applied.
564
26
ACS Paragon Plus Environment