How Biotransformation Influences Toxicokinetics of Azole Fungicides in the Aquatic Invertebrate Gammarus pulex

Metabolism of a fungicide propiconazole by Cunninghamella elegans ATCC36112

The metabolic breakdown of propiconazole by fungi was examined, and it was found that the microbial model (Cunninghamella elegans ATCC36112) efficiently degrades the triazole fungicide propiconazole through the action of cytochrome P450. This enzyme primarily facilitates the oxidation and hydrolysis processes involved in phase I metabolism. We observed major metabolites indicating hydroxylation/oxidation of propyl groups of propiconazole. Around 98% of propiconazole underwent degradation within a span of 3 days post-treatment, leading to the accumulation of five metabolites (M1-M5). The experiments started with a preliminary identification of propiconazole and its metabolites using GC-MS. The identified metabolites were then separated and identified by in-depth analysis using preparative UHPLC and MS/MS. The metabolites of propiconazole are M1 (CGA-118245), M2(CGA-118244), M3(CGA-136735), M4(GB-XLIII-42-1), and M5(SYN-542636). To further investigate the role of key enzymes in potential fungi, we treated the culture medium with piperonyl butoxide (PB) and methimazole (MZ), and then examined the kinetic responses of propiconazole and its metabolites. The results indicated a significant reduction in the metabolism rate of propiconazole in the medium treated with PB, while methimazole showed weaker inhibitory effects on the metabolism of propiconazole in the fungus C. elegans.

Fate of synthetic chemicals in the agronomic insect pest Spodoptera littoralis: experimental feeding-contact assay and toxicokinetic model

Insecticides prevent or reduce insect crop damage, maintaining crop quality and quantity. Physiological traits, such as an insect's feeding behavior, influence the way insecticides are absorbed and processed in the body (toxicokinetics), which can be exploited to improve species selectivity. To fully understand the uptake of insecticides, it is essential to study their total uptake and toxicokinetics independent of their toxic effects on insects. We studied the toxicokinetics (TK) of insecticidally inactive test compounds incorporating agro-like structural motifs in larvae of the Egyptian cotton leafworm (Spodoptera littoralis, Lepidoptera), and their distribution across all biological matrices, using laboratory experiments and modeling. We measured Spodoptera larval behavior and temporal changes of whole-body concentrations of test compounds during feeding on treated soybean leaf disks and throughout a subsequent depuration period. Differences in the distribution of the total quantities of compounds were found between the biological matrices leaf, larva, and feces. Rate constants for uptake and elimination of test compounds were derived by calibrating a toxicokinetic model to the whole-body concentrations. Uptake and elimination rate constants depended on the physicochemical properties of the test compounds. Increasing hydrophobicity increased the bioaccumulation potential of test compounds. Incomplete quantities in larval matrices indicated that some compounds may undergo biotransformation. As fecal excretion was a major elimination pathway, the variable time of release and number of feces pellets led to a high variability in the body burden. We provide quantitative models to predict the toxicokinetics and bioaccumulation potential of inactive insecticide analogs (parent compounds) in Spodoptera.

Immunoassay for the detection of cyproconazole in foods: From hapten synthesis to the establishment of a gold immunochromatographic assay

Cyproconazole (CPZ) is extensively used in agricultural production. However, its overuse can lead to high residue problems in crops. Existing detection methods are still dominated by instrumental methods and the development of rapid, sensitive field detection remains a challenge. In this study, we designed a novel hapten synthetic pathway and prepared a monoclonal antibody (mAb) that could specifically recognize CPZ with high sensitivity (half inhibition rate was 0.27 ng/mL). From this, a gold immunochromatographic assay (GICA) for the detection of CPZ was established by combining the mAb with gold nanoparticles, with limits of detection in rice, tomatoes and grapes of 0.02 mg/kg, 0.01 mg/kg and 0.05 mg/kg, respectively. The spiked recoveries ranged from 86.5 % to 115.1 %, and the results showed that the GICA was not significantly different from detection using LC-MS/MS. Therefore, we have successfully developed a GICA method for the reliable in situ, rapid and sensitive detection of CPZ.

Multiomics Sequencing and AlphaFold2 Analysis of the Stereoselective Behavior of Mefentrifluconazole for Bioactivity Improvement and Risk Reduction

As the first isopropanol chiral triazole fungicide, mefentrifluconazole has broad prospects for application. In this study, the stereoselective stability, bioactivity, fate, and biotoxicity were systematically investigated. Our results indicated that the stability of mefentrifluconazole enantiomers differed between environmental media, and they were stable in water and sediment in the dark. The bactericidal activity of R-mefentrifluconazole against the four target pathogens was 4.6-43 times higher than that of S-mefentrifluconazole. In the water-sediment system, S-mefentrifluconazole dissipated faster than R-mefentrifluconazole in water; however, its accumulation capacity was higher than that of R-mefentrifluconazole in sediment and zebrafish. S-Mefentrifluconazole induced more differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) in zebrafish than did R-mefentrifluconazole. Multiomics sequencing results showed that S-mefentrifluconazole enhanced the antioxidant, detoxification, immune, and metabolic functions of zebrafish by interacting with related proteins. Based on AlphaFold2 modeling and molecular docking, mefentrifluconazole enantiomers had different binding modes with key target proteins in pathogens and zebrafish, which may be the main reason for the stereoselective differences in bioactivity and biotoxicity. Based on its excellent bioactivity and low biotoxicity, the R-enantiomer can be developed to improve the bioactivity and reduce the risk of mefentrifluconazole.

Elucidating the spatial distribution of organic contaminants and their biotransformation products in amphipod tissue by MALDI- and DESI-MS-imaging

The application of mass spectrometry imaging (MSI) is a promising tool to analyze the spatial distribution of organic contaminants in organisms and thereby improve the understanding of toxicokinetic and toxicodynamic processes. MSI is a common method in medical research but has been rarely applied in environmental science. In the present study, the suitability of MSI to assess the spatial distribution of organic contaminants and their biotransformation products (BTPs) in the aquatic invertebrate key species Gammarus pulex was studied. Gammarids were exposed to a mixture of common organic contaminants (carbamazepine, citalopram, cyprodinil, efavirenz, fluopyram and terbutryn). The distribution of the parent compounds and their BTPs in the organisms was analyzed by two MSI methods (MALDI- and DESI-HRMSI) after cryo-sectioning, and by LC-HRMS/MS after dissection into different organ compartments. The spatial distribution of contaminats in gammarid tissue could be successfully analyzed by the different analytical methods. The intestinal system was identified as the main site of biotransformation, possibly due to the presence of biotransforming enzymes. LC-HRMS/MS was more sensitive and provided higher confidence in BTP identification due to chromatographic separation and MS/MS. DESI was found to be the more sensitive MSI method for the analyzed contaminants, whereas additional biomarkers were found using MALDI. The results demonstrate the suitability of MSI for investigations on the spatial distribution of accumulated organic contaminants. However, both MSI methods required high exposure concentrations. Further improvements of ionization methods would be needed to address environmentally relevant concentrations.

