Examining Uncertainty in In Vitro-In Vivo Extrapolation Applied in Fish Bioconcentration Models

Many human pharmaceuticals are weak inhibitors of the cytochrome P450 system in rainbow trout (Oncorhynchus mykiss) liver S9 fractions

Introduction

Pharmaceutical residues are widely detected in aquatic environment and can be taken up by nontarget species such as fish. The cytochromes P450 (CYP) represent an important detoxification mechanism in fish, like in humans. In the present study, we assessed the correlation of the substrate selectivities of rainbow trout CYP1A and CYP3A homologues with those of human, through determination of the half-maximal inhibitory concentrations (IC50) of a total sixteen human pharmaceuticals toward CYP1A-like ethoxyresorufin O-deethylase (EROD) and CYP3A-like 7-benzyloxy-4-trifluoromethylcoumarin O-debenzylase (BFCOD) in rainbow trout (Oncorhynchus mykiss) liver S9 fractions (RT-S9).

Methods

The inhibitory impacts (IC50) of atomoxetine, atorvastatin, azelastine, bimatoprost, clomethiazole, clozapine, desloratadine, disulfiram, esomeprazole, felbinac, flecainide, orphenadrine, prazosin, quetiapine, sulpiride, and zolmitriptan toward the EROD and BFCOD activities in RT-S9 were determined using the IC50 shift assay, capable of identifying time-dependent inhibitors (TDI). Additionally, the nonspecific binding of the test pharmaceuticals to RT-S9 was assessed using equilibrium dialysis.

Results

Most test pharmaceuticals were moderate to weak inhibitors of both EROD and BFCOD activity in RT-S9, even if most are noninhibitors of human CYP1A or CYP3A. Only bimatoprost, clomethiazole, felbinac, sulpiride, and zolmitriptan did not inhibit either activity in RT-S9. EROD inhibition was generally stronger than that of BFCOD and some substances (atomoxetine, flecainide, and prazosin) inhibited selectively only EROD activity. The strongest EROD inhibition was detected with azelastine and esomeprazole (unbound IC50 of 3.8 ± 0.5 µM and 3.0 ± 0.8 µM, respectively). None of the test substances were TDIs of BFCOD, but esomeprazole was a TDI of EROD. Apart from clomethiazole and disulfiram, the nonspecific binding of the test pharmaceuticals to the RT-S9 was extensive (unbound fractions <0.5) and correlated well (R 2 = 0.7135) with their water-octanol distribution coefficients.

Discussion

The results indicate that the P450 interactions in RT-S9 cannot be explicitly predicted based on human data, but the in vitro data reported herein can shed light on the substrate selectivity of rainbow trout CYP1A1 and CYP3A27 in comparison to their human homologues. The IC50 concentrations are however many orders of magnitude higher than average environmental concentrations of pharmaceuticals. The time-dependent EROD inhibition by esomeprazole could warrant further research to evaluate its possible interlinkages with hepatotoxic impacts on fish.

A Comparison of In Vitro Metabolic Clearance of Various Regulatory Fish Species Using Hepatic S9 Fractions

Bioaccumulation predictions can be substantially improved by combining in vitro metabolic rate measurements derived from rainbow trout hepatocytes and/or hepatic S9 fractions with quantitative structure-activity relationship (QSAR) modeling approaches. Compared with in vivo testing guidelines Organisation for Economic Co-operation and Development (OECD) 305 and Office of Chemical Safety and Pollution Prevention (OCSPP; an office of the US Environmental Protection Agency) 850.1730, the recently adopted OECD test guidelines 319A and 319B are in vitro approaches that have the potential to provide a time- and cost-efficient, humane solution, reducing animal use while addressing uncertainties in bioaccumulation across species. The present study compares the hepatic clearance of the S9 subcellular fraction of rainbow trout, bluegill, common carp, fathead minnow, and largemouth bass, discerning potential differences in metabolism between different warm- and cold-water species. With refinements to the in vitro metabolic S9 assay for high-throughput analysis, we measured in vitro clearance rates of seven chemicals crossing multiple classes of chemistry and modes of action. We confirmed that data from rainbow trout liver S9 fraction metabolic rates can be utilized to predict rainbow trout bioconcentration factors using an in vitro to in vivo extrapolation model, as intended in the OECD 319B applicability domain per the bioaccumulation prediction. Also, we determined that OECD 319B can be applied to other species, modified according to their habitat, adaptations to feeding behavior, and environmental conditions (e.g., temperature). Once toxicokinetics for each species is better understood and appropriate models are developed, this method can be an excellent tool to determine hepatic clearance and potential bioaccumulation across species. The present study could be leveraged prior to or in place of initiating in vivo bioconcentration studies, thus optimizing selection of appropriate fish species. Environ Toxicol Chem 2024;43:1390-1405. © 2024 SETAC.

Organ-specific biotransformation in salmonids: Insight into intrinsic enzyme activity and biotransformation of three micropollutants

