The Scanning Cryomacroscope - A Device Prototype for the Study of Cryopreservation

Thermomechanical stress analyses of nanowarming-assisted recovery from cryopreservation by vitrification in human heart and rat heart models

This study investigates thermomechanical stress in cryopreservation by vitrification of the heart, while exploring the effects of nanowarming-assisted recovery from cryogenic storage. This study expands upon a recently published study, combining experimental investigation and thermal analysis of cryopreservation on a rat heart model. Specifically, this study focuses on scenarios with variable concentrations of silica-coated iron-oxide nanoparticles (sIONPs), while accounting for loading limitations associated with the heart physiology, as well as the properties of cryoprotective agent (CPA) solution and the geometry of the container. Results of this study suggest that variable sIONP concentration based on the heart physiology will elevate mechanical stresses when compared with the mathematically simplified, uniform distribution case. The most dangerous part of rewarming is below glass transition and at the onset of nanowarming past the glass transition temperature on the way for organ recovery from cryogenic storage. Throughout rewarming, regions that rewarm faster, such as the chambers of the heart (higher sIONP concentration), undergo compressive stresses, while the slower rewarming regions, such as the heart myocardium (low sIONP concentration), undergo tension. Being a brittle material, the vitrified organ is expected to fail under tension in lower stresses than in compression. Unfortunately, the location and magnitude of the maximum stress in the investigated cases varied, while general rules were not identified. This investigation demonstrates the need to tailor the thermal protocol of heart cryopreservation on a case-by-case basis, since the location, orientation, magnitude, and instant at which the maximum mechanical stress is found cannot be predicted a priori. While thermomechanical stress poses a significant risk to organ integrity, careful design of the thermal protocol can be instrumental in reducing the likelihood of structural damage, while taking full advantage of the benefits of nanowarming.

Large surface deformation due to thermo-mechanical effects during cryopreservation by vitrification - mathematical model and experimental validation

This study presents a simplified thermal-fluids (TF) mathematical model to analyze large surface deformations in cryoprotective agents (CPA) during cryopreservation by vitrification. The CPA deforms during vitrification due to material flow caused by the combined effects of thermal gradients within the domain, thermal contraction due to temperature, and exponential increase in the viscosity of the CPA as it is cooled towards glass transition. While it is well understood that vitrification is associated with thermo-mechanical stress, which might lead to structural damage, those large deformations can lead to stress concentration, further intensifying the probability to structural failure. The results of the TF model are experimentally validated by means of cryomacroscopy on a cuvette containing 7.05M dimethyl sulfoxide (DMSO) as a representative CPA. The TF model presented in this study is a simplified version of a previously presented thermo-mechanics (TM) model, where the TM model is set to solve the coupled heat transfer, fluid mechanics and solid mechanics problems, while the TF model omits further deformations in the solid state. It is demonstrated in this study that the TF model alone is sufficient to capture large-body deformations during vitrification. However, the TF model alone cannot be used to estimate mechanical stresses, which become significant only when the deformation rates become so small that the deformed body practically behaves as an amorphous solid. This study demonstrates the high sensitivity of deformation predictions to variation in material properties, chief among which are the variations of density and viscosity with temperature. Finally, this study includes a discussion on the possibility of turning on and off the TF and TM models in respective parts of the domain, in order to solve the multiphysics problem in a computationally cost-effective manner.

Mathematical modeling of surface deformation during vitrification

Mathematical modeling of surface deformation during cryopreservation by vitrification is presented in this study. The specific problem under consideration is of a cryoprotective agent (CPA) solution vitrifying in a vial, following previously obtained cryomacroscopy observations. A multiphysics solution is proposed in this study, combining coupled effects associated with heat transfer, fluid mechanics, and solid mechanics. Consistent with previous investigations, this study demonstrates that surface deformation is the result of material flow, which is the combined outcome of temperature gradients developed during the inward cooling process, the tendency of the material to change its volume with temperature, and the exponential increase in material viscosity with the decreasing temperature. During this process, the behavior of the CPA changes from liquid to a solid-like amorphous material, where the arrested flow in the vitrified state results in mechanical stresses. Results of this study show a good qualitative agreement of surface deformation with previously obtained experimental data, and support prior investigations to explain fracture tendencies propagating from the deformed surface. Results of this study also highlight the effect of heat convection in the CPA at the early stage of cooling.

Thermal conductivity of cryoprotective agents loaded with nanoparticles, with application to recovery of preserved tissues and organs from cryogenic storage

