Mathematical modeling of cryoprotectant addition and removal for the cryopreservation of engineered or natural tissues

Rational synthesis of total damage during cryoprotectant equilibration: modelling and experimental validation of osmomechanical, temperature, and cytotoxic damage in sea urchin (Paracentrotus lividus) oocytes

Sea urchins (e.g., Paracentrotus lividus) are important for both aquaculture and as model species. Despite their importance, biobanking of urchin oocytes by cryopreservation is currently not possible. Optimized cryoprotectant loading may enable novel vitrification methods and thus successful cryopreservation of oocytes. One method for determining an optimized loading protocol uses membrane characteristics and models of damage, namely osmomechanical damage, temperature damage (e.g., chill injury) and cytotoxicity. Here we present and experimentally evaluate existing and novel models of these damage modalities as a function of time and temperature. In osmomechanical damage experiments, oocytes were exposed for 2 to 30 minutes in hypertonic NaCl or sucrose supplemented seawater or in hypotonic diluted seawater. In temperature damage experiments, oocytes were exposed to 1.7 °C, 10 °C, or 20 °C for 2 to 90 minutes. Cytotoxicity was investigated by exposing oocytes to solutions of Me2SO for 2 to 30 minutes. We identified a time-dependent osmotic damage model, a temperature-dependent damage model, and a temperature and time-dependent cytotoxicity model. We combined these models to estimate total damage during a cryoprotectant loading protocol and determined the optimal loading protocol for any given goal intracellular cryoprotectant concentration. Given our fitted models, we find sea urchin oocytes can only be loaded to 13% Me2SO v/v with about 50% survival. This synthesis of multiple damage modalities is the first of its kind and enables a novel approach to modelling cryoprotectant equilibration survival for cells in general.

Crystallisation Degree Analysis during Cryopreservation of Biological Tissue Applying Interval Arithmetic

This paper presents the numerical modelling of heat transfer and changes proceeding in the homogeneous sample, caused by the crystallisation phenomenon during cryopreservation by vitrification. Heat transfer was simulated in a microfluidic system in which the working fluid flowed in micro-channels. The analysed process included single-phase flow during warming, and two-phase flow during cooling. In the model under consideration, interval parameters were assumed. The base of the mathematical model is given by the Fourier equation, with a heat source including the degree of ice crystallisation. The formulated problem has been solved using the interval version of the finite difference method, with the rules of the directed interval arithmetic. The fourth order Runge-Kutta algorithm has been applied to determine the degree of crystallisation. In the final part of this paper, examples of numerical computations are presented.

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.

Cryopreservation of NK and T Cells Without DMSO for Adoptive Cell-Based Immunotherapy

Dimethylsufoxide (DMSO) being universally used as a cryoprotectant in clinical adoptive cell-therapy settings to treat hematological malignancies and solid tumors is a growing concern, largely due to its broad toxicities. Its use has been associated with significant clinical side effects-cardiovascular, neurological, gastrointestinal, and allergic-in patients receiving infusions of cell-therapy products. DMSO has also been associated with altered expression of natural killer (NK) and T-cell markers and their in vivo function, not to mention difficulties in scaling up DMSO-based cryoprotectants, which introduce manufacturing challenges for autologous and allogeneic cellular therapies, including chimeric antigen receptor (CAR)-T and CAR-NK cell therapies. Interest in developing alternatives to DMSO has resulted in the evaluation of a variety of sugars, proteins, polymers, amino acids, and other small molecules and osmolytes as well as modalities to efficiently enable cellular uptake of these cryoprotectants. However, the DMSO-free cryopreservation of NK and T cells remains difficult. They represent heterogeneous cell populations that are sensitive to freezing and thawing. As a result, clinical use of cryopreserved cell-therapy products has not moved past the use of DMSO. Here, we present the state of the art in the development and use of cryopreservation options that do not contain DMSO toward clinical solutions to enable the global deployment of safer adoptively transferred cell-based therapies.

Numerical Study of Heat and Mass Transfer during Cryopreservation Process with Application of Directed Interval Arithmetic