Aquatic ecosystems continue to be threatened by chemical pollution. To what extent organisms are able to cope with chemical exposure depends on their ability to display mechanisms of defense across different organs. Among these mechanisms, biotransformation processes represent key physiological responses that facilitate detoxification and reduce the bioaccumulation potential of chemicals. Biotransformation does not only depend on the ability of different organs to display biotransformation enzymes but also on the affinity of chemicals towards these enzymes. In the present study, we explored the ability of different organs and of two freshwater fish to support biotransformation processes through the determination of in vitro phase I and II biotransformation enzyme activity, and their role in supporting intrinsic clearance and the formation of biotransformation products. Three environmentally relevant pollutants were evaluated: the polycyclic aromatic hydrocarbon (PAH) pyrene (as recommended by the OECD 319b test guideline), the fungicide azoxystrobin, and the pharmaceutical propranolol. Comparative studies using S9 sub-cellular fractions derived from the liver, intestine, gills, and brain of brown trout (Salmo trutta) and rainbow trout (Oncorhynchus mykiss) revealed significant phase I and II enzyme activity in all organs. However, organ- and species-specific differences were found. In brown trout, significant extrahepatic biotransformation was observed for pyrene but not for azoxystrobin and propranolol. In rainbow trout, the brain appeared to biotransform azoxystrobin. In this same species, propranolol appeared to be biotransformed by the intestine and gills. Biotransformation products could be detected only from hepatic biotransformation, and their profiles and formation rates displayed species-specific patterns and occurred at different magnitudes. Altogether, our findings further contribute to the current understanding of organ-specific biotransformation capacity, beyond the expression and activity of enzymes, and its dependence on specific enzyme-chemical interactions to support mechanisms of defense against exposure.

Bioconcentration Assessment in Fish Based on In Vitro Intrinsic Clearance: Predictivity of an Empirical Model Compared to In Vitro-In Vivo Extrapolation Models

To estimate the bioconcentration factor (BCF), the in vitro intrinsic clearance (CLIN VITRO,INT) from rainbow trout liver S9 fractions (RT-S9) can be applied to in vitro-in vivo extrapolation (IVIVE) models, yet uncertainties remain in model parameterization. An alternative model approach is evaluated: a regression model was built in the form log BCF = a × log Kow + b × log CLIN VITRO,INT. The coefficients a and b were fitted based on a training set of 40 chemicals. A high robustness of the coefficients and good accuracy of BCF prediction were found on independent datasets of neutral organic chemicals (measured log Kow 3.3-6.2). BCF predictions were similar to or in better agreement with in vivo BCFs compared to IVIVE models (2.4- to 2.9- vs 2.8- to 3.6-fold misprediction) for training and test sets. Species-matched models (trout, carp) did not result in improvements. This study presents the largest dataset on CLIN VITRO,INT and BCFs to assess predictivity of the RT-S9 assay. The robustness of the regression statistics on different datasets and the high statistical weight of the CLIN VITRO,INT term illustrate the predictive power of the RT-S9 assay as an important step toward regulatory acceptance to replace animal experiments.

Kidney cell culture Macaca fascicularis as a candidate for vaccine development and in vitro model

BACKGROUND

Cell culture is the proliferation of a cell population in vitro by isolating from the original tissue or growing from existing ones. One essential source is the monkey kidney cell cultures which have an essential role in biomedical study. This is due to the significant homology between the human and macaque genomes making these useful for cultivating human viruses, especially enteroviruses, and growing vaccines.

METHODS

This study developed cell cultures derived from the kidney of Macaca fascicularis (Mf) and validated its gene expression.

RESULTS

The primary cultures were successfully subcultured up to six passages, grew as monolayers, and exhibited epithelial-like morphology. The cultured cells remained heterogeneous in phenotype and they expressed CD155 and CD46 as viral receptors, cell morphology (CD24, endosialin, and vWF), proliferation, also apoptosis markers (Ki67 and p53).

CONCLUSIONS

These results indicated that the cell cultures can be used as in vitro model cells for vaccine development and bioactive compound.

Miniaturized method for the quantification of persistent organic pollutants and their metabolites in HepG2 cells: assessment of their biotransformation

Biotransformation can greatly influence the accumulation and, subsequently, toxicity of substances in living beings. Although traditionally these studies to quantify metabolization of a compound have been carried out with in vivo species, currently, in vitro test methods with very different cell lines are being developed for their evaluation. However, this is still a very limited field due to multiple variables of a very diverse nature. So, an increasing number of analytical chemists are working with cells or other similar biological samples of very small size. This makes it necessary to address the development of analytical methods that allow determining their concentration both inside the cells and in their exposure medium. The aim of this study is to develop a set of analytical methodologies for the quantification of polycyclic aromatic hydrocarbons, PAHs (phenanthrene, PHE), and polybrominated diphenyl ethers, PBDEs (2,2',4,4'-tetrabromodiphenyl ether, BDE-47), and their major metabolites in cells and their exposure medium. Analytical methodologies, based on miniaturized ultrasound probe-assisted extraction, gas chromatography-mass spectrometry-microelectron capture detector (GC-MS-µECD), and liquid chromatography-fluorescence detector (LC-FL) determination techniques, have been optimized and then applied to a biotransformation study in HepG2 at 48 h of exposure. Significant concentrations of the major metabolites of PHE (1-OH, 2-OH, 3-OH, 4-OH-, and 9-OH-PHE) and BDE-47 (5-MeO-, 5-OH-, and 3-OH-BDE-47) were detected and quantified inside the cells and in the exposure medium. These results provide a new method for determination and improve information on the metabolization ratios for a better knowledge of the metabolic pathways and their toxicity.