Comparative insight of pesticide transformations between river and wetland systems

The widespread use of pesticides threatens the environment and ecosystems. Despite the positive effects of plant protection products, pesticides also have unexpected negative effects on nontarget organisms. The microbial biodegradation of pesticides is one of the major pathways for reducing their risks at aquatic ecosystems. The objective of this study was to compare the biodegradability of pesticides in simulated wetland and river systems. Parallel experiments were conducted with 17 pesticides based on the OECD 309 guidelines. A comprehensive analytical method, such as target screening combined with suspect and non-target screening, was performed to evaluate the biodegradation via identification of transformation products (TPs) using LC-HRMS. As evidence of biodegradation, we identified 97 TPs for 15 pesticides. Metolachlor and dimethenamid had 23 and 16 TPs, respectively, including Phase II glutathione conjugates. The analysis of 16S rRNA sequences for microbials characterized operational taxonomic units. Rheinheimera and Flavobacterium, which have the potential for glutathione S-transferase, were dominant in wetland systems. Estimation of toxicity, biodegradability, and hydrophobicity using QSAR prediction indicated lower environmental risks of detected TPs. We conclude that the wetland system is more favorable for pesticide degradation and risk mitigation mainly attributed to the abundance and variety of the microbial communities.

Speed it up: How temperature drives toxicokinetics of organic contaminants in freshwater amphipods

The acceleration of global climate change draws increasing attention towards interactive effects of temperature and organic contaminants. Many studies reported a higher sensitivity of aquatic invertebrates towards contaminant exposure with increasing or fluctuating temperatures. The hypothesis of this study was that the higher sensitivity of invertebrates is associated with the changes of toxicokinetic processes that determine internal concentrations of contaminants and consequently toxic effects. Therefore, the influence of temperature on toxicokinetic processes and the underlying mechanisms were studied in two key amphipod species (Gammarus pulex and Hyalella azteca). Bioconcentration experiments were carried out at four different temperatures with a mixture of 12 exposure relevant polar organic contaminants. Tissue and medium samples were taken in regular intervals and analysed by online solid-phase extraction liquid chromatography high-resolution tandem mass spectrometry. Subsequently, toxicokinetic rates were modelled and analysed in dependence of the exposure temperature using the Arrhenius equation. An exponential relationship between toxicokinetic rates versus temperature was observed and could be well depicted by applying the Arrhenius equation. Due to a similar Arrhenius temperature of uptake and elimination rates, the bioconcentration factors of the contaminants were generally constant across the temperature range. Furthermore, the Arrhenius temperature of the toxicokinetic rates and respiration was mostly similar. However, in some cases (citalopram, cyprodinil), the bioconcentration factor appeared to be temperature dependent, which could potentially be explained by the influence of temperature on active uptake mechanisms or biotransformation. The observed temperature effects on toxicokinetics may be particularly relevant in non-equilibrated systems, such as exposure peaks in summer as exemplified by the exposure modelling of a field measured pesticide peak where the internal concentrations increased by up to fourfold along the temperature gradient. The results provide novel insights into the mechanisms of chemical uptake, biotransformation and elimination in different climate scenarios and can improve environmental risk assessment.

The prochloraz chronic exposure to Daphnia magna derived in biochemical alterations of F0 generation daphnids and malformed F1 progeny

In the present study, D. magna individuals were exposed to several sublethal prochloraz concentrations (87, 130, 170, 230 and 380 μg/L) for 21 days according to; the previous acute toxicity results. The fungicide effects on reproduction, survival, individual size, and growth population rate were evaluated after an exposure of 21 days, and no changes were observed compared to the control group. On the other hand, F1 generation neonates were collected and their external morphology evaluated; to estimate if the fungicide concentrations used induced effects during oogenesis and; embryogenesis processes. Neonates from parents which were previously exposed to 170 μg/L and higher concentrations were malformed since 16-d of exposure onward. All animals presented the same malformation: asymmetrical shell morphology and separated valves that did not cover the complete animal body regardless of the; concentration. The biochemical parameters tested in the broodstock were cholesterol, triglycerides, glucose and LDH activity. At the end of the chronic exposure experiment, cholesterol and triglycerides remained unaltered while glucose and the LDH enzyme levels increased significantly. The results of the present work showed a direct effect of; prochloraz on D. magna individual growth, along with mobilization of some; biochemical intermediate metabolism. A daphnid stress response as a result of the fungicide presence in the medium could be an explanation for the metabolic disorders. On the other hand, the F1 malformed neonates found in the present study suggested an effect of prochloraz among different daphnid generations and more studies would be necessary in this field.

Toxicokinetics of Two Oxathiapiprolin Enantiomers in Rats and Their Stereoselective Interaction with Oxysterol Binding Protein

Oxathiapiprolin is a chiral fungicide, and it can affect the metabolism of the cholesterol compounds by inhibiting oxysterol binding protein (OSBP) to exert its fungicidal effect. The application of oxathiapiprolin in agriculture is widespread, and its residue in the environment is a threat to both human and animal health. The two oxathiapiprolin enantiomers differ in their fungicidal activity, biotoxicity, and degradation by environmental forces. However, their biotoxicity has not been reported in animals. The toxicokinetics of a pesticide should be a crucial component for the evaluation of its toxicity in vivo. In this study, we investigated the absorption, bioavailability, tissue distribution, and excretion of the two oxathiapiprolin enantiomers in rats to verify their toxicokinetic process in animals. An ultrahigh-performance liquid chromatography triple quadrupole tandem mass spectrometry (UHPLC-QQQ/MS) method was established to quantify the two oxathiapiprolin enantiomers in vivo. The two oxathiapiprolin enantiomers were found to have approximately the same absorption rate and bioavailability, and both were excreted mainly in the feces. The half-life of R-(-)-oxathiapiprolin was nearly twice that of S-(+)-oxathiapiprolin. R-(-)-oxathiapiprolin also had greater distribution than S-(+)-oxathiapiprolin in the liver, lungs, heart, spleen, kidneys, stomach, large intestine, small intestine, brain, and pancreas, supporting the notion that R-(-)-oxathiapiprolin could better bind with OSBP. The stereoselectivity of S-(+)-oxathiapiprolin in these tissues may be responsible for it being readily metabolized in vivo. The molecular docking technique was subsequently used to verify the more superior binding between R-(-)-oxathiapiprolin and OSBP compared with the binding between S-(+)-oxathiapiprolin and OSBP. The findings of this study could provide more reliable data for determining the toxicokinetics of a single enantiomer of oxathiapiprolin in animals, thereby providing some theoretical basis for its subsequent toxicological study.

