Phase Separation in a Lipid/Cholesterol System: Comparison of Coarse-Grained and United-Atom Simulations

Photo-Responsive Control of Adsorption and Structure Formation at the Air-Water Interface with Arylazopyrazoles

Smart interfaces that are responsive to external triggers such as light are of great interest for the development of responsive or adaptive materials and interfaces. Using alkyl-arylazopyrazole butyl sulfonate surfactants (alkyl-AAP) that can undergo E/Z photoisomerization when irradiated with green (E) and UV (Z) lights, we demonstrate through a combination of experiments and computer simulations that there can be surprisingly large changes in surface tension and in the molecular structure and order at air-water interfaces. Surface tensiometry, vibrational sum-frequency generation (SFG) spectroscopy, and neutron reflectometry (NR) are applied to the study of custom-synthesized AAP surfactants with octyl- and H-terminal groups at air-water interfaces as a function of their bulk concentration and E/Z configuration. Upon photoswitching, a drastic influence of the alkyl chain on both the surface activity and the responsiveness of interfacial surfactants is revealed from changes in the surface tension, γ, where the largest changes in γ are observed for octyl-AAP (Δγ ∼ 23 mN/m) in contrast to H-AAP with Δγ < 10 mN/m. Results from vibrational SFG spectroscopy and NR show that the interfacial composition and the molecular order of the surfactants drastically change with E/Z photoisomerization and surface coverage. Indeed, from analysis of the S-O (head group) and C-H vibrational bands (hydrophobic tail), a qualitative analysis of orientational and structural changes of interfacial AAP surfactants is provided. The experiments are complemented by resolution of thermodynamic parameters such as equilibrium constants from ultra-coarse-grained simulations, which also capture details like island formation and interaction parameters of interfacial molecules. Here, the interparticle interaction ("stickiness") and the interaction with the surface are adjusted, closely reflecting experimental conditions.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Frame-guided assembly from a theoretical perspective

The molecular self-assembly of various structures such as micelles and vesicles has been the subject of comprehensive studies. Recently, a new approach to design these structures, the frame-guided assembly, has been developed to progress towards fabrics of predefined shape and size, following an initially provided frame of guiding elements. Here we study frame-guided assembly into a two-dimensional membrane via computer simulations, based on a single-bead coarse grained surfactant model in continuous space. In agreement with the experiment the assembly process already starts for surfactant concentrations below the critical micelle concentration. Furthermore, upon increasing temperature the formation process gets more delocalized. Additionally, the assembly process of the resulting membrane plane is modelled by a lattice gas model. It displays a similar phenomenology but additionally allows the derivation of analytical mean-field predictions. In this way a fundamental understanding of frame-guided assembly can be gained.

Capturing Membrane Phase Separation by Dual Resolution Molecular Dynamics Simulations

Understanding the lateral organization in plasma membranes remains an open problem and is of great interest to many researchers. Model membranes consisting of coexisting domains are commonly used as simplified models of plasma membranes. The coarse-grained (CG) Martini force field has successfully captured spontaneous separation of ternary membranes into a liquid-disordered and a liquid-ordered domain. With all-atom (AA) models, however, phase separation is much harder to achieve due to the slow underlying dynamics. To remedy this problem, here, we apply the virtual site (VS) hybrid method on a ternary membrane composed of saturated lipids, unsaturated lipids, and cholesterol to investigate the phase separation. The VS scheme couples the two membrane leaflets at CG and AA resolution. We found that the rapid phase separation reached by the CG leaflet can accelerate and guide this process in the AA leaflet.

Nonconverged Constraints Cause Artificial Temperature Gradients in Lipid Bilayer Simulations

