Water replacement hypothesis in atomic details: effect of trehalose on the structure of single dehydrated POPC bilayers.

We present molecular dynamics (MD) simulations to study the plausibility of the water replacement hypothesis (WRH) from the viewpoint of structural chemistry. A total of 256 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) lipids were modeled for 400 ns at 11.7 or 5.4 waters/lipid. To obtain a single dehydrated bilayer relevant to the WRH, simulations were performed in the NP(xy)h(z)T ensemble with h(z) > 8 nm, allowing interactions between lipids in the membrane plane and preventing interactions between neighboring membranes via periodic boundary conditions. This setup resulted in a stable single bilayer in (or near) the gel state. Trehalose caused a concentration-dependent increase of the area per lipid (APL) accompanied by fluidizing the bilayer core. This mechanism has been suggested by the WRH. However, dehydrated bilayers in the presence of trehalose were not structurally identical to fully hydrated bilayers. The headgroup vector was in a more parallel orientation in dehydrated bilayers with respect to the bilayer plane and maintained this orientation in the presence of trehalose in spite of APL increase. The total dipole potential changed sign in dehydrated bilayers and remained slightly positive in the presence of trehalose. The model of a dehydrated bilayer presented here allows the study of the mechanisms of membrane protection against desiccation by different compounds.