Environmental behavior of the chiral fungicide epoxiconazole in earthworm-soil system: Enantioselective enrichment, degradation kinetics, chiral metabolite identification, and biotransformation mechanism

The environmental impact of the chiral fungicide epoxiconazole and its chiral transformation products (TPs) on non-target organisms and the environment has become a significant concern due to its widespread use in agricultural practice. Enantioselectivity studies of parent contaminants cannot adequately assess the complexity of its chiral TPs in the environment. This study aimed to investigate the environmental behavior of epoxiconazole in an earthworm-soil system. 2S,3R-(-)-epoxiconazole was preferentially enriched in earthworms during the accumulation phase (p < 0.05), but no enantioselectivity was observed during the elimination phase. One methoxylated and four hydroxylated chiral TPs were identified in soil, earthworm, and excrement. The epoxy ring hydroxylated TP and methoxylated TP of epoxiconazole were discovered for the first time in the environment. The chemically specific enantioselectivity with enantiomer fraction (EF) > 0.8 was observed for the TPs in different matrices. The CYP450 monooxygenase of earthworm was significant activated. In vitro enzyme metabolism experiments (earthworm microsomes and recombinant CYP450 enzymes CYP2A6, CYP 2C9, and CYP 3A4) were carried out to further explain the biotransformation mechanism of epoxiconazole in earthworm. This study provides new evidence of enantiomeric biotransformation of chiral fungicide epoxiconazole in the earthworm-soil system and could provide valuable insights into their environmental risk assessment.

Systematic Underestimation of Pesticide Burden for Invertebrates under Field Conditions: Comparing the Influence of Dietary Uptake and Aquatic Exposure Dynamics

Pesticides used in agriculture can end up in nearby streams and can have a negative impact on nontarget organisms such as aquatic invertebrates. During registration, bioaccumulation potential is often investigated using laboratory tests only. Recent studies showed that the magnitude of bioaccumulation in the field substantially differs from laboratory conditions. To investigate this discrepancy, we conducted a field bioaccumulation study in a stream known to receive pollutant loadings from agriculture. Our work incorporates measurements of stream pesticide concentrations at high temporal resolution (every 20 min), as well as sediment, leaves, and caged gammarid analyses (every 2-24 h) over several weeks. Of 49 investigated pesticides, 14 were detected in gammarids with highly variable concentrations of up to 140 ± 28 ng/g_ww. Toxicokinetic modeling using laboratory-derived uptake and depuration rate constants for azoxystrobin, cyprodinil, and fluopyram showed that despite the highly resolved water concentrations measured, the pesticide burden on gammarids remains underestimated by a factor of 1.9 ± 0.1 to 31 ± 3.0, with the highest underestimations occurring after rain events. Including dietary uptake from polluted detritus leaves and sediment in the model explained this underestimation only to a minor proportion. However, suspended solids analyzed during rain events had high pesticide concentrations, and uptake from them could partially explain the underestimation after rain events. Additional comparison between the measured and modeled data showed that the pesticide depuration in gammarids is slower in the field. This observation suggests that several unknown mechanisms may play a role, including lowered enzyme expression and mixture effects. Thus, it is important to conduct such retrospective risk assessments based on field investigations and adapt the registration accordingly.

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.

Monitoring of pharmaceuticals in aquatic biota (Procambarus clarkii) of the Doñana National Park (Spain)

In this work the presence of different pharmaceuticals at Doñana National Park (Spain) and their main entry sources (input source or entry points) have been stated over the 2011-2016 years period. Twenty-three selected pharmaceuticals (corresponding to eight therapeutic families) were evaluated in crayfish and water samples from Doñana National Park (Spain) (six sampling points selected in order to cover different possible pollution sources into and surrounding the Park). The multiresidue determination was carried out using enzymatic-microwave assisted extraction prior to high performance liquid chromatography mass spectrometry detection. Sulphonamides (sulfadiazine, sulfamerazine, sulfamethazine, and sulfamethoxazole); trimethoprim, an antibiotic that is frequently co-administered with sulfamethoxazole; amphenicols (chloramphenicol, florfenicol and thiamphenicol); fluoroquinolones (ciprofloxacin, enrofloxacin, flumequine, danofloxacin, gatifloxacin, norfloxacin, marbofloxacin and grepafloxacin); penicillins (amoxicillin); tetracyclines (chlortetracycline and oxytetracycline); non-steroidal anti-inflammatory drugs (salicylic acid and ibuprofen); beta-blocker drugs (atenolol); and antiepileptics (carbamazepine) were analysed. Ciprofloxacin, ibuprofen, salicylic acid, flumequine, and carbamazepine were detected and/or quantified at some of the selected sampling points. A clear ecotoxicological risk to the ecosystem was demonstrated from the occurrence of ciprofloxacin in samples obtained after the punctual and massive presence of people inside the Park. Furthermore, flumequine and carbamazepine have been detected in Procambarus clarkii specimens in concentrations around 30 ng g^-1 and 14 ng g^-1, respectively, and their occurrence in the specimens could indicate the persistence of the discharge sources. The main source of pharmaceuticals into the Park might be the livestock farming activities, and the influence of urban wastewaters from surrounding villages does not seem to be very important.

Biotransformation Capacity of Zebrafish (Danio rerio) Early Life Stages: Functionality of the Mercapturic Acid Pathway

