Characterization of the Glass Transition of Water Predicted by Molecular Dynamics Simulations Using Nonpolarizable Intermolecular Potentials

Compositional relaxation on the approach to the glass transition in a model trehalose solution

Molecular dynamics simulation was used to study the temperature dependence of the mutual diffusion coefficient D_{m} and the intermediate scattering function of equilibrium and metastable aqueous solutions of the cryoprotectant molecule trehalose at very low (2.2 and 9wt.%) and very high (80 and 95wt.%) concentrations. The simulations were conducted over a range of temperatures approaching the glass transition temperature T_{g} for each concentration. Similar to a recent observation made on a glass-forming model polydisperse colloidal suspension [Hannam et al., Phys. Rev. E 96, 022609 (2017)2470-004510.1103/PhysRevE.96.022609], we confirmed by a set of independent computations that D_{m} is responsible for the long-time decay of the intermediate scattering function. We observed that D_{m} decreased on the approach to the glass transition temperature, resulting in an extremely slow long-time decay in the intermediate scattering function that culminated in the arrest of compositional fluctuations and a plateau in the intermediate scattering function at T_{g}. In both cases, crystallization requires a change in the composition of the solution, which is a process controlled by D_{m}. This transport coefficient can either increase or decrease as solidification is approached, because it depends on a product of thermodynamic and mobility factors. Our observations show that in both cases, for the glass-forming liquids, it is observed to decrease, while for a previously studied monodisperse colloidal suspension which crystallizes easily, it increases. The similarity in the behavior of these two very different glass-forming systems (the polydisperse colloidal suspension and the sugar solution) shows the importance of the mutual diffusion coefficient to our understanding of vitrification and suggests a possible distinction between between glass-forming and crystallizing solutions.

Molecular mechanism of the synergistic effects of vitrification solutions on the stability of phospholipid bilayers

The vitrification solutions used in the cryopreservation of biological samples aim to minimize the deleterious formation of ice by dehydrating cells and promoting the formation of the glassy state of water. They contain a mixture of different cryoprotective agents (CPAs) in water, typically polyhydroxylated alcohols and/or dimethyl sulfoxide (DMSO), which can damage cell membranes. Molecular dynamics simulations have been used to investigate the behavior of pure DPPC, pure DOPC, and mixed DOPC-β-sitosterol bilayers solvated in a vitrification solution containing glycerol, ethylene glycol, and DMSO at concentrations that approximate the widely used plant vitrification solution 2. As in the case of solutions containing a single CPA, the vitrification solution causes the bilayer to thin and become disordered, and pores form in the case of some bilayers. Importantly, the degree of thinning is, however, substantially reduced compared to solutions of DMSO containing the same total CPA concentration. The reduction in the damage done to the bilayers is a result of the ability of the polyhydroxylated species (especially glycerol) to form hydrogen bonds to the lipid and sterol molecules of the bilayer. A decrease in the amount of DMSO in the vitrification solution with a corresponding increase in the amount of glycerol or ethylene glycol diminishes further its damaging effect due to increased hydrogen bonding of the polyol species to the bilayer headgroups. These findings rationalize, to our knowledge for the first time, the synergistic effects of combining different CPAs, and form the basis for the optimization of vitrification solutions.