Deriving in vivo biotransformation rate constants and metabolite parent concentration factor/stable metabolite factor from bioaccumulation and bioconcentration experiments: An illustration with worm accumulation data

Organophosphate esters and their metabolites in silver pomfret (Pampus argenteus) of the Vietnamese coastal areas: Spatial-temporal distribution and exposure risk

A large number of studies on organophosphate esters (tri-OPEs) in marine organisms have not assessed the simultaneous occurrence of tri-OPEs and their metabolites (di-OPEs) in these species. This research investigated the concentration and geographical distribution of 15 tri-OPEs and 7 di-OPEs in 172 samples of Pampus argenteus that were collected annually from 2021 to 2023 at three distinct locations along the Vietnamese coast. As a result, tri-OPEs and di-OPEs were detected in numerous fish samples, indicating their widespread spatial and temporal occurrence in marine fish and pointing out the importance of monitoring their levels. The tri-OPEs and di-OPEs ranged within 2.1-38.9 ng g-1 dry weight (dw) and 3.2-263.4 ng g-1 dw, respectively. The mean concentrations of tri-OPEs ranged from 0.4 (TIPrP) to 5.4 ng g-1 dw (TBOEP), with TBOEP and TEHP having the highest mean values. In addition, the profiles of tri-OPEs in fish exhibited a descending order: Σalkyl OPEs > ΣCl-alkyl OPEs > Σaryl OPEs. The di-OPEs, namely BEHP and DMP, had the highest mean levels, measuring 33.4 ng g-1 dw and 23.8 ng g-1 dw, respectively. Furthermore, there have been significant findings of strong positive correlations between di-OPEs and tri-OPE pairs (p < 0.05). It is worth noting that there is a noticeable difference in the composition of tri-OPEs between the North and other regions. Despite these findings, the presence of OPE-contaminated fish did not pose any health risks to Vietnam's coastal population.

First insights into the bioaccumulation, biotransformation and trophic transfer of typical tetrabromobisphenol A (TBBPA) analogues along a simulated aquatic food chain

Tetrabromobisphenol A (TBBPA) analogues have been investigated for their prevalent occurrence in environments and potential hazardous effects to humans and wildlife; however, there is still limited knowledge regarding their toxicokinetics and trophic transfer in aquatic food chains. Using a developed toxicokinetic model framework, we quantified the bioaccumulation, biotransformation and trophic transfer of tetrabromobisphenol S (TBBPS) and tetrabromobisphenol A di(allyl ether) (TBBPA-DAE) during trophic transfer from brine shrimp (Artemia salina) to zebrafish (Danio rerio). The results showed that the two TBBPA analogues could be readily accumulated by brine shrimp, and the estimated bioconcentration factor (BCF) value of TBBPS (5.68 L kg-1 ww) was higher than that of TBBPA-DAE (1.04 L kg-1 ww). The assimilation efficiency (AE) of TBBPA-DAE in zebrafish fed brine shrimp was calculated to be 16.3%, resulting in a low whole-body biomagnification factor (BMF) in fish (0.684 g g-1 ww). Based on the transformation products screened using ultra-high-performance liquid chromatograph-high resolution mass spectrometry (UPLC-HRMS), oxidative debromination and hydrolysis were identified as the major transformation pathways of TBBPS, while the biotransformation of TBBPA-DAE mainly took place through ether bond breaking and phase-II metabolism. Lower accumulation of TBBPA as a metabolite than its parent chemical was observed in both brine shrimp and zebrafish, with metabolite parent concentration factors (MPCFs) < 1. The investigated BCFs for shrimp of the two TBBPA analogues were only 3.77 × 10-10 - 5.59 × 10-3 times of the theoretical Kshrimp-water based on the polyparameter linear free energy relationships (pp-LFERs) model, and the BMF of TBBPA-DAE for fish was 0.299 times of the predicted Kshrimp-fish. Overall, these results indicated the potential of the trophic transfer in bioaccumulation of specific TBBPA analogues in higher trophic-level aquatic organisms and pointed out biotransformation as an important mechanism in regulating their bioaccumulation processes. ENVIRONMENTAL IMPLICATION: The internal concentration of a pollutant in the body determines its toxicity to organisms, while bioaccumulation and trophic transfer play important roles in elucidating its risks to ecosystems. Tetrabromobisphenol A (TBBPA) analogues have been extensively investigated for their adverse effects on humans and wildlife; however, there is still limited knowledge regarding their toxicokinetics and trophic transfer in aquatic food chains. This study investigated the bioaccumulation, biotransformation and trophic transfer of TBBPS and TBBPA-DAE in a simulated di-trophic food chain. This state-of-art study will provide a reference for further research on this kind of emerging pollutant in aquatic environments.