In the present paper, numerical modelling of heat and mass transfer proceeding in a two-dimensional axially symmetrical articular cartilage sample subjected to a cryopreservation process is presented. In the model under consideration, interval parameters were assumed. The heat transfer process is described using the Fourier interval equation, while the cryoprotectant transport (DMSO) across the cell membrane is analyzed using a two-parameter model taking into account the simulation of the water volume in the chondrocytes and the change in DMSO concentration over time. The liquidus tracking (LT) protocol introduced by Pegg et al. was used to model the cryopreservation process. This procedure divides the heating and cooling phases into eight and seven steps, respectively, allowing precise regulation of temperature and cryoprotectant (CPA) concentration of bathing solutions. This protocol protects chondrocytes from ice crystal, osmotic stress, and electrolyte damage. The obtained interval concentrations of cryoprotectant in chondrocytes were compared with previous simulations obtained using the deterministic model and they are mostly in agreement with the simulation data.

Vitrification of particulated articular cartilage via calculated protocols

Preserving viable articular cartilage is a promising approach to address the shortage of graft tissue and enable the clinical repair of articular cartilage defects in articulating joints, such as the knee, ankle, and hip. In this study, we developed two 2-step, dual-temperature, multicryoprotectant loading protocols to cryopreserve particulated articular cartilage (cubes ~1 mm3 in size) using a mathematical approach, and we experimentally measured chondrocyte viability, metabolic activity, cell migration, and matrix productivity after implementing the designed loading protocols, vitrification, and warming. We demonstrated that porcine and human articular cartilage cubes can be successfully vitrified and rewarmed, maintaining high cell viability and excellent cellular function. The vitrified particulated articular cartilage was stored for a period of 6 months with no significant deterioration in chondrocyte viability and functionality. Our approach enables high-quality long-term storage of viable articular cartilage that can alleviate the shortage of grafts for use in clinically repairing articular cartilage defects.

Numerical Modeling of Heat and Mass Transfer during Cryopreservation Using Interval Analysis

In the paper, the numerical analysis of heat and mass transfer proceeding in an axially symmetrical articular cartilage sample subjected to the cryopreservation process is presented. In particular, a two-dimensional (axially symmetrical) model with imprecisely defined parameters is considered. The base of the heat transfer model is given by the interval Fourier equation and supplemented by initial boundary conditions. The phenomenon of cryoprotectant transport (Me2SO) through the extracellular matrix is described by the interval mass transfer equation. The liquidus-tracking (LT) method is used to control the temperature, which avoids the formation of ice regardless of the cooling and warming rates. In the LT process, the temperature decreases/increases gradually during addition/removal of the cryoprotectant, and the articular cartilage remains on or above the liquidus line so that no ice forms, independent of the cooling/warming rate. The discussed problem is solved using the interval finite difference method with the rules of directed interval arithmetic. Examples of numerical computations are presented in the final part of the paper. The obtained results of the numerical simulation are compared with the experimental results, realized for deterministically defined parameters.

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.

Robust numerical solution for the two parameter solute solvent transport model

Prediction of solute and solvent transport in cells is central to developing and testing cryopreservation protocols. As we show here, however, the models used can be difficult to accurately numerically integrate in some key cases, and thus are a challenge to implement when determining the time dependent cell state during cryoprotectant equilibration and cooling. Exact solution techniques exist for overcoming this problem, but their implementation is also challenging: inversion of a nonlinear function is required that negates much of the utility of the approach. This communication describes a simple approach for more robust numerical integration that can be implemented using any numerical differential equation solver, and can facilitate arbitrarily accurate solutions to transport models without the complication of inversion formulae or complicated numerical integration schemes. Further, a simple relevant example of red blood cell equilibration with 40% glycerol is presented with comments on extending the approach to other settings.

Differential Evolution for the Optimization of DMSO-Free Cryoprotectants: Influence of Control Parameters

