Generic parameterization for a pharmacokinetic model to predict Cd concentrations in several tissues of different fish species

Dynamics of copper regulation in a marine clam Sinonovacula constricta at the organ level: Insight from a physiologically based pharmacokinetic model

Copper (Cu) is a common pollutant in estuaries and has received considerable attention worldwide. To gain an insight into the physiological mechanisms of waterborne Cu absorption, tissue distribution, storage, metabolism, and excretion in an estuarine razor clam Sinonovacula constricta, we developed a physiologically-based pharmacokinetic model based on prolonged Cu exposure with two exposure treatments. The tissues of S. constricta were divided into four parts: blood, digestive gland, gill, and other tissues. Our results showed that the waterborne Cu entered and exchanged with the gills and digestive gland, whereas digestive gland and other tissues were the main storage sites for Cu. Gills of S. constricta were able to maintain their Cu concentrations under both exposure treatments. Additionally, the gills exhibited a remarkable ability to remove Cu from water, with a transfer rate constant of 1.73 d-1 from the gills to water, while restricting its transfer from the blood with a transfer rate constant of 0.0131 d-1 from blood to gills. These results highlighted the crucial role of gills in regulating Cu levels in S. constricta as well as the detoxification and maintenance of metal homeostasis. Cu uptake rate constant in gill from waterborne was similar to that of digestive gland (0.294 vs. 0.364 L g-1 d-1), thus water entering the digestive tract was considered as another route of waterborne Cu absorption in bivalves. A significant amount of Cu in the blood was transferred to the digestive glands. These two factors explained the relatively higher Cu accumulation in the digestive glands than in other tissues in clams. The findings of this study enhanced our understanding of the homeostatic regulation and transportation mechanisms in marine bivalves.

Co-exposure effects of lead and TiO2 nanoparticles in primary kidney cell culture from the freshwater fish Hoplias malabaricus

This study evaluated the effects of Lead (Pb) and titanium dioxide nanoparticles (TiO₂ NPs) alone or in combination in anterior kidney macrophages of the freshwater fish Hoplias malabaricus, naïve or stimulated with 1ng.mL^-1 lipopolysaccharide (LPS). Pb (1×10^-5 to 1×10^-1mg.mL^-1) or TiO₂ NPs (1.5×10^-6 to 1.5×10^-2mg.mL^-1) reduced cell viability despite LPS stimulation, especially Pb 10^-1mg.mL^-1. In combination, lower concentrations of NPs intensified Pb-induced cell viability reduction while higher concentrations restored the cell viability independently of LPS stimulation. Basal and LPS- induced NO production was reduced by both TiO2 NPs and Pb isolated. The combination of both xenobiotics avoided this reduction of NO production by the isolated compounds at lower concentrations but the protective effect was lost as the concentrations increased. None xenobiotic increase DNA fragmentation. Therefore, at specific conditions, TiO₂ NPs may have a protective effect over Pb toxicity, may also provide additional toxicity at higher concentrations.

An Arsenic-76 radiotracer to study the routes of assimilation, hemolymph distribution, and tissue inventories in the bioindicator organism Pomacea canaliculata

The aim of this work was to study the absorption, distribution through the hemolymph, and bioaccumulation of arsenic by the freshwater Pomacea canaliculata using a short-lived tracer (⁷⁶As, t_1/2: 1.07 d) with high specific activity. Arsenic travels mainly dissolved in the plasma of the snail's hemolymph. This element is transferred from the hemolymph to the tissues (87%) 4 h after the inoculation of 50 μL of a 0.04 g/L of ⁷⁶As radiotracer solution, being the digestive gland, kidney, and head-foot the main places of arsenical inventories. Snails exhibited a rapid arsenic accumulation response in a wide range of concentrations (from 1 to 1000 μg/L) of the metalloid dissolved in water and in a concentration-dependent manner. Also, snails incorporated As from the digestive system when they received a single safe dose of ~2 μg of ⁷⁶As inoculated in a fish food pellet. The (semi) physiologically based toxicokinetic model developed in this study is based on anatomical and physiological parameters (blood flow, irrigation, tissue volume and other). Together, these findings make P. canaliculata an excellent sentinel organism to evaluate freshwater bodies naturally contaminated with As.

Physiologically Based Pharmacokinetic Model for the Biotransportation of Arsenic in Marine Medaka (Oryzias melastigma)

