Combining biomarkers of susceptibility with toxicodynamic models: a novel approach for assessing modifying factors.

Biomarkers of susceptibility are indicative of an individual's capacity to react to particular exposures, whereas toxicodynamic models elucidate the correlation between exposure and response. In this article, a novel methodology is presented for the evaluation of modifying factors in the field of toxicology. The proposed approach involves the integration of biomarkers of susceptibility with toxicodynamic models. Through the integration of these two methodologies, scientists are able to gain a more comprehensive understanding of the impact of modifying factors, such as genetic polymorphisms or epigenetic profiles, on an individual's reaction to toxic substances. This methodology has the potential to facilitate a more thorough evaluation of the hazards linked to the contact with combinations of chemicals and the cumulative effects of such exposures. The utilization of biomarkers in the evaluation of exposure for risk assessment is progressively incorporating the examination of susceptibility factors alongside exposure factors. This may involve the identification of a particular genetic polymorphism for a metabolic enzyme. The integration of ecotoxicological tests with models is crucial for achieving a comprehensive assessment. This approach exhibits the potential to enhance our comprehension of disease causation and facilitate the identification of populations that may exhibit an elevated susceptibility to disease.

Systemic effects of settleable atmospheric particulate matter (SePM) on swamp ghost crab Ucides cordatus.

Metallurgical activities are a significant source of settleable atmospheric particulate matter (SePM). The material is exposed to wind action, leading to its deposition throughout terrestrial and aquatic ecosystems, thus promoting contamination by metals and metalloids. However, knowledge of the impacts on biota is scarce. In aquatic coastal zones, evaluating hemolymph in invertebrates makes it possible to have insights into the pre-pathogenic effects and health status of organisms. Our study aimed to evaluate bioaccumulation and the sublethal effects of SePM on the mangrove crab Ucides cordatus by assessing biomarkers of cito-genotoxicity in the hemolymph. Organisms underwent a 30-day experiment with four treatments: control; 0.01 g.L-1, 0.1 g.L-1, 1 g.L-1 of SePM, with hemolymph sampled at 2, 7, 15, and 30 days of exposure to assess lipid peroxidation (LPO), DNA damage (strand break), cholinesterase (ChE) and lysosomal membrane stability (LMS). The results revealed metals' bioaccumulation in soft tissues (Al, Fe+, Fe++, Cu, Zr, Nb) and dose-time-dependent responses for LPO, DNA strand break, ChE, and LMS. Significant correlation was found between LPO and Cu (tissue), reduced LMS and Al and Fe (tissue), and Cu, Zn, Ag, and Bi in water. Hemolymph was related to the toxicokinetic and toxicodynamic of metals and metalloids from SePM in Ucides cordatus. New toxicological evidence was obtained to shed light on the impacts of SePM on the ecological status of coastal zones.

Modeling temperature-dependent life-cycle toxicity of thiamethoxam in Chironomus riparius using a DEB-TKTD model.

The neonicotinoid insecticide thiamethoxam (TMX) is widely used to protect crops against insect pests. Despite some desirable properties such as its low toxicity to birds and mammals, concerns have been raised about its toxicity to non-target arthropods, including freshwater insects like chironomids. Whereas multiple studies have investigated chronic effects of neonicotinoids in chironomid larvae at standardized laboratory conditions, a better understanding of their chronic toxicity under variable temperatures and exposure is needed for coherent extrapolation from the laboratory to the field. Here, we developed a quantitative mechanistic effect model for Chironomus riparius, to simulate the species' life history under dynamic temperatures and exposure concentrations of TMX. Laboratory experiments at four different temperatures (12, 15, 20, 23 °C) and TMX concentrations between 4 and 51 µg/L were used to calibrate the model. Observed concentration-dependent effects of TMX in C. riparius included slower growth, later emergence, and higher mortality rates with increasing concentrations. Furthermore, besides a typical accelerating effect on the organisms' growth and development, higher temperatures further increased the effects associated with TMX. With some data-informed modeling decisions, most prominently the inclusion of a size dependence that makes larger animals more sensitive to TMX, the model was parametrized to convincingly reproduce the data. Experiments at both a constant (20 °C) and a dynamically increasing temperature (15-23 °C) with pulsed exposure were used to validate the model. Finally, the model was used to simulate realistic exposure conditions using two reference exposure scenarios measured in Missouri and Nebraska, utilizing a moving time window (MTW) and either a constant temperature (20 °C) or the measured temperature profiles belonging to each respective scenario. Minimum exposure multiplication factors leading to a 10% effect (EP10) in the survival at pupation, i.e., the most sensitive endpoint found in this study, were 25.67 and 21.87 for the Missouri scenario and 38.58 and 44.64 for the Nebraska scenario, when using the respective temperature assumptions. While the results illustrate that the use of real temperature scenarios does not systematically modify the EPx in the same direction (making it either more or less conservative when used as a risk indicator), the advantage of this approach is that it increases the realism and thus reduces the uncertainty associated with the model predictions.

