Predicting concentrations of organic chemicals in fish by using toxicokinetic models

Ecotoxicity of phenol and cresols to aquatic organisms: A review

With the development of industrial production and continuous demand for chemicals, a large volume of wastewater containing phenols was discharged into the aquatic environment. Moreover, chemical leakage further increased the emission of phenols into aquatic systems. Phenol and its methylated derivative (cresols) were selected due to their extensive use in industry and ecotoxicity to freshwater and marine organisms. This review focused on the ecotoxicity of phenol and m-, o-, and p-cresol on aquatic systems. The mechanism of action of phenols was also discussed. The aim of this literature review was to summarise the knowledge of the behaviour, and toxicity on marine and freshwater organisms, of phenols as well as to try to select a series of sensitive biomarkers suitable for ecotoxicological assessment and environmental monitoring in aquatic environments.

Statistical relationship between metabolic decomposition and chemical uptake predicts bioconcentration factor data for diverse chemical exposures

BACKGROUND

A challenge of in vitro to in vivo extrapolation (IVIVE) is to predict the physical state of organisms exposed to chemicals in the environment from in vitro exposure assay data. Although toxicokinetic modeling approaches promise to bridge in vitro screening data with in vivo effects, they are often encumbered by a need for redesign or re-parameterization when applied to different tissues or chemicals.

RESULTS

We demonstrate a parameterization of reverse toxicokinetic (rTK) models developed for the adult zebrafish (Danio rerio) based upon particle swarm optimizations (PSO) of the chemical uptake and degradation rates that predict bioconcentration factors (BCF) for a broad range of chemicals. PSO reveals a relationship between chemical uptake and decomposition parameter values that predicts chemical-specific BCF values with moderate statistical agreement to a limited yet diverse chemical dataset, and all without a need to retrain the model to new data.

CONCLUSIONS

The presented model requires only the octanol-water partitioning ratio to predict BCFs to a fidelity consistent with existing QSAR models. This success begs re-evaluation of the modeling assumptions; specifically, it suggests that chemical uptake into arterial blood may be limited by transport across gill membranes (diffusion) rather than by counter-current flow between gill lamellae (convection). Therefore, more detailed molecular modeling of aquatic respiration may further improve predictive accuracy of the rTK approach.

Reliability of In Vitro Methods Used to Measure Intrinsic Clearance of Hydrophobic Organic Chemicals by Rainbow Trout: Results of an International Ring Trial

In vitro assays are widely employed to obtain intrinsic clearance estimates used in toxicokinetic modeling efforts. However, the reliability of these methods is seldom reported. Here we describe the results of an international ring trial designed to evaluate two in vitro assays used to measure intrinsic clearance in rainbow trout. An important application of these assays is to predict the effect of biotransformation on chemical bioaccumulation. Six laboratories performed substrate depletion experiments with cyclohexyl salicylate, fenthion, 4-n-nonylphenol, deltamethrin, methoxychlor, and pyrene using cryopreserved hepatocytes and liver S9 fractions from trout. Variability within and among laboratories was characterized as the percent coefficient of variation (CV) in measured in vitro intrinsic clearance rates (CLIN VITRO, INT; ml/h/mg protein or 106 cells) for each chemical and test system. Mean intralaboratory CVs for each test chemical averaged 18.9% for hepatocytes and 14.1% for S9 fractions, whereas interlaboratory CVs (all chemicals and all tests) averaged 30.1% for hepatocytes and 22.4% for S9 fractions. When CLIN VITRO, INT values were extrapolated to in vivo intrinsic clearance estimates (CLIN VIVO, INT; l/d/kg fish), both assays yielded similar levels of activity (<4-fold difference for all chemicals). Hepatic clearance rates (CLH; l/d/kg fish) calculated using data from both assays exhibited even better agreement. These findings show that both assays are highly reliable and suggest that either may be used to inform chemical bioaccumulation assessments for fish. This study highlights several issues related to the demonstration of assay reliability and may provide a template for evaluating other in vitro biotransformation assays.

In vitro or not in vitro: a short journey through a long history

The aim of ecotoxicology is to study toxic effects on constituents of ecosystems, with the protection goal being populations and communities rather than individual organisms. In this ecosystem perspective, the use of in vitro methodologies measuring cellular and subcellular endpoints at a first glance appears to be odd. Nevertheless, more recently in vitro approaches gained momentum in ecotoxicology. In this article, we will discuss important application domains of in vitro methods in ecotoxicology. One area is the use of in vitro assays to replace, reduce, and refine (3R) in vivo tests. Research in this field has focused mainly on the use of in vitro cytotoxicity assays with fish cells as non-animal alternative to the in vivo lethality test with fish and on in vitro biotransformation assays as part of an alternative testing strategy for bioaccumulation testing with fish. Lessons learned from this research include the importance of a critical evaluation of the sensitivity, specificity and exposure conditions of in vitro assays, as well as the availability of appropriate in vitro-in vivo extrapolation models. In addition to this classical 3R application, other application domains of in vitro assays in ecotoxicology include the screening and prioritization of chemical hazards, the categorization of chemicals according to their modes of action and the provision of mechanistic information for the pathway-based prediction of adverse outcomes. The applications discussed in this essay may highlight the potential of in vitro technologies to enhance the environmental hazard assessment of single chemicals and complex mixtures at a reduced need of animal testing.

Importance of Toxicokinetics to Assess the Utility of Zebrafish Larvae as Model for Psychoactive Drug Screening Using Meta-Chlorophenylpiperazine (mCPP) as Example

The number of new psychoactive substances (NPS) increases rapidly, harming society and fuelling the need for alternative testing strategies. These should allow the ever-increasing number of drugs to be tested more effectively for their toxicity and psychoactive effects. One proposed strategy is to complement rodent models with zebrafish (Danio rerio) larvae. Yet, our understanding of the toxicokinetics in this model, owing to the waterborne drug exposure and the distinct physiology of the fish, is incomplete. We here explore the toxicokinetics and behavioral effects of an NPS, meta-chlorophenylpiperazine (mCPP), in zebrafish larvae. Uptake kinetics of mCPP, supported by toxicokinetic modeling, strongly suggested the existence of active transport processes. Internal distribution showed a dominant accumulation in the eye, implying that in zebrafish, like in mammals, melanin could serve as a binding site for basic drugs. We confirmed this by demonstrating significantly lower drug accumulation in two types of hypo-pigmented fish. Comparison of the elimination kinetics between mCPP and previously characterized cocaine demonstrated that drug affinities to melanin in zebrafish vary depending on the structure of the test compound. As expected from mCPP-elicited responses in rodents and humans, zebrafish larvae displayed hypoactive behavior. However, significant differences were seen between zebrafish and rodents with regard to the concentration-dependency of the behavioral response and the comparability of tissue levels, corroborating the need to consider the organism-internal distribution of the chemical to allow appropriate dose modeling while evaluating effects and concordance between zebrafish and mammals. Our results highlight commonalities and differences of mammalian versus the fish model in need of further exploration.

