Hydrophobic mismatch-induced clustering as a primer for protein sorting in the secretory pathway

Super-resolution microscopy to study membrane nanodomains and transport mechanisms in the plasma membrane

Biological membranes are complex, heterogeneous, and dynamic systems that play roles in the compartmentalization and protection of cells from the environment. It is still a challenge to elucidate kinetics and real-time transport routes for molecules through biological membranes in live cells. Currently, by developing and employing super-resolution microscopy; increasing evidence indicates channels and transporter nano-organization and dynamics within membranes play an important role in these regulatory mechanisms. Here we review recent advances and discuss the major advantages and disadvantages of using super-resolution microscopy to investigate protein organization and transport within plasma membranes.

Energetics of the Transmembrane Peptide Sorting by Hydrophobic Mismatch

Hydrophobic mismatch between a lipid membrane and embedded transmembrane peptides or proteins plays a role in their lateral localization and function. Earlier studies have resolved numerous mechanisms through which the peptides and membrane proteins adapt to mismatch, yet the energetics of lateral sorting due to hydrophobic mismatch have remained elusive due to the lack of suitable computational or experimental protocols. Here, we pioneer a molecular dynamics simulation approach to study the sorting of peptides along a membrane thickness gradient. Peptides of different lengths tilt and diffuse along the membrane to eliminate mismatch with a rate directly proportional to the magnitude of mismatch. We extract the 2-dimensional free energy profiles as a function of local thickness and peptide orientation, revealing the relative contributions of sorting and tilting, and suggesting their thermally accessible regimes. Our approach can readily be applied to study other membrane systems of biological interest where hydrophobic mismatch, or membrane thickness in general, plays a role.

Impacts of Hydrophobic Mismatch on Antimicrobial Peptide Efficacy and Bilayer Permeabilization

Antimicrobial resistance continues to be a major threat to world health, with the continued emergence of resistant bacterial strains. Antimicrobial peptides have emerged as an attractive option for the development of novel antimicrobial compounds in part due to their ubiquity in nature and the general lack of resistance development to this class of molecules. In this work, we analyzed the antimicrobial peptide C18G and several truncated forms for efficacy and the underlying mechanistic effects of the sequence truncation. The peptides were screened for antimicrobial efficacy against several standard laboratory strains, and further analyzed using fluorescence spectroscopy to evaluate binding to model lipid membranes and bilayer disruption. The results show a clear correlation between the length of the peptide and the antimicrobial efficacy. Furthermore, there is a correlation between peptide length and the hydrophobic thickness of the bilayer, indicating that hydrophobic mismatch is likely a contributing factor to the loss of efficacy in shorter peptides.

Spontaneous local membrane curvature induced by transmembrane proteins

The (local) curvature of cellular membranes acts as a driving force for the targeting of membrane-associated proteins to specific membrane domains, as well as a sorting mechanism for transmembrane proteins, e.g., by accumulation in regions of matching spontaneous curvature. The latter measure was previously experimentally employed to study the curvature induced by the potassium channel KvAP and by aquaporin AQP0. However, the direction of the reported spontaneous curvature levels as well as the molecular driving forces governing the membrane curvature induced by these integral transmembrane proteins could not be addressed experimentally. Here, using both coarse-grained and atomistic molecular dynamics (MD) simulations, we report induced spontaneous curvature values for the homologous potassium channel Kv 1.2/2.1 Chimera (KvChim) and AQP0 embedded in unrestrained lipid bicelles that are in very good agreement with experiment. Importantly, the direction of curvature could be directly assessed from our simulations: KvChim induces a strong positive membrane curvature (≈0.036 nm^-1) whereas AQP0 causes a comparably small negative curvature (≈-0.019 nm^-1). Analyses of protein-lipid interactions within the bicelle revealed that the potassium channel shapes the surrounding membrane via structural determinants. Differences in shape of the protein-lipid interface of the voltage-gating domains between the extracellular and cytosolic membrane leaflets induce membrane stress and thereby promote a protein-proximal membrane curvature. In contrast, the water pore AQP0 displayed a high structural stability and an only faint effect on the surrounding membrane environment that is connected to its wedge-like shape.