An amended in vitro-in vivo extrapolation model that accounts for first pass clearance effects on chemical bioaccumulation in fish

Measured rates of in vitro intrinsic clearance for fish may be extrapolated to the whole animal as a means of estimating a whole-body biotransformation rate constant (k_B; d^-1). This estimate of k_B can then be used as an input to existing bioaccumulation prediction models. Most in vitro-in vivo extrapolation/bioaccumulation (IVIVE/B) modeling efforts to date have focused on predicting the chemical bioconcentration in fish (aqueous only exposure), with less attention paid to dietary exposures. Following dietary uptake, biotransformation in the gut lumen, intestinal epithelia, and liver can reduce chemical accumulation; however, current IVIVE/B models do not consider these first pass clearance effects on dietary uptake. Here we present an amended IVIVE/B model that accounts for first pass clearance. The model is then used to examine how biotransformation in the liver and intestinal epithelia (alone or combined) may impact chemical accumulation that occurs during dietary exposure. First pass clearance by the liver can greatly reduce dietary uptake of contaminants, but these effects are only apparent at rapid rates of in vitro biotransformation (first order depletion rate constant k_DEP ≥ 10 h^-1). The impact of first pass clearance becomes more pronounced when biotransformation in the intestinal epithelia is included in the model. Modelled results suggest that biotransformation in the liver and intestinal epithelia cannot entirely explain reduced dietary uptake reported in several in vivo bioaccumulation tests. This unexplained reduction in dietary uptake is attributed to chemical degradation in the gut lumen. These findings underscore the need for research to directly investigate luminal biotransformation in fish.

Trophodynamic of endocrine disrupting compounds in the aquatic food webs: Association with hydrophobicity and biota metabolic rate

Increasing concentration of endocrine disrupting compounds (EDCs) are released into the aquatic environment, resulting in irreversible effects on the endocrine and reproductive systems of biota. How the liver enzymes affect metabolic rate of these compounds and thus their structure-related trophic transfer in aquatic food webs remains largely unknown. In this study, the concentrations of seven common EDCs were measured in 15 species of fish, 7 invertebrate species and plankton collected from Liuxi River to Pearl River, South China. The mean ΣEDC concentrations generally were found to increase as follows: plankton (29.59 ng g^-1 dw) < invertebrate species (50.69 ng g^-1 dw) < fish (122.56 ng g^-1 dw), with 4-nonylphenol (4-NP) and bisphenol S (BPS) as the predominant components. Trophic magnification factors (TMFs) values were >1.0 ranged from 1.30 (BPS) to 4.07 (4-NP), indicating trophic magnification potential. Measurement of metabolism and activities of microsomal CYP450 enzymes were performed in the fish liver microsomes of Hypophthalmichthys molitrix ([TL] = 2.27), Cirrhinus mrigala (TL = 3.87) and Odontamblyopus rubicundus (TL = 4.73). TMFs were significantly negatively correlated with the obtained in vitro biotransformation clearance rates (CL _{in vitro}) of EDCs and CYP450 enzymes activities. A multiple linear regression model indicated that biotransformation clearance is a more powerful predictor for TMFs than the hydrophobicity (Kow) to drive changes in the studied aquatic food web trophodynamics.

Investigating the bioaccumulation potential of anionic organic compounds using a permanent rainbow trout liver cell line