Biotransformation of Benzo[ a]pyrene by Three Rainbow Trout ( Onchorhynchus mykiss) Cell Lines and Extrapolation To Derive a Fish Bioconcentration Factor

Permanent fish cell lines constitute a promising complement or substitute for fish in the environmental risk assessment of chemicals. We demonstrate the potential of a set of cell lines originating from rainbow trout ( Oncorhynchus mykiss) to aid in the prediction of chemical bioaccumulation in fish, using benzo[ a]pyrene (BaP) as a model chemical. We selected three cell lines from different tissues to more fully account for whole-body biotransformation in vivo: the RTL-W1 cell line, representing the liver as major site of biotransformation, and the RTgill-W1 (gill) and RTgutGC (intestine) cell lines, as important environment-organism interfaces, which likely influence chemical uptake. All three cell lines were found to effectively biotransform BaP. However, rates of in vitro clearance differed, with the RTL-W1 cell line being most efficient, followed by RTgutGC. Co-exposures with α-naphthoflavone as potent inhibitor of biotransformation, assessment of CYP1A catalytic activity, and the progression of cellular toxicity upon prolonged BaP exposure revealed that BaP is handled differently in the RTgill-W1 compared to the other two cell lines. Application of the cell-line-derived in vitro clearance rates into a physiology-based toxicokinetic model predicted a BaP bioconcentration factor (BCF) of 909-1057 compared to 920 reported for rainbow trout in vivo.

Tissue distribution, bioaccumulation characteristics and health risk of antibiotics in cultured fish from a typical aquaculture area

The negative impacts of residual antibiotics in the environment on ecosystem and human health are big concerns. However, little information is available on the antibiotic bioaccumulation in aquaculture farms. In this study, the bioaccumulative potentials of 21 antibiotics in the plasma, bile, liver and muscle of cultured fish from a typical aquaculture area were systematically investigated.

RESULTS

indicated that antibiotic distribution in the cultured fish was mainly influenced by species and specific substances. The mean values of log bioaccumulation factors (Log BAFs) for the detected antibiotics were in the range of 0.43-3.70, 0.36-4.75, -0.31-4.48, and 0.23-4.33 in the fish plasma, bile, liver and muscle tissues, respectively. For grass carp, both ciprofloxacin and enrofloxacin showed high transportability from the plasma to the muscle and liver. The correlations of various antibiotic concentrations between the plasma and the other three tissues indicated that the concentrations of ciprofloxacin and enrofloxacin in the fish tissues could be predicted by their concentrations in the plasma. Based on the calculated hazard quotients, human health risk evaluation of antibiotic exposure by fish consumption indicated that the consumption of these cultured fish posed low risks to human health.

The influence of chemical degradation during dietary exposures to fish on biomagnification factors and bioaccumulation factors

The chemical dietary absorption efficiency (E_D) quantifies the amount of chemical absorbed by an organism relative to the amount of chemical an organism is exposed to following ingestion. In particular, E_D can influence the extent of bioaccumulation and biomagnification for hydrophobic chemicals. A new E_D model is developed to quantify chemical process rates in the gastrointestinal tract (GIT). The new model is calibrated with critically evaluated measured E_D values (n = 250) for 80 hydrophobic persistent chemicals. The new E_D model is subsequently used to estimate chemical reaction rate constants (k_R) assumed to occur in the lumen of the GIT from experimental dietary exposure tests (n = 255) for 165 chemicals. The new k_R estimates are corroborated with k_R estimates for the same chemicals from the same data derived previously by other methods. The roles of k_R and the biotransformation rate constant (k_B) on biomagnification factors (BMFs) determined under laboratory test conditions and on BMFs and bioaccumulation factors (BAFs) in the environment are examined with the new model. In this regard, differences in lab and field BMFs are highlighted. Recommendations to address uncertainty in E_D and k_R data are provided.

Evaluating parameter availability for physiologically based pharmacokinetic (PBPK) modeling of perfluorooctanoic acid (PFOA) in zebrafish

Physiologically based pharmacokinetic (PBPK) models are considered useful tools to describe the absorption, distribution, metabolism and excretion of xenobiotics. For accurate predictions, PBPK models require species-specific and compound-specific parameters. Zebrafish are considered an appropriate vertebrate model for investigating the toxicity of a wide variety of compounds. However, no specific mechanistic model exists for the pharmacokinetics of perfluoroalkyl acids (PFAAs) in zebrafish, despite growing concern about this class of ubiquitous environmental contaminants. The purpose of this study was to evaluate the current state of knowledge for the parameters that would be needed to construct such a model for zebrafish. We chose perfluorooctanoic acid (PFOA) as a model PFAA with greater data availability. We have updated a previous PBPK model for rainbow trout to simulate PFOA fate in zebrafish following waterborne exposure. For the first time, the model considers hepatobiliary circulation. In order to evaluate the availability of parameters to implement this model, we performed an extensive literature review to find zebrafish-specific parameters. As in previous approaches, we broadened our search to include mammalian and other fish studies when zebrafish-specific data were lacking. Based on the method used to measure or estimate parameters, or based on their species-specific origin, we scored and ranked the quality of available parameters. These scores were then used in Monte Carlo and partial rank correlation analyses to identify the most critical data gaps. The liver, where fatty acid binding proteins (FABPs) and plasma proteins are considered, represented the best model-data agreement. Lack of agreement in other tissues suggest better parameters are needed. The results of our study highlight the lack of zebrafish-specific parameters. Based on sensitivity and uncertainty analysis, parameters associated with PFAA-protein interactions and passive diffusion need further refinement to enable development of predictive models for these emerging chemicals in zebrafish.

Can poly-parameter linear-free energy relationships (pp-LFERs) improve modelling bioaccumulation in fish?