Elastic deformations mediate interaction of the raft boundary with membrane inclusions leading to their effective lateral sorting

Liquid-ordered lipid domains represent a lateral inhomogeneity in cellular membranes. These domains have elastic and physicochemical properties different from those of the surrounding membrane. In particular, their thickness exceeds that of the disordered membrane. Thus, elastic deformations arise at the domain boundary in order to compensate for the thickness mismatch. In equilibrium, the deformations lead to an incomplete register of monolayer ordered domains: the elastic energy is minimal if domains in opposing monolayers lie on the top of each other, and their boundaries are laterally shifted by about 3 nm. This configuration introduces a region, composed of one ordered and one disordered monolayers, with an intermediate bilayer thickness. Besides, a jump in a local monolayer curvature takes place in this intermediate region, concentrating here most of the elastic stress. This region can participate in a lateral sorting of membrane inclusions by offering them an optimal bilayer thickness and local curvature conditions. In the present study, we consider the sorting of deformable lipid inclusions, undeformable peripheral and deeply incorporated peptide inclusions, and undeformable transmembrane inclusions of different molecular geometry. With rare exceptions, all types of inclusions have an affinity to the ordered domain boundary as compared to the bulk phases. The optimal lateral distribution of inclusions allows relaxing the elastic stress at the boundary of domains.

A distinct function of the retinoblastoma protein in the control of lipid composition identified by lipidomic profiling

Here, by combining lipidomics with transcriptome analysis, we demonstrate that Rb depletion in mouse embryonic fibroblastss induces significant alterations in their lipid composition. We discovered that Rb depletion induced increase in lysophosphatidylserine, diacylglycerol (DAG), fatty acid (FA), acylcarnitine, phosphatidylcholine (PC), arachidonoyl ethanolamine, and decrease in phosphatidylglycerol, monoacylglycerol, without change in total lipid per protein levels. Analysis of the acyl chain composition of DAG, PC and phosphatidylserine revealed increase of saturated and mono-unsaturated acyl chains with specific carbon chain length. Consistently, we observed that Rb depletion increased the levels of fatty acids with the corresponding carbon chain length and number of carbon-carbon double bondssuch as myristic acid (14:0), palmitic acid (16:0), stearic acid (18:0) and all forms of FA 18:1. Microarray analysis revealed that Rb depletion induced significant upregulation of enzymes involved in elongation and desaturation of fatty acids. Among these, we found that elongation of long chain fatty acid family member 6 (Elovl6) and stearoyl-CoA desaturase 1 (Scd1) are the most robustly controlled by Rb possibly through E2F and sterol regulatory element-binding protein transcription factors. Depletion of Elovl6 or Scd1 significantly suppressed colony formation, sphere formation and xenograft tumor growth of Rb-deficient tumor cells. Suppression of self-renewal by the SCD1 inhibitor was rescued upon supplementation of the mono-unsaturated fatty acids generated by this enzyme. This study suggests a novel role for Rb in suppressing the malignant progression of tumors by controlling the lipid composition.

Localisation Microscopy of Breast Epithelial ErbB-2 Receptors and Gap Junctions: Trafficking after γ-Irradiation, Neuregulin-1β, and Trastuzumab Application

In cancer, vulnerable breast epithelium malignance tendency correlates with number and activation of ErbB receptor tyrosine kinases. In the presented work, we observe ErbB receptors activated by irradiation-induced DNA injury or neuregulin-1β application, or alternatively, attenuated by a therapeutic antibody using high resolution fluorescence localization microscopy. The gap junction turnover coinciding with ErbB receptor activation and co-transport is simultaneously recorded. DNA injury caused by 4 Gray of 6 MeV photon γ-irradiation or alternatively neuregulin-1β application mobilized ErbB receptors in a nucleograde fashion—a process attenuated by trastuzumab antibody application. This was accompanied by increased receptor density, indicating packing into transport units. Factors mobilizing ErbB receptors also mobilized plasma membrane resident gap junction channels. The time course of ErbB receptor activation and gap junction mobilization recapitulates the time course of non-homologous end-joining DNA repair. We explain our findings under terms of DNA injury-induced membrane receptor tyrosine kinase activation and retrograde trafficking. In addition, we interpret the phenomenon of retrograde co-trafficking of gap junction connexons stimulated by ErbB receptor activation.

