Toxicokinetic Modeling of Per- and Polyfluoroalkyl Substance Concentrations within Developing Zebrafish (Danio rerio) Populations.

Per- and polyfluoroalkyl substances (PFAS) are pervasive environmental contaminants, and their relative stability and high bioaccumulation potential create a challenging risk assessment problem. Zebrafish (Danio rerio) data, in principle, can be synthesized within a quantitative adverse outcome pathway (qAOP) framework to link molecular activity with individual or population level hazards. However, even as qAOP models are still in their infancy, there is a need to link internal dose and toxicity endpoints in a more rigorous way to further not only qAOP models but adverse outcome pathway frameworks in general. We address this problem by suggesting refinements to the current state of toxicokinetic modeling for the early development zebrafish exposed to PFAS up to 120 h post-fertilization. Our approach describes two key physiological transformation phenomena of the developing zebrafish: dynamic volume of an individual and dynamic hatching of a population. We then explore two different modeling strategies to describe the mass transfer, with one strategy relying on classical kinetic rates and the other incorporating mechanisms of membrane transport and adsorption/binding potential. Moving forward, we discuss the challenges of extending this model in both timeframe and chemical class, in conjunction with providing a conceptual framework for its integration with ongoing qAOP modeling efforts.

Multiple cryoprotectant toxicity model for vitrification solution optimization.

Vitrification is a promising cryopreservation technique for complex specimens such as tissues and organs. However, it is challenging to identify mixtures of cryoprotectants (CPAs) that prevent ice formation without exerting excessive toxicity. In this work, we developed a multi-CPA toxicity model that predicts the toxicity kinetics of mixtures containing five of the most common CPAs used in the field (glycerol, dimethyl sulfoxide (DMSO), propylene glycol, ethylene glycol, and formamide). The model accounts for specific toxicity, non-specific toxicity, and interactions between CPAs. The proposed model shows reasonable agreement with training data for single and binary CPA solutions, as well as ternary CPA solution validation data. Sloppy model analysis was used to examine the model parameters that were most important for predictions, providing clues about mechanisms of toxicity. This analysis revealed that the model terms for non-specific toxicity were particularly important, especially the non-specific toxicity of propylene glycol, as well as model terms for specific toxicity of formamide and interactions between formamide and glycerol. To demonstrate the potential for model-based design of vitrification methods, we paired the multi-CPA toxicity model with a published vitrification/devitrification model to identify vitrifiable CPA mixtures that are predicted to have minimal toxicity. The resulting optimized vitrification solution composition was a mixture of 7.4 molal glycerol, 1.4 molal DMSO, and 2.4 molal formamide. This demonstrates the potential for mathematical optimization of vitrification solution composition and sets the stage for future studies to optimize the complete vitrification process, including CPA mixture composition and CPA addition and removal methods.

General tissue mass transfer model for cryopreservation applications.

Successful cryopreservation of complex specimens, such as tissues and organs, would greatly benefit both the medical and scientific research fields. Vitrification is one of the most promising techniques for complex specimen cryopreservation, but toxicity remains a major challenge because of the high concentration of cryoprotectants (CPAs) needed to vitrify. Our group has approached this problem using mathematical optimization to design less toxic CPA equilibration methods for cells. To extend this approach to tissues, an appropriate mass transfer model is required. Fick's law is commonly used, but this simple modeling framework does not account for the complexity of mass transfer in tissues, such as the effects of fixed charges, tissue size changes, and the interplay between cell membrane transport and transport through the extracellular fluid. Here, we propose a general model for mass transfer in tissues that accounts for all of these phenomena. To create this model, we augmented a previously published acellular model of mass transfer in articular cartilage to account for the effects of cells. We show that the model can accurately predict changes in CPA concentration and tissue size for both articular cartilage and pancreatic islets, tissue types with vastly different properties.

Rapid quantification of multi-cryoprotectant toxicity using an automated liquid handling method.