Toxicodynamic of combined mycotoxins: MicroRNAs and acute-phase proteins as diagnostic biomarkers.

Mycotoxins, ubiquitous contaminants in food, present a global threat to human health and well-being. Mitigation efforts, such as the implementation of sound agricultural practices, thorough food processing, and the advancement of mycotoxin control technologies, have been instrumental in reducing mycotoxin exposure and associated toxicity. To comprehensively assess mycotoxins and their toxicodynamic implications, the deployment of effective and predictive strategies is imperative. Understanding the manner of action, transformation, and cumulative toxic effects of mycotoxins, moreover, their interactions with food matrices can be gleaned through gene expression and transcriptome analyses at cellular and molecular levels. MicroRNAs (miRNAs) govern the expression of target genes and enzymes that play pivotal roles in physiological, pathological, and toxicological responses, whereas acute phase proteins (APPs) exert regulatory control over the metabolism of therapeutic agents, both endogenously and posttranscriptionally. Consequently, this review aims to consolidate current knowledge concerning the regulatory role of miRNAs in the initiation of toxicological pathways by mycotoxins and explores the potential of APPs as biomarkers following mycotoxin exposure. The findings of this research highlight the potential utility of miRNAs and APPs as indicators for the detection and management of mycotoxins in food through biological processes. These markers offer promising avenues for enhancing the safety and quality of food products.

Comparing Sensitivity of Different Bee Species to Pesticides: A TKTD modeling approach.

Risk assessment for bees is mainly based on data for honey bees; however, risk assessment is intended to protect all bee species. This raises the question of whether data for honey bees are a good proxy for other bee species. This issue is not new and has resulted in several publications in which the sensitivity of bee species is compared based on the values of the 48-h median lethal dose (LD50) from acute test results. When this approach is used, observed differences in sensitivity may result both from differences in kinetics and from inherent differences in species sensitivity. In addition, the physiology of the bee, like its overall size, the size of the honey stomach (for acute oral tests), and the physical appearance (for acute contact tests) also influences the sensitivity of the bee. The recently introduced Toxicokinetic-Toxicodynamic (TKTD) model that was developed for the interpretation of honey bee tests (Bee General Uniform Threshold Model for Survival [BeeGUTS]) could integrate the results of acute oral tests, acute contact tests, and chronic tests within one consistent framework. We show that the BeeGUTS model can be calibrated and validated for other bee species and also that the honey bee is among the more sensitive bee species. In addition, we found that differences in sensitivity between species are smaller than previously published comparisons based on 48-h LD50 values. The time-dependency of the LD50 and the specifics of the bee physiology are the main causes of the wider variation found in the published literature. Environ Toxicol Chem 2024;43:1431-1441. © 2024 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Mechanistic analysis of the sub chronic toxicity of La and Gd in Daphnia magna based on TKTD modelling and synchrotron X-ray fluorescence imaging.