Clustered localization of EGFRvIII in glioblastoma cells as detected by high precision localization microscopy

For receptor tyrosine kinases supramolecular organization on the cell membrane is critical for their function. Super-resolution fluorescence microscopy techniques have offered new opportunities for the analysis of single receptor localization. Here, we analysed the cluster formation of the epidermal growth factor receptor variant III (EGFRvIII), a deletion variant which is expressed in glioblastoma. The constitutively activated variant EGFRvIII is expressed in cells with an egfr gene amplification and is thought to enhance the tumorigenic potential especially of glioblastoma cells. Due to the lack of an adequate model system, it is still unclear how endogenous EGFRvIII expression alters cellular signalling and if it is organized in clusters like the wild type receptor. We have recently described the establishment of two pairs of iso-genetic cell lines (BS153 and DKMG), displaying endogenous EGFRvIII expression or not. Using these cell lines we investigated single receptor localization of EGFRvIII by high precision localization microscopy. Cluster analysis revealed that EGFRvIII is present in clusters on the surface of the cells, with about 60% or even more receptor molecules being assembled in clusters of approximately 100 nm in diameter whereby the cluster definition was iteratively determined. The signal to signal distance may indicate dimer formation while signal quantification indicates 1 × 10⁶-5 × 10⁶ EGFRvIII molecules per cell. Altogether, these data give unique insights into the membrane surface localization of EGFRvIII in glioblastoma cells. These insights will help to unveil the function of this tumour associated receptor variant which might lead to a better understanding of glioblastoma and therefore could lead to improved therapy approaches.

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

We used coarse-grained molecular dynamics simulations to examine the effects of transmembrane α-helical WALP peptides on the behavior of four-component lipid mixtures. These mixtures contain a high-melting temperature (high-Tm) lipid, a nanodomain-inducing low-Tm lipid, a macrodomain-inducing low-Tm lipid and cholesterol to model the outer leaflet of cell plasma membranes. In a series of simulations, we incrementally replace the nanodomain-inducing low-Tm lipid by the macrodomain-inducing low-Tm lipid and measure how lipid and phase properties are altered by the addition of WALPs of different length. Regardless of the ratio of the two low-Tm lipids, shorter WALPs increase domain size and all WALPs increase domain alignment between the two leaflets. These effects are smallest for the longest WALP tested, and increase with increasing WALP concentration. Thus, our simulations explain the experimental observation that WALPs induce macroscopic domains in otherwise nanodomain-forming lipid-only mixtures (unpublished). Since the cell plasma membrane contains a large fraction of transmembrane proteins, these findings link the behavior of lipid-only model membranes in vitro to phase behavior in vivo.

Phospholipid profiling identifies acyl chain elongation as a ubiquitous trait and potential target for the treatment of lung squamous cell carcinoma

Lung cancer is the leading cause of cancer death. Beyond first line treatment, few therapeutic options are available, particularly for squamous cell carcinoma (SCC). Here, we have explored the phospholipidomes of 30 human SCCs and found that they almost invariably (in 96.7% of cases) contain phospholipids with longer acyl chains compared to matched normal tissues. This trait was confirmed using in situ 2D-imaging MS on tissue sections and by phospholipidomics of tumor and normal lung tissue of the L-IkkαKA/KA mouse model of lung SCC. In both human and mouse, the increase in acyl chain length in cancer tissue was accompanied by significant changes in the expression of acyl chain elongases (ELOVLs). Functional screening of differentially expressed ELOVLs by selective gene knockdown in SCC cell lines followed by phospholipidomics revealed ELOVL6 as the main elongation enzyme responsible for acyl chain elongation in cancer cells. Interestingly, inhibition of ELOVL6 drastically reduced colony formation of multiple SCC cell lines in vitro and significantly attenuated their growth as xenografts in vivo in mouse models. These findings identify acyl chain elongation as one of the most common traits of lung SCC discovered so far and pinpoint ELOVL6 as a novel potential target for cancer intervention.