Organophosphate diesters (DAPs) and hydroxylated organophosphate flame retardants (HO-OPFRs) as biomarkers of OPFR contamination in a typical freshwater food chain

Organophosphate flame retardants (OPFRs) can rapidly biotransform into two types of metabolites in biota: (1) organophosphate diesters (DAPs) and (2) hydroxylated OPFRs (HO-OPFRs). Therefore, the levels of parent OPFRs alone are not sufficient to indicate OPFR pollution in biological organisms. This study analyzed 12 OPFR metabolites, including 6 DAPs and 6 HO-OPFRs, in a typical freshwater food chain consisted of crucian carp, catfish, mud carp, snakehead, and oriental river prawn. The total concentrations of OPFR metabolites were comparable to those of parent OPFRs, and ranged from 0.65 to 17 ng/g ww. Bis(2-butoxyethyl) 3'-hydroxy-2-butoxyethyl phosphate (14%-77%), di-n-butyl phosphate (DNBP) (6.7%-24%), bis(1-chloro-2-propyl) phosphate (BCIPP) (0.7%-35%), and 1-hydroxy-2-propyl bis(1-chloro-2-propyl) phosphate (BCIPHIPP) (6.0%-24%) were the major OPFR metabolites. Various aquatic species exhibited significant differences in their OPFR metabolite/parent ratios (MPR) (p < 0.05), indicating varying biotransformation potentials of different organisms for various OPFRs. The growth-independent accumulation of tri-n-butyl phosphate (TNBP), tris(chloro-2-propyl) phosphate (TCIPP), triphenyl phosphate, and 2-ethylhexyl diphenyl phosphate in mud carps could be explained by their biotransformation potential. A significant negative correlation was found between the concentration of bis(2-butoxyethyl) phosphate and δ15N values (p < 0.05), with a calculated trophic magnification factor (TMF) of 0.66. Significant positive correlations were observed between BCIPP and TCIPP (R2 = 0.25, p < 0.05), as well as between DNBP and TNBP (R2 = 0.30, p < 0.01), implying that these two DAPs could be used as biomarkers to quantitatively assess TCIPP and TNBP contamination in wild aquatic organisms.

Bioaccumulation and Trophic Transfer of Organophosphate Flame Retardants and Their Metabolites in the Estuarine Food Web of the Pearl River, China

The accumulation and trophodynamics of organophosphate flame retardants (OPFRs) and their metabolites were investigated in the estuarine food web of the Pearl River, China. The mean ∑OPFR concentration among the investigated species increased in the following order: fish [431 ± 346 ng/g lipid weight (lw)] < snail (1310 ± 621 ng/g lw) < shrimp (1581 ± 1134 ng/g lw) < crab (1744 ± 1397 ng/g lw). The di-alkyl phosphates (DAPs) of di-(n-butyl) phosphate (DNBP), bis(2-butoxyethyl) phosphate (BBOEP), and diphenyl phosphate (DPHP) were the most abundant metabolites, with concentrations same as or even higher than their corresponding parent compounds. The log bioaccumulation factors for most OPFRs were lower than 3.70, and significant biomagnification was only found for trisphenyl phosphate [TPHP, with the trophic magnification factors (TMFs) > 1]. The TMFs of OPFRs, except for TPHP and tributyl phosphate had a positive correlation with lipophilicity (log K_OW, p ≤ 0.05) and a negative correlation with the biotransformation rate (log K_M, p ≤ 0.05). The mean TMF > 1 was observed for all of the OPFR metabolites based on the bootstrap regression method. The "pseudo-biomagnification" of OPFR metabolites might be attributed to the biotransformation of OPFRs in organisms at high trophic levels.

