Non-ideal solution thermodynamics of cytoplasm

Temperature Dependence of Membrane Permeability Parameters for Five Cell Types Using Nonideal Thermodynamic Assumptions to Mathematically Model Cryopreservation Protocols

Cryopreservation is the process of preserving biological matter at subzero temperatures for long-term storage. During cryopreservation, cells are susceptible to various injuries that can be mitigated by controlling the cooling and warming profiles and cryoprotective agent (CPA) addition and removal procedures. Mathematical modeling of the changing cell volume at different temperatures can greatly reduce the experiments needed to optimize cryopreservation protocols. Such mathematical modeling requires as inputs the cell membrane permeabilities to water and CPA and the osmotically inactive fraction of the cell. Since the intra- and extracellular solutions are generally thermodynamically nonideal, our group has been incorporating the osmotic virial equation to model the solution thermodynamics that underlie the cell volume change equations, adding the second and third osmotic virial coefficients of the grouped intracellular solute to the cell osmotic parameters that must be measured. In our previous work, we reported methods to obtain cell osmotic parameters at room temperature by fitting experimental cell volume kinetic data with equations that incorporated nonideal solution thermodynamics assumptions. Since the relevant cell volume excursions occur at different temperatures, the temperature dependence of the osmotic parameters plays an important role. In this work, we present a new two-part fitting method to obtain five cell-type-specific parameters (water permeability, dimethyl sulfoxide permeability, osmotically inactive fraction, and the second and third osmotic virial coefficients of the intracellular solution) from experimental measurements of equilibrium cell volume and cell volume as a function of time at room temperature and 0 °C for five cell types, namely, human umbilical vein endothelial cells (HUVECs), H9c2 rat myoblasts, porcine corneal endothelial cells (PCECs), the Jurkat T-lymphocyte cell line, and human cerebral microvascular endothelial cells (hCMECs/D3 cell line). The fitting method in this work is based on both equilibrium and kinetic cell volume data, enabling us to solve some technical challenges and expand our previously reported measurement technique to 0 °C. Finally, we use the measured parameters to model the cell volume changes for a HUVEC cryopreservation protocol to demonstrate the impact of the nonideal thermodynamic assumptions on predicting the changing cell volume during freezing and thawing.

Predicting freezing points of ternary salt solutions with the multisolute osmotic virial equation

Previously, the multisolute osmotic virial equation with the combining rules of Elliott et al. has been shown to make accurate predictions for multisolute solutions with only single-solute osmotic virial coefficients as inputs. The original combining rules take the form of an arithmetic average for the second-order mixed coefficients and a geometric average for the third-order mixed coefficients. Recently, we derived generalized combining rules from a first principles solution theory, where all mixed coefficients could be expressed as arithmetic averages of suitable binary coefficients. In this work, we empirically extended the new model to account for electrolyte effects, including solute dissociation, and demonstrated its usefulness for calculating the properties of multielectrolyte solutions. First, the osmotic virial coefficients of 31 common salts in water were tabulated based on the available freezing point depression (FPD) data. This was achieved by polynomial fitting, where the degree of the polynomial was determined using a special criterion that accounts for the confidence intervals of the coefficients. Then, the multisolute model was used to predict the FPD of 11 ternary electrolyte solutions. Furthermore, models with the new combining rules and the original combining rules of Elliott et al. were compared using both mole fraction and molality as concentration units. We find that the mole-fraction-based model with the new combining rules performs the best and that the results agree well with independent experimental measurements with an all-system root-mean-square error of 0.24 osmoles/kg (0.45 °C) and close to zero mean bias for the entire dataset (371 data points).

Permeability and Osmotic Parameters of Human Umbilical Vein Endothelial Cells and H9C2 Cells under Non-ideal Thermodynamic Assumptions: A Novel Iterative Fitting Method

Cryopreservation is the use of very low subzero temperatures to preserve cells and tissues for later use. This is achieved by controlled cooling in the presence of cryoprotectants that moderate the amount of ice formed. Mathematical modeling of the cryopreservation process is a useful tool to investigate the different variables that affect the results of this process. The changing cell volume during cryopreservation can be modeled using cell membrane water and cryoprotectant permeabilities and the osmotically inactive fraction of the intracellular contents. These three cell-specific parameters have been found previously for different cell types under ideal and dilute assumptions, but biological solutions at subzero temperatures are far from ideal and dilute, especially when cryoprotectants are included. In this work, the osmotic virial equation is used to model the changing cell volume under non-ideal assumptions, and the intracellular environment is described using the grouped solute, which consists of all impermeant intracellular solutes grouped together, leading to two additional cell-specific parameters, the second and third osmotic virial coefficients of the grouped solute. Herein, we present a novel fitting method to efficiently determine these five cell-specific parameters by fitting kinetic cell volume data under non-ideal assumptions and report the results of applying this method to obtain the parameters for two cell types: human umbilical vein endothelial cells and H9C2 rat myoblasts.

Gibbsian Surface Thermodynamics