Permanent rainbow trout (Oncorhynchus mykiss) cell lines represent potential in vitro alternatives to experiments with fish. We here developed a method to assess the bioaccumulation potential of anionic organic compounds in fish, using the rainbow trout liver-derived RTL-W1 cell line. Based on the availability of high quality in vivo bioconcentration (BCF) and biomagnification (BMF) data and the substances' charge state at physiological pH, four anionic compounds were selected: pentachlorophenol (PCP), diclofenac (DCF), tecloftalam (TT) and benzotriazol-tert-butyl-hydroxyl-phenyl propanoic acid (BHPP). The fish cell line acute toxicity assay (OECD TG249) was used to derive effective concentrations 50 % and non-toxic exposure concentrations to determine exposure concentrations for bioaccumulation experiments. Bioaccumulation experiments were performed over 48 h with a total of six time points, at which cell, medium and plastic fractions were sampled and measured using high resolution tandem mass spectrometry after online solid phase extraction. Observed cell internal concentrations were over-predicted by KOW-derived predictions while pH-dependent octanol-water partitioning (DOW) and membrane lipid-water partitioning (DMLW) gave better predictions of cell internal concentrations. Measured medium and cell internal concentrations at steady state were used to calculate RTL-W1-based BCF, which were compared to DOW- or DMLW-based model approaches and in vivo data. With the exception of PCP, the cell-derived BCF best compared to DOW-based model predictions, which were higher than predictions based on DMLW. All methods predicted the in vivo BCF for diclofenac well. For PCP, the cell-derived BCF was lowest although all BCF predictions underestimated the in vivo BCF by ≥ 1 order of magnitude. The RTL-W1 cells, and all other prediction methods, largely overestimated in vivo BMF, which were available for PCP, TT and BHPP. We conclude that the RTL-W1 cell line can supplement BCF predictions for anionic compounds. For BMF estimations, however, in vitro-in vivo extrapolations need adaptation or a multiple cell line approach.

Bioaccumulation and Biotransformation of BDE-47 Using Zebrafish Eleutheroembryos (Danio rerio)

Polybrominated diphenyl ethers (PBDEs) are well-known endocrine disrupting chemicals identified as organic persistent pollutants. Their metabolites OH-BDE and MeO-BDE have been reported to be potentially more toxic than the postulated precursor PBDEs. One of the most predominant congeners of PBDEs in the environment is BDE-47, due to its high presence in industrially used mixtures. In the present study, the bioaccumulation and biotransformation of BDE-47 into its major metabolites is evaluated using zebrafish (Danio rerio) eleutheroembryos adapting a previously developed alternative method to bioconcentration official guideline Organisation for Economic Co-ordination and Development 305, which reduces the animal suffering, time, and cost. For the simultaneous determination of BDE-47 and its metabolites in larvae and exposure medium, and considering the polarity difference of the analytes and the small sample size, the development of a validated analytical method is a step to ensure quality results. In the present study, an ultrasound-assisted extraction followed by a solid phase extraction dispersive clean-up step and gas chromatography-mass spectrometry-microelectron capture detector (GC-MS-μECD) with a previous derivatization process was optimized and validated. Bioconcentration factors (BCFs) were calculated using a first-order one-compartment toxicokinetic model. The profiles found show rapid absorption in the first hours of larval development and great bioaccumulative capacity, finding BCFs of 7294 ± 899 and 36 363 ± 5702 at nominal concentrations of 10 and 1 μg L^-1 , respectively. Metabolization studies show increasing concentrations of the metabolites BDE-28, 2'-OH-BDE-28, and 5-MeO-BDE-47 throughout the exposure time. The results obtained show the feasibility of the method for bioaccumulation and open up the possibility of metabolic studies with zebrafish eleutheroembryos, which is a very underdeveloped field without official testing or regulation. Environ Toxicol Chem 2023;42:835-845. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Models Used to Predict Chemical Bioaccumulation in Fish from in Vitro Biotransformation Rates Require Accurate Estimates of Blood-Water Partitioning and Chemical Volume of Distribution

Methods for extrapolating measured in vitro intrinsic clearance to a whole-body biotransformation rate constant (kB ) have been developed to support modeled bioaccumulation assessments for fish. The inclusion of extrapolated kB values into existing bioaccumulation models improves the prediction of chemical bioconcentration factors (BCFs), but there remains a tendency for these methods to overestimate BCFs relative to measured values. Therefore, a need exists to evaluate the extrapolation procedure to assess potential sources of error in predicted kB values. We examined how three different approaches (empirically based, composition based, and polyparameter linear free energy relationships [ppLFERs]) used to predict chemical partitioning in vitro (liver S9 system; KS9W ), in blood (KBW ), and in whole fish tissues (KFW ) impact the prediction of a chemical's hepatic clearance binding term (fU ) and apparent volume of distribution (VD ), both of which factor into the calculation of kB and the BCF. Each approach yielded different KS9W , KBW , and KFW values, but resulted in fU values that were of similar magnitude and remained relatively constant at log octanol-water partition ratios (KOW ) greater than 4. This is because KBW and KS9W values predicted by any given approach exhibit a similar dependence on log KOW (i.e., regression slope), which results in a cancelation of "errors" when fU is calculated. In contrast, differences in KBW values predicted by the three approaches translate to differences in VD , and by extension kB and the BCF, which become most apparent at log KOW greater than 6. There is a need to collect KBW and VD data for hydrophobic chemicals in fish that can be used to evaluate and improve existing partitioning prediction approaches in extrapolation models for fish. Environ Toxicol Chem 2023;42:33-45. © 2022 SETAC.

Predicting Hepatic Clearance of Psychotropic Drugs in Isolated Perfused Fish Livers Using a Combination of Two In Vitro Assays