Determination of In Vivo Biotransformation Kinetics Using Early-Time Biota Concentrations

Technical challenges have hampered the characterization of biotransformation kinetics-a critical link in understanding and predicting the toxicokinetics and ecotoxicology of organic compounds. A shortcut approach to characterize the in vivo biotransformation rate constant (k_M ) with incomplete pathway or metabolite details was proposed. The value of k_M can be derived as 2tln1fPC(t)) , with f_PC (t) being the molar equivalent fraction of the parent compound (PC) at an early time t in both constant exposure and decay source chemical uptake scenarios. The approximation-based k_M values agreed well with k_M values derived from rigorous fitting or toxicokinetic modeling (n = 42, root mean square error = 0.30) with accuracy exceeding those of typical toxicokinetic or partitioning models. The method is accurate when sampling time is adequately resolved (i.e., t < ln(2)/k_M ) but will likely produce biased k_M values with improper time-averaging. The approximate equation yields consistent theoretical expectations for fast and slow biotransformation reactions and is fully compatible with standard bioaccumulation and toxicity testing protocols. The simplification strategy circumvents statistical complications and numerical issues inherent in regressing or modeling the toxicokinetics of multimetabolite systems and may be adapted to similar problems at other physiological scales or ecotoxicological contexts. The method can help advance interspecies comparison of chemical metabolism and support the development of in vitro-in vivo extrapolations and in silico models needed for building next-generation ecological and health risk-assessment practices. Environ Toxicol Chem 2022;41:148-158. © 2021 SETAC.

The added value of Bayesian inference for estimating biotransformation rates of organic contaminants in aquatic invertebrates

Toxicokinetic (TK) models refer to the process of contaminant bioaccumulation as a balance between rate of uptake from different sources (e.g., water or food), and rate of elimination via different processes such as excretion, growth and/or biotransformation. Biotransformation can considerably modify the fate of chemicals in an organism, especially their bioavailability, residence time, and toxicity. Invertebrate models generally neglect this process as they assume a low metabolic activity. However, some species such as Gammarus sp. amphipods are able to metabolize a vast range of organic compounds. Some recent TK models include biotransformation, but they prove limited for estimating related parameters by giving negative values and/or large uncertainties for biotransformation rate(s). Here we propose a generic TK model accounting for biotransformation using a Bayesian framework for simultaneously estimating the parameters. We illustrated the added value of our method by fitting this generic TK model to 22 published datasets of several benthic invertebrate species exposed to different chemicals. All parameters are estimated simultaneously for all datasets and showed narrow estimates. Furthermore, the median model predictions and their 95% credibility intervals showed that the model confidently fitted the data. In most cases the uncertainties around biotransformation rate(s) were reduced in comparison to the original studies. From a methodology standpoint, this paper reflects that Bayesian inference has real added value for simultaneously estimating all TK parameters for parent chemicals and their metabolite(s) based on all available data, while accounting for different types of data and the correlation between parameters. Bayesian inference was able to overcome the limits of previous methods, since no parameters were fixed and no irrelevant negative values were obtained. Moreover, the 95% credibility intervals around model predictions, which are core uncertainties for Environmental Risk Assessment, were easily acquired.

A Reduced Model for Bioconcentration and Biotransformation of Neutral Organic Compounds in Midge