The release of lanthanides (Ln) into the environment has increased in recent decades due to their expanding applications in society. Studying their toxicity in aquatic ecosystems is urgent and challenging, with contradictory evidence presented in the literature. This study compared the biodistribution of La and Gd in Daphnia magna exposed to sub-chronic conditions and developed the first Toxicokinetic-Toxicodynamic (TKTD) model for these lanthanides with this model crustacean. D. magna were initially exposed for 7 days to concentrations close to the LC50 of La (2.10 mg L-1) and Gd (1.70 mg L-1). After exposure, half of the live daphnids were introduced in a clean media to allow depuration over 24 h, while the other organisms were directly prepared for synchrotron imaging measurements. Synchrotron X-ray fluorescence analysis revealed that metal distribution in the organisms was similar for both La and Gd, predominantly localized in the intestinal tract, even after the depuration process. These results indicate that ingested metal can adversely affect organisms under sub-chronic exposure conditions, highlighting the importance of using nominal concentrations as a more suitable indicator of metal bioavailability for risk assessment. The General Unified Threshold Model of Survival (GUTS) TKTD framework, in its reduced form (GUTS-RED), was developed for La and Gd using dissolved and nominal concentrations. D. magna were exposed for 7 days to concentrations from 0.5 to 5 mg L-1 of La or Gd and mortality monitored daily. The mechanistic model revealed a faster toxicokinetics for La than Gd and a higher toxicity for Gd than La in the organism. This study confirmed, despite similar chemical properties, the variation in both toxicity and toxicokinetics between these two metals.

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.

Pharmacokinetic-toxicodynamic Modeling to Elucidate the Involvement of Dorsal Root Ganglion Neuron in Oxaliplatin-induced Peripheral Neuropathy.

BACKGROUND/AIM: Oxaliplatin (L-OHP)-induced peripheral neuropathy (OIPN) limits L-OHP dosage due to nerve cell damage in the dorsal root ganglion (DRG) caused by platinum (Pt). Despite various recommended approaches for OIPN management, no effective approach has been established. The aim of this study was to evaluate Pt distribution into DRG after repeat administrations of L-OHP in rats and to develop a pharmacokinetic-toxicodynamic (PK-TD) model using Pt concentrations in DRG to predict neuropathy severity. MATERIALS AND METHODS: Male Wistar rats were administered L-OHP (3, 5, or 8 mg/kg i.v.) once weekly. Blood and DRG samples were collected following L-OHP administration. For toxicodynamic (TD) study, OIPN was evaluated using the von Frey test. Plasma and DRG Pt concentrations and thresholds values in von Frey test were used for PK-TD modeling using Phoenix WinNonlin® version 8.3 software. RESULTS: Pt concentration in the DRG increased with repeated administration of L-OHP in a dose-dependent manner, indicating Pt accumulation in DRG following multiple administrations. The PK-TD model, consisting of an indirect response model and a transit compartment model with the DRG compartment, adequately described the temporal changes in OIPN with reliable TD parameters (≤36.4% with coefficient of variation). The maximum drug inhibition model could be employed to establish a quantitative correlation between the Pt content present in the DRG and the toxic potency resulting in OIPN. CONCLUSION: The utility of the PK-TD model for predicting neuropathy outcomes was established, suggesting that models composed of the DRG compartment contribute to determining an optimal dosing strategy for reducing OIPN.

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.

Recalling the reported toxicity assessment of deoxynivalenol, mitigating strategies and its toxicity mechanisms: Comprehensive review.

