Concentration dependence of the cell membrane permeability to cryoprotectant and water and implications for design of methods for post-thaw washing of human erythrocytes

Cryopreservation of organoids: Strategies, innovation, and future prospects

Organoid technology has demonstrated unique advantages in multidisciplinary fields such as disease research, tumor drug sensitivity, clinical immunity, drug toxicology, and regenerative medicine. It will become the most promising research tool in translational research. However, the long preparation time of organoids and the lack of high-quality cryopreservation methods limit the further application of organoids. Although the high-quality cryopreservation of small-volume biological samples such as cells and embryos has been successfully achieved, the existing cryopreservation methods for organoids still face many bottlenecks. In recent years, with the development of materials science, cryobiology, and interdisciplinary research, many new materials and methods have been applied to cryopreservation. Several new cryopreservation methods have emerged, such as cryoprotectants (CPAs) of natural origin, ice-controlled biomaterials, and rapid rewarming methods. The introduction of these technologies has expanded the research scope of cryopreservation of organoids, provided new approaches and methods for cryopreservation of organoids, and is expected to break through the current technical bottleneck of cryopreservation of organoids. This paper reviews the progress of cryopreservation of organoids in recent years from three aspects: damage factors of cryopreservation of organoids, new protective agents and loading methods, and new technologies of cryopreservation and rewarming.

Silymarin mediated osmotic responses and damage in HepG2 cell suspensions and monolayers

Maintenance of cells within a volume range compatible with their functional integrity is a critical determinant of cell survival after cryopreservation, and quantifying this osmotically induced damage is a part of the rational design of improved cryopreservation protocols. The degree that cells tolerate osmotic stress significantly impacts applicable cryoprotocols, but there has been little research on the time dependence of this osmotic stress. Additionally, the flavonoid silymarin has been shown to be hepatoprotective. Therefore, here we test the hypotheses that osmotic damage is time-dependent and that flavonoid inclusion reduces osmotic damage. In our first experiment, cells were exposed to a series of anisosmotic solutions of graded hypo- and hypertonicity for 10-40 min, resulting in a conclusion that osmotically induced damage is time dependent. In the next experiment, adherent cells preincubated with silymarin at the concentration of 10-4 mol/L and 10-5 mol/L showed a significant increase in cell proliferation and metabolic activity after osmotic stress compared to untreated matched controls. For instance, when adherent cells preincubated with 10-5 mol/L silymarin were tested, resistance to osmotic damage and a significant increase (15%) in membrane integrity was observed in hypo-osmotic media and a 22% increase in hyperosmotic conditions. Similarly, significant protection from osmotic damage was observed in suspended HepG2 cells in the presence of silymarin. Our study concludes that osmotic damage is time dependent, and the addition of silymarin leads to elevated resistance to osmotic stress and a potential increase in the cryosurvival of HepG2 cells.

An impedance based microfluidic sensor for evaluation of individual red blood cell solute permeability

-In this paper, we investigate a microfluidic based sensing device for cell membrane permeability measurements in real time with applications in rapid assessment of red blood cell (RBC) quality at the individual cell level. The microfluidic chip was designed with unique abilities to line up the RBCs in the centerline of the microchannel using positive dielectrophoresis (p-DEP) forces, rapid mixing of RBCs with various media (e.g. containing permeating or nonpermeating solutes) injected from different inlets to achieve high mixing efficiency. The chip detects the impedance values of the RBCs within 0.19 s from the start of mixing with other media, at ten electrodes along the length of the channel and enables time series measurements of volume change of individual cell caused by cell osmosis in anisosmotic fluids over a 0.8 s postmixing timespan. This technique enables estimating water permeability of individual cell accurately. Here we first present confirmation of a linear voltage-diameter relationship in polystyrene bead standards. Next, we show that under equilibrium conditions, the voltage-volume relationship in rat red blood cells (RBCs) is linear, corresponding to previously published Boyle van 't Hoff plots. Using rat cells as a model for human, we present the first measurement of water permeability in individual red blood cells and confirm that these data align with previously published population level values for human RBC. Finally, we present preliminary evidence for possible application of our device to identify individual RBCs infected with Plasmodium falciparum malaria parasites. Future developments using this device will address the use of whole blood with non-homogenous cell populations, a task currently performed by clinical Coulter counters.

Comparison of dilute and nondilute osmotic equilibrium models for erythrocytes