A bioconcentration factor (BCF) database and a toxicokinetic model considering only biota-water partitioning and biotransformation were constructed for neutral organic chemicals in midge. The database contained quality-reviewed BCF and toxicokinetic data with variability constrained to within 0.5 to 1 log unit. Diverse conditions in exposure duration, flow set-up, substrate presence, temperature, and taxonomic classification did not translate into substantial variability in BCF, uptake rate constant (k₁ ), or depuration rate constant (k_T ), and no systematic bias was observed in BCFs derived in unlabeled versus radiolabeled studies. Substance-specific biotransformation rate constants k_M were derived by difference between the calculated biota-water partitioning coefficient (K_BW ) and experimental BCF for developing a midge biotransformation model. Experimental midge BCF was modeled as BCF = K_BW /(1 + k_M/ k₂ ) with log k_M (k_M in h^-1 ) = -0.37 log K_OW - 0.06T (in K) + 18.87 (root mean square error [RMSE] = 0.60), log k₁ (k₁ in L kg_wet.wt^-1  h^-1 ) = -0.0747 W (body weight in mg_wet.wt ) + 2.35 (RMSE = 0.48). The K_BW value was estimated using midge biochemical composition and established polyparameter linear free energy relationships, and the diffusive elimination rate constant (k₂ ) was computed as k₂  = k₁ /K_BW. The BCF model predicted >85% of BCFs that associated with neutral organic compounds (log K_OW  = 1.46 - 7.75) to within 1 log-unit error margin and had comparable accuracy similar to amphipod or fish models. A number of outliers and critical limitations of the k_M model were identified and examined, and they largely reflected the inherent limitation of difference-derived k_M , the lack of chemical diversity, and inadequate temperature variation in existing data. Future modeling efforts can benefit from more BCF and toxicokinetic observations of BCF on structurally diverse chemicals for model training, validation, and diagnosis. Environ Toxicol Chem 2021;40:57-71. © 2020 SETAC.

Tissue-specific bioaccumulation, metabolism and excretion of tris (2-ethylhexyl) phosphate (TEHP) in rare minnow (Gobiocyprisrarus)

Tris (2-ethylhexyl) phosphate (TEHP) is one of the most commonly used organophosphorus flame retardant (OPFR) analogues and is commonly detected in surface water and sediments. Limited information is available about the metabolic pathway or metabolite formation related to TEHP in fish. In this study, rare minnows (Gobiocyprisrarus) were exposed to TEHP in static water for 30 d to investigate the bioaccumulation and metabolite distribution in the fish muscle, liver, kidney, gill, GI-tract, ovary and testis. Based on the estimated k_up,_parent and k_d,parent values, the bioconcentration factors (BCF_parent) of TEHP in fish tissues were calculated in the order of kidney > ovary ≈ liver ≈ testis > gill ≈ GI-tract > muscle; this finding was consistent with the results of our previous study on other alkyl-substituted OPFRs. In addition, this study identified the metabolic profiles of TEHP in the liver. TEHP was oxidatively metabolized by the fish to a dealkylated metabolite (di 2-ethylhexyl phosphate; DEHP) and hydroxylated TEHP (OH-TEHP). OH-TEHP further underwent extensive phase II metabolism to yield glucuronic acid conjugates. DEHP was mainly distributed in rare minnow in the following order: liver > GI-tract > kidney ≫ other tissues. However, the metabolite showed lower accumulation potential in fish tissues than TEHP, with metabolite parent concentration factors (MPCFs) for DEHP of less than 0.1 in all the investigated tissues. The BCF_parent values of TEHP in various fish tissues were only 9.0 × 10^-3-7.2 × 10^-4 times its estimated tissue-water partition coefficient (K_tissue-water) values based on tissue lipid, protein and water contents, which indicated the significance of biotransformation in reducing the bioaccumulation potential of TEHP in fish. The toxicokinetic data in the present study help in understanding the tissue-specific bioaccumulation and metabolism pathways of TEHP in fish and highlight the importance of toxicology research on TEHP metabolites in aquatic organisms.

Multicompartmental Toxicokinetic Modeling of Discrete Dietary and Continuous Waterborne Uptake of Two Polycyclic Aromatic Hydrocarbons by Zebrafish Danio rerio