This study presents the influence of control parameters including population (NP) size, mutation factor (F), crossover (Cr), and four types of differential evolution (DE) algorithms including random, best, local-to-best, and local-to-best with self-adaptive (SA) modification for the purpose of optimizing the compositions of dimethylsufloxide (DMSO)-free cryoprotectants. Post-thaw recovery of Jurkat cells cryopreserved with two DMSO-free cryoprotectants at a cooling rate of 1 °C/min displayed a nonlinear, four-dimensional structure with multiple saddle nodes, which was a suitable training model to tune the control parameters and select the most appropriate type of differential evolution algorithm. Self-adaptive modification presented better performance in terms of optimization accuracy and sensitivity of mutation factor and crossover among the four different types of algorithms tested. Specifically, the classical type of differential evolution algorithm exhibited a wide acceptance to mutation factor and crossover. The optimization performance is more sensitive to mutation than crossover and the optimization accuracy is proportional to the population size. Increasing population size also reduces the sensitivity of the algorithm to the value of the mutation factor and crossover. The analysis of optimization accuracy and convergence speed suggests larger population size with F > 0.7 and Cr > 0.3 are well suited for use with cryopreservation optimization purposes. The tuned differential evolution algorithm is validated through finding global maximums of other two DMSO-free cryoprotectant formulation datasets. The results of these studies can be used to help more efficiently determine the optimal composition of multicomponent DMSO-free cryoprotectants in the future.

Using engineering models to shorten cryoprotectant loading time for the vitrification of articular cartilage

Osteochondral allograft transplantation can treat full thickness cartilage and bone lesions in the knee and other joints, but the lack of widespread articular cartilage banking limits the quantity of cartilage available for size and contour matching. To address the limited availability of cartilage, vitrification can be used to store harvested joint tissues indefinitely. Our group's reported vitrification protocol [Biomaterials 33 (2012) 6061-6068] takes 9.5 h to load cryoprotectants into intact articular cartilage on bone and achieves high cell viability, but further optimization is needed to shorten this protocol for clinical use. Herein, we use engineering models to calculate the spatial and temporal distributions of cryoprotectant concentration, solution vitrifiability, and freezing point for each step of the 9.5-h protocol. We then incorporate the following major design choices for developing a new shorter protocol: (i) all cryoprotectant loading solution concentrations are reduced, (ii) glycerol is removed as a cryoprotectant, and (iii) an equilibration step is introduced to flatten the final cryoprotectant concentration profiles. We also use a new criterion-the spatially and temporally resolved prediction of solution vitrifiability-to assess whether a protocol will be successful instead of requiring that each cryoprotectant individually reaches a certain concentration. A total cryoprotectant loading time of 7 h is targeted, and our new 7-h protocol is predicted to achieve a level of vitrifiability comparable to the proven 9.5-h protocol throughout the cartilage thickness.

Vitrification of in vitro-derived bovine embryos: targeting enhancement of quality by refining technology and standardising procedures

Bovine invitro fertilisation technology has been widely exploited in commercial settings. The majority of invitro-derived cattle embryos are transferred into recipient cows as recently collected (i.e. 'fresh') embryos due to the lack of a reliable cryopreservation method that results in favourable pregnancy rates following transfer of thawed embryos. This is a primary reason for the poor industry uptake of this extreme temperature freezing process. Numerous investigations into vitrification have revealed the importance of rapid cooling and warming rates, enhancing embryo viability after cryopreservation compared with conventional slow freezing. Those studies spawned a considerable assortment of cryovessels and diversity of procedures, delivering variable rates of success, which makes performing vitrification consistently a practical challenge. Hence, further research is required in order to both optimise and standardise vitrification methodology and to design a cryovessel that enables direct transfer of vitrified embryos to recipients after warming. In parallel with improvements in vitrification, it is important to continue to raise the quality of invitro-derived cattle embryos through modifications in laboratory culture techniques. The twin goals of methodology refinement and standardisation, leading to embryo quality enhancement, are each imperative if invitro fertilisation technology is to be adopted in the field.

Foundations of modeling in cryobiology-II: Heat and mass transport in bulk and at cell membrane and ice-liquid interfaces

