Molecular dynamics simulation study of the effect of glycerol dialkyl glycerol tetraether hydroxylation on membrane thermostability

Abstract Archaeal tetraether membrane lipids span the whole membrane width and present two C 40 isoprenoid chains bound by two glycerol groups (or one glycerol and calditol). These lipids confer stability and maintain the membrane fluidity in mesophile to extremophile environments, making them very attractive for biotechnological applications. The isoprenoid lipid composition in archaeal membranes varies with temperature, which has placed these lipids in the focus of paleo-climatological studies for over a decade. Non-hydroxylated isoprenoid archaeal lipids are typically used as paleo-thermometry proxies, but recently identified hydroxylated (OH) derivatives have also been proposed as temperature proxies. The relative abundance of hydroxylated lipids increases at lower temperatures, but the physiological function of the OH moiety remains unknown. Here we present molecular dynamics simulations of membranes formed by the acyclic glycerol-dialkyl-glycerol-tetraether caldarchaeol (GDGT-0), the most widespread archaeal core lipid, and its mono-hydroxylated variant (OH-GDGT-0) to better understand the physico-chemical properties conferred to the membrane by this additional moiety. The molecular dynamics simulations indicate that the additional OH group forms hydrogen bonds mainly with the sugar moieties of neighbouring lipids and with water molecules, effectively increasing the size of the polar headgroups. The hydroxylation also introduces local disorder that propagates along the entire alkyl chains, resulting in a slightly more fluid membrane. These changes would help to maintain trans-membrane transport in cold environments, explaining why the relative abundance of hydroxylated Archaea lipids increases at lower temperatures. The in silico approach aids to understand the underlying physiological mechanisms behind the hydroxylated lipid based paleo-thermometer recently proposed.