Energetics of hydrophobic matching in lipid-protein interactions

Lipid chain length modulates the activity of transmembrane proteins by mismatch between the hydrophobic span of the protein and that of the lipid membrane. Relative binding affinities of lipids with different chain lengths are used to estimate the excess free energy of lipid-protein interaction that arises from hydrophobic mismatch. For a wide range of integral proteins and peptides, the energy cost is much less than the elastic penalty of fully stretching or compressing the lipid chains to achieve complete hydrophobic matching. The chain length dependences of the free energies of lipid association are described by a model that combines elastic chain extension with a free energy term that depends linearly on the extent of residual mismatch. The excess free energy densities involved lie in the region of 0.5-2.0 k_BT.nm⁻². Values of this size could arise from exposure of hydrophobic groups to polar portions of the lipid or protein, but not directly to water, or alternatively from changes in tilt of the transmembrane helices that are energetically comparable to those activating mechanosensitive channels. The influence of hydrophobic mismatch on dimerization of transmembrane helices and their transfer between lipid vesicles, and on shifts in chain-melting transitions of lipid bilayers by incorporated proteins, is analyzed by using the same thermodynamic model. Segmental order parameters confirm that elastic lipid chain distortions are insufficient to compensate fully for the mismatch, but the dependence on chain length with tryptophan-anchored peptides requires that the free energy density of hydrophobic mismatch should increase with increasing extent of mismatch.     

Lipid-protein interactions

Intrinsic membrane proteins are solvated by a shell of lipid molecules interacting with the membrane-penetrating surface of the protein; these lipid molecules are referred to as annular lipids. Lipid molecules are also found bound between transmembrane α-helices; these are referred to as non-annular lipids. Annular lipid binding constants depend on fatty acyl chain length, but the dependence is less than expected from models based on distortion of the lipid bilayer alone. This suggests that hydrophobic matching between a membrane protein and the surrounding lipid bilayer involves some distortion of the transmembrane α-helical bundle found in most membrane proteins, explaining the importance of bilayer thickness for membrane protein function. Annular lipid binding constants also depend on the structure of the polar headgroup region of the lipid, and hotspots for binding anionic lipids have been detected on some membrane proteins; binding of anionic lipid molecules to these hotspots can be functionally important. Binding of anionic lipids to non-annular sites on membrane proteins such as the potassium channel KcsA can also be important for function. It is argued that the packing preferences of the membrane-spanning α-helices in a membrane protein result in a structure that matches nicely with that of the surrounding lipid bilayer, so that lipid and protein can meet without either having to change very much.     

How protein transmembrane segments sense the lipid environment

Integral membrane proteins have central roles in a vast number of vital cellular processes. A structural feature that most membrane proteins have in common is the presence of one or more alpha-helices with which they traverse the lipid bilayer. Because of the interaction with the surrounding lipids, the organization of these transmembrane helices will be sensitive to lipid properties like lateral packing, hydrophobic thickness, and headgroup charge. The helices may adapt to the lipids in different ways, which in turn can influence the structure and function of the intact membrane protein. In this review, we will focus on how the lipid environment influences two specific properties of transmembrane segments: their lateral association and their tilt with respect to the bilayer normal.     

Functional competition within a membrane: Lipid recognition vs. transmembrane helix oligomerization

Binding of specific lipids to large, polytopic membrane proteins is well described, and it is clear that such lipids are crucial for protein stability and activity. In contrast, binding of defined lipid species to individual transmembrane helices and regulation of transmembrane helix monomer–oligomer equilibria by binding of distinct lipids is a concept, which has emerged only lately. Lipids bind to single-span membrane proteins, both in the juxta-membrane region as well as in the hydrophobic membrane core. While some interactions counteract transmembrane helix oligomerization, in other cases lipid binding appears to enhance oligomerization. As reversible oligomerization is involved in activation of many membrane proteins, binding of defined lipids to single-span transmembrane proteins might be a mechanism to regulate and/or fine-tune the protein activity. But how could lipid binding trigger the activity of a protein? How can binding of a single lipid molecule to a transmembrane helix affect the structure of a transmembrane helix oligomer, and consequently its signaling state? These questions are discussed in the present article based on recent results obtained with simple, single-span transmembrane proteins. This article is part of a Special Issue entitled: Lipid–protein interactions.     