Zebrafish (Danio rerio) early life stages offer a versatile model system to study the efficacy and safety of drugs or other chemicals with regard to human and environmental health. This is because, aside from the well-characterized genome of zebrafish and the availability of a broad range of experimental and computational research tools, they are exceptionally well suited for high-throughput approaches. Yet, one important pharmacokinetic aspect is thus far only poorly understood in zebrafish embryo and early larvae: their biotransformation capacity. Especially, biotransformation of electrophilic compounds is a critical pathway because they easily react with nucleophile molecules, such as DNA or proteins, potentially inducing adverse health effects. To combat such adverse effects, conjugation reactions with glutathione and further processing within the mercapturic acid pathway have evolved. We here explore the functionality of this pathway in zebrafish early life stages using a reference substrate (1-chloro-2,4-dinitrobenzene, CDNB). With this work, we show that zebrafish embryos can biotransform CDNB to the respective glutathione conjugate as early as 4 h postfertilization. At all examined life stages, the glutathione conjugate is further biotransformed to the last metabolite of the mercapturic acid pathway, the mercapturate, which is slowly excreted. Being able to biotransform electrophiles within the mercapturic acid pathway shows that zebrafish early life stages possess the potential to process xenobiotic compounds through glutathione conjugation and the formation of mercapturates. The presence of this chemical biotransformation and clearance route in zebrafish early life stages supports the application of this model in toxicology and chemical hazard assessment.

Grandmother's pesticide exposure revealed bi-generational effects in Daphnia magna

Man-made chemicals are a significant contributor to the ongoing deterioration of numerous ecosystems. Currently, risk assessment of these chemicals is based on observations in a single generation of animals, despite potential adverse intergenerational effects. Here, we investigate the effect of the fungicide prochloraz across three generations of Daphnia magna. We studied both the effects of continuous exposure over all generations and the effects of first-generation (F0) exposure on two subsequent generations. Effects at different levels of biological organization from genome-wide gene expression, whole organism metabolite levels, CYP enzyme activity and key phenotypic effects, such as reproduction, were monitored. Acclimation to prochloraz was found after continuous exposure. Following F0-exposure, embryonically exposed F1-offspring showed no significant effects. However, in the potentially germline exposed F2 animals, several parameters differed significantly from controls. A direct association between these F2 effects and the toxic mode of action of prochloraz was found, showing that chemicals can be harmful not only to the directly exposed generation, but also to prenatally exposed generations and in that way effects may even appear to skip a generation. This implies that current risk assessment practices are neglecting an important aspect of toxicity, such as delayed effects across generations due to a time gap between chemical exposure and emergence of effects.

Pesticides: formulants, distribution pathways and effects on human health - a review

Pesticides are commonly used in agriculture to enhance crop production and control pests. Therefore, pesticide residues can persist in the environment and agricultural crops. Although modern formulations are relatively safe to non-target species, numerous theoretical and experimental data demonstrate that pesticide residues can produce long-term negative effects on the health of humans and animals and stability of ecosystems. Of particular interest are molecular mechanisms that mediate the start of a cascade of adverse effects. This is a review of the latest literature data on the effects and consequences of contamination of agricultural crops by pesticide residues. In addition, we address the issue of implicit risks associated with pesticide formulations. The effects of pesticides are considered in the context of the Adverse Outcome Pathway concept.

Biochemical and molecular responses of the Mediterranean mussel (Mytilus galloprovincialis) to short-term exposure to three commonly prescribed drugs

Pharmaceuticals represent a group of emerging contaminants. The short-term effect (3 and 7 days) of warfarin (1 and 10 mg L^-1), dexamethasone (0.392 and 3.92 mg L^-1) and imidazole (0.013 and 0.13 mg L^-1) exposure was evaluated on mussels (Mytilus galloprovincialis). Total antioxidant status, glutathione reductase, glutathione peroxidase (GPx) and superoxide dismutase enzyme activities, and the expression of genes involved in the xenobiotic response (ATP binding cassette subfamily B member 1 (abcb1) and several nuclear receptor family J (nr1j) isoforms), were evaluated. All nr1j isoforms are suggested to be the xenobiotic receptor orthologs of the NR1I family. All drugs increased GPx activity and altered the expression of particular nr1j isoforms. Dexamethasone exposure also decreased abcb1 expression. These findings raised some concerns regarding the release of these pharmaceuticals into the aquatic environment. Thus, further studies might be needed to perform an accurate environmental risk assessment of these 3 poorly studied drugs.

Biodegradation of typical azole fungicides in activated sludge under aerobic conditions

Widespread use of azole fungicides and low removal efficiency in wastewater treatment plants (WWTPs) have led to the elevated concentration of azole fungicides in receiving environment. However, there was limited research about the removal mechanism of azole fungicides in the biological treatment of WWTPs. Imidazole fungicide climbazole and triazole fungicide fluconazole were selected to investigate the biodegradation mechanism of azole fungicides in activated sludge under aerobic conditions. Climbazole was found to be adsorbed to solid sludge and resulted in quick biodegradation. The degradation of climbazole in the aerobic activated sludge system was fitted well by the first-order kinetic model with a half-life of 5.3 days, while fluconazole tended to stay in liquid and had only about 30% of loss within 77 days incubation. Ten biotransformation products of climbazole were identified by high resolution mass spectrometry using suspect and non-target screening method. But no biodegradation products of fluconazole were identified due to its limited removal. The possible biodegradation pathways for climbazole were proposed based on the products identification and pathway prediction system, and involves oxidative dehalogenation, side chain oxidation and azole ring loss. The findings from this study suggest that it should be a concern for the persistence of fluconazole in the environment.

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.

Emerging pharmaceuticals and industrial chemicals in Nakdong River, Korea: Identification, quantitative monitoring, and prioritization

The extensive development and use of new anthropogenic chemicals have inevitably led to their presence in aquatic environments. Surface waters affected by sewage effluents have been exposed to these new substances. In the present study, the occurrence of anthropogenic substances, including pharmaceuticals and industrial chemicals, was investigated in one of the major rivers in Korea, the Nakdong River. Furthermore, seasonal variations in their content were determined via annual monitoring. Through the suspect and non-target screening (SNTS) technique, 58 substances were newly identified in the river and integrated in the quantitative monitoring practice. The results revealed that niflumic acid and melamine exhibited the highest median concentrations, i.e., 320 ng/L and 11,000 ng/L, respectively. The results associated with seasonal change revealed that the concentration of a considerable number of substances increased in winter when the flow rate was low. Conversely, some substances exhibited high concentrations in summer (e.g., polyethylene glycol) and spring (e.g., niflumic acid). This was attributed to the seasonal changes in the consumption, prescriptions, or the application of alternative substances. These changes were also reflected by the risk quotient (RQ) values calculated from the concentration and toxicity values. Pharmaceuticals such as telmisartan and carbamazepine and industrial chemicals such as organophosphorus flame retardants (OPFRs) and melamine accounted for approximately 90% of the total RQ. Major substances prioritized using the production of the RQ value and the detection frequency included OPFRs and telmisartan. It is recommended that these results be reflected in future water quality monitoring plans.

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.