Rigidity of transmembrane proteins determines their cluster shape

Protein aggregation in cell membrane is vital for the majority of biological functions. Recent experimental results suggest that transmembrane domains of proteins such as α-helices and β-sheets have different structural rigidities. We use molecular dynamics simulation of a coarse-grained model of protein-embedded lipid membranes to investigate the mechanisms of protein clustering. For a variety of protein concentrations, our simulations under thermal equilibrium conditions reveal that the structural rigidity of transmembrane domains dramatically affects interactions and changes the shape of the cluster. We have observed stable large aggregates even in the absence of hydrophobic mismatch, which has been previously proposed as the mechanism of protein aggregation. According to our results, semiflexible proteins aggregate to form two-dimensional clusters, while rigid proteins, by contrast, form one-dimensional string-like structures. By assuming two probable scenarios for the formation of a two-dimensional triangular structure, we calculate the lipid density around protein clusters and find that the difference in lipid distribution around rigid and semiflexible proteins determines the one- or two-dimensional nature of aggregates. It is found that lipids move faster around semiflexible proteins than rigid ones. The aggregation mechanism suggested in this paper can be tested by current state-of-the-art experimental facilities.

Effects of protein crowding on membrane systems

Cellular membranes are typically decorated with a plethora of embedded and adsorbed macromolecules, e.g. proteins, that participate in numerous vital processes. With typical surface densities of 30,000 proteins per μm(2) cellular membranes are indeed crowded places that leave only few nanometers of private space for individual proteins. Here, we review recent advances in our understanding of protein crowding in membrane systems. We first give a brief overview on state-of-the-art approaches in experiment and simulation that are frequently used to study crowded membranes. After that, we review how crowding can affect diffusive transport of proteins and lipids in membrane systems. Next, we discuss lipid and protein sorting in crowded membrane systems, including effects like protein cluster formation, phase segregation, and lipid droplet formation. Subsequently, we highlight recent progress in uncovering crowding-induced conformational changes of membranes, e.g. membrane budding and vesicle formation. Finally, we give a short outlook on potential future developments in the field of crowded membrane systems. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

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.

Influenza A virus hemagglutinin and neuraminidase mutually accelerate their apical targeting through clustering of lipid rafts

In polarized epithelial cells, influenza A virus hemagglutinin (HA) and neuraminidase (NA) are intrinsically associated with lipid rafts and target the apical plasma membrane for viral assembly and budding. Previous studies have indicated that the transmembrane domain (TMD) and cytoplasmic tail (CT) of HA and NA are required for association with lipid rafts, but the raft dependencies of their apical targeting are controversial. Here, we show that coexpression of HA with NA accelerated their apical targeting through accumulation in lipid rafts. HA was targeted to the apical plasma membrane even when expressed alone, but the kinetics was much slower than that of HA in infected cells. Coexpression experiments revealed that apical targeting of HA and NA was accelerated by their coexpression. The apical targeting of HA was also accelerated by coexpression with M1 but not M2. The mutations in the outer leaflet of the TMD and the deletion of the CT in HA and NA that reduced their association with lipid rafts abolished the acceleration of their apical transport, indicating that the lipid raft association is essential for efficient apical trafficking of HA and NA. An in situ proximity ligation assay (PLA) revealed that HA and NA were accumulated and clustered in the cytoplasmic compartments only when both were associated with lipid rafts. Analysis with mutant viruses containing nonraft HA/NA confirmed these findings. We further analyzed lipid raft markers by in situ PLA and suggest a possible mechanism of the accelerated apical transport of HA and NA via clustering of lipid rafts.