A wide range of studies have characterized different types of biosorbent, with regard to their interactions with chemicals. This has resulted in the development of poly-parameter linear free energy relationships (pp-LFERs) for the estimation of partitioning of neutral organic compounds to biological phases (e.g., storage lipids, phospholipids and serum albumins). The aims of this study were to explore and evaluate the influence of implementing pp-LFERs both into a one-compartment fish model and a multi-compartment physiologically based toxicokinetic (PBTK) fish model and the associated implications for chemical risk assessment. For this purpose, fish was used as reference biota, due to their important role in aquatic food chains and dietary exposure to humans. The bioconcentration factor (BCF) was utilized as the evaluation metric. Overall, our results indicated that models incorporating pp-LFERs (R² = 0.75) slightly outperformed the single parameter (sp) LFERs approach in the one-compartmental fish model (R² = 0.72). A pronounced enhancement was achieved for compounds with log K_OW between 4 and 5 with increased R² from 0.52 to 0.71. The minimal improvement was caused by the overestimation of lipid contribution and underestimation of protein contribution by the sp-approach, which cancelled each other out. Meanwhile, a greater improvement was observed for multi-compartmental PBTK models with consideration of metabolism, making all predictions fall within a factor of 10 compared with measured data. For screening purposes, the K_OW-based (sp-LFERs) approach should be sufficient to quantify the main partitioning characteristics. Further developments are required for the consideration of ionization and more accurate quantification of biotransformation in biota.

Assessing the impact of benzo[a]pyrene with the in vitro fish gut model: An integrated approach for eco-genotoxicological studies

In vitro models are emerging tools for reducing reliance on traditional toxicity tests, especially in areas where information is sparse. For studies of fish, this is especially important for extrahepatic organs, such as the intestine, which, until recently, have been largely overlooked in favour of the liver or gill. Considering the importance of dietary uptake of contaminants, the rainbow trout (Oncorhynchus mykiss) intestine-derived cell line RTgutGC was cultured, to test its suitability as a high-throughput in vitro model. Benzo[a]pyrene (B[a]P) is an important contaminant and a model polycyclic aromatic hydrocarbon (PAH). Over 48 h exposure, a range of endpoints and xenobiotic metabolism rates were examined at three different pH levels indicative of the in vitro (pH 7.5) and in vivo mid-gut (pH 7.7) and hind-gut (pH 7.4) regions as a function of time. These endpoints included (i) cell viability: acid phosphatase (APH) and lactate dehydrogenase (LDH) assays; (ii) glucose uptake; (iii) cytochrome P450 enzyme activity: 7-ethoxyresoorufin-O-deethylase (EROD) assay; (iv) glutathione transferase (GST) activity; (v) genotoxic damage determined using the comet assay. Absence of cell viability loss, in parallel with decrease in the parent compound (B[a]P) in the medium and its subsequent increase in the cells suggested active sequestration, biotransformation, and removal of this representative PAH. With respect to genotoxic response, significant differences were observed at both the sampling times and the two highest concentrations of B[a]P. No significant differences were observed for the different pH conditions. Overall, this in vitro xenobiotic metabolism system appears to be a robust model, providing a basis for further development to evaluate metabolic and toxicological potential of contaminants without use of animals.

Calculating the first-order kinetics of three coupled, reversible processes

Calculation of the concentration-time profile in four consecutive, well-mixed compartments that are connected by diffusional transport is a frequently occurring problem for chemists and engineers. Mathematically this is equivalent to many other problems such as the concentration profiles of a parent compound and its three consecutive reaction products resulting from reversible first-order kinetics. Here we present an analytical solution to this problem implemented in a Microsoft Excel spreadsheet (available for download) and we discuss various examples of how this simple-to-use tool can be applied to very different scenarios from various fields of science.

Aquatic ecotoxicological models and their applicability in Arctic regions

Dose-response modeling is one of the most important steps of ecological risk assessment. It requires concentration-effects relationships for the species under consideration. There are very limited studies and experimental data available for the Arctic aquatic species. Lack of toxicity data hinders obtaining dose-response relationships for lethal (LC50 values), sub-lethal and carcinogenic effects. Gaps in toxicity data could be filled using a variety of in-silico ecotoxicological methods. This paper reviews the suitability of such methods for the Arctic scenario. Mechanistic approaches like toxicokinetic and toxicodynamic analysis are found to be better suited for interspecies extrapolation than statistical methods, such as Quantitative Structure-Activity Relationships/Quantitative Structure Activity-Activity Relationship, ICE, and other empirical models, such as Haber's law and Ostwald's equation. A novel approach is proposed where the effects of the toxicant exposure are quantified based on the probability of cellular damage and metabolites interactions. This approach recommends modeling cellular damage using a toxicodynamic model and physiology or metabolites interactions using a toxicokinetic model. Together, these models provide more reliable estimates of toxicity in the Arctic aquatic species, which will assist in conducting ecological risk assessment of Arctic environment.

A toxicokinetic model for fish including multiphase sorption features

In the scientific field of physiologically based toxicokinetic modeling the complexity of the model used depends on the complexity of the problem to be handled, leading to a broad range of existing models from simple 1-box models to complex multicompartment models. Most of these models work with lumped parameters, for example, an uptake efficiency parameter that can only be obtained with a fit of experimental data. The authors' goal was a model that is completely based on well-defined physiological and physicochemical parameters. Lumped parameters fitted on training data sets would limit the model's applicability. This would enable a new view on process understanding for uptake, distribution, and elimination procedures. Eventual goals are a better localization of chemicals within the organism itself, and to set the stage for future extensions toward ionic compounds and active transport across membranes. The model evaluation reported in the present study has shown that uptake, clearance, and bioaccumulation data for nonpolar chemicals are well predicted. Environ Toxicol Chem 2017;36:1538-1546. © 2016 SETAC.

Toxicodynamic modeling of zebrafish larvae to metals using stochastic death and individual tolerance models: comparisons of model assumptions, parameter sensitivity and predictive performance

Process-based toxicodynamic (TD) models are playing an increasing role in predicting chemical toxicity to aquatic organism. Stochastic death (SD) and individual tolerance distribution (IT) are two often used assumptions in TD models which could lead to different consequences for risk assessment of chemicals. Here, using the toxicity data of single (Cu, Zn, Cd, and Pb) and their binary metal mixtures on survival of zebrafish larvae, we assessed the parameter sensitivity and evaluated the predictive performance of SD and IT models. The sensitivity analysis indicated the parameters related to toxicodynamics such as k _k and threshold, had a great influence on the SD model's output and α had a great influence on the IT model's output. The predicted survival probability was highly sensitive to the assumptions of SD or IT models, and the SD model explained toxicity of single metal and binary metal mixtures better than IT model. Our results suggested that SD model is more suitable in assessing the metal toxicity to zebrafish larvae. Moreover, different combinations of laboratory metal-specific and species-specific experiments with SD and IT models need further study for better understanding and predicting toxic effects for different metals and organisms.