In the present study, we developed a multicompartmental toxicokinetic model for two polycyclic aromatic hydrocarbons (phenanthrene and anthracene) in their deuterated form (PAHs-d₁₀) in zebrafish considering continuous waterborne uptake and discrete dietary uptake. We quantified the bioconcentration, bioaccumulation, and depuration of these two PAHs-d₁₀ in zebrafish, and then estimated the kinetic parameters by fitting the model into the experimental data. The experimental and fitting results both showed that there was a peak concentration in each compartment of zebrafish after every dietary uptake, while the peak value depended on the ingestion amount of the PAH-d₁₀ and varied among different compartments. The PAH-d₁₀ amount in the blood reached 20-27% of the total amount bioaccumulated in zebrafish at steady-state, followed by skin (20-26%), and fillet (16-22%). The rank of PAH-d₁₀ steady-state concentrations in each compartment showed inconsistency with its lipid contents, indicating that the distribution of the PAHs-d₁₀ in zebrafish was not merely affected by the lipid content in each compartment, but also affected by their kinetics and biotransformation. This study suggests that discrete dietary uptake caused by intermittent food ingestion significantly affects the bioaccumulation of PAHs in fish. Further studies are needed to investigate such effect on other toxicants that are more resistant to biotransformation than PAHs in fish.

Occurrence, biomagnification and maternal transfer of legacy and emerging organophosphorus flame retardants and plasticizers in water snake from an e-waste site

Water snake and small common carp samples collected from a Chinese pond polluted with electronic waste (e-waste) were analyzed for organophosphorus flame retardants (PFRs), PFR metabolites, and plasticizers to investigate their occurrence, biomagnification, and maternal transfer in ovoviviparous species. Mean concentrations of total PFRs, PFR metabolites, and plasticizers were 2.2-16, 1.3-2.8 and 151-1320 ng/g wet weight (ww), respectively in analyzed organisms. Metabolites of PFRs were found in the same order of magnitude as or even higher than their parent compounds, indicating the importance of monitoring metabolites to evaluate the internal exposure of PFRs in organisms. Biomagnification factors (BMFs) were below 1 for all targeted chemicals and negatively correlated with metabolite/parent ratios (MPRs), suggesting a biodilution driven by metabolism. The lipid normalized concentrations were lower in eggs than in muscle for most of targeted chemicals. The maternal transfer potential was significantly and positively correlated with log K_OW (p < 0.05) when log K_OW was below 6.

Toxicokinetic patterns, metabolites formation and distribution in various tissues of the Chinese rare minnow (Gobiocypris rarus) exposed to tri(2‑butoxyethyl) phosphate (TBOEP) and tri-n-butyl phosphate (TNBP)

Alkylated organophosphate esters (alkyl-OPEs) are widely used and extensively detected in aquatic organisms. This work investigated the tissue-specific toxicokinetics of two common alkyl-OPEs, tri(2‑butoxyethyl) phosphate (TBOEP) and tri‑n‑butyl phosphate (TNBP) in Chinese rare minnow (Gobiocypris rarus) through a 50 day uptake and depuration experiment. The tissue-specific bioconcentration factor (BCF) values for the two alkyl-OPEs ranged from 1 to 30 L/kg wet weight (ww), with the kidney and ovary as the tissues with the highest accumulation. The tissue BCFs only exhibited a significant correlation with lipid contents only in storage tissues (i.e., muscle, brain, ovary and testis), indicating that lipids might not be the major contributor to tissue distribution of TBOEP and TNBP. However, the contribution of blood perfusion and active transport to tissue-specific OPE accumulation needs to be further investigated. Lower accumulation of metabolites than parent chemicals was observed, with metabolite parent concentration factors (MPCFs) <1. Di-alkyl phosphate (DAP), bis(2‑butoxyethyl) phosphate (BBOEP) and di(n-butyl) phosphate (DNBP) were the most abundantly formed metabolites of TBOEP and TNBP in various tissues, followed by the monohydroxylated OPEs (OH-OPEs). However, bis(2‑butoxyethyl) hydroxyethyl phosphate (BBOEHEP), was detected at much lower levels in the tissues. All the investigated metabolites showed high production rates (k_{prod,metabolites}) in the fish liver, followed by the GI tract and the kidney, indicating the importance of the hepatobiliary and urinary systems in eliminating the metabolites. Our study suggested that metabolism plays an important role in eliminating these two alkyl-OPEs in rare minnow and results in different tissue distribution mechanisms for metabolites and their compounds.

