Predicting the Metal Mixture Toxicity with a Toxicokinetic-Toxicodynamic Model Considering the Time-Dependent Adverse Outcome Pathways

Chemicals mainly exist in ecosystems as mixtures, and understanding and predicting their effects are major challenges in ecotoxicology. While the adverse outcome pathway (AOP) and toxicokinetic-toxicodynamic (TK-TD) models show promise as mechanistic approaches in chemical risk assessment, there is still a lack of methodology to incorporate the AOP into a TK-TD model. Here, we describe a novel approach that integrates the AOP and TK-TD models to predict mixture toxicity using metal mixtures (specifically Cd-Cu) as a case study. We preliminarily constructed an AOP of the metal mixture through temporal transcriptome analysis together with confirmatory bioassays. The AOP revealed that prolonged exposure time activated more key events and adverse outcomes, indicating different modes of action over time. We selected a potential key event as a proxy for damage and used it as a measurable parameter to replace the theoretical parameter (scaled damage) in the TK-TD model. This refined model, which connects molecular responses to organism outcomes, effectively predicts Cd-Cu mixture toxicity over time and can be extended to other metal mixtures and even multicomponent mixtures. Overall, our results contribute to a better understanding of metal mixture toxicity and provide insights for integrating the AOP and TK-TD models to improve risk assessment for chemical mixtures.

Investigating the Impact of Humic Acid on Copper Accumulation in Sinonovacula constricta Using a Toxicokinetic-Toxicodynamic Model

In aquatic ecosystems, the interaction between heavy metals and dissolved organic carbon (DOC) plays a pivotal role in modifying the bioavailability of these metals. This study, employing a toxicokinetic-toxicodynamic model, delves into the interactive effects of humic acid (HA), a significant component of DOC, on the bioaccumulation and toxicity of copper (Cu) in the estuarine economic bivalve Sinonovacula constricta. Utilizing the stable isotope 65Cu as a tracer, we evaluated Cu uptake in S. constricta under varied DOC concentrations in a controlled laboratory setting. Our findings reveal that at DOC concentrations below 3.05 mg L-1, the bioavailability of Cu is reduced due to shifts in the speciation distribution of Cu, resulting in decreased bioaccumulation within S. constricta. Conversely, at DOC levels exceeding 3.05 mg L-1, the formation of colloidal Cu-HA complexes allows its entry into the bivalves' digestive system. Moreover, toxicity assays demonstrate an increase in S. constricta survival rates with higher DOC concentrations, suggesting a protective effect of DOC against Cu toxicity. The integration of accumulation and toxicity data infers that Cu-HA complexes, when ingested via the digestive tract, exhibit lower toxicity compared to Cu directly assimilated from the water phase. These findings emphasize the need to consider environmental DOC levels in assessing Cu pollution risks and provide insights for managing heavy metal toxicity in estuarine aquaculture.

Application of biotic ligand and toxicokinetic-toxicodynamic modeling to predict the accumulation and toxicity of metal mixtures to zebrafish larvae

Predicting the accumulation and toxicity of mixtures of metals to aquatic organisms is a key challenge in ecotoxicological studies. In this study, the accumulation and toxicity of mixed essential (Cu) and nonessential (Cd and Pb) metals in zebrafish larvae exposed to a binary mixture of these elements at environmentally relevant concentrations were predicted using a refined toxicokinetic (TK)-toxicodynamic (TD) model aided with biotic ligand model (BLM) and toxic equivalent factor (TEF) approach. Competitive inhibition and non-competitive interaction/inhibition were observed in bio-uptake. Both Pb and Cd behaved as competitive inhibitors of Cu uptake at high Cu concentrations (>0.1 μM). By contrast, Cu uptake was independent of Cd or Pb when the Cu concentrations were below 10(-7) M. Furthermore, low concentrations of Cu had an adiaphorous effect on Cd or Pb uptake. Cd uptake was inhibited by Pb, and the Pb uptake rates consistently decreased in the presence of Cd. The accumulation processes of Cd-Pb, Cu-Cd, and Cu-Pb were accurately predicted by the BLM-aided TK models. The traditional TD model could successfully predict the toxicity of Cd-Pb mixtures, but not those of Cu-Cd or Cu-Pb mixtures. The revised TD model, which considered the possible different killing rates (Kk) above or below the threshold, offered better prediction for the toxicity of Cu-Cd or Cu-Pb mixtures. The overall findings may be of key significance in understanding and predicting metal uptake, accumulation, and toxicity in binary or multiple metal exposure scenarios.