Toxicokinetic models and related tools in environmental risk assessment of chemicals

Environmental risk assessment of chemicals for the protection of ecosystems integrity is a key regulatory and scientific research field which is undergoing constant development in modelling approaches and harmonisation with human risk assessment. This review focuses on state-of-the-art toxicokinetic tools and models that have been applied to terrestrial and aquatic species relevant to environmental risk assessment of chemicals. Both empirical and mechanistic toxicokinetic models are discussed using the results of extensive literature searches together with tools and software for their calibration and an overview of applications in environmental risk assessment. These include simple tools such as one-compartment models, multi-compartment models to physiologically-based toxicokinetic (PBTK) models, mostly available for aquatic species such as fish species and a number of chemical classes including plant protection products, metals, persistent organic pollutants, nanoparticles. Data gaps and further research needs are highlighted.

Pharmaceutical Metabolism in Fish: Using a 3-D Hepatic In Vitro Model to Assess Clearance

At high internal doses, pharmaceuticals have the potential for inducing biological/pharmacological effects in fish. One particular concern for the environment is their potential to bioaccumulate and reach pharmacological levels; the study of these implications for environmental risk assessment has therefore gained increasing attention. To avoid unnecessary testing on animals, in vitro methods for assessment of xenobiotic metabolism could aid in the ecotoxicological evaluation. Here we report the use of a 3-D in vitro liver organoid culture system (spheroids) derived from rainbow trout to measure the metabolism of seven pharmaceuticals using a substrate depletion assay. Of the pharmaceuticals tested, propranolol, diclofenac and phenylbutazone were metabolised by trout liver spheroids; atenolol, metoprolol, diazepam and carbamazepine were not. Substrate depletion kinetics data was used to estimate intrinsic hepatic clearance by this spheroid model, which was similar for diclofenac and approximately 5 fold higher for propranolol when compared to trout liver microsomal fraction (S9) data. These results suggest that liver spheroids could be used as a relevant and metabolically competent in vitro model with which to measure the biotransformation of pharmaceuticals in fish; and propranolol acts as a reproducible positive control.

Numerical evaluation of bioaccumulation and depuration kinetics of PAHs in Mytilus galloprovincialis

Polycyclic aromatic hydrocarbons (PAHs) are important organic pollutants in the aquatic environment due to their persistence and bioaccumulation potential both in organisms and in sediments. Benzo(a)anthracene (BaA) and phenanthrene (PHE), which are in the priority pollutant list of the U.S. EPA (Environmental Protection Agency), are selected as model compounds of the present study. Bioaccumulation and depuration experiments with local Mediterranean mussel species, Mytilus galloprovincialis were used as the basis of the study. Mussels were selected as bioindicator organisms due to their broad geographic distribution, immobility and low enzyme activity. Bioaccumulation and depuration kinetics of selected PAHs in Mytilus galloprovincialis were described using first order kinetic equations in a three compartment model. The compartments were defined as: (1) biota (mussel), (2) surrounding environment (seawater), and (3) algae (Phaeodactylum tricornutum) as food source of the mussels. Experimental study had been performed for three different concentrations. Middle concentration of the experimental data was used as the model input in order to represent other high and low concentrations of selected PAHs. Correlations of the experiment and model data revealed that they are in good agreement. Accumulation and depuration trend of PAHs in mussels regarding also the durations can be estimated effectively with the present study. Thus, this study can be evaluated as a supportive tool for risk assessment in addition to monitoring studies.

From the Cover: Zebrafish Larvae Are Insensitive to Stimulation by Cocaine: Importance of Exposure Route and Toxicokinetics

Zebrafish (Danio rerio) larvae have been suggested as vertebrate model to complement or even replace mammals for rapidly assessing behavioral effects of psychoactive drugs. Yet, divergent responses have been reported in mammals and fish despite the conservation of many drug targets. Cocaine, eg, acts as stimulant in mammals but no such response has been documented for zebrafish larvae. We hypothesized that differences in exposure routes (inhalation or injection in mammals vs waterborne in fish) may be a reason for differences in behavioral responses. We characterized cocaine toxicokinetics by liquid chromatography-mass spectrometry and found its rapid uptake into larvae. We used Matrix-assisted laser desorption ionization-mass spectrometry imaging for the first time to characterize internal distribution of cocaine in zebrafish larvae. Surprisingly, eyes accumulated the highest amount of cocaine and retained most of it even after 48 h depuration. We attribute this to trapping by pigment melanin, a thus far little explored mechanism that may also be relevant for other basic drugs. Cocaine also reached the brain but with levels similar to those in trunk indicating simple passive diffusion as means of distribution which was supported by toxicokinetic models. Although brain levels covered those known to cause hyperactivity in mammals, only hypoactivity (decreased locomotion) was recorded in zebrafish larvae. Our results therefore point to cocaine's anesthetic properties as the dominant mechanism of interaction in the fish: upon entry through the fish skin and gills, it first acts on peripheral nerves rapidly overriding any potential stimulatory response in the brain.

Tralopyril bioconcentration and effects on the gill proteome of the Mediterranean mussel Mytilus galloprovincialis

Antifouling (AF) systems are used worldwide as one of the most cost-effective ways of protecting submerged structures against heavy biofouling. The emergence of environmentally friendly AF biocides requires knowledge on their environmental fate and toxicity. In this study we measured the bioconcentration of the emerging AF biocide tralopyril (TP) in the Mediterranean mussel Mytilus galloprovincialis and investigated the effects of TP on the mussel gill proteome following acute (2days) and chronic (30days) exposure, as well as after a 10-day depuration period. The experiments were carried out with 1μg/L TP; blank and solvent (5×10(-5)% DMSO) controls were also included. Proteomics analysis was performed by mass spectrometry-based multidimensional protein identification technology (MudPIT). Differentially expressed proteins were identified using a label-free approach based on spectral counts and G-test. Our results show that TP is rapidly accumulated by mussels at concentrations up to 362ng/g dw (whole tissues), reaching steady-state condition within 13days. Ten days of depuration resulted in 80% elimination of accumulated TP from the organism, suggesting that a complete elimination could be reached with longer depuration times. In total, 46 proteins were found to be regulated in the different exposure scenarios. Interestingly, not only TP but also DMSO alone significantly modulated the protein expression in mussel gills following acute and chronic exposure. Both compounds regulated proteins involved in bioenergetics, immune system, active efflux and oxidative stress, often in the opposite way. Alterations of several proteins, notably several cytoskeletal ones, were still observed after the depuration period. These may reflect either the continuing chemical effect due to incomplete elimination or an onset of recovery processes in the mussel gills. Our study shows that exposure of adult mussels to sublethal TP concentration results in the bioconcentration of this biocide in the tissues and modulates the expression of several proteins that may intervene in important metabolic pathways.

