Line Tension Assists Membrane Permeation at the Transition Temperature in Mixed-Phase Lipid Bilayers

The umbrella sampling method has been used to evaluate the free energy profile for a large permeant moving through a lipid bilayer, represented using a coarse-grained simulation model, at and below its gel-fluid transition temperature. At the lipid transition temperature, determined to be 302 K for the MARTINI 2.0 model of DPPC, the permeation barrier for passage through an enclosed fluid domain embedded in a patch of gel was significantly lower than that for passage through a fluid stripe domain. In contrast, permeation through a fluid domain in a stripe geometry produced a free energy profile nearly identical to that of a gel-free fluid bilayer. In both cases, insertion of the permeant into a fluid domain coexisting with the gel phase led to a shift in phase composition, as lipids transitioned from fluid to gel to accommodate the area occupied by the permeant. In the case of the enclosed fluid domain, this transition produced a decrease in the length of the fluid-gel interface as the approximately circular fluid domain shrank. The observed decrease in the apparent permeation barrier, combined with an approximation for the change in interfacial length, enabled estimation of the interfacial line tension to be between 10 and 13 pN for this model. The permeation barrier was shown to drop even further in simulations performed at temperatures below the transition temperature. The results suggest a mechanism to explain the experimentally observed anomalous peak in the temperature-dependent permeability of lipid bilayers near their transition temperatures. The contribution of this mechanism toward the permeability of a gel phase containing a thermal distribution of fluid-phase domains is estimated using a simple statistical thermodynamic model.