Mycotoxins frequently contaminate a variety of food items, posing significant concerns for both food safety and public health. The adverse consequences linked to poisoning from these substances encompass symptoms such as vomiting, loss of appetite, diarrhea, the potential for cancer development, impairments to the immune system, disruptions in neuroendocrine function, genetic damage, and, in severe cases, fatality. The deoxynivalenol (DON) raises significant concerns for both food safety and human health, particularly due to its potential harm to vital organs in the body. It is one of the most prevalent fungal contaminants found in edible items used by humans and animals globally. The presence of harmful mycotoxins, including DON, in food has caused widespread worry. Altered versions of DON have arisen as possible risks to the environment and well-being, as they exhibit a greater propensity to revert back to the original mycotoxins. This can result in the buildup of mycotoxins in both animals and humans, underscoring the pressing requirement for additional investigation into the adverse consequences of these modified mycotoxins. Furthermore, due to the lack of sufficient safety data, accurately evaluating the risk posed by modified mycotoxins remains challenging. Our review study delves into conjugated forms of DON, exploring its structure, toxicity, control strategies, and a novel animal model for assessing its toxicity. Various toxicities, such as acute, sub-acute, chronic, and cellular, are proposed as potential mechanisms contributing to the toxicity of conjugated forms of DON. Additionally, the study offers an overview of DON's toxicity mechanisms and discusses its widespread presence worldwide. A thorough exploration of the health risk evaluation associated with conjugated form of DON is also provided in this discussion.

Extrapolating Metal (Cu, Ni, Zn) Toxicity from Individuals to Populations Across Daphnia Species Using Mechanistic Models: The Roles of Uncertainty Propagation and Combined Physiological Modes of Action.

Mechanistic effect modeling is a promising tool to improve the ecological realism of environmental risk assessment. An open question for the mechanistic modeling of metal toxicity is whether the same physiological mode of action (PMoA) could be assumed for closely related species. The implications of various modeling choices, such as the use of parameter point estimates and assumption of simplistic toxicodynamic models, are largely unexplored. We conducted life-table experiments with Daphnia longispina, Daphnia magna, and Daphnia pulex exposed to the single metals Cu, Ni, and Zn, and calibrated toxicokinetic-toxicodynamic (TKTD) models based on dynamic energy budget theory. We developed TKTD models with single and combined PMoAs to compare their goodness-of-fit and predicted population-level sensitivity. We identified the PMoA reproduction efficiency as most probable in all species for Ni and Zn, but not for Cu, and found that combined-PMoA models predicted higher population-level sensitivity than single-PMoA models, which was related to the predicted individual-level sensitivity, rather than to mechanistic differences between models. Using point estimates of parameters, instead of sampling from the probability distributions of parameters, could also lead to differences in the predicted population-level sensitivity. According to model predictions, apical chronic endpoints (cumulative reproduction, survival) are conservative for single-metal population effects across metals and species. We conclude that the assumption of an identical PMoA for different species of Daphnia could be justified for Ni and Zn, but not for Cu. Single-PMoA models are more appropriate than combined-PMoA models from a model selection perspective, but propagation of the associated uncertainty should be considered. More accurate predictions of effects at low concentrations may nevertheless motivate the use of combined-PMoA models. Environ Toxicol Chem 2024;43:338-358. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Development of a mechanistic model for analyzing avian reproduction data for pesticide risk assessment.

Mechanistic effect models are increasingly recommended as tools for refining evaluations of risk from exposure to pesticides. In the context of bird and mammal risk assessments, DEB-TKTD models have been recommended for characterizing sublethal effects at lower tiers. However, there are currently no such models. Currently, chronic, multi-generational studies are performed to characterize potential effects of pesticides on avian reproduction, but it is has not been established to what extent results from these studies can inform effect models. Here, a standard Dynamic Energy Budget (DEB) model was extended to account for the avian toxicity endpoints observed in regulatory studies. We linked this new implementation to a toxicological module to capture observed pesticide effects on reproduction via a decreased efficiency of egg production. We analysed ten reproduction studies with five different pesticides conducted with the mallard (Anas platyrhynchos) and the northern bobwhite (Colinus virginianus). The new model implementation accurately distinguished between effects on egg production from direct mechanism of toxicity and from food avoidance. Due to the specific nature of regulatory studies, model applicability for risk refinement is currently limited. We provide suggestions for next steps in model development.

hb or not hb - When and why accounting for background mortality in toxicological survival models matters?