Successful cryopreservation requires the addition of cryoprotective agents (CPAs). The addition of permeating CPAs, such as glycerol, is associated with some risk to the cells and tissues. These risks are both related to the CPA themselves (CPA toxicity) and to the volume response of the cell (osmotic damage). To minimize the potential for damage during cryopreservation, mathematical models are often employed to understand the interactions between protocols and cell volume responses. In the literature, this volume response is usually captured using ideal and dilute approximations of chemical potential and osmolality, an approach that has been called into question for cells in high concentrations of CPAs. To address this, the relevance of non-ideal and non-dilute models has been explored in a number of cell types in the presence of permeating CPAs. However, it has not been explored in erythrocytes, which have a cytosolic hemoglobin content of more than 20% by volume and are cryopreserved in 40% glycerol. Because hemoglobin has been suggested to be a highly non-ideal solute, if the non-ideal and non-dilute transport model is relevant to any cells, it should be relevant to erythrocytes. Here we investigate the use, and accuracy, of both the dilute and non-dilute models in predicting cell volume changes during CPA equilibration in erythrocytes, and demonstrate that using published values for the non-ideal and non-dilute model, applied to erythrocytes, leads to model predictions inconsistent with experimental data, whereas dilute approximations align well with experimental data.

Exogenous Melatonin Ameliorates the Negative Effect of Osmotic Stress in Human and Bovine Ovarian Stromal Cells

Ovarian tissue cryopreservation transplantation (OTCT) is the most flexible option to preserve fertility in women and children with cancer. However, OTCT is associated with follicle loss and an accompanying short lifespan of the grafts. Cryopreservation-induced damage could be due to cryoprotective agent (CPA) toxicity and osmotic shock. Therefore, one way to avoid this damage is to maintain the cell volume within osmotic tolerance limits (OTLs). Here, we aimed to determine, for the first time, the OTLs of ovarian stromal cells (OSCs) and their relationship with reactive oxygen species (ROS) and mitochondrial respiratory chain activity (MRCA) of OSCs. We evaluated the effect of an optimal dose of melatonin on OTLs, viability, MRCA, ROS and total antioxidant capacity (TAC) of both human and bovine OSCs in plated and suspended cells. The OTLs of OSCs were between 200 and 375 mOsm/kg in bovine and between 150 and 500 mOsm/kg in human. Melatonin expands OTLs of OSCs. Furthermore, melatonin significantly reduced ROS and improved TAC, MRCA and viability. Due to the narrow osmotic window of OSCs, it is important to optimize the current protocols of OTCT to maintain enough alive stromal cells, which are necessary for follicle development and graft longevity. The addition of melatonin is a promising strategy for improved cryopreservation media.

Effect of cryoprotectant concentration on bovine oocyte permeability and comparison of two membrane permeability modelling approaches

The plasma membrane permeability to water and cryoprotectant (CPA) significantly impacts vitrification efficiency of bovine oocytes. Our study was designed to determine the concentration-dependent permeability characteristics for immature (GV) and mature (MII) bovine oocytes in the presence of ethylene glycol (EG) and dimethyl sulphoxide (Me₂SO), and to compare two different modeling approaches: the two parameter (2P) model and a nondilute transport model. Membrane permeability parameters were determined by consecutively exposing oocytes to increasing concentrations of Me₂SO or EG. Higher water permeability was observed for MII oocytes than GV oocytes in the presence of both Me₂SO and EG, and in all cases the water permeability was observed to decrease as CPA concentration increased. At high CPA concentrations, the CPA permeability was similar for Me₂SO and EG, for both MII and GV oocytes, but at low concentrations the EG permeability of GV oocytes was substantially higher. Predictions of cell volume changes during CPA addition and removal indicate that accounting for the concentration dependence of permeability only has a modest effect, but there were substantial differences between the 2P model and the nondilute model during CPA removal, which may have implications for design of improved methods for bovine oocyte vitrification.

Insights into the model of non-perfect osmometer cells for cryopreservation: A parametric sweep analysis

Recently, a mathematical model able to describe the non-perfect osmotic behavior of cells during cryopreservation was proposed. The model improves the two-parameter formalism typically adopted in cryopreservation literature by allowing the transmembrane permeation of ions/salt, through the temporary opening of mechanosensitive channels whenever membrane stretching occurs: cells can reach an equilibrium volume different from the initial one, when isotonic conditions are re-established after contacting with impermeant or permeant solutes, such as sucrose or a cryoprotectant agent like dimethyl sulfoxide, respectively. Although the model was conceived as a conservative development of the two-parameter formalism to avoid over-parameterization, a complex picture of the system emerges. To describe this, first an appropriate non-dimensional version of the model equations is derived. Then, a parametric sweep analysis is performed and discussed to highlight the features of the novel model in comparison with the two-parameter formalism: the conditions by which the first reduces to the second are identified. Only equilibrium equations with impermeant sucrose may be analytically derived from the model: their validity is here extended much more than originally assumed. When permeant dimethyl sulfoxide comes into play, the temporary opening of mechanosensitive channels is difficult to predict and prevents the derivation of the equilibrium equations: in this case, a numerical integration of system dynamics up to steady state is required to determine the cell volume at equilibrium. In conclusion, cell volume at equilibrium depends on the position of the temporal window of mechanosensitive channels opening, which, in general, is a complex function of model parameters and operating conditions.