Modeling coupled heat and mass transport in biological systems is critical to the understanding of cryobiology. In Part I of this series we derived the transport equation and presented a general thermodynamic derivation of the critical components needed to use the transport equation in cryobiology. Here we refine to more cryobiologically relevant instances of a double free-boundary problem with multiple species. In particular, we present the derivation of appropriate mass and heat transport constitutive equations for a system consisting of a cell or tissue with a free external boundary, surrounded by liquid media with an encroaching free solidification front. This model consists of two parts-namely, transport in the "bulk phases" away from boundaries, and interfacial transport. Here we derive the bulk and interfacial mass, energy, and momentum balance equations and present a simplification of transport within membranes to jump conditions across them. We establish the governing equations for this cell/liquid/solid system whose solution in the case of a ternary mixture is explored in Part III of this series.

Numerical solution of inward solidification of a dilute ternary solution towards a semi-permeable spherical cell

Modeling a cell's response to encroaching ice has informed the development of cryopreservation protocols for four decades. It has been well documented that knowledge of the cellular state as a function of media and cooling rate faciliate informed cryopreservation protocol design and explain mechanisms of damage. However, previous efforts have neglected the interaction between solutes and the encroaching ice front and their effects on the cell state. To address this, here we examine the cryobiologically relevant setting of a spherically-symmetric model of a biological cell separated by a ternary fluid mixture from an encroaching solid-liquid interface. The cell and liquid regions contain cell membrane impermeable intracellular and extracellular salts, respectively, a cell membrane permeable solute commonly used in cryopreservation protocols known as a cryoprotective agent (CPA), and water as a membrane permeable solvent. As cooling and solidification proceed the extracellular chemical environment evolves and leads to mass transport across the cell membrane. Consequently, both the solidification front and the cell membrane are free boundaries whose dynamics are coupled through transport processes in the solid, liquid and cell regions. We describe a numerical procedure to solve this coupled free-boundary problem based on a domain transformation and method of lines approach. We also investigate how the thermal and chemical states inside the cell are influenced by different cooling protocols at the external boundary. Finally, we observe that the previously unaccounted-for partial solute rejection at the advancing solid-liquid interface increases the CPA and salt concentrations in the extracellular liquid as a function of the interface speed and segregation coefficients, suggesting that previous model predictions of the cell state during cryopreservation were inaccurate.

Advances in the slow freezing cryopreservation of microencapsulated cells

Over the past few decades, the use of cell microencapsulation technology has been promoted for a wide range of applications as sustained drug delivery systems or as cells containing biosystems for regenerative medicine. However, difficulty in their preservation and storage has limited their availability to healthcare centers. Because the preservation in cryogenic temperatures poses many biological and biophysical challenges and that the technology has not been well understood, the slow cooling cryopreservation, which is the most used technique worldwide, has not given full measure of its full potential application yet. This review will discuss the different steps that should be understood and taken into account to preserve microencapsulated cells by slow freezing in a successful and simple manner. Moreover, it will review the slow freezing preservation of alginate-based microencapsulated cells and discuss some recommendations that the research community may pursue to optimize the preservation of microencapsulated cells, enabling the therapy translate from bench to the clinic.

Droplet based vitrification for cell aggregates: Numerical analysis

Cell aggregates represent the main format of cells existing in vivo and have been widely used as tissue and disease models in vitro. Nevertheless, the preservation of cell aggregates while maintaining their functionalities for off-the-shelf applications is still challenging. Among various preservation methods, droplet-based vitrification exhibits superior advantages for the cryopreservation of cell aggregates; however, the physical mechanisms underlying droplet-based vitrification of cell aggregate using this method remain elusive. To address this issue, we proposed a voronoi model to construct two-dimensional geometric morphologies of cell aggregates and established a coupled physical model to describe the diffusion, heat transfer and crystallization processes during vitrification. Based on these models, we performed a numerical study on the variation and distribution of cryoprotectant (CPA) concentration, temperature and crystallization in cell aggregates during droplet-based vitrification. The results show that although cell membrane is not an obvious barrier in heat transfer, it affects the diffusion of CPA remarkably as a biologic film and thus the following crystallization in cell aggregates. The effective protection of CPA during vitrification occurs during the initial stage of CPA diffusion, thus a longer CPA loading time does not necessarily lead to significant decrease in crystallization, but rather may induce more toxicity to cells.