IMPORTANCE

Lipid rafts serve as sites for viral entry, particle assembly, and budding, leading to efficient viral replication. The influenza A virus utilizes lipid rafts for apical plasma membrane targeting and particle budding. The hemagglutinin (HA) and neuraminidase (NA) of influenza virus, key players for particle assembly, contain determinants for apical sorting and lipid raft association. However, it remains to be elucidated how lipid rafts contribute to the apical trafficking and budding. We investigated the relation of lipid raft association of HA and NA to the efficiency of apical trafficking. We show that coexpression of HA and NA induces their accumulation in lipid rafts and accelerates their apical targeting, and we suggest that the accelerated apical transport likely occurs by clustering of lipid rafts at the TGN. This finding provides the first evidence that two different raft-associated viral proteins induce lipid raft clustering, thereby accelerating apical trafficking of the viral proteins.

An investigation of the effect of membrane curvature on transmembrane-domain dependent protein sorting in lipid bilayers

Sorting of membrane proteins within the secretory pathway of eukaryotic cells is a complex process involving discrete sorting signals as well as physico-chemical properties of the transmembrane domain (TMD). Previous work demonstrated that tail-anchored (TA) protein sorting at the interface between the Endoplasmic Reticulum (ER) and the Golgi complex is exquisitely dependent on the length and hydrophobicity of the transmembrane domain, and suggested that an imbalance between TMD length and bilayer thickness (hydrophobic mismatch) could drive long TMD-containing proteins into curved membrane domains, including ER exit sites, with consequent export of the mismatched protein out of the ER. Here, we tested a possible role of curvature in TMD-dependent sorting in a model system consisting of Giant Unilamellar Vesicles (GUVs) from which narrow membrane tubes were pulled by micromanipulation. Fluorescent TA proteins differing in TMD length were incorporated into GUVs of uniform lipid composition or made of total ER lipids, and TMD-dependent sorting and diffusion, as well as the bending rigidity of bilayers made of microsomal lipids, were investigated. Long and short TMD-containing constructs were inserted with similar orientation, diffused equally rapidly in GUVs and in tubes pulled from GUVs, and no difference in their final distribution between planar and curved regions was detected. These results indicate that curvature alone is not sufficient to drive TMD-dependent sorting at the ER-Golgi interface, and set the basis for the investigation of the additional factors that must be required.

Reduced lateral mobility of lipids and proteins in crowded membranes

Coarse-grained molecular dynamics simulations of the E. coli outer membrane proteins FhuA, LamB, NanC, OmpA and OmpF in a POPE/POPG (3∶1) bilayer were performed to characterise the diffusive nature of each component of the membrane. At small observation times (<10 ns) particle vibrations dominate phospholipid diffusion elevating the calculated values from the longer time-scale bulk value (>50 ns) of 8.5×10⁻⁷ cm² s⁻¹. The phospholipid diffusion around each protein was found to vary based on distance from protein. An asymmetry in the diffusion of annular lipids in the inner and outer leaflets was observed and correlated with an asymmetry in charged residues in the vicinity of the inner and outer leaflet head-groups. Protein rotational and translational diffusion were also found to vary with observation time and were inversely correlated with the radius of gyration of the protein in the plane of the bilayer. As the concentration of protein within the bilayer was increased, the overall mobility of the membrane decreased reflected in reduced lipid diffusion coefficients for both lipid and protein components. The increase in protein concentration also resulted in a decrease in the anomalous diffusion exponent α of the lipid. Formation of extended clusters and networks of proteins led to compartmentalisation of lipids in extreme cases.