Decisions in Environmental Risk Assessment (ERA) about impacts of chemical compounds on different species are based on critical effect indicators such as the 50% lethal concentration (LC50). Regulatory documents recommend concentration-response (or concentration-effect) model fitting on standard toxicity test data to get LC50 values. However, toxicokinetic-toxicodynamic (TKTD) models proved their efficiency to better exploit toxicity test data, at Tier-2 but also at Tier-1, delivering time-independent indicators. In particular, LC50 values can be obtained from the reduced General Unified Threshold model of Survival (GUTS-RED) with both variants, Stochastic Death and Individual Tolerance, that include parameter hb, the background mortality. Estimating hb during the fitting process or not depends on studies and fitting habits, while it may strongly influence the other GUTS-RED parameters, and consequently the LC50 estimate. We hypothesized that estimating hb from all data in all replicates over time should provide more precise LC50 estimates. We then explored how estimating hb impacted: (i) GUTS-RED model parameters; (ii) goodness-of-fit criteria (fitting plot, posterior predictive check, parameter correlations); (iii) LC50 accuracy and precision. We finally show that estimating hb does not impact the LC50 precision while providing more accurate and precise GUTS parameter estimates. Hence, estimating hb would lead to a more protective ERA.

Behavioral toxicity, histopathological alterations and oxidative stress in Tubifex tubifex exposed to aromatic carboxylic acids- acetic acid and benzoic acid: A comparative time-dependent toxicity assessment.

This study evaluated Acetic acid (AA) and Benzoic acid's (BA) acute and sublethal toxicity by observing mortality, behavioral responses, and changes in the levels of oxidative stress enzymes in Tubifex tubifex. Exposure-induced changes in antioxidant activity (Catalase, Superoxide dismutase), oxidative stress (Malondialdehyde concentrations), and histopathological alterations in the tubificid worms were also noted across exposure intervals. The 96 h LC50 values of AA and BA to T. tubifex were 74.99 and 37.15 mg/l, respectively. Severity in behavioral alterations (including increased mucus production, wrinkling, and reduction in clumping) and autotomy showed concentration-dependent trends for both toxicants. Although histopathological effects also showed marked degeneration in the alimentary and integumentary systems in highest exposure groups (worms exposed to 14.99 mg/l for AA and 7.42 mg/l for BA) for both toxicants. Antioxidant enzymes (catalase and superoxide dismutase) also showed a marked increase of up to 8-fold and 10-fold for the highest exposure group of AA and BA respectively. While species sensitivity distribution analysis revealed T. tubifex as most sensitive to AA and BA compared to other freshwater vertebrates and invertebrates, General Unified Threshold model of Survival (GUTS) predicted individual tolerance effects (GUTS-IT), with slower potential for toxicodynamic recovery, as a more likely pathway for population mortality. Study findings demonstrate BA with greater potential for ecological effects compared to AA within 24 h of exposure. Furthermore, ecological risks to critical detritus feeders like T. tubifex may have severe implications for ecosystem services and nutrient availability within freshwater habitats.

Modelling the effects of the pyrethroid insecticide cypermethrin on the life cycle of the soil dwelling annelid Enchytraeus crypticus, an original experimental design to calibrate a DEB-TKTD model.

The Dynamic Energy Budget theory (DEB) enables ecotoxicologists to model the effects of chemical stressors on organism life cycles through the coupling of toxicokinetic-toxicodynamic (TK-TD) models. While good progress has been made in the application of DEB-TKTD models for aquatic organisms, applications for soil fauna are scarce, due to the lack of dedicated experimental designs suitable for collecting the required time series effect data. Enchytraeids (Annelida: Clitellata) are model organisms in soil ecology and ecotoxicology. They are recognised as indicators of biological activity in soil, and chemical stress in terrestrial ecosystems. Despite this, the application of DEB-TKTD models to investigate the impact of chemicals has not yet been tested on this family. Here we assessed the impact of the pyrethroid insecticide cypermethrin on the life cycle of Enchytraeus crypticus. We developed an original experimental design to collect the data required for the calibration of a DEB-TKTD model for this species. E. crypticus presented a slow initial growth phase that has been successfully simulated with the addition of a size-dependent food limitation for juveniles in the DEB model. The DEB-TKTD model simulations successfully agreed with the data for all endpoints and treatments over time. The highlighted physiological mode of action (pMoA) for cypermethrin was an increase of the growth energy cost. The threshold for effects on survival was estimated at 73.14 mg kg- 1, and the threshold for effects on energy budget (i.e., sublethal effects) at 19.21 mg kg- 1. This study demonstrates that DEB-TKTD models can be successfully applied to E. crypticus as a representative soil species, and may improve the ecological risk assessment for terrestrial ecosystems, and our mechanistic understanding of chemical effects on non-target species.