The effect of additive compounds on glycerol-induced damage to human chondrocytes

High concentrations of cryoprotective agents are required for cryopreservation techniques such as vitrification. Glycerol is a common cryoprotective agent used in cryopreservation protocols but this agent is toxic at high concentrations. This work is an attempt to mitigate the toxic effects of high concentrations of glycerol on intact chondrocytes in human knee articular cartilage from total knee arthroplasty patients by simultaneous exposure to glycerol and a variety of additive compounds. The resulting cell viability in the cartilage samples as measured by membrane integrity staining showed that, in at least one concentration or in combination, all of the tested additive compounds (tetramethylpyrazine, ascorbic acid, chondroitin sulphate, glucosamine sulphate) were able to reduce the deleterious effects of glycerol exposure when examination of membrane integrity took place on a delayed time frame. The use of additive compounds to reduce cryoprotectant toxicity in articular cartilage may help improve cell recovery after cryopreservation.

Anatomical study: comparing the human, sheep and pig knee meniscus

BACKGROUND

Animal models are commonly used in investigating new treatment options for knee joint injuries including injuries to the meniscus. The reliability and applicability of these models to replicate findings in humans depends on determining the most suitable animal proxy. Therefore, this study was designed to compare the wet weight, volume and dimensions of the human meniscus with two commonly used animal models: sheep and pig.

METHODS

Human menisci (n = 6 pairs) were obtained from the knee joints of cadaveric male donors. Sheep menisci (n = 6 pairs) and pig menisci (n = 22 pairs) were obtained from the stifle joints of adult sheep and pigs. Meniscal wet weight, volume and dimensions of the body were measured and compared among the species. Anatomical dimensions included circumference, width, peripheral height, articular height and superior articular length.

RESULTS

The circumference of human menisci (lateral: 84.0 mm, medial: 88.7 mm) was significantly longer than that of sheep (lateral: 50.0 mm, medial: 55.5 mm) and pig (lateral: 66.8 mm, medial: 64.9 mm). The majority of the remaining dimensions of the medial and all of the remaining dimensions of the lateral menisci in sheep showed no statistical difference in comparison to the human menisci. The meniscal weight in pig was significantly larger (lateral: 6.4 g, medial: 5.0 g) than the human (lateral: 4.9 g, medial: 4.4 g) and sheep (lateral: 2.5 g, medial: 2.2 g). Porcine meniscal volume (lateral: 6.5 ml, medial: 5.1 ml) was also larger than the human (lateral: 5.0 ml, medial: 4.5 ml) and sheep (lateral: 2.3 ml, medial: 2.2 ml) menisci. The dimensions measured in the pig meniscus were generally larger than human menisci with statistically significant differences in most categories.

CONCLUSION

Sheep meniscal dimensions more closely matched human meniscal dimensions than the pig meniscal dimensions. This information may help guide the choice of an animal proxy in meniscal research.

Cryoprotectant kinetic analysis of a human articular cartilage vitrification protocol

We recently published a protocol to vitrify human articular cartilage and a method of cryoprotectant removal in preparation for transplantation. The current study's goal was to perform a cryoprotectant kinetic analysis and theoretically shorten the procedure used to vitrify human articular cartilage. First, the loading of the cryoprotectants was modeled using Fick's law of diffusion, and this information was used to predict the kinetics of cryoprotectant efflux after the cartilage sample had been warmed. We hypothesized that diffusion coefficients obtained from the permeation of individual cryoprotectants into porcine articular cartilage could be used to provide a reasonable prediction of the cryoprotectant loading and of the combined cryoprotectant efflux from vitrified human articular cartilage. We tested this hypothesis with experimental efflux measurements. Osteochondral dowels from three patients were vitrified, and after warming, the articular cartilage was immersed in 3 mL X-VIVO at 4 °C in two consecutive solutions, each for 24 h, with the solution osmolality recorded at various times. Measured equilibrium values agreed with theoretical values within a maximum of 15% for all three samples. The results showed that diffusion coefficients for individual cryoprotectants determined from experiments with 2-mm thick porcine cartilage can be used to approximate the rate of efflux of the combined cryoprotectants from vitrified human articular cartilage of similar thickness. Finally, Fick's law of diffusion was used in a computational optimization to shorten the protocol with the constraint of maintaining the theoretical minimum cryoprotectant concentration needed to achieve vitrification. The learning provided by this study will enable future improvements in tissue vitrification.

