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

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.

Loading equine oocytes with cryoprotective agents captured with a finite element method model

Cryopreservation can be used to store equine oocytes for extended periods so that they can be used in artificial reproduction technologies at a desired time point. It requires use of cryoprotective agents (CPAs) to protect the oocytes against freezing injury. The intracellular introduction of CPAs, however, may cause irreversible osmotic damage. The response of cells exposed to CPA solutions is governed by the permeability of the cellular membrane towards water and the CPAs. In this study, a mathematical mass transport model describing the permeation of water and CPAs across an oocyte membrane was used to simulate oocyte volume responses and concomitant intracellular CPA concentrations during the exposure of oocytes to CPA solutions. The results of the analytical simulations were subsequently used to develop a phenomenological finite element method (FEM) continuum model to capture the response of oocytes exposed to CPA solutions with spatial information. FEM simulations were used to depict spatial differences in CPA concentration during CPA permeation, namely at locations near the membrane surface and towards the middle of the cell, and to capture corresponding changes in deformation and hydrostatic pressure. FEM simulations of the multiple processes occurring during CPA loading of oocytes are a valuable tool to increase our understanding of the mechanisms underlying cryopreservation outcome.

Investigating the Solubility and Activity of a Novel Class of Ice Recrystallization Inhibitors

O-aryl-β-d-glucosides and N-alkyl-d-gluconamides are two classes of effective ice recrystallization inhibitors (IRIs), however their solubilities limit their use in cryopreservation applications. Herein, we have synthesized and assessed phosphonate analogues of small-molecule IRIs as a method to improve their chemical and physical properties. Four sodium phosphonate compounds 4–7 were synthesized and exhibited high solubilities greater than 200 mM. Their IRI activity was evaluated using the splat cooling assay and only the sodium phosphonate derivatives of α-methyl-d-glucoside (5-Na) and N-octyl-d-gluconamide (7-Na) exhibited an IC50 value less than 30 mM. It was found that the addition of a polar sodium phosphonate group to the alkyl gluconamide (1) and aryl glucoside (2) structure decreased its IRI activity, indicating the importance of a delicate hydrophobic/hydrophilic balance within these compounds. The evaluation of various cation-phosphonate pairs was studied and revealed the IRI activity of ammonium and its ability to modulate the IRI activity of its paired anion. A preliminary cytotoxicity study was also performed in a HepG2 cell line and phosphonate analogues were found to have relatively low cytotoxicity. As such, we present phosphonate small-molecule carbohydrates as a biocompatible novel class of IRIs with high solubilities and moderate-to-high IRI activities.

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.