The toxicity of arsenic (As) targets specific tissues of organisms, while the biotransportation of As among the tissues of fish remains poorly understood. In the present study, radiotracer techniques followed by a physiologically based pharmacokinetic (PBPK) modeling were applied to simulate the biotransportation (absorption, distribution, and elimination) of ⁷³As(V) and biotransformation of As(V) in the marine medaka Oryzias melastigma after waterborne As exposure. Fish were simulated by a six-compartment model by assuming that blood was the intermediate exchange among different compartments (gill, intestine, liver, head, and carcass). Modeling suggested that intestine and gill were the uptake, exchange, as well as elimination sites of waterborne As, while carcass and head were the main storage sites. Intestine played a vital role in the metabolism of As(V) by biotransforming inorganic As into arsenobetaine (AsB), possibly because of the important role of gut microbiota. The correlation between the PBPK model constants and the As speciation (e.g., AsB %, inorganic As %, and methylated As %) indicated that AsB tended to be stored in the tissues rather than being depurated, while inorganic and methylated As were more easily transferred from tissues to the blood and eliminated. Modeling simulation coupling with biotransformation for the first time demonstrated that the fish intestine was the main metabolic site, and synthesis of AsB as mediated by the microbiota in the intestine contributed to the high As bioaccumulation in marine fish.

Applications of dynamic models in predicting the bioaccumulation, transport and toxicity of trace metals in aquatic organisms

This review evaluates the three dynamic models (biokinetic model: BK, physiologically based pharmacokinetic model: PBPK, and toxicokinetic-toxicodynamic model: TKTD) in our understanding of the key questions in metal ecotoxicology in aquatic systems, i.e., bioaccumulation, transport and toxicity. All the models rely on the first-order kinetics principle of metal uptake and elimination. The BK model basically treats organisms as a single compartment, and is both physiologically and geochemically based. With a good understanding of each kinetic parameter, bioaccumulation of metals in any aquatic organisms can be studied holistically and mechanistically. Modeling efforts are not merely restrained from the prediction of metal accumulation in the tissues, but instead provide the direction of the key processes that need to be addressed. PBPK is more physiologically based since it mainly addresses the transportation, transformation and distribution of metals in the organisms. It can be treated conceptually as a multi-compartmental kinetic model, whereas the physiology is driving the development of any good PBPK model which is no generic for aquatic animals and contaminants. There are now increasingly applications of the PBPK modeling specifically in metal studies, which reveal many important processes that are impossible to be teased out by direct experimental measurements without adequate modeling. TKTD models further focus on metal toxicity in addition to metal bioaccumulation. The TK part links exposure and bioaccumulation, while the TD part links bioaccumulation and toxic effects. The separation of TK and TD makes it possible to model processes, e.g., toxicity modification by environmental factors, interaction between different metals, at both the toxicokinetic and toxicodynamic levels. TKTD models provide a framework for making full use of metal toxicity data, and thus provide more information for environmental risk assessments. Overall, the three models reviewed here will continue to provide guiding principles in our further studies of metal bioaccumulation and toxicity in aquatic organisms.

Physiologically based toxicokinetic and toxicodynamic (PBTK-TD) modelling of Cd and Pb exposure in adult zebrafish Danio rerio: Accumulation and toxicity

Accurately predicting the accumulation and toxicity of metals in organisms is a challenging work in ecotoxicology. Here, we developed and validated a physiologically based toxicokinetic and toxicodynamic (PBTK-TD) model for adult zebrafish exposed to Cd and Pb. The model included the gill, liver, intestine, gonad, carcass, and brain, which were linked by blood circulation in the PBTK process and by dynamic relationships between the target organ concentrations and mortality in the TD process. Results showed that the PBTK sub-model can accurately describe and predict the uptake, distribution and disposition kinetics of Cd and Pb in zebrafish. The exchange rates and the accumulation of the metals in the organs were significantly different. For Cd, the highest exchange rate was between blood and liver, and the greatest accumulation of Cd occurred in the liver. For Pb, the greatest accumulation occurred in the gill. The TD sub-model further indicated that metal concentrations in the gill may effectively act as more suitable indicator of Cd and Pb toxic effect than whole body or other organs. The proposed PBTK-TD model is helpful to understanding the fundamental processes by which zebrafish regulate the uptake and disposition of metal and to quantitatively predicting metal toxicity.

