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

Screening of pesticide distributions in foods of animal origin: a matrix-based approach for biotransfer factor modeling of grazing mammals

Pesticides are detected in a wide array of foods of animal origin, posing a threat to food safety and human health. Here, to facilitate the management of pesticide residues in livestock products, we proposed a screening model for simulating pesticide biotransfer factors (BTFs, ratio of steady-state pesticide concentration between animal tissues and feed) in mammalian bodies. The proposed model was developed based on simple matrix algebra using first-order kinetics. Simulation of over 700 pesticides in common food products derived from cattle and sheep indicated that pesticide biotransfer is a balancing process between uptake efficiency from the gut lumen and the uptake-elimination ratio (ratio of the overall uptake and elimination rate constants) in individual tissues. Furthermore, we parameterized the developed BTF model using the octanol-water partition coefficient. The simulated pesticide BTFs could be categorized into three lipophilicity phases, namely, the lipophilicity-boosting (log K_OW < 4), lipophilicity-balancing (4 ≤ log K_OW ≤ 8), and lipophilicity-limiting (log K_OW > 8) phases, and the simulated general trends of BTFs were consistent with the experimental data. Although the metabolic rate constants of pesticides in mammalian bodies warrant further evaluation, the model proposed here can assist in the risk assessment and regulatory management of pesticide residues in foods of animal origin.

A perspective on the role of fugacity and activity for evaluating the PBT properties of organic chemicals and providing a multi-media synoptic indicator of environmental contamination

The hazard and risk from organic chemicals present in the environment are routinely evaluated using P (persistence), B (bioaccumulation) and T (toxicity) criteria. We present a perspective on how models based on the equilibrium criteria of fugacity and chemical activity can contribute to all three evaluations, thus providing a consistent and seamless assessment process. Persistence and its closely related, but fundamentally different property, chemical residence time can be determined from degradation half-lives (typically obtained from monitoring data) and multi-media mass balance models describing chemical fate and transport in real or evaluative environments. Bioaccumulation is driven by equilibrium partitioning processes that can be estimated from fugacity models treating uptake by respiration and diet in single organisms and in food webs, most commonly for aquatic systems for which confirmatory monitoring data can be obtained. Biomagnification is readily evaluated both experimentally and as a fugacity (or activity) increase from prey to predator. Toxicokinetic evaluations of chemical uptake and disposition in a variety of organisms, including humans, can be obtained using fugacity- and physiologically-based pharmacokinetic models. Toxicodynamic processes of chemical interactions with organisms leading to adverse toxic outcomes are less obviously amenable to fugacity assessment with the notable exception of baseline toxicity (narcosis), the most common mode of action of commercial chemicals. It is shown that a range of fugacities can be defined and correlated for specific species thus enabling direct comparison of fugacities estimated or monitored in the environment with fugacities resulting in baseline toxicity. Finally, it is shown that fugacity and activity can serve as multi-media metrics of chemical contamination in ecosystems, thus enabling identification of specific regions and species at greatest risk in a pictorial format that is readily understandable to scientific and regulatory communities and to the general public.

The forgotten role of toxicodynamics: How habitat quality alters the mite, Oppia nitens, susceptibility to zinc, independent of toxicokinetics

Soil habitat quality is thought to influence metal toxicity via changes in speciation and thereby toxicokinetics. Here, we assessed the toxicokinetic and toxicodynamic effects of habitat quality on mite, Oppia nitens when exposed to zinc (Zn) contaminated soils. Forty-seven soils were ranked into three habitat qualities; high, medium, and low based on biological reproduction of Folsomia candida, Enchytraeus crypticus, and Elymus lanceolatus. From the 47 soils, eighteen soils (comprising of six soils from each habitat quality) were randomly selected and dosed with field relevant concentrations of Zn. Mite survival and reproduction were assessed after 28 days. Total Zn, bioaccessible Zn, Zn bioavailability, Zn body burden, lactate dehydrogenase activity (LDH) and glucose-6-phosphate dehydrogenase (G6PDH) activities of the mites were determined. Zinc toxicity and potency were much less in the high compared to low quality soils and the mites in the high habitat quality soils tolerated higher zinc body burdens (2040 ± 130 μg/g b.w) than the lower habitat quality (1180 ± 310 μg/g b.w). Lower LDH activity (20 ± 2 μU mg^-1) in the high quality soils compared to lower quality soils (50 ± 8 μU mg^-1) suggested that there was less stress in the high habitat quality mites. Despite changes in speciation across habitat qualities, bioavailability of zinc was similar (∼20%) irrespective of habitat quality. Our results suggest that the influence of soil properties on survival is modulated by toxicodynamics rather than toxicokinetics. Restoring habitat quality may be more important for soil invertebrate protection than metal concentration at contaminated sites.

Bioconcentration, bioaccumulation, biomagnification and trophic magnification: a modelling perspective

We present a modelling perspective on quantifying metrics of bio-uptake of organic chemicals in fish. The models can be in concentration, partition ratio, rate constant (CKk) format or fugacity, Z and D value (fZD) format that are shown to be exactly equivalent, each having it merits. For most purposes a simple, parameter-parsimonious one compartment steady-state model containing some 13 parameters is adequate for obtaining an appreciation of the uptake equilibria and kinetics for scientific and regulatory purposes. Such a model is first applied to the bioaccumulation of a series of hypothetical, non-biotransforming chemicals with log K_OW (octanol-water partition ratio) values of 4 to 8 in 10 g fish ranging in lipid contents to deduce wet-weight and lipid normalized concentrations, bioaccumulation and biomagnification factors. The sensitivity of biomagnification factors to relative lipid contents is discussed. Second, a hypothetical 5 species linear food chain is simulated to evaluate trophic magnification factors (TMFs) showing the critical roles of K_OW and biotransformation rate. It is shown that lipid normalization of concentrations is most insightful for less hydrophobic chemicals (log K_OW < 5) when bio-uptake is largely controlled by respiratory intake and equilibrium (equi-fugacity) is approached. For more hydrophobic chemicals when dietary uptake kinetics dominate, wet weight concentrations and BMFs are more insightful. Finally, a preferred strategy is proposed to advance the science of bioaccumulation using a combination of well-designed ecosystem monitoring, laboratory determinations and modelling to confirm that the perceived state of the science contained in the models is consistent with observations.

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.