Molecular dynamics (MD) simulations have become an indispensable tool to investigate phase separation in model membrane systems. In particular, simulations based on coarse-grained (CG) models have found widespread use due to their increased computational efficiency, allowing for simulations of multicomponent lipid bilayers undergoing phase separation into liquid-ordered and liquid-disordered domains. Here, we show that a significant temperature difference between molecule types can artificially arise in CG MD membrane simulations with the standard Martini simulation parameters in GROMACS. In particular, the linear constraint solver (LINCS) algorithm does not converge with its default settings, resulting in serious temperature differences between molecules in a time step-dependent manner. We demonstrate that the underlying reason for this behavior is the presence of highly constrained moieties, such as cholesterol. Their presence can critically impact numerous structural and dynamic membrane properties obtained from such simulations. Furthermore, any preference of these molecules toward a certain membrane phase can lead to spatial temperature gradients, which can amplify the degree of phase separation or even induce it in compositions that would otherwise mix well. We systematically investigated the effect of the integration time step and LINCS settings on membrane properties. Our data show that for cholesterol-containing membranes, a time step of 20 fs should be combined with at least lincs_iter = 2 and lincs_order = 12, while using a time step of 30 fs requires at least lincs_iter = 3 and lincs_order = 12 to bring the temperature differences to a level where they do not perturb central membrane properties. Moreover, we show that in cases where stricter LINCS settings are computationally too demanding, coupling the lipids in multiple groups to the temperature bath offers a practical workaround to the problem, although the validity of this approach should be further verified. Finally, we show that similar temperature gradients can also emerge in atomistic simulations using the CHARMM force field in combination with settings that allow for a 5 fs integration step.

Role of polyunsaturated phospholipids in liquid-ordered and liquid-disordered phases

Polyunsaturated phospholipids play a strong repulsive role in the liquid-disordered phase but a weak role in the liquid-ordered phase.

Bilayer compositional asymmetry influences the nanoscopic to macroscopic phase domain size transition

The eukaryotic plasma membrane (PM) exhibits lipid mixing heterogeneities known as lipid rafts. These lipid rafts, the result of liquid-liquid phase separation, can be modeled by coexisting liquid ordered (Lo) and liquid disordered (Ld) domains. Four-lipid component systems with a high-melting lipid, a nanodomain-inducing low-melting lipid, a macrodomain-inducing low-melting lipid, and cholesterol (chol) can give rise to domains of different sizes. These four-component systems have been characterized in experiments, yet there are few studies that model the asymmetric distribution of lipids actually found in the PM. We used molecular dynamics (MD) simulations to analyze the transition from nanoscopic to macroscopic domains in symmetric and in asymmetric model membranes. Using coarse-grained MD simulations, we found that asymmetry promotes macroscopic domain growth in a case where symmetric systems exhibit nanoscopic domains. Also, macroscopic domain formation in symmetric systems is highly dependent on registration of like phases in the cytoplasmic and exoplasmic leaflets. Using united-atom MD simulations, we found that symmetric Lo domains are only slightly more ordered than asymmetric Lo domains. We also found that large Lo domains in our asymmetric systems induce a slight chain ordering in the apposed cytoplasmic regions. The chol fractions of phase-separated Lo and Ld domains of the exoplasmic leaflet were unchanged whether the system was symmetric or asymmetric.

To Hop or not to Hop: Exceptions in the FCS Diffusion Law

Diffusion obstacles in membranes have not been directly visualized because of fast membrane dynamics and the occurrence of subresolution molecular complexes. To understand the obstacle characteristics, mobility-based methods are often used as an indirect way of assessing the membrane structure. Molecular movement in biological plasma membranes is often characterized by anomalous diffusion, but the exact underlying mechanisms are still elusive. Imaging total internal reflection fluorescence correlation spectroscopy (ITIR-FCS) is a well-established mobility-based method that provides spatially resolved diffusion coefficient maps and is combined with FCS diffusion law analysis to examine subresolution membrane organization. In recent years, although FCS diffusion law analysis has been instrumental in providing new insights into the membrane structure below the optical diffraction limit, there are certain exceptions and anomalies that require further clarification. To this end, we correlate the membrane structural features imaged by atomic force microscopy (AFM) with the dynamics measured using ITIR-FCS. We perform ITIR-FCS measurements on supported lipid bilayers (SLBs) of various lipid compositions to characterize the anomalous diffusion of lipid molecules in distinct obstacle configurations, along with the high-resolution imaging of the membrane structures with AFM. Furthermore, we validate our experimental results by performing simulations on image grids with experimentally determined obstacle configurations. This study demonstrates that FCS diffusion law analysis is a powerful tool to determine membrane heterogeneities implied from dynamics measurements. Our results corroborate the commonly accepted interpretations of imaging FCS diffusion law analysis, and we show that exceptions happen when domains reach the percolation threshold in a biphasic membrane and a network of domains behaves rather like a meshwork, resulting in hop diffusion.

Dual Resolution Membrane Simulations Using Virtual Sites

