Molecular simulation of the effect of cholesterol on lipid-mediated protein-protein interactions

Experiments and molecular simulations have shown that the hydrophobic mismatch between proteins and membranes contributes significantly to lipid-mediated protein-protein interactions. In this article, we discuss the effect of cholesterol on lipid-mediated protein-protein interactions as function of hydrophobic mismatch, protein diameter and protein cluster size, lipid tail length, and temperature. To do so, we study a mesoscopic model of a hydrated bilayer containing lipids and cholesterol in which proteins are embedded, with a hybrid dissipative particle dynamics-Monte Carlo method. We propose a mechanism by which cholesterol affects protein interactions: protein-induced, cholesterol-enriched, or cholesterol-depleted lipid shells surrounding the proteins affect the lipid-mediated protein-protein interactions. Our calculations of the potential of mean force between proteins and protein clusters show that the addition of cholesterol dramatically reduces repulsive lipid-mediated interactions between proteins (protein clusters) with positive mismatch, but does not affect attractive interactions between proteins with negative mismatch. Cholesterol has only a modest effect on the repulsive interactions between proteins with different mismatch.     

Wetting and capillary condensation as means of protein organization in membranes.

Wetting and capillary condensation are thermodynamic phenomena in which the special affinity of interfaces to a thermodynamic phase, relative to the stable bulk phase, leads to the stabilization of a wetting phase at the interfaces. Wetting and capillary condensation are here proposed as mechanisms that in membranes may serve to induce special lipid phases in between integral membrane proteins leading to long-range lipid-mediated joining forces acting between the proteins and hence providing a means of protein organization. The consequences of wetting in terms of protein aggregation and protein clustering are derived both within a simple phenomenological theory as well as within a concrete calculation on a microscopic model of lipid-protein interactions that accounts for the lipid bilayer phase equilibria and direct lipid-protein interactions governed by hydrophobic matching between the lipid bilayer hydrophobic thickness and the length of the hydrophobic membrane domain. The theoretical results are expected to be relevant for optimizing the experimental conditions required for forming protein aggregates and regular protein arrays in membranes.     

Exploring the membrane proteome—Challenges and analytical strategies

The analysis of proteins in biological membranes forms a major challenge in proteomics. Despite continuous improvements and the development of more sensitive analytical methods, the analysis of membrane proteins has always been hampered by their hydrophobic properties and relatively low abundance. In this review, we describe recent successful strategies that have led to in-depth analyses of the membrane proteome. To facilitate membrane proteome analysis, it is essential that biochemical enrichment procedures are combined with special analytical workflows that are all optimized to cope with hydrophobic polypeptides. These include techniques for protein solubilization, and also well-matched developments in protein separation and protein digestion procedures. Finally, we discuss approaches to target membrane–protein complexes and lipid–protein interactions, as such approaches offer unique insights into function and architecture of cellular membranes.     

Lipids do influence protein function-the hydrophobic matching hypothesis revisited

A topical review of the current state of lipid-protein interactions is given with focus on the physical interactions between lipids and integral proteins in lipid-bilayer membranes. The concepts of hydrophobic matching and curvature stress are revisited in light of recent data obtained from experimental and theoretical studies which demonstrate that not only do integral proteins perturb the lipids, but the physical state of the lipids does also actively influence protein function. The case of the trans-membrane water-channel protein aquaporin GlpF from E. coli imbedded in lipid-bilayer membranes is discussed in some detail. Numerical data obtained from Molecular Dynamics simulations show on the one side that the lipid bilayer adapts to the channel by a hydrophobic matching condition which reflects the propensity of the lipid molecules for forming curved structures. On the other side, it is demonstrated that the transport function of the channel is modulated by the matching condition and/or the curvature stress in a lipid-specific manner.     

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.     

Characterizing the binding of annexin V to a lipid bilayer using molecular dynamics simulations