Toxicity Minimized Cryoprotectant Addition and Removal Procedures for Adherent Endothelial Cells

Ice-free cryopreservation, known as vitrification, is an appealing approach for banking of adherent cells and tissues because it prevents dissociation and morphological damage that may result from ice crystal formation. However, current vitrification methods are often limited by the cytotoxicity of the concentrated cryoprotective agent (CPA) solutions that are required to suppress ice formation. Recently, we described a mathematical strategy for identifying minimally toxic CPA equilibration procedures based on the minimization of a toxicity cost function. Here we provide direct experimental support for the feasibility of these methods when applied to adherent endothelial cells. We first developed a concentration- and temperature-dependent toxicity cost function by exposing the cells to a range of glycerol concentrations at 21°C and 37°C, and fitting the resulting viability data to a first order cell death model. This cost function was then numerically minimized in our state constrained optimization routine to determine addition and removal procedures for 17 molal (mol/kg water) glycerol solutions. Using these predicted optimal procedures, we obtained 81% recovery after exposure to vitrification solutions, as well as successful vitrification with the relatively slow cooling and warming rates of 50°C/min and 130°C/min. In comparison, conventional multistep CPA equilibration procedures resulted in much lower cell yields of about 10%. Our results demonstrate the potential for rational design of minimally toxic vitrification procedures and pave the way for extension of our optimization approach to other adherent cell types as well as more complex systems such as tissues and organs.

Fatigue damage to pig erythrocytes during repeated swelling and shrinkage

During the removal of cryoprotectants from cryopreserved-thawed blood with the dialysis-based or dilution-filtration method, due to the change in the extracellular osmolality, erythrocytes usually undergo repeated swelling and shrinkage. However, the erythrocyte fatigue damage induced by this repeated volume change has not yet been studied. In this work, by successively loading hypotonic and hypertonic solutions, we mimicked the repeated swelling and shrinkage of pig erythrocytes and then examined the effect of the number of cycle loops on the steady-state volume and the mortality of the pig erythrocytes. The results suggest that because of cell leakage in the swelling process, the steady-state volume of the pig erythrocytes after one cycle is smaller than the volume before the cycle, even though the cell performs a self-protective regulatory procedure. If the number of cycle loops is increased, the repeated swelling and shrinkage will cause a continuous decrease in the steady-state volume, and the ability of the pig erythrocytes to resist osmotic damage will decrease; as a result, the mortality of the pig erythrocytes increases as the number of cycle loops increases. The viability of the cells is also affected by the hypotonic and isotonic processing times: a short processing time may contribute to a decrease in the mortality of the pig erythrocytes. This work is of significance to optimizing the process of removing cryoprotectants.

Advances in cell encapsulation technology and its application in drug delivery

INTRODUCTION

Cell encapsulation technology has improved enormously since it was proposed 50 years ago. The advantages offered over other alternative systems, such as the prevention of repetitive drug administration, have triggered the use of this technology in multiple therapeutic applications.

AREAS COVERED

In this article, improvements in cell encapsulation technology and strategies to overcome the drawbacks that prevent its use in the clinic have been summarized and discussed. Different studies and clinical trials that have been performed in several therapeutic applications have also been described.

EXPERT OPINION

The authors believe that the future translation of this technology from bench to bedside requires the optimization of diverse aspects: i) biosafety, controlling and monitoring cell viability; ii) biocompatibility, reducing pericapsular fibrotic growth and hypoxia suffered by the graft; iii) control over drug delivery; iv) and the final scale up. On the other hand, an area that deserves more attention is the cryopreservation of encapsulated cells as this will facilitate the arrival of these biosystems to the clinic.

A model for predicting the permeation of dimethyl sulfoxide into articular cartilage, and its application to the liquidus-tracking method