YidC, an assembly site for polytopic Escherichia coli membrane proteins located in immediate proximity to the SecYE translocon and lipids

Like its mitochondrial homolog Oxa1p, the inner membrane protein YidC of Escherichia coli is involved in the integration of membrane proteins. We have analyzed individual insertion steps of the polytopic E. coli membrane protein MtlA targeted as ribosome-nascent chain complexes to inner membrane vesicles. YidC can accommodate at least the first two transmembrane segments of MtlA at the protein lipid interface and retain them even though the length of the nascent chain would amply allow insertion into membrane lipids. An even longer insertion intermediate of MtlA is described that still has the first transmembrane helix bound to YidC while the third contacts SecE and YidC during integration. Our findings suggest that YidC forms a contiguous integration unit with the SecYE translocon and functions as an assembly site for polytopic membrane proteins mediating the formation of helix bundles prior to their release into the membrane lipids.     

Lipid-protein interactions in biological membranes: a structural perspective

Lipid molecules bound to membrane proteins are resolved in some high-resolution structures of membrane proteins. An analysis of these structures provides a framework within which to analyse the nature of lipid-protein interactions within membranes. Membrane proteins are surrounded by a shell or annulus of lipid molecules, equivalent to the solvent layer surrounding a water-soluble protein. The lipid bilayer extends right up to the membrane protein, with a uniform thickness around the protein. The surface of a membrane protein contains many shallow grooves and protrusions to which the fatty acyl chains of the surrounding lipids conform to provide tight packing into the membrane. An individual lipid molecule will remain in the annular shell around a protein for only a short period of time. Binding to the annular shell shows relatively little structural specificity. As well as the annular lipid, there is evidence for other lipid molecules bound between the transmembrane α-helices of the protein; these lipids are referred to as non-annular lipids. The average thickness of the hydrophobic domain of a membrane protein is about 29 Å, with a few proteins having significantly smaller or greater thicknesses than the average. Hydrophobic mismatch between a membrane protein and the surrounding lipid bilayer generally leads to only small changes in membrane thickness. Possible adaptations in the protein to minimise mismatch include tilting of the helices and rotation of side chains at the ends of the helices. Packing of transmembrane α-helices is dependent on the chain length of the surrounding phospholipids. The function of membrane proteins is dependent on the thickness of the surrounding lipid bilayer, sometimes on the presence of specific, usually anionic, phospholipids, and sometimes on the phase of the phospholipid.     

Lipid-protein interactions in biological membranes: a structural perspective

Lipid molecules bound to membrane proteins are resolved in some high-resolution structures of membrane proteins. An analysis of these structures provides a framework within which to analyse the nature of lipid-protein interactions within membranes. Membrane proteins are surrounded by a shell or annulus of lipid molecules, equivalent to the solvent layer surrounding a water-soluble protein. The lipid bilayer extends right up to the membrane protein, with a uniform thickness around the protein. The surface of a membrane protein contains many shallow grooves and protrusions to which the fatty acyl chains of the surrounding lipids conform to provide tight packing into the membrane. An individual lipid molecule will remain in the annular shell around a protein for only a short period of time. Binding to the annular shell shows relatively little structural specificity. As well as the annular lipid, there is evidence for other lipid molecules bound between the transmembrane alpha-helices of the protein; these lipids are referred to as non-annular lipids. The average thickness of the hydrophobic domain of a membrane protein is about 29 A, with a few proteins having significantly smaller or greater thicknesses than the average. Hydrophobic mismatch between a membrane protein and the surrounding lipid bilayer generally leads to only small changes in membrane thickness. Possible adaptations in the protein to minimise mismatch include tilting of the helices and rotation of side chains at the ends of the helices. Packing of transmembrane alpha-helices is dependent on the chain length of the surrounding phospholipids. The function of membrane proteins is dependent on the thickness of the surrounding lipid bilayer, sometimes on the presence of specific, usually anionic, phospholipids, and sometimes on the phase of the phospholipid.     

Sphingomyelin chain length influences the distribution of GPI-anchored proteins in rafts in supported lipid bilayers