The objective of this study is to provide thermal conductivity data for CPA-based nanofluids for the benefit of the analyses of cryopreservation by vitrification. Thermal conductivity measurements were conducted using a hot-wire technique on an experimentation platform of the cryomacroscope, to correlate measurements with observed physical effects such as crystallization and fracturing. Tested materials in this study include the CPA cocktails M22, VS55, DP6, and DP6+sucrose. Nanofluids in this study include the above CPA cocktails as base solutions, when mixed with either iron-oxide nanoparticles (IONP) or silica-coated iron-oxide nanoparticles (sIONP). Results of this study demonstrated the addition of sIONP to any of the CPA cocktails tested did not significantly affect its thermal conductivity, its tendency to vitrify or, conversely, its tendency to form rewarming phase crystallization (RPC). Fractures were observed with cryomacroscopy at the onset of rewarming for DP6+sIONP under carefully controlled rewarming conditions without RF activation, despite the inherent opacity of the sIONP solutions. It is likely that using RF heating in order to accelerate rewarming while unifying the temperature distribution would prevent fracture and RPC. However, sIONP were not activated in this study, as the RF heating mechanism would interfere with thermal conductivity measurements. The addition of IONP to DP6 appears to hinder the tendency of the CPA to vitrify, which is a detrimental effect. But it is unlikely that uncoated nanoparticle solutions will be used in practical applications.

Analysis of polarized-light effects in glass-promoting solutions with applications to cryopreservation and organ banking

This study presents experimental results and an analysis approach for polarized light effects associated with thermomechanical stress during cooling of glass promoting solutions, with applications to cryopreservation and tissue banking in a process known as vitrification. Polarized light means have been previously integrated into the cryomacroscope-a visualization device to detect physical effects associated with cryopreservation success, such as crystallization, fracture formation, and contamination. The experimental study concerns vitrification in a cuvette, which is a rectangular container. Polarized light modeling in the cuvette is based on subdividing the tridimensional (3D) domain into a series of planar (2D) problems, for which a mathematical solution is available in the literature. The current analysis is based on tracking the accumulated changes in light polarization and magnitude, as it passes through the sequence of planar problems. Results of this study show qualitative agreement in light intensity history and distribution between experimental data and simulated results. The simulated results help explaining differences between 2D and 3D effects in photoelasticity, most notably, the counterintuitive observation that high stress areas may correlate with low light intensity regions based on the particular experimental conditions. Finally, it is suggested that polarized-light analysis must always be accompanied by thermomechanical stress modeling in order to explain 3D effects.

Polarized light scanning cryomacroscopy, part II: Thermal modeling and analysis of experimental observations

This study aims at developing thermal analysis tools and explaining experimental observations made by means of polarized-light cryomacroscopy (Part I). Thermal modeling is based on finite elements analysis (FEA), where two model parameters are extracted from thermal measurements: (i) the overall heat transfer coefficient between the cuvette and the cooling chamber, and (ii) the effective thermal conductivity within the cryoprotective agent (CPA) at the upper part of the cryogenic temperature range. The effective thermal conductivity takes into account enhanced heat transfer due to convection currents within the CPA, creating the so-called Bénard cells. Comparison of experimental results with simulation data indicates that the uncertainty in simulations due to the propagation of uncertainty in measured physical properties exceeds the uncertainty in experimental measurements, which validates the modeling approach. It is shown in this study that while a cavity may form in the upper-center portion of the vitrified CPA, it has very little effect on estimating the temperature distribution within the domain. This cavity is driven by thermal contraction of the CPA, with the upper-center of the domain transitioning to glass last. Finally, it is demonstrated in this study that additional stresses may develop within the glass transition temperature range due to nonlinear behavior of the thermal expansion coefficient. This effect is reported here for the first time in the context of cryobiology, using the capabilities of polarized-light cryomacroscopy.

Polarized light scanning cryomacroscopy, part I: Experimental apparatus and observations of vitrification, crystallization, and photoelasticity effects

