Potential of the zebrafish (Danio rerio) embryo test to discriminate between chemicals of similar molecular structure-a study with valproic acid and 14 of its analogues.

Valproic acid is a frequently used antiepileptic drug and known pediatric hepatotoxic agent. In search of pharmaceuticals with increased effectiveness and reduced toxicity, analogue chemicals came into focus. So far, toxicity and teratogenicity data of drugs and metabolites have usually been collected from mammalian model systems such as mice and rats. However, in an attempt to reduce mammalian testing while maintaining the reliability of toxicity testing of new industrial chemicals and drugs, alternative test methods are being developed. To this end, the potential of the zebrafish (Danio rerio) embryo to discriminate between valproic acid and 14 analogues was investigated by exposing zebrafish embryos for 120 h post fertilization in the extended version of the fish embryo acute toxicity test (FET; OECD TG 236), and analyzing liver histology to evaluate the correlation of liver effects and the molecular structure of each compound. Although histological evaluation of zebrafish liver did not identify steatosis as the prominent adverse effect typical in human and mice, the structure-activity relationship (SAR) derived was comparable not only to human HepG2 cells, but also to available in vivo mouse and rat data. Thus, there is evidence that zebrafish embryos might serve as a tool to bridge the gap between subcellular, cell-based systems and vertebrate models.

Multistate models of developmental toxicity: Application to valproic acid-induced malformations in the zebrafish embryo.

For the determination of acute toxicity of chemicals in zebrafish (Danio rerio) embryos, the OECD test guideline 236, relative to the Fish Embryo Toxicity Test (FET), stipulates a dose-response analysis of four lethal core endpoints and a quantitative characterization of abnormalities including their time-dependency. Routinely, the data are analyzed at the different observation times separately. However, observations at a given time strongly depend on the previous effects and should be analyzed jointly with them. To solve this problem, we developed multistate models for occurrence of developmental malformations and live events in zebrafish embryos exposed to eight concentrations of valproic acid (VPA) the first five days of life. Observations were recorded daily per embryo. We statistically infer on model structure and parameters using a numerical Bayesian framework. Hatching probability rate changed with time and we compared five forms of its time-dependence; a constant rate, a piecewise constant rate with a fixed hatching time at 48 h post fertilization, a piecewise constant rate with a variable hatching time, as well as a Hill and Gaussian form. A piecewise constant function of time adequately described the hatching data. The other transition rates were conditioned on the embryo body concentration of VPA, obtained using a physiologically-based pharmacokinetic model. VPA impacted mostly the malformation probability rate in hatched and non-hatched embryos. Malformation reversion probability rates were lowered by VPA. Direct mortality was low at the concentrations tested, but increased linearly with internal concentration. The model makes full use of data and gives a finer grain analysis of the teratogenic effects of VPA in zebrafish than the OECD-prescribed approach. We discuss the use of the model for obtaining toxicological reference values suitable for inter-species extrapolation. A general result is that complex multistate models can be efficiently evaluated numerically.

Neurotoxic effects in zebrafish embryos by valproic acid and nine of its analogues: the fish-mouse connection?

Since teratogenicity testing in mammals is a particular challenge from an animal welfare perspective, there is a great need for the development of alternative test systems. In this context, the zebrafish (Danio rerio) embryo has received increasing attention as a non-protected embryonic vertebrate in vivo model. The predictive power of zebrafish embryos for general vertebrate teratogenicity strongly depends on the correlation between fish and mammals with respect to both overall general toxicity and more specific endpoints indicative of certain modes-of-action. The present study was designed to analyze the correlation between (1) effects of valproic acid and nine of its analogues in zebrafish embryos and (2) their known neurodevelopmental effects in mice. To this end, zebrafish embryos exposed for 120 h in an extended version of the acute fish embryo toxicity test (FET; OECD TG 236) were analyzed with respect to an extended list of sublethal endpoints. Particular care was given to endpoints putatively related to neurodevelopmental toxicity, namely jitter/tremor, deformation of sensory organs (eyes) and craniofacial deformation, which might correlate to neural tube defects caused by valproic acid in mammals. A standard evaluation of lethal (LC according to OECD TG 236) and sublethal toxicity (EC) merely indicated that four out of ten compounds tested in zebrafish correlate with positive results in mouse in vivo studies. A detailed assessment of more specific effects, however, namely, jitter/tremor, small eyes and craniofacial deformation, resulted in a correspondence of 75% with in vivo mouse data. A refinement of endpoint analysis from an integration of all observations into one LCx or ECx data (as foreseen by current ecotoxicology-driven OECD guidelines) to a differential evaluation of endpoints specific of selected modes-of-action thus increases significantly the predictive power of the zebrafish embryo model for mammalian teratogenicity. However, for some of the endpoints observed, e.g., scoliosis, lordosis, pectoral fin deformation and lack of movement, further experiments are required for the identification of underlying modes-of-action and an unambiguous interpretation of their predictive power for mammalian toxicity.