Identification of transformation products to characterize the ability of a natural wetland to degrade synthetic organic pollutants

Natural wetlands have been recognized as a natural reactor for degradation and elimination of environmental pollutants. The Upo Wetland, the largest inland wetland in Korea, is mainly surrounded by agricultural lands and it is susceptible to contamination from excess nutrient loads and synthetic organic contaminants (SOCs) (e.g., pesticides). The aim of this study was to identify major SOCs in the wetland and evaluate their degradation. We used high resolution mass spectrometry (HRMS) with a two-step analysis approach (i.e., 1^st analysis for target measurement along with suspect and non-target screening (SNTS) and 2^nd analysis for complimentary suspect screening) to identify and quantify the transformation products (TPs) of the identified parent SOCs. Quantitative analysis of 30 targets, mainly including pesticides, showed that fungicides were the major SOCs detected in the wetland, accounting for about 50% of the composition ratio of the total SOCs quantified. Orysastrobin occurred at the highest mean concentration (>700 ng/L), followed by two other fungicides, carbendazim and tricyclazole. The first analysis (SNTS) tentatively identified 39 TPs (30 by suspect, 9 by non-target screening) of 14 parent pesticides. Additionally, the second analysis (complimentary suspect screening) identified 9 more TPs. Among the 48 total TPs identified, 7 were confirmed with reference standards. The identification of the remaining TPs had a high confidence level (e.g., level 2 or 3). Regarding transport though the wetland, most TPs showed greater peak area ratios (i.e., the relative portion of chromatographic area of the TPs to the parent compound) at the outlet point of the wetland compared to the inlet point. The risk quotient, which was calculated using the concentrations of parent compounds, decreased toward the outlet, demonstrating the degradation capacity of the wetland. The estimates for biodegradability, hydrophobicity, and toxicity by an in-silico quantitative structure-activity relationship (QSAR) model indicated a lower half-life, lower logD_OW, and greater effect concentration for most TPs compared to the parent compounds. Based on these results, we conclude that natural wetlands play a role as an eco-friendly reactor for degrading SOCs to form numerous TPs that are lower risk than the parent compounds.

Retrospective screening of high-resolution mass spectrometry archived digital samples can improve environmental risk assessment of emerging contaminants: A case study on antifungal azoles

Environmental risk assessment associated with aquatic and terrestrial contamination is mostly based on predicted or measured environmental concentrations of a limited list of chemicals in a restricted number of environmental compartments. High resolution mass spectrometry (HRMS) can provide a more comprehensive picture of exposure to harmful chemicals, particularly through the retrospective analysis of digitally stored HRMS data. Using this methodology, our study characterized the contamination of various environmental compartments including 154 surface water, 46 urban effluent, 67 sediment, 15 soil, 34 groundwater, 24 biofilm, 41 gammarid and 49 fish samples at 95 sites widely distributed over the Swiss Plateau. As a proof-of-concept, we focused our investigation on antifungal azoles, a class of chemicals of emerging concern due to their endocrine disrupting effects on aquatic organisms and humans. Our results demonstrated the occurrence of antifungal azoles and some of their (bio)transformation products in all the analyzed compartments (0.1-100 ng/L or ng/g d.w.). Comparison of actual and predicted concentrations showed the partial suitability of level 1 fugacity modelling in predicting the exposure to azoles. Risk quotient calculations additionally revealed risk of exposure especially if some of the investigated rivers and streams are used for drinking water production. The case study clearly shows that the retrospective analysis of HRMS/MS data can improve the current knowledge on exposure and the related risks to chemicals of emerging concern and can be effectively employed in the future for such purposes.

Linking Morphology, Toxicokinetic, and Toxicodynamic Traits of Aquatic Invertebrates to Pyrethroid Sensitivity

Pyrethroid insecticides are known to be highly toxic to most aquatic nontarget organisms, but little is known about the mechanisms causing some species to be highly sensitive while others are hardly affected by the pyrethroids. The aim of the present study was to measure the sensitivity (EC₅₀-values) of 10 aquatic invertebrates toward a 24 h pulse of the pyrethroid cypermethrin and subsequently test if the difference in sensitivity could be explained by measured morphological and physiological traits and modeled toxicokinetic (TK) and toxicodynamic (TD) parameters. Large differences were observed for the measured uptake and elimination kinetics, with bioconcentration factors (BCFs) ranging from 53 to 2337 at the end of the exposure. Similarly, large differences were observed for the TDs, and EC₅₀-values after 168 h varied 120-fold. Modeling the whole organism cypermethrin concentrations indicated compartmentation into a sorbed fraction and two internal fractions: a bioavailable and non-bioavailable internal fraction. Strong correlations between surface/volume area and the TK parameters (sorption and uptake rate constants and the resulting BCF) were found, but none of the TK parameters correlated with sensitivity. The only parameter consistently correlating with sensitivity across all species was the killing rate constant of the GUTS-RED-SD model (the reduced general unified threshold models of survival assuming stochastic death), indicating that sensitivity toward cypermethrin is more related to the TD parameters than to TK parameters.

Characterizing biotransformation products and pathways of the flame retardant triphenyl phosphate in Daphnia magna using non-target screening

Triphenyl phosphate (TPHP), one of the organophosphate flame retardants, has been widely used in manufacturing, thereby causing a gradual increase in TPHP concentrations in aquatic environments. However, the information on the biotransformation mechanism of TPHP in invertebrates is lacking. The study identified the biotransformation products of TPHP in Daphnia magna, which showed particularly high toxicity in aquatic organisms, and determined the rates of depuration. Daphnia magna, a standard species for toxicity studies, was exposed to triphenyl phosphate and transferred to the pure medium. The biotransformation products of TPHP and its depuration rates were determined by liquid chromatography-high resolution mass spectrometry. Nine biotransformation products (five in the positive mode and four in the negative mode) of triphenyl phosphate were identified in D. magna. Based on the depuration ratio, the major biotransformation mechanism is estimated to be cysteine conjugation and sulfation. Certain biotransformation products (diphenyl phosphate, hydroxylated triphenyl phosphate, and thiol triphenyl phosphate) might induce toxicity in biota. The results could be used to predict main biotransformation processes and toxic products of organophosphate flame retardants in aquatic invertebrates.

