Tristearin bilayers: structure of the aqueous interface and stability in the presence of surfactants

We report results of atomistic molecular dynamics simulations of an industrially-relevant, exemplar triacylglycerol (TAG), namely tristearin (TS), under aqueous conditions, at different temperatures and in the presence of an anionic surfactant, sodium dodecylbenzene sulphonate (SDBS). We predict the TS bilayers to be stable and in a gel phase at temperatures of 350 K and below. At 370 K the lipid bilayer was able to melt, but does not feature a stable liquid-crystalline phase bilayer at this elevated temperature. We also predict the structural characteristics of TS bilayers in the presence of SDBS molecules under aqueous conditions, where surfactant molecules are found to spontaneously insert into the TS bilayers. We model TS bilayers containing different amounts of SDBS, with the presence of SDBS imparting only a moderate effect on the structure of the system. Our study represents the first step in applying atomistic molecular dynamics simulations to the investigation of TAG-aqueous interfaces. Our results suggest that the CHARMM36 force-field appears suitable for the simulation of such systems, although the phase behaviour of the system may be shifted to lower temperatures than is the case for the actual system. Our findings provide a foundation for further simulation studies of the TS-aqueous interface.