A fugacity-based toxicokinetic model for narcotic organic chemicals in fish

A novel dynamic fugacity-based model is described, developed, and tested that simulates the uptake of narcotic organic chemicals in fish from water as occurs in aquatic bioconcentration and toxicity tests. The physiologically based toxicokinetic model treats the time course of chemical distribution in 4 compartments (tissue groups) in the fish, including the liver, in which biotransformation may occur. In addition to calculating bioconcentration and toxicokinetics, 5 possible toxic endpoints are defined corresponding to chemical concentration, fugacity, or activity reaching a critical value that causes 50% mortality. The mathematical description of multicompartment uptake is simplified by expressing the equations in the fugacity format. The model is parameterized and tested against reported empirical data for the bioconcentration of pentachloroethane in rainbow trout and for uptake and mortality from aquatic exposures to naphthalene and 1,2,4-trichlorobenzene in fathead minnows. Model performance is evaluated, and it is concluded that with suitable parameterization it has potential for application for assessment of both bioconcentration and toxicity expressed as median lethal concentrations, critical body residues, and chemical activity as a function of time to death.

Variation in predicted internal concentrations in relation to PBPK model complexity for rainbow trout

The present study is motivated by the increasing demand to consider internal partitioning into tissues instead of exposure concentrations for the environmental toxicity assessment. To this end, physiologically based pharmacokinetic (PBPK) models can be applied. We evaluated the variation in accuracy of PBPK model outcomes depending on tissue constituents modeled as sorptive phases and chemical distribution tendencies addressed by molecular descriptors. The model performance was examined using data from 150 experiments for 28 chemicals collected from US EPA databases. The simplest PBPK model is based on the "Kow-lipid content" approach as being traditional for environmental toxicology. The most elaborated one considers five biological sorptive phases (polar and non-polar lipids, water, albumin and the remaining proteins) and makes use of LSER (linear solvation energy relationship) parameters to describe the compound partitioning behavior. The "Kow-lipid content"-based PBPK model shows more than one order of magnitude difference in predicted and measured values for 37% of the studied exposure experiments while for the most elaborated model this happens only for 7%. It is shown that further improvements could be achieved by introducing corrections for metabolic biotransformation and compound transmission hindrance through a cellular membrane. The analysis of the interface distribution tendencies shows that polar tissue constituents, namely water, polar lipids and proteins, play an important role in the accumulation behavior of polar compounds with H-bond donating functional groups. For compounds without H-bond donating fragments preferable accumulation phases are storage lipids and water depending on compound polarity.

Cross-Species Extrapolation of Uptake and Disposition of Neutral Organic Chemicals in Fish Using a Multispecies Physiologically-Based Toxicokinetic Model Framework

The potential to bioconcentrate is generally considered to be an unwanted property of a substance. Consequently, chemical legislation, including the European REACH regulations, requires the chemical industry to provide bioconcentration data for chemicals that are produced or imported at volumes exceeding 100 tons per annum or if there is a concern that a substance is persistent, bioaccumulative, and toxic. For the filling of the existing data gap for chemicals produced or imported at levels that are below this stipulated volume, without the need for additional animal experiments, physiologically-based toxicokinetic (PBTK) models can be used to predict whole-body and tissue concentrations of neutral organic chemicals in fish. PBTK models have been developed for many different fish species with promising results. In this study, we developed PBTK models for zebrafish (Danio rerio) and roach (Rutilus rutilus) and combined them with existing models for rainbow trout (Onchorhynchus mykiss), lake trout (Salvelinus namaycush), and fathead minnow (Pimephales promelas). The resulting multispecies model framework allows for cross-species extrapolation of the bioaccumulative potential of neutral organic compounds. Predictions were compared with experimental data and were accurate for most substances. Our model can be used for probabilistic risk assessment of chemical bioaccumulation, with particular emphasis on cross-species evaluations.

Bioconcentration, metabolism and effects of diphenhydramine on behavioral and biochemical markers in crucian carp (Carassius auratus)

Diphenhydramine (DPH), an antihistamine used to alleviate human allergies, is widespread in aquatic environments. However, little is known about the biochemical and behavioral effects of DPH on non-target aquatic animals. In the present study, the tissue distribution, bioconcentration, metabolism, biochemical and behavioral effects were investigated in crucian carp (Carassius auratus) exposed to various concentrations of DPH (0.84, 4.23, 21.7 and 112 μg L(-1)) for 7d. DPH can accumulate in crucian carp, and high concentrations have been observed in the liver and brain with maximum bioconcentration factors of 148 and 81.6, respectively. A portion of the absorbed DPH was metabolized by the crucian carp to N-demethyl DPH and N,N-didemethyl DPH via N-demethylation. Direct fluorimetric assay was employed to assess metabolic activity, while oxidative stress and neurotransmission biomarkers were determined by Diagnostic Reagent Kits. DPH was found to increase hepatic 7-ethoxyresorufin O-deethylase activity in crucian carp with maximal induction of 119%. Concerning the oxidative stress status, DPH significantly inhibited superoxide dismutase (SOD, 37-58%) and glutathione S-transferase (GST, 43-65%) activities and led to a significant increase in malondialdehyde (MDA, 67-140%) levels and catalase (CAT, 38-143%) and glutathione peroxidase (GPx, 39-189%) activities in fish liver. Brain acetylcholinesterase activity was also induced in DPH-exposed crucian carp with maximal induction of 174%. In addition, shoaling was significantly enhanced, while swimming activity and feeding rates were markedly suppressed at DPH concentrations equal to or higher than 21.7 μg L(-1). Furthermore, significant correlations were found between oxidative stress biomarkers (SOD, CAT, GPx, GST and MDA) and behavioral parameters. Collectively, our results confirmed that DPH can accumulate and be metabolized in fish and exert a negative effect at different levels of biological organization.

Advancing In Vitro-In Vivo Extrapolations of Mechanism-Specific Toxicity Data Through Toxicokinetic Modeling