Coupling River Concentration Simulations with a Toxicokinetic Model Effectively Predicts the Internal Concentrations of Wastewater-Derived Micropollutants in Field Gammarids

Although the exposure assessment of wastewater-derived micropollutants via chemical, bioanalytical, and modeling methods in environmental compartments is becoming more frequent, the whole-body burden (i.e., internal concentrations) in nontarget organisms is rarely assessed. An understanding of the internal concentration fluctuation is especially important when exploring the mechanistic linkage between exposure and effects. In this study, we coupled a simple river model with a first-order toxicokinetic (TK) model to predict the concentrations of wastewater-derived micropollutants in freshwater invertebrates (Gammarus spp.). We applied Monte Carlo simulations and conducted laboratory experiments to account for the uncertain input data and the lack of uptake/depuration rate constants required for the TK model. The internal concentrations in field gammarids were predicted well, and the estimates varied only by a factor of 0.1-1.9. Fast equilibrium may also be assumed such that bioconcentration factors (BCFs) are used together with the daily river dilution patterns to predict internal concentrations. While this assumption is suitable for compounds observed in our experiment to reach the steady state within 48 h in gammarids, the model overpredicted the concentrations of substances that reach this condition after longer periods. Nevertheless, this approach provides conservative estimates and simplifies the coupling of models as BCFs are slightly more accessible than the rate constants. However, if one is interested in a more detailed exposure information (e.g., peak concentration and the whole-body burden recovery after a spill), then the nonsteady-state formulation should be employed.

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.

Enantioselective degradation and bioaccumulation of sediment-associated fipronil in Lumbriculus variegatus: Toxicokinetic analysis

Enantioselective degradation and biotransformation are critical processes affecting the bioaccumulation and toxicity of chiral pesticides in the environment. In the present study, enantioselective uptake, biotransformation and elimination of a current use pesticide, fipronil in a benthic invertebrate, Lumbriculus variegatus were assessed using a sediment bioaccumulation test. Toxicokinetic models were constructed to quantitatively describe kinetic processes of fipronil enantiomers. The degradation of fipronil in sediment significantly affected chemical uptake, thus degradation kinetic model was incorporated into toxicokinetic modeling. It was shown that S-(+)-fipronil degraded faster than R-(-)-fipronil in sediment, with dissipation rate constants being 0.090 ± 0.008 and 0.023 ± 0.006 1/d, respectively. As a result, R-(-)-enantiomer preferentially accumulated in sediment over time. Similarly, higher concentrations of R-(-)-fipronil were detected in L. variegatus compared with S-(+)-fipronil. Toxicokinetic modeling showed R-(-)-fipronil had larger uptake and elimination rate coefficients and apparent maximum reaction rate, but a smaller apparent half-saturation constant than S-(+)-fipronil. Preferential uptake of R-(-)-fipronil from sediment to L. variegatus was the main reason for greater R-(-)-fipronil concentrations in organism. Biotransformation of fipronil in L. variegatus was also enantioselective, yet it played fewer roles on enantioselective bioaccumulation than uptake. Overall, our findings highlight the importance of selective degradation, uptake and biotransformation of sediment-associated fipronil on its enantioselective bioaccumulation in benthic invertebrates, which helps to improve the accuracy for assessing aquatic toxicity of the chiral pesticide. CAPSULE: Enantioselective bioaccumulation of sediment-associated fipronil in Lumbriculus variegatus was quantitatively explained by selective degradation, uptake, biotransformation and elimination parameters using a combination of degradation and toxicokinetic modeling.

The synergistic potential of azole fungicides does not directly correlate to the inhibition of cytochrome P450 activity in aquatic invertebrates

The ability of azole fungicides to inhibit cytochrome P450 dependent metabolism is proposed to be the main mechanism for their synergizing effect on pyrethroid insecticide toxicity in aquatic invertebrates. This study investigates the correlation between inhibition strength and synergistic potential of azole fungicides in the crustacean Daphnia magna and the insect larvae Chironomus riparius. Inhibition strength was measured in vivo toward the cytochrome P450 catalysed conversion of 7-ethoxycoumarin to 7-hydroxycoumarin (ECOD). Synergistic potentials were determined as the ratio between predicted and observed toxicity of mixtures based on the model of concentration addition (CA) and independent action (IA). Azoles (n = 9-11) enhanced the toxicity of α-cypermethrin in D. magna (Synergy ratios CA: 0.8 - 16; IA: 1.1 - 22) and inhibited cytochrome P450 activity by different degrees (IC50: 0.0023 - 36 μM for D. magna and 0.08 - 24 μM for C. riparius). Inhibition strengths were strongly correlated in the two organisms (r: 0.937 p: 0.019 for triazoles and r: 0.903 p: 0.097 for imidazoles). Lipophilicity governed the inhibition strength of triazoles in both species (r > 0.9, p < 0.05). No correlation was observed between inhibition strengths and synergistic potentials. Several reasons for the apparent lack of correlation were discussed.

Bioaccumulation, Biotransformation, and Synergistic Effects of Binary Fungicide Mixtures in Hyalella azteca and Gammarus pulex: How Different/Similar are the Two Species?

Aquatic organisms are consistently exposed to a mixture of micropollutants that can bioaccumulate, undergo biotransformation, and may exert mixture effects. However, little is known on the underlying mechanisms and species-specificity. Herein we investigated bioaccumulation, biotransformation and synergistic effects of azole (i.e., prochloraz) and strobilurin (i.e., azoxystrobin) fungicides in the two aquatic invertebrate species, Hyalella azteca and Gammarus pulex. Bioaccumulation of azoxystrobin was similar, whereas bioaccumulation of prochloraz was slightly different in the two species but was still significantly below the REACH criteria for bioaccumulative substances. Similar biotransformation patterns were observed in both species, and only a few unique biotransformation reactions were detected in H. azteca such as malonyl-glucose and taurine conjugation. Toxicokinetic modeling additionally indicated that biotransformation is a more important elimination pathway in H. azteca. In mixtures, no-observed-adverse-effect levels of prochloraz decreased the LC₅₀s of azoxystrobin in both species which correlated well with increased internal azoxystrobin concentrations. This synergistic effect is partly due to the inhibition of cytochrome P450 monooxygenases by prochloraz which subsequently triggered the reduced biotransformation of azoxystrobin (lower by five folds in H. azteca). The largely similar responses in both species suggest that the easier-to-cultivate H. azteca is a promising representative of invertebrates for toxicity testing.