Long-term storage of articular cartilage (AC) has excited great interest due to the practical surgical significance of this tissue. The liquidus-tracking (LT) method developed by Pegg et al. (2006) [29] for vitreous preservation of AC achieved reasonable survival of post-warming chondrocytes in situ, but the design of the entire procedure was more dependent on trial and error. Mathematical modeling would help to better understand the LT process, and thereby make possible improvements to attain higher cell survival. Mass transfer plays a dominant role in the LT process. In the present study, a diffusion model based on the free-volume theory and the Flory-Huggins thermodynamics theory was developed to predict the permeation of dimethyl sulfoxide (Me2SO) into AC. A comparison between the predicted mean concentration of Me2SO in the AC disc and the experimental data over wide temperature and concentration ranges [-30 to 37 °C, 10 to 64.5% (w/w)] shows that the developed model can accurately describe the permeation of Me2SO into AC [coefficient of determination (R(2)): 0.951-1.000, mean relative error (MRE): 0.8-12.8%]. With this model, the spatial and temporal distribution of Me2SO in the AC disc during a loading/unloading process can be obtained. Application of the model to Pegg et al.'s LT procedure revealed that the liquidus line is virtually not followed for the center part of the AC disc. The presently developed model will be a useful tool in the analysis and design of the LT method for vitreous preservation of AC.

Cryopreservation of articular cartilage

Cryopreservation has numerous practical applications in medicine, biotechnology, agriculture, forestry, aquaculture and biodiversity conservation, with huge potentials for biological cell and tissue banking. A specific tissue of interest for cryopreservation is the articular cartilage of the human knee joint for two major reasons: (1) clinically, there exists an untapped potential for cryopreserved cartilage to be used in surgical repair/reconstruction/replacement of injured joints because of the limited availability of fresh donor tissue and, (2) scientifically, successful cryopreservation of cartilage, an avascular tissue with only one cell type, is considered a stepping stone for transition from biobanking cell suspensions and small tissue slices to larger and more complicated tissues. For more than 50years, a great deal of effort has been directed toward understanding and overcoming the challenges of cartilage preservation. In this article, we focus mainly on studies that led to the finding that vitrification is an appropriate approach toward successful preservation of cartilage. This is followed by a review of the studies on the main challenges of vitrification, i.e. toxicity and diffusion, and the novel approaches to overcome these challenges such as liquidus tracking, diffusion modeling, and cryoprotective agent cocktails, which have resulted in the recent advancements in the field.

A numerical study on distributions during cryoprotectant loading caused by laminar flow in a microchannel

In this work, we conduct a computational study on the loading of cryoprotective agents into cells in preparation for cryopreservation. The advantages of microfluidics in cryopreserving cells include control of fluid flow parameters for reliable cryoprotectant loading and reproducible streamlined processing of samples. A 0.25 m long, three inlet T-junction microchannel serves as an idealized environment for this process. The flow field and concentration distribution are determined from a computational fluid dynamics study and cells are tracked as inert particles in a Lagrangian frame. These particles are not confined to streamlines but can migrate laterally due to the Segre-Sildeberg effect for particles in a shear flow. During this tracking, the local concentration field surrounding the cell is monitored. This data are used as input into the Kedem-Katchalsky equations to numerically study passive solute transport across the cell membrane. As a result of the laminar flow, each cell has a unique pathline in the flow field resulting in different residence times and a unique external concentration field along its path. However, in most previous studies, the effect of a spatially varying concentration field on the transport across the cell membrane is ignored. The dynamics of this process are investigated for a population of cells released from the inlet. Using dimensional analysis, we find a governing parameter α, which is the ratio of the time scale for membrane transport to the average residence time in the channel. For [Formula: see text], cryoprotectant loading is completed to within 5% of the target concentration for all of the cells. However, for [Formula: see text], we find the population of cells does not achieve complete loading and there is a distribution of intracellular cryoprotective agent concentration amongst the population. Further increasing α beyond a value of 2 leads to negligible cryoprotectant loading. These simulations on populations of cells may lead to improved microfluidic cryopreservation protocols where more consistent cryoprotective agent loading and freezing can be achieved, thus increasing cell survival.