International legislation, such as the European REACH regulation (registration, evaluation, authorization, and restriction of chemicals), mandates the assessment of potential risks of an ever-growing number of chemicals to the environment and human health. Although this legislation is considered one of the most important investments in consumer safety ever, the downside is that the current testing strategies within REACH rely on extensive animal testing. To address the ethical conflicts arising from these increased testing requirements, decision-makers, such as the European Chemicals Agency (ECHA), are committed to Russel and Burch's 3R principle (i.e., reduction, replacement, refinement) by demanding that animal experiments should be substituted with appropriate alternatives whenever possible. A potential solution of this dilemma might be the application of in vitro bioassays to estimate toxic effects using cells or cellular components instead of whole organisms. Although such assays are particularly useful to assess potential mechanisms of toxic action, scientists require appropriate methods to extrapolate results from the in vitro level to the situation in vivo. Toxicokinetic models are a straightforward means of bridging this gap. The present chapter describes different available options for in vitro-in vivo extrapolation (IVIVE) of mechanism-specific effects focused on fish species and also reviews the implications of confounding factors during the conduction of in vitro bioassays and their influence on the optimal choice of different dose metrics.

Potential Biotechnological Strategies for the Cleanup of Heavy Metals and Metalloids

Global mechanization, urbanization, and various natural processes have led to the increased release of toxic compounds into the biosphere. These hazardous toxic pollutants include a variety of organic and inorganic compounds, which pose a serious threat to the ecosystem. The contamination of soil and water are the major environmental concerns in the present scenario. This leads to a greater need for remediation of contaminated soils and water with suitable approaches and mechanisms. The conventional remediation of contaminated sites commonly involves the physical removal of contaminants, and their disposition. Physical remediation strategies are expensive, non-specific and often make the soil unsuitable for agriculture and other uses by disturbing the microenvironment. Owing to these concerns, there has been increased interest in eco-friendly and sustainable approaches such as bioremediation, phytoremediation and rhizoremediation for the cleanup of contaminated sites. This review lays particular emphasis on biotechnological approaches and strategies for heavy metal and metalloid containment removal from the environment, highlighting the advances and implications of bioremediation and phytoremediation as well as their utilization in cleaning-up toxic pollutants from contaminated environments.

A physiologically based toxicokinetic (PBTK) model for moderately hydrophobic organic chemicals in the European eel (Anguilla anguilla)

The European eel (Anguilla anguilla) is a facultatively catadromous fish species with a complex life cycle. Its current population status is alarming: recruitment has decreased drastically since the 1980s and its stock is still considered to be outside safe biological limits. Although there is no consensus on the reasons for this situation, it is currently thought to have resulted from a combination of different stressors, including anthropogenic contaminants. To deepen our understanding of the processes leading to the accumulation of lipophilic organic contaminants in yellow eels (i.e. the feeding, continental growth stage), we developed a physiologically based toxicokinetic model using our own data and values from the literature. Such models can predict the uptake and distribution of water-borne organic chemicals in the whole fish and in different tissues at any time during exposure. The predictive power of the model was tested against experimental data for six chemicals with n-octanol-water partitioning coefficient (log Kow) values ranging from 2.13-4.29. Model performance was excellent, with a root mean squared error of 0.28 log units. This model has the potential to help identify suitable habitats for restocking under eel management plans.

Modeling bioaccumulation and biomagnification of nonylphenol and its ethoxylates in estuarine-marine food chains

There are several studies on bioaccumulation and biomagnification of nonylphenol (NP) and its ethoxylates (NPEOs), but their toxico-kinetic mechanisms remain unclear. In the present investigation, we explored the accumulation of NP and NPEOs in estuarine-marine food chains with a bioaccumulation model comprising five trophic levels. Using this model, we estimated uptake and elimination rate constants for NPEOs based on the organisms' weight and lipid content and the chemicals' Kow. Further, we calculated accumulation factors for NP and NPEOs, including biota-sediment accumulation factors (BSAF) and biomagnification factors (BMF), and compared these to independent field measurements collected in the Western Scheldt estuary in The Netherlands and field data reported in the literature. The estimated BSAF values for NP and total NPEOs were below 1 for all trophic levels. The estimated BMF values were around 1 for all trophic levels except for the highest level (carnivorous mammals and birds). For this trophic level, the estimated BMF value varied between 0.1 and 2.4, depending on the biotransformation capacity. For all trophic levels, except primary producers, the accumulation estimates that accounted for biotransformation of NPEOs into NP were closer to the field data than model estimates that did not include biotransformation, indicating that NP formation by biotransformation of NPEOs might occur in organisms.

Toxicology across scales: Cell population growth in vitro predicts reduced fish growth

Environmental risk assessment of chemicals is essential but often relies on ethically controversial and expensive methods. We show that tests using cell cultures, combined with modeling of toxicological effects, can replace tests with juvenile fish. Hundreds of thousands of fish at this developmental stage are annually used to assess the influence of chemicals on growth. Juveniles are more sensitive than adult fish, and their growth can affect their chances to survive and reproduce. Thus, to reduce the number of fish used for such tests, we propose a method that can quantitatively predict chemical impact on fish growth based on in vitro data. Our model predicts reduced fish growth in two fish species in excellent agreement with measured in vivo data of two pesticides. This promising step toward alternatives to fish toxicity testing is simple, inexpensive, and fast and only requires in vitro data for model calibration.

Tissue-specific bioaccumulation of human and veterinary antibiotics in bile, plasma, liver and muscle tissues of wild fish from a highly urbanized region

We investigated the bioaccumulation of antibiotics in bile, plasma, liver and muscle tissues of wild fish from four rivers in the Pearl River Delta region. In total, 12 antibiotics were present in at least one type of fish tissues from nine wild fish species in the four rivers. The mean values of log bioaccumulation factors (log BAFs) for the detected antibiotics in fish bile, plasma, liver, and muscle tissues were at the range of 2.06-4.08, 1.85-3.47, 1.41-3.51, and 0.48-2.70, respectively. As the digestion tissues, fish bile, plasma, and liver showed strong bioaccumulation ability for some antibiotics, indicating a different bioaccumulation pattern from hydrophobic organic contaminants. Human health risk assessment based on potential fish consumption indicates that these antibiotics do not appear to pose an appreciable risk to human health. To the best of our knowledge, this is first report of bioaccumulation patterns of antibiotics in wild fish bile and plasma.

Development and application of the adverse outcome pathway framework for understanding and predicting chronic toxicity: I. Challenges and research needs in ecotoxicology

