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.     

Is the protein/lipid hydrophobic matching principle relevant to membrane organization and functions?

Biological membranes are complex and well-organized multimolecular assemblies composed of a wide variety of protein and lipid molecular species. If such a diversity in protein and lipid polar headgroup structures may easily be related to a large panel of functions, the wide dispersion in acyl chain length and structure which the lipids display is more difficult to understand. It is not required for maintaining bilayer assembly and fluidity. Direct information on the lateral distribution of these various molecular species, on their potential specificity for interaction between themselves and with proteins and on the functional implications of these interactions is also still lacking. Because hydrophobic interactions play a major role in stabilizing membrane structures, we suggest considering the problem from the point of view of the matching of the hydrophobic surface of proteins by the acyl chains of the lipids. After a brief introduction to the hydrophobic matching principle, we will present experimental results which demonstrate the predictive power of the current theories and then, we will introduce the new and important concept of protein/lipid sorting in membranes. Finally, we will show how the hydrophobic matching condition may play a key role in the membrane organization and function.     

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.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.     

Effect of nonbilayer lipids on membrane binding and insertion of the catalytic domain of leader peptidase

Biological membranes contain a substantial amount of "nonbilayer lipids", which have a tendency to form nonlamellar phases. In this study the hypothesis was tested that the presence of nonbilayer lipids in a membrane, due to their overall small headgroup, results in a lower packing density in the headgroup region, which might facilitate the interfacial insertion of proteins. Using the catalytic domain of leader peptidase (delta2-75) from Escherichia coli as a model protein, we studied the lipid class dependence of its insertion and binding. In both lipid monolayers and vesicles, the membrane binding of (catalytically active) delta2-75 was much higher for the nonbilayer lipid DOPE compared to the bilayer lipid DOPC. For the nonbilayer lipids DOG and MGDG a similar effect was observed as for DOPE, strongly suggesting that no specific interactions are involved but that the small headgroups create hydrophobic interfacial insertion sites. On the basis of the results of the monolayer experiments, calculations were performed to estimate the space between the lipid headgroups accessible to the protein. We estimate a maximal size of the insertion sites of 15 +/- 7 A2/lipid molecule for DOPE, relative to DOPC. The size of the insertion sites decreases with an increase in headgroup size. These results show that nonbilayer lipids stimulate the membrane insertion of delta2-75 and support the idea that such lipids create insertion sites by reducing the packing density at the membrane-water interface. It is suggested that PE in the bacterial membrane facilitates membrane insertion of the catalytic domain of leader peptidase, allowing the protein to reach the cleavage site in preproteins.     

Lipid bilayer deformation and the free energy of interaction of a Kv channel gating-modifier toxin

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein.     

The electrostatics of lipid surfaces

Charged lipids constitute a substantial fraction of all membrane lipids. Their charges vary in quantity and distribution within their headgroup regions. In long range interactions, their charges' value and electrostatic potential in the vicinity of the membrane surface can be approximated by the Guy-Chapman theory. This theory treats the interface as a charged structureless plain surrounded by uniform environments. However, if one considers intermolecular interactions, such assumptions need to be revised. The interface is in reality a thick region containing the residual charges of lipid headgroups. Their arrangement depends on the type of lipid present in the membrane. The variety of lipids and their biological functions suggests that charge distribution determines the extent and type of interaction with surface associated molecules. Numerous examples show that protein behavior at the lipid bilayer surface is determined by the type of lipid present, indicating protein specificity towards certain surface locations and local properties (determined by lipid composition) of a particular type. Such specificity is achieved by a combination of electrostatic, hydrophobic and enthropic effects. Comparing lipid biological activity, it can be stated that residual charge distribution is one of the factors of intermolecular recognition leading to the specific interaction of lipid molecules and selected proteins in various processes, particularly those involved with signal transduction pathways. Such specificity enables a variety of processes occurring simultaneously on the same membrane surface to function without cross-reaction interference.     

ESR and monolayer study of the localization of coenzyme Q10 in artificial membranes

The data obtained from the ESR experiments show a complex, depth dependent effect of CoQ10 on the lipid molecules mobility in the bilayer. These effects depend both on its concentration and the temperature. CoQ10 disturbs not only the hydrophobic core of the membrane but also the region close to the hydrophilic headgroups of phospholipids. Both these effects could be explained by the fact that the high hydrophobicity of CoQ10 causes the molecules to position itself in the interior of the bilayer, but at the same time its water seeking headgroup is located close to the region of the polar headgrops of membrane lipids. The presence of CoQ10 in the hydrophobic core has further implications on the properties of membrane intrinsic domain. Results of monolayer experiments indicate that CoQ10 may form aggregates when mixed with PC molecules in the lipid hydrocarbon chain-length dependent manner. CoQ10 is not fully miscible with DMPC or DPPC but it is well miscible with the long-chain DSPC molecules. Our suggestion is that CoQ10 when present in long-chain phospholipid bilayer, interacts with saturated fatty acyl-chains and adapt the structure which allows such interactions: either parallel to the saturated acyl chains or "pseudo-ring" conformation resembling sterol structure.     