In vitro biotransformation assays with primary trout hepatocytes (RT-HEP) or liver subcellular fractions (RT-S9) have been proposed as valuable tools to help scientists and regulators better understand the toxicokinetics of chemicals. While both assays have been applied successfully to a diversity of neutral organic chemicals, only the RT-S9 assay has been applied to a large number of ionizable organic chemicals. Here, a combination of an in vitro biotransformation assay with RT-HEP with an active transport assay based on the permanent rainbow trout liver cell line RTL-W1 was used to qualitatively predict the potential hepatic clearance of nine psychotropic drugs with various degrees of ionization. Predictions were compared with rates of clearance measured in isolated perfused rainbow trout livers, and the importance of active transport was verified in the presence of the active transport inhibitor cyclosporin A. For the first time, it was demonstrated that a combination of biotransformation and active transport assays is powerful for the prediction of rates of hepatic clearance of ionizable chemicals. Ultimately, it is expected that this approach will allow for use of fewer animals while at the same time improving our confidence in the use of data from in vitro assays in chemical risk assessment.

Bioaccumulation Screening of Neutral Hydrophobic Organic Chemicals in Air-Breathing Organisms Using In Vitro Rat Liver S9 Biotransformation Assays

To advance methods for bioaccumulation assessment of organic substances in air-breathing organisms, the present study developed an in vitro approach for screening neutral hydrophobic organic substances for their bioaccumulation potential in air-breathing organisms consisting of (1) depletion assays for chemicals in rat liver S9 subcellular fractions, (2) in vitro-in vivo extrapolation, and (3) whole-organism bioaccumulation modeling to assess the biomagnification potential of neutral organic substances in the rat. Testing of the in vitro method on 14 test chemicals of potentially biomagnifying substances showed that the bioassays could be conducted with a high level of reproducibility and that in vitro-derived elimination rate constants were in good agreement with in vivo-determined elimination rate constants in the rat. Exploring the potential of the in vitro approach for screening organic chemicals for bioaccumulation in air-breathing organisms indicated that chemical substances that exhibit a depletion rate constant in the S9 in vitro bioassay ≥0.3 h^-1 are not expected to biomagnify in rats independent of their octanol-water partitioning coefficient (K_OW ) or octanol-air partitioning coefficient (K_OA ). The high level of reproducibility achieved in the test, combined with the good agreement between in vitro-derived and in vivo-determined depuration rates, suggests that the in vitro approach in combination with a K_OA - and K_OW -based screening approach has good potential for screening chemicals in commerce for their bioaccumulation potential in air-breathing organisms in a cost-effective and expedient manner, especially if the bioassay can be automated. Environ Toxicol Chem 2022;41:2565-2579. © 2022 SETAC.

In vitro biotransformation assays using fish liver cells: Comparing rainbow trout and carp hepatocytes

Biotransformation assays using primary hepatocytes from rainbow trout, Oncorhynchus mykiss, were validated as a reliable in vitro tool to predict in vivo bioconcentration factors (BCF) of chemicals in fish. Given the pronounced interspecies differences of chemical biotransformation, the present study aimed to compare biotransformation rate values and BCF predictions obtained with hepatocytes from the cold-water species, rainbow trout, to data obtained with hepatocytes of the warm-water species, common carp (Cyprinus carpio). In a first step, we adapted the protocol for the trout hepatocyte assay, including the cryopreservation method, to carp hepatocytes. The successful adaptation serves as proof of principle that the in vitro hepatocyte biotransformation assays can be technically transferred across fish species. In a second step, we compared the in vitro intrinsic clearance rates (CL_{in vitro, int}) of two model xenobiotics, benzo[a]pyrene (BaP) and methoxychlor (MXC), in trout and carp hepatocytes. The in vitro data were used to predict in vivo biotransformation rate constants (k_B) and BCFs, which were then compared to measured in vivo k_B and BCF values. The CL_{in vitro, int} values of BaP and MXC did not differ significantly between trout and carp hepatocytes, but the predicted BCF values were significantly higher in trout than in carp. In contrast, the measured in vivo BCF values did not differ significantly between the two species. A possible explanation of this discrepancy is that the existing in vitro-in vivo prediction models are parameterized only for trout but not for carp. Therefore, future research needs to develop species-specific extrapolation models.

Validation of Methods for in Vitro-in Vivo Extrapolation Using Hepatic Clearance Measurements in Isolated Perfused Fish Livers

In vitro biotransformation assays using hepatocytes or liver subcellular fractions, combined with in vitro-in vivo extrapolation (IVIVE) models, have been proposed as an alternative to live fish bioconcentration studies. The uncertainty associated with IVIVE approaches to date has been attributed to assay protocols, model assumptions, or variability of in vivo data. An isolated perfused trout liver model that measures hepatic clearance has been proposed for validating IVIVE predictions in the absence of other confounding factors. Here, we investigated the hepatic clearances of five chemicals (pyrene, phenanthrene, 4-n-nonlyphenol, deltamethrin, and methoxychlor) in this model and compared measured rates to values predicted from published in vitro intrinsic clearances for validation of IVIVE models. Additionally, we varied protein concentrations in perfusates to test binding assumptions of these models. We found that measured and predicted hepatic clearances were in very good agreement (root mean squared error 16.8 mL h^-1 g^-1) across three levels of protein binding and across a more diverse chemical space than previously studied within this system. Our results show that current IVIVE methods can reliably predict in vivo clearance rates and indicate that discrepancies from measured bioconcentration factors might be driven by other processes, such as extrahepatic biotransformation, etc., and help streamline optimization efforts to the processes that truly matter.