Glycosyl-phosphatidylinositol (GPI)-anchored proteins are enriched in cholesterol- and sphingolipid-rich lipid rafts within the membrane. Rafts are known to have roles in cellular organization and function, but little is understood about the factors controlling the distribution of proteins in rafts. We have used atomic force microscopy to directly visualize proteins in supported lipid bilayers composed of equimolar sphingomyelin, dioleoyl-sn-glycero-3-phosphocholine and cholesterol. The transmembrane anchored angiotensin converting enzyme (TM-ACE) was excluded from the liquid ordered raft domains. Replacement of the transmembrane and cytoplasmic domains of TM-ACE with a GPI anchor (GPI-ACE) promoted the association of the protein with rafts in the bilayers formed with brain sphingomyelin (mainly C18:0). Association with the rafts did not occur if the shorter chain egg sphingomyelin (mainly C16:0) was used. The distribution of GPI-anchored proteins in supported lipid bilayers was investigated further using membrane dipeptidase (MDP) whose GPI anchor contains distearoyl phosphatidylinositol. MDP was also excluded from rafts when egg sphingomyelin was used but associated with raft domains formed using brain sphingomyelin. The effect of sphingomyelin chain length on the distribution of GPI-anchored proteins in rafts was verified using synthetic palmitoyl or stearoyl sphingomyelin. Both GPI-ACE and MDP only associated with the longer chain stearoyl sphingomyelin rafts. These data obtained using supported lipid bilayers provide the first direct evidence that the nature of the membrane-anchoring domain influences the association of a protein with lipid rafts and that acyl chain length hydrophobic mismatch influences the distribution of GPI-anchored proteins in rafts.     

Characterization of a membrane protein folding motif, the Ser zipper, using designed peptides

Polar residues play important roles in the association of transmembrane helices and the stabilities of membrane proteins. Although a single Ser residue in a transmembrane helix is unable to mediate a strong association of the helices, the cooperative interactions of two or more appropriately placed serine hydroxyl groups per helix has been hypothesized to allow formation of a "serine zipper" that can stabilize transmembrane helix association. In particular, a heptad repeat Sera Xxx Xxx Leud Xxx Xxx Xxx (Xxx is a hydrophobic amino acid) appears in both antiparallel helical pairs of polytopic membrane proteins as well as the parallel helical dimerization motif found in the murine erythropoietin receptor. To examine the intrinsic conformational preferences of this motif independent of its context within a larger protein, we synthesized a peptide containing three copies of a SeraLeud heptad motif. Computational results are consistent with the designed peptide adopting either a parallel or antiparallel structure, and conformational search calculations yield the parallel dimer as the lowest energy configuration, which is also significantly more stable than the parallel trimer. Analytical ultracentrifugation indicated that the peptide exists in a monomer-dimer equilibrium in dodecylphosphocholine micelles. Thiol disulfide interchange studies showed a preference for forming parallel dimers in micelles. In phospholipid vesicles, only the parallel dimer was formed. The stability of the SerZip peptide was studied in vesicles prepared from phosphatidylcholine (PC) lipids of different chain length: POPC (C16:0C18:1 PC) and DLPC (C12:0PC). The stability was greater in POPC, which has a good match between the length of the hydrophobic region of the peptide and the bilayer length. Finally, mutation to Ala of the Ser residues in the SerZip motif gave rise to a relatively small decrease in the stability of the dimer, indicating that packing interactions rather than hydrogen-bonding provided the primary driving force for association.     

Electron spin resonance in membrane research: protein–lipid interactions from challenging beginnings to state of the art

Conventional electron paramagnetic resonance (EPR) spectra of lipids that are spin-labelled close to the terminal methyl end of the acyl chains are able to resolve the lipids directly contacting the protein from those in the fluid bilayer regions of the membrane. This allows determination of both the stoichiometry of lipid–protein interaction (i.e., number of lipid sites at the protein perimeter) and the selectivity of the protein for different lipid species (i.e., association constants relative to the background lipid). Spin-label EPR data are summarised for 20 or more different transmembrane peptides and proteins, and 7 distinct species of lipids. Lineshape simulations of the two-component conventional spin-label EPR spectra allow estimation of the rate at which protein-associated lipids exchange with those in the bulk fluid regions of the membrane. For lipids that do not display a selectivity for the protein, the intrinsic off-rates for exchange are in the region of 10 MHz: less than 10× slower than the rates of diffusive exchange in fluid lipid membranes. Lipids with an affinity for the protein, relative to the background lipid, have off-rates for leaving the protein that are correspondingly slower. Non-linear EPR, which depends on saturation of the spectrum at high radiation intensities, is optimally sensitive to dynamics on the timescale of spin-lattice relaxation, i.e., the microsecond regime. Both progressive saturation and saturation transfer EPR experiments provide definitive evidence that lipids at the protein interface are exchanging on this timescale. The sensitivity of non-linear EPR to low frequencies of spin exchange also allows the location of spin-labelled membrane protein residues relative to those of spin-labelled lipids, in double-labelling experiments.     