Structural Thermodynamics of myr-Src(2–19) Binding to Phospholipid Membranes

Many proteins are anchored to lipid bilayer membranes through a combination of hydrophobic and electrostatic interactions. In the case of the membrane-bound nonreceptor tyrosine kinase Src from Rous sarcoma virus, these interactions are mediated by an N-terminal myristoyl chain and an adjacent cluster of six basic amino-acid residues, respectively. In contrast with the acyl modifications of other lipid-anchored proteins, the myristoyl chain of Src does not match the host lipid bilayer in terms of chain conformation and dynamics, which is attributed to a tradeoff between hydrophobic burial of the myristoyl chain and repulsion of the peptidic moiety from the phospholipid headgroup region. Here, we combine thermodynamic information obtained from isothermal titration calorimetry with structural data derived from ²H, ¹³C, and ³¹P solid-state nuclear magnetic resonance spectroscopy to decipher the hydrophobic and electrostatic contributions governing the interactions of a myristoylated Src peptide with zwitterionic and anionic membranes made from lauroyl (C12:0) or myristoyl (C14:0) lipids. Although the latter are expected to enable better hydrophobic matching, the Src peptide partitions more avidly into the shorter-chain lipid analog because this does not require the myristoyl chain to stretch extensively to avoid unfavorable peptide/headgroup interactions. Moreover, we find that Coulombic and intrinsic contributions to membrane binding are not additive, because the presence of anionic lipids enhances membrane binding more strongly than would be expected on the basis of simple Coulombic attraction.     

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.     

Interaction of an Amphiphilic Peptide with a Phospholipid Bilayer Surface by Molecular Dynamics Simulation Study

Abstract

Corticotropin-releasing factor (CRF) is the principal neuroregulator of adrenocorticotropic hormone (ACTH) secretion. Previous experiments have demonstrated that CRF binds avidly to the surface of single egg phosphatidylcholine vesicles and its amphiphilic secondary structure might play an important role in the function. In this study, the interaction of the residues 13–41 in human CRF with the surface of a DOPC bilayer was investigated by molecular dynamics (MD) simulation in order to understand the role of the membrane surface in the formation of the amphiphilic α helix as well as to determine the effects of the peptide on the lipid bilayer. The model used included 60 DOPC molecules, 1 helical peptide (CRF_13–41) on the bilayer surface, and explicit waters of solvation in the lipid polar head group regions, together with constant-volume periodic boundary conditions in three dimensions. The MD simulation was carried out for 510 ps. In addition, CRF_13–41, initially in a helical form, was simulated in vacuo as a control. The results indicate that while it was completely unstable in vacuo, the peptide helical form was generally maintained on the bilayer surface, but with distortions near the terminal ends. The peptide was confined to the bilayer headgroup/water region, similar to that reported from neutron diffraction measurement of tripeptides bound to the phosphatidylcholine bilayer surface (Ref 1). The amphiphilicity of the peptide matched that of the bilayer headgroup environment, with the hydrophilic side oriented toward water and the hydrophobic side making contact with the bilayer hydrocarbon core. These results support the hypothesis that the amphiphilic environment of a membrane surface is important in the induction of peptide amphiphilic α-helical secondary structure. Two major effects of the peptide on the lipids were found: the first CH₂ segment in the lipid chains was significantly disordered and the lipid headgroup distribution was broadened towards the water region.     

Identification of the hydrophobic thickness of a membrane protein using fluorescence spectroscopy: studies with the mechanosensitive channel MscL

The hydrophobic thickness of a membrane protein is an important parameter, defining how the protein sits within the hydrocarbon core of the lipid bilayer that surrounds it in a membrane. Here we show that Trp scanning mutagenesis combined with fluorescence spectroscopy can be used to define the hydrophobic thickness of a membrane protein. The mechanosensitive channel of large conductance (MscL) contains two transmembrane alpha-helices, of which the second (TM2) is lipid-exposed. The region of TM2 that spans the hydrocarbon core of the bilayer when MscL is reconstituted into bilayers of dioleoylphosphatidylcholine runs from Leu-69 to Leu-92, giving a hydrophobic thickness of ca. 25 A. The results obtained using Trp scanning mutagenesis were confirmed using Cys residues labeled with the N-methyl-amino-7-nitroben-2-oxa-1,3-diazole [NBD] group; both fluorescence emission maxima and fluorescence lifetimes for the NBD group are sensitive to solvent dielectric constant over the range (2-40) thought to span the lipid headgroup region of a lipid bilayer. Changing phospholipid fatty acyl chain lengths from C14 and C24 results in no significant change for the fluorescence of the interfacial residues, suggesting very efficient hydrophobic matching between the protein and the surrounding lipid bilayer.     

Manipulating the folding of membrane proteins: using the bilayer to our advantage

