Evaluating coarse-grained MARTINI force-fields for capturing the ripple phase of lipid membranes

Abstract

Phospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (P{beta}') membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained (CG) resolution. Here, with an aggregate ~2.5 s all-atom and ~122 s CG MD simulations, we systematically assess the ability of six CG MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid and water force-field (FF) variants, parametrized to capture the DPPC gel and fluid phases, for their ability to capture the P{beta}' phase, and compared observations with those from an all-atom FF. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (> 2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated P{beta}' and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the P{beta}' phase, similar to previous studies, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically-trapped structures caused by finite size effects rather than being representative of true P{beta}' phases. We showed that even a MARTINI FF variant that could capture the tilted L{beta}' gel phase, a prerequisite for stabilizing the P{beta}' phase, could not capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs (including MARTINI3) may require specific re-parametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=114 SRC="FIGDIR/small/408674v2_ufig1.gif" ALT="Figure 1"> View larger version (66K): org.highwire.dtl.DTLVardef@36a8fcorg.highwire.dtl.DTLVardef@90511eorg.highwire.dtl.DTLVardef@bdb3f1org.highwire.dtl.DTLVardef@70bb9a_HPS_FORMAT_FIGEXP M_FIG C_FIG