Cryomacroscopy is an effective means to observe physical events affecting cryopreservation success in large-size specimens. The current study aims at integrating polarized-light in the study of large-size cryopreservation, using the scanning cryomacroscope as a development platform. Results of this study demonstrate polarized light as a visualization enhancement means, including the following effects: contaminants in the CPA solution, crystallization, fracture formation, thermal contraction, and solute precipitation. In addition, photoelasticity effects are used to demonstrate the development of residual stresses and the potential for stress relaxation above the glass transition temperature. Furthermore, this study suggests that the ability to periodically switch between non-polarized light and polarized light is an essential feature of investigation. When using polarized light for example, a dark region may represent a free-of-stress and free-of-crystals material, or fully crystallized material, which may potentially experience mechanical stress; switching to a non-polarized light would help to distinguish between the different cases. The analysis of thermo-mechanical stress in cryopreservation is essentially based on four key elements: identification of physical events, knowledge of physical properties, thermal analysis of the specimen, and description of the mechanical behavior of the cryopreserved material (also known as the constitutive law). With the above knowledge, one can investigate the conditions to preserve structural integrity. While the current study aims at identification of physical events, critical knowledge on physical properties and mechanical behavior has already been developed in previous studies. The companion manuscript (Part II) aims at providing means for thermal analysis in the specimen, which will serve as the basis for a multi-scale analysis of thermo-mechanical stress in large-size specimens.

Thermal conductivity of the cryoprotective cocktail DP6 in cryogenic temperatures, in the presence and absence of synthetic ice modulators

The thermal conductivity of the cryoprotective agent (CPA) cocktail DP6 in combination with synthetic ice modulators (SIMs) is measured in this study, using a transient hot-wire method. DP6 is a mixture of 3 M dimethyl sulfoxide (DMSO) and 3 M propylene glycol, which received significant attention in the cryobiology community in recent years. Tested SIMs include 6% 1,3Cyclohexanediol, 6% 2,3Butanediol, and 12% PEG400 (percentage by volume). This study integrates the scanning cryomacroscope for visual verification of crystallization and vitrification events. It is demonstrated that the thermal conductivity of the vitrifying CPA cocktail decreases monotonically with the decreasing temperature down to -180 °C. By contrast, the thermal conductivity of the crystalline material increases with decreasing temperature in the same temperature range. Results of this study demonstrate that the thermal conductivity may vary by three fold between the amorphous and crystalline phases of DP6 below the glass transition temperature of DP6 (Tg = -119 °C). The selected SIMs demonstrate the ability to inhibit crystallization in DP6, even at subcritical cooling rates. An additional ice suppression capability is observed by the Euro-Collins as a vehicle solution, disproportionate to its volume ratio in the cocktail. The implication of the observed thermal conductivity differences between the amorphous and crystalline phases of the same cocktail on cryopreservation simulations is significant in some cases and must be taken into account in thermal analyses of cryopreservation protocols.

Large Thermal Conductivity Differences between the Crystalline and Vitrified States of DMSO with Applications to Cryopreservation

Thermal conductivity of dimethyl-sulfoxide (DMSO) solution is measured in this study using a transient hot wire technique, where DMSO is a key ingredient in many cryoprotective agent (CPA) cocktails. Characterization of thermal properties of cryoprotective agents is essential to the analysis of cryopreservation processes, either when evaluating experimental data or for the design of new protocols. Also presented are reference measurements of thermal conductivity for pure water ice and glycerol. The thermal conductivity measurement setup is integrated into the experimentation stage of a scanning cryomacroscope apparatus, which facilitates the correlation of measured data with visualization of physical events. Thermal conductivity measurements were conducted for a DMSO concentration range of 2M and 10M, in a temperature range of -180°C and 25°C. Vitrified samples showed decreased thermal conductivity with decreasing temperature, while crystalline samples showed increased thermal conductivity with decreasing temperature. These different behaviors result in up to a tenfold difference in thermal conductivity at -180°C. Such dramatic differences can drastically impact heat transfer during cryopreservation and their quantification is therefore critical to cryobiology.

On the Effects of Thermal History on the Development and Relaxation of Thermo-Mechanical Stress in Cryopreservation

This study investigates the effects of the thermal protocol on the development and relaxation of thermo-mechanical stress in cryopreservation by means of glass formation, also known as vitrification. The cryopreserved medium is modeled as a homogeneous viscoelastic domain, constrained within either a stiff cylindrical container or a highly compliant bag. Annealing effects during the cooling phase of the cryopreservation protocol are analyzed. Results demonstrate that an intermediate temperature-hold period can significantly reduce the maximum tensile stress, thereby decreasing the potential for structural damage. It is also demonstrated that annealing at temperatures close to glass transition significantly weakens the dependency of thermo-mechanical stress on the cooling rate. Furthermore, a slower initial rewarming rate after cryogenic storage may drastically reduce the maximum tensile stress in the material, which supports previous experimental observations on the likelihood of fracture at this stage. This study discusses the dependency of the various stress components on the storage temperature. Finally, it is demonstrated that the stiffness of the container wall can affect the location of maximum stress, with implications on the development of cryopreservation protocols.