Molecular genetic and biochemical approaches for defining lipid-dependent membrane protein folding

We provide an overview of lipid-dependent polytopic membrane protein folding and topogenesis. Lipid dependence of this process was determined by employing Escherichia coli cells in which specific lipids can be eliminated, substituted, tightly titrated or controlled temporally during membrane protein synthesis and assembly. The secondary transport protein lactose permease (LacY) was used to establish general principles underlying the molecular basis of lipid-dependent effects on protein domain folding, protein transmembrane domain (TM) orientation, and function. These principles were then extended to several other secondary transport proteins of E. coli. The methods used to follow proper conformational organization of protein domains and the topological organization of protein TMs in whole cells and membranes are described. The proper folding of an extramembrane domain of LacY that is crucial for energy dependent uphill transport function depends on specific lipids acting as non-protein molecular chaperones. Correct TM topogenesis is dependent on charge interactions between the cytoplasmic surface of membrane proteins and a proper balance of the membrane surface net charge defined by the lipid head groups. Short-range interactions between the nascent protein chain and the translocon are necessary but not sufficient for establishment of final topology. After release from the translocon short-range interactions between lipid head groups and the nascent protein chain, partitioning of protein hydrophobic domains into the membrane bilayer, and long-range interactions within the protein thermodynamically drive final membrane protein organization. Given the diversity of membrane lipid compositions throughout nature, it is tempting to speculate that during the course of evolution the physical and chemical properties of proteins and lipids have co-evolved in the context of the lipid environment of membrane systems in which both are mutually dependent on each other for functional organization of proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.     

Association of spin-labeled lipids with beta-barrel proteins from the outer membrane of Escherichia coli

The interaction of spin-labeled lipids with beta-barrel transmembrane proteins has been studied by the electron spin resonance (ESR) methods developed for alpha-helical integral proteins. The outer membrane protein OmpA and the ferrichrome-iron receptor FhuA from the outer membrane of Escherichia coli were reconstituted in bilayers of dimyristoylphosphatidylglycerol. The ESR spectra from phosphatidylglycerol spin labeled on the 14-C atom of the sn-2 chain contain a second component from motionally restricted lipids contacting the intramembranous surface of the beta-barrel, in addition to that from the fluid bilayer lipids. The stoichiometry of motionally restricted lipids, 11 and 32 lipids/monomer for OmpA and FhuA, respectively, is constant irrespective of the total lipid/protein ratio. It is proportional to the number of transmembrane beta-strands, eight for OmpA and 22 for FhuA, and correlates reasonably well with the intramembranous perimeter of the protein. Spin-labeled lipids with different polar headgroups display a differential selectivity of interaction with the two proteins. The more pronounced pattern of lipid selectivity for FhuA than for OmpA correlates with the preponderance of positively charged residues facing the lipids in the extensions of the beta-sheet and shorter interconnecting loops on the extracellular side of FhuA.     

Protein-lipid interactions with Fusobacterium nucleatum major outer membrane protein FomA: spin-label EPR and polarized infrared spectroscopy

FomA, the major outer membrane protein of Fusobacterium nucleatum, was expressed and purified in Escherichia coli and reconstituted from detergent in bilayer membranes of phosphatidylcholines with chain lengths from C(12:0) to C(17:0). The conformation and orientation of membrane-incorporated FomA were determined from polarized, attenuated total reflection, infrared (IR) spectroscopy, and lipid-protein interactions with FomA were characterized by using electron paramagnetic resonance (EPR) spectroscopy of spin-labeled lipids. Approximately 190 residues of membranous FomA are estimated to be in a beta-sheet configuration from IR band fitting, which is consistent with a 14-strand transmembrane beta-barrel structure. IR dichroism of FomA indicates that the beta-strands are tilted by approximately 45 degrees relative to the sheet/barrel axis and that the order parameter of the latter displays a discontinuity corresponding to hydrophobic matching with fluid C(13:0) lipid chains. The stoichiometry ( N b = 23 lipids/monomer) of lipid-protein interaction from EPR demonstrates that FomA is not trimeric in membranes of diC(14:0) phosphatidylcholine and is consistent with a monomeric beta-barrel of 14-16 strands. The pronounced selectivity of interaction found with anionic spin-labeled lipids places basic residues of the protein in the vicinity of the polar-apolar membrane interfaces, consistent with current topology models. Comparison with similar data from the 8- to 22-stranded E. coli outer membrane proteins, OmpA, OmpG, and FhuA, supports the above conclusions.     

