Predict the glass transition temperature of glycerol-water binary cryoprotectant by molecular dynamic simulation

Vitrification is proposed to be the best way for the cryopreservation of organs. The glass transition temperature (T_g) of vitrification solutions is a critical parameter of fundamental importance for cryopreservation by vitrification. The instruments that can detect the thermodynamic, mechanical and dielectric changes of a substance may be used to determine the glass transition temperature. T_g is usually measured by using differential scanning calorimetry (DSC). In this study, the T_g of the glycerol-aqueous solution (60%, wt/%) was determined by isothermal-isobaric molecular dynamic simulation (NPT-MD). The software package Discover in Material Studio with the Polymer Consortium Force Field (PCFF) was used for the simulation. The state parameters of heat capacity at constant pressure (C_p), density (ρ), amorphous cell volume (V_cell) and specific volume (V_specific) and radial distribution function (rdf) were obtained by NPT-MD in the temperature range of 90–270 K. These parameters showed a discontinuity at a specific temperature in the plot of state parameter versus temperature. The temperature at the discontinuity is taken as the simulated T_g value for glycerol–water binary solution. The T_g values determined by simulation method were compared with the values in the literatures. The simulation values of T_g (160.06–167.51 K) agree well with the DSC results (163.60–167.10 K) and the DMA results (159.00 K). We drew the conclusion that molecular dynamic simulation (MDS) is a potential method for investigating the glass transition temperature (T_g) of glycerol–water binary cryoprotectants and may be used for other vitrification solutions.