Higher Risks of Copper Toxicity in Turbid Waters: Quantifying the Bioavailability of Particle-Bound Metals to Set Site-Specific Water Quality Criteria.

In coastal waters, particulate metals constitute a substantial fraction of the total metals; however, the prevalent water quality criteria are primarily based on dissolved metals, seemingly neglecting the contribution of particulate metals. Here we developed a method to quantify the toxicity risk of particulate metals, and proposed a way to calculate modifying factors (MFs) for setting site-specific criteria in turbid waters. Specifically, we used a side-by-side experimental design to study copper (Cu) bioaccumulation and toxicity in an estuarine clam, Potamocorbula laevis, under the exposure to "dissolved only" and "dissolved + particulate" 65Cu. A toxicokinetic-toxicodynamic model (TK-TD) was used to quantify the processes of Cu uptake, ingestion, assimilation, egestion, and elimination, and to relate mortality risk to tissue Cu. We find that particulate Cu contributes 40-67% of the Cu bioaccumulation when the suspended particulate matter (SPM) ranges from 12 to 229 mg L-1. The Cu-bearing SPM also increases the sensitivity of organisms to internalized Cu by decreasing the internal threshold concentration (CIT) from 141 to 76.8 µg g-1. MFs were derived based on the TK-TD model to consider the contribution of particulate Cu (in the studied SPM range) for increasing Cu bioaccumulation (MF = 1.3-2.2) and toxicity (MF = 2.3-3.9). Water quality criteria derived from dissolved metal exposure need to be lowered by dividing by an MF to provide adequate protection. Overall, the method we developed provides a scientifically sound framework to manage the risks of metals in turbid waters.

Liquid-Phase Microextraction Approaches for Preconcentration and Analysis of Chiral Compounds: A Review on Current Advances.

Chirality is a critical issue in pharmaceutics, forensic chemistry, therapeutic drug monitoring, doping control, toxicology, or environmental investigations as enantiomers of a chiral compound can exhibit different activities, i.e., one enantiomer can have the desired effect while the other one can be inactive or even toxic. To monitor enantioselective metabolism or toxicokinetic/toxicodynamic mechanisms in extremely low content in biological or environmental matrices, sample preparation is vital. The present review describes current status of development of liquid-phase microextraction approaches such as hollow fiber liquid-phase microextraction (HF-LPME), electromembrane extraction (EME), dispersive liquid-liquid microextraction (DLLME), and supramolecular solvent-based microextraction (SSME), used for sample preparation of enantiomers/chiral compounds. The advantages and limitations of the above techniques are discussed. Attention is also focused on chiral separation approaches commonly applied to study the stereo-selective metabolism or toxicokinetic/toxicodynamic mechanisms of enantiomers in the biological and environmental samples.

Effect of dissolved organic matter on the bioavailability and toxicity of cadmium in zebrafish larvae: Determination based on toxicokinetic-toxicodynamic processes.