Sphingomyelin chain length influences the distribution of GPI-anchored proteins in rafts in supported lipid bilayers

Glycosyl-phosphatidylinositol (GPI)-anchored proteins are enriched in cholesterol- and sphingolipid-rich lipid rafts within the membrane. Rafts are known to have roles in cellular organization and function, but little is understood about the factors controlling the distribution of proteins in rafts. We have used atomic force microscopy to directly visualize proteins in supported lipid bilayers composed of equimolar sphingomyelin, dioleoyl-sn-glycero-3-phosphocholine and cholesterol. The transmembrane anchored angiotensin converting enzyme (TM-ACE) was excluded from the liquid ordered raft domains. Replacement of the transmembrane and cytoplasmic domains of TM-ACE with a GPI anchor (GPI-ACE) promoted the association of the protein with rafts in the bilayers formed with brain sphingomyelin (mainly C18:0). Association with the rafts did not occur if the shorter chain egg sphingomyelin (mainly C16:0) was used. The distribution of GPI-anchored proteins in supported lipid bilayers was investigated further using membrane dipeptidase (MDP) whose GPI anchor contains distearoyl phosphatidylinositol. MDP was also excluded from rafts when egg sphingomyelin was used but associated with raft domains formed using brain sphingomyelin. The effect of sphingomyelin chain length on the distribution of GPI-anchored proteins in rafts was verified using synthetic palmitoyl or stearoyl sphingomyelin. Both GPI-ACE and MDP only associated with the longer chain stearoyl sphingomyelin rafts. These data obtained using supported lipid bilayers provide the first direct evidence that the nature of the membrane-anchoring domain influences the association of a protein with lipid rafts and that acyl chain length hydrophobic mismatch influences the distribution of GPI-anchored proteins in rafts.     

Stoichiometry of lipid interactions with transmembrane proteins--Deduced from the 3D structures

The stoichiometry of the first shell of lipids interacting with a transmembrane protein is defined operationally by the population of spin-labeled lipid chains whose motion is restricted directly by the protein. Interaction stoichiometries have been determined experimentally for a wide range of alpha-helical integral membrane proteins by using spin-label ESR spectroscopy. Here, we determine the spatially defined number of first-shell lipids at the hydrophobic perimeter of integral membrane proteins whose 3D structure has been determined by X-ray crystallography and lipid-protein interactions characterized by spin-labeling. Molecular modeling is used to build a single shell of lipids surrounding transmembrane structures derived from the PDB. Constrained energy optimization of the protein-lipid assemblies is performed by molecular mechanics. For relatively small proteins (up to 7-12 transmembrane helices), the geometrical first shell corresponds to that defined experimentally by perturbation of the lipid-chain dynamics. For larger, multi-subunit alpha-helical proteins, the lipids perturbed directly by the protein may either exceed or be less in number than those that can be accommodated at the intramembranous perimeter. In these latter cases, the motionally restricted spin-labeled lipids can be augmented by intercalation, or can correspond to a specific subpopulation at the protein interface, respectively. For monomeric beta-barrel proteins, the geometrical lipid stoichiometry corresponds to that determined from lipid mobility for a 22-stranded barrel, but fewer lipids are motionally restricted than can be accommodated around an eight-stranded barrel. Deviations from the geometrical first shell, in the beta-barrel case, are for the smaller protein with a highly curved barrel.     

Alpha-synuclein association with phosphatidylglycerol probed by lipid spin labels

