Ternary mixtures of saturated and unsaturated phospholipids and cholesterol constitute a well-known model system to study raft formation in membranes. This phenomenon is, e.g., observed in cell membranes. Here, coarse-grained (CG) and microscopic united-atom (UA) simulations are performed to investigate the phase separation of a membrane bilayer for the ternary system of saturated 16:0 (DPPC) and unsaturated 18:2 (DUPC) phospholipids and cholesterol (CHOL). The results of a 9 μs UA simulation demonstrate the onset of phase separation and can be compared with properties of the corresponding CG system. While sharing the structural features of phase separation in the CG model, the onset of demixing for the UA model is 40 times slower. This factor can be rationalized by the different short-time diffusion constants. Various system properties such as order parameters and lipid-CHOL and CHOL-CHOL interactions are analyzed and compared between the UA and CG models. From the structural perspective, the phase separation process appears to be rather similar for both models, which illustrates once more the power of using CG approaches.

The lipid mixture of DPPC (saturated lipid)/DUPC (unsaturated lipid)/CHOL (cholesterol) is studied with respect to its ability to form liquid-ordered and liquid-disordered phases. We employ coarse-grained simulations with MARTINI force field. All three components are systematically modified in order to explore the relevant molecular properties, leading to phase separation. Specifically, we show that the DPPC/DUPC/CHOL system unmixes due to enthalpic DPPC-DPPC and DPPC-CHOL interactions. The phase separation remains unchanged, except for the formation of a gel phase at long times after decreasing the conformational degrees of freedom of the unsaturated DUPC. In contrast, the phase separation can be suppressed by softening the DPPC chains. In an attempt to mimic the ordering and unmixing effect of CHOL the latter is replaced by a stiff and shortened DPPC-like lipid. One still observes phase separation, suggesting that it is mainly the rigid and planar structure of CHOL which is important for raft formation. Addition of an extra bead to the head of CHOL has no notable impact on the phase separation of the system, supporting the irrelevance of the Umbrella model for the phase separation. Reduction of the conformational entropy of CHOL by stiffening its last bead results in a significant increase of the order of the DPPC/CHOL domain. This suggests that the conformational entropy of CHOL is important to prohibit the gelation process. The interleaflet interactions as mediated by the terminal molecular groups seem to have a strong impact on the possibility of a subsequent gelation process after phase separation.

The annexins are a family of Ca2+-regulated phospholipid binding proteins that are involved in membrane domain organization and membrane trafficking. Although they are widely studied and crystal structures are available for several soluble annexins their mode of membrane association has never been studied at the molecular level. Here we obtained molecular information on the annexin-membrane interaction that could serve as paradigm for the peripheral membrane association of cytosolic proteins by Molecular Dynamics simulations. We analyzed systems containing the monomeric annexin A2 (AnxA2), a membrane with negatively charged phosphatidylserine (POPS) lipids as well as Ca2+ ions. On the atomic level we identify the AnxA2 orientations and the respective residues which display the strongest interaction with Ca2+ ions and the membrane. The simulation results fully agree with earlier experimental findings concerning the positioning of bound Ca2+ ions. Furthermore, we identify for the first time a significant interaction between lysine residues of the protein and POPS lipids that occurs independently of Ca2+ suggesting that AnxA2-membrane interactions can also occur in a low Ca2+ environment. Finally, by varying Ca2+ concentrations and lipid composition in our simulations we observe a calcium-induced negative curvature of the membrane as well as an AnxA2-induced lipid ordering.

Based on all-atom Molecular Dynamics (MD) simulations of a lipid bilayer we present a systematic mapping on a 2D lattice model. Keeping the lipid type and the chain order parameter as key variables we derive a free energy functional, containing the enthalpic interaction of adjacent lipids as well as the tail entropy. The functional form of both functions is explicitly determined for saturated and polyunsaturated lipids. By studying the lattice model via Monte Carlo simulations it is possible to reproduce the temperature dependence of the distribution of order parameters of the pure lipids, including the prediction of the gel transition. Furthermore, application to a mixture of saturated and polyunsaturated lipids yields the correct phase separation behavior at lower temperatures with a simulation time reduced by approximately 7 orders of magnitude as compared to the corresponding MD simulations. Even the time-dependence of the de-mixing is reproduced on a semi-quantitative level. Due to the generality of the approach we envisage a large number of further applications, ranging from modeling larger sets of lipids, sterols, and solvent proteins to predicting nucleation barriers for the melting of lipids. Particularly, from the properties of the 2D lattice model one can directly read off the enthalpy and entropy change of the 1,2-dipalmitoyl-sn-glycero-3-phosphocholine gel-to-liquid transition in excellent agreement with experimental and MD results.

Graphical Abstract