All-atomistic (AA) and coarse-grain (CG) simulations have been successfully applied to investigate a broad range of biomolecular processes. However, the accessible time and length scales of AA simulation are limited and the specific molecular details of CG simulation are simplified. Here, we propose a virtual site (VS) based hybrid scheme that can concurrently couple AA and CG resolutions in a single membrane simulation, mitigating the shortcomings of either representation. With some adjustments to make the AA and CG force fields compatible, we demonstrate that lipid bilayer properties are well kept in our hybrid approach. Our VS hybrid method was also applied to simulate a small lipid vesicle, with the inner leaflet and interior solvent represented in AA, and the outer leaflet together with exterior solvent at the CG level. Our multiscale method opens the way to investigate biomembrane properties at increased computational efficiency, in particular applications involving large solvent filled regions.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Aescin-Cholesterol Complexes in DMPC Model Membranes: A DSC and Temperature-Dependent Scattering Study

The saponin aescin, a mixture of triterpenoid saponins, is obtained from the seeds of the horse chestnut tree Aesculus hippocastanum. The β-form employed in this study is haemolytically active. The haemolytic activity results from the ability of aescin to form strong complexes with cholesterol in the red blood cell membrane. In this study, we provide a structural analysis on the complex formation of aescin and cholesterol when embedded in a phospholipid model membrane formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). In this work, the temperatures investigated extend from DMPC’s L_β′ to its L_α phase in dependence of different amounts of the saponin (0–6 mol% for calorimetric and 0–1 mol% for structural analyses) and the steroid (1–10 mol%). At these aescin contents model membranes are conserved in the form of small unilamellar vesicles (SUVs) and major overall structural modifications are avoided. Additionally, interactions between aescin and cholesterol can be studied for both phase states of the lipid, the gel and the fluid state. From calorimetric experiments by differential scanning calorimetry (DSC), it could be shown that both, the steroid and the saponin content, have a significant impact on the cooperative phase transition behaviour of the DMPC molecules. In addition, it becomes clearly visible that the entire phase behaviour is dominated by phase separation which indeed also depends on the complexes formed between aescin and cholesterol. We show by various methods that the addition of cholesterol alters the impact of aescin on structural parameters ranging from the acyl chain correlation to vesicle-vesicle interactions. While the specific saponin-phospholipid interaction is reduced, addition of cholesterol leads to deformation of SUVs. The analyses of the structures formed were performed by wide-angle X-ray scattering (WAXS), small-angle X-ray scattering (SAXS), and small-angle neutron scattering (SANS).

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Structure of lateral heterogeneities in a coarse-grained model for multicomponent membranes

We study the lateral domain structure in a coarse-grained molecular model for multicomponent lipid bilayers by semi-grandcanonical Monte Carlo simulations. The membranes are filled with liquid ordered (lo) domains surrounded by a liquid disordered (ld) matrix. Depending on the membrane composition and temperature, we identify different morphological regimes: one regime (I) where the lo domains are small and relatively compact, and two regimes (II, II') where they are larger and often interconnected. In the latter two regimes, the ld matrix forms a network of disordered trenches separating the lo domains, with a relatively high content of interdigitated line defects. Since such defects are also a structural element of the modulated ripple phase in one component membranes, we argue that the regimes II, II' may be amorphous equivalents of the ripple phase in multicomponent membranes. We also analyze the local structure and provide evidence that the domains in regime I are stabilized by a monolayer curvature mechanism postulated in earlier work [S. Meinhardt et al., PNAS, 2013, 110, 4476].

ForConX: A forcefield conversion tool based on XML

The force field conversion from one MD program to another one is exhausting and error-prone. Although single conversion tools from one MD program to another exist not every combination and both directions of conversion are available for the favorite MD programs Amber, Charmm, Dl-Poly, Gromacs, and Lammps. We present here a general tool for the force field conversion on the basis of an XML document. The force field is converted to and from this XML structure facilitating the implementation of new MD programs for the conversion. Furthermore, the XML structure is human readable and can be manipulated before continuing the conversion. We report, as testcases, the conversions of topologies for acetonitrile, dimethylformamide, and 1-ethyl-3-methylimidazolium trifluoromethanesulfonate comprising also Urey-Bradley and Ryckaert-Bellemans potentials. © 2017 Wiley Periodicals, Inc.

Regimes of Complex Lipid Bilayer Phases Induced by Cholesterol Concentration in MD Simulation

Cholesterol is essential to the formation of phase-separated lipid domains in membranes. Lipid domains can exist in different thermodynamic phases depending on the molecular composition and play significant roles in determining structure and function of membrane proteins. We investigate the role of cholesterol in the structure and dynamics of ternary lipid mixtures displaying phase separation using molecular dynamics simulations, employing a physiologically relevant span of cholesterol concentration. We find that cholesterol can induce formation of three regimes of phase behavior: 1) miscible liquid-disordered bulk, 2) phase-separated, domain-registered coexistence of liquid-disordered and liquid-ordered domains, and 3) phase-separated, domain-antiregistered coexistence of liquid-disordered and newly identified nanoscopic gel domains composed of cholesterol threads we name "cholesterolic gel" domains. These findings are validated and discussed in the context of current experimental knowledge, models of cholesterol spatial distributions, and models of ternary lipid-mixture phase separation.