Quantitative study of unsaturated transport of glycerol through aquaglyceroporin that has high affinity for glycerol

The structures of several aquaglyceroporins have been resolved to atomic resolution showing two or more glycerols bound inside a channel and confirming a glycerol-facilitator's affinity for its substrate glycerol. However, the kinetics data of glycerol transport experiments all point to unsaturated transport that is characteristic of low substrate affinity in terms of the Michaelis-Menten kinetics. In this article, we present an in silico-in vitro research focused on AQP3, one of the human aquaglyceroporins that is natively expressed in the abundantly available erythrocytes. We conducted 2.1 μs in silico simulations of AQP3 embedded in a model erythrocyte membrane with intracellular-extracellular asymmetries in leaflet lipid compositions and compartment salt ions. From the equilibrium molecular dynamics (MD) simulations, we elucidated the mechanism of glycerol transport at high substrate concentrations. From the steered MD simulations, we computed the Gibbs free-energy profile throughout the AQP3 channel. From the free-energy profile, we quantified the kinetics of glycerol transport that is unsaturated due to glycerol-glycerol interactions mediated by AQP3 resulting in the concerted movement of two glycerol molecules for the transport of one glycerol molecule across the cell membrane. We conducted in vitro experiments on glycerol uptake into human erythrocytes for a wide range of substrate concentrations and various temperatures. The experimental data quantitatively validated our theoretical-computational conclusions on the unsaturated glycerol transport through AQP3 that has high affinity for glycerol.

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.

Administration of glycerol-based formulations in sheep results in similar ovulation rate to eCG but red blood cell indices may be affected

BACKGROUND

The objective of this study was to investigate the metabolic and osmotic effects of different doses of glycerol or a glycerol - propylene glycol mixture in Sarda sheep with the aim to identify those able to beneficially modify ewe's metabolic status without harmful changes in red blood cell (RBC) indices. Thereafter, the selected doses were tested for their effects on ewe's ovarian activity during an induced follicular phase and compared to the effects of a hormonal treatment with equine chorionic gonadotrophin (eCG).

RESULTS

Glycerol was administered alone (G groups: 90% glycerol and 10% water; % v/v) or in combination with propylene glycol (M groups: 70% glycerol, 20% propylene glycol, 10% water; % v/v). Treatments were formulated to provide 100, 75, 50 and 25% of the amount of energy supplied in previous experiments. Obtained results showed that the formulations G75 and M75 (22.5 and 18.2% on DM basis, respectively) induce metabolic changes comparable to those induced by M100. The latter dose has been already evaluated for its effects on sheep metabolism and reproductive performance. However, with these high doses, plasma osmolality increased significantly, and RBC indices showed significant alterations. The low dose groups (G25 and M25, 8.6 and 6.9% on DM basis, respectively) did not show any alterations in plasma osmolality and RBC indices, but the metabolic milieu differed markedly from that of M100. Between the medium dose groups, M50 (12.9% on DM basis) showed a more comparable milieu to M100 than G50 (15.9% on DM basis) and no RBC alterations. Therefore, M75, G75 and M50 doses were tested for their effect on ovarian functions and proved to be equally effective as eCG.

CONCLUSION

The results of the present study evidenced an alteration of RBC indices, and possibly of their functions, as a side effect of glycerol administration at high doses in the diet of ewes. Therefore, protocols foreseeing the administration of glycerol should be tested for their effects on RBC indices and functions. In general terms, the medium dose of the glucogenic mixture (12.9% of dietary DM on offer) should be preferred.

Implications of variability in cell membrane permeability for design of methods to remove glycerol from frozen-thawed erythrocytes

