Quantifying Acyl Chain Interdigitation in Simulated Bilayers via Direct Transbilayer Interactions

In a lipid bilayer, the interactions between the lipid hydrocarbon chains from opposing leaflets can influence membrane properties. These interactions include the phenomenon of interdigitation, in which an acyl chain of one leaflet extends past the bilayer midplane and into the opposing leaflet. While static interdigitation is well understood in gel-phase bilayers from X-ray diffraction measurements, much less is known about dynamic interdigitation in fluid phases. In this regard, atomistic molecular dynamics simulations can provide mechanistic information on interleaflet interactions that can be used to generate experimentally testable hypotheses. To address limitations of existing computational methodologies that provide results that are either indirect or averaged over time and space, here we introduce three novel ways of quantifying the extent of chain interdigitation. Our protocols include the analysis of instantaneous interactions at the level of individual carbon atoms, thus providing temporal and spatial resolution for a more nuanced picture of dynamic interdigitation. We compare the methods on bilayers composed of lipids with an equal total number of carbon atoms, but different mismatches between the <i>sn</i>-1 and <i>sn</i>-2 chain lengths. We find that these metrics, which are based on freely available software packages and are easy to implement, provide complementary details that help characterize various features of lipid-lipid contacts at the bilayer midplane. The new frameworks thus allow for a deeper look at fundamental molecular mechanisms underlying bilayer structure and dynamics and present a valuable expansion of the membrane biophysics toolkit.