Effects of Transmembrane α-Helix Length and Concentration on Phase Behavior in Four-Component Lipid Mixtures: A Molecular Dynamics Study

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.

A Trifunctional Linker for Palmitoylation and Peptide and Protein Localization in Biological Membranes

Attachment of lipophilic groups is an important post-translational modification of proteins, which involves the coupling of one or more anchors such as fatty acids, isoprenoids, phospholipids, or glycosylphosphatidyl inositols. To study its impact on the membrane partitioning of hydrophobic peptides or proteins, we designed a tyrosine-based trifunctional linker. The linker allows the facile incorporation of two different functionalities at a cysteine residue in a single step. We determined the effect of the lipid modification on the membrane partitioning of the synthetic α-helical model peptide WALP with or without here and in all cases below; palmitoyl groups in giant unilamellar vesicles that contain a liquid-ordered (Lo ) and liquid-disordered (Ld ) phase. Introduction of two palmitoyl groups did not alter the localization of the membrane peptides, nor did the membrane thickness or lipid composition. In all cases, the peptide was retained in the Ld phase. These data demonstrate that the Lo domain in model membranes is highly unfavorable for a single membrane-spanning peptide.

Inserting Small Molecules across Membrane Mixtures: Insight from the Potential of Mean Force

Small solutes have been shown to alter the lateral organization of cell membranes and reconstituted phospholipid bilayers; however, the mechanisms by which these changes happen are still largely unknown. Traditionally, both experiment and simulation studies have been restricted to testing only a few compounds at a time, failing to identify general molecular descriptors or chemical properties that would allow extrapolating beyond the subset of considered solutes. In this work, we probe the competing energetics of inserting a solute in different membrane environments by means of the potential of mean force. We show that these calculations can be used as a computationally efficient proxy to establish whether a solute will stabilize or destabilize domain phase separation. Combined with umbrella-sampling simulations and coarse-grained molecular dynamics simulations, we are able to screen solutes across a wide range of chemistries and polarities. Our results indicate that for the system under consideration, preferential partitioning and therefore effectiveness in altering membrane phase separation are strictly linked to the location of insertion in the bilayer (i.e., midplane or interface). Our approach represents a fast and simple tool for obtaining structural and thermodynamic insight into the partitioning of small molecules between lipid domains and its relation to phase separation, ultimately providing a platform for identifying the key determinants of this process.

Membrane Bending Moduli of Coexisting Liquid Phases Containing Transmembrane Peptide

A number of highly curved membranes in vivo, such as epithelial cell microvilli, have the relatively high sphingolipid content associated with "raft-like" composition. Given the much lower bending energy measured for bilayers with "nonraft" low sphingomyelin and low cholesterol content, observing high curvature for presumably more rigid compositions seems counterintuitive. To understand this behavior, we measured membrane rigidity by fluctuation analysis of giant unilamellar vesicles. We found that including a transmembrane helical GWALP peptide increases the membrane bending modulus of the liquid-disordered (Ld) phase. We observed this increase at both low-cholesterol fraction and higher, more physiological cholesterol fraction. We find that simplified, commonly used Ld and liquid-ordered (Lo) phases are not representative of those that coexist. When Ld and Lo phases coexist, GWALP peptide favors the Ld phase with a partition coefficient of 3-10 depending on mixture composition. In model membranes at high cholesterol fractions, Ld phases with GWALP have greater bending moduli than the Lo phase that would coexist.

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.

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.

Partitioning of nanoscale particles on a heterogeneous multicomponent lipid bilayer

Partitioning of nanoparticles into different lipid phases of a cell membrane is regulated by the physical properties of both the membrane and nanoparticles.

Line Tension Controls Liquid-Disordered + Liquid-Ordered Domain Size Transition in Lipid Bilayers

To better understand animal cell plasma membranes, we studied simplified models, namely four-component lipid bilayer mixtures. Here we describe the domain size transition in the region of coexisting liquid-disordered (Ld) + liquid-ordered (Lo) phases. This transition occurs abruptly in composition space with domains increasing in size by two orders of magnitude, from tens of nanometers to microns. We measured the line tension between coexisting Ld and Lo domains close to the domain size transition for a variety of lipid mixtures, finding that in every case the transition occurs at a line tension of ∼0.3 pN. A computational model incorporating line tension and dipole repulsion indicated that even small changes in line tension can result in domains growing in size by several orders of magnitude, consistent with experimental observations. We find that other properties of the coexisting Ld and Lo phases do not change significantly in the vicinity of the abrupt domain size transition.