Physiologically based toxicokinetics (PBTK) models for pharmaceuticals and personal care products in wild common carp (Cyprinus carpio)

Pharmaceutical and personal care products (PPCPs) are regarded as "pseudo-persistent" pollutants due to their being continuously loaded into the aquatic environment. Physiologically based toxicokinetics (PBTK) models that can quantitatively describe absorption, distribution, metabolism and excretion processes of chemicals in biota are of importance to predict internal exposure (e.g. doses at specific target tissues/organs) from external exposure concentrations. In this study, PBTK models with up to six compartments including brain, liver, kidney, gills, richly perfused tissues and poorly perfused tissues, were developed for predicting internal distribution of 10 PPCPs in wild common carp (Cyprinus carpio). The PBTK predicted concentrations were close to the measured ones, with deviations less than 1 log unit for most of PPCPs. Sensitivity analysis showed that various partition coefficients of the chemicals exerted significant influence on model outputs.

Development of a PBPK Model for Silver Accumulation in Chub Infected with Acanthocephalan Parasites

Simultaneous presence of metals and parasites in fish might lead to potential risks to human health. Parasites might influence metal accumulation and disturb detoxification in fish, thereby affecting biomarkers of fish responses as well as metal biomagnification in humans. It is, therefore, of importance to take into account parasite infection when investigating metal accumulation in fish. However, mechanisms of metal accumulation and distribution in fish-parasite systems are not integrated into current approaches. The present study proposes a new physiologically based pharmacokinetic model for mechanistic simulation of metal partitioning between intestinal parasites and their hosts. As a particular case, Ag accumulation in the system of chub Squalius cephalus and the acanthocephalan Pomphorhynchus tereticollis was investigated. As a novelty, fish cardiac output and organ-specific blood flow distribution were incorporated in our model. This approach distinguishes the current model from the ones developed previously. It also facilitates model extrapolation and application to varying conditions. In general, the model explained Ag accumulation in the system well, especially in chub gill, storage (including skin, muscle, and carcass), and liver. The highest concentration of Ag was found in the liver. The accumulation of Ag in the storage, liver, and gill compartments followed a similar pattern, i.e., increasing during the exposure and decreasing during the depuration. The model also generated this observed trend. However, the model had a weaker performance for simulating Ag accumulation in the intestine and the kidney. Silver accumulation in these organs was less evident with considerable variations.

Selenium induces the demethylation of mercury in marine fish

The antagonistic effect of selenium (Se) on mercury (Hg) toxicity has been known for decades. Earlier studies mainly focused on Hg-Se interaction based on biokinetics and bioaccumulation, but the influences of Se on in vivo biotransformation of methylmercury (MeHg) have not been well understood. We conducted a 42-day exposure study to investigate the dynamic changes of MeHg and its primary degradation product - inorganic mercury (IHg) - in different organs of black seabream (Acanthopagrus schlegeli) exposed to different dietary Se levels. A physiologically based pharmacokinetic (PBPK) model was then developed to describe the biotransformation and disposition of MeHg under the influence of Se. Our results demonstrated that Se significantly increased the transformation from MeHg into IHg, thereby decreasing the accumulation of MeHg. The simulation further showed that the intestine was the major site for demethylation, with an estimated rate 1.5-fold higher in high Se treatment than in low Se treatment. However, the hepatic demethylation rate was extremely low and comparable between the two treatments (0.012-0.015 d^-1). These results strongly suggested that the intestine instead of the commonly assumed liver was the major site for Hg-Se interaction. Furthermore, Se did not show significant influences on the distribution and elimination of MeHg, but promoted the uptake and elimination of the generated IHg from demethylation. Therefore, Se-induced demethylation especially in the intestine played an important role in mitigating the MeHg accumulation. This study provided new sight to elucidate the Hg-Se interaction in fish.

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.

Development and Validation of a Biodynamic Model for Mechanistically Predicting Metal Accumulation in Fish-Parasite Systems