Predicting chemical ecotoxicity by learning latent space chemical representations

In silico prediction of chemical ecotoxicity (HC₅₀) represents an important complement to improve in vivo and in vitro toxicological assessment of manufactured chemicals. Recent application of machine learning models to predict chemical HC₅₀ yields variable prediction performance that depends on effectively learning chemical representations from high-dimension data. To improve HC₅₀ prediction performance, we developed an autoencoder model by learning latent space chemical embeddings. This novel approach achieved state-of-the-art prediction performance of HC₅₀ with R² of 0.668 ± 0.003 and mean absolute error (MAE) of 0.572 ± 0.001, and outperformed other dimension reduction methods including principal component analysis (PCA) (R² = 0.601 ± 0.031 and MAE = 0.629 ± 0.005), kernel PCA (R² = 0.631 ± 0.008 and MAE = 0.625 ± 0.006), and uniform manifold approximation and projection dimensionality reduction (R² = 0.400 ± 0.008 and MAE = 0.801 ± 0.002). A simple linear layer with chemical embeddings learned from the autoencoder model performed better than random forest (R² = 0.663 ± 0.007 and MAE = 0.591 ± 0.008), fully connected neural network (R² = 0.614 ± 0.016 and MAE = 0.610 ± 0.008), least absolute shrinkage and selection operator (R² = 0.617 ± 0.037 and MAE = 0.619 ± 0.007), and ridge regression (R² = 0.638 ± 0.007 and MAE = 0.613 ± 0.005) using unlearned raw input features. Our results highlighted the usefulness of learning latent chemical representations, and our autoencoder model provides an alternative approach for robust HC₅₀ prediction.

New insights on in vitro biotransformation of anticoagulant rodenticides in fish

The assessment of the bioaccumulation potential of chemicals is an essential and mandatory part of their regulatory environmental risk and hazard assessment. So far, in vitro data on fish metabolism is rarely available for biocidal active substances such as anticoagulant rodenticides. In this case study we present in vitro biotransformation rates of eight biocidal and one pharmaceutical anticoagulants in rainbow trout (Oncorhynchus mykiss) liver subcellular S9 fraction (RT-S9) determined following the Organisation for Economic Co-operation and Development test guideline 319B method at two different incubation temperatures (i.e., 12 ± 1 °C and 23 ± 2 °C). Furthermore, we address challenges associated with the usability and interpretation of in vitro data to support the decision making within the regulatory bioaccumulation assessment in bridging the gap between in silico methods and in vivo studies. According to our results, four of the tested substances (i.e., chlorophacinone, coumatetralyl, bromadiolone, and difenacoum) exhibited significant intrinsic clearance (p < .001) in the RT-S9 assay. Overall, the observed metabolism was (very) slow and clearance rates were temperature-dependent. Whether the determined in vitro biotransformation rate had a substantial influence on the predicted bioconcentration factor during extrapolation was subject to the lipophilicity of the test substance. Further improvements of existing concepts are needed to overcome uncertainties in the prediction of bioconcentration factors for chemicals such as anticoagulants.

In vitro-in vivo biotransformation and phase I metabolite profiling of benzo[a]pyrene in Gulf killifish (Fundulus grandis) populations with different exposure histories

Chronic exposure to pollution may lead populations to display evolutionary adaptations associated with cellular and physiological mechanisms of defense against xenobiotics. This could result in differences in the way individuals of the same species, but inhabiting different areas, cope with chemical exposure. In the present study, we explore two Gulf killifish (Fundulus grandis) populations with different exposure histories for potential differences in the biotransformation of benzo[a]pyrene (BaP), and conduct a comparative evaluation of in vitro and in vivo approaches to describe the applicability of new approach methodologies (NAMs) for biotransformation assessments. Pollution-adapted and non-adapted F. grandis were subjected to intraperitoneal (IP) injections of BaP in time-course exposures, prior to measurements of CYP biotransformation activity, BaP liver concentrations, and the identification and quantification of phase I metabolites. Additionally, substrate depletion bioassays using liver S9 fractions were employed for measurements of intrinsic hepatic clearance and to evaluate the production of metabolites in vitro. Pollution-adapted F. grandis presented significantly lower CYP1A activity and intrinsic clearance rates that were 3 to 4 times lower than non-adapted fish. The metabolite profiling of BaP showed the presence of 1‑hydroxy-benzo[a]pyrene in both the in vitro and in vivo approaches but with no significant population differences. Contrarily, 9‑hydroxy-benzo[a]pyrene and benzo[a]pyrene-4,5-dihydrodiol, only identified through the in vivo approach, presented higher concentrations in the bile of pollution-adapted fish relative to non-adapted individuals. These observations further the understanding of the evolutionary adaptation of F. grandis inhabiting heavily polluted environments in the Houston Ship Channel, TX, USA, and highlight the need to consider the evolutionary history of populations of interest during the implementation of NAMs.

Reduced biotransformation of polycyclic aromatic hydrocarbons (PAHs) in pollution-adapted Gulf killifish (Fundulus grandis)

Anthropogenic pollution represents a significant source of selection, potentially leading to the emergence of evolutionary adaptations in chronically exposed organisms. A recent example of this scenario corresponds to Gulf killifish (Fundulus grandis) populations inhabiting the Houston Ship Channel (HSC), Texas, USA, which have been documented to have adapted to this heavily contaminated environment. Although not fully elucidated, one particularly important aspect of their adaptation involves the reduced inducibility of the aryl hydrocarbon receptor (AhR) and, potentially, the alteration of major biotransformation pathways. In the present study, we employed a modified Organization for Economic Cooperation and Development (OECD) 319-B test guideline to explore population and sex-related differences in the hepatic biotransformation of six polycyclic aromatic hydrocarbons (PAHs) in F. grandis populations with different exposure histories. Pollution-adapted F. grandis showed significantly lower hepatic clearance of PAHs than non-adapted fish, especially for high molecular weight PAHs (chrysene, benzo[k]fluoranthene, and benzo[a]pyrene), with pollution-adapted females presenting the lowest clearance. The characterization of different phase I biotransformation enzymes revealed that the basal activity of CYP1A, fundamental in the biotransformation of PAHs, was significantly lower in pollution-adapted fish, especially in females, which showed the lowest activity. Contrarily, basal CYP2C9-like activity was significantly higher in pollution-adapted fish. These results demonstrate the importance of exposure and evolutionary histories in shaping organisms' responses to pollution and provide significant evidence of sex-specific biotransformation differences in F. grandis populations.

Bioaccumulation of current-use herbicides in fish from a global biodiversity hotspot: Lake St Lucia, South Africa

Agricultural pesticides threaten aquatic systems and biodiversity at a global scale, but limited information is available on the accumulation of current-use herbicides in tissues of aquatic organisms. Here, we examine the potential exposure and accumulation of currently used herbicides in two species of fish from Lake St Lucia, a global biodiversity hotspot located in South Africa. Muscle tissue samples were analysed for 11 widely used multi-residue and phenoxy-acid herbicides. Herbicide residues were detected in all samples analysed, with total concentrations ranging from 44.3 to 238 ng g^-1 (Clarias gariepinus) and 72.2-291 ng g^-1 dw (Oreochromis mossambicus). The most prominent herbicides detected included the two phenoxy-acid herbicides, MCPA (17.6 ± 12 ng g^-1) and 2,4-D (28.9 ± 16 ng g^-1), along with acetochlor (15.4 ± 5.8 ng g^-1), atrazine (12.7 ± 7.1 ng g^-1) and terbuthylazine (12.4 ± 12 ng g^-1). Results indicate that fish at Lake St Lucia accumulate a complex mixture of herbicides, some previously unreported in tissue, highlighting the potential threat that agricultural runoff may pose to conservation areas. However, assessing the impact of herbicide accumulation on wild fish populations is difficult at present and urgent toxicological data are needed to better understand chronic exposure effects on aquatic organisms. A preliminary human health risk assessment indicated minimal risk associated with the consumption of local fish, but potential additive and synergistic effects of contaminant mixtures remain unknown.

Bioavailability and toxicity of silver nanoparticles: Determination based on toxicokinetic-toxicodynamic processes

Determining the bioavailability and toxicity mechanism of silver nanoparticles (AgNPs) is challenging as Ag⁺ is continuously released by external or internal AgNP dissolution in the actual exposure system (regardless of the laboratory or the natural environment). Here a novel pulsed-gradient Ag⁺ (AgNO₃) exposure was conducted with zebrafish (Danio rerio) larvae to simulate dissolved gradient concentrations of Ag⁺ from polyvinylpyrrolidone (PVP)-coated AgNPs. The accumulation and toxicity of the pulsed-gradient Ag⁺ (AgNO₃) and, in the meantime, the released Ag⁺ from PVP-AgNPs were predicted using a toxicokinetic-toxicodynamic (TK-TD) model with obtained Ag⁺ parameters. In order to further understand the possible mechanism of PVP-AgNP releasing Ag⁺ in the body, subcellular fractions (S9) of zebrafish were also used to incubate with AgNPs in vitro to mimic the realistic in vivo scenarios. In the TK process, in vivo analysis showed that AgNPs released around twice as many Ag⁺ into the body than were detected with a single Ag⁺ pulse-exposure system; this was supported by evidence that subcellular S9 fractions might cause the PVP-AgNPs to lose the capping agent and favor Ag⁺ release. In the TD process, toxicity (survival rate) was predicted by the total bodily Ag(I) concentration, suggesting that AgNP toxicity in larvae was mainly due to gradually released Ag⁺ rather than AgNPs themselves. This study helps clarify the role of Ag⁺ in AgNP toxicity and offers a novel framework by which to investigate the toxicity of metal nanoparticles and corresponding metal ions in biological systems.