Phospholipid Chain Interactions with Cholesterol Drive Domain Formation in Lipid Membranes

Cholesterol is a key component of eukaryotic membranes, but its role in cellular biology in general and in lipid rafts in particular remains controversial. Model membranes are used extensively to determine the phase behavior of ternary mixtures of cholesterol, a saturated lipid, and an unsaturated lipid with liquid-ordered and liquid-disordered phase coexistence. Despite many different experiments that determine lipid-phase diagrams, we lack an understanding of the molecular-level driving forces for liquid phase coexistence in bilayers with cholesterol. Here, we use atomistic molecular dynamics computer simulations to address the driving forces for phase coexistence in ternary lipid mixtures. Domain formation is directly observed in a long-timescale simulation of a mixture of 1,2-distearoyl-sn-glycero-3-phosphocholine, unsaturated 1,2-dilinoleoyl-sn-glycero-3-phosphocholine, and cholesterol. Free-energy calculations for the exchange of the saturated and unsaturated lipids between the ordered and disordered phases give insight into the mixing behavior. We show that a large energetic contribution to domain formation is favorable enthalpic interactions of the saturated lipid in the ordered phase. This favorable energy for forming an ordered, cholesterol-rich phase is opposed by a large unfavorable entropy. Martini coarse-grained simulations capture the unfavorable free energy of mixing but do not reproduce the entropic contribution because of the reduced representation of the phospholipid tails. Phospholipid tails and their degree of unsaturation are key energetic contributors to lipid phase separation.

Critical size dependence of domain formation observed in coarse-grained simulations of bilayers composed of ternary lipid mixtures

Model cellular membranes are known to form micro- and macroscale lipid domains dependent on molecular composition. The formation of macroscopic lipid domains by lipid mixtures has been the subject of many simulation investigations. We present a critical study of system size impact on lipid domain phase separation into liquid-ordered and liquid-disordered macroscale domains in ternary lipid mixtures. In the popular di-C16:0 PC:di-C18:2 PC:cholesterol at 35:35:30 ratio mixture, we find systems with a minimum of 1480 lipids to be necessary for the formation of macroscopic phase separated domains and systems of 10 000 lipids to achieve structurally converged conformations similar to the thermodynamic limit. To understand these results and predict the behavior of any mixture forming two phases, we develop and investigate an analytical Flory-Huggins model which is recursively validated using simulation and experimental data. We find that micro- and macroscale domains can coexist in ternary mixtures. Additionally, we analyze the distributions of specific lipid-lipid interactions in each phase, characterizing domain structures proposed based on past experimental studies. These findings offer guidance in selecting appropriate system sizes for the study of phase separations and provide new insights into the nature of domain structure for a popular ternary lipid mixture.

Exploring the free energy gain of phase separation via Markov state modeling

The gain of free energy upon unmixing is determined via application of Markov state modeling (MSM), using an Ising model with a fixed number of up- and down-spins. MSM yields reasonable estimates of the free energies. However, a closer look reveals significant differences that point to residual non-Markovian effects. These non-Markovian effects are rather unexpected since the typical criteria to study the quality of Markovianity indicate complete Markovian behavior. We identify the sparse connectivity between different Markov states as a likely reason for the observed bias. By studying a simple five state model system, we can analytically elucidate different sources of the bias and thus explain the different deviations that were observed for the Ising model. Based on this insight, we can modify the determination of the count matrix in the MSM approach. In this way, the estimation of the free energy is significantly improved.