Integrate mechanistic evidence from new approach methodologies (NAMs) into a read-across assessment to characterise trends in shared mode of action.

Read-across approaches often remain inconclusive as they do not provide sufficient evidence on a common mode of action across the category members. This read-across case study on thirteen, structurally similar, branched aliphatic carboxylic acids investigates the concept of using human-based new approach methods, such as in vitro and in silico models, to demonstrate biological similarity. Five out of the thirteen analogues have preclinical in vivo studies. Three out of them induced lipid accumulation or hypertrophy in preclinical studies with repeated exposure, which leads to the read-across hypothesis that the analogues can potentially induce hepatic steatosis. To confirm the selection of analogues, the expression patterns of the induced differentially expressed genes (DEGs) were analysed in a human liver model. With increasing dose, the expression pattern within the tested analogues got more similar, which serves as a first indication of a common mode of action and suggests differences in the potency of the analogues. Hepatic steatosis is a well-known adverse outcome, for which over 55 adverse outcome pathways have been identified. The resulting adverse outcome pathway (AOP) network, comprised a total 43 MIEs/KEs and enabled the design of an in vitro testing battery. From the AOP network, ten MIEs, early and late KEs were tested to systematically investigate a common mode of action among the grouped compounds. The targeted testing of AOP specific MIE/KEs shows that biological activity in the category decreases with side chain length. A similar trend was evident in measuring liver alterations in zebra fish embryos. However, activation of single MIEs or early KEs at in vivo relevant doses did not necessarily progress to the late KE "lipid accumulation". KEs not related to the read-across hypothesis, testing for example general mitochondrial stress responses in liver cells, showed no trend or biological similarity. Testing scope is a key issue in the design of in vitro test batteries. The Dempster-Shafer decision theory predicted those analogues with in vivo reference data correctly using one human liver model or the CALUX reporter assays. The case study shows that the read-across hypothesis is the key element to designing the testing strategy. In the case of a good mechanistic understanding, an AOP facilitates the selection of reliable human in vitro models to demonstrate a common mode of action. Testing DEGs, MIEs and early KEs served to show biological similarity, whereas the late KEs become important for confirmation, as progression from MIEs to AO is not always guaranteed.

Development of a generic zebrafish embryo PBPK model and application to the developmental toxicity assessment of valproic acid analogs.

In order to better explain, predict, or extrapolate to humans the developmental toxicity effects of chemicals to zebrafish (Danio rerio) embryos, we developed a physiologically-based pharmacokinetic (PBPK) model designed to predict organ concentrations of neutral or ionizable chemicals, up to 120 h post-fertilization. Chemicals' distribution is modeled in the cells, lysosomes, and mitochondria of ten organs of the embryo. The model's partition coefficients are calculated with sub-models using physicochemical properties of the chemicals of interest. The model accounts for organ growth and changes in metabolic clearance with time. We compared ab initio model predictions to data obtained on culture medium and embryo concentrations of valproic acid (VPA) and nine analogs during continuous dosing under the OECD test guideline 236. We further improved the predictions by estimating metabolic clearance and partition coefficients from the data by Bayesian calibration. We also assessed the performance of the model at reproducing data published by Brox et al. (2016) on VPA and 16 other chemicals. We finally compared dose-response relationships calculated for mortality and malformations on the basis of predicted whole embryo concentrations versus those based on nominal water concentrations. The use of target organ concentrations substantially shifted the magnitude of dose-response parameters and the relative toxicity ranking of chemicals studied.