Internal Concentrations in Gammarids Reveal Increased Risk of Organic Micropollutants in Wastewater-Impacted Streams

Internal concentrations link external exposure to the potential effect, as they reflect what the organisms actually take up and experience physiologically. In this study, we investigated whether frequently detected risk-driving substances in water were found in the exposed organisms and if they are classified the same based on the whole body internal concentrations. Field gammarids were collected upstream and downstream of ten wastewater treatment plants in mixed land use catchments. The sampling was conducted in autumn and winter, during low flow conditions when diffuse agricultural input was reduced. The field study was complemented with laboratory and flume experiments to determine the bioaccumulation potentials of selected substances. For 32 substances, apparent bioaccumulation factors in gammarids were determined for the first time. With a sensitive multiresidue method based on online-solid phase extraction followed by liquid chromatography coupled to high resolution mass spectrometry, we detected 63 (semi-) polar organic substances in the field gammarids, showing higher concentrations downstream than upstream. Interestingly, neonicotinoids, which are particularly toxic toward invertebrates, were frequently detected and were further determined as major contributors to the toxic pressure based on the toxic unit approach integrating internal concentration and toxic potency. The total toxic pressure based on internal concentrations was substantially higher compared to when external concentrations were used. Thus, internal concentrations may add more value to the current environmental risk assessment that is typically based solely on external exposure.

Stable Isotope Labeling-Assisted Metabolite Probing for Emerging Contaminants in Plants

Biotransformation is a notable modulator of the fate, bioaccumulation, and toxicity of contaminants in the environment. However, it is often formidable to identify unknown biotransformation products in the absence of reference standards, and this analytical challenge is particularly true for contaminants of emerging concern (CECs) that are mostly polar molecules without characteristic structures (e.g., Cl and Br) and in complex matrices such as plants. In this study, using the fibrate drug gemfibrozil as a model CEC and Arabidopsis thaliana as a model plant, we developed and demonstrated a novel analytical framework coupling deuterium stable isotope labeling with high-resolution mass spectrometry (SILAMS) in identifying plant biotransformation products. When exposed in A. thaliana cells, gemfibrozil was quickly taken up into the cells and extensively metabolized. The use of nonlabeled and deuterated gemfibrozil at a 3:1 ratio created unique diagnostic patterns in mass spectra, enabling the identification of 11 novel phase II amino acid/peptide conjugates. Similarity in mass fragmentation patterns and chromatographic behaviors was then employed to establish the probable structures. Two major metabolites were further confirmed as glutamate and glutamine conjugates using authentic standards. Most of the identified conjugates were also detected in the whole A. thaliana plant. Therefore, SILAMS offers unique advantages by excluding false matrix positives and helping discern unknown metabolites, including polar conjugates with endogenous biomolecules, with a high degree of confidence. This novel framework may be readily applied to other CECs for high-throughput metabolite screening in plants to improve our understanding of their food safety and human health risks and potential deleterious effects on other species living on plants.

A review of the pharmaceutical exposome in aquatic fauna

Pharmaceuticals have been considered 'contaminants of emerging concern' for more than 20 years. In that time, many laboratory studies have sought to identify hazard and assess risk in the aquatic environment, whilst field studies have searched for targeted candidates and occurrence trends using advanced analytical techniques. However, a lack of a systematic approach to the detection and quantification of pharmaceuticals has provided a fragmented literature of serendipitous approaches. Evaluation of the extent of the risk for the plethora of human and veterinary pharmaceuticals available requires the reliable measurement of trace levels of contaminants across different environmental compartments (water, sediment, biota - of which biota has been largely neglected). The focus on pharmaceutical concentrations in surface waters and other exposure media have therefore limited both the characterisation of the exposome in aquatic wildlife and the understanding of cause and effect relationships. Here, we compile the current analytical approaches and available occurrence and accumulation data in biota to review the current state of research in the field. Our analysis provides evidence in support of the 'Matthew Effect' and raises critical questions about the use of targeted analyte lists for biomonitoring. We provide six recommendations to stimulate and improve future research avenues.

Screening of Pesticide and Biocide Patterns As Risk Drivers in Sediments of Major European River Mouths: Ubiquitous or River Basin-Specific Contamination?

Pesticides and biocides (PaB) are ubiquitously present in aquatic ecosystems due to their widespread application and have been detected in rivers at concentrations that may cause distress to aquatic life. Many of these compounds accumulate in sediments acting as long-term source for aquatic ecosystems. However, data on sediment contamination with current-use PaB in Europe are scarce. Thus, in this study, we elucidated PaB patterns and associated risks in sediments of seven major European rivers focusing on their last stretch as an integrative sink of particles transported by these rivers. Sediments were extracted with pressurized liquid extraction (PLE) using a broad-spectrum method recovering many compound classes with a wide range of physicochemical properties. Altogether 126 compounds were analyzed and 81 of them were detected with LC-HRMS and GC-NCI-MS/MS at least in one of the sediments. The highest number of compounds was detected (59) in River Elbe sediments close to Cuxhaven with outstanding concentrations ranging from 0.8 to 1691 mg/g organic carbon. Multivariate analysis identified a cluster with 3 ubiquitous compounds (cyhalothrin, carbendazim, fenpropimorph) and three clusters of chemicals with higher variability within and between rivers. Risk assessment indicates an acute toxic risk to benthic crustaceans at all investigated sites with the pyrethroids tefluthrin and cyfluthrin together with the fungicide carbendazim as the main drivers. Risks to algae were driven at most sites almost exclusively by photosynthesis inhibitors with estuary-specific herbicide mixtures, while in the rivers Po and Gironde cell division inhibitors played an important role at some sites. Mixtures of specific concern have been defined and suggested for integration in future monitoring programs.

Binding Specificity Determines the Cytochrome P450 3A4 Mediated Enantioselective Metabolism of Metconazole