Biological membranes are selective barriers which control the entry/exit of molecules to/from the interior of a cell. They are composed of a lipid bilayer in which are embedded many membrane proteins. Whilst the individual components of membranes are relatively well characterised, the lateral organization and dynamics of the membrane remain less well understood. The lateral mobility of constituent membrane species affects many processes, including how quickly proteins complexes form and protein recruitment occurs, how quickly lipids can be modified/lysed, and the formation of disordered and ordered microdomains. Biological membranes can contain as much as 50% protein. The dynamics of these crowded environments differ greatly from the sparsely populated membranes often studied in silico or in vitro. We use molecular dynamics computer simulations to quantify how mobility within the membrane decreases as the protein concentration increases. We calculate a baseline diffusion of both lipids and selected bacterial outer membrane proteins in the simplest of systems, namely a single protein in a large lipid bilayer patch. In this case diffusion can be correlated with the size of the protein. We observe how proteins affect the mobility of adjacent lipids. As the protein concentration within our systems increases we show that diffusion of both the proteins and lipids is reduced.

Hydrophobic mismatch and lipid sorting near OmpA in mixed bilayers: atomistic and coarse-grained simulations

To understand the effects of lipid composition on membrane protein function in a mixture as complex as a biomembrane, one must know whether the lipid composition local to the protein differs from the mean lipid composition. In this study, we simulated the transmembrane domain of a β-barrel protein, OmpA, in mixtures of lipids of different tail lengths under conditions of negative hydrophobic mismatch, i.e., local bilayer thinning. We modeled the influence of OmpA on the local lipid composition both at a coarse-grained (CG) resolution using conventional molecular dynamics, and at an atomistic resolution within the semi-grand canonical ensemble using mutation moves to rapidly approach an equilibrium lateral distribution of lipids. Moderate enrichment of the shorter tail component (either DDPC in DDPC/DMPC mixtures or DMPC in DMPC/DSPC mixtures) extending 2-3 nm away from the protein surface was observed with both the atomistic and CG models. The similarity in trends suggests that the more computationally economical CG models capture the essential features of lipid sorting induced by hydrophobic mismatch.

The MARVEL transmembrane motif of occludin mediates oligomerization and targeting to the basolateral surface in epithelia

Occludin (Ocln), a MARVEL-motif-containing protein, is found in all tight junctions. MARVEL motifs are comprised of four transmembrane helices associated with the localization to or formation of diverse membrane subdomains by interacting with the proximal lipid environment. The functions of the Ocln MARVEL motif are unknown. Bioinformatics sequence- and structure-based analyses demonstrated that the MARVEL domain of Ocln family proteins has distinct evolutionarily conserved sequence features that are consistent with its basolateral membrane localization. Live-cell microscopy, fluorescence resonance energy transfer (FRET) and bimolecular fluorescence complementation (BiFC) were used to analyze the intracellular distribution and self-association of fluorescent-protein-tagged full-length human Ocln or the Ocln MARVEL motif excluding the cytosolic C- and N-termini (amino acids 60-269, FP-MARVEL-Ocln). FP-MARVEL-Ocln efficiently arrived at the plasma membrane (PM) and was sorted to the basolateral PM in filter-grown polarized MDCK cells. A series of conserved aromatic amino acids within the MARVEL domain were found to be associated with Ocln dimerization using BiFC. FP-MARVEL-Ocln inhibited membrane pore growth during Triton-X-100-induced solubilization and was shown to increase the membrane-ordered state using Laurdan, a lipid dye. These data demonstrate that the Ocln MARVEL domain mediates self-association and correct sorting to the basolateral membrane.

Protein sorting and membrane-mediated interactions

Sorting of membrane proteins is of vital importance for living cells. Indeed, roughly one-third of a eukaryotic cell's proteome consists of peripheral and transmembrane proteins. These need to be properly distributed and dynamically maintained at distinct locations in the compartmentalized cell, and one may wonder how proteins determine where, when, and how to travel to reach a specific organelle. While specific binary interactions between proteins have been invoked in explaining the trafficking and sorting processes, a more active role of lipids in this context has become visible in recent years. In particular, membrane-mediated interactions have been suggested to serve as a robust physicochemical mechanism to facilitate protein sorting. Here, we will review some recent insights into these aspects.

Aggregation of model membrane proteins, modulated by hydrophobic mismatch, membrane curvature, and protein class