Gibbsian composite-system thermodynamics is the framework governing the equilibrium of composite systems, including systems that at equilibrium have more than one value of pressure because of the action of surface tension, semipermeable membranes, or fields, and thus cannot be treated as simple systems. J. W. Gibbs's paper that lays out composite-system thermodynamics, "On the Equilibrium of Heterogeneous Substances", communicated in two parts in 1876 and 1878, is widely regarded as one of the most important pieces of scientific literature of its century. Many scientists adopted and stressed the importance of Gibbsian thermodynamics. In 1960, H. B. Callen wrote a textbook that made Gibbsian composite-system thermodynamics more accessible to thermodynamicists. However, Callen's book left out Gibbs's work on curved fluid interfaces and did not treat the complicated nonideal systems of interest to today's thermodynamicists. In this Feature Article, I have attempted to convey in a comprehensive manner the framework of Gibbsian composite-system thermodynamics including in detail the treatment of systems with interface effects and with nonideal, multicomponent phases. This work lays out the relationships between important equilibrium equations including the following: the Gibbs-Duhem equation, the Gibbs adsorption equation, the Young-Laplace equation, the Young equation, the Cassie-Baxter equation, the Wenzel equation, the Kelvin equation, the Gibbs-Thompson equation, and the Ostwald-Freundlich equation, including nonideal and multicomponent forms. Equations of state that are often useful for Gibbsian composite-system thermodynamics are reviewed including adsorption isotherms and our own work on two semiempirical equations of state: the Elliott et al. form of the osmotic virial equation and the Shardt-Elliott-Connors-Wright equation for the temperature and composition dependence of surface tension. I summarize the work of our group developing Gibbisan composite-system thermodynamics including new equations for such things as the curvature-induced depression of the eutectic temperature or the removal of azeotropes by nanoscale fluid interface curvature. Gibbsian composite-system thermodynamics has broad applications in biotechnology, nanostructured materials, surface textures and coatings, microfluidics, nanoscience, atmospheric and environmental physics, among others.

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.

Measuring the effects of macromolecular crowding on antibody function with biolayer interferometry

Biotherapeutic proteins are commonly dosed at high concentrations into the blood, which is an inherently complex, crowded solution with substantial protein content. The effects of macromolecular crowding may lead to an appreciable level of non-specific hetero-association in this physiological environment. Therefore, developing a method to characterize the diverse consequences of non-specific interactions between proteins under such non-ideal, crowded conditions, which deviate substantially from those commonly employed for in vitro characterization, is vital to achieving a more complete picture of antibody function in a biological context. In this study, we investigated non-specific interactions between human serum albumin (HSA) and two monoclonal antibodies (mAbs) by static light scattering and determined these interactions are both ionic strength-dependent and mAb-dependent. Using biolayer interferometry (BLI), we assessed the effect of HSA on antigen binding by mAbs, demonstrating that these non-specific interactions have a functional impact on mAb:antigen interactions, particularly at low ionic strength. While this effect is mitigated at physiological ionic strength, our in vitro data support the notion that HSA in the blood may lead to non-specific interactions with mAbs in vivo, with a potential impact on their interactions with antigen. Furthermore, the BLI method offers a high-throughput advantage compared to orthogonal techniques such as analytical ultracentrifugation and is amenable to a greater variety of solution conditions compared to nuclear magnetic resonance spectroscopy. Our study demonstrates that BLI is a viable technology for examining the impact of non-specific interactions on specific biologically relevant interactions, providing a direct method to assess binding events in crowded conditions.

Improved Cryopreservation of Human Umbilical Vein Endothelial Cells: A Systematic Approach

Cryopreservation of human umbilical vein endothelial cells (HUVECs) facilitated their commercial availability for use in vascular biology, tissue engineering and drug delivery research; however, the key variables in HUVEC cryopreservation have not been comprehensively studied. HUVECs are typically cryopreserved by cooling at 1 °C/min in the presence of 10% dimethyl sulfoxide (DMSO). We applied interrupted slow cooling (graded freezing) and interrupted rapid cooling with a hold time (two-step freezing) to identify where in the cooling process cryoinjury to HUVECs occurs. We found that linear cooling at 1 °C/min resulted in higher membrane integrities than linear cooling at 0.2 °C/min or nonlinear two-step freezing. DMSO addition procedures and compositions were also investigated. By combining hydroxyethyl starch with DMSO, HUVEC viability after cryopreservation was improved compared to measured viabilities of commercially available cryopreserved HUVECs and viabilities for HUVEC cryopreservation studies reported in the literature. Furthermore, HUVECs cryopreserved using our improved procedure showed high tube forming capability in a post-thaw angiogenesis assay, a standard indicator of endothelial cell function. As well as presenting superior cryopreservation procedures for HUVECs, the methods developed here can serve as a model to optimize the cryopreservation of other cells.

Foundations of modeling in cryobiology-I: concentration, Gibbs energy, and chemical potential relationships

Mathematical modeling plays an enormously important role in understanding the behavior of cells, tissues, and organs undergoing cryopreservation. Uses of these models range from explanation of phenomena, exploration of potential theories of damage or success, development of equipment, and refinement of optimal cryopreservation/cryoablation strategies. Over the last half century there has been a considerable amount of work in bio-heat and mass-transport, and these models and theories have been readily and repeatedly applied to cryobiology with much success. However, there are significant gaps between experimental and theoretical results that suggest missing links in models. One source for these potential gaps is that cryobiology is at the intersection of several very challenging aspects of transport theory: it couples multi-component, moving boundary, multiphase solutions that interact through a semipermeable elastic membrane with multicomponent solutions in a second time-varying domain, during a two-hundred Kelvin temperature change with multi-molar concentration gradients and multi-atmosphere pressure changes. In order to better identify potential sources of error, and to point to future directions in modeling and experimental research, we present a three part series to build from first principles a theory of coupled heat and mass transport in cryobiological systems accounting for all of these effects. The hope of this series is that by presenting and justifying all steps, conclusions may be made about the importance of key assumptions, perhaps pointing to areas of future research or model development, but importantly, lending weight to standard simplification arguments that are often made in heat and mass transport. In this first part, we review concentration variable relationships, their impact on choices for Gibbs energy models, and their impact on chemical potentials.