Alpha-synuclein is a small presynaptic protein, which is linked to the development of Parkinson's disease. Alpha-synuclein partitions between cytosolic and vesicle-bound states, where membrane binding is accompanied by the formation of an amphipathic helix in the N-terminal section of the otherwise unstructured protein. The impact on alpha-synuclein of binding to vesicle-like liposomes has been studied extensively, but far less is known about the impact of alpha-synuclein on the membrane. The interactions of alpha-synuclein with phosphatidylglycerol membranes are studied here by using spin-labeled lipid species and electron spin resonance (ESR) spectroscopy to allow a detailed analysis of the effect on the membrane lipids. Membrane association of alpha-synuclein perturbs the ESR spectra of spin-labeled lipids in bilayers of phosphatidylglycerol but not of phosphatidylcholine. The interaction is inhibited at high ionic strength. The segmental motion is hindered at all positions of spin labeling in the phosphatidylglycerol sn-2 chain, while still preserving the chain flexibility gradient characteristic of fluid phospholipid membranes. Direct motional restriction of the lipid chains, resulting from penetration of the protein into the hydrophobic interior of the membrane, is not observed. Saturation occurs at a protein/lipid ratio corresponding to approximately 36 lipids/protein added. Alpha-synuclein exhibits a selectivity of interaction with different phospholipid spin labels when bound to phosphatidylglycerol membranes in the following order: stearic acid > cardiolipin > phosphatidylcholine > phosphatidylglycerol approximately phosphatidylethanolamine > phosphatidic acid approximately phosphatidylserine > N-acyl phosphatidylethanolamine > diglyceride. Accordingly, membrane-bound alpha-synuclein associates at the interfacial region of the bilayer where it may favor a local concentration of certain phospholipids.     

Life at the border: Adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water‐soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co‐translational and post‐translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large‐scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen‐bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region—between double bonds and carbonyl groups of lipids. These “midpolar” regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side‐chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein–lipid binding.     

Different membrane anchoring positions of tryptophan and lysine in synthetic transmembrane alpha-helical peptides.

Specific interactions of membrane proteins with the membrane interfacial region potentially define protein position with respect to the lipid environment. We investigated the proposed roles of tryptophan and lysine side chains as "anchoring" residues of transmembrane proteins. Model systems were employed, consisting of phosphatidylcholine lipids and hydrophobic alpha-helical peptides, flanked either by tryptophans or lysines. Peptides were incorporated in bilayers of different thickness, and effects on lipid structure were analyzed. Induction of nonbilayer phases and also increases in bilayer thickness were observed that could be explained by a tendency of Trp as well as Lys residues to maintain interactions with the interfacial region. However, effects of the two peptides were remarkably different, indicating affinities of Trp and Lys for different sites at the interface. Our data support a model in which the Trp side chain has a specific affinity for a well defined site near the lipid carbonyl region, while the lysine side chain prefers to be located closer to the aqueous phase, near the lipid phosphate group. The information obtained in this study may further our understanding of the architecture of transmembrane proteins and may prove useful for refining prediction methods for transmembrane segments.     

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.     

Protein–lipid interactions studied with designed transmembrane peptides: role of hydrophobic matching and interfacial anchoring (Review)

Biological membranes are characterized by a heterogeneous composition, which is not only manifested in the wide variety of their components, but also in aspects like the lateral organization, topology, and conformation of proteins and lipids. In bringing about the correct membrane structure, protein–lipid interactions can be expected to play a prominent role. The extent of hydrophobic matching between transmembrane protein segments and lipids potentially constitutes a versatile director of membrane organization, because a tendency to avoid hydrophobic mismatch could result in compensating adaptations such as tilt of the transmembrane segment or segregation into distinct domains. Also, interfacial interactions between lipid headgroups and the aromatic and charged residues that typically flank transmembrane domains may act as an organizing element. In this review, we discuss the numerous model studies that have systematically explored the influence of hydrophobic matching and interfacial anchoring on membrane structure. Designed peptides consisting of a polyleucine or polyleucine/alanine hydrophobic stretch, which is flanked on both sides by tryptophan or lysine residues, reflect the general layout of transmembrane protein segments. It is shown for phosphatidylcholine bilayers and for other model membranes that these peptides adapt a transmembrane topology without extensive peptide or lipid adaptations under conditions of hydrophobic matching, but that significant rearrangements can result from hydrophobic mismatch. Moreover, these effects depend on the nature of the flanking residues, implying a modulation of the mismatch response by interfacial interactions of the flanking residues. The implications of these model studies for the organization of biomembranes are discussed in the context of recent experiments with more complex systems.