Trehalose-glycerol mixtures are known to be effective in the long time preservation of proteins. However, the microscopic mechanism of their effective preservation abilities remains unclear. In this article we present a molecular dynamics simulation study of the short time, less than 1 ns, dynamics of four trehalose-glycerol mixtures at temperatures below the glass transition temperature. We found that a mixture of 5% glycerol and 95% trehalose has the most suppressed short time dynamics (fast dynamics). This result agrees with the experimental analysis of the mean-square displacement of the hydrogen atoms, as measured via neutron scattering, and correlates with the experimentally observed enhancement of the stability of some enzymes at this particular concentration. Our microscopic analysis suggests that the formation of a robust intermolecular hydrogen bonding network is most effective at this concentration and is the main mechanism for the suppression of the fast dynamics.

We explore possible molecular mechanisms behind the coupling of protein and solvent dynamics using atomistic molecular-dynamics simulations. For this purpose, we analyze the model protein lysozyme in glycerol, a well-known protein-preserving agent. We find that the dynamics of the hydrogen bond network between the solvent molecules in the first shell and the surface residues of the protein controls the structural relaxation (dynamics) of the whole protein. Specifically, we find a power-law relationship between the relaxation time of the aforementioned hydrogen bond network and the structural relaxation time of the protein obtained from the incoherent intermediate scattering function. We demonstrate that the relationship between the dynamics of the hydrogen bonds and the dynamics of the protein appears also in the dynamic transition temperature of the protein. A study of the dynamics of glycerol as a function of the distance from the surface of the protein indicates that the viscosity seen by the protein is not the one of the bulk solvent. The presence of the protein suppresses the dynamics of the surrounding solvent. This implies that the protein sees an effective viscosity higher than the one of the bulk solvent. We also found significant differences in the dynamics of surface and core residues of the protein. The former is found to follow the dynamics of the solvent more closely than the latter. These results allowed us to propose a molecular mechanism for the coupling of the solvent-protein dynamics.

The traditional approach used to predict the ability of a glassy matrix to maximally preserve the activity of a protein solute is the glass transition temperature (T(g)) of the glass. Recently it has been shown that the addition of a low T(g) diluent (glycerol) can rigidify the structure of a high T(g) glassy matrix in binary glycerol-trehalose glasses. The optimal density of glycerol in trehalose minimizes the average mean square displacements of non-exchangeable protons in the glass samples. The amount of glycerol added to a trehalose glass coincides with the maximal recovery of biological activity in a separate study using similar binary glass samples. In this study, we use molecular dynamics (MD) simulations to investigate the dynamics of a hydrated protein encased in glycerol, unary trehalose and binary glycerol-trehalose glasses. We have found that we are able to reproduce the rigidification of the glycerol-trehalose glassy matrix and that there is a direct correlation between bulk glass dynamics and the extent of atomic fluctuation of protein atoms. The detailed microscopic picture that emerges is that protein dynamics are suppressed mainly by inertia of the bulk glass and to a lesser extent specific interactions at the protein-solvent interface. Thus, the inertia of the glassy matrix may be an influential factor in the determination of pharmaceutically relevant formulations.