Physical mechanisms of micro- and nanodomain formation in multicomponent lipid membranes

This article summarizes a variety of physical mechanisms proposed in the literature, which can generate micro- and nanodomains in multicomponent lipid bilayers and biomembranes. It mainly focusses on lipid-driven mechanisms that do not involve direct protein-protein interactions. Specifically, it considers (i) equilibrium mechanisms based on lipid-lipid phase separation such as critical cluster formation close to critical points, and multiple domain formation in curved geometries, (ii) equilibrium mechanisms that stabilize two-dimensional microemulsions, such as the effect of linactants and the effect of curvature-composition coupling in bilayers and monolayers, and (iii) non-equilibrium mechanisms induced by the interaction of a biomembrane with the cellular environment, such as membrane recycling and the pinning effects of the cytoplasm. Theoretical predictions are discussed together with simulations and experiments. The presentation is guided by the theory of phase transitions and critical phenomena, and the appendix summarizes the mathematical background in a concise way within the framework of the Ginzburg-Landau theory. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Coarse-grained simulations of polyelectrolyte complexes: MARTINI models for poly(styrene sulfonate) and poly(diallyldimethylammonium)

We present simulations of aqueous polyelectrolyte complexes with new MARTINI models for the charged polymers poly(styrene sulfonate) and poly(diallyldimethylammonium). Our coarse-grained polyelectrolyte models allow us to study large length and long time scales with regard to chemical details and thermodynamic properties. The results are compared to the outcomes of previous atomistic molecular dynamics simulations and verify that electrostatic properties are reproduced by our MARTINI coarse-grained approach with reasonable accuracy. Structural similarity between the atomistic and the coarse-grained results is indicated by a comparison between the pair radial distribution functions and the cumulative number of surrounding particles. Our coarse-grained models are able to quantitatively reproduce previous findings like the correct charge compensation mechanism and a reduced dielectric constant of water. These results can be interpreted as the underlying reason for the stability of polyelectrolyte multilayers and complexes and validate the robustness of the proposed models.

Lipid Raft Size and Lipid Mobility in Non-raft Domains Increase during Aging and Are Exacerbated in APP/PS1 Mice Model of Alzheimer's Disease. Predictions from an Agent-Based Mathematical Model

A connection between lipid rafts and Alzheimer's disease has been studied during the last decades. Mathematical modeling approaches have recently been used to correlate the effects of lipid composition changes in the physicochemical properties of raft-like membranes. Here we propose an agent based model to assess the effect of lipid changes in lipid rafts on the evolution and progression of Alzheimer's disease using lipid profile data obtained in an established model of familial Alzheimer's disease. We have observed that lipid raft size and lipid mobility in non-raft domains are two main factors that increase during age and are accelerated in the transgenic Alzheimer's disease mouse model. The consequences of these changes are discussed in the context of neurotoxic amyloid β production. Our agent based model predicts that increasing sterols (mainly cholesterol) and long-chain polyunsaturated fatty acids (LCPUFA) (mainly DHA, docosahexaenoic acid) proportions in the membrane composition might delay the onset and progression of the disease.

Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms

This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Membrane-mediated regulation of the intrinsically disordered CD3ϵ cytoplasmic tail of the TCR

The regulation of T-cell-mediated immune responses depends on the phosphorylation of immunoreceptor tyrosine-based activation motifs (ITAMs) on T-cell receptors. Although many details of the signaling cascades are well understood, the initial mechanism and regulation of ITAM phosphorylation remains unknown. We used molecular dynamics simulations to study the influence of different compositions of lipid bilayers on the membrane association of the CD3ϵ cytoplasmic tails of the T-cell receptors. Our results show that binding of CD3ϵ to membranes is modulated by both the presence of negatively charged lipids and the lipid order of the membrane. Free-energy calculations reveal that the protein-membrane interaction is favored by the presence of nearby basic residues and the ITAM tyrosines. Phosphorylation minimizes membrane association, rendering the ITAM motif more accessible to binding partners. In systems mimicking biological membranes, the CD3ϵ chain localization is modulated by different facilitator lipids (e.g., gangliosides or phosphoinositols), revealing a plausible regulatory effect on activation through the regulation of lipid composition in cell membranes.