Aggregation of transmembrane proteins is important for many biological processes, such as protein sorting and cell signaling, and also for in vitro processes such as two-dimensional crystallization. We have used large-scale simulations to study the lateral organization and dynamics of lipid bilayers containing multiple inserted proteins. Using coarse-grained molecular dynamics simulations, we have studied model membranes comprising ∼7000 lipids and 16 identical copies of model cylindrical proteins of either α-helical or β-barrel types. Through variation of the lipid tail length and hence the degree of hydrophobic mismatch, our simulations display levels of protein aggregation ranging from negligible to extensive. The nature and extent of aggregation are shown to be influenced by membrane curvature and the shape or orientation of the protein. Interestingly, a model β-barrel protein aggregates to form one-dimensional strings within the bilayer plane, whereas a model α-helical bundle forms two-dimensional clusters. Overall, it is clear that the nature and extent of membrane protein aggregation is dependent on several aspects of the proteins and lipids, including hydrophobic mismatch, protein class and shape, and membrane curvature.

Lateral sorting in model membranes by cholesterol-mediated hydrophobic matching

Theoretical studies predict hydrophobic matching between transmembrane domains of proteins and bilayer lipids to be a physical mechanism by which membranes laterally self-organize. We now experimentally study the direct consequences of mismatching of transmembrane peptides of different length with bilayers of different thicknesses at the molecular level. In both model membranes and simulations we show that cholesterol critically constrains structural adaptations at the peptide-lipid interface under mismatch. These constraints translate into a sorting potential and lead to selective lateral segregation of peptides and lipids according to their hydrophobic length.

Dynamic structure formation of peripheral membrane proteins

Using coarse-grained membrane simulations we show here that peripheral membrane proteins can form a multitude of higher-order structures due to membrane-mediated interactions. Peripheral membrane proteins characteristically perturb the lipid bilayer in their vicinity which supports the formation of protein assemblies not only within the same but surprisingly also across opposing leaflets of a bilayer. In addition, we also observed the formation of lipid-protein domains on heteregeneous membranes. The clustering ability of proteins, as quantified via the potential of mean force, is enhanced when radius and hydrophobic penetration depth of the proteins increases. Based on our data, we propose that membrane-mediated cluster formation of peripheral proteins supports protein assembly in vivo and hence may play a pivotal role in the formation of templates for signaling cascades and in the emergence of transport intermediates in the secretory pathway.

Eukaryotic cells are subdivided into a variety of compartments by membranes, i.e. by lipid bilayers into which a multitude of proteins are embedded. About 30% of all protein species in a cell are associated with membranes to perform vital functions, e.g. in signaling and transport pathways. A plethora of membrane-associated proteins, so-called peripheral membrane proteins, penetrate only one monolayer whereas transmembrane proteins span the entire thickness of a lipid bilayer. Despite the crucial role of peripheral membrane proteins in many cellular processes, little is known about the general principles of how they can form higher-order structures that eventually are biologically active units. Using a simulation approach we show and quantify here the influence of membrane-mediated interactions between peripheral membrane proteins. These interactions most likely play a pivotal role in forming higher-order structures in living cells. In particular, template formation in signaling cascades and coat assembly during the formation of transport intermediates may strongly benefit from the phenomena described here.

From Coarse Grained to Atomistic: A Serial Multiscale Approach to Membrane Protein Simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and the interactions of membrane protein/lipid complexes at a reduced level of representation, allowing longer and larger simulations. We describe a fragment-based protocol for converting membrane simulation systems, comprising a membrane protein embedded in a phospholipid bilayer, from coarse-grained to atomistic resolution, for further refinement and analysis via atomistic simulations. Overall, this provides a method for generating an accurate and well equilibrated membrane protein/lipid complex. We exemplify the protocol using the acid-sensing/amiloride-sensitive ion channel protein (ASIC) channel protein, a trimeric integral membrane protein. The method is further evaluated using a test set of 10 different membrane proteins of differing size and complexity. Simulations are assessed in terms of protein conformational drift, lipid/protein interactions, and lipid dynamics.