Cytochrome P450 3A4 (CYP3A4) is a promiscuous enzyme, mediating the biotransformations of ∼50% of clinically used drugs, many of which are chiral molecules. Probing the interactions between CYP3A4 and chiral chemicals is thus essential for the elucidation of molecular mechanisms of enantioselective metabolism. We developed a stepwise-restrained-molecular-dynamics (MD) method to model human CYP3A4 in a complex with cis-metconazole (MEZ) isomers and performed conventional MD simulations with a total simulation time of 2.2 μs to probe the molecular interactions. Our current study, which employs a combined experimental and theoretical approach, reports for the first time on the distinct conformational changes of CYP3A4 that are induced by the enantioselective binding of cis-MEZ enantiomers. CYP3A4 preferably metabolizes cis-RS MEZ over the cis-SR isomer, with the resultant enantiomer fraction for cis-MEZ increasing rapidly from 0.5 to 0.82. cis-RS MEZ adopts a more extended structure in the active pocket with its Cl atom exposed to the solvent, whereas cis-SR MEZ sits within the hydrophobic core of the active pocket. Free-energy-perturbation calculations indicate that unfavorable van der Waals interactions between the cis-MEZ isomers and the CYP3A4 binding pocket predominantly contribute to their binding-affinity differences. These results demonstrate that binding specificity determines the cytochrome P450 3A4 mediated enantioselective metabolism of cis-MEZ.

Can Toxicokinetic and Toxicodynamic Modeling Be Used to Understand and Predict Synergistic Interactions between Chemicals?

Some chemicals are known to enhance the effect of other chemicals beyond what can be predicted with standard mixture models, such as concentration addition and independent action. These chemicals are called synergists. Up until now, no models exist that can predict the joint effect of mixtures including synergists. The aim of the present study is to develop a mechanistic toxicokinetic (TK) and toxicodynamic (TD) model for the synergistic mixture of the azole fungicide, propiconazole (the synergist), and the insecticide, α-cypermethrin, on the mortality of the crustacean Daphnia magna. The study tests the hypothesis that the mechanism of synergy is the azole decreasing the biotransformation rate of α-cypermethrin and validates the predictive ability of the model on another azole with a different potency: prochloraz. The study showed that the synergistic potential of azoles could be explained by their effect on the biotransformation rate but that this effect could only partly be explained by the effect of the two azoles on cytochrome P450 activity, measured on D. magna in vivo. TKTD models of interacting mixtures seem to be a promising tool to test mechanisms of interactions between chemicals. Their predictive ability is, however, still uncertain.

Mechanistic Understanding of the Synergistic Potential of Azole Fungicides in the Aquatic Invertebrate Gammarus pulex

Azole fungicides are known inhibitors of the important enzyme class cytochrome P450 monooxygenases (CYPs), thereby influencing the detoxification of co-occurring substances via biotransformation. This synergism in mixtures containing an azole has mostly been studied by effect measurements, while the underlying mechanism has been less well investigated. In this study, six azole fungicides (cyproconazole, epoxiconazole, ketoconazole, prochloraz, propiconazole, and tebuconazole) were selected to investigate their synergistic potential and their CYP inhibition strength in the aquatic invertebrate Gammarus pulex. The strobilurin fungicide azoxystrobin was chosen as co-occurring substrate, and the synergistic potential was measured in terms of internal concentrations of azoxystrobin and associated biotransformation products (BTPs). Azoxystrobin is biotransformed by various reactions, and 18 BTPs were identified. By measuring internal concentrations of azoxystrobin and its BTPs with high-resolution tandem mass spectrometry in the presence and absence of azole fungicides followed by toxicokinetic modeling, we showed that the inhibition of CYP-catalyzed biotransformation reactions indeed played a role for the observed synergism. However, synergism was only observed for prochloraz at environmentally realistic concentrations. Increased uptake rate constants, an increase in the total internal concentration of azoxystrobin and its BTPs, in vivo assays for measuring CYP activities, and G. pulex video-tracking suggested that the 2-fold increase in bioaccumulation, and, thereby, the raised toxicity of azoxystrobin in the presence of prochloraz is not only caused by inhibited biotransformation but even more by increased azoxystrobin uptake induced by hyperactivity.

Uptake, biotransformation and elimination of selected pharmaceuticals in a freshwater invertebrate measured using liquid chromatography tandem mass spectrometry

Methods were developed to assess uptake and elimination kinetics in Gammarus pulex of nine pharmaceuticals (sulfamethazine, carbamazepine, diazepam, temazepam, trimethoprim, warfarin, metoprolol, nifedipine and propranolol) using targeted LC-MS/MS to determine bioconcentration factors (BCFs) using a 96 h toxicokinetic exposure and depuration period. The derived BCFs for these pharmaceuticals did not trigger any regulatory thresholds and ranged from 0 to 73 L kg^-1 (sulfamethazine showed no bioconcentration). Metabolism of chemicals can affect accurate BCF determination through parameterisation of the kinetic models. The added selectivity of LC-MS/MS allowed us to develop confirmatory methods to monitor the biotransformation of propranolol, carbamazepine and diazepam in G. pulex. Varying concentrations of the biotransformed products; 4-hydroxypropranolol sulphate, carbamazepine-10,11-epoxide, nordiazepam, oxazepam and temazepam were measured following exposure of the precursor compounds. For diazepam, the biotransformation product nordiazepam was present at higher concentrations than the parent compound at 94 ng g^-1 dw. Overall, the results indicate that pharmaceutical accumulation is low in these freshwater amphipods, which can potentially be explained by the rapid biotransformation and excretion.

Direct Conjugation of Emerging Contaminants in Arabidopsis: Indication for an Overlooked Risk in Plants?

Agricultural use of treated wastewater, biosolids, and animal wastes introduces a multitude of contaminants of emerging concerns (CECs) into the soil-plant system. The potential for food crops to accumulate CECs depends largely on their metabolism in plants, which at present is poorly understood. Here, we evaluated the metabolism of naproxen and ibuprofen, two of the most-used human drugs from the Profen family, in Arabidopsis thaliana cells and the Arabidopsis plant. The complementary use of high-resolution mass spectrometry and ¹⁴C labeling allowed the characterization of both free and conjugated metabolites, as well as nonextractable residues. Naproxen and ibuprofen, in their parent form, were conjugated quickly and directly with glutamic acid and glutamine, and further with peptides, in A. thaliana cells. For example, after 120 h, the metabolites of naproxen accounted for >90% of the extractable chemical mass, while the intact parent itself was negligible. The structures of glutamate and glutamine conjugates were confirmed using synthesized standards and further verified in whole plants. Amino acid conjugates may easily deconjugate, releasing the parent molecule. This finding highlights the possibility that the bioactivity of such CECs may be effectively preserved through direct conjugation, a previously overlooked risk. Many other CECs are also carboxylic acids, such as the profens. Therefore, direct conjugation may be a common route for plant metabolism of these CECs, making it imperative to consider conjugates when assessing their risks.