Because of different reported effects of parasitism on the accumulation of metals in fish, it is important to consider parasites while interpreting bioaccumulation data from biomonitoring programmes. Accordingly, the first step is to take parasitism into consideration when simulating metal bioaccumulation in the fish host under laboratory conditions. In the present study, the accumulation of metals in fish-parasite systems was simulated by a one-compartment toxicokinetic model and compared to uninfected conspecifics. As such, metal accumulation in fish was assumed to result from a balance of different uptake and loss processes depending on the infection status. The uptake by parasites was considered an efflux from the fish host, similar to elimination. Physiological rate constants for the uninfected fish were parameterised based on the covalent index and the species weight while the parameterisation for the infected fish was carried out based on the reported effects of parasites on the uptake kinetics of the fish host. The model was then validated for the system of the chub Squalius cephalus and the acanthocephalan Pomphorhynchus tereticollis following 36-day exposure to waterborne Pb. The dissolved concentration of Pb in the exposure tank water fluctuated during the exposure, ranging from 40 to 120 μg/L. Generally, the present study shows that the one-compartment model can be an effective method for simulating the accumulation of metals in fish, taking into account effects of parasitism. In particular, the predicted concentrations of Cu, Fe, Zn, and Pb in the uninfected chub as well as in the infected chub and the acanthocephalans were within one order of magnitude of the measurements. The variation in the absorption efficiency and the elimination rate constant of the uninfected chub resulted in variations of about one order of magnitude in the predicted concentrations of Pb. Inclusion of further assumptions for simulating metal accumulation in the infected chub led to variations of around two orders of magnitude in the predictions. Therefore, further research is required to reduce uncertainty while characterising and parameterising the model for infected fish.

Pollutant levels in discarded fish species by Spanish trawlers operating in the Great Sole Bank and the Atlantic coast of the Iberian Peninsula

Organic and inorganic pollutant levels were determined for the most discarded species from trawlers operating in Great Sole and Spanish coastal fishing grounds. Results for heavy metals indicated that Cd can reach values higher than legal limits for some species and tissues, while Hg and Pb concentrations are below established values. No significant variation was noticed with fishing grounds, but both season influences in the case of Pb and interspecies variation for Hg and Cd have been detected. Valorization recommendations could be therefore established according to the levels found in the different species.

Modulatory effect of nano TiO₂on Pb in Hoplias malabaricus trophically exposed

This study investigated the hepatic and neural effects of TiO₂ nanoparticle and Pb in Hoplias malabaricus trophically exposed. The alanine transaminase activity was altered at the high dose of exposed group to Pb and at the lowest doses of co-exposed groups. It may reflect the hepatic effects of TiO₂ on Pb toxicity, but the aspatate transaminase activity was not altered. The decreased injury index observed at the highest dose of co-exposed group compared to TiO₂ may be related to the increased energy demand and can explain the more pronounced toxic effects observed in this group. The liver authomethallography revealed the metals presence at high dose groups. Serotonin concentration increased at the Pb lowest dose and at the highest dose of co-exposed group compare to control. Most importantly, when associated the contaminants were able to interact and altered some biomarkers. However, further studies, about action mechanisms of this co-exposure are needed.

Parameter uncertainty in modeling bioaccumulation factors of fish

We quantified the uncertainty due to biota-related parameters in estimated bioaccumulation factors (BAFs) of persistent organic pollutants for fish through Monte Carlo simulations. For this purpose, the bioaccumulation model OMEGA (Optimal Modeling for EcotoxicoloGical Applications) was parameterized based on data from the existing literature, analysis of allometric data, and maximum likelihood estimation. Lipid contents, fractions of food assimilated, the allometric rate exponent, normalized food intakes, respiration and growth dilution rates, and partial mass transfer resistances in water and lipid layers were included as uncertain parameters. The uncertainty in partial resistances was particularly important in the estimation of the rate constants for chemical intake from water by fish. Uncertainties in the fractions of food assimilated and partial water layer resistances from and to food were particularly important in the estimation of the rate constants of chemical intake from food. The uncertainty in the model outcomes for the bioaccumulation factors for fish was a factor of 10 (ratio of 95th and fifth percentile estimates), which was mainly caused by the uncertainty in the lipid fraction. For chemicals with a K(OW) of 10(3) to 10(6), the uncertainty in the lipid contents of fish accounted for more than 50% of the uncertainty in the estimated bioaccumulation factor. For chemicals with a high K(OW) (10(7) and higher), the fractions of food assimilated and partial resistances also contributed to uncertainty in the estimated bioaccumulation factor (up to 60%). A case study showed that uncertainty in estimated BAF for nonpersistent substances can be dominated by uncertainty in the rate constants for metabolic transformation.