Annexins play critical roles in membrane organization, membrane trafficking and vesicle transport. The family members share the ability to bind to membranes with high affinities, but the interactions between annexins and membranes remain unclear. Here, using long‐time molecular dynamics simulations, we provide detailed information for the binding of an annexin V trimer to a POPC/POPS lipid bilayer. Calcium ions function as bridges between several negatively charged residues of annexin V and the oxygen atoms of lipids. The preferred calcium‐bridges are those formed via the carboxyl oxygen atoms of POPS lipids. H‐bonds and hydrophobic interactions formed by several critical residues have also been observed in the annexin‐membrane interface. The annexin‐membrane binding causes small changes of annexin trimer structures, while has significant effects on lipid bilayer structures. The lipid bilayer shows a bent shape and forms a concave region in the annexin‐membrane interaction interface, which provides an atomic‐level evidence to support the view that annexins could disturb the stability of lipids and bend membranes. This study provides insights into the commonly occurring PS‐dependent and calcium‐dependent binding of proteins to membranes. Proteins 2014; 82:312–322. © 2013 Wiley Periodicals, Inc.     

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.     

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.     

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.     

Lipid-protein interactions and the function of the Ca2+-ATPase of sarcoplasmic reticulum

Regardless of the nature of the protein constituents of membranes, the molecular arrangement of lipids interacting with them must satisfy hydrophobic, ionic, and steric requirements. Biological membranes have a great diversity of lipid constituents, and this diversity might have functional roles. It has been proposed, for example, that the hydrophobic regions of membrane proteins are stabilized in the membrane through interactions with lipids able to adopt configurations other than the bilayer structure. Progress in understanding at the molecular level how lipid-protein interactions control the properties of membrane proteins has been hindered by the lack of information concerning the structure of the hydrophobic regions of membrane proteins. Nevertheless, there are many examples in the literature describing how changes in the lipid environment affect physical and biochemical properties of membrane proteins. From these studies, discussed in this review, an overall picture of how lipids and proteins interact in membranes is beginning to emerge.     

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.     

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.     

Folding and assembly of beta-barrel membrane proteins

Beta-barrel membrane proteins occur in the outer membranes of Gram-negative bacteria, mitochondria and chloroplasts. The membrane-spanning sequences of beta-barrel membrane proteins are less hydrophobic than those of alpha-helical membrane proteins, which is probably the main reason why completely different folding and membrane assembly pathways have evolved for these two classes of membrane proteins. Some beta-barrel membrane proteins can be spontaneously refolded into lipid bilayer model membranes in vitro. They may also have this ability in vivo although lipid and protein chaperones likely assist with their assembly in appropriate target membranes. This review summarizes recent work on the thermodynamic stability and the mechanism of membrane insertion of beta-barrel membrane proteins in lipid model and biological membranes. How lipid compositions affect folding and assembly of beta-barrel membrane proteins is also reviewed. The stability of these proteins in membranes is not as large as previously thought (<10 kcal/mol) and is modulated by elastic forces of the lipid bilayer. Detailed kinetic studies indicate that beta-barrel membrane proteins fold in distinct steps with several intermediates that can be characterized in vitro. Formation of the barrel is synchronized with membrane insertion and all beta-hairpins insert simultaneously in a concerted pathway.     

Contributions of Gaussian curvature and nonconstant lipid volume to protein deformation of lipid bilayers

An elastic model for membrane deformations induced by integral membrane proteins is presented. An earlier theory is extended to account for nonvanishing saddle splay modulus within lipid monolayers and perturbations to lipid volume proximal to the protein. Analytical results are derived for the deformation profile surrounding a single cylindrical protein inclusion, which compare favorably to coarse-grained simulations over a range of protein sizes. Numerical results for multi-protein systems indicate that membrane-mediated interactions between inclusions are strongly affected by Gaussian curvature and display nonpairwise additivity. Implications for the aggregation of proteins are discussed.     

Interactions between misfolded protein oligomers and membranes: A central topic in neurodegenerative diseases?