Computer simulations suggest a key role of membranous nanodomains in biliary lipid secretion

The bile fluid contains various lipids that are secreted at the canalicular membrane of hepatocytes. As the secretion mechanism is still a matter of debate and a direct experimental observation of the secretion process is not possible so far, we used a mathematical model to simulate the extraction of the major bile lipids cholesterol, phosphatidylcholine and sphingomyelin from the outer leaflet of the canalicular membrane. Lipid diffusion was modeled as random movement on a triangular lattice governed by next-neighbor interaction energies. Phase separation in liquid-ordered and liquid-disordered domains was modeled by assigning two alternative ordering states to each lipid species and minimization of next-neighbor ordering energies. Parameterization of the model was performed such that experimentally determined diffusion rates and phases in ternary lipid mixtures of model membranes were correctly recapitulated. The model describes the spontaneous formation of nanodomains in the external leaflet of the canalicular membrane in a time window between 0.1 ms to 10 ms at varying lipid proportions. The extraction of lipid patches from the bile salt soluble nanodomain into the bile reproduced observed biliary phospholipid compositions for a physiologi-cal membrane composition. Comparing the outcome of model simulations with available experi-mental observations clearly favors the extraction of tiny membrane patches composed of about 100–400 lipids as the likely mechanism of biliary lipid secretion.

Formation of the bile is one of the central functions of the liver. The bile fluid aids in the digestion of edible fats and removal of drugs and toxins from the body. The bile fluid is mainly composed of bile salts (BS), phosphatidylcholine (PC) and cholesterol (CH) in a fairly fixed proportion that prevents liver impairment by gallstone formation or cholestasis. During bile formation, BS are actively pumped out of the hepatocyte into the extracellular space where they extract PC and CH from the canalicular membrane. This extraction process bears the risk for the canalicular membrane to be destructed. Hence, only a certain fraction of the membrane should be accessible to the solubilizing activity of BS. We have developed a mathematical model that describes the temporal formation of CH-enriched ordered and PC-enriched disordered nanodomains in the canalicular membrane. Model simulations reveal that the disordered nanodomains exhibit a composition of PC and CH similar to that also found in the bile. From this finding and the good concordance of model simulations with experimental data we conclude that PC and CH are mainly secreted into the bile from the disordered nanodomain. Our work adds a new layer of physiological importance to the spontaneous formation of lipid domains in biological membranes.

Multiscale modeling of four-component lipid mixtures: domain composition, size, alignment, and properties of the phase interface

Simplified lipid mixtures are often used to model the complex behavior of the cell plasma membrane. Indeed, as few as four components-a high-melting lipid, a nandomain-inducing low-melting lipid, a macrodomain-inducing low-melting lipid, and cholesterol (chol)-can give rise to a wide range of domain sizes and patterns that are highly sensitive to lipid compositions. Although these systems are studied extensively with experiments, the molecular-level details governing their phase behavior are not yet known. We address this issue by using molecular dynamics simulations to analyze how phase separation evolves in a four-component system as it transitions from small domains to large domains. To do so, we fix concentrations of the high-melting lipid 16:0,16:0-phosphatidylcholine (DPPC) and chol, and incrementally replace the nanodomain-inducing low-melting lipid 16:0,18:2-PC (PUPC) by the macrodomain-inducing low-melting lipid 18:2,18:2-PC (DUPC). Coarse-grained simulations of this four-component system reveal that lipid demixing increases as the amount of DUPC increases. Additionally, we find that domain size and interleaflet alignment change sharply over a narrow range of replacement of PUPC by DUPC, indicating that intraleaflet and interleaflet behaviors are coupled. Corresponding united atom simulations show that only lipids within ∼2 nm of the phase interface are significantly perturbed regardless of domain composition or size. Thus, whereas the fraction of interface-perturbed lipids is negligible for large domains, it is significant for smaller ones. Together, these results reveal characteristic traits of bilayer thermodynamic behavior in four-component mixtures, and provide a baseline for investigation of the effects of proteins and other lipids on membrane phase properties.

Computer simulations of phase separation in lipid bilayers and monolayers