The presence of dissolved organic matter (DOM) strongly influences the bioavailability of metals in aquatic environments; however, the association between the binding activities and the concentrations of DOM compositions is not well documented, leading to uncertainties in metal toxicity assessment. We creatively quantify the mitigation and acceleration effects of DOM compositions on cadmium (Cd) bioaccumulation and toxicity in zebrafish larvae using abiotic ligand (ABLs) and biotic ligand (BLs) in a toxicokinetic-toxicodynamic (TK-TD) model. The BL-TK-TD model could accurately predict the protective effect of fulvic acid while overestimating the complexing capacity of citric acid. The model also could successfully simulate the protective effects of native DOM in most cases from 32 natural water bodies in China. The observed LC50 values of Cd showed a peak effect for the native DOM fraction comprising hydrophilic acidic contents (3.55 ± 0.44 mg L - 1) in natural water from 32 sites. The BL-TK-TD model provides practically useful information to identify the effect of different DOM compositions on metal bioavailability and toxicity in aquatic environments and guides future water management policies aimed at controlling aquatic heavy metal pollution.

Elucidating the Differences in Metal Toxicity by Quantitative Adverse Outcome Pathways.

Numerous studies have reported that the toxicity differences among metals are widespread; however, little is known about the mechanism of differences in metal toxicity to aquatic organisms due to the lack of quantitative understanding of their adverse outcome pathway. Here, we investigated the effects of Cd and Cu on bioaccumulation, gene expression, physiological responses, and apical effects in zebrafish larvae. RNA sequencing was conducted to provide supplementary mechanistic information for the effects of Cd and Cu exposure. On this basis, we proposed a quantitative adverse outcome pathway (qAOP) suitable for metal risk assessment of aquatic organisms. Our work provides a mechanistic explanation for the differences in metal toxicity where the strong bioaccumulation of Cu enables the newly accumulated Cu to reach the threshold that causes different adverse effects faster than Cd in zebrafish larvae, resulting in a higher toxicity of Cu than that of Cd. Furthermore, we proposed a parameter CIT/BCF (the ratio of internal threshold concentration and bioaccumulation factor) that helps to understand the toxicity differences by combining the information of bioaccumulation and internal threshold of adverse effects. This work demonstrated that qAOP is an effective quantitative tool for understanding the toxicity mechanism and highlight the importance of toxicokinetics and toxicodynamics at different biological levels in determining the metal toxicity.

Development of physiologically-based toxicokinetic-toxicodynamic (PBTK-TD) model for 4-nonylphenol (4-NP) reflecting physiological changes according to age in males: Application as a new risk assessment tool with a focus on toxicodynamics.

Environmental exposure to 4-nonylphenol (4-NP) is extensive, and studies related to human risk assessment must continue. Especially, prediction of toxicodynamics (TDs) related to reproductive toxicity in males is very important in risk-level assessment and management of 4-NP. This study aimed to develop a physiologically-based-toxicokinetic-toxicodynamic (PBTK-TD) model that added a TD prostate model to the previously reported 4-n-nonylphenol (4-n-NP) physiologically-based-pharmacokinetic (PBPK) model. Modeling was performed under the assumption of similar TKs between 4-n-NP and 4-NP because TK experiments on 4-NP, a random-mixture, are practically difficult. This study was very important to quantitatively predict the TKs and TDs of 4-NP by age at exposure using an advanced PBTK-TD model that reflected physiological-changes according to age. TD-modeling was performed based on the reported toxic effects of 4-NP on RWPE-1 cells, a human-prostate-epithelial-cell-line. Through a meta-analysis of reported human physiological data, body weight, tissue volume, and blood flow rate patterns according to age were mathematically modeled. These relationships were reflected in the PBTK-TD model for 4-NP so that the 4-NP TK and TD changes according to age and their differences could be confirmed. Differences in TK and TD parameters of 4-NP at various ages were not large, within 3.61-fold. Point-of-departure (POD) and reference-doses for each age estimated using the model varied as 426.37-795.24 and 42.64-79.52 µg/kg/day, but the differences (in POD or reference doses between ages) were not large, at less than 1.87-times. The PBTK-TD model simulation predicted that even a 100-fold 4-NP PODman dose would not have large toxicity to the prostate. With a focus on TDs, the predicted maximum possible exposure of 4-NP was as high as 6.06-23.60 mg/kg/day. Several toxicity-related values estimated by the dose-response curve were higher than those calculated, depending upon the PK or TK, which would be useful as a new exposure limit for prostate toxicity of 4-NP.