In North America, red blood cells (RBCs) are currently cryopreserved in a solution of 40% glycerol. While glycerol is not inherently toxic to humans, it must be removed prior to transfusion to prevent intravascular osmotic hemolysis. The current deglycerolization procedure requires about 45 min per RBC unit. We previously presented predictions suggesting that glycerol could be safely removed from RBCs in less than 1 min. However, experimental evaluation of these methods resulted in much higher hemolysis than expected. Here we extend our previous study by considering both concentration-dependence of permeability and variability in permeability values in the mathematical optimization algorithm. To establish a model for the concentration dependence of glycerol permeability, we combined literature data with new measurements of permeability in the presence of 40% glycerol. To account for cell-dependent variability we scaled the concentration-dependent permeability model to define a permeability range for optimization. Methods designed using a range extending to 50% of the model-predicted glycerol permeability had a duration of less than 3 min and resulted in hemolysis ranging from 34% to 83%; hemolysis values were highly dependent on the blood donor. Extending the permeability range to 5% of the model-predicted value yielded a 30 min method that resulted in an average hemolysis of 12%. Our results suggest high variability in the glycerol permeability between donors and within a population of cells from the same donor. Such variability has broad implications for design of methods for equilibration of cells with cryoprotectants.

A microfluidic approach for synchronous and nondestructive study of the permeability of multiple oocytes

Investigation of oocyte membrane permeability plays a crucial role in fertility preservation, reproductive medicine, and reproductive pharmacology. However, the commonly used methods have disadvantages such as high time consumption, low efficiency, and cumbersome data processing. In addition, the developmental potential of oocytes after measurement has not been fully validated in previous studies. Moreover, oocytes can only maintain their best status in vitro within a very limited time. To address these limitations, we developed a novel multichannel microfluidic chip with newly designed micropillars that provide feasible and repeatable oocyte capture. The osmotic responses of three oocytes at different or the same cryoprotectant (CPA) concentrations were measured simultaneously, which greatly improved the measurement efficiency. Importantly, the CPA concentration dependence of mouse oocyte membrane permeability was found. Moreover, a neural network algorithm was employed to improve the efficiency and accuracy of data processing. Furthermore, analysis of fertilization and embryo transfer after perfusion indicated that the microfluidic approach does not damage the developmental potential of oocytes. In brief, we report a new method based on a multichannel microfluidic chip that enables synchronous and nondestructive measurement of the permeability of multiple oocytes.

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.

Factors Affecting the Membrane Permeability Barrier Function of Cells during Preservation Technologies

Cellular membranes are exposed to extreme conditions during the processing steps involved in cryopreservation (and freeze-drying) of cells. The first processing step involves adding protective agents. Exposing cells to protective agents causes fluxes of both water and solutes (i.e., permeating cryoprotective agents) across the cellular membrane, resulting in cell volume changes and possibly osmotic stress. In addition, protective molecules may interact with lipids, which may lead to membrane structural changes and permeabilization. After loading with protective agents, subsequent freezing exposes cells to severe osmotic and mechanical stresses, caused by extra and/or intracellular ice formation and a drastically increased solute concentration in the unfrozen fraction. Furthermore, cellular membranes undergo thermotropic and lyotropic phase transitions during cooling and freezing, which drastically alter the membrane permeability and its barrier function. In this article, it is shown that membrane permeability to water and solutes is dependent on the temperature, medium osmolality, types of solutes present, cell hydration level, and absence or presence of ice. Freezing most drastically alters the membrane permeability barrier function, which is reflected as a change in the activation energy for water transport. In addition, membranes become temporarily leaky during freezing-induced fluid-to-gel membrane phase transitions, resulting in the uptake of impermeable solutes.

Single-channel permeability and glycerol affinity of human aquaglyceroporin AQP3

For its fundamental relevance, transport of water and glycerol across the erythrocyte membrane has long been investigated before and after the discovery of aquaporins (AQPs), the membrane proteins responsible for water and glycerol transport. AQP1 is abundantly expressed in the human erythrocyte for maintaining its hydrohomeostasis where AQP3 is also expressed (at a level ~30-folds lower than AQP1) facilitating glycerol transport. This research is focused on two of the remaining questions: How permeable is AQP3 to water? What is the glycerol-AQP3 affinity under near-physiological conditions? Through atomistic modelling and large-scale simulations, we found that AQP3 is two to three times more permeable to water than AQP1 and that the glycerol-AQP3 affinity is approximately 500/M. Using these computed values along with the data from the latest literature on AQP1 and on erythrocyte proteomics, we estimated the water and glycerol transport rates across the membrane of an entire erythrocyte. We used these rates to predict the time courses of erythrocyte swelling-shrinking in response to inward and outward osmotic gradients. Experimentally, we monitored the time course of human erythrocytes when subject to an osmotic or glycerol gradient with light scattering in a stopped-flow spectrometer. We observed close agreement between the experimentally measured and the computationally predicted time courses of erythrocytes, which corroborated our computational conclusions on the AQP3 water-permeability and the glycerol-AQP3 affinity.