The deposition of amyloid material has been associated with many different diseases. Although these diseases are very diverse the amyloid material share many common features such as cross-β-sheet structure of the backbone of the proteins deposited. Another common feature of the aggregation process for a wide variety of proteins is the presence of prefibrillar oligomers. These oligomers are linked to the cytotoxicity occurring during the aggregation of proteins. These prefibrillar oligomers interact extensively with lipid membranes and in some cases leads to destabilization of lipid membranes. This interaction is however highly dependent on the nature of both the oligomer and the lipids. Anionic lipids are often required for interaction with the lipid membrane while increased exposure of hydrophobic patches from highly dynamic protein oligomers are structural determinants of cytotoxicity of the oligomers. To explore the oligomer lipid interaction in detail the interaction between oligomers of α-synuclein and the 4th fasciclin-1 domain of TGFBIp with lipid membranes will be examined here. For both proteins the dynamic species are the ones causing membrane destabilization and the membrane interaction is primarily seen when the lipid membranes contain anionic lipids. Hence the dynamic nature of oligomers with exposed hydrophobic patches alongside the presence of anionic lipids could be essential for the cytotoxicity observed for prefibrillar oligomers in general. This article is part of a Special Issue entitled: Lipid–protein interactions.     

Direct measurement of the mechanical properties of lipid phases in supported bilayers

Biological membranes define not only the cell boundaries but any compartment within the cell. To some extent, the functionality of membranes is related to the elastic properties of the lipid bilayer and the mechanical and hydrophobic matching with functional membrane proteins. Supported lipid bilayers (SLBs) are valid biomimetic systems for the study of membrane biophysical properties. Here, we acquired high-resolution topographic and quantitative mechanics data of phase-separated SLBs using a recent atomic force microscopy (AFM) imaging mode based on force measurements. This technique allows us to quantitatively map at high resolution the mechanical differences of lipid phases at different loading forces. We have applied this approach to evaluate the contribution of the underlying hard support in the determination of the elastic properties of SLBs and to determine the adequate indentation range for obtaining reliable elastic moduli values. At ~200 pN, elastic forces dominated the force-indentation response and the sample deformation was <20% of the bilayer thickness, at which the contribution of the support was found to be negligible. The obtained Young's modulus (E) of 19.3 MPa and 28.1 MPa allowed us to estimate the area stretch modulus (k(A)) as 106 pN/nm and 199 pN/nm and the bending stiffness (k(c)) as 18 k(B)T and 57 k(B)T for the liquid and gel phases, respectively.     

Lateral pressure profile, spontaneous curvature frustration, and the incorporation and conformation of proteins in membranes

Lipid-protein interactions are an important determinant of the stability and function of integral and transmembrane proteins. In addition to local interactions at the lipid-protein interface, global interactions such as the distribution of internal lateral pressure may also influence protein conformation. It is shown here that the effects of the membrane lateral pressure profile on the conformation or insertion of proteins in membranes are equivalent to the elastic response to the frustrated spontaneous curvature, c_o, of the component lipid monolayer leaflets. The chemical potential of the protein in the membrane is predicted to depend linearly on the spontaneous curvature of the lipid leaflets, just as does the contribution of the protein to the elastic bending energy of the lipid, and to be independent of the hydrophobic tension, γ_phob, at the lipid-water interface. Analysis of the dependence of protein partitioning or conformational transitions on spontaneous curvature of the constituent lipids gives an experimental estimate for the cross-sectional intramembrane shape of the protein or its difference between conformations. Values in the region of 50-110 Å² are estimated for the effective cross-sectional shape changes on the insertion and conductance transitions of alamethicin, and on the activation of CTP:phosphocholine cytidylyltransferase or rhodopsin in lipid membranes. Much larger values are estimated for the mechanosensitive channel, MscL. Values for the change in intramembrane shape may also be used, together with determinations of lipid relative association constants, to estimate contributions of direct lipid-protein interactions to the lateral pressure experienced by the protein. Changes in chemical potential ∼12 kJ mol⁻¹ can be estimated for radial changes of 1 Å in a protein of diameter 40 Å.