Bioconcentration model for non-ionic, polar, and ionizable organic compounds in amphipod

The present study presents a bioconcentration model for non-ionic, polar, and ionizable organic compounds in amphipod based on first-order kinetics. Uptake rate constant k₁ is modeled as logk1=10.81logKOW + 0.15 (root mean square error [RMSE] = 0.52). Biotransformation rate constant k_M is estimated using an existing polyparameter linear free energy relationship model. Respiratory elimination k₂ is calculated as modeled k₁ over theoretical biota-water partition coefficient K_biow considering the contributions of lipid, protein, carbohydrate, and water. With negligible contributions of growth and egestion over a typical amphipod bioconcentration experiment, the bioconcentration factor (BCF) is modeled as k₁ /(k_M + k₂ ) (RMSE = 0.68). The proposed model performs well for non-ionic organic compounds (log K_OW range = 3.3-7.62) within 1 log-unit error margin. Approximately 12% of the BCFs are underpredicted for polar and ionizable compounds. However, >50% of the estimated k₂ values are found to exceed the total depuration rate constants. Analyses suggest that these excessive k₂ values and underpredicted BCFs reflect underestimation in K_biow , which may be improved by incorporating exoskeleton as a relevant partitioning component and refining the membrane-water partitioning model. The immediate needs to build up high-quality experimental k_M values, explore the sorptive role of exoskeleton, and investigate the prevalence of k₂ overestimation in other bioconcentration models are also identified. The resulting BCF model can support, within its limitations, the ecotoxicological and risk assessment of emerging polar and ionizable organic contaminants in aquatic environments and advance the science of invertebrate bioaccumulation. Environ Toxicol Chem 2018;37:1378-1386. © 2018 SETAC.

Accumulation and distribution of organophosphate flame retardants (PFRs) and their di-alkyl phosphates (DAPs) metabolites in different freshwater fish from locations around Beijing, China

Organophosphate flame retardants (PFRs) can be rapidly metabolized in the body, and recent studies have shown that the di-alkyl phosphates (DAPs) are important metabolites. The accumulation and distribution of 8 PFRs and their 4 DAPs metabolites were first investigated in whole-body samples and various tissues of three freshwater fish species (topmouth gudgeon, crucian carp and loach) with different feeding habits from locations around Beijing, China. Concentrations of ΣPFRs in whole-body samples across all sampling locations ranged from 264.7 to 1973 ng g^-1 lipid weight (lw), while all the paired DAP metabolites were detected in the total range from 35.3 to 510 ng g^-1 lw. The calculated log bioconcentration factors (BCFs) of PFRs in whole fish were correlated with their log K_OW (P < 0.05). The metabolite/parent ratios (MPRs) of ΣDAPs were calculated and ranged from 0.10 to 1.12 in whole-fish of all species. The MPRs of BBOEP/TBOEP were the highest. With respect to their distribution in different tissues, both the parent PFRs and metabolites were found at relatively higher levels in the liver than in other tissues (muscle, intestine, kidney and ovary), which was markedly different from those observed in avian species in previous studies. The accumulation of PFRs and DAPs in various tissues was not significantly correlated with the lipid content. The highest PFRs level in the liver may be related to the active hepatic accumulation processes. Meanwhile, the MPRs for all 4 pairs were the highest in the kidney relative to the other tissues. To the best of our knowledge, this is first study of DAPs in wild animals, and our study may improve the understanding of the accumulation and metabolism of PFRs in the body.