Studying phase coexistence in lipid bilayers and monolayers is important for understanding lipid-lipid interactions underlying lateral organization in biological membranes. Computer simulations follow experimental approaches and use model lipid mixtures of simplified composition. Atomistic simulations give detailed information on the specificity of intermolecular interactions, while coarse-grained simulations achieve large time and length scales and provide a bridge towards state-of-the-art experimental techniques. Computer simulations allow characterizing the structure and composition of domains during phase transformations at Angstrom and picosecond resolution, and bring new insights into phase behavior of lipid membranes.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Lipid clustering correlates with membrane curvature as revealed by molecular simulations of complex lipid bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.

Cell membranes play important roles in vivo both in shielding the cell interior from the surrounding environment and in cell function through lipid components of the membrane having roles in controlling protein function, cell signaling etc. We employ molecular dynamics simulations to explore the behavior of biologically realistic membrane models. Our simulations reveal nano-domain clustering of the glycolipid GM3 and to a lesser extent of the anionic lipid phosphatidylinositol 4,5-bisphophate (PIP₂). When including transmembrane proteins we are able to observe preferential interactions of known regulatory lipids (e.g. GM3, PIP₂ and cholesterol) with the proteins. Membrane curvature is shown to be coupled to the local lipid composition, suggestive of a link between lipid nano-domains and membrane geometry.

Methodologies for the analysis of instantaneous lipid diffusion in MD simulations of large membrane systems

Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically comprise millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library.

Cholesterol affects C₆₀ translocation across lipid bilayers

Cholesterol plays an important role in regulating the structural properties of phospholipid membranes and further influences the permeability of molecules and nanoparticles. However, nanoparticles' translocation across phospholipid membranes in the presence of cholesterol on the molecular scale is rarely studied. Here, we performed coarse-grained molecular dynamics simulations to probe the translocation of C60, one of the most popular nanoparticles, across dipalmitoylphosphatidylcholine bilayers with different concentrations of cholesterol molecules (0-50 mol%). The results reveal that the presence of cholesterol molecules induces lower area per lipid, larger bilayer thickness, and more ordered orientation of lipid tails. The higher the concentration of cholesterol molecules, the more significant is the condensing effect of lipid bilayer as just mentioned. Besides, dynamic processes, free energy profiles and permeability coefficients further indicate that the permeability of C60 decreases with increasing cholesterol concentration, which can be explained by the condensation effect and reduced free volume. Our researches provide an explicit description of the impact of cholesterol on C60 translocation across lipid bilayers.

Key molecular requirements for raft formation in lipid/cholesterol membranes

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.

Pairing of cholesterol with oxidized phospholipid species in lipid bilayers

We claim that (1) cholesterol protects bilayers from disruption caused by lipid oxidation by sequestering conical shaped oxidized lipid species such as 1-palmitoyl-2-azelaoyl-sn-glycero-3-phosphocholine (PZPC) away from phospholipid, because cholesterol and the oxidized lipid have complementary shapes and (2) mixtures of cholesterol and oxidized lipids can self-assemble into bilayers much like lysolipid–cholesterol mixtures. The evidence for bilayer protection comes from molecular dynamics (MD) simulations and dynamic light scattering (DLS) measurements. Unimodal size distributions of extruded vesicles (LUVETs) made up of a mixture of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and PZPC containing high amounts of PZPC are only obtained when cholesterol is present in high concentrations. In simulations, bilayers containing high amounts of PZPC become porous, unless cholesterol is also present. The protective effect of cholesterol on oxidized lipids has been observed previously using electron paramagnetic resonance (EPR) and electron microscopy imaging of vesicles. The evidence for the pairing of cholesterol and PZPC comes mainly from correlated 2-D density and thickness plots from simulations, which show that these two molecules co-localize in bilayers. Further evidence that the two molecules can cohabitate comes from self-assembly simulations, where we show that cholesterol-oxidized lipid mixtures can form lamellar phases at specific concentrations, reminiscent of lysolipid–cholesterol mixtures. The additivity of the packing parameters of cholesterol and PZPC explains their cohabitation in a planar bilayer. Oxidized lipids are ubiquitously present in significant amounts in high- and low-density lipoprotein (HDL and LDL) particles, diseased tissues, and in model phospholipid mixtures containing polyunsaturated lipids. Therefore, our hypothesis has important consequences for cellular cholesterol trafficking; diseases related to oxidized lipids, and to biophysical studies of phase behaviour of cholesterol-containing phospholipid mixtures.