The folding mechanisms of integral membrane proteins have largely eluded detailed study. This is owing to the inherent difficulties in folding these hydrophobic proteins in vitro, which, in turn, reflects the often apparently insurmountable problem of mimicking the natural membrane bilayer with lipid or detergent mixtures. There is, however, a large body of information on lipid properties and, in particular, on phosphatidylcholine and phosphatidylethanolamine lipids, which are common to many biological membranes. We have exploited this knowledge to develop efficient in vitro lipid-bilayer folding systems for the membrane protein, bacteriorhodopsin. Furthermore, we have shown that a rate-limiting apoprotein folding step and the overall folding efficiency appear to be controlled by particular properties of the lipid bilayer. The properties of interest are the stored curvature elastic energy within the bilayer, and the lateral pressure that the lipid chains exert on the their neighbouring folding proteins. These are generic properties of the bilayer that can be achieved with simple mixtures of biological lipids, and are not specific to the lipids studied here. These bilayer properties also seem to be important in modulating the function of several membrane proteins, as well as the function of membranes in vivo. Thus, it seems likely that careful manipulations of lipid properties will shed light on the forces that drive membrane protein folding, and will aid the development of bilayer folding systems for other membrane proteins.     

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.     

Thylakoid membrane landscape in the sixties: a tribute to Andrew Benson

Prior to the 1960s, the model for the molecular structure of cell membranes consisted of a lipid bilayer held in place by a thin film of electrostatically-associated protein stretched over the bilayer surface: (the Danielli–Davson–Robertson “unit membrane” model). Andrew Benson, an expert in the lipids of chloroplast thylakoid membranes, questioned the relevance of the unit membrane model for biological membranes, especially for thylakoid membranes, instead of emphasizing evidence in favour of hydrophobic interactions of membrane lipids within complementary hydrophobic regions of membrane-spanning proteins. With Elliot Weier, Benson postulated a remarkable subunit lipoprotein monolayer model for thylakoids. Following the advent of freeze fracture microscopy and the fluid lipid-protein mosaic model by Singer and Nicolson, the subunits, membrane-spanning integral proteins, span a dynamic lipid bilayer. Now that high resolution X-ray structures of photosystems I and II are being revealed, the seminal contribution of Andrew Benson can be appreciated.     

Structural transitions in short-chain lipid assemblies studied by (31)P-NMR spectroscopy

The self-assembled supramolecular structures of diacylphosphatidylcholine (diC(n)PC), diacylphosphatidylethanolamine (diC(n)PE), diacylphosphatidyglycerol (diC(n)PG), and diacylphosphatidylserine (diC(n)PS) were investigated by (31)P nuclear magnetic resonance (NMR) spectroscopy as a function of the hydrophobic acyl chain length. Short-chain homologs of these lipids formed micelles, and longer-chain homologs formed bilayers. The shortest acyl chain lengths that supported bilayer structures depended on the headgroup of the lipids. They increased in the order PE (C(6)) < PC (C(9)) < or = PS (C(9) or C(10)) < PG (C(11) or C(12)). This order correlated with the effective headgroup area, which is a function of the physical size, charge, hydration, and hydrogen-bonding capacity of the four headgroups. Electrostatic screening of the headgroup charge with NaCl reduced the effective headgroup area of PS and PG and thereby decreased the micelle-to-bilayer transition of these lipid classes to shorter chain lengths. The experimentally determined supramolecular structures were compared to the assembly states predicted by packing constraints that were calculated from the hydrocarbon-chain volume and effective headgroup area of each lipid. The model accurately predicted the chain-length threshold for bilayer formation if the relative displacement of the acyl chains of the phospholipid were taken into account. The model also predicted cylindrical rather than spherical micelles for all four diacylphospholipid classes and the (31)P-NMR spectra provided evidence for a tubular network that appeared as an intermediate phase at the micelle-to-bilayer transition. The free energy of micellization per methylene group was independent of the structure of the supramolecular assembly, but was -0.95 kJ/mol (-0.23 kcal/mol) for the PGs compared to -2.5 kJ/mol (-0.60 kcal/mol) for the PCs. The integral membrane protein OmpA did not change the bilayer structure of thin (diC(10)PC) bilayers.     

CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

Protein modulation of lipids, and vice-versa, in membranes

This review describes: (i) perturbations of the membrane lipids that are induced by integral membrane proteins, and reciprocally, (ii) the effects that the lipids may have on the function of membrane-associated proteins. Topics of the first category that are covered include: stoichiometry and selectivity of the first shell of lipids associated at the intramembranous perimeter of transmembrane proteins; the chain configuration and exchange rates of the first-shell lipids; the effects of transmembrane peptides on transbilayer movement of lipids (flip-flop); the effects of membrane proteins on lipid polymorphism and formation of non-lamellar phases; and the effects of hydrophobic mismatch on lipid chain configuration, phase stability and selectivity of lipid-protein association. Topics of the second category are: the influence of lipid selectivity on conformational changes in the protein; the effects of elastic fluctuations of the lipid bilayer on protein insertion and orientation in membranes; the effects of hydrophobic matching on intramembrane protein-protein association; and the effects of intrinsic lipid curvature on membrane integration, oligomer formation and activity of membrane proteins.     

CGDB: a database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.