Cryopreservation in a vitrified state has vast potential for long-term storage of tissues and organs that may be damaged by ice formation. However, the toxicity imparted by the high concentration of cryoprotectants (CPAs) required to vitrify these specimens remains a hurdle. To address this challenge, we previously developed a mathematical approach to design less toxic CPA equilibration methods based on the minimization of a toxicity cost function. This approach was used to design improved methods for equilibration of bovine pulmonary artery endothelial cells (BPAEC) with glycerol. To fully capitalize on the toxicity cost function approach, it is critical to describe the toxicity kinetics of additional CPAs, including multi-CPA mixtures that are commonly used for vitrification. In this work, we used automated liquid handling to characterize the toxicity kinetics of five of the most common CPAs (glycerol, dimethyl sulfoxide (DMSO), propylene glycol, ethylene glycol, and formamide), along with their binary and ternary mixtures for BPAEC. In doing so, we developed experimental methods that can be used to determine toxicity kinetics more quickly and accurately. Our results highlight some common CPA toxicity trends, including the relatively low toxicity of ethylene glycol and a general increase in toxicity as the CPA concentration increases. Our results also suggest potential new approaches to reduce toxicity, including a surprising toxicity neutralization effect of glycerol on formamide. In the future, this dataset will serve as the basis to expand our CPA toxicity model, enabling application of the toxicity cost function approach to vitrification solutions containing multiple CPAs.

Osmotic response during kidney perfusion with cryoprotectant in isotonic or hypotonic vehicle solution.

Organ cryopreservation would revolutionize transplantation by overcoming the shelf-life limitations of conventional organ storage. To prepare an organ for cryopreservation, it is first perfused with cryoprotectants (CPAs). These chemicals can enable vitrification during cooling, preventing ice damage. However, CPAs can also cause toxicity and osmotic damage. It is a major challenge to find the optimal balance between protecting the cells from ice and avoiding CPA-induced damage. In this study, we examined the organ perfusion process to shed light on phenomena relevant to cryopreservation protocol design, including changes in organ size and vascular resistance. In particular, we compared perfusion of kidneys (porcine and human) with CPA in either hypotonic or isotonic vehicle solution. Our results demonstrate that CPA perfusion causes kidney mass changes consistent with the shrink-swell response observed in cells. This response was observed when the kidneys were relatively fresh, but disappeared after prolonged warm and/or cold ischemia. Perfusion with CPA in a hypotonic vehicle solution led to a significant increase in vascular resistance, suggesting reduced capillary diameter due to cell swelling. This could be reversed by switching to perfusion with CPA in isotonic vehicle solution. Hypotonic vehicle solution did not cause notable osmotic damage, as evidenced by low levels of lactate dehydrogenase (LDH) in the effluent, and it did not have a statistically significant effect on the delivery of CPA into the kidney, as assessed by computed tomography (CT). Overall, our results show that CPA vehicle solution tonicity affects organ size and vascular resistance, which may have important implications for cryopreservation protocol design.

Mathematical Modeling of Protectant Transport in Tissues.

Mass transfer of protectant chemicals is a fundamental aspect of cryopreservation and freeze-drying protocols. As such, mass transfer modeling is useful for design of preservation methods. Cell membrane transport modeling has been successfully used to guide design of preservation methods for isolated cells. For tissues, though, there are several mass transfer modeling challenges that arise from phenomena associated with cells being embedded in a tissue matrix. Both cells and the tissue matrix form a barrier to the free diffusion of water and protective chemicals. Notably, the extracellular space becomes important to model. The response of cells embedded in the tissue is dependent on the state of the extracellular space which varies both spatially and temporally. Transport in the extracellular space can also lead to changes in tissue size. In this chapter, we describe various mass transfer models that can be used to describe transport phenomena occurring during loading of tissues with protective molecules for cryopreservation applications. Assumptions and simplifications that limit the applicability of each of these models are discussed.