Development of a method for assessing the accumulation and metabolization of antidepressant drugs in zebrafish (Danio rerio) eleutheroembryos

Antidepressant drugs are widely used for the treatment of common mental or other psychiatric disorders such as depression, which affect about 121 million people worldwide. This widespread use has contributed to the input of these pharmaceuticals and their metabolites into the environment. The aim of this work was to develop an analytical method to quantify the most widely used antidepressant drugs, selective serotonin reuptake inhibitors (SSRI), and their main metabolites in the environment. For this, a new and reliable miniaturized extraction method based on dispersive SPE cleanup procedure for extraction of SSRI followed by derivatization with n-heptafluorobutyrylimidazole, and detection by GC-MS was developed. The methodology, including a first-order one-compartment model, was then applied to a bioconcentration study in zebrafish (Danio rerio) eleutheroembryos. The results showed low bioaccumulation of these compounds; however, a biotransformation evidence of the parent compounds into their metabolites was observed after 6 h of exposure. These results indicate the need to integrate metabolic transformation rates to fully model and understand the bioaccumulation patterns of SSRI and their metabolites.

Relevance of desorption kinetics and permeability for in vitro-based predictions of hepatic clearance in fish

The impact of desorption kinetics and permeation kinetics on in vitro-based predictions of in vivo hepatic blood clearances is investigated in the present study. Most commonly, possible limitations due to slow desorption of chemicals from albumin or slow permeation of chemicals through cellular membranes are not considered when in vivo clearances are predicted from in vitro biotransformation rate constants. To evaluate whether the most commonly used extrapolation models might thus overlook important kinetic limitations, we compare predictions of in vivo clearance that explicitly consider desorption and permeation kinetics with predictions of in vivo clearance that neglect these aspects. Our results show that strong limitations due to slow permeation kinetics are possible depending on the assumed permeability value. While permeability values estimated with a mechanistic approach are fast enough to avoid significant limitations, other experimentally derived permeability values lead to dramatically decreased in vivo clearance predictions. These latter values lead to unrealistically low in vivo biotransformation estimates. Furthermore, we also evaluated the implications of desorption kinetics using experimentally determined desorption rate constants. These evaluations show that slow desorption kinetics are unlikely to limit in vivo clearance.

Could chemical exposure and bioconcentration in fish be affected by slow binding kinetics in blood?

The possible implications of slow binding kinetics on respiratory uptake, bioconcentration and exposure of chemicals were evaluated in the present study. Most physiological and chemical information needed for such an evaluation is already known from the literature or can be estimated. However, data for binding kinetics of chemicals in fish plasma have not been reported in the literature yet. In the first part of this study, we therefore experimentally investigated the plasma binding kinetics for ten chemicals, including pollutants like polycyclic aromatic hydrocarbons and a pesticide. The determined desorption rate constants were in the range of 0.4 s-1 to 0.1 s-1. In the second part of this study, we present a comparative modeling analysis of generic predictions with binding kinetics of different velocities. For doing so, a model that explicitly represents binding kinetics in blood was developed and applied for different hypothetical scenarios. The evaluation showed that slow sorption kinetics only limits respiratory uptake and thus influences the levels of bioaccumulation for extreme and, by that, rather unlikely parameter combinations (i.e. for strongly sorbing chemicals with very slow binding kinetics). It can therefore be assumed that limitations on respiratory uptake due to slow binding kinetics in blood are rather unlikely for most chemicals.

Prediction of Unbound Fractions for in Vitro-in Vivo Extrapolation of Biotransformation Data

For in vitro-in vivo extrapolation of biotransformation data, the different sorptive environments in vitro and in vivo need to be considered. The most common approach for doing so is using the ratio of unbound fractions in vitro and in vivo. In the literature, several algorithms for prediction of these unbound fractions are available. In this study, we present a theoretical evaluation of the most commonly used algorithms for prediction of unbound fractions in S9 assays and blood and compare prediction results with empirical values from the literature. The results of this analysis prove a good performance of "composition-based" algorithms, i.e. algorithms that represent the inhomogeneous composition of in vitro assay and in vivo system and describe sorption to the individual components (lipids, proteins, water) in the same way. For strongly sorbing chemicals, these algorithms yield constant values for the ratio of unbound fractions in vitro and in vivo. This is mechanistically plausible, because in these cases, the chemicals are mostly bound, and the ratio of unbound fractions is determined by the volume ratio of sorbing components in both phases.