To elucidate the effects of chemicals on populations of different species in the environment, efficient testing and modeling approaches are needed that consider multiple stressors and allow reliable extrapolation of responses across species. An adverse outcome pathway (AOP) is a concept that provides a framework for organizing knowledge about the progression of toxicity events across scales of biological organization that lead to adverse outcomes relevant for risk assessment. In this paper, we focus on exploring how the AOP concept can be used to guide research aimed at improving both our understanding of chronic toxicity, including delayed toxicity as well as epigenetic and transgenerational effects of chemicals, and our ability to predict adverse outcomes. A better understanding of the influence of subtle toxicity on individual and population fitness would support a broader integration of sublethal endpoints into risk assessment frameworks. Detailed mechanistic knowledge would facilitate the development of alternative testing methods as well as help prioritize higher tier toxicity testing. We argue that targeted development of AOPs supports both of these aspects by promoting the elucidation of molecular mechanisms and their contribution to relevant toxicity outcomes across biological scales. We further discuss information requirements and challenges in application of AOPs for chemical- and site-specific risk assessment and for extrapolation across species. We provide recommendations for potential extension of the AOP framework to incorporate information on exposure, toxicokinetics and situation-specific ecological contexts, and discuss common interfaces that can be employed to couple AOPs with computational modeling approaches and with evolutionary life history theory. The extended AOP framework can serve as a venue for integration of knowledge derived from various sources, including empirical data as well as molecular, quantitative and evolutionary-based models describing species responses to toxicants. This will allow a more efficient application of AOP knowledge for quantitative chemical- and site-specific risk assessment as well as for extrapolation across species in the future.

Relationships between exposure and dose in aquatic toxicity tests for organic chemicals

There is continuing debate about the merits of exposure-based toxicity metrics such as median lethal concentration (LC50) versus organism-based metrics such as critical body residue (CBR) as indicators of chemical toxicity to aquatic organisms. To demonstrate relationships and differences between these 2 metrics, the authors applied a simple one-compartment toxicokinetic mass-balance model for water-exposed fish for a series of hypothetical organic chemicals exhibiting baseline narcotic toxicity. The authors also considered the influence of several toxicity-modifying factors. The results showed that the results of standard toxicity tests, such as the LC50, are strongly influenced by several modifying factors, including chemical and organism characteristics such as hydrophobicity, body size, lipid content, metabolic biotransformation, and exposure durations. Consequently, reported LC50s may not represent consistent dose surrogates and may be inappropriate for comparing the relative toxicity of chemicals. For comparisons of toxicity between chemicals, it is preferable to employ a delivered dose metric, such as the CBR. Reproducible toxicity data for a specific combination of chemical, exposure conditions, and organism can be obtained only if the extent of approach to steady state is known. Suggestions are made for revisions in test protocols, including the use of models in advance of empirical testing, to improve the efficiency and effectiveness of tests and reduce the confounding influences of toxicity-modifying factors, especially exposure duration and metabolic biotransformation. This will assist in linking empirical measurements of LC50s and CBRs, 2 different but related indicators of aquatic toxicity, and thereby improve understanding of the large existing database of aquatic toxicity test results.

Physiologically-based toxicokinetic models help identifying the key factors affecting contaminant uptake during flood events

As a consequence of global climate change, we will be likely facing an increasing frequency and intensity of flood events. Thus, the ecotoxicological relevance of sediment re-suspension is of growing concern. It is vital to understand contaminant uptake from suspended sediments and relate it to effects in aquatic biota. Here we report on a computational study that utilizes a physiologically based toxicokinetic model to predict uptake, metabolism and excretion of sediment-borne pyrene in rainbow trout (Oncorhynchus mykiss). To this end, data from two experimental studies were compared with the model predictions: (a) batch re-suspension experiments with constant concentration of suspended particulate matter at two different temperatures (12 and 24°C), and (b) simulated flood events in an annular flume. The model predicted both the final concentrations and the kinetics of 1-hydroxypyrene secretion into the gall bladder of exposed rainbow trout well. We were able to show that exhaustive exercise during exposure in simulated flood events can lead to increased levels of biliary metabolites and identified cardiac output and effective respiratory volume as the two most important factors for contaminant uptake. The results of our study clearly demonstrate the relevance and the necessity to investigate uptake of contaminants from suspended sediments under realistic exposure scenarios.

QSARs for aquatic toxicity: celebrating, extending and displaying the pioneering contributions of Ferguson, Konemann and Veith

Significant advances were made in the development of quantitative structure-activity relationships (QSARs) relating molecular structure to aquatic toxicity by three studies over 30 years ago by Ferguson in 1939, Konemann in 1981, and Veith and colleagues in 1983. We revisit the original concepts and data from these studies and review these contributions from the bases of current perspectives on the hypothesized mechanism of baseline narcotic toxicity and the underlying thermodynamic and kinetic aspects. The relationships between LC50, octanol-water partition coefficient, aqueous solubility, chemical activity and chemical volume fraction in lipid phases are outlined including kinetic influences on measured toxicities. These relationships provide a compelling and plausible explanation of the success of these and other QSARs for aquatic toxicity. Suggestions are made for further advances in these QSARs to improve assessments of toxicity by baseline narcotic toxicity and selective modes of action, especially using emerging quantum chemical computational capabilities.

Measured and modeled toxicokinetics in cultured fish cells and application to in vitro-in vivo toxicity extrapolation

Effect concentrations in the toxicity assessment of chemicals with fish and fish cells are generally based on external exposure concentrations. External concentrations as dose metrics, may, however, hamper interpretation and extrapolation of toxicological effects because it is the internal concentration that gives rise to the biological effective dose. Thus, we need to understand the relationship between the external and internal concentrations of chemicals. The objectives of this study were to: (i) elucidate the time-course of the concentration of chemicals with a wide range of physicochemical properties in the compartments of an in vitro test system, (ii) derive a predictive model for toxicokinetics in the in vitro test system, (iii) test the hypothesis that internal effect concentrations in fish (in vivo) and fish cell lines (in vitro) correlate, and (iv) develop a quantitative in vitro to in vivo toxicity extrapolation method for fish acute toxicity. To achieve these goals, time-dependent amounts of organic chemicals were measured in medium, cells (RTgill-W1) and the plastic of exposure wells. Then, the relation between uptake, elimination rate constants, and log KOW was investigated for cells in order to develop a toxicokinetic model. This model was used to predict internal effect concentrations in cells, which were compared with internal effect concentrations in fish gills predicted by a Physiologically Based Toxicokinetic model. Our model could predict concentrations of non-volatile organic chemicals with log KOW between 0.5 and 7 in cells. The correlation of the log ratio of internal effect concentrations in fish gills and the fish gill cell line with the log KOW was significant (r>0.85, p = 0.0008, F-test). This ratio can be predicted from the log KOW of the chemical (77% of variance explained), comprising a promising model to predict lethal effects on fish based on in vitro data.

Understanding receptor-mediated effects in rainbow trout: in vitro-in vivo extrapolation using physiologically based toxicokinetic models

The European REACH regulation requires the use of animal experimentation to assess the risk of industrial chemicals. However, the 3R principle (reduction, replacement, refinement) demands the use of suitable alternative test methods. Many dossiers submitted for the authorization of chemicals have attempted to provide the required data without performing new experiments, relying heavily on in silico methods; in vitro assays were scarcely used. We propose a methodology that uses physiologically based toxicokinetic (PBTK) models to extrapolate in vitro data to the in vivo level. We collected experimental results for in vitro and in vivo ethoxyresorufin-O-deethylase and vitellogenin induction following chemical exposure and compared those results with model predictions. We found that the predictive power of aqueous chemical concentrations was limited; median effect concentrations (EC50s) based on internal concentrations in fish correlated better with in vitro EC50s. Our data show that in vitro assays could offer a substitute for fish studies when combined with PBTK models.

A physiologically based toxicokinetic model for the zebrafish Danio rerio

Zebrafish (Danio rerio) is a widely used model for toxicological studies, in particular those related to investigations on endocrine disruption. The development and regulatory use of in vivo and in vitro tests based on this species can be enhanced by toxicokinetic modeling. For this reason, we propose a physiologically based toxicokinetic (PBTK) model for zebrafish describing the uptake and disposition of organic chemicals. The model is based on literature data on zebrafish, other cyprinidae and other fish families, new experimental physiological information (volumes, lipids and water contents) obtained from zebrafish, and chemical-specific parameters predicted by generic models. The relevance of available models predicting the latter parameters was evaluated with respect to gill uptake and partition coefficients in zebrafish. This evaluation benefited from the fact that the influence of confounding factors such as body weight and temperature on ventilation rate was included in our model. The predictions for six chemicals (65 data points) yielded by our PBTK model were compared to available toxicokinetics data for zebrafish and 88% of them were within a factor of 5 of the corresponding experimental values. Sensitivity analysis highlighted that the 1-octanol/water partition coefficient, the metabolism rate, and all the parameters that enable the prediction of assimilation efficiency and partitioning of chemicals need to be precisely determined in order to allow an effective toxicokinetic modeling.

Highly time-variable exposure to chemicals--toward an assessment strategy

Organisms in the environment experience fluctuating, pulsed, or intermittent exposure to pollutants. Accounting for effects of such exposures is an important challenge for environmental risk assessment, particularly given the simplified design of standard ecotoxicity tests. Dynamic simulation using toxicokinetic-toxicodynamic (TK-TD) models describes the processes that link exposure with effects in an organism and provides a basis for extrapolation to a range of exposure scenarios. In so doing, TK-TD modeling makes the risk assessment more robust and aids use and interpretation of experimental data. Toxicokinetic-toxicodynamic models are well-developed for predicting survival of individual organisms and are increasingly applied to sublethal endpoints. In the latter case particularly, linkage to individual-based models (IBMs) allows extrapolation to population level as well as accounting for differences in effects of toxicant exposure at different stages in the life cycle. Extrapolation between species remains an important constraint because there is currently no systematic understanding of species traits that cause differences in the relevant processes. Toxicokinetic-toxicodynamic models allow interrogation of exposure profiles to determine intrinsic toxicity potential rather than using absolute maximum concentrations or time-weighted averages as surrogates. A decision scheme is proposed to guide selection of risk assessment approaches using dose extrapolation based on Haber's Law, TK-TD models, and/or IBMs depending on the nature of toxic effect and timing in relation to life history.

A list of fish species that are potentially exposed to pesticides in edge-of-field water bodies in the European Union--a first step towards identifying vulnerable representatives for risk assessment

Surrogate species are used in standard toxicity tests for the environmental risk assessment of chemicals. Test results are then extrapolated to the situation in the field, which is often associated with a large degree of uncertainty. Since a vulnerable species in the field is not only characterised by its intrinsic sensitivity to a stressor but also by its potential for exposure and its population resilience, the identification of focal species based on these three components of vulnerability is needed for a more ecologically relevant risk assessment. This study listed European fish species that are susceptible to pesticide exposure in the field and thus achieved the first step towards identifying focal species for the risk assessment of pesticides for fish in Europe. A step-wise filtering approach was applied to list freshwater fish species that are native to Europe and widespread in the European Union, which inhabit streams, ditches or ponds in agricultural landscapes and therefore, are at an elevated risk of being exposed to pesticides. Out of the 579 fish species occurring in European freshwater, 27 species met the filtering criteria. The resulting list was verified based on monitoring studies that were conducted in agricultural landscapes over the past 20 years. Focal fish species that can be used for a more ecologically relevant environmental risk assessment of pesticides in Europe can be identified from the produced list of species by further assessing their ecological (life history and dispersal characteristics) and intrinsic sensitivities.

A method to predict and understand fish survival under dynamic chemical stress using standard ecotoxicity data

The authors present a method to predict fish survival under exposure to fluctuating concentrations and repeated pulses of a chemical stressor. The method is based on toxicokinetic-toxicodynamic modeling using the general unified threshold model of survival (GUTS) and calibrated using raw data from standard fish acute toxicity tests. The model was validated by predicting fry survival in a fish early life stage test. Application of the model was demonstrated by using Forum for Co-ordination of Pesticide Fate Models and Their Use surface water (FOCUS-SW) exposure patterns as model input and predicting the survival of fish over 485 d. Exposure patterns were also multiplied by factors of five and 10 to achieve higher exposure concentrations for fish survival predictions. Furthermore, the authors quantified how far the exposure profiles were below the onset of mortality by finding the corresponding exposure multiplication factor for each scenario. The authors calculated organism recovery times as additional characteristic of toxicity as well as number of peaks, interval length between peaks, and mean duration as additional characteristics of the exposure pattern. The authors also calculated which of the exposure patterns had the smallest and largest inherent potential toxicity. Sensitivity of the model to parameter changes depends on the exposure pattern and differs between GUTS individual tolerance and GUTS stochastic death. Possible